基于離散滑坡位移條件下油氣管道受力分析*

李明

(濟南大學；中國石油天然氣管道工程有限公司)

0 引 言

我國地質條件復雜多樣，尤其在西南各省山區和西北黃土高原地區，雨季來臨時滑坡地質災害較為常見。油氣長輸管道具有距離長、沿線分輸點多和覆蓋面廣的特點，兼顧避繞人口密集區和城鎮規劃區的需求，在特殊情況下受場地條件限制時，不可避免穿越滑坡體，從而受到滑坡地質災害的威脅。據統計，中緬油氣管道在2014—2018年期間發生滑坡22起[1]，西氣東輸管道沿線統計滑坡155處，澀寧蘭輸氣管道沿線統計滑坡13處[2-3]，因此長輸油氣管道受滑坡體影響的區段很多。當坡體穩定性遭到破壞而發生滑動時，管道在坡體滑動位移的作用下輕則會產生變形，重則會造成斷裂、火災和爆炸。為保障油氣管道安全，需要對穿越滑坡區段油氣管道的受力狀況進行分析。

以往的研究成果中，陶勇寅等[4]將管道簡化為彈性地基梁，考慮了管土之間的剪切應力，并給出了彈性地基梁的撓曲線微分方程，該方程可用于埋地管道彎矩、剪力和應力的計算；鄧道明等[5]在滑坡地質災害對管道破壞機理研究的基礎上，進一步發展了滑坡作用下管道的內力和變形計算公式；劉慧[6]在前人分析的基礎上，參考斷層錯動的相關研究，給出了管道在軸向滑坡、橫向滑坡及圓弧滑坡作用下的解析方法；王聯偉等[7]基于土彈簧模型對管道滑坡力學影響因素進行了分析，討論了縱向滑坡和橫向滑坡作用下管道受力情況。

目前大多數研究都是基于固定滑坡位移條件下的管道應力分析，考慮到滑坡體不同部位的滑動位移并不相同，采用固定滑坡位移對管道進行應力分析所得計算結果還不夠精確，因此需要首先對滑坡體沿管道敷設方向各點的實際位移進行求解，然后再對管道相應位置施加離散滑坡位移，發現管道受力狀況與滑坡體的滑動規律更為相符，并采用該方法分別對穿越縱向滑坡體管道和穿越橫向滑坡體管道進行了受力分析。

1 模型建立

1.1 建立滑坡體模型

建立滑坡體模型用于求解沿管道敷設方向各點的滑動位移。滑坡體土壤本構關系選用Mohr-Coulomb模型，該模型假設坡體剪應力達到材料抗剪強度極限時即會發生破壞，可以用作模擬滑坡體土壤材料，其抗剪強度與其所受正應力呈線性關系[8]，表達式為：

τ=c+σtanφ

(1)

式中：τ為土壤材料的剪切強度，kPa；σ為切點的正應力，kPa；c為黏聚力，kPa；φ為內摩擦角，(°)。

采用有限元強度折減法模擬坡體的失穩滑動過程，將坡體所能提供的最大抗剪強度與外載荷在坡體內所產生的實際剪應力之比定義為強度折減系數，通過調整強度折減系數來降低巖土強度，直至坡體在外載荷作用下發生失穩滑動[9-11]。對于Mohr-Coulomb模型，強度折減系數可以表示為：

cm=c/Fr

(2)

φm=arctan(tanφ/Fr)

(3)

式中：cm為強度折減之后的黏聚力，kPa；φm為強度折減之后的內摩擦角，(°)；Fr為強度折減系數。

坡體是否發生失穩滑動的判定準則有3類[12-14]：①以有限元計算是否收斂作為判定準則；②以等效塑性應變自坡腳至坡頂是否發生貫通作為判定準則;③以坡面節點位移是否發生突變作為判定準則。其中坡體是否發生塑形貫通具有直觀、易判斷的特點，因此這里將其作為坡體是否發生失穩滑動的判定準則，然后提取沿管道敷設方向各點的滑動位移作為離散滑坡位移。建立滑坡體模型，坡頂長a，坡角α，滑坡體坡面長L，模型高h，坡體基礎長b，滑坡體模型如圖1所示。

圖1 滑坡體模型Fig.1 Sliding mass model

1.2 建立管土相互作用模型

建立管土相互作用模型用于求解離散滑坡位移條件下管道受力狀況。管道材料具有非線性特點，其應力應變曲線可以根據GB/T 50470—2017中Ramberg-Osgood方程進行擬合[15]。管道本體采用PIPE20管單元模擬，該單元是具有拉壓、彎曲和扭轉性能的單軸單元，每個節點有6個自由度，即沿節點坐標系X、Y、Z方向的平動自由度和繞X、Y、Z軸的轉動自由度，具有塑形、蠕變和膨脹功能[16]。土壤本體采用COMBIN39土彈簧單元模型。GB/T 50470—2017中列出了三軸土彈簧的計算公式，包含管道軸向、側向、垂直向上和垂直向下土彈簧[17-21]。管土相互作用模型考慮了材料和幾何大變形的非線性特征，可以實現離散滑坡位移條件下管道受力規律及失效因素分析[22-25]。

