李 霖， 蘭才富， 李 虎

(1.國家管網(wǎng)集團(tuán)公司西南管道南寧輸油氣分公司， 南寧 530022； 2.西南石油大學(xué) 石油與天然氣工程學(xué)院，成都 610500； 3.四川省地質(zhì)工程勘察院集團(tuán)有限公司， 成都 610072)

在長輸油氣管道快速建設(shè)過程中會穿越眾多地質(zhì)災(zāi)害區(qū)域，其中也不可避免地穿越巖溶沉陷區(qū)[1]。據(jù)統(tǒng)計中國可溶巖分布面積約達(dá)到國土面積的1/3以上，是世界上巖溶最發(fā)育的國家之一，其中四川、云南、廣西等地分布最為廣泛[2]?？λ固貛r溶是易溶巖，主要隱藏在地下，由于其地下水發(fā)育，水溶性溶解巖石在水的作用下發(fā)生溶蝕現(xiàn)象，造成一定區(qū)域的空洞，在相應(yīng)的環(huán)境作用或外部荷載作用下造成上部巖土體發(fā)生沉降或坍塌，從而對埋地長輸油氣管道造成嚴(yán)重破壞，危機(jī)附近居民和環(huán)境[3-4]。

國外對喀斯特巖溶地區(qū)的研究始于19世紀(jì)末[5]，主要以對喀斯特地貌的評價為主。Magdalene等[6]和Gao等[7]運(yùn)用G1S技術(shù)對地質(zhì)條件、覆土厚度、地表水系等一系列因素進(jìn)行了量化分析并以此建立了相應(yīng)的評價模型。Zisman[8]運(yùn)用SPT參數(shù)且將地質(zhì)特征作為巖溶塌陷區(qū)危險評價因子，對佛羅里達(dá)州巖溶區(qū)進(jìn)行了評價。Katarina等[9]運(yùn)用加權(quán)回歸法、K-Function法等統(tǒng)計算法在空間上分析了人為因素和自然因素對巖溶塌陷區(qū)的影響。國內(nèi)學(xué)者針對塌陷區(qū)長輸油氣管道的研究主要為解析法和數(shù)值分析法。高惠瑛、馮啟民[10-11]在研究塌陷區(qū)使管道幾何大變形的基礎(chǔ)上將軸力與彎矩作用考慮到管道塑性狀態(tài)下，建立了塌陷區(qū)管道失效判據(jù)。張鵬等[12]在分析管道不同角度穿越沉陷區(qū)時得到了管道最大應(yīng)力的時變性規(guī)律以及管道最大應(yīng)力位置。游東潘[13]運(yùn)用ABAQUS有限元軟件對巖溶塌陷區(qū)管道的力學(xué)行為做了分析，提出了穿越巖溶地區(qū)管道建設(shè)的相關(guān)建議。王惠[14]運(yùn)用數(shù)值模擬分析方法對沉陷盆地的輸油管道力學(xué)行為進(jìn)行分析得出了相應(yīng)因素的敏感性順序。以上研究均未直接指出油氣管道在不同巖溶沉降長度下管道的最大位移和相應(yīng)狀態(tài)下的應(yīng)力分布狀態(tài)，也未指出巖溶地區(qū)塌陷后最大懸空長度和在一定沉陷長度下管道的預(yù)警位移。且對于長輸油氣管道穿越巖溶地帶時，由于不同的輸送內(nèi)壓可能會使得不同長度沉降下管道的最大位移不同。在一定沉降長度下，管道是否達(dá)到管道允許的最大位移將是初步判斷管道安全性的重要參數(shù)，且按90%屈服應(yīng)力所對應(yīng)的預(yù)警位移可作為管道安全運(yùn)營的預(yù)警參數(shù)，保證管道安全。因此，研究巖溶地區(qū)不同內(nèi)壓下管道最大位移和預(yù)警位移以及管道最大懸空長度對管道安全運(yùn)行具有重要作用。

針對以上問題，以西南某穿越巖溶地區(qū)埋地輸氣管道為例，運(yùn)用ABAQUS有限元軟件建立管-土三維模型，分析巖溶地區(qū)在不同輸送內(nèi)壓時不同沉陷長度下管道的最大位移和管道預(yù)警位移，得到其最危險位置，且考慮懸空狀態(tài)下管道最大懸空長度，以期為穿越巖溶地區(qū)管道的安全評價和及時治理提供參考。

1 穿越巖溶地區(qū)埋地管道變形特征

