胡生寶，毛 建，曾躍輝

(1.中石油管道有限責任公司 西氣東輸分公司銀川管理處，寧夏 銀川 750001；2.中石油管道有限責任公司 西氣東輸分公司，上海 200122；3.彭州華潤燃氣有限公司，四川 彭州 611930)

輸氣管道滑坡災害影響因素敏感性分析

胡生寶1，毛 建2，曾躍輝3

(1.中石油管道有限責任公司 西氣東輸分公司銀川管理處，寧夏 銀川 750001；2.中石油管道有限責任公司 西氣東輸分公司，上海 200122；3.彭州華潤燃氣有限公司，四川 彭州 611930)

為保證輸氣管道在滑坡災害下能夠安全運行，需對輸氣管道滑坡災害進行影響因素敏感性分析，以識別出管道應力與影響因素之間的敏感狀況.利用ABAQUS有限元軟件建立輸氣管道滑坡災害數(shù)值分析模型，計算出管道的最大Mises應力值；應用敏感性系數(shù)分析法分析了輸氣管道滑坡災害的主要影響因素，得出主要影響因素與管道應力之間的變化規(guī)律.計算結(jié)果表明，管道在橫穿滑坡區(qū)域時，除管道壁厚與管道應力呈反向變化關(guān)系外，滑坡寬度、內(nèi)摩擦角度、土體黏聚力、管道埋深和管內(nèi)壓力都與管道應力呈正向變化關(guān)系；管道壁厚為影響管道應力的關(guān)鍵因素，其余各因素敏感性從強到弱的排列順序為，滑坡寬度>土體黏聚力>內(nèi)摩擦角>管內(nèi)壓力>管道埋深.

輸氣管道；滑坡；敏感性分析；ABAQUS；數(shù)值模擬

0 引 言

隨著經(jīng)濟的快速發(fā)展，國內(nèi)長距離輸送天然氣管道(以下簡稱管道)已形成橫貫東西，縱貫南北的天然氣輸氣管網(wǎng)，但地形地貌的多樣性也使建成的管道常受到各種地質(zhì)災害的影響，其中滑坡就是影響較嚴重的災害之一[1-2].目前，相關(guān)專家對滑坡災害下管道應力與其影響因素進行了研究，且取得了一定的研究成果[3-6].但以上研究均未涉及對管道應力影響因素的敏感性分析.鑒于此，本研究采用ABAQUS有限元軟件對滑坡災害下管道的受力情況進行模擬，綜合考慮影響滑坡災害下管道受力的各主要因素，并分析出各因素的敏感性順序，識別出關(guān)鍵因素，以期為滑坡災害下輸氣管道的風險評估、防災減災措施的制定提供相關(guān)依據(jù).

1 管道滑坡災害有限元模型

1.1 土體滑動平面模型的建立

目前，管道穿越滑坡體的方式有3種，即橫穿、縱穿與斜穿，而橫向穿越是最常見的滑坡體穿越方式，在該種穿越方式下管道受力最大，最容易失效[7].因此，本研究重點分析管道在滑坡體中橫向穿越的情況.通常，管道在橫穿滑坡體時，因滑動土體摩擦力的作用，管道不僅會因土體下滑彎曲變形而產(chǎn)生彎曲應力，還會導致滑坡體兩端的管道受到拉應力甚至剪切破壞.管道橫穿滑坡體的示意圖如圖1所示.

圖1管道橫穿滑坡體示意圖

本研究以實際工程為例，通過建立ABAQUS有限元模型對相關(guān)問題進行分析.模型中，管道管材為X80鋼，管道直徑取1 016 mm，管道壁厚取17.5 mm，管道埋設深度為2 m，管道的設計運行壓力為8 MPa.管道穿越滑坡體的主要參數(shù)如表1所示.

表1 管道穿越滑坡體的主要參數(shù)

