关键词：停泵水锤 计算 防护措施

停泵水锤是水锤现象中的一种，是指水泵机组因突然断电或其他原因而造成的开阀状态下突然停车时，在水泵及管路系统中，因流速突然变化而引起的一系列急剧的压力交替升降的水力冲击现象。一般情况下停泵水锤最为严重，其对泵房和管路的安全有极大的威胁，国内有几座水泵房曾发生停泵水锤而导致泵房淹没或管路破裂的重大事故。

停泵水锤值的大小与泵房中水泵和输水管路的具体情况有关。在泵房和输水管路设计时应考虑可能发生的水锤情况，并采取相应的防范措施避免水锤的发生，或将水锤的影响控制在允许范围内。我院在综合国内外关于水锤的最新科研成果并结合多年工程实践的经验，以特征线法为基础开发了水锤计算程序。这一程序可较好地模拟各种工况条件下水泵及输水管路系统的水锤状况，为高扬程长距离输水工程提供设计依据。

1、停泵水锤的计算原理

停泵水锤的计算有多种方法:图解法、数解法和电算法。其基本原理是按照弹性水柱理论，建立水锤过程的运动方程和连续方程，这两个方程是双曲线族偏微分方程。

运动方程式为：

连续方程式为：

式中：H ——管中某点的水头

V——管内流速

a——水锤波传播速度

x——管路中某点坐标

g——重力加速度

t——时间

f——管路摩阻系数

D——管径

通过简化求解得到水锤分析计算的最重要的基础方程：

H－H₀＝F(t－x／a)＋F(t＋x／a) (3)

V－V₀＝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之间的相互关系。为了方便计算机的计算，将上述方程组变换为水头平衡方程和转速改变方程，即成事故停泵时水泵的两个边界条件方程式：

F₁＝PM-BQv＋H_n(β²＋v²)(A₀＋A₁x)－ΔH₀v²／(τ²)＝0　　　(5)

F₂＝(β²＋v²)(B₀＋B₁x)＋m₀－C₃(β₀-β)＝0　　　　(6)

式中　β——N／N_n(实际转速／额定转速)

v——Q／Q_n(实际流量／额定流量)

通过上述两式的联立，采用牛顿——莱福生迭代公式，可以解出v和β的近似数值。

将水泵的全面性能曲线改造为仅与转速和流速有关的全面性能曲线，以便计算机在解方程时取值，即：

WH(x)＝h／(β²＋v²) (7)

WM(x)＝m／(β²＋v²) (8)

式中　h——H／H_n(实际扬程／额定扬程)

m——M／M_n(实际转矩／额定转矩)

现行的水锤计算方法就是基于上述原理。

2、几种边界条件下水锤工况的模拟结果

根据我国南方某城市取水泵房的水泵及输水管线的实际情况，采用计算机程序模拟水锤情况如下（均按最低枯水位计算）：

基本情况：

水泵机组：Q_n＝5 000 m³／h，H_n＝55 m，N_n＝741 r／min，N_s＝132.4，GD²＝874.7 kg·m²，M_n＝932.72 kg·m，近期单台运行，远期两台运行。

输水管线：DN＝1 400 mm，L＝5 750 m，几何扬程：35 m（近期），45 m(远期)。

泵房和输水管线如图1所示：

2.1　假设为有阀管路停泵水锤

①普通止回阀

普通止回阀管路停泵水锤计算结果如表1所示。

表1　普通止回阀管路停泵水锤计算结果

运行条件	泵出口处最大压力值(kPa)	当加大水泵机组转动惯量适当的 倍数时，泵出口处最大压力值(kPa)
两台水泵并联运行	1536(156.7 m)	830(84.68 m)
一台水泵运行	892(91.05 m)	617(63.0 m)

表1中所列数据为假设水泵出口处的流速为零时阀门即刻关闭所产生的水锤压力值。实际工况中，阀门关闭总要一段时间，因此实际水锤值将与表中所列数据有出入。根据计算机模拟结果，如果在此条件下适当增加水泵机组的转动惯量可以将水锤压力值明显降低。