摘 要：針對水錘引起的管道風(fēng)險問題，考慮閥門初始開度、閥門關(guān)閉歷時的不確定性，推導(dǎo)了簡單管道極限水錘壓力概率分布函數(shù)。相對極限水錘出現(xiàn)在0.8～1.2之間為大概率事件，宜選取確定性分析得到的極限水錘壓力的1.2倍進行結(jié)構(gòu)設(shè)計。運用蒙特卡羅法計算得到兩種設(shè)計風(fēng)險值分別為12 285.00、588.44，對比說明了有必要基于水錘的不確定性對管道進行風(fēng)險評估和設(shè)計校核。

關(guān)鍵詞：閥門初始開度;不確定性;風(fēng)險;蒙特卡羅法;極限水錘

中圖分類號：TV732.4+1 ? 文獻標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2020.04.024

Abstract：In order to solve the risk in a simple pipeline caused by water hammer， firstly， the uncertainty of valve shutdown and the initial valve opening were taken into account. The paper deducted the probability distribution function of limit water hammer of simple pipe. It shows that the value of relative limit water hammer pressure is great probability in 0.8-1.2. It should take 1.2 times of the deterministic analysis results as a guideline of design. Then using the Monte Carlo method to calculate the two kinds of design risk values of 12 285.00 and 588.44 respectively. Contrasting the two situations， it illustrates the risk assessment and design checking based on water hammer uncertainty is necessary.

Key words： water hammer; uncertainty; risk; method of Monte Carlo; limit water hammer

影響壓力管道結(jié)構(gòu)安全的因素眾多，如水文氣象條件、地形地質(zhì)狀況、荷載大小及材料性能等，其中水錘現(xiàn)象是導(dǎo)致壓力管道失效的主要原因[1]。閥門快速關(guān)閉時管路中液體流速急劇改變，造成壓力顯著、反復(fù)、迅速變化的現(xiàn)象叫水錘現(xiàn)象[2]。水錘現(xiàn)象發(fā)生時，壓力的反復(fù)變化和瞬時升高會造成管道損壞，嚴(yán)重時其升高值將超過管壁的極限承載力，發(fā)生管道炸裂等事故[1-3]。水錘壓力與閥門關(guān)閉歷時、閥門初始開度、初始流量有關(guān)，而這些變量都是隨機變量，因此水錘壓力也是隨機變量。張靜[4]對比了多種參數(shù)的不確定性，指出水錘是引起管道風(fēng)險的主要不確定性變量，其他變量可按確定性變量考慮。郭文鑄等[5]采用隨機分析方法考慮電站負荷、水庫水位和關(guān)閉歷時的不確定性，得到了簡單管系統(tǒng)極限水錘壓力服從對數(shù)正態(tài)分布。張芹芬[6]考慮水庫水位、導(dǎo)葉初始開度，推導(dǎo)了簡單管系統(tǒng)極限水擊壓力分布。水錘壓力隨機分析已經(jīng)取得一定成果，但多數(shù)研究并沒有考慮閥門初始開度和關(guān)閉歷時不確定性隨機變量的影響。筆者考慮閥門關(guān)閉歷時、初始開度的不確定性，通過理論分析推導(dǎo)了庫-管道-閥門簡單系統(tǒng)極限水錘壓力的概率分布，采用舍取法產(chǎn)生服從極限水錘壓力的隨機數(shù)，用蒙特卡羅法對管道失效概率進行了分析。

1 極限水錘統(tǒng)計特性

1.1 閥門初始開度統(tǒng)計特征

由于閥門控制組具有復(fù)雜性、人為操作的不確定性，或者閥門經(jīng)過長時間運行材料老化，出現(xiàn)卡死情況，無法保證每次閥門初始開度，因此初始開度具有不確定性。通過統(tǒng)計92個閥門初始開度的實測資料，繪制了頻數(shù)直方圖，見圖1。

一般情況下根據(jù)頻率或頻數(shù)直方圖可假設(shè)可能樣本分布。Dn為離散直方圖頻數(shù)Fn（xk）（xk為開度值）與其擬合連續(xù)函數(shù)相對應(yīng)開度的頻數(shù)值之差的最大絕對值，即Dn=max|Fn（xk）-F（xk）|，比較Dn與Dn，1-α的大小，Dn，0.90≈1.23/n0.5，Dn，0.95≈1.36/n0.5，Dn，0.99≈1.63/n0.5（n為樣本數(shù)量）。檢驗樣本是否服從對數(shù)正態(tài)分布：當(dāng)Dn

水庫斷面面積遠遠大于壓力管道斷面面積，且水錘現(xiàn)象持續(xù)時間短，故水庫水位變化可忽略不計，H0為常數(shù)。

2.2 產(chǎn)生隨機數(shù)

若計算所得概率密度函數(shù)不是基本的分布，則用公式直接產(chǎn)生服從該概率密度函數(shù)的隨機數(shù)具有一定的困難。采用舍取法可以產(chǎn)生服從任意分布的隨機數(shù)，舍取法[11-14]產(chǎn)生隨機數(shù)的基本思想為：產(chǎn)生2個相應(yīng)范圍內(nèi)均勻分布隨機數(shù)xi、yi。xi、yi正交并形成笛卡兒平面坐標(biāo)系，已知概率密度函數(shù)f（x），比較yi與f（xi）的大小，若yi >f（xi），則隨機數(shù)xi不服從此分布，丟棄xi;若yi

3 算 例

3.1 確定性計算

