魏良帥， 陳國輝， 王正平， 王佃明

(1.中國地質(zhì)科學(xué)院探礦工藝研究所,四川 成都 611734； 2.中國地質(zhì)調(diào)查局地質(zhì)災(zāi)害防治技術(shù)中心,四川 成都 611734； 3.四川省地質(zhì)工程勘察院，四川 成都 610072； 4.中國石油化工股份有限公司中原油田普光分公司，四川 達州 635002)

油氣管道經(jīng)常不可避免地穿越一些地形地質(zhì)條件復(fù)雜區(qū)，而這些地區(qū)往往又是地質(zhì)災(zāi)害的高發(fā)區(qū)，不同類型的地質(zhì)災(zāi)害以不同的方式威脅著管道的安全運營。管道地質(zhì)災(zāi)害的發(fā)生往往會導(dǎo)致管道輸送介質(zhì)的泄漏，帶來巨大的財產(chǎn)損失與環(huán)境污染。據(jù)統(tǒng)計：美國天然氣管道事故的8.5%是由地質(zhì)災(zāi)害導(dǎo)致的；西歐管道地質(zhì)災(zāi)害造成的事故率為7%；加拿大管道地質(zhì)災(zāi)害造成的事故占12%；南美洲安第斯山區(qū)的Andean天然氣管道，地質(zhì)災(zāi)害導(dǎo)致的事故率高達50%[1]。因此，針對管道沿線地質(zhì)災(zāi)害進行調(diào)查與風(fēng)險評價，對保障油氣管道安全運營，減少地質(zhì)災(zāi)害對管道的危害具有重要意義。

1 管道沿線地質(zhì)災(zāi)害發(fā)育概況

普光氣田位于四川省達州市宣漢縣普光鎮(zhèn)，屬超深、高含硫、高壓、復(fù)雜山地氣田。普光氣田是中國目前發(fā)現(xiàn)的最大規(guī)模海相整裝高含硫氣田，截至目前已探明天然氣地質(zhì)儲量4122億m3。根據(jù)中國石油天然氣行業(yè)氣藏分類標(biāo)準(zhǔn)，屬于特大氣藏。為有效開發(fā)利用高含硫天然氣資源，國務(wù)院將“川氣東送”工程列為“十一五”國家重大工程，將普光氣田作為工程的主供氣源。作為氣田內(nèi)集輸氣工程的天然氣管線共計64.49 km。管道沿線地形起伏大，地形坡度一般25°以上，相對高差達800 m。管道穿越區(qū)屬亞熱帶濕潤季風(fēng)氣候，年均降水量達1213.5 mm；管道沿線地質(zhì)環(huán)境條件整體較為復(fù)雜，為地質(zhì)災(zāi)害的發(fā)育創(chuàng)造了有利的條件，也為管道的安全運營帶來了隱患[2]。

根據(jù)2017年汛期管道沿線地災(zāi)排查資料顯示[2]，此次排查共發(fā)現(xiàn)各類地質(zhì)災(zāi)害隱患點63處，類型涉及滑坡、崩塌(危巖)、水毀(坡面水毀、河溝道水毀、臺田地水毀)、地面塌陷和潛在不穩(wěn)定斜坡共7種災(zāi)害類型，其中滑13坡處，不穩(wěn)定斜坡22處，崩塌17處，水毀9處，地面塌陷1處(圖1)。

圖1 管道沿線地質(zhì)災(zāi)害發(fā)育分布圖

2 地質(zhì)災(zāi)害風(fēng)險評價模型及評價指標(biāo)體系

2.1 地質(zhì)災(zāi)害風(fēng)險評價模型

通常情況下，風(fēng)險評價是指對不良結(jié)果或事件發(fā)生的概率進行定量描述的過程。亦或說，風(fēng)險評價是對一特定時期內(nèi)健康、安全、環(huán)境、生態(tài)等受損傷的可能性大小及可能程度作出的評估過程[3-4]。對區(qū)域管道地質(zhì)災(zāi)害風(fēng)險而言，風(fēng)險評價可以定義為：是對地質(zhì)災(zāi)害發(fā)生的危險性和災(zāi)害損失(包括管道因地質(zhì)災(zāi)害破壞泄漏造成的后果損失)的可能性做出的綜合性分析評價。它是地質(zhì)災(zāi)害危險性、易損性的組合，災(zāi)害的易損不僅包括災(zāi)害本身造成的損失，還包括管道因地質(zhì)災(zāi)害破壞泄漏造成的后果損失。因此區(qū)域管道地質(zhì)災(zāi)害風(fēng)險要綜合考慮危險性與易損性兩方面的影響[4]：

