長輸油氣管道河溝段水毀災害氣象風險評價研究

薛俊卓，詹 輝，張國良，何 開，鄧清祿

(1.中國地質大學(武漢) 工程學院，湖北 武漢 430074;2.江西省天然氣管道有限公司，江西 南昌 330096)

對于淺埋鋼質輸油氣管道，自然災害是主要風險源之一[1]。國內外相關統計分析表明地質災害對于管道的危害性不可忽視[1-3]。

在各類地質災害中，滑坡和崩塌對管道的危害最大，但數量相對少；而數量多、發育廣泛的地質災害是水毀[3-6]，尤其是穿越河道或沖溝(統稱河溝)段的長輸油氣管道易受暴雨洪水沖刷的危害，造成河溝段水毀災害的危害性最大。河溝段水毀災害的危害性主要表現在，暴雨時，引發的洪水沖刷侵蝕河床與岸坡，可能造成順河溝或穿河溝敷設的管道出現露管或漂管，并伴隨帶來管道損傷、彎折等后果[7-8]。管道穿越中小型河流，以開挖直埋的敷設方式為主，其埋深較淺(一般小于2 m)。中小型河流實測水文資料通常不足[9]，缺乏專門的勘察，通常未設置針對性的防治暴雨洪水沖刷措施，因此開展暴雨狀況下管道河溝段水毀災害氣象風險評價具有重要意義。

對管道地質災害風險評價方面的研究，既有整體風險管控方法的研究[10]，也有專門針對某類型地質災害風險評價的研究，如對滑坡、崩塌、泥石流等高危害性地質災害風險評價方面的研究成果較多[3,11]，且模糊邏輯法[12]、耦合協調[13]等復雜數學方法也有所應用。而對埋地管道水毀災害方面的研究則相對較少，主要有以下方面：①河溝段水毀災害中埋地管道的力學特性[14]；②水毀災害風險的現場調查分析[6-7]；③應用知識或數據驅動模型來評價水毀災害的風險，如層次分析法[15]、灰色關聯可拓理論[16]、GIS因子疊加[5]等；④水毀災害發生后埋地管道水工保護措施的重新設計[8]、基于投入—產出的效率評價[17]等災后反思研究；⑤考慮多因素的半定量風險概率指數評價[18]等。雖然這些研究取得了一定的成果，但在埋地管道河溝段水毀水文氣象風險評價方面未見專門的研究。因此，本文探討了適合長輸油氣埋地管道河溝段水毀災害氣象風險評價的暴雨洪水和河溝道沖刷深度等關鍵指標的計算方法，在此基礎上設計了風險半定量評價指標體系，并對江西省天然氣管網一期工程進行了實例應用與分析，可為長輸油氣管道河溝段水毀災害氣象風險評價提供技術路線，對管道水毀災害進行預測評價。

1 河溝道沖刷深度計算方法

山地丘陵地區河溝道的沖淤變化主要是通過降雨匯流作用產生的，具有季節性和突發性。估算暴雨洪水狀態下河溝道的沖刷深度，是河溝段管道風險評價的關鍵指標。河溝道沖刷深度與洪峰時水流速、水深、河床物質等有關，其計算方法研究是河溝段管道水毀災害氣象風險評價的重要研究內容。

1. 1 暴雨洪水計算

在我國大部分地區，洪水主要由暴雨形成，本文采用根據降水量匯流計算暴雨洪水的方法(即水文氣象法)[9]，該計算過程需要綜合考慮產流計算、匯流計算以及洪峰流量與水深、流速的計算等幾個部分。

1.1.1 產流計算

一次降雨過程可分成若干間隔相等的時段Δt，則降雨時程可計為1×Δt、2×Δt、…、t×Δt，如3 h、6 h、24 h。某時程t×Δt內降雨的面雨量P面t×Δt可表示為

P面t×Δt=P點·t×Δt·α

(1)

式中：P面t×Δt為t×Δt時程內降雨的面雨量(mm)；P點t×Δt為t×Δt時程內降雨的點雨量(mm);α為點面折算系數，與暴雨時程、流域面積有關。

某時程t×Δt內的凈面雨量ht×Δt，等于毛雨量扣除下滲損失量。結合江西省水文局的暴雨統計資料，采用蓄滿產流原理的降雨徑流分析模型，可表示如下：

ht×Δt=(R總t×Δt)-min(R總t×Δt，3×fc)

(2)

R總t×Δt=f(P面t×Δt+Pa)

(3)

式中：ht×Δt為t×Δt時程內的凈面雨量(mm)；fc為穩定下滲率(mm/h)；R總t×Δt為t×Δt時程內降雨的徑流總量(mm)；Pα為雨前土壤含水量(mm)。

