基于特征线法的液控止回蝶阀消除水锤仿真研究.docxVIP

基于特征线法的液控止回蝶阀消除水锤仿真研究 1 应用案例及仿真的背景 黄河虞门污水处理厂二级水泵。在管道安全方面，采用液体控制器停止回归军阀，避免因事故中断引起的水回流、水锤侵蚀和高速逆转，确保机组和管道的安全运行。该蝶阀由铁岭阀门厂生产,型号为HD741-16 DN600(见图1)。该蝶阀的工作原理是:在水泵起动前,由于流道出口淹没在出水池水位以下,为防止出水池水体倒流,蝶阀处于关闭状态;在水泵起动过程中,随着机组转速的上升,水泵出口水位逐渐上升,流道中的空气从蝶阀前的空气阀排出,经过很短的造压过程,蝶阀开始开启;在水泵正常停泵过程中,待蝶阀快速关至一定角度后机组断电,液控蝶阀继续关闭直至全关,切断水流。本研究结合该供水工程的实际,在液控止回蝶阀安全防护下,对停泵暂态过程进行计算机仿真研究,并通过现场测试进行比较,以期得到液控止回蝶阀防护水锤的最优关闭过程,从而为类似供水工程的建设、运行和开展这一课题的深层次研究提供技术支持。 禹门口二级泵站位于山西、陕西两省交界处,属于高扬程(H=105 m)、大流量(Q=1.6 m3/s)的抽黄灌溉工程,设计灌溉面积3.33万hm2,安装10台机组,单机单管,明敷预应力混凝土管道,内径1.0 m,长550 m。总流量14.4 m3/s,地形扬程98.8 m,河水平均含沙量达50 kg/m3。选用甘肃工业大学设计、阳泉水泵厂生产的单级双吸离心泵,水泵型号为32SH-7,其额定参数如表1所列。 2 现场试验设备和测量方法 2.1 实时量测 基于该提水泵站的实际情况,水泵断电情况下,对泵及液控止回蝶阀水力过渡过程有关参数进行实时量测。内容包括:(1)突然断电后,液控止回蝶阀前后、起坡处和管路变坡处压力随时间变化规律的实时量测;(2)机组断电后其转速随时间的变化规律;(3)实时量测断电后主管流量随时间的变化规律。 2.2 试验设备和测量方法 2.2.1 压力和时间检测 通过球型阀门将主管与ZQ-Y型压力传感器互联,通过有线通信方式传至计算机采集系统记录相关的数据,求出相应的压力与时间过程。为校核,在试验现场的每个压力测点上,均有专人进行目测。 2.2.2 锚点两侧导向系统 利用测速电机进行顺势转速量测,将发电机用支架固定在主轴旁,用一干电池组供给测速发电机励磁,测速电机的输出电压与所测转速成正比例输出,求得水泵的瞬态转速和对应的时间过程。 2.2.3 主管流量的计算 分支管测流法,包括支管、闸阀及差压变速器(输出标准直流信号4~20 mA),送至计算机可实时记录支管流量随时间的变化规律。通过流体力学理论“并联管端的水头损失相等”特点,可求得下列关系式:Q=kq,计算便可求得主管流量。式中,q为分支管测出的瞬态流量(m3/s);Q为主管流量(m3/s);k为测流装置的系统流量倍数。现场试验开始前,测流装置在水泵稳定工况下利用超声波流量仪进行精确标定。该供水工程的水力过渡工程的实测结果见表2,限于篇幅表中仅给出最大最小压力包络线数据。 3 理论计算与实验结果表明 3.1 特征线法的改进 水锤图解法对于布置复杂和摩阻损失占比重较大的长管道系统,计算精度较差。美国学者斯利特和怀利提出的特征线法,可考虑管道的摩阻损失及水锤方程的其它次要项,其计算精度高、速度快、收敛性好。本研究在基于对该法深刻理解的基础上,建立了考虑管路摩阻时压力管路水锤的基本方程,原理如下。 3.1.1 求一个管道断面的流量和流量 设i-1、i和i+1为管道上3个相邻的断面,断面间的距离为Δx。如果在t时刻i-1和i+1处的水头和流量已知,则根据下列正、负水锤特征方程联立求解,即可求得t+Δt时刻的管中i断面的水头Hi和流量Qi,即 +a时:Qi=Cp?CaHi?a时:Qi=Cn?CaHi(1) 式中,Cp=Qi-1+CaHi-1-CfQi-1|Qi-1|;Cn=Qi+1-CaHi+1-CfQi+1|Qi+1|;Ca=gAa,a为水锤波波速;Cf=fΔt2DA,f为管子摩阻系数;A、D分别为管子过流断面积和管径。 求整条管道上不同时刻的流量和水头,将管道分为若干长为Δx的管段,再求出波速a,通过距离Δx所需时间Δt=Δx/a,当在Δt开始时各断面的流量和水头已知时,即可利用正、负特征方程求出相邻两管段中间断面在Δt时刻末的流量和水头,依次类推,可求出整条管道各断面在Δt时刻末的流量和水头值,然后再根据Δt末求得的已知值求出2Δt时刻末的未知值,依次再求出3Δt时刻的未知值,直至计算到所规定的时间或者精度为止。 本输入管路中装有加压泵,在泵前后的管路断面水头和流量以及泵工作变量的确定,除利用特征方程外,尚需根据泵的特性列出附加方程才能求得,在装有泵的节点处,要求的变量有10个,即管路断面和泵出口断面上的Qpi-1,Hpi-1,Qpi,H