當管道縱向穿越滑坡體時，管道可分為坡底固定土體段、縱向穿越滑坡體段和坡頂固定土體段三部分，其中縱向穿越滑坡體段長度Lp與滑坡體模型中坡面長度L相等，這樣滑坡體坡面各點滑動位移與相應位置管道所受位移載荷就可以一一對應，實現在管道上施加離散滑坡位移過程。當管道橫向穿越滑坡體時，管道可分為左側固定土體段、橫向穿越滑坡體段和右側固定土體段三部分，其中橫向穿越滑坡體段長度為Wp，由于滑坡體上部、中部和下部位移分布并不相同，所以管道橫向穿越滑坡體不同部位時需要分別進行受力分析。管道縱向和橫向穿越滑坡體幾何模型分別如圖2和圖3所示。

圖2 埋地管道縱向穿越滑坡體幾何模型Fig.2 Geometric model for buried pipeline to cross sliding mass vertically

圖3 埋地管道橫向穿越滑坡體幾何模型Fig.3 Geometric model for buried pipeline to cross sliding mass horizontally

2 實例分析

2.1 基礎數據

某一均勻土質滑坡體，L=100 m，α=45°，h=100 m，a=80 m，b=200 m，土壤容重γ=20 kN/m3，c=12.38 kPa，φ=20°，求解坡體滑動時坡面各點水平、豎向和總位移。某輸氣管道穿越該滑坡體，管道規格為D1219×27.5 mm，材質為X80，管道內壓在0～10 MPa之間，管道運行時與下溝回填時的溫差在0～50 ℃之間變化，管頂埋深1.2 m，當管道縱向穿越滑坡體時Lp=100 m，當管道橫向穿越滑坡體時Wp=100 m，分別對管道縱向和橫向穿越該滑坡體時進行受力分析。根據土壤參數計算土彈簧參數，結果如表1所示。

表1 三軸土彈簧參數Table 1 Parameters of triaxial soil spring

2.2 離散滑坡位移

當坡體發生滑動時，滑動面近似為圓弧形，自坡腳至坡頂沿管道敷設方向建立一條路徑，該路徑上土壤位移即為管周離散滑坡位移，其長度為100 m。根據有限元計算結果，沿該路徑的豎向滑動位移絕對值逐漸增大，近似呈直線增長，說明滑坡體頂部主要發生豎向移動，而在滑坡體底部豎向位移較小；沿該路徑的水平位移絕對值呈現先增長后降低的趨勢，說明坡腳滑動面下部土體水平位移較小，而滑動面上部土體水平位移較大；滑坡體總位移呈現先增長后略有降低的趨勢，說明滑動面上部土體總位移變化不大，位移最大點出現在坡腳滑動面附近。滑坡體位移分布如圖4所示。滑坡體沿管道敷設方向位移的變化曲線如圖5所示。

圖4 滑坡位移分布云圖Fig.4 Cloud chart for landslide displacement distribution

圖5 滑坡體沿管道敷設方向位移的變化曲線Fig.5 Displacement variation curve of sliding mass in pipelining direction

2.3 縱向穿越滑坡體管道受力分析

當管道縱向穿越滑坡體時，將上述離散滑坡位移施加到對應位置的管道上，在管道幾何模型中該區段管道長度為200～300 m。根據有限元計算結果，管道Mises應力極值點有4個，分別處于坡腳左側水平段、坡腳彎管、滑動面附近和坡頂附近；管道最大Mises應力點在坡腳彎管處，主要由于該點不但承受了上部土體的滑動位移載荷，而且附近土體的滑動剪切作用還在該點產生了較大的彎矩；滑動面附近存在一處Mises應力極值點，主要由于土體的滑動剪切作用使管道產生了較強的應力和較大的彎曲變形；坡腳左側水平段管道也有一處應力極值點，主要由于坡腳彎管在土壤載荷作用下產生一定的位移，帶動坡腳左側水平段管道也產生較大的Mises應力和彎矩；坡頂附近管道存在一處Mises應力極值點，主要由于該處豎向位移載荷較大，管道在豎向載荷和重力的共同作用下產生了較大的拉伸應力和彎矩。管道坡腳Mises應力分布和坡頂附近Mises應力分布如圖6所示。

圖6 管道Mises應力分布極值點Fig.6 Extreme point of Mises stress distribution in pipeline