巖溶地區(qū)由于地下水的影響形成土洞，土洞上部土體沉陷使管道發(fā)生彎曲，甚至發(fā)生破壞，此時管道受到自身重力之外還受到上部土體重力影響。根據(jù)現(xiàn)有文獻(xiàn)對巖溶塌陷區(qū)管道力學(xué)行為的影響，將管道看作是受均布荷載的Winkler地基梁模型[15-16]，在此模型基礎(chǔ)上計算得到管道位移為

y=

(1)

式中：參數(shù)λ與彎矩M的表達(dá)式為

(2)

(3)

式中：y為管道撓度,m;W為管道最大沉降值,m;M為管道沉降邊緣點(diǎn)彎矩,N·m;l為最大沉降點(diǎn)到邊緣點(diǎn)的距離,m;EI為管道的彎曲剛度;x為管道上一點(diǎn)到沉降邊緣點(diǎn)的距離,m;k為彈性地基系數(shù)。

通過管道撓度方程和邊界點(diǎn)彎矩可以得到受彎曲影響下管道上任意一點(diǎn)應(yīng)力值，即

(4)

根據(jù)式(1)～式(4)可在一定精度上計算管道的力學(xué)響應(yīng)。

2 巖溶沉陷區(qū)數(shù)值分析模型

2.1 有限元模型

本模型以西南某穿越巖溶地區(qū)輸氣管道為例，該線所用輸送管道為X80管材，經(jīng)勘探管道平均埋深2 m，外徑1 016 mm，壁厚15.3 mm，設(shè)計運(yùn)輸內(nèi)壓為10 MPa，管材屈服應(yīng)力為555 MPa。采用ABAQUS有限元軟件進(jìn)行模擬分析，土體采用摩爾-庫侖模型，材料參數(shù)如表1所示，查閱文獻(xiàn)[17]得到管道應(yīng)力應(yīng)變數(shù)據(jù)作為本次模擬的管道力學(xué)參數(shù)，如圖1所示。計算模擬長度為10～100 m時管道不同內(nèi)壓下的最大位移狀態(tài)和在管道發(fā)生懸空時不同內(nèi)壓下的極限懸空長度，將90%的屈服應(yīng)力作為管道的預(yù)警應(yīng)力[18]，分析預(yù)警應(yīng)力時不同沉陷長度下管道的預(yù)警位移。

表1 材料參數(shù)

圖1 x80管材應(yīng)力-應(yīng)變曲線

圖2 巖溶沉陷模型

查閱相關(guān)資料，由于管道受影響的非沉陷區(qū)域長度是沉降區(qū)域長度的0.49倍[19]，且模擬沉陷區(qū)長度為10～100 m，故選取模型長度為200 m，土體橫截面簡化為5 m×5 m。將模型走向和橫向部分分別施加對稱約束，將模型底部未沉陷區(qū)域設(shè)置為固定約束，如圖2所示。管-土之間采用罰函數(shù)約束其接觸行為，管-土間摩擦系數(shù)取0.5。由于在外界荷載作用下短距離的沉陷區(qū)域仍然可能使得管道發(fā)生破壞，且?guī)r溶沉降為漏斗式，故在沉降區(qū)域土體上施加二次拋物線位移函數(shù)以模擬在自然沉降或受外界荷載作用下管道達(dá)到極限位移，其函數(shù)為

(5)

式中：y為沉陷區(qū)各點(diǎn)位移值,m;x為沉陷區(qū)內(nèi)一點(diǎn)到沉陷邊緣的距離,m;c為沉陷下沉值,m;l為沉陷區(qū)長度,m。

2.2 模型驗(yàn)證

本次建模所采取的邊界約束條件與文獻(xiàn)[19]相似，故僅對模型的網(wǎng)格敏感性進(jìn)行分析，以期能使模擬結(jié)果更加精確。以沉陷區(qū)域長度10 m為例，將非沉陷區(qū)網(wǎng)格尺寸取為2 m，沉陷區(qū)網(wǎng)格尺寸取2、1.5、1、0.5 m作為敏感性分析尺寸，劃分網(wǎng)格總體數(shù)量分別為16 796、18 116、20 656、28 376個，分別施加最大位移為0.2 m的位移二次函數(shù)，結(jié)果如圖3所示。

圖3 模型網(wǎng)格驗(yàn)證

由圖3可知，網(wǎng)格尺寸為0.5、1、1.5、2 m時管道最大應(yīng)力分別為571.04、570.46、564.34、563.18 MPa，可知當(dāng)網(wǎng)格尺寸小于1 m時，其管道最大應(yīng)力相差極小，相對誤差值在5%之內(nèi)，認(rèn)為結(jié)果可行可信，故取網(wǎng)格尺寸為1 m。

3 沉陷區(qū)管道響應(yīng)分析