根據地區水文氣象資料，可直接獲得fc、Pa值；根據資料，統計實測的R總t×Δt值和P面t×Δt+Pa值，可擬合出兩者間二次經驗關系函數f。

1.1.2 匯流計算

匯流計算方法有等流時線法、單位線法、經驗公式法等。通過對比各方法的特點，瞬時單位線法、水科院推理公式法以其適用性廣、結果精準的特點，常被作為水利水電工程防洪設計的依據。實踐研究表明，對于有較大流域的河流，采用瞬時單位線法計算設計洪水更合理，對于小匯水面積的河溝，通常采用推理公式法計算設計洪水更合理[9]?？紤]到長輸管道穿越的河溝流域面積相差較大，本文針對不同流域面積的河溝采用不同的匯流計算方法。

(1) 瞬時單位線法。采用瞬時單位線法計算暴雨洪水流量的公式為[9]

(4)

(5)

Qt×Δt=h1×Δt·q(1×Δt,(t+1)×Δt)+h2×Δt·q(2×Δt,t×Δt)+…+ht×Δt·q(t×Δt,1×Δt)

(6)

上式中：u(0,t)為無因次單位線；q(t,Δt)為時段單位暴雨洪水流量(km3/h)；Qt×Δt為t×Δt時程內暴雨洪水流量(m3/s);F為至計算斷面的流域面積(km2)；Γ為伽馬函數；n、k為瞬時單位線模式參數，由各流域水文資料獲得，n根據本流域資料總結的n～F關系表取值，k由本流域資料總結的含F、主河道加權平均坡降j(‰)、n、經驗系數等參數的指數函數經驗公式計算得到；其他符號含義同前。

暴雨洪水流量的時程曲線呈現先遞增后遞減的趨勢，遞減達到Qt×Δt=0時，表示地面暴雨洪水徑流終止。

(2) 推理公式法。推理公式法計算洪峰流量的公式為

(7)

(8)

式中：h(rank)t×Δt為各時程凈面雨量按大小排列；L為主河道長度(km)；j為主河道加權平均坡降(‰);τ為匯流時間(s)；m為地區推理公式參數，根據本域流資料總結的含L、j、經驗系數等參數的經驗公式計算得到。

求解匯流時間τ，即聯立求解公式(7)=(8)，得到t=τ的值，此時Qτ·Δt=Qt·Δt為地面洪峰流量。

1.1.3 洪峰流量與水深、流速的計算

所有時程的洪水流量Qt×Δt之最大值即為洪峰流量Qm。河流斷面根據實際情況設為矩形、梯形等形態，記河流水面寬為b(m)，斷面邊坡系數為m(平距與垂距之比)，根據謝才公式，可反解出洪峰流量時的水深h(m)，其計算公式為

(9)

流速v(m/s)的計算公式為

v=Q/ω

(10)

1. 2 河溝道沖刷深度的計算

開挖施工敷設的管道通常埋深淺，管道易受到暴雨沖刷的影響。在管道河溝段沖刷深度計算方法方面，已有研究探究了將墩、臺跨度擴展為整個河寬的一般、局部沖刷深度的計算公式[3]，但該方法在江西省天然氣管道多點試算結果與野外經驗值偏差較大；而水利部《堤防工程設計規范》(GB 50286—2013)[19]中的順壩及平順護岸沖刷深度計算模型，與管道河溝段施工后平整岸坡或修建護坡的物理模型一致[6]，且試算結果與野外經驗值較為吻合，局部沖刷深度的計算公式為

(11)

2 長輸油氣管道河溝段水毀災害氣象風險評價方法

長輸油氣管道輸送介質具有高壓、易燃易爆的特點，故管道地質災害風險評價不僅要考慮地質災害的易發性、易損性，還要結合管道失效后果[3-4,19]。管道河溝段可能修建有護岸、護底等針對水毀災害的工程防治措施，這些工程防護措施能在多大程度上阻止自然災害的發生，也應納入評價指標體系。

本文建立的評價指標體系包括地質災害的自然條件加工程防治措施下的地質災害易發性P、管道易損性C和管道失效后果E三個部分。首先，分別對地質災害易發性、管道易損性建立多個評價指標，每個指標以(0,1]上的數值來表示其可能產生風險的程度，即風險概率指數；然后對管道失效后果E按照公式計算歸一化管道失效后果損失指數；最終將易發性、易損性、失效后果以概率代數、模糊邏輯等方式運算，得到五級劃分的氣象風險等級。

