非均勻沉降下沿海天然氣埋地管道安全運(yùn)行分析

摘 要 埋地天然氣管道運(yùn)行安全受地基非均勻沉降現(xiàn)象影響嚴(yán)重。為確保管道安全運(yùn)行，通過對沿海地區(qū)埋地天然氣管道進(jìn)行應(yīng)力監(jiān)測，獲取管道在地基非均勻沉降影響下的環(huán)向、軸向附加應(yīng)力，并結(jié)合管道工作內(nèi)壓進(jìn)行管道應(yīng)力校核。基于ANSYS軟件建立埋地天然氣管道有限元模型，對比分析應(yīng)力監(jiān)測實(shí)際值與數(shù)值模擬結(jié)果，驗(yàn)證了應(yīng)力監(jiān)測數(shù)據(jù)的有效性和有限元模型的可行性。進(jìn)一步通過有限元數(shù)值模擬探究管道內(nèi)徑、壁厚和埋土參數(shù)對管道應(yīng)力的影響。數(shù)值模擬結(jié)果表明：在管道敷設(shè)階段，通過減小管道內(nèi)徑、適當(dāng)增加管道壁厚和選擇較硬的埋土可以減小地基非均勻沉降對管道的影響。研究成果可為受地基非均勻沉降影響的埋地天然氣管道運(yùn)行安全提供參考。

關(guān)鍵詞 非均勻沉降 埋地天然氣管道 應(yīng)力分析 安全運(yùn)行

中圖分類號 TE973" "文獻(xiàn)標(biāo)識碼 A" "文章編號 0254?6094（2023）01?0076?07

我國沿海地區(qū)的天然氣站場多建在海濱、湖濱及河流沿岸等地，地質(zhì)條件多為強(qiáng)度低、壓縮性較高的軟土地基，因其流變性和不均勻性易產(chǎn)生非均勻沉降，造成場地運(yùn)動。根據(jù)歐洲天然氣管道事故數(shù)據(jù)組織（EGIG）的統(tǒng)計分析，由于場地運(yùn)動造成的埋地天然氣管道事故占管道事故總數(shù)的16%[1]。

不少學(xué)者通過有限元數(shù)值模擬研究埋地管道的場地運(yùn)動，TRIFONOV O V和CHERNIY V P研究了活動斷層作用下埋地管道的應(yīng)力應(yīng)變情況[2];WANG Y等研究隧道開挖造成的土體隆起及沉降對管道的影響[3];趙新威等通過建立土-埋地管道在動力接觸條件下的相互作用模型，探究管道在地震載荷作用下的失效過程[4];馬樹峰等研究高陡邊坡輸氣管道受坡高、管徑、內(nèi)壓和壁厚對管道應(yīng)力和變形的影響[5];張羽翀等提出坡體蠕滑作用下管道初始應(yīng)力的簡化思路[6];呂彥偉總結(jié)了管道在土體非均勻沉降作用下的受力變形特征，預(yù)測管道破壞的臨界條件[7]。還有學(xué)者通過監(jiān)測技術(shù)對埋地管道進(jìn)行研究，劉玉卿等利用振弦式應(yīng)變計監(jiān)測管道軸向應(yīng)力，建立管道軸向檢測應(yīng)力的預(yù)警模型[8];WANG F等針對多年凍土地區(qū)長輸管道建立了一套全面的現(xiàn)場監(jiān)測系統(tǒng)，闡明不同地質(zhì)條件下管道與凍土退化之間的相互作用[9]。以上學(xué)者分別通過有限元數(shù)值分析和應(yīng)力監(jiān)測對埋地管道地場運(yùn)動進(jìn)行研究，對地基非均勻沉降的地場運(yùn)動有指導(dǎo)意義，但有限元的結(jié)果難以在實(shí)際工程中得到驗(yàn)證，應(yīng)力監(jiān)測結(jié)果也存在一定的誤差，無法校核其準(zhǔn)確性。

筆者選取華南地區(qū)某分輸站為研究對象，該分輸站處于沿海地區(qū)，站場建成至今非均勻沉降現(xiàn)象嚴(yán)重，通過對管道進(jìn)行應(yīng)力監(jiān)測，判斷非均勻沉降對管道安全運(yùn)行的影響，基于有限元數(shù)值模擬校核管道應(yīng)力，并探究地基非均勻沉降作用下管道內(nèi)徑、壁厚和埋土性質(zhì)對管道應(yīng)力的影響。

