基于非线性建模的水电机组动态特性仿真研究_魏加达.pdf

12第 46 卷 第 02 期2023 年 02 月Vol.46No.02Feb.2023水 电 站 机 电 技 术Mechanical&ElectricalTechniqueofHydropowerStation1 引言水轮机调节系统作为受水机电三大因素耦合作用的复杂非线性系统，是保证水轮发电机组安全稳定运行的重要控制设备1。搭建水轮机调节系统模型，进行仿真计算可为发电机组参数设计、导叶关闭规律整定、调速系统参数辨识以及模式切换等提供依据，并将有利于改善大小波动过渡过程品质以及一次调频动态响应特性等2-4。当前，一方面，大型机组突增突减负荷现象频繁出现，水电机组大波动过渡过程动态特性对电网安全稳定的影响日益显现；另一方面，新能源建设的蓬勃发展、多类型电源互相协调的发展结构对水电机组一次调频性能提出了新的挑战，因此建立一个准确可行的数学模型对机组动态特性加以分析研究具有重大意义。由于水力因素的不稳定性，特别是受水流惯性、流量时变性因素影响，以及水轮机功率调节非线性、发电机非线性特性以及电液随动系统中主配压阀输出饱和非线性、主配压阀死区非线性和主接力器速度限幅、主接力器输出饱和非线性等众多因素影响，使得想建立一个反映水轮机调节系统动态特性的精确数学模型变得极为复杂与困难，需研究者进行大量有针对性的细节研究与处理。此前，调节系统相关建模仿真应用多为静态特性研究、参数辨识等5-8，缺乏对于机组动态特性的研究；在利用调节系统整体模型进行研究时，多采用理想水轮机模型与随动系统一阶模型9-11，未针对随动系统进行专门深入研究。本文综合考虑调节系统各组成部分的非线性因素，建立了一个高精度反映调节系统真实动态特性的非线性模型。基于所建立的模型以大型水利枢纽 TG 电站为例，进行了甩负荷大波动过渡过程与一次调频动态响应过程仿真，证明了文中所建立的非线性模型的能够高精度、准确反映机组的甩负荷以及一次调频动态特性，从而为机组甩负荷过渡过程品质及一次调频性能优化研究提供了重要技术手段。2 水轮机调节系统非线性建模2.1 调节器模型调节器数学模型的建立与其控制规律息息相关，目前应用广泛的并联 PID 控制数学模型如图 1所示。图 1 中 f 为机组频率；Cf为频率给定；死区为频率死区；Cy为开度给定；KP1为开度模式下比例增收稿日期：2022-09-13作者简介：魏加达（1996-），男，助理工程师，从事水电厂自动装置诊断分析工作。基于非线性建模的水电机组动态特性仿真研究魏加达，王卫玉，张 培（湖南五凌电力科技有限公司，湖南 长沙 410000）摘 要：考虑到大型机组突增突减负荷的频繁出现以及新能源不稳定性对水电机组一次调频带来的挑战，建立能准确反映机组动态性能的仿真模型显得尤为重要。水轮机调节系统具有非最小相位、强非线性特点，需对其包含的大量非线性因素进行针对性研究和处理。针对水轮机调节系统各组成部分，综合考虑各因素建立了相应数学模型，重点建立了随动系统精细化模型、水轮机外特性数学模型，最终建立了水轮机调节系统整体非线性数学模型。基于建立的模型，进行甩负荷及一次调频仿真，通过与真机数据进行对比分析，证明了模型可高精度、准确反映机组甩负荷及一次调频动态特性，为水电机组动态特性研究提供了新的技术手段。关键词：水轮机调节系统；非线性建模；仿真；甩负荷；一次调频中图分类号：TK730文献标识码：A文章编号：1672-5387（2023）02-0012-06DOI：10.13599/ki.11-5130.202302.00413第 02 期魏加达，等：基于非线性建模的水电机组动态特性仿真研究益；KI1为开度模式下积分增益；Kd1为开度模式下微分增益；Td1为开度模式下微分衰减时间常数；bp为永态转差系数；死区为开度死区；p 为机组功率；Cp为功率给定；死区为功率死区；KP2为功率模式下比例增益；KI2为功率模式下积分增益；Kd2为功率模式下微分增益；Td2为功率模式下微分衰减时间常数；ep为调差率。开度模式传递函数为：G1（s）=Y（s）X（s）=-Kds2+Kps+KibpKds2+（bpKp+1）s+bpKi （1）功率模式传递函数为：G2（s）=Y（s）X（s）=-Kds2+Kps+Kieps （2）2.2 随动系统精细化模型典型的调速器电液随动系统主要由综合放大器、电气机械/液压转换装置、主配压阀和主接力器等部分构成。电液随动系统包含诸多非线性因素，譬如液压随动系统中的主配压阀死区非线性、主配压阀输出饱和非线性和主接力器速度限幅、主接力器输出饱和非线性等，对系统动静态特性存在较大影响。大量的仿真调试经验表明，采用简化的一阶线性模型进行仿真分析得出的参数与实际相差较远。由此造成了一系列问题，比如为保证系统的动态稳定性，只能大幅度调小综合放大器的放大倍数或调大内环反馈增益，结果使得难以保证系统静态精度，且减缓了随动系统的响应速度，导致随动系统输出滞后于控制器输出，恶化了动态品质。因此，在进行调节系统建模时，必须将随动系统一系列非线性因素加以考虑，建立精细化模型。综合考虑各类死区、饱和限幅和开关机时间分段调整等非线性因素，建立两段开启、三段关闭的随动系统非线性精细化数学模型，如图 2 所示。图 2 随动系统二段开启、三段关闭模型图 2 中 KA为综合放大器放大倍数，1Ty1s+1为第一级液压放大环节的传递函数；限幅为主配压阀限幅；死区为主配压阀搭叠量造成的死区，Ty为接力器反应时间常数；限幅为接力器速度对应的限幅；限幅为接力器行程限幅。对于不同的开关闭特性相应模型，分别以不同开启速度 1/To1、1/To2，不同的关闭速度 1/Tc1、1/Tc2、1/Tc3等对接力器开关闭速度进行整定，实现两段开启与三段关闭。2.3 外特性水轮机模型水轮机外特性一般由水轮机模型综合特性曲线与飞逸特性曲线描述，是水轮机内部流动规律的外在表现，可以表示水轮机在不同工况下的能量转换特性。两种曲线均可以通过模型试验或现场试验获取。由水轮机单位参数公式、出力公式及力矩公式知：n11=nD1 H （3）Q11=QD2 H （4）P=9.81QH （5）M=P=P2n/60 （6）由上述式（3）（6）可推出单位力矩公式：M11=609.812Q11n11 （7）因此通过读取模型综合特性上不同单位流量、不同等开度线下的效率值则可计算得到水轮机的力矩特性，流量特性则可直接读取得到。但模型试验得到的模型综合特性曲线只涵盖水轮机运行的高效区，不包括空载、开机等过渡工况区域，因此通常采用特征曲线外延的方法来补全水轮机全工域特性，以此得到完整的水轮机流量特性与力矩特性，来表征水轮机动态特性。本文基于模型综合特性曲线与飞逸特性曲线等外特性曲线，建立外特性水轮机模型来描述水轮机流量特性与力矩特性随转速、导叶图 1 并联 PID 调速器控制逻辑框图14第 46 卷水 电 站 机 电 技 术开度的动态变化关系如图 3 所示。图 3 外特性水轮机模型图 3 中 y 为随动系统输出的导叶开度；n 为机组转速；n11为单位转速；n11r为额定单位转速；q 为机组流量；q11为单位流量；h 为水头；mt为机组输出力矩；mt11为额定单位力矩。在进行水轮机全工域特性补全时，主要以飞逸特性曲线为限制线，以模型综合特性曲线获取的流量、力矩数值为基础，采用恰当的数值拟合方法来补全小开度区域的数据。最终得到的水轮机流量特性与力矩特性如图 4、图 5 所示。2.4 引水系统模型由于引水系统结构多样化和水流状况的复杂性，引水系统蕴含着丰富的非线性动态特性，建立引水系统精确模型存在较大困难，且与水机电耦合求解时过于复杂。调节系统建模时大多依据引水管道的长度来选择采用刚性或弹性水击模型，本文引水系统采用弹性水击模型进行描述，其传递函数为：G（hq）=-2hweTr2s-e-Vr2seTr2s+e-Tr2s （8）式中：hw=Tw/Tr，为水轮机管道特性系数；TW为水流惯性时间常数，s；Tr为水击波相长，s；q 和 h 分别为水轮机流量和水轮机工作水头标幺值的拉氏变换；s 为拉普拉斯算子。2.5 发电机模型发电机作为机械能转化为动能的核心设备，其性能同时受多种因素影响，包括发电机本身作为旋转刚体具有的机械惯性、发电机功率与关系、励磁系统以及 d、q 轴绕组阻尼等。电站水轮机调节系统动态性能仅受发电机转速变化过程影响，将发电机作为旋转刚体考虑即可满足仿真计算精度要求，采用一阶模型描述，传递函数为：G（s）=1Tas+en （9）式中：Ta为机组惯性时间常数；en=eg-ex，其中ex为水轮机力矩对转速的偏导数，eg为发电机力矩对转速的偏导数；en计算时只考虑负载与电网的影响，水轮机的 ex已反映在水轮机模型中；s 为拉氏算子。（a）单位流量等开度线（b）单位流量曲面图 4 水轮机流量特性（a）单位力矩等开度线（b）单位力矩曲面图 5 水轮机力矩特性15第 02 期魏加达，等：基于非线性建模的水电机组动态特性仿真研究3 机组动态特性仿真以某大型水利枢纽电站为例，基于上述建立的调节系统非线性模型，进行机组甩负荷、一次调频等动态特性仿真。根据提供的资料，电站水轮机采用混流式水轮机，机组的相关参数如表 1 所示。表 1 机组参数表名称单位大小转轮出口直径 Dm9.44最大水头 Hmaxm113.0额定水头 Hrm85.0最小水头 Hminm65.0额定流量 Qrm3/s941.27额定出力 PrMW700额定转速 nrr/min75.0水流惯性时间常数 Tws2.0发电机转动惯量 GD2t m2450 0003.1 机组甩负荷动态过程仿真水轮机甩负荷作为水轮发电机组在日常运行中常见的大波动过渡过程，机组甩负荷过程中会造成最大转速上升值和最大压力上升值等参数的急剧变化，若设计不当，将引发安全事故，严重威胁电站安全稳定运行，因此对于机组甩负荷工况加以研究对电站安全运行具有重要指导意义。根据电站提供的实测甩负荷数据，分别验证甩小负荷 175 MW（25%Pr）和大负荷 756 MW（108%Pr）。设定负荷相应的接力器行程与导叶开度百分数，调整仿真初始工况出力稳定后，设定适当的仿真时长，进行甩负荷。待仿真结束，读取仿真计算结果，将得到的仿真结果绘制曲线与实测曲线对比，如图 6、图 7所示。1：14322：3：4：图 6 机组甩 175 MW 负荷1：14322：3：4：图 7 机组甩 756 MW 负荷由图 6、图 7 中转速变化曲线对比知：仿真转速上升速度略快于实测转速上升，达到转速上升峰值的时间早，仿真峰值 105.2 s/139.3 Hz 略小于实测峰值 106.6 s/140.8 Hz，峰值误差最大不超过 1.1%，仿真转速与实测转速下降过程趋势基本一致，吻合得较好。由图 6、图 7 中导叶开度变化曲线对比可知：当机组甩 756 MW 负荷时，两者导叶关闭规律吻合较好；当机组甩 175 MW 负荷时，导叶开度变化吻合性较差。误差产生的原因可能是缺乏小开度的飞逸特性曲线，采用大开度外延的小开度飞逸与实际情况有一定的偏差，插补的制动区力矩特性与实际值存在差异导致转速峰值与转速下降速度存在差异。通过上述机组甩大、小负荷过程分析可知，所建立的模型能够较为准确反映机组甩负荷过渡过程中的动态特性。3.2 一次调频动态过程仿真根据电站提供的真机一次调频试验数据，分别进行开度模式和功率模式下的一次调频仿真验证。理论上，在开度模式下，主要考虑开度变化曲线的吻合度；在功率模式下，主要考虑功率变化曲线的吻合度。依据仿真计算结果绘制导叶开度变化曲线与功率响应过程曲线对比图，具体分析如下所述。3.2.1 开度模式下一次调频仿真（1）一次调频过程 1仿真条件为：H=102.6 m，P=603.3 MW,接力器行程百分数为 47.50%，调速器参数 kp=3.5、ki=0.95、kd=0，永态转差系数Bp=0.04，频率扰动f=+0.15 Hz，开度和功率响应过程对比如图 8 所示。16第 46 卷水 电 站 机 电 技 术1：2：21(a)导叶开度变化