小波动状况下简单式调压室水位波动分析.docVIP

小波动状况下简单式调压室水位波动分析 摘要：计算小波动时，由于调压室水位的波动情况复杂性，需要调用程序去解决，而本文根据其基本方程推导出了解析法，避免了调用繁琐的程序。当水电站发生小波动时，应用泰勒展开式将引水道和调压室的动量方程的摩阻损失线性简化，通过求解常微分方程就可得出流量瞬时改变或流量随时间改变时调压室中的水位波动方程，进而可分析水位波动振幅的衰减和验证调压室的稳定性。 关键词：小波动；调压室；水位波动；振幅 1 前 言 水电站在正常运行时，可能发生较小的突然性负荷变化，水轮机引用流量的改变将引起调压室水位的波动。而这微小的波动可能一经此扰动不会再回到原来的平衡位置，也会引起大量的弃水，同时调压室水位振幅的衰减也将影响到机组的出力，因此小波动状况下的调压室的水位波动也应受到重视。该状况下水位波动的变化能很好的反应过渡过程，在大多数过渡过程中小波动的稳定均能满足大波动的稳定，因此小波动过程中调压室水位的衰减振荡程度能很好的反应其调压室的稳定性，为研究小波动对机组的影响也有重要的意义[1,2]。 目前，小波动状况下调压室水位的波动大多是采用编程程序或专业的软件进行分析，需要调用繁琐的程序或建立复杂的模型。为此其解析公式的推求有着重要的意义，本文采用从简单调压室的基本方程出发，根据小波动状况下流速变化较小的特点，通过对动量方程中摩阻损失简化为速度的线性函数的处理得出小波动状况下调压室的水位波动计算公式。 作者简介：史晓阳(1987–)，男，汉族，河北人，河海大学硕士研究生，研究方向为水利水电工程，E-mail:shixiaoyanghb@163.com 2 小波动下简单式式调压室的水位波动计算 2.1基本方程的建立 如图1，调压室水位波动的基本方程[3]： 连续性方程：fv=F+Q (1) 动量方程： (2) 图1 引水系统示意图 Fig.1 schematic diagram of conduit system 式中：f为引水道的断面面积；v为引水道中的流速；F为调压室的断面面积；Z为以水库静水位为基准的调压室涌波，向下为正，向上为负；L为引水道的长度；为引水道中的水头损失；其中与v2成比例。 2.2 线性振荡近似简化 在动量方程中由于与v2成比例，所以在计算中出现了非线性的水头损失。但对于小波动中的较小波动，可以使用线性函数来简化M=的取值[4]。如图所示 图2 摩阻、流速关系图 Fig.2 relation schema of fraction drag and flow velocity 在较小的水位波动中，其流速v和最终流速v0相差不大，因此可以采用泰勒展开式的前两项近似计算： M(v)=M0+M′ (v0)(v-v0)=M0-kx (3) 2.3 流量瞬时变化方程的求解 引入变量x和y，令x=v-v0，y=M0-Z= -Z。将x，y和式(3)带入式(1)，(2)得到流量瞬时改变而产生的波动方程： F=fx (4) =-y-kx (5) 其中k=2，初值的计算：在t=0时，初始值x1=v1-v0 ，y1=M0-M1= - 将(5)式两边对t求导，然后将(4)带可得： (6) 对于该常系数微分方程，可令x=， (7) 其中 将x带入(5)即可求出函数y。 (1)当n1时，稳定后流速比较大： ,两根均为实数，所以： x= (8) y=+ (9) (2)当n1时，稳定后流速比较小时： ,两根为共轭复根，所以 (10) (11) (3) 当n=1时，所以 x= (12) y=