1 工程概況

筆者對某分輸站出站管道展開研究，站內(nèi)地基主要由人工填土、粉質(zhì)黏土、砂土構(gòu)成，站外地基主要由淤泥質(zhì)土、淤泥、黏土構(gòu)成，站內(nèi)出現(xiàn)多處沉降，站外較為平整。管道采用專用材質(zhì)L450，內(nèi)徑1 000 mm，壁厚22 mm，最大埋深4 450 mm，通過對管道站場周期性（月）測試記錄，管道運(yùn)行過程中工作壓力從2021年1月的6.75 MPa逐步調(diào)升至2021年10月的7.61 MPa，出站管道處地基沉降從21.2 mm增至78.6 mm，出站管道實(shí)際工況記錄見表1，這種非均勻沉降現(xiàn)象會導(dǎo)致出站的埋地管道因受載不均而發(fā)生變形，管道應(yīng)力隨之發(fā)生變化，應(yīng)力的變化實(shí)際是埋地管道與管周土壤相互作用的結(jié)果[10]，應(yīng)力過大會導(dǎo)致管道開裂，發(fā)生泄漏，危害天然氣分輸站運(yùn)行安全，影響下游天然氣供應(yīng)。因此對分輸站管道進(jìn)行應(yīng)力監(jiān)測與分析對管道安全意義重大。

2 管道應(yīng)力監(jiān)測

為了避免埋地管道受地基非均勻沉降影響過大而存在安全隱患，需要掌握埋地管道安全狀態(tài)，通過對管道開展應(yīng)力監(jiān)測，可反映管道受非均勻沉降的應(yīng)力水平，從而評判管道的安全并采取相應(yīng)對策措施。應(yīng)力監(jiān)測是通過多種設(shè)備和儀器在管道施工、管道安裝和管道運(yùn)行過程對管道應(yīng)力變化實(shí)現(xiàn)可視化處理的方法[11]。筆者使用電阻應(yīng)變片法對分輸站的進(jìn)站管道進(jìn)行周期性（月）應(yīng)力監(jiān)測。

2.1 監(jiān)測點(diǎn)設(shè)置

管道應(yīng)力監(jiān)測中選擇管道應(yīng)力較大且集中的部位，結(jié)合管道走向、結(jié)構(gòu)及焊接位置等因素確定測點(diǎn)。根據(jù)分輸站管道情況，在進(jìn)站埋地管道上彎頭處設(shè)置了4個測點(diǎn)，如圖1、2所示，每個測點(diǎn)沿管道的環(huán)向和軸向設(shè)置3組應(yīng)變片進(jìn)行管道應(yīng)變監(jiān)測，應(yīng)變片型號為BE120?3BA兩向直角應(yīng)變片，通過電纜信號線將東華DH3815靜態(tài)應(yīng)變測試儀與應(yīng)變片連接，以2 Hz的頻率采集管道因地基非均勻沉降而產(chǎn)生的應(yīng)變數(shù)據(jù)，為方便后續(xù)數(shù)據(jù)分析，取對應(yīng)測點(diǎn)3組應(yīng)變片采集的數(shù)據(jù)的平均值為測試值。

2.2 監(jiān)測數(shù)據(jù)處理

因管道有內(nèi)壓，在連接好監(jiān)測系統(tǒng)后，需要進(jìn)行系統(tǒng)調(diào)零平衡，系統(tǒng)平衡后對管道進(jìn)行周期性（月）監(jiān)測，記錄監(jiān)測管道的應(yīng)變值，通過應(yīng)力-應(yīng)變關(guān)系公式計算系統(tǒng)監(jiān)測產(chǎn)生的環(huán)向應(yīng)力σ和軸向應(yīng)力σ：

σ=（ε+με） （1）

σ=（ε+με）（2）

式中 E——管道用材的彈性模量，200 GPa;

ε——所測環(huán)向應(yīng)變;

ε——所測軸向應(yīng)變;

μ——管道用材的泊松比，0.30。