風(fēng)險(R)=危險性(H)×易損性(V)

本次風(fēng)險評價采用該模型作為區(qū)域管道地質(zhì)災(zāi)害風(fēng)險評價模型。

地質(zhì)災(zāi)害危險性、易損性評價采用基于GIS多因子綜合疊加模型[5-9]，其計算公式如下：

H=ωihi

(1)

V=ωjvj

(2)

式中：H——地質(zhì)災(zāi)害危險值；ωi——第i個危險性評價指標(biāo)的權(quán)重；hi——第i個危險性評價指標(biāo)的狀態(tài)值；V——管道地質(zhì)災(zāi)害的易損值；ωj——第j個易損性評價指標(biāo)的權(quán)重；vj——第j個易損性評價指標(biāo)的狀態(tài)值。

本次評價以管道為中心線，管道兩側(cè)1 km范圍作為評價區(qū)，采用30 m×30 m的規(guī)則柵格格網(wǎng)單元作為此次管道地質(zhì)災(zāi)害風(fēng)險評價的最小單元。

2.2 地質(zhì)災(zāi)害風(fēng)險評價指標(biāo)體系

地質(zhì)災(zāi)害危險性評價指標(biāo)從地質(zhì)環(huán)境及誘發(fā)因素角度盡量全面地考慮影響地質(zhì)災(zāi)害發(fā)生的各種影響因素，從宏觀上主要分為兩類：基本因素和影響(誘發(fā))因素?；疽蛩厥侵复_定地質(zhì)災(zāi)害發(fā)生背景的基本地質(zhì)因素，如地形、地層巖性、地質(zhì)構(gòu)造等。影響因素指影響和誘發(fā)地質(zhì)災(zāi)害向不利方向演化的各種外動力和人類活動因素。根據(jù)野外調(diào)查及室內(nèi)分析研究，綜合考慮現(xiàn)有數(shù)據(jù)與GIS技術(shù)功能的局限性等因素，采用兩兩對比法對指標(biāo)篩選優(yōu)化，最終選取坡度、高差、匯流面積、災(zāi)害點密度、斜坡結(jié)構(gòu)類型、坡面形態(tài)等6個因子構(gòu)建管道沿線地質(zhì)災(zāi)害危險性評價指標(biāo)體系(圖2)。

圖2 區(qū)域管道地質(zhì)災(zāi)害危險性評價指標(biāo)體系

區(qū)域管道地質(zhì)災(zāi)害易損性不僅僅包括地質(zhì)災(zāi)害本身對承載體的破壞，還包括管道在遭受地質(zhì)災(zāi)害破壞后，導(dǎo)致的管道破裂造成油氣泄漏對周圍人口、道路、環(huán)境等的破壞。分析輸氣管道泄漏后對周圍人口、道路、環(huán)境等造成的影響，綜合考慮現(xiàn)有數(shù)據(jù)與GIS技術(shù)功能的局限性等條件，最終構(gòu)建的區(qū)域管道地質(zhì)災(zāi)害易損性評價指標(biāo)體系如圖3所示。

圖3 區(qū)域管道地質(zhì)災(zāi)害易損性評價指標(biāo)體系

2.3 評價指標(biāo)權(quán)重

本次指標(biāo)權(quán)重的確定采用貢獻率模型，分析管道沿線地質(zhì)災(zāi)害發(fā)育特征及易損程度，對評價指標(biāo)在每個災(zāi)害的影響程度進行打分(1分為評價指標(biāo)對災(zāi)害的影響小或易損小，2分影響中等或易損中等，3分影響大或易損大)，然后匯總各災(zāi)害對每種指標(biāo)的貢獻總值，通過歸一化處理，最終確定的各因子的權(quán)重(表1、表2)[10]。

(3)

(4)

式中：rij——第j災(zāi)害樣本中第i個評價指標(biāo)分值；n——災(zāi)害樣本的個數(shù)；ri——第i個評價指標(biāo)的貢獻總值；m——評價指標(biāo)的個數(shù)；Si——第i個評價指標(biāo)的權(quán)重。

表1 地質(zhì)災(zāi)害危險性評價指標(biāo)權(quán)重

表2 地質(zhì)災(zāi)害易損性評價指標(biāo)權(quán)重

3 評價指標(biāo)的量化與分級

3.1 危險性評價指標(biāo)的量化與分級

在建立評價指標(biāo)體系后，進一步要確定他們的分級，一般情況下，指標(biāo)的分級應(yīng)與危險性分級相對應(yīng)。地質(zhì)災(zāi)害危險性分為高危險、較高危險區(qū)、中等危險區(qū)、較低危險區(qū)、低危險區(qū)五級，分別賦值為5、4、3、2、1，相應(yīng)地將評價指標(biāo)也用五級來描述。危險性評價指標(biāo)分級賦值標(biāo)準(zhǔn)如表3所示，分級評價結(jié)果參見圖4。

