冲击式水轮机甩负荷水击反分析与控制方法研究
北极星电力网技术频道  
作者:电力论文3
2007/12/24 18:41:16
【摘要】通过对冲击式水轮机甩负荷水击问题的反分析研究,结合PID控制理论,总结出一套如何确定阀门关闭规律的方法。这对于冲击式水轮机水电站的稳定工作具有相当大的实用价值。图9幅。
【关键词】冲击式水轮机水击反分析控制
1前言
在水力发电站的运行中,不可避免地要涉及到水轮机因各种工作情况而进行的负荷变更,电能的生产与其它部门的生产不同(电能的生产、分配和消费是在同一个时间内进行的),故其不能大量地储存。电力系统中发电厂供电量的多少,决定于用户的需要;因此,从发电厂到用户的各个环节中,任何一个环节的故障或其运行方式和用电方式的改变均将影响整个电力系统中供电和用电的平衡性。这种影响在水轮机上的体现就是各种工况的甩负荷,而在水电站的压力引水系统中的体现则为水击现象。如何通过各种甩负引起的水击测量值,反分析出阀门关闭的规律就变得更加具有实际意义,即可以通过阀门的关闭规律来控制水击,这一问题称为水击问题的反分析。
2 水击反分析理论
在给定水击压力允许峰值以后,确定的相关流速初始值,应满足数学模型:
求:V(x,0)=Vo。其中,g为重力加速度;H为水头;y为速度;f为摩擦因子;D为管道直径;α为管线与水平面的倾角。
2.1数值解法
我们构造出一类判别泛函,在假定水锤正问题存在唯一稳定解的条件下,根据判断函数的特性,可确定反问题的解的存在范围和具体求解方法。
定义1.设,水锤正问题存在唯一稳定解,而且Vo在某一物理范围,内变化,即Vm≤Vo≤VM,则称为反问题的允许流速初始值。
定义2.设Hc为给定的水击允许峰值,在(1)的允许流速初始值类中任取一个Vo,其对应的水击最大峰值为HV(Vo),则反问题(1)的判别函数为:
注意到HV(Vo)是Vo的单调增连续函数,很容易证明有下面这个定理:
定理:设,是(1)的允许流速初值类,如果一切,均有μv<0,则系统绝对安全;如果一切,均有μv>0,则系统绝对不安全。除了以上两种情况外,(1)有唯一解。从流体力学的知识可以知道,流速与流量之间为正比例关系,通过这一关系可以确定管道截面上流量的变化关系。
2.2阀门边界处理
通过以上关系的确定,可以反分析出阀门因流量变化而引起的阀门开度变化。这里我们只以冲击式水轮机的关闭方式为例子说明。由于冲击式水轮机的转速与负荷之间没有直接的联系,因此在处理时只需要考虑喷嘴阀门的开度[1]。
3控制原理
在一个水力系统中,影响水击峰值的先天因素为:管路参数(管长、管截面积),管路的连接方式,部件参数(如阀门参数),操作方式(如阀门关闭方式)及压力波速。后天因素为流速初值,操作时间(如阀门关闭时间)。其中,流速初值、压力波速和操作时间属于可控制因素,其它因素在设计安装时已经确定,一般在运行时不可控制。在后边3个因素中,操作时间是由机械性能或操作技能所决定(个别操作时间,如阀门关闭时间与别的因素有关),所以真正需要探讨其控制原理的因素只有一个:流速的初值。实现对流速初值控制的理论基本上是求其反问题,原则上这类反问题为一切水轮机甩负荷水击的控制问题提供了理论基础。
4 PID控制理论
由于阀门关闭过程的分析在现有理论基础上均采用线性分析理论进行研究。PID(Proportional,Integral,Derivative)控制方法作为一种广泛应用于结构时频域内的控制理论,同样可以应用于水轮机甩负荷水击的控制。PID控制的传递方程为:
其中,KP、KI、KD分别为相应的控制参数;s表示时间。
通过控制传递方程的参数调整以达到对最后水击压力峰值的控制,详细的PID控制理论可参阅其它文献。这样通过PID控制的水击压力峰值来寻求阀门关闭的规律,也就为水轮机甩负荷操作提供了理论依据。
5数据分析
实测数据采集地点为云南大华水电站。电站设计装机为2台卧式冲击水轮机,装机容量2×7MW。工程是由位于花坪河和扎局清河的首部取水建筑物、2条引水渠(分别引花坪河河水和扎局清河河水)、压力前池(容积5×104m3)、压力钢管及厂区枢纽等组成。
5.1工况1
2台机同时甩100负荷时2号机阀门处的水击压力峰值:
图12台水轮机同时甩满负荷时2号机阀门处的压力时程
相应PID控制数值仿真计算压力峰值:
图22台水轮机同时甩满负荷时2号机阀门处的数值仿真压力时程
5.2工况2
2台机同时甩50负荷时2号机阀门处的实测压力峰值:
图32台水轮机同时甩50荷时2号机阀门处的压力时程
相应PID控制数值仿真计算压力峰值:
图42台水轮机同时甩50负荷时2号机阀门处的数值仿真压力时程
5.3工况3
2台机同时甩25负荷时2号机阀门处的实测压力峰值:
图52台水轮机同时甩25时2号机阀门处的压力时程
相应PID控制数值仿真计算压力峰值:
图62台水轮机同时甩25负荷时2号机阀门处的数值仿真压力时程
5.4数据说明
HD控制仿真压力峰值的计算是以常规的特征线法为计算基础,在计算过程中通过PID各个参数的控制以达到计算结果与实测压力峰值相当的水平,一般XP参数通常取为50~100。对于冲击式水轮机而言,水轮机的负荷完全通过喷针开度来调节,这样通过数值仿真而得到的水力峰值就可以作为阀门开度反分析的依据,从而得到阀门开度的变化规律。依照这原理,我们得到了阀门的开度控制曲线:
图7工况1阀门开度曲线圈
图8工况2阀门开度曲线图
通过以上阀门开度曲线(图中:1表示完全关闭;0表示完全打开)可以看出,在正常的仿真反分析中,阀门开度变化经常会在某一位置做定幅振动,这在实际操作过程中是根本无法采用的,因此我们建议在关闭过程中采用间歇方式,即在振荡处保持恒定状态。而阀门为何不能最终关闭,其主要原因是我们在建立模型时没有考虑阀门自身的摩擦阻力。