停泵水锤计算及其防护措施
停泵水锤是水锤现象中的一种，是指水泵机组因突然断电或其他原因而造成的开阀状态下突然停车时，在水泵及管路系统中，因流速突然变化而引起的一系列急剧的压力交替升降的水力冲击现象。一般情况下停泵水锤最为严重，其对泵房和管路的安全有极大的威胁，国有几座水泵房曾发生停泵水锤而导致泵房淹没或管路破裂的重大事故。
停泵水锤值的大小与泵房中水泵和输水管路的具体情况有关。在泵房和输水管路设计时应考虑可能发生的水锤情况，并采取相应的防措施避免水锤的发生，或将水锤的影响控制在允许围。我院在综合国外关于水锤的最新科研成果并结合多年工程实践的经验，以特征线法为基础开发了水锤计算程序。这一程序可较好地模拟各种工况条件下水泵及输水管路系统的水锤状况，为高扬程长距离输水工程提供设计依据。 1 停泵水锤的计算原理
停泵水锤的计算有多种方法:图解法、数解法和电算法。其基本原理是按照弹性水柱理论，建立水锤过程的运动方程和连续方程，这两个方程是双曲线族偏微分方程。 运动方程式为：
连续方程式为：
式中 H ——管中某点的水头 V——管流速 a——水锤波传播速度 x——管路中某点坐标 g——重力加速度 t——时间 f——管路摩阻系数 D——管径
通过简化求解得到水锤分析计算的最重要的基础方程： H-H0＝F(t-x/a)+F(t+x/a) (3)
V-V0＝g/a×F(t-x/a)-g/a×F(t+x/a) (4) 式中 F(t-x/a)——直接波 F(t+x/a)——反射波
在波动学中，直接波和反射波的传播在坐标轴(H，V)中的表现形式为射线，即特征线。它表示管路中某两点处在水锤过程中各自相应时刻的水头H与流速V之间的相互关系。为了方便计算机的计算，将上述方程组变
换为水头平衡方程和转速改变方程，即成事故停泵时水泵的两个边界条件方程式：
F1＝PM-BQv+Hn(β2+v2)(A0+A1x)-ΔH0v2/(τ2)＝0 (5) F2＝(β2+v2)(B0+B1x)+m0-C3(β0-β)＝0 (6) 式中 β——N/Nn(实际转速/额定转速) v——Q/Qn(实际流量/额定流量)
通过上述两式的联立，采用牛顿—莱福生迭代公式，可以解出v和β的近似数值。
将水泵的全面性能曲线改造为仅与转速和流速有关的全面性能曲线，以便计算机在解方程时取值，即： WH(x)＝h/(β2+v2) (7) WM(x)＝m/(β2+v2) (8) 式中 h——H/Hn(实际扬程/额定扬程) m——M/Mn(实际转矩/额定转矩) 现行的水锤计算方法就是基于上述原理。 2 几种边界条件下水锤工况的模拟结果
根据我国南方某城市取水泵房的水泵及输水管线的实际情况，采用计算机程序模拟水锤情况如下（均按最低枯水位计算）：
基本情况：
水泵机组：Qn=5 000 m3/h,Hn=55 m,Nn=741 r/min,Ns=132.4,GD2=874.7 kg．m2，Mn=932.72 kg．m，近期单台运行，远期两台运行。
输水管线：DN=1 400 mm,L=5 750 m，几何扬程：35 m（近期），45 m(远期)。
泵房和输水管线如图1所示：
2.1 假设为有阀管路停泵水锤 ① 普通止回阀
普通止回阀管路停泵水锤计算结果如表1所示。
表1 普通止回阀管路停泵水锤计算结果 运行 条件 泵出口处最大 当加大水泵机组转动惯量适当的 压力值(kPa) 倍数时，泵出口处最大压力值(kPa) 两台水泵 1 536 830(84.68 m) 并联运行 (156.7 m) 一台水 892 617(63.0 m) 泵运行 (91.05 m) 表1中所列数据为假设水泵出口处的流速为零时阀门即刻关闭所产生的水锤压力值。实际工况中，阀门关闭总要一段时间，因此实际水锤值将与表中所列数据有出入。根据计算机模拟结果，如果在此条件下适当增加水泵机组的转动惯量可以将水锤压力值明显降低。 ② 缓闭止回阀
缓闭止回阀管路停泵水锤计算结果如表2所示。
表2 缓闭止回阀管路停泵水锤计算结果 最佳阀门关 运行 条件 最佳阀门关 水泵最大倒 水泵最大 闭条件下泵 闭组合条件 转速度比 倒流量比 出口处最大 (关阀时间和 (β/βn) (Q/Qn) 压力值(kPa) 阀门关闭程度) 快关：3 s，80% 0.06 并联运行 (88.5 m) 一台水 655 0.03 泵运行 (66.8 m) -0.15 慢关：21 s，20% -0.23 慢关：24 s，20% 快关：3 s，80% 两台水泵 867 经过计算机模拟，当关阀时间和快慢组合与最佳模拟条件不同时，泵前最大压力值都将有所增加。因此一个装有两阶段关闭阀门的输水系统，其阀门的操作过程应经过计算确定，并应在试运行中调整。此种设备定货
停泵水锤的计算方法