閥門初始開度τ0=1，管道長度L=100 m，水庫水位H0=40 m，閥門有效關(guān)閉歷時Tc0=10 s，管道初始流速V0=4 m/s，波速a=1 200 m/s。ρτ0>1，最大水錘壓力為極限水錘壓力。確定性計算得極限水錘壓力Hm0=86.16 m。

3.2 極限水錘壓力概率密度

考慮龍貝格求解Hm概率得相對極限水錘壓力β=Hm/Hm0（Hm0為確定性計算極限水錘壓力;Hm為水錘壓力）。相對極限水錘壓力概率密度函數(shù)如圖2所示。

采用不確定性分析，相對確定性分析失效概率大大降低，僅為確定性分析的1/30。計算1 000次到5 000次相對誤差后者大概為前者的12倍，計算到500 000次相對誤差均很小，說明計算的失效概率較大時，隨計算次數(shù)增多，誤差收斂速度快，計算較小失效概率不具備此特點，但當(dāng)計算次數(shù)足夠多時，計算結(jié)果精度能滿足要求。蒙特卡羅法計算的管道失效概率為0.048 1，經(jīng)濟總量為255 400元，風(fēng)險值P=失效概率×經(jīng)濟總量=12 285.00?？紤]相對極限水錘壓力的1.2倍進行結(jié)構(gòu)設(shè)計，失效概率為0.001 632，經(jīng)濟總量為367 776元，風(fēng)險值為588.44，遠遠小于確定性的風(fēng)險值。當(dāng)計算次數(shù)為1 000 000時，P和Pf計算時間分別為7.22 s和9.01 s，證明在現(xiàn)代技術(shù)條件下蒙特卡羅法能夠得到很好地應(yīng)用，且具有較高精度。

4 結(jié) 論

（1）考慮閥門初始開度和關(guān)閉歷時的不確定性，通過理論推導(dǎo)，得到了極限水錘壓力概率密度函數(shù)，計算得到其極限水錘壓力為確定性分析結(jié)果的0.6～1.4倍，故有必要在實際工程中考慮閥門開度和閥門關(guān)閉歷時的不確定性，但其值出現(xiàn)在確定性計算結(jié)果的0.8～1.2倍之間為大概率事件。為了確保供水系統(tǒng)正常運行，宜選取0.2倍極限水錘壓力設(shè)計值作為安全裕度。

（2）確定性分析和不確定性分析得到的風(fēng)險值分別為12 285.00和588.44，對比說明了考慮水錘不確定性進行風(fēng)險分析的必要性。計算次數(shù)為1 000 000時，兩者計算時間分別為7.22 s和9.01 s，且計算結(jié)果相對誤差可忽略，說明蒙特卡羅法具有較高的精度且在現(xiàn)代計算機技術(shù)條件下具有可行性。

參考文獻：

[1] 劉梅清，梁興，劉志勇.長管道事故停泵水錘現(xiàn)場測試與信號分析[J].排灌機械工程學(xué)報.2012，30（3）：249-253.

[2] 金錐.停泵水錘及其防護[M].北京：中國建筑工業(yè)出版社，2003：79-113.

[3] 劉華，鞠小明，陳家遠.供水管道系統(tǒng)中的水力過渡過程研究[J].四川大學(xué)學(xué)報，1999，10（1）：56-59.

[4] 張靜.水電站壓力鋼管的風(fēng)險分析與計算[D].烏魯木齊：新疆農(nóng)業(yè)大學(xué)，2015：32-34.

[5] 郭文鑄，索麗生.簡單管水擊隨機分析[J].水力發(fā)電學(xué)報，1996，15（4）：72-81.

[6] 張芹芬.水電站水擊與調(diào)壓室涌浪的隨機分析[D].南京：河海大學(xué)，1998：78-81.

[7] 楊惠蓮，張濤.誤差理論與數(shù)據(jù)處理[M].天津：天津大學(xué)出版社，1992：56-65.

[8] 王長新.施工導(dǎo)流風(fēng)險分析及計算[D].南京：河海大學(xué)，1995：23-26.

[9] 趙國藩.工程結(jié)構(gòu)可靠性理論與應(yīng)用[M].大連：大連理工大學(xué)出版社，1996：76-88.

[10] 中華人民共和國建設(shè)部，國家質(zhì)量監(jiān)督檢驗檢疫總局.給水排水工程管道結(jié)構(gòu)設(shè)計規(guī)范：GB 50332—2002[S].北京：中國建筑工業(yè)出版社，2002：17-19.

[11] 楊綠峰，張偉，韓曉鳳.水電站壓力鋼管整體安全評估方法研究[J].水力發(fā)電學(xué)報，2011，30（5）：149-156.

[12] 袁子厚，陳明祥，何英明.水電站中明管可靠度校準(zhǔn)分析[J].重慶建筑大學(xué)學(xué)報，2007，3（5）：75-77.

[13] 劉敬敏，楊綠峰，張偉.等基于系統(tǒng)可靠度的明鋼管整體安全性研究[J].水利水運工程學(xué)報，2013，6（4）：74-80.

[14] 王菊，王和明，徐海龍，等.任意分布隨機數(shù)的FPGA實現(xiàn)[J].火力與指揮控制，2015，40（4）：173-175.

[15] 郭強.泵站出水管風(fēng)險分析[D].烏魯木齊：新疆農(nóng)業(yè)大學(xué)，2016：23-25.

【責(zé)任編輯 張華巖】