何禮元1 賈彥琨2 劉琨2 童驍2 梁伍1

1.北京中地華安環(huán)境工程有限公司；2.國家管網(wǎng)集團西氣東輸南昌輸氣分公司 

 

摘要：以某長(cháng)輸管道水毀地質(zhì)災害為研究對象，采用Logistic回歸模型分別完成地質(zhì)災害易發(fā)性定量評價(jià)和影響條件地質(zhì)災害易發(fā)性定量評價(jià)，再綜合兩個(gè)定量評價(jià)結果獲得最終地質(zhì)災害易發(fā)性定量評價(jià)結果，據此將研究區域地質(zhì)災害易發(fā)程度劃分為高易發(fā)區、中易發(fā)區、低易發(fā)區和非易發(fā)區，對油氣管道防災減災有一定指導意義。

關(guān)鍵詞：油氣管道；地質(zhì)災害；水毀；邏輯回歸模型；易發(fā)性預測分區

 

我國長(cháng)輸油氣管道沿線(xiàn)地質(zhì)環(huán)境和自然條件十分復雜，部分地段地質(zhì)災害嚴重，對油氣管道建設及運營(yíng)構成較大制約和影響。根據近10年的統計，由地質(zhì)災害造成的油氣管道運行失效約占11%，僅次于第三方施工損壞和打孔盜油破壞，嚴重威脅著(zhù)管道的運行安全。

本文以某長(cháng)輸管道沿線(xiàn)為研究區，依據其年度地質(zhì)災害專(zhuān)業(yè)化排查出的地質(zhì)災害點(diǎn)及相關(guān)數據，根據野外調查數據分析地質(zhì)災害的發(fā)育條件、分布規律及致災因素，選取有效的評價(jià)指標。采用邏輯回歸模型構建易發(fā)性評價(jià)模型，預測災害易發(fā)性概率。根據R方檢驗和顯著(zhù)性檢驗等結果判斷模型精度和可靠性。最后對地質(zhì)災害易發(fā)區域進(jìn)行劃分，為油氣管道沿線(xiàn)的減災防災及規劃管理提供一定程度的參考。

1  Logistic回歸模型介紹

Logistic函數由比利時(shí)數學(xué)家Pierre François Verhulst在1938年提出[1-3]，形式如下：

其中，N（t）表示t時(shí)刻個(gè)體數，r表示內在增長(cháng)率，K表示承載力，化簡(jiǎn)為：

當因變量為二值變量時(shí)，0表示事件未發(fā)生；1表示事件發(fā)生。定義P為事件發(fā)生的概率，未發(fā)生即為1-P，把P看作為自變量Xi的線(xiàn)性函數。P對于Xi在P=0或P=1附近是不敏感且緩慢的，非線(xiàn)性程度較高，于是尋求一個(gè)關(guān)于P的函數，使得它在P=0或P=1附近變化幅度較大，引入Logistic函數變換，解得邏輯回歸模型方程為：

式中：P為研究區地質(zhì)災害發(fā)生概率；e為邏輯回歸方程自然常數；b0為邏輯回歸方程截距（邏輯回歸方程的一個(gè)常數項）；bi（i＝1,2,3…,n）為邏輯回歸模型中對應的評價(jià)因子X(jué)i（i＝1,2,3…,n）的邏輯回歸系數。

2  擬合分析

選取的指標有：高程X1、管道走向X2、積水影響長(cháng)度X3、地表徑流與管道關(guān)系X4、地形地貌X5、管道埋深X6。將數據帶入Logistic回歸模型進(jìn)行擬合，其中地表徑流與管道關(guān)系定義為分類(lèi)變量，0表示否，1表示是。分類(lèi)截止值設定為0.5。擬合后預測正確度統計如表 1所示，結果可知Logistic預測正確率高達81.1%，預測效果較好。

表 1 擬合優(yōu)度統計量

Cox & Snell R方和Nagelkerke R方檢驗結果參數表示模型的變量對因變量的解釋能力，取值為0～1，越接近1解釋能力越好。如表 2所示檢驗結果均大于0.5，說(shuō)明具有較好的擬合程度。

表 2 擬合結果檢驗

擬合結果中自變量系數如表 3所示，B項為每個(gè)自變量所對應系數，S.E項為標準差，Sig項為顯著(zhù)性試驗。在易發(fā)性評價(jià)過(guò)程中，邏輯回歸模型運算結果中的B值為各評價(jià)因子對災害發(fā)生的貢獻度，其值越大說(shuō)明該因子對災害發(fā)生的影響越大。通過(guò)比較Wald和Sig值來(lái)判斷各評價(jià)指標的顯著(zhù)性，Wald值越大或Sig值越小，顯著(zhù)性就越高。根據顯著(zhù)性檢驗結果，可以看出地表徑流與管道關(guān)系、管道埋深等因素對災害具有較高的敏感性。

表 3 自變量系數及結果檢驗

將常數、X1－X6自變量系數帶入式（3），q＝2.592－0.009X1＋0.006X2－0.002X3－1.359X4＋0.127X5＋0.326X6。最終得出水毀災害易發(fā)性概率計算模型如下：