根據(jù)所建立的模型，由于需要計算到管道達(dá)到屈服的最大位移，故在沉陷區(qū)施加從0.2到1.2 m，增量為0.05 m的位移，分別計算無壓、5 MPa、10 MPa時管道在不同沉降長度下達(dá)到管道屈服應(yīng)力和預(yù)警應(yīng)力時的力學(xué)響應(yīng)分析，并得到管道屈服和預(yù)警時的最大位移和預(yù)警位移。

3.1 無壓時管道響應(yīng)分析

為分析管道在無壓時受沉陷作用的力學(xué)行為變化，得到管道最大位移和預(yù)警位移以及最危險位置，由有限元計算得到沉降長度為10 m的土體位移云圖如圖4所示，管道位移云圖如圖5所示，各長度下管道達(dá)到屈服時應(yīng)力如圖6所示，管道最大位移如圖7所示。其達(dá)到屈服時最大位移、管道預(yù)警位移及危險位置見表2。

圖4 無壓時沉降10 m土體位移云圖

圖5 無壓時沉降10 m管道位移云圖

圖6 無壓時不同沉陷長度下管道應(yīng)力曲線

圖7 無壓時不同沉陷長度下管道位移曲線

由圖6可知，當(dāng)沉降長度小于等于25 m時，管道在跨中位置達(dá)到屈服，此時兩端最大應(yīng)力為393.75 MPa，這是由于在短距離沉降時管道應(yīng)力主要來自于受到沉降的彎曲作用，而在兩端位置管道應(yīng)力則隨著沉降長度的增長而不斷增加，兩端受到的剪力作用不斷加強(qiáng)。當(dāng)沉降長度大于等于30 m時，管道則在兩端位置率先達(dá)到屈服強(qiáng)度，沉陷長度30 m時跨中最大應(yīng)力為395.86 MPa，管道主要收到兩端土體位移的剪切作用，跨中位置管道應(yīng)力則隨著沉陷長度的增長而不斷減小。

由圖7可知，當(dāng)管道達(dá)到屈服時，其最小極限位移為沉降長度10 m時達(dá)到112.05 mm，最大極限位移為沉降長度100 m時達(dá)到1 035.25 mm，且隨著沉降長度的增加管道能承受的最大位移也在不斷增加，且基本成線性增長。

表2 無壓時沉降管道危險位移

由表2可知，根據(jù)預(yù)警位移與最大位移對比，沉陷長度為10 m時兩者相差較小，預(yù)警位移為107.81 mm，達(dá)到預(yù)警級別，差值僅有4.24 mm，隨著長度的增加差值也逐漸增加，長度達(dá)到100 m時達(dá)到以預(yù)警級別的預(yù)警位移為906.92 mm，差值達(dá)到128.33 mm。

3.2 5 MPa內(nèi)壓下管道響應(yīng)分析

在5 MPa的運(yùn)輸內(nèi)壓下管道的力學(xué)行為變化以及預(yù)警位移結(jié)果如下，不同長度下管道達(dá)到屈服時應(yīng)力如圖8所示，管道屈服時達(dá)到最大位移如圖9所示。其達(dá)到屈服時最大位移、管道預(yù)警位移及危險位置見表3。

圖8 5 MPa時不同沉陷長度下管道應(yīng)力曲線

圖9 5 MPa時不同沉陷長度下管道位移曲線

表3 5 MPa時沉陷管道危險位移

由圖8可知，在5 MPa內(nèi)壓下管道沉陷長度小于等于20 m時，管道在沉陷段跨中率先達(dá)到555 MPa屈服強(qiáng)度，此時沉陷20 m時管道兩端應(yīng)力為405.42 MPa，管道沿線應(yīng)力成山字形分布。當(dāng)沉陷長度大于等于25 m時，管道則在兩端位置達(dá)到屈服強(qiáng)度，沉陷25 m時管道跨中最大應(yīng)力為458.54 MPa，且隨著沉陷長度的增加其跨中位置最大應(yīng)力逐漸減小。

由圖9可知，管道沿線位移趨勢與無壓時管道位移趨勢相同，其最小極限位移為沉陷長度10 m時達(dá)到103.97 mm，最大極限位移為沉陷長度100 m時達(dá)到958.82 mm。

由表3可知，管道沉陷長度在20 m以下時危險位置位移跨中，沉陷長度超過25 m則危險位置在兩端，比無壓時危險位置最早出現(xiàn)在兩端的長度短，且預(yù)警位移與屈服位移均比無壓時小，預(yù)警位移與最大位移之差隨著沉陷長度的增大逐漸增大，最小為沉陷10 m時，預(yù)警位移為98.11 mm，相差5.86 mm。當(dāng)沉陷長度為100 m時最大位移比無壓時少76.5mm，減少了7.39%。

