基于EFAST的空气罐水锤防护效果的全局敏感性分析

冉红; 蒋劲; 廖志芳; 张白云; 彭宇; Hong RAN; Jin JIANG; Zhi-fang LIAO; Bai-yun ZHANG; Yu PENG

中国农村水利水电 ›› 2021 ›› (8) : 216-220.

供水工程

基于EFAST的空气罐水锤防护效果的全局敏感性分析

作者信息 +

Global Sensitivity Analysis of Air Vessel for Water Hammer Protection Based on EFAST Method

Author information +

稿件信息 +

摘要

分析探讨长距离有压输水系统中立式空气罐的各参数变化对其水锤防护能力的影响。以某供水工程为例，采用扩展傅里叶幅度检验（EFAST）法，把事故停泵工况中的最大正压力、最低负压和水泵最大倒转转速作为空气罐水锤防护的结果，分析各参数基于基准值±50%变幅时，对空气罐水锤防护效果的全局敏感性。结果表明：空气罐影响最大水锤压力的主要参数为连接管直径、空气罐体积和气液比；影响最低负压的主要参数为连接管直径、出口局部阻力系数、气液比以及连接管长度；影响水泵最大倒转转速的主要参数依次为连接管直径、气液比和空气罐体积。其中，空气罐与主管道之间的连接管直径是影响最大压力、最低负压和水泵最大倒转转速最敏感的参数。

Abstract

The influence of various parameter changes of the vertical air vessel in the long-distance water conveyance system on its water hammer protection ability is analyzed and discussed. By taking a water supply project as an example， the Extended Fourier Amplitude Test （EFAST） method is used. The maximum positive pressure， the minimum negative pressure and the maximum reversing speed of the pump in the accidental pumping condition are taken as the results of the air vessel water hammer protection. An analysis of the parameter is based on the global sensitivity of the air vessel water hammer protection effect when the reference value is ±50%. The results show that the main parameters that affect the maximum water hammer pressure of the air vessel are the diameter of the connecting pipe， the volume of the air vessel and the gas-liquid ratio； the main parameters that affect the lowest negative pressure are the diameter of the connecting pipe， the outlet local resistance coefficient， the gas-liquid ratio and the length of the connecting pipe. The main parameters that affect the maximum reversing speed of the pump are the diameter of the connecting pipe， the gas-liquid ratio and the volume of the air vessel in order. Among them， the diameter of the connecting pipe between the air vessel and the main pipe is the most sensitive parameter that affects the maximum pressure， the minimum negative pressure and the maximum reverse speed of the pump.

导出引用

冉红 , 蒋劲 , 廖志芳 , 张白云 , 彭宇. 基于EFAST的空气罐水锤防护效果的全局敏感性分析[J].中国农村水利水电, 2021(8): 216-220

Hong RAN , Jin JIANG , Zhi-fang LIAO , Bai-yun ZHANG , Yu PENG. Global Sensitivity Analysis of Air Vessel for Water Hammer Protection Based on EFAST Method[J].China Rural Water and Hydropower, 2021(8): 216-220

长距离有压输水工程通常具有输水流量大、地形起伏多变、扬程高和管路系统复杂等特点。当管路系统局部水力边界条件发生变化，如启泵或停泵、开阀或关阀和进出口水位变化等，都会诱导管路发生剧烈的压力波动，产生巨大的水锤压力，严重影响了管路系统的正常运行^［1］。通常，水锤防护措施主要有：单向调压塔、双向调压塔、空气阀、空气罐、超压泄压阀及联合防护等^［2，3］。空气罐相对于其他防护措施具有良好的水锤调控能力、施工方便和安全稳定性更高等特点。目前，工程上经常应用的空气罐结构类型主要有水气自然分离式、隔膜式和囊式3种^［4］。影响空气罐水锤防护效果的主要参数有空气罐体积、气液比、阻抗孔尺寸、安装方式和连接管参数等。国内外学者对这些参数的优化做了大量研究。邓安利等^［5］在保证空气罐体积相同和罐内水体不放空条件下，得出罐内初始气体体积越大，水锤防护效果更好。龚娟等^［6］分析了阻抗孔尺寸大小对空气罐水锤防护效果的影响，在保证管道沿线不出现负压时，尽量选取较小的阻抗孔能有效降低管道压力波动。魏振荣^［4］等分析了空气罐和主管道之间的连接管直径；在一定范围内直径越大，空气罐水锤防护效果越好，且存在一个最优的连接管直径。孙强^［7］等对比分析了空气罐连接管带有旁通管、带有节流阀的旁通管和带有独立流入与流出短管的三种类型空气罐的水锤防护效果，并利用序列二次规划法对空气罐参数进行了优化。此外，大多数都是空气罐与空气阀、单向调压塔、水泵出口阀门等联合防护水锤的研究。然而，由于空气罐模型参数较多、参数间多为非线性关系且实际工程复杂多变，给模型参数的拟定带来了一定的难度。合适的空气罐尺寸不仅能够达到预期的水锤防护效果，还能够降低工程造价节约成本。因此，开展对空气罐模型参数的敏感性研究，确定各个参数对水锤防护效果的影响程度，进而对影响程度较大的参数进行重点拟定，忽略次要因素，提高空气罐的设计效率。