研究發(fā)現(xiàn)，土體滑動對管道的影響可通過在土體端面施加位移載荷的方式代替，并且當管道后方土體長度大于管徑5倍時，其導致管道的極限彎曲影響已不明顯[8].根據(jù)管道直徑大小，本研究所建模型通過在管道后方6 m的土體上施加位移載荷來模擬滑動土體對管道的影響.為更加貼近滑坡實際工況，通常認為非滑動區(qū)的寬度應大于滑動區(qū)寬度的1/4[8].因此，取整個坡體寬度為40 m，其中滑動土體寬度為20 m，非滑動土體寬10 m，土體厚度為8 m，建立如圖2所示的土體滑動平面示意圖.

圖2土體滑動平面示意圖

1.2 管土相互作用模型選擇

目前，常用的管道與土體相互作用模型可分為彈性地基梁模型、土彈簧模型和非線性接觸模型3種.本研究在對這3種相互作用模型的優(yōu)缺點[9]和使用條件進行對比的基礎上，根據(jù)管道滑坡災害模擬的實際情況選擇非線性接觸模型.

1.3 邊界條件設定及網(wǎng)格劃分

1)邊界條件.在滑坡災害下，管道因受到滑動土體的影響，既可能發(fā)生豎直方向上的位移，也可能發(fā)生水平方向上的位移，但管道兩端的非滑動區(qū)土體仍保持穩(wěn)定，限定了管道所有方向的位移和轉(zhuǎn)動.據(jù)此，可將管道兩端面設定為固定約束.

2)網(wǎng)絡劃分.網(wǎng)格劃分的質(zhì)量對管道滑坡災害模擬結(jié)果有著重要影響，網(wǎng)格劃分越細，模擬結(jié)果也與真實情況越接近，但計算成本也就越高.因此，應該從計算精度和計算成本兩方面綜合考慮網(wǎng)格的劃分，本研究管道和土體網(wǎng)格劃分的分別如圖3所示.

(a)管道網(wǎng)格劃分

(b)土體網(wǎng)格劃分

圖3模型對象網(wǎng)格劃分示意圖

1.4 管道應力計算

在滑動區(qū)寬度為20 m的條件下，通過ABAQUS有限元分析軟件計算出管道的Mises應力云圖如圖4所示.

圖4滑坡災害下管道的Mises應力云圖

從圖4可知，管道的最大Mises應力值為466.2 MPa，其位于滑動區(qū)與非滑動區(qū)的過渡段區(qū)域.且由計算出的管道最大Mises應力值小于管道的屈服極限499.5 MPa可知，該輸氣管道目前處于安全運行狀態(tài)[10].

2 滑坡災害下管道應力敏感性分析

管道應力敏感性分析，是指通過改變一個或者多個影響因素，分析管道應力的變化情況，以確定出各因素對于滑坡災害下管道應力之間的影響關(guān)系[11].

2.1 敏感性因素選擇

滑坡災害下的管道應力影響因素可以分成3種類型：第一種是滑坡土體性質(zhì)對管道應力的影響，如滑坡寬度、內(nèi)摩擦角度、土體黏聚力等；第二種是管道自身因素對管道應力的影響，如管道壁厚、管道埋深等；第三種是管內(nèi)介質(zhì)情況對管道應力的影響，如管內(nèi)壓力.

2.2 敏感性分析方法

由于影響滑坡災害下管道應力的因素較多，為達到簡化模擬過程并定量分析影響因素與管道應力間數(shù)值關(guān)系的目的，本研究選擇單因素分析法中的敏感性系數(shù)分析法對影響因素進行敏感性分析.敏感性系數(shù)分析法每次只改變一個影響因素，其他因素保持恒定即參考值，觀察模型響應值f(x)的變化，其中參考值的設定情況如表2所示.敏感性系數(shù)分析法的計算公式[11]為，

(1)

(2)

(3)

式中，ηi為因素對模型的敏感性系數(shù)；xi為因素各級水平值；xi0為因素設定的參考值；f(xi)為因素為第i級水平時模型的響應值；f(xi0)為因素在參考值時模型的響應值.

|ηi|值的大小表示因素敏感性的強弱.|ηi|值越大，表示影響因素的敏感性越強，則其對目標的影響越大；反之，該影響因素的敏感性越弱，則其對目標的影響越小.

表2 敏感性分析的參考值

2.3 管道應力敏感性