監(jiān)測過程中，管道工作內(nèi)壓也有升降，內(nèi)壓變動也會導(dǎo)致管道應(yīng)力變化，管道工作內(nèi)壓產(chǎn)生的基本應(yīng)力主要有環(huán)向應(yīng)力σ和軸向應(yīng)力σ：

σ=（3）

σ=（4）

式中 D——管道外徑，mm;

p——管道工作內(nèi)壓，MPa;

δ——管道壁厚，mm。

此時需要減去壓力變化對應(yīng)的應(yīng)力值才是地基非均勻沉降造成的管道附加應(yīng)力。

通過10次應(yīng)變監(jiān)測的數(shù)據(jù)，計算管道4個測點(diǎn)因地基非均勻沉降導(dǎo)致的附加環(huán)、軸向應(yīng)力大小，應(yīng)力變化曲線結(jié)果圖3所示。

受地基非均勻沉降影響，4個測點(diǎn)變化趨勢大致相同，但呈現(xiàn)應(yīng)力狀態(tài)有所不同：

a. 1#測點(diǎn)位于彎頭上端，受地基非沉降的作用，直管段受土壓力的影響，對彎頭上方造成拉伸，產(chǎn)生拉應(yīng)力（應(yīng)力數(shù)值為正），最大應(yīng)力為環(huán)向拉應(yīng)力29.14 MPa;

b. 3#測點(diǎn)位于彎頭下端，與1#測點(diǎn)上下對稱，該測點(diǎn)區(qū)域受下方地基支撐力和直管段的擠壓產(chǎn)生壓應(yīng)力（應(yīng)力數(shù)值為負(fù)），最大應(yīng)力為環(huán)向壓應(yīng)力22.63 MPa;

c. 2#、4#測點(diǎn)位于彎頭左右兩側(cè)，該測點(diǎn)區(qū)域受側(cè)向土壓力的影響，產(chǎn)生的應(yīng)力狀態(tài)在拉-壓應(yīng)力波動，但應(yīng)力水平較小，后逐漸轉(zhuǎn)為拉應(yīng)力，最大應(yīng)力為2#測點(diǎn)的環(huán)向拉應(yīng)力22.69 MPa。

測試結(jié)果表明非均勻沉降下，埋地管道受載的復(fù)雜性、應(yīng)力狀態(tài)的多變性。彎頭是管道的薄弱部位，受地基非均勻沉降造成的影響最大，地基非均勻沉降容易使彎頭前后出現(xiàn)變形，造成應(yīng)力過大。

為判斷埋地管道在地基非均勻沉降下的安全狀態(tài)，需要校核管道應(yīng)力。根據(jù)氣體輸送和分配管道系統(tǒng)標(biāo)準(zhǔn)[12]要求，由壓力、活載荷和靜載荷以及風(fēng)和地震等耦合荷載在管道中產(chǎn)生的軸向應(yīng)力之和不超過管道材料最小屈服強(qiáng)度的90%，工程案例中管道載荷主要由管道內(nèi)壓的永久載荷和地基沉降附加在管道上的載荷，管道軸向應(yīng)力校核：

σ=σ+σ≤0.9σ（5）

式中 σ——附加載荷產(chǎn)生的管道軸向應(yīng)力，MPa;

σ——管道材料最小屈服強(qiáng)度，MPa。

管道軸向應(yīng)力校核結(jié)果如圖4所示，4個測點(diǎn)中，最大軸向應(yīng)力為81.67 MPa，小于材料最小屈服強(qiáng)度的90%（450×0.9=405 MPa），故判定管道屬于安全運(yùn)行狀態(tài)。

3 管道應(yīng)力有限元分析

應(yīng)力監(jiān)測過程中可能會遇到應(yīng)變片失效、設(shè)備數(shù)據(jù)采集時外部環(huán)境干擾的影響，為了驗(yàn)證應(yīng)力監(jiān)測數(shù)據(jù)的準(zhǔn)確性，通過有限元數(shù)值模擬，分析地基非均勻沉降下埋地管道的應(yīng)力，并通過有限元分析研究影響管道應(yīng)力的因素。

3.1 有限元模型

建立模型與實(shí)物1：1的管道幾何模型，如圖5、6所示，幾何模型參數(shù)見表2、3。