敏感性分析方法一般分为局部敏感性分析和全局敏感性分析。局部敏感性分析为考虑模型单个输入变量对模型输出结果的直接影响，未考虑参数间互相影响；全局敏感性分析则同时考虑多个输入变量及变量间的相互作用对模型输出结果的影响，多用于输入变量较多的非线性模型。当前，常用的全局敏感性分析方法主要有Morris参数筛选法、Sobo法和傅立叶幅度敏感性分析法（FAST，及EFAST）等^［8］。EFAST法是SAILTELI^［9］等结合Sobol法和FAST法提出的一种新的全局敏感性分析方法，通常具有计算量大、稳健和获取信息多等特点。目前，大多数研究集中在空气罐参数对水锤防护影响的局部敏感性分析，对空气罐参数的全局敏感性分析尚未见报道。因此，本次研究采用计算量大、可靠性更高的EFAST法，以某工程为例，分析空气罐的体积、气液比，空气罐与主管道之间的连接管直径、长度、进出局部阻力系数对水锤防护效果的影响，为空气罐参数的选取和设计提供科学依据。

1 EFAST方法

水锤计算一般采用一维非恒定流连续性方程和动量方程联立求解。在有压输水管道系统中，通常借助于特征线法，将水锤基本微分方程离散为有限差分方程，进而求解微分方程的数值解。水锤基本微分方程为：