3  易發(fā)性[4]概率預測

研究樣本為某長(cháng)輸管道年度地質(zhì)災害專(zhuān)業(yè)化排查數據，地質(zhì)災害隱患點(diǎn)主要包括河溝道水毀、坡面水毀、臺田地水毀等。野外調查時(shí)數據源采集參數包括里程、經(jīng)度、緯度、高程、災害點(diǎn)類(lèi)型、土質(zhì)類(lèi)型、植被類(lèi)型、管道走向、積水影響長(cháng)度、地表徑流與管道關(guān)系、地形地貌、管道狀況、管道埋深等。如表 4所示，選取高程X1、管道走向X2、積水影響長(cháng)度X3、地表徑流與管道關(guān)系X4、地形地貌X5、管道埋深X6六個(gè)指標數據進(jìn)行概率預測。

表 4 管道沿線(xiàn)現場(chǎng)原始數據（以10組數據為例）

將原始數據在表格中進(jìn)行整合，輸入IBM SPSS Statistics中，經(jīng)計算后將類(lèi)別數據輸出為量化數據，將SPSS輸出結果帶入易發(fā)性概率計算公式中，計算結果如表 5所示。

表 5 原始數據SPSS輸出

4  地質(zhì)災害易發(fā)區[5-6]劃分

按照“區內相似、區際差異”的原則，即在同一類(lèi)型的區內，地質(zhì)環(huán)境背景條件、主要誘發(fā)條件和地質(zhì)災害發(fā)育特征應基本相似，而不同類(lèi)型的區內，則具有明顯的差異性。依據定量評價(jià)結果，借助GIS（地理信息系統）技術(shù)，用電子信息化手段合并易發(fā)性相同或者相近的概率單元格，對地質(zhì)災害易發(fā)性程度進(jìn)行等級分區，分區主要特征如表 6所示。

表 6 地質(zhì)災害易發(fā)程度分區主要特征

5  總結與建議

（1）通過(guò)Logistic回歸模型預測油氣管道水毀地質(zhì)災害易發(fā)性概率效果較好，根據預測結果對易發(fā)程度分區，為管道地質(zhì)災害日常管理劃出了重點(diǎn)，有利于管道管理單位將有限的資源合理利用，立足于治早治小，對油氣管道水毀地質(zhì)災害防災減災有一定指導意義。

（2）由于地質(zhì)災害的發(fā)育是地質(zhì)環(huán)境變化和影響因素綜合作用的結果，是不斷變化和發(fā)展的，在同一或臨近區域的影響因素相似性高。因此本文介紹的關(guān)于Logistic回歸模型水毀災害易發(fā)性概率計算結果適用于長(cháng)輸油氣管道的易發(fā)性分區。同時(shí)，該方法所需采集的樣本較多，難以在日常管理中實(shí)施，主要依靠承包商單位周期性地質(zhì)災害專(zhuān)業(yè)化排查成果進(jìn)行易發(fā)性預測分區。

（3）地質(zhì)災害的調查、防治將是一項長(cháng)期的工作，下一步將結合區域監測預警信息，減少現場(chǎng)數據采集工作量，實(shí)現短周期地質(zhì)災害易發(fā)性分區，為長(cháng)輸油氣管道的安全運營(yíng)提供長(cháng)期保障。

 

參考文獻：

[1]田春山，劉希林，汪佳.基于 CF 和 Logistic 回歸模型的廣東省地質(zhì)災害易發(fā)性評價(jià)[J].水文地質(zhì)工程地質(zhì)，2016(11)：154-170.

[2]鐘畝鋒.基于確定性系數與Logistic回歸模型的地質(zhì)災害易發(fā)性評價(jià)[J].世界有色金屬，2020(01)：186-190.

[3]樊芷吟，茍曉峰，秦明月，范強，余建樂(lè )，趙建軍.基于信息量模型與 Logistic 回歸模型耦合的地質(zhì)災害易發(fā)性評價(jià) [J].工程地質(zhì)學(xué)報，2018(04)：340-347.

[4]鐘蔭乾.地質(zhì)災害易發(fā)性評價(jià)[J].湖北地礦，2002，16(4)：81-86.

[5]尚慧，倪萬(wàn)魁，程花.斜坡單元劃分在彭陽(yáng)縣地質(zhì)災害危險性區劃的應用[J].中國水土保持，2011，3(2)：48-50.

[6]苗雨.基于GIS和Logistic模型的地質(zhì)災害危險性區劃研究[D].西安：長(cháng)安大學(xué)，2010：1-59.

作者簡(jiǎn)介：何禮元，1995年生，2018年畢業(yè)于安徽理工大學(xué)地質(zhì)工程專(zhuān)業(yè)，大學(xué)本科學(xué)歷，助理工程師，主要從事長(cháng)輸油氣管道地質(zhì)災害勘查、評估、監測預警工作。聯(lián)系方式：13135548788，757279399@qq.com。