3.3 10 MPa內(nèi)壓下管道響應(yīng)分析

在設(shè)計內(nèi)壓下管道的力學(xué)行為變化以及管道預(yù)警位移變化如下，不同長度下管道達(dá)到屈服時應(yīng)力如圖10所示，管道達(dá)到屈服時的位移如圖11所示。其達(dá)到屈服時最大位移和預(yù)警位移及危險位置見表4。

圖10 10 MPa時不同沉陷長度下管道應(yīng)力曲線

圖11 10 MPa時不同沉陷長度下管道位移曲線

由圖10可知，在10 MPa設(shè)計內(nèi)壓下沉陷長度小于等于15 m時，管道在沉陷段跨中位置率先達(dá)到屈服，此時管道跨端位置應(yīng)力為394.95 MPa，管道沿線應(yīng)力成山字形分布。當(dāng)沉陷長度大于等于20 m時，管道則在兩端位置達(dá)到屈服強(qiáng)度，此時沉陷25m時管道跨中位置應(yīng)力為485.21 MPa，其應(yīng)力變化趨勢與無壓和5 MPa工況相同。

由圖11可知，管道沿線位移趨勢與無壓時管道位移趨勢相同，其最小極限位移為沉陷長度10 m時達(dá)到95.99 mm，最大極限位移為沉陷長度100 m時達(dá)到826.78 mm。

表4 10 MPa時沉陷管道危險位移

由表4可知，管道沉陷100 m時預(yù)警位移為684.42 mm，與最大位移相差，危險位置在兩端，沉陷長度小于15 m時危險位置位于跨中，超過20 m則危險位置位于兩端處。對比表2、3、4可知，管道內(nèi)壓越大管道沉陷危險位移越小，且管道內(nèi)壓越大，危險位置出現(xiàn)在兩端時的沉陷長度越短。故在發(fā)生溶巖沉陷時及時停輸斷氣不僅可以防止泄露也可以加大管道的危險位移，在一定程度上防止管道發(fā)生二次危害，且根據(jù)90%屈服應(yīng)力所得到的管道預(yù)警位移可在一定程度上為管道安全預(yù)警做出判斷，提高管道安全性。

4 懸空時管道響應(yīng)分析

管道在穿越巖溶地帶時，若巖溶空洞發(fā)生突然性的塌陷，則使得管道會發(fā)生突發(fā)性懸空，管道在不同內(nèi)壓下最大的懸空長度則是對管道安全性進(jìn)行初步判斷的重要參數(shù)。因此，將分別計算不同內(nèi)壓下管道能達(dá)到的最大懸空長度，根據(jù)相關(guān)文獻(xiàn)[20-21]建立相關(guān)的有限元模型，分析在懸空長度從50 m到600 m時管道的力學(xué)響應(yīng)變化，不同懸空長度下管道最大應(yīng)力-應(yīng)變變化曲線如圖12所示，最大位移變化曲線如圖13所示。

圖12 管道最大應(yīng)力-應(yīng)變曲線

圖13 管道最大位移曲線

由圖12可知，管道在內(nèi)壓為10、5 MPa和無壓狀態(tài)下達(dá)到屈服時的最大懸空長度分別約為280、340、380 m，而按照0.9倍屈服極限的預(yù)警位移則最大懸空長度分別約為190、250、270 m，可知管道最大應(yīng)力隨著懸空長度的增加先迅速增加，達(dá)到屈服應(yīng)力后開始區(qū)域平穩(wěn)，而管道應(yīng)變則隨著懸空長度的增加而增加，且基本成線性增長趨勢。由圖13可知，管道在不同內(nèi)壓下最大位移隨著懸空長度的增加差值逐漸增大，且隨著懸空長度的增加管道最大位移成線性增長。

5 結(jié)論

通過有限元數(shù)值分析，分析了在不同內(nèi)壓下管道不同沉陷長度時的最大位移和預(yù)警位移以及不同內(nèi)壓下管道的最大懸空距離，得出以下結(jié)論：

1)當(dāng)沉陷長度較小，管道跨中達(dá)到屈服時，隨著沉陷長度的增加管道兩端最大應(yīng)力逐漸增加，當(dāng)管道兩端達(dá)到屈服時，管道跨中位置最大應(yīng)力隨著沉陷長度的增加而不斷減小，且管道應(yīng)力呈“山”字形分布。

2)管道內(nèi)壓越大，危險位置出現(xiàn)在兩端時的沉陷長度越短，不同沉陷長度下管道最大位移和預(yù)警位移越小。

3)管道發(fā)生懸空時，內(nèi)壓越大其管道臨界懸空長度越短，在設(shè)計內(nèi)壓下管道最大懸空長度約為280 m。