表3 地質(zhì)災(zāi)害危險性評價指標(biāo)分級標(biāo)準(zhǔn)

3.2 易損性評價指標(biāo)的量化與分級

易損性評價指標(biāo)同樣相應(yīng)地將評價指標(biāo)也用五級來分級描述，分級賦值標(biāo)準(zhǔn)如表4所示。分級評價結(jié)果參見圖5。

表4 地質(zhì)災(zāi)害易損性評價指標(biāo)分級標(biāo)準(zhǔn)

圖4 管道地質(zhì)災(zāi)害危險分區(qū)

圖5 管道地質(zhì)災(zāi)害易損分區(qū)

4 管道地質(zhì)災(zāi)害風(fēng)險分區(qū)

4.1 管道地質(zhì)災(zāi)害風(fēng)險評價

基于GIS多因子綜合疊加模型，利用柵格計算器將獲取到的地質(zhì)災(zāi)害危險性、易損性各指標(biāo)與其對應(yīng)的權(quán)重進行加權(quán)計算，然后綜合兩方面的評價成果，最終完成區(qū)域管道地質(zhì)災(zāi)害風(fēng)險值計算。

合理地確定地質(zhì)災(zāi)害風(fēng)險等級界線值是風(fēng)險區(qū)劃的關(guān)鍵環(huán)節(jié)之一，一般采用自然斷點法和等間距法[11-12]，本次采用前者，將管道沿線地質(zhì)災(zāi)害風(fēng)險劃分為地質(zhì)災(zāi)害高風(fēng)險區(qū)、較高風(fēng)險區(qū)、中等風(fēng)險區(qū)、較低風(fēng)險區(qū)和低風(fēng)險區(qū)五個不同等級[13-14],各單元確定的地質(zhì)災(zāi)害風(fēng)險性等級標(biāo)準(zhǔn)參見表5，最終得到管道沿線地質(zhì)災(zāi)害風(fēng)險分區(qū)圖(圖6)。

表5 地質(zhì)災(zāi)害危險等級分區(qū)

圖6 管道地質(zhì)災(zāi)害風(fēng)險分區(qū)、分段

4.2 管道地質(zhì)災(zāi)害風(fēng)險分段

根據(jù)地質(zhì)災(zāi)害風(fēng)險評價分區(qū)圖對管道進行地質(zhì)災(zāi)害風(fēng)險分段，將管道劃分成了50個風(fēng)險等級段，其中：高風(fēng)險段2段，較高風(fēng)險段14段，中風(fēng)險段20段，較低風(fēng)險段14段(表6)。

表6 管道沿線地質(zhì)災(zāi)害風(fēng)險等級分段評價統(tǒng)計

4.3 評價結(jié)果分析

將現(xiàn)有地質(zhì)災(zāi)害點圖層與各管道地質(zhì)災(zāi)害危險性分區(qū)圖疊加，結(jié)果顯示災(zāi)害大多落在中等～高風(fēng)險區(qū)段。同時，將評價分區(qū)分段結(jié)果與野外調(diào)查分段結(jié)果對比，評價結(jié)果與調(diào)查評價結(jié)果相符(圖6)。

5 結(jié)論

(1)區(qū)域管道地質(zhì)災(zāi)害風(fēng)險是地質(zhì)災(zāi)害危險性、易損性的組合，災(zāi)害的易損不僅包括災(zāi)害本身造成的損失，還包括管道因地質(zhì)災(zāi)害破壞泄漏造成的后果損失。

(2)采用貢獻率模型作為確權(quán)方法，利用管道沿線地質(zhì)災(zāi)害作為統(tǒng)計樣本點，有效避免了人為因素的影響。

(3)完成了普光氣田地面集輸氣管線地質(zhì)災(zāi)害風(fēng)險評價，將管道沿線地質(zhì)災(zāi)害風(fēng)險劃分為地質(zhì)災(zāi)害高風(fēng)險區(qū)、較高風(fēng)險區(qū)、中等風(fēng)險區(qū)、較低風(fēng)險區(qū)和低風(fēng)險區(qū)五個不同等級,并根據(jù)地質(zhì)災(zāi)害風(fēng)險評價分區(qū)圖對管道進行地質(zhì)災(zāi)害風(fēng)險分段，將管道劃分成了53個風(fēng)險等級段，其中：高風(fēng)險段2段，較高風(fēng)險段14段，中風(fēng)險段20段，較低風(fēng)險段14段。