中国水能蕴藏量1万kW以上的河流有300多条,水能资源丰富程度居世界第一.水资源总量约2.8万亿m³,约占全世界水资源总量的6%,可开发量为3.78亿kW.和世界其它国家相比中国水能利用情况处于较低水平.随着煤炭等不可再生资源的过度开采,水能这种绿色、可持续的发电能源必将得到大力开发,大型水电站也必将相继投入运行,对运行人员的有效培训也将进一步加强,这就需要建立更加完善的模型来建立水电站的仿真培训系统^[1].同时,水电站孤网运行,水电厂和风电、火电的联调也需要建立能够更好地模拟大波动过程的仿真模型.

本文通过分析水电站调节系统存在的非线性因素,对相关部分分别建立非线性模型,将各部分组合,得到完整的调节系统非线性模型.在Matlab/Simulink软件中建立其模型,验证模型的准确性.

水轮机调节系统包括调速器和调节对象,如图 1所示.电厂常用的调速器是微机调速器,调节对象主要包括水轮机及其引水系统和发电机.

为了简化分析,以混流式水轮机为例.水轮机动态特性指调节过程中水轮机力矩M_t、流量Q随导叶开度α、水头H和转速n变化的特性^[2],即

图 1

Fig. 1

图 1 水轮机调节系统结构图 Fig. 1 Turbine regulating system structure

非线性水轮机模型中考虑了在大波动过程中机组变化比较大,传递函数不能恒定为常数的特点.水轮机力矩和流量用积分形式表达,即

机组执行机构为数字控制,由综合放大环节、电液伺服环节(比例阀)、配压阀以及主接力器等构成.建立更加符合现场实际的模型需注意:

(1) 考虑现场实际中的非线性情况:转速、导叶开度均有实际的物理意义和取值范围;并加入了限速、限幅的饱和非线性环节.

(2) 考虑到实际情况中,开启和关闭时接力器的机械动作特点不同,开启和关闭时接力器反应时间常数也不同.所以建立液压随动系统模型时要考虑这一点,开启和关闭时选择不同的接力器反应时间常数.

在常用机械液压随动系统基础上,结合以上两点,建立的液压随动系统模型如图 2所示.其中:u为输入变量;t为时间;Kp、Ki、Kd均为调节器参数;To、Tc分别为开启和关闭的时间常数;T₂为惯性时间常数;P_max、Pmin为速度的上限和下限;s为复变量.

图 2

Fig. 2

图 2 液压随动系统模型 Fig. 2 Hydraulic servo system model

水轮机是一个复杂的时变非线性系统,目前还没有公认的表达式可描述它的流量特性和力矩特性.本文根据水力机械的主要参数及模型单位参数间的关系,搭建非线性水轮机模型,其中无法用数学模型表达的关系就借助水轮机转轮特性曲线查表得到^{[3, 4]},则有

在仿真分析中,变量均以相对偏差值表示,但是在水轮机特性表中变量的形式是相对值或全量值,所以在查表Q₁₁=f(n₁₁,α)和表η=f(n₁₁,α)之前要将相对偏差值转换为全值量^{[5， 6]}.变量对应关系如表 1所示.根据式(3)和表 1可在Matlab/Simulink软件中得到水轮机非线性仿真模型,如图 3所示.其中:h₀、n₀、y₀分别为水头、转速、接力器行程的初始值.通过输入导叶开度和单位转速查表得到相应的水轮机流量和机组效率.

图 3

Fig. 3

图 3 流量特性和力矩特性模型 Fig. 3 Flow characteristic and torque characteristic model

表 1中,n_r、Q_r、H_r、M_r分别为额定工况下的转速、流量、水头和力矩;Q₀、H₀、M₀、Y₀分别为初始工况下的流量、水头、力矩和接力器行程;Y、Y_m分别为接力器行程和接力器行程的最大值;x、q、h、m、y分别为转速、流量、水头、力矩、接力器行程的偏差相对值.

表 1(Table 1)

表 1 变量对应关系 Tab.1 Correspondence between variables 参数名称 偏差相对值 全量值 转速 x=n-n0 nr n=n0+xnr 流量 q=Qt-Q0 Qr Q=Q0+qQr 水头 h=H-H0 Hr H=H0+hHr 力矩 m=Mt-M0 Mr Mt=M0+mMr 接力器行程 y=Y-Y0 Ym Y=Y0+yYm

表 1 变量对应关系 Tab.1 Correspondence between variables

3 水轮机调节系统非线性仿真模型

水电站的调速器都是选用PID调节模块,但是考虑到现场频率的反馈信号与发电机出口的频率信号相比存在延迟和滤波,信号在传输过程中有一些变化,所以结合某水电站机组的调速器模型在反馈信号中加入一阶微分环节,模拟这种影响.得到的水轮机调节系统非线性模型如图 4所示.

图 4

Fig. 4

图 4 水轮机调节系统非线性模型 Fig. 4 Nonlinear turbine regulating system model

选定云南省某水电站机组作为对象进行仿真.在水头为169 m,导叶开度为92%,带有功功率为73%时进行甩负荷仿真实验,控制器的参数设置为:主环PID参数中调节器参数分别为Kp=4,Ki=1,Kd=0.1,永态转差系数b_p=0.040 4;副环PID参数中调节器参数分别为K_p=8.38,Ki=0,Kd=0,接力器动作时间T_yT_o=18.74 s;被控对象模型参数中水流惯性时间常数T_w=1.6,发电机惯性时间常数Ta=10,发电机负荷自调节系数en=1.甩负荷仿真实验结果如图 5所示.

图 5

Fig. 5

图 5 甩负荷仿真实验结果 Fig. 5 Load rejection simulation results

甩负荷过程中系统波动较大,这是由于甩负荷时要求频率保持在额定值(50 Hz),所以发电机转速和频率的变化幅度较大,运用非线性模型能够更好地模拟这一过程.但是由于本文建立的非线性模型用到的二维查表法是一线性插值方法,其精度不高,因此该模型还需要进一步优化、完善.

对水电站调节系统存在的非线性因素进行了分析.首先,水轮机是调节系统非线性的最大来源,所以本文在水轮机综合特性曲线中水轮机特性关系的基础上建立了水轮机非线性模型;其次,由于信号反馈中存在惯性和延迟,阀门开启和关闭等动作存在死区和限幅,因此增加了相应的非线性环节.

在Matlab/Simulink软件中搭建仿真模型,各部分的传递函数和参数设定参考了云南某水电站650 MW机组的实际情况.通过仿真实验证明了模型的准确性,为今后建立大机组水轮机调节系统准确、可靠模型提供了参考.