對模型施加約束條件，整體考慮重力載荷，對埋土底部和四周施加固定約束，埋土頂部中部設(shè)置豎直向下位移（圖7中綠色部分），表征非均勻沉降，位移參數(shù)與實(shí)際所測數(shù)據(jù)一致，管道內(nèi)壓與實(shí)測工況一致，管道模型依據(jù)Ramberg?Osgood本構(gòu)模型，地基埋土本構(gòu)模型依據(jù)Mohr?Coulomb準(zhǔn)則[13]。管道和土體均采用Solid186四面體單元，建立管道與土體之間的面接觸，管道為目標(biāo)幾何體，埋土為接觸幾何體，將模型分別劃分為139 214個單元和248 029個節(jié)點(diǎn)，如圖7所示。

3.2 有限元數(shù)值結(jié)果對比

通過ANSYS有限元軟件對管土模型進(jìn)行數(shù)值模擬，以第10次應(yīng)力監(jiān)測實(shí)際工況為例，有限元模擬管道應(yīng)力云圖如圖8所示，實(shí)測彎頭段模擬應(yīng)力（Von?Mises等效應(yīng)力）最大值189.9 MPa、最小值128.67 MPa、均值為156.24 MPa，根據(jù)GB/T 9711—2017《石油天然氣工業(yè)管線輸送系統(tǒng)用鋼管》[14]，該最大應(yīng)力在安全范圍以內(nèi)。

根據(jù)Von?Mises屈服準(zhǔn)則[15]，計算管道Von?Mises等效應(yīng)力σ：

σ=（6）

其中σ、σ、σ表示3個方向（軸、環(huán)、徑）主應(yīng)力。

由于管道滿足薄壁特征，徑向應(yīng)力遠(yuǎn)小于其他兩個方向的應(yīng)力，可忽略。將應(yīng)力監(jiān)測測出的管道附加環(huán)軸應(yīng)力與管道內(nèi)壓導(dǎo)致的環(huán)軸應(yīng)力組合計算Von?Mises等效應(yīng)力，計算出管道應(yīng)力結(jié)果與有限元模擬應(yīng)力結(jié)果對比如圖9所示。

實(shí)測應(yīng)力變化曲線與模擬應(yīng)力變化曲線，在趨勢上總體一致，模擬結(jié)果與實(shí)測結(jié)果相對誤差最大6.24%，表明應(yīng)力監(jiān)測可以反映管道應(yīng)力狀態(tài)的可靠性，同時該管土模型也能準(zhǔn)確模擬管道在地基非均勻沉降作用下的應(yīng)力狀態(tài)。

3.3 管道應(yīng)力影響因素

通過應(yīng)力監(jiān)測和有限元分析可知地基非均勻沉降會對管道造成影響，會使管道產(chǎn)生不可忽視的附加應(yīng)力。為了減小地基非均勻沉降的影響，從管道本身和地基埋土兩方面考慮，通過有限元模型探究管道受地基非均勻沉降作用的影響因素。參考近階段管道實(shí)際工況數(shù)據(jù)，使模型接近實(shí)際情況，設(shè)置管道內(nèi)壓為7.11 MPa和地基局部沉降70 mm。選擇管道內(nèi)徑、管道壁厚、埋土力學(xué)參數(shù)3個因素進(jìn)行研究。

3.3.1 管道內(nèi)徑

根據(jù)GB/T 9711—2017對管道規(guī)格的要求，選擇內(nèi)徑960、1 000、1 100、1 200、1 300 mm5種規(guī)格的管道進(jìn)行管道應(yīng)力數(shù)值模擬，由圖10可見，管道內(nèi)徑越大，地基非均勻沉降造成的影響越大，管道應(yīng)力越大，在選用大口徑管道輸送時，應(yīng)加強(qiáng)對埋地管道地基的防護(hù)。

3.3.2 管道壁厚

管道壁厚是影響管道剛度和強(qiáng)度的主要因素[16]，依據(jù)GB/T 9711—2017對管道規(guī)格的要求，選擇管道壁厚17、22、27、32、37 mm5種規(guī)格的管道進(jìn)行管道應(yīng)力數(shù)值模擬，由圖11可見，隨著管道壁厚增加，地基非均勻沉降造成的影響越小，管道應(yīng)力越小，影響趨勢隨著管道壁厚增加有所放緩。在實(shí)際埋地管道建設(shè)中，可適當(dāng)增加管道壁厚可有效緩解地基非均勻沉降的影響，但管道壁厚過厚，減緩效果并不明顯同時加大建設(shè)成本。