分別對縱向穿越滑坡體管道施加離散滑坡位移和固定位移，比較兩者有限元計算結果，前者管道最大Mises應力值比后者增大約50%，因此對管道施加離散滑坡位移進行受力分析的結論更加安全；雖然后者在坡腳彎管同樣也出現了應力極值點，但是在坡腳左側水平段、滑動面附近和坡頂附近管段均未出現應力極值點，主要由于對管道施加固定滑坡位移不能反映坡腳附近的滑動剪切作用及由此產生的影響，也不能反映坡頂豎向載荷比坡腳豎向載荷大對管道所產生的作用。

施加離散和固定滑坡位移條件下管道應力與彎矩的變化曲線如圖7所示。由圖7可看出：比較施加離散和固定位移條件下管道的彎矩分布，前者管道最大彎矩值比后者大1倍多，坡腳彎管兩側的彎矩最大，坡腳左側水平段、滑動面附近和坡頂附近也存在彎矩極值點，說明坡腳彎管變形最劇烈，這在一定程度上也緩解了滑坡位移對管道的作用；彎管左側水平段管道存在一處彎矩極值點，是由于坡腳彎管變形帶動該點變形引起的；滑動面附件存在彎矩極值點，說明管道在剪切作用下產生了明顯的局部彎曲變形；坡頂附近彎矩極值點是由于豎向滑坡位移分布不均使管道局部發生變形引起的。

圖7 離散和固定滑坡位移條件下管道應力與彎矩的變化曲線Fig.7 Variation curve of pipeline stress and bending moment under discrete and fixed landslide displacement

施加離散和固定滑坡位移條件下管道位移的變化曲線如圖8所示。由圖8可以看出：對穿越滑坡體管道施加離散滑坡位移時，管道在滑動面附近產生的水平位移最大，在坡腳彎管和坡頂彎管處產生的豎向位移較大，穿越滑坡體段管道總位移量相對于其兩側彎管位移量較小，沿坡頂固定土體段敷設的管道總位移量逐漸減小；對穿越滑坡體管道施加固定位移載荷時，管道產生的水平位移很小而豎向位移較大，穿越滑坡體段管道的豎向位移最大且基本不變，滑坡體兩端固定土體內敷設管道的位移則逐漸減小。比較施加固定和離散位移條件下管道產生位移的最大值，前者比后者增大約17%，因此對穿越滑坡體管道施加離散滑坡位移，雖然管道產生的位移相對較小，但是在局部土壤剪切作用下管道的最大應力卻較高；對穿越滑坡體管道施加固定滑坡位移，雖然管道產生的位移載荷較大，但是管道位移載荷分布較為均勻，不能體現局部土壤剪切作用下管道的受力狀況，管道最大應力反而較小。

圖8 離散和固定滑坡位移條件下管道位移的變化曲線Fig.8 Variation curve of pipeline displacement under discrete and fixed landslide displacement

內壓和溫差變化對縱向穿越滑坡體管道Mises應力影響如圖9所示。

由圖9可以看出：對穿越滑坡體管道施加離散滑坡位移，管道內壓p由0逐漸增大到10 MPa，坡腳左側水平段、坡腳彎管和滑動面附近管道Mises應力變化不大，坡頂附近Mises應力略有增加，增幅僅為7.3%，其余地段管道Mises應力均有明顯增加，所以提高管道內壓對坡腳左側水平段、坡腳彎管和滑動面附近管道影響不大，因此對于縱向穿越滑坡體管道，其最大Mises應力主要由滑坡位移載荷引起；當管道溫差ΔT由0增大到50 ℃時，由于管道沿線Mises應力分布基本不變，所以管道溫差變化對縱向穿越滑坡體管道應力分布狀況基本無影響。

圖9 內壓和溫差變化對縱向穿越滑坡體管道Mises應力的影響Fig.9 Effect of internal pressure and temperature difference variation on Mises stress of pipeline crossing sliding mass vertically

2.4 橫向穿越滑坡體管道受力分析

當油氣管道橫向穿越滑坡體時，由于滑坡體上、中、下三個部位水平位移、豎向位移和總位移均不相同，所以橫向穿越滑坡體上、中、下三個部位的管道受力狀況也各不相同。滑動面以上區域沿坡頂至坡腳方向，水平滑坡位移逐漸增大，豎向滑坡位移卻逐漸減小，總的滑坡位移量略有增加，如表2所示。

表2 滑坡體上、中、下部位的位移量 m