\frac{\partial V}{\partial t} + V \frac{\partial V}{\partial x} + g \frac{\partial H}{\partial x} + \frac{f' V |V|}{2 D} = 0

（1）

\frac{\partial H}{\partial t} + V \frac{\partial H}{\partial x} - V s i n θ + \frac{c^{2}}{g} \frac{\partial V}{\partial x} = 0

（2）

式中：H为管道测压管水头值m；V为管道横断面平均流速，m/s；t为时间，s；g为重力加速度，m/s²；

x

为管道轴线上的坐标；

θ

为管轴线与水平线的夹角，°；D为管道直径，m；c为水锤波速，m/s；

f'

为沿程阻力系数。

空气罐是内部充满高压气体的金属密封容器，其上层为气体，下层为液体；通常为了方便安装、运行维护和安全管理，将其安装在泵房内。当水泵机组事故停机时，管道中压力降低，空气罐内上层气体迅速膨胀，下层液体在气体压力作用下，通过空气罐与主管道之间的连接管迅速补给到主管道，防止管道出现严重负压；反之，防止管道出现过大正压。如图1所示，假设空气罐内气体符合等温绝热的气体状态方程，忽略气体与罐壁摩擦力。根据液体的连续性方程、水头平衡方程及管道水锤相容性方程，确定空气罐的边界条件如下。

图1 空气罐示意图

Fig.1 Schematic diagram of air vessel

Full size|PPT slide

流量连续性方程为：

Q_{p 1} = Q_{s t} + Q_{p 2}

（3）

水头平衡方程为：

H_{P} = H_{s t} + \frac{(P - P_{0})}{γ} + f Q_{s t} |Q_{s t}|

（4）

罐内水位与流量方程：

A_{s t} \frac{d H_{s t}}{d t} = Q_{s t}

（5）

等温绝热气体状态方程：

p V^{n} = C

（6）

式中：

Q_{p 1}

，

Q_{p 2}

分别为空气罐上、下游流量，m³/s；

Q_{s t}

为流入空气罐的流量，m³/s；

H_{P}

为空气罐与主管连接处压力，m；

H_{s t}

为空气罐内液位高度，m；

p

为罐内气体绝对压力，Pa；

p_{0}

为空气罐外部大气压，Pa；

γ

为水的容重，N/m³；

A_{s t}

为空气罐截面积，m²；

f

为空气罐与主管道之间的水力损失系数；

V

为罐内气体体积，m³；n为等温绝热指数，取值为1~1.4，一般取为1.2。

全局敏感性分析主要考虑了模型输入变量的变化方差以及变量之间的交互作用对模型输出结果的总体影响，能够定性和定量的分析模型输入变量对输出结果的贡献率^［10］。EFASTF法是一种基于模型方差分解的全局敏感性分析方法，该方法整合了FAST方法和Sobol方法的优点，通过对模型输出结果进行傅里叶分析，得到参数及参数之间相互作用所引起的模型结果方差；其一阶敏感性指数反映了单一参数对模型输出结果的贡献率；总敏感性指数则反映了参数直接和间接对模型输出结果的贡献率；二者的差值可以表示参数之间相互作用的程度，差值越大，表明参数之间的相互影响越大^［11］。模型

Y = f (x_{1}, x_{2}, \dots, x_{n})

的总方差可以分解：

V (Y) = \sum_{i = 1}^{n} V_{i} + \sum_{i < j \leq n}^{n} V_{i j} + \dots + \sum_{i < j \leq n}^{n} V_{(1,2, \dots, n)}

（9）

S_{i} = \frac{V_{i}}{V}

（10）

S_{T i} = \frac{V - V_{- i}}{V}

（11）

式中：

V (Y)

为模型总方差；

V_{i}

为

x_{i}

自身变化引起的模型方差；

V_{i j}

为

x_{i}

通过

x_{j}

作用贡献的方差；

V_{(1,2, \dots, n)}

为

x_{i}

通过余下n-1个参数相互作用贡献的方差；

S_{i}

为

x_{i}

的一阶敏感性指数；

V_{- i}

为不包括参数

x_{i}

的所有参数方差之和；

S_{T i}

为

x_{i}

的总敏感性指数。

2 基本资料

2.1 工程基本概况

参数的敏感性分析通常需要设置模拟情景。本次研究以某供水工程为例，输水管线全长2 593 m，干管为DN1300的球墨铸铁管，采用离心泵进行有压输水。泵站为一管三机，三用一备，额定流量为1.415 m³/s，额定扬程为610 m，额定转速为 1 000 r/min，配套电机功率为9 962 kW，机组转动惯量为1 800 kg·m²。

图2 输水管线纵断面图

Fig.2 The elevation profile of the pipeline

Full size|PPT slide

2.2 空气罐参数选取

本文的水锤防护计算选择应用非常广泛的流体管网分析软件PIPENET。应用PIPENET建模，对该工程的水锤防护进行仿真计算。通常空气罐和主管之间由连接管、阀门和一些附属设施组成。局部阻力系数主要考虑连接管的弯头、阀门及进出口断面突然扩散和收缩，取值参考《水力计算手册》^［12］和前人研究^［7］。

经过多次试算，确定防护措施为：泵后阀门采用液控缓闭蝶阀，发生事故停机时，阀门3.5 s快关75%开度，26.5 s慢关，全管线只用空气罐防，且安装在水泵出口。当空气罐总体积为39.25 m³，其中液体体积为总体积的50%，与主管道之间的连接管直径为DN400的钢管，长度为10 m，进口阻力系数为5，出口阻力系数为3.5时，水锤防护计算结果为：水锤最大压力765.90 m，出现在阀门出口；最小压力为-1.41 m，出现在桩号2+596处；机组最大倒转速为-1 176.5 r/min。以上计算值均满足《泵站设计规范》（GB 50265-2010）^［13］中关于水锤防护最大压力和水泵倒转等相关要求，且管道全线未出现水柱断裂。

图3 管道压力包络线

Fig.3 Pipeline pressure envelope

Full size|PPT slide

用以上计算采用的空气罐体积、气液比以及与主管道之间的连接管参数，作为空气罐参数敏感性分析的基准值。由于基于空气罐模型参数的敏感性分析较少，结合其他领域相关敏感性研究^{［11，14］}，将空气罐各参数增减变化幅调整为50%，参数具体取值范围如表1，此时各参数变化值在合理区间范围内，且模拟结果差异性明显，可进一步做空气罐的敏感性分析。

表1 空气罐模型参数取值

Tab.1 Air Vessel model parameter value

参数名	参数基准值	取值区间
空气罐体积V/m³	39.25	19.6~58.88
液体体积占总体积比率	0.5	0.25~0.75
连接管长度/m	10	5~15
连接管进口阻力系数	5	2.5~7.5
出口阻力系数	3.5	1.75~5.25
连接管管直径/mm	400	200~600

2.3 敏感性分析方案

将表1中各参数作为水锤计算的输入数据，保持阀门关闭规律和其他参数不变。应用非常广泛的流体管网分析软件PIPENET进行水锤数值仿真计算，以最大压力、最小压力和水泵机组最大倒转速作为空气罐水锤防护效果的输出结果。应用EFAST法分析空气罐各参数对其水锤防护效果的影响。

全局敏感性分析借助欧盟委员会联和研究中心提供的 simlab（Version 2.2.1）的EFAST模块完成。此软件是基于Monte Carlo方法来研究模型参数的交互作用及敏感性。具体步骤如下：

（1）选取合适的空气罐各参数作为基准值，并以基准值上下变幅50%，确定参数取值范围，同时保持其他参数不变。

（2）在Simlab中输入参数范围，且参数服从均匀分布。

（3）利用Monte Carlo方法对参数进行随机采样。EFAST法需要每个参数采样次数至少大于65次，才能进行敏感性分析。此次，每个参数的采样次数505次，共采样3 030 组。

（4）将采样的参数依次对应输入空气罐模型，调用PIPENET进行仿真计算，读取计算结果的最大压力、最小压力及水泵的最大倒转速。

（5）将计算结果整理为Simlab可识别的文件格式，将其作为外部文件，输入Simlab对参数进行敏感性分析。

3 结果分析与讨论

3.1 敏感性分析结果

空气罐体积、气液比以及与主管道之间的连接管直径、长度和局部阻力系数对空气罐水锤防护影响的一阶敏感性指数如表2所示和全局敏感性指数如表3所示。

表2 空气罐参数的一阶敏感性指数

Tab.2 The first-order sensitivity index of air Vessel parameters

参数	输水系统最大压力	输水系统最小压力	水泵的最大倒转速
空气罐体积	0.183 4	0.018 9	0.039 8
气液比	0.171 1	0.027 9	0.034 1
连接管直径	0.624 8	0.610 6	0.581 6
连接管长度	0.000 8	0.014 6	0.000 4
出口局部阻力系数	0.001 4	0.039 0	0.013 1
进口局部阻力系数	0.013 2	0.000 4	0.007 5

表3 空气罐参数的全局敏感性指数

Tab.3 Global sensitivity index of air Vessel parameters

参数	输水系统最大压力	输水系统最小压力	水泵的最大倒转速
空气罐体积	0.208 2	0.139 5	0.298 8
气液比	0.196 9	0.149 6	0.315 7
连接管直径	0.668 1	0.885 3	0.702 4
连接管长度	0.002 7	0.122 3	0.024 1
出口局部阻力系数	0.007 5	0.229 7	0.007 8
进口局部阻力系数	0.025 2	0.077 5	0.013 8

目前，对敏感性分析划定尚无统一的标准，据DEJONG^［15］在EFAST法中界定的一阶敏感指数S_i >0.05，总敏感指数ST_i >0.1为敏感性指数的取值标准，当敏感性指数小于界定的值时，不予以讨论。

3.2 参数敏感性分析

根据图4所示，在一定范围内，连接管直径、空气罐体积和气液比的一阶敏感性指数超过了0.05，全局敏感性指数超过了0.1，是影响空气罐防护正压水锤的主要参数，且都是直接效应占主导；其中，连接管直径是最敏感的参数，可解释其变化方差的66.81%。在长距离输水工程中，应用空气罐防护正压水锤时，应首先考虑空气罐连接管直径、体积、气液比；且3个参数之间的相互作用不强，通常可以单独对3个参数做优化，降低空气罐造价。

图4 影响最大压力的参数敏感性分析

Fig.4 Sensitivity analysis of parameters affecting the maximum pressure

Full size|PPT slide

根据图5所示，在一定范围内，只有连接管直径的一阶敏感性指数超过了0.05，但连接管直径、出口局部阻力系数、气液比、空气罐体积和连接管长度的全局敏感性指数均超过了0.1。空气罐参数对于最低负压的影响顺序依次为：连接管直径、出口局部阻力系数、气液比、空气罐体积、连接管长度。但其主要的参数是连接管直径，可解释其变化方差的88.53%；其直接影响和与其他参数的相互作用都是最强的。在长距离输水工程中应用空气罐防护负压时，可以尽量通过减小出口局部阻力系数、连接管长度，在合理范围内选取较大的连接管直径，或适当增加的空气罐体积和气液比，来达到需要的防护标准。

图5 影响最小压力的参数敏感性分析

Fig.5 Sensitivity analysis of parameters affecting the maximum pressure

Full size|PPT slide

根据图6所示，在一定范围内，只有连接管直径的一阶敏感性指数超过了0.05，但连接管直径、空气罐体积、气液比的全局敏感性指数均超过了0.1。水泵最大倒转速对连接管直径最敏感，可解释其变化方差的70.24%。而空气罐体积和气液比主要是通过相互作用影响水泵的倒转。在实际工程应用中，在有效范围内，空气罐体积较小，气液比越高，水泵最大倒转越小。

图6 影响水泵倒转的参数敏感性分析

Fig.6 Sensitivity analysis of parameters affecting pump inversion

Full size|PPT slide

4 结论

本次研究以供水工为例，全管线水锤主防护为空气罐，借助可靠性较高的EFAST方法对空气罐的体积、气液比，与主管道之间连接管长度、直径和进出口局部阻力系数共6个参数进行了敏感性分析。得到的主要结论如下：

（1）对于系统最大压力，空气罐与主管道之间的连接管直径是最敏感的参数，其次是空气罐体积和气液比，且都是直接效应占主主导，而连接管长度和进出口阻力系数影响对其影响相对较小。

（2）对于系统最低负压，空气罐与主管道之间的连接管直径是最敏感的参数。然而，出口局部阻力系数、气液比、空气罐体积和连接管长度的参数交互作用不能忽视，其中出口局部阻力系数的交互作用最强。进口局部阻力系数对最低负压的影响不大。

（3）对于水泵最大倒转，空气罐与主管道之间的连接管直径是最敏感的参数，其次是空气罐体积和气液比，连接管长度和进出口局部阻力系数影响不大。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	WYLIE E B， STREETER V L， SUO L S. Fluid transient in systems［M］.New York：McGraw-Hill Press，1993. 本文引用 [1]

2	李琨，吴建华，刘亚明，等.空气罐对泵站水锤的防护效果研究［J］.人民长江，2020，51（2）：200-204. 本文引用 [1]

3	马世波，张健．长距离输水工程停泵水锤防护措施研究［J］．人民长江，2009，40（1）：85-86. 本文引用 [1]

4	魏振荣.空气罐的水锤防护特性及其工程应用［J］.东北水利水电，2017，35（6）：67-70. 本文引用 [2]

5	邓安利，蒋劲，兰刚，等.长距离输水工程停泵水锤的空气罐防护特性［J］.武汉大学学报（工学版），2015，48（3）：402-406. 本文引用 [1]

6	龚娟，张健，俞晓东.高扬程输水系统空气罐阻抗孔尺寸优化［J］.水电能源科学，2013，31（5）：166-169. 本文引用 [1]

7	SUN Q， WU Y B， XU Y，et al. Optimal sizing of an air vessel in a long-distance water-supply pumping system using the SQP method［J］.Journal of Pipeline Systems Engineering and Practice， 2016，7（3）：161－174. 本文引用 [2]

8	任启伟，陈洋波，舒晓娟.基于Extend FAST方法的新安江模型参数全局敏感性分析［J］.中山大学学报（自然科学版），2010，49（3）：127-134. 本文引用 [1]

9	SALTELLI A， STAR ANTOLA， CHAN P S. A quantitative model-in-dependent method for global sensitivity analysis of model output［J］. Technometrics， 1999，41（1）：39-56. 本文引用 [1]

10	SALTELLI A， CHAN K， SCOTT E M.Sensitivity Analysis［M］.New York：Wiley，2000：58-66．本文引用 [1]

11	崔金涛，邵光成，林洁，等.基于EFAST的CROPGRO-Tomato模型参数全局敏感性分析［J］.农业机械学报，2020，51（1）：237-244. 本文引用 [2]

12	李炜.水力计算手册［M］.北京：中国水利水电出版社，2007：9-12. 本文引用 [1]

13	中国华人民共和国水利部. 泵站设计规范：GB50265-2010 ［S］.北京：中国计划出版社，2011：49. 本文引用 [1]

14	霍世璐，段磊，刘明明，等.VG模型参数对灌溉补给系数的敏感性分析［J］.水利水电技术， 2020，51（7）：177-185. 本文引用 [1]

15	兴安，卓志清，赵云泽，等.基于EFAST的不同生产水平下 WOFOST模型参数敏感性分析［J］.农业机械学报，2020，51（2）：161-171. 本文引用 [1]

PDF(1379 KB)

访问

引用

详细情况

段落导航

收稿日期	出版日期
2020-11-09	2021-08-15
发布日期
2021-09-24

选择文件类型/文献管理软件名称

选择包含的内容