在河溝道沖刷強度易發性評價指標方面，按照輸氣管道設計規范[21]，管道埋深不應小于0.8 m，將管道埋深減去沖刷深度記為h(m)，當h小于0.8 m時管道埋深不足，且h越小風險概率越高，故風險概率指數設為h的函數；而在護岸、護底的有效性評價指標方面，單指標的風險概率指數則難以定量計算；同樣，管道易損性的3個單指標風險概率指數也難以定量計算。針對此種情況，設置每個評價指標所屬的幾種狀態，根據專家經驗，對某單評價指標εi(ε1～6為P1～3和C1～3)不同狀態，如河溝道沖刷強度評價指標P1中P11(h≤0)，P12(0

熵權法是根據指標變異性的大小，即信息熵來確定指標的客觀權重，信息熵越小則提供的信息量越多，更利于評價區分，其重要性越大。某評價指標εi對應評價點的指標值記為xiq，如試算的河溝道沖刷強度P1，第1、2個評價點，則有x11=0.459，x12=0.522；管道敷設方式C1，x41=2(表示溝埋穿越)，x42=3(表示沿河凹岸敷設),…。n個評價點、k個評價指標之間(如3個管道易損性指標之間)，評價指標εi(Pi或Ci)的熵權wi計算過程如下：

xiq的歸一化數據為

(12)

評價指標值比重為

(13)

(14)

評價指標熵權為

(15)

評價指標風險概率指數εij′的計算公式為

(16)

管道河溝段水毀災害易發性、易損性評價指標(Pi和Ci)在不同狀態(Pij和Cij)下的風險概率指數εij′，見表1。

表1 管道河溝段水毀災害易發性、易損性評價指標在 不同狀態下的風險概率指數

災害事件是多項不利因素共同作用造成的，故某評價點管道水毀災害的風險概率(不含失效后果)P(R)，可采用單指標概率指數代數積的方法來計算。地質災害易發性概率Pi、管道易損性概率Ci取值，采用在該指標所屬的某一狀態對應的風險概率指數。P(R)的計算公式為

(17)

管道失效后果E的計算公式[20]為

E=PH×SP×DI×RC

(18)

其中：PH為產品危害系數，根據管道輸送的產品類別(如柴油、原油、汽油、天然氣)，取值為5～10；SP為泄漏系數，計算輸氣管道10 min或輸油管道1 h的泄漏量，并劃分等級，取值為1～5；DI為擴散系數，根據管道周圍巖土體密實程度、與流動水系距離等關系，取值為1.5～5；RC為受體系數，根據管道所在土地用途(如工業區、農村等)以及環境敏感性，取值為0.5～4.9。

管道歸一化失效后果E′的計算公式為

E′=(E-3.75)/1 446.25

(19)

將風險評價等級劃分為五級即低、較低、中、較高、高，不同風險評價等級的確定以風險概率為主、后果損失為輔，后果損失偏低者，按風險概率分級，后果損失偏高者，則其風險等級應當比概率分級有所提高。本文采用模糊邏輯關系，建立風險等級的模糊邏輯關系圖，以此確定風險評價等級，見圖1。

3 實例應用與分析

江西省天然氣管網一期工程全長約1 020 km，為環鄱陽湖8個地市及一些縣區供氣，沿線地形地貌條件復雜多變，降水充沛，該地區多年平均降水量1 638 mm，且年際變化較大，洪水頻發。在管道路線上選取19處開挖施工的一定規模中小河溝道作為評價點，使用GIS軟件，并結合quick bird、worldview衛星高分辨率遙感影像以及江西省測繪局制作的1∶50 000矢量地形圖中的等高線及河流線，提取主河道、流域面積、主河道加權平均坡降等流域參數。由于該輸氣管道投產時間不長，故本文根據江西省水文局暴雨統計資料，以概率統計方法推求N年一遇暴雨降水量，以此對管道河溝段水毀災害氣象風險進行評價。

3. 1 暴雨洪水和河溝道沖刷深度的計算

某時段Δt內設計頻率為a%(1/a%年一遇)的暴雨面雨量P面Δt(a%)計算公式為

P面Δt(a%)=PΔt(年最大)·CvΔt·Kp(a%)·α

(20)

式中:PΔt為年最大Δt時段(年最大)暴雨點雨量均值，如年1 h、3 h最大暴雨點雨量均值，由當地暴雨統計資料相應等值線圖讀??；CvΔt為Δt時段年最大暴雨變差系數，由當地暴雨統計資料相應等值線圖讀??；Kp(a%)為a%所對應的偏差系數Cs=3.5Cv(Cv為變差系數)的皮爾遜Ⅲ型曲線值，即：

(21)

式中：f(x,β,α)為伽馬分布概率密度函數，當β=1時，其為

(22)

由江西省水文局暴雨統計資料可知，該管道經過地區1 h年最大暴雨點雨量均值范圍為40～50 mm，1 h年最大暴雨變差系數點為0.35～5；24 h年最大暴雨點雨量均值范圍為100～140 mm，24 h年最大暴雨變差系數為0.4～0.65。根據暴雨的特征，分別計算時段為1 h、3 h、6 h和24 h，設計頻率分別為2%(50年一遇)、1%(100年一遇)的面暴雨量。在此基礎上，根據對水文調節最不利的雨型時程分配表，設計每個時程的降雨量。

根據江西省天然氣管網一期工程，針對該長輸管道穿越的河溝道流域面積的不同，對于面積大于50 km2的流域采用瞬時單位線法計算并設計洪峰流量，小于或等于50 km2的流域采用推理公式法計算設計洪峰流量，最后根據野外實調的河道斷面資料，計算河溝道沖刷深度。

3. 2 風險評價結果與討論

各評價點的信息根據野外調查記錄、測繪資料、地質資料，按照上述方法進行計算與評價。其中，該管道管徑為508 mm，輸送天然氣壓力為6.3 MPa，按規范推薦資料[20]，后果損失指數中PH均為10，SP統一取值為4，DI統一取值為5。50年一遇暴雨下水毀災害的風險概率P(R)、風險等級和歸一化管道失效后果E′以及100年一遇暴雨下水毀災害的風險概率的計算結果，見表2。

由表2可以看出：

(1) 該評價計算結果與專家經驗判斷結果基本一致。

(2) 在計算過程中，有7處評價點，在50年一遇暴雨條件下管道埋深減去沖刷深度的值處于0～0.8 m之間，但多數評價點處已修建護岸護底等防治措施，工程質量優良，故總風險概率并不高；而金溪縣徐家鄉管道穿越山間河流(19號評價點)，其現狀已出現漿砌石溝底硬化損毀現象，下蝕、溯源侵蝕發育(見圖2和圖3)，經計算該評價點在50年一遇設計暴雨洪水下，洪峰流量為42.2 m3/s，沖刷深度為0.2 m，沖刷后管道埋深為0.3 m，其風險概率P(R)為0.149 9，歸一化管道失效后果E′為0.12，風險等級較高，因此為了降低風險，建議重修溝底硬化以及在管道下游5 m外加修淤積壩，以降低局部坡降，提高護底措施的有效性。

表2 江西省天然氣管網一期工程河溝道水毀災害氣象風險評價結果

圖2 19號評價點較高風險現狀 Fig.2 Situation of the high risk point No.19

圖3 19號評價點較高風險現狀剖面圖Fig.3 Profile map of the high risk point No.19

(3) 在100年一遇暴雨條件下，相比50年一遇暴雨條件，其風險概率P(R)有升高的趨勢，且有2處評價點的風險概率升高明顯。其中，升高最明顯的10號評價點位于安義縣石鼻鎮，管道穿越河流，河床主要由漂石、砂卵石組成，管道上已修建蓋板涵防護，兩岸已有擋墻護岸，但其現狀已出現左側墻基礎掏蝕，右岸部分垮塌現象(見圖4)，經計算，該評價點在50年一遇設計暴雨洪水下其風險概率為0.011 6，在100年一遇設計暴雨洪水下其風險概率為0.016 7，且調查期間護岸被沖毀，護底工程尚完整，經預測暴雨時河流侵蝕能力強，建議汛期應加強巡查，并根據氣象風險適時加強水工保護措施的強度。

圖4 10號評價點中風險現狀Fig.4 Situation of the medium risk point No.10

4 結論與建議

本文探討了適合于長輸油氣管道河溝段暴雨洪水和河溝道沖刷深度計算的方法，并設計了半定量評價指標體系，實現了長輸管道河溝段水毀災害的氣象風險評價。運用此方法，利用推求的 50年、100年一遇的暴雨降水量，對江西省天然氣管網一期工程的19處評價點的水毀災害進行氣象風險評價，其評價結果與實際情況較為符合，說明該方法可為長輸油氣管道等河溝道水毀災害的氣象風險評價提供技術路線參考。

需要注意的是，以概率方式推求N年一遇暴雨降雨量有一定的局限性，汛期時可采用暴雨降雨量實測值或短期預報值，并根據現場實際情況對其他評價指標狀態進行動態修改，從而得到短期預測氣象風險。另外，評價指標的風險概率指數與等級劃分判據等可以隨著實踐認知而優化，以使風險評價結果更加合理。