在分析管道應力敏感性時，本研究分別以表1和表2的數(shù)據(jù)建立管道滑坡災害有限元模擬模型.計算出管道最大Mises應力值，再分別以參考值的±10%、±20%、±30%、±40%、±50%為變化量改變模型的影響因素，通過ABAQUS有限元軟件計算出管道的最大Mises應力值.計算結(jié)果如表3所示，根據(jù)計算結(jié)果繪制的影響因素敏感性分析圖如圖5所示.

表3 管道最大Mises應力值匯總表

圖5滑坡災害下管道應力影響因素敏感性分析圖

根據(jù)圖5可知，在影響滑坡災害下管道應力的眾多因素中，除管道壁厚與滑坡災害下管道應力呈反向變化關(guān)系外，其余各影響因素都呈正向變化關(guān)系.

在得出各影響因素與管道間的應力變化規(guī)律后，將表3內(nèi)的計算結(jié)果帶入式(1)～式(3)即可求出在相同變化量下各影響因素的敏感性系數(shù).各影響因素的敏感系數(shù)如表4所示，其對應的敏感系數(shù)關(guān)系圖如圖6所示.

表4 不同影響因素的敏感系數(shù)表

2.4 敏感性因素分析與結(jié)果

2.4.1 滑坡寬度.

隨著滑坡寬度的增加，管道應力對滑坡寬度的敏感度系數(shù)減小，其敏感度系數(shù)的變化范圍為0.62～1.05.表明隨著滑坡寬度增大，滑坡寬度對管道應力的影響程度減弱，且減弱效果明顯.

圖6不同影響因素的敏感系數(shù)變化關(guān)系圖

2.4.2 內(nèi)摩擦角度.

隨著滑動土體內(nèi)摩擦角度的增加，管道應力對內(nèi)摩擦角度的敏感度系數(shù)減小，其敏感度系數(shù)的變化范圍為0.40～0.63.表明隨著土體內(nèi)摩擦角度增大，土體內(nèi)摩擦角度對管道應力的影響程度減弱，但減弱效果不太明顯.

2.4.3 土體黏聚力.

隨著滑動土體黏聚力的增加，管道應力對土體黏聚力的敏感度系數(shù)減小，其敏感度系數(shù)的變化范圍為0.25～0.73.表明隨著土體黏聚力增大，土體黏聚力對管道應力的影響程度減弱，且減弱效果明顯.

2.4.4 管道壁厚.

隨著管道壁厚的增加，管道應力對管道壁厚的敏感度系數(shù)的絕對值在減小，其敏感度系數(shù)的變化范圍為-0.57～-1.16.表明隨著管道壁厚增大，管道壁厚對管道應力的影響程度減弱，且減弱效果明顯.

2.4.5 管道埋深.

隨著管道埋深的增加，管道應力對管道埋深的敏感度系數(shù)幾乎保持不變，其敏感度系數(shù)的變化范圍為0.15～0.22.表明隨著管道埋深增大，管道埋深對管道應力的影響程度保持不變.

2.4.6 管內(nèi)壓力.

隨著管內(nèi)壓力的增加，管道應力對管內(nèi)壓力的敏感度系數(shù)幾乎保持不變，其敏感度系數(shù)的變化范圍為0.26～0.31.表明隨著管內(nèi)壓力增大，管內(nèi)壓力對管道應力的影響程度保持不變.

2.4.7 結(jié) 論.

根據(jù)敏感度系數(shù)的計算結(jié)果可以得出各個因素敏感性從強到弱的排列順序為：管道壁厚>滑坡寬度>土體黏聚力>內(nèi)摩擦角>管內(nèi)壓力>管道埋深.因此，當管道處于設計施工階段時，如果管道需從滑坡區(qū)或潛在不穩(wěn)定斜坡區(qū)穿過，一方面可通過適當增加管道壁厚以提高管道承受滑坡災害的能力；另一方面可采用輕質(zhì)回填材料對管道進行回填.當管道處于生產(chǎn)運營階段時，可以在滿足輸氣量的前提下適當降低管道的運行壓力.

3 結(jié) 語

