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

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

中圖分類號 TE973" "文獻標識碼 A" "文章編號 0254?6094（2023）01?0076?07

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

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

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

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，通過對管道站場周期性（月）測試記錄，管道運行過程中工作壓力從2021年1月的6.75 MPa逐步調(diào)升至2021年10月的7.61 MPa，出站管道處地基沉降從21.2 mm增至78.6 mm，出站管道實際工況記錄見表1，這種非均勻沉降現(xiàn)象會導致出站的埋地管道因受載不均而發(fā)生變形，管道應力隨之發(fā)生變化，應力的變化實際是埋地管道與管周土壤相互作用的結(jié)果[10]，應力過大會導致管道開裂，發(fā)生泄漏，危害天然氣分輸站運行安全，影響下游天然氣供應。因此對分輸站管道進行應力監(jiān)測與分析對管道安全意義重大。

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

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

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

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

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

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

σ=（ε+με） （1）

σ=（ε+με）（2）

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

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

ε——所測軸向應變;

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

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

σ=（3）

σ=（4）

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

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

δ——管道壁厚，mm。

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

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

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

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

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

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

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

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

σ=σ+σ≤0.9σ（5）

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

σ——管道材料最小屈服強度，MPa。

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

3 管道應力有限元分析

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

3.1 有限元模型

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

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

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

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

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

σ=（6）

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

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

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

3.3 管道應力影響因素

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

3.3.1 管道內(nèi)徑

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

3.3.2 管道壁厚

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

3.3.3 埋土參數(shù)

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

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

4 結(jié)論

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

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

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

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

參 考 文 獻

[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].石油機械，2014，42（3）：104-109.

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

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

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

[8]" "劉玉卿，余志峰，佟雷，等.基于軸向應力監(jiān)測數(shù)據(jù)的管道應力狀態(tài)預警模型[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è)科技大學，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強度準則的巷道圍巖數(shù)值模擬對比分析[J].礦業(yè)研究與開發(fā)，2019，39（12）：28-34.

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

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

[16]" "康遜.天然氣管道在非均勻沉降狀態(tài)下的應力測試與分析[D].廣州：華南理工大學，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。

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