3.3.3 埋土參數(shù)

埋土的力學(xué)性能參數(shù)包括：彈性模量、內(nèi)摩擦角和黏聚力，都是用來評價土的承載力、穩(wěn)定性的參數(shù)。有學(xué)者在做埋土參數(shù)對管道應(yīng)力影響研究時，往往都只考慮土的彈性模量，或者將土的各參數(shù)分開研究，這樣缺少實(shí)際指導(dǎo)意義。根據(jù)《工程地質(zhì)手冊（第四版）》[17]對土的力學(xué)性能參數(shù)的經(jīng)驗(yàn)數(shù)據(jù)，選擇5種粉質(zhì)黏土研究埋土參數(shù)對管道應(yīng)力的影響，埋土參數(shù)見表4。

由圖12可見，埋土的抗剪強(qiáng)度越大、抗壓縮強(qiáng)度越大，地基穩(wěn)定性越高，非均勻沉降對管道的影響越小。在實(shí)際工程中，選擇硬質(zhì)的埋土能夠較好緩解非均勻沉降的影響。

4 結(jié)論

4.1 通過監(jiān)測分輸站出站管道應(yīng)力可知，地基非均勻沉降會對管道產(chǎn)生影響，造成管道產(chǎn)生的拉、壓應(yīng)力，出站管道埋地彎頭處4個測點(diǎn)中最大拉應(yīng)力為彎頭上端的環(huán)向拉應(yīng)力29.14 MPa，最大壓應(yīng)力為彎頭下端的環(huán)向壓應(yīng)力22.63 MPa。

4.2 通過管道軸向應(yīng)力校核結(jié)果，最大管道軸向應(yīng)力為81.67 MPa，滿足設(shè)計規(guī)范中的強(qiáng)度要求，地基非均勻沉降對管道有一定影響，現(xiàn)階段不會對管道安全運(yùn)行造成危害，但仍需加強(qiáng)監(jiān)測管理。

4.3 通過進(jìn)行非均勻沉降埋地管道有限元分析，管道應(yīng)力數(shù)值模擬結(jié)果與管道應(yīng)力監(jiān)測結(jié)果對比，相對誤差最大6.24%，驗(yàn)證了應(yīng)力監(jiān)測結(jié)果的準(zhǔn)確性和可信度，針對存在地基非均勻沉降的埋地管道，實(shí)施管道應(yīng)力監(jiān)測能有效反映管道應(yīng)力狀態(tài)，可用于評估管道安全，同時建立的有限元管土模型可準(zhǔn)確反映非均勻沉降作用下管道的應(yīng)力狀態(tài)，可為管道運(yùn)行管理提供參考依據(jù)。

4.4 通過有限元分析研究埋地管道應(yīng)力影響因素，在沿海地區(qū)敷設(shè)輸氣管道時，為了減小地基非均勻沉降作用對埋地管道的影響，可采取減小管道內(nèi)徑、適當(dāng)增大管道壁厚、填鋪較硬的埋土等措施。

參 考 文 獻(xiàn)

[1]" "European Gas Pipeline Incident Data Group.11th

EGIG report [R/OL].2020?12?17[2022?02?20].http：//www.egig.eu/.

[2]" TRIFONOV O V，CHERNIY V P.A semi?analytical approach to a nonlinear stress?strain analysis of buried steel pipelines crossing active faults[J].Soil Dynamics and Earthquake Engineering，2010，30（11）：1298-1308.

[3]" WANG Y，WANG Q，ZHANG K Y.An Analytical Model for Pipe?Soil?Tunneling Interaction[J].Procedia Engineering，2011，14：3127-3135.

[4]" "趙新威，曾祥國，姚安林，等.地震載荷作用下埋地輸氣管道的數(shù)值模擬[J].石油機(jī)械，2014，42（3）：104-109.

[5]" "馬樹鋒，王勃，蘇海平，等.高陡邊坡輸氣管道的受力特性分析[J].石油機(jī)械，2020，48（6）：143-148.