本研究顯示，除管道壁厚與滑坡災害下輸氣管道應力呈反向變化關(guān)系外，滑坡寬度、內(nèi)摩擦角度、土體黏聚力、管道埋深和管內(nèi)壓力都與滑坡災害下輸氣管道應力呈正向變化關(guān)系.輸氣管道壁厚對滑坡災害下管道應力的影響最大，各因素敏感性從強到弱的排列順序為：管道壁厚>滑坡寬度>土體黏聚力>內(nèi)摩擦角>管內(nèi)壓力>管道埋深.本研究認為，將敏感性分析引入到了滑坡災害下管道的應力分析中，在一定程度上可定量分析各因素與輸氣管道應力之間的規(guī)律，此為制定滑坡災害管道的防治措施提供了一定的借鑒依據(jù).

[1]姚志祥，龍軍．格拉管道一次橫向滑坡斷裂跑油搶修實踐[J]．管道技術(shù)與設備，2004，12(1)：31-32．

[2]帥健，王曉霖，左尚志．地質(zhì)災害作用下管道的破壞行為與防護對策[J]．焊管，2008，31(5)：9-15，93．

[3]林冬，雷宇，許可方，等．橫向滑坡對管道的影響試驗[J]．石油學報，2011，32(4)：728-732．

[4]焦中良，谷海威，郭杰，等．滑坡條件下埋地管道的應力分析[J]．煤氣與熱力，2014，34(12)：11-16．

[5]練章富，李風雷．滑坡帶埋地管道力學強度分析[J]．西南石油大學學報(自然科學版)，2014，36(2)：165-170．

[6]郝建斌，劉建平，荊宏遠，等．橫穿狀態(tài)下滑坡對管道推力的計算[J]．石油學報，2012，33(6)：1093-1097．

[7]吳銳，梅永貴，鄧清祿，等．滑坡作用下輸氣管道受力分析[J]．建筑科學與工程學報，2014，31(3)：105-111．

[8]張伯君．山體滑坡區(qū)域內(nèi)長輸埋地油氣管道強度研究[D]．杭州：浙江大學，2013．

[9]劉慧．滑坡作用下埋地管線反應分析[D]．大連：大連理工大學，2008．

[10]中國石油天然氣集團公司.輸氣管道工程設計規(guī)范：GB50251-2015[S]．北京：中國計劃出版社，2015.

[11]蔡毅，邢巖，胡丹．敏感性分析綜述[J]．北京師范大學學報(自然科學版)，2008，44(1)：9-16．

SensitivityAnalysisofInfluencingFactorsofLandslideinGasPipeline

HUShengbao1,MAOJian2,ZENGYuehui3

(1.Yinchuan Management Office, West-East Gas Pipeline Company，Ningxia 750001, China; 2.West-East Gas Pipeline Company, China Petroleum Pipeline Co., Ltd., Shanghai 200122, China; 3.Pengzhou China Resources Gas Co., Ltd., Pengzhou 610500, China)

In order to ensure the safe operation of the gas pipeline in the landslide disaster,it is necessary to analyze the sensitivity of the influencing factors of the landslide hazard of the gas pipeline to identify the sensitive relationship between the pipeline stress and the influencing factors.The numerical analysis model of landslide hazard is established by using ABAQUS finite element software to calculate the maximum Mises stress value of the pipeline.The main factors influencing the landslide hazard of gas pipeline are analyzed by sensitivity coefficient analysis,and the analysis identifies the main influencing factors and the law of change between pipeline stresses.The results show that the landslide width,internal friction angle, soil cohesion,pipeline depth and pipe pressure are all positive with the pipeline stress when the pipe crosses the landslide,while the pipe wall thickness is inversely related to the pipe stress.The relationship between the thickness of the pipe and the wall is the key factor affecting the stress of the pipeline.The sensitivity of the other factors rank as the following:landslide width,the soil cohesion,the internal friction angle,the pipe pressure and the pipeline depth.

gas pipeline;landslide;sensitivity analysis;ABAQUS;numerical simulation

TE973；P642.22

A

1004-5422(2017)04-0427-05

2017-07-06.

胡生寶(1975 — )，男，工程師，從事油氣管道完整性相關(guān)技術(shù)研究.