當管道橫向穿越滑坡體不同部位時，需要根據穿越部位實際滑坡位移量對管道進行加載，雖然滑坡體上部總位移量相對較小，但是管道Mises應力卻最大，滑坡體下部總位移量相對較小，但管道Mises應力卻較小，主要由于管道本身具有一定的自重，豎向載荷對其影響更為明顯，因此管道橫向穿越滑坡體時易從滑坡體下部通過，雖然滑坡體下部水平位移和總位移量均較大，但管道Mises應力卻較低。管道Mises應力最大值處于滑坡體與固定土體之間的交界面處，管道彎矩最大值處于交界面附近兩側，主要由于交界面兩側土體發生了較為明顯的位移錯動，導致穿越交界面處管道在土體剪切作用下產生較大的應力，交界面兩側管道同時也產生了較大的彎矩。橫向穿越不同部位滑坡體管道應力和彎矩的變化曲線如圖10所示。

圖10 橫向穿越不同部位滑坡體管道應力和彎矩的變化曲線Fig.10 Variation curves of stress and bending moment of pipeline crossing different parts of sliding mass horizontally

當管道橫向穿越滑坡體下部時，隨著管道內壓逐漸增大，管道Mises應力也逐漸加大，但是在滑坡體與固定土體交界面處，因滑坡位移產生的管道Mises應力占比卻逐漸減小，例如案例中管道內壓由0提高到10 MPa，因滑坡位移產生的管道Mises應力占管道總Mises應力比值由81%降低到9%，滑坡位移對管道應力的影響逐漸變得不明顯，所以當管道橫向穿越滑坡體時，高壓管道受到滑坡位移載荷的影響較小，而低壓管道受到的影響則較大。管道橫向穿越滑坡體下部，隨著管道溫差逐漸增大，管道Mises應力也逐漸加大，尤其在滑坡體與固定土體交界面處，因滑坡位移產生的管道Mises應力占比也逐漸增大，例如案例中溫差由0提高到50 ℃時，因滑坡位移產生的管道Mises應力占管道總Mises應力比值由2%增加到12%，所以當管道橫向穿越滑坡體時，溫差變化較大的管道受到滑坡位移載荷的影響更為明顯。內壓或溫差變化對橫向穿越滑坡體下部管道Mises應力影響如圖11所示。

圖11 內壓和溫差變化對橫向穿越滑坡體下部管道Mises應力的影響Fig.11 Effect of internal pressure and temperature difference variation on Mises stress of pipeline crossing lower part of sliding mass horizontally

3 結 論

(1)基于離散滑坡位移條件下油氣管道受力分析主要通過兩步來實現，第一步采用有限元法求解滑坡體各點位移，通過逐漸降低坡體強度來實現坡體的滑動過程，當管道縱向穿越滑坡體時需要沿管道敷設方向提取坡體各點的滑坡位移，當管道橫向穿越滑坡體時，需根據穿越坡體的不同部位來提取相應的滑動位移；第二步采用有限元建立管土作用模型，當管道縱向穿越滑坡體時，將沿管道敷設方向的離散滑坡位移載荷施加到對應位置的管道上，當管道橫向穿越滑坡體時，將根據不同穿越部位所提取的滑坡位移施加到管道上，從而實現離散滑坡位移載荷作用下管道受力分析。采用該種方法對管道進行受力分析，計算結論可使管道更加安全。

(2)對于縱向穿越滑坡體油氣管道施加離散滑坡位移，所得管道應力分布云圖能夠反映局部土體剪切作用對管道的影響，在坡腳左側水平段、坡腳彎管、滑動面附近和坡頂附近共計4處產生了應力極值點，管道受力狀況與滑坡體的滑動規律也更為相符，具有較好的適用性。縱向穿越滑坡體油氣管道Mises應力最大值和彎矩最大值均位于坡腳彎管附近，坡腳左側水平段和滑動面附近管道Mises應力值和彎矩值大致相當，坡頂附近管道Mises應力值和彎矩值進一步減小，因此坡腳彎管是縱向穿越滑坡體管道的防護重點。當內壓逐漸增大時，縱向穿越滑坡體管道Mises應力極值點變化不大，其余段管道Mises應力均隨著內壓的增大而增大；溫差變化對縱向穿越滑坡體油氣管道受力基本無影響。

(3)當油氣管道橫向穿越滑坡體時，從滑坡體下部位置通過時管道Mises應力值較低，主要由于滑坡體坡腳附近的水平位移較大，而坡頂附近的豎向位移較大，同時考慮管道自身重力的影響，滑坡豎向位移載荷比水平位移載荷對管道的影響更明顯。橫向穿越滑坡體油氣管道Mises應力極值點和彎矩極值點均位于固定土體和滑動土體的交界面處，呈對稱分布，因此穿越交界面處的管道為防護重點。當油氣管道橫向穿越滑坡體時，低壓管道相對于高壓管道受滑坡位移的影響更為明顯，而溫差變化較大的管道比溫差變化較小的管道更易受到滑坡位移的影響。