[6]" "張羽翀，李龍，黃鑫，等.凍土區(qū)坡體蠕滑作用下管道應(yīng)力數(shù)值計算[J].油氣儲運(yùn)，2021，40（2）：185-191.

[7]" "呂彥偉.土體非均勻沉降作用下埋地含缺陷管道破壞規(guī)律研究[D].淮南：安徽理工大學(xué)，2021.

[8]" "劉玉卿，余志峰，佟雷，等.基于軸向應(yīng)力監(jiān)測數(shù)據(jù)的管道應(yīng)力狀態(tài)預(yù)警模型[J].石油機(jī)械，2018，46（6）：105-109.

[9]" "WANG F，LI G Y，MA W，et al.Pipelinepermafrost interaction monitoring system along the ChinaRussia crude oil pipeline[J].Engineering Geology，2019，254：113-125.

[10]" "吳敬龍.非均勻地基條件下地下綜合管廊受力特性研究[D].長沙：中南林業(yè)科技大學(xué)，2019.

[11]" HU X L，TAN F L，TANG H M，et al.In?situ monitoring platform and preliminary analysis of monitoring data of Majiagou landslide with stabilizing piles[J].Engineering Geology，2017，228：323-336.

[12] ASME.Gas transmission and distribution piping systems：ASME B31.8?2018[S].New York：American Society of Machanical Engineering Press，2018.

[13]" " 范磊，王衛(wèi)軍，袁超，等.基于Mohr?Coulomb與Drucker?Prager強(qiáng)度準(zhǔn)則的巷道圍巖數(shù)值模擬對比分析[J].礦業(yè)研究與開發(fā)，2019，39（12）：28-34.

[14]" "中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局，中國國家標(biāo)準(zhǔn)化管理委員會.石油天然氣工業(yè)管線輸送系統(tǒng)用鋼管：GB/T 9711—2017[S].北京：中國標(biāo)準(zhǔn)出版社，2017.

[15]" "林俊丞.天然氣分輸站不均勻沉降變形管線應(yīng)力計算[D].青島：中國石油大學(xué)（華東），2018.

[16]" "康遜.天然氣管道在非均勻沉降狀態(tài)下的應(yīng)力測試與分析[D].廣州：華南理工大學(xué)，2017.

[17]" "常士驃，張?zhí)K民.工程地質(zhì)手冊[M].4版.北京：中國建筑工業(yè)出版社，2007.

（收稿日期：2022-03-09，修回日期：2022-04-01）

Analysis of Safe Operation of Buried Coastal Natural Gas

Pipelines under Non?uniform Settlement

MA Xiao?ming， ZHAN Di

（School of Mechanical and Automotive Engineering， South China University of Technology）

Abstract" "The non?uniform settlement of foundations seriously influences operation safety of the buried natural gas pipeline. For purpose of ensuring the safe operation of the pipeline， the stress monitoring of the buried natural gas pipeline in coastal areas was carried out to obtain additional circumferential and axial stress of the pipeline influenced by non?uniform settlement of the foundation， and to check the pipeline stress in combination with working internal pressure of the pipeline. In addition， based on ANSYS software， the finite element model of the buried natural gas pipeline was established and the actual stress monitoring value and the numerical simulation results were contrastively analyzed to verify validity of the stress monitoring data and feasibility of the finite element model; meanwhile， influence of the pipeline’s inner diameter， wall thickness and buried soil parameters on the pipeline stress were further explored through finite element numerical simulation to show that， in the pipeline laying stage， reducing pipeline’s inner diameter， appropriately increasing its wall thickness and selecting harder buried soils can decrease the influence of the foundation’s non?uniform settlement on the pipeline. The research results provide the reference for safe operation of the buried natural gas pipelines affected by the foundation’s non?uniform settlement.

Key words" "non?uniform settlement， buried natural gas pipeline， stress analysis， safe operation

作者簡介：馬小明（1962-），副教授，從事安全檢測與失效分析、事故調(diào)查與分析、液化天然氣安全技術(shù)與裝備等研究工作，mars_scut@163.com。

引用本文：馬小明，詹迪.非均勻沉降下沿海天然氣埋地管道安全運(yùn)行分析[J].化工機(jī)械，2023，50（1）：76-82.