伍穎 李都 陳朗

摘要：由于我國西南地區(qū)山脈連綿，存在大量高落差埋地輸氣管道，并且隨著服役年限延長也會出現(xiàn)不同程度的腐蝕。同時山區(qū)高落差埋地含腐蝕管道在地震作用下受力復(fù)雜，而目前難以通過實地監(jiān)測或者在試驗中完成其地震響應(yīng)研究。為此，基于西南地區(qū)某一典型高落差埋地X80管道的實際工況，建立了地震作用下高落差埋地含腐蝕X80管道有限元模型，探討了腐蝕深度、腐蝕寬度和腐蝕長度對高落差埋地含腐蝕管道地震動力響應(yīng)的影響規(guī)律。研究結(jié)果表明：在地震波和內(nèi)壓的加載下，腐蝕深度的增加會使最大等效應(yīng)力呈線性向上增加趨勢，而內(nèi)壓和敷設(shè)角度會影響管道發(fā)生腐蝕后的初始應(yīng)力值；最大等效應(yīng)力隨腐蝕長度的增大而增大，當(dāng)敷設(shè)角度為60°時，腐蝕長度超過0.3D（D為管道外徑）后，敷設(shè)角度對管道最大等效應(yīng)力的影響超過了腐蝕長度；最大等效應(yīng)力隨腐蝕寬度的增大而減??；通過參數(shù)敏感性分析得出，腐蝕深度對地震作用下高落差埋地管道最大等效應(yīng)力影響最大（占比為0.71），其次為腐蝕長度（占比為0.27），腐蝕寬度影響最?。ㄕ急葹?.02）。在天然氣管道設(shè)計和施工階段，應(yīng)避免敷設(shè)角度大于45°，對于在役管段出現(xiàn)腐蝕處應(yīng)重點監(jiān)測其腐蝕深度。所得結(jié)論可為山區(qū)長輸管道的設(shè)計、運營安全和工程管理提供參考。

關(guān)鍵詞：埋地輸氣管道；高落差；腐蝕缺陷；地震響應(yīng)；有限元

In mountainous Southwest China， there are a large number of high drop buried gas pipelines， which have corroded to different degrees with the service life.Moreover， these corroded pipelines were subject to complex forces under the action of earthquakes；and unfortunately， their seismic responses can hardly be obtained by field monitoring or tests.Considering the actual working conditions of a typical high drop buried X80 pipeline in Southwest China， a finite element model of high drop buried corroded X80 pipeline under seismic action was built， and the influences of corrosion depth， corrosion width and corrosion length on the seismic dynamic response of high drop buried corroded pipeline were identified.The study results show that under the loading of seismic wave and internal pressure， the increase of corrosion depth makes the maximum equivalent stress increase linearly upward， while the internal pressure and laying angle affects the initial stress value of the pipeline after corrosion.The maximum equivalent stress increases with the increase of corrosion length.When the laying angle is 60° and the corrosion length exceeds 0.3D， the influence of the laying angle on the maximum equivalent stress of the pipeline exceeds that of the corrosion length.The maximum equivalent stress decreases with the increase of corrosion width.The parameter sensitivity analysis shows that the corrosion depth has the greatest influence on the maximum equivalent stress of high drop buried pipeline under seismic action， accounting for 0.71， followed by corrosion length， accounting for 0.27， and the corrosion width has the least influence， accounting for 0.02.Therefore， in the design and construction stage of natural gas pipeline， the laying angle should be avoided to be greater than 45°， and for the corrosion area of inservice pipeline， monitoring should focus on the corrosion depth.The conclusions provide reference for the design， operation safety and project management of longdistance pipelines in mountainous areas.

buried gas pipeline；high drop；corrosion defect；seismic response；finite element

0 引 言

我國地勢復(fù)雜，西高東低，山地面積占全國國土面積的70％。在天然氣管道規(guī)劃、設(shè)計及施工建設(shè)過程中，大部分長輸天然氣管道都不可避免地要穿越地面障礙物，但因隧道穿越施工難度大、成本高，又因高陡坡地帶山體底面積較大，不宜采用隧道穿越或繞行的方式［1］，一般采用沿坡埋地敷設(shè)方式，從而形成眾多高落差（大落差）埋地管道。高落差埋地管道具有坡度陡、落差大及彎管多的特點，且隨著服役年限的延長及內(nèi)部介質(zhì)和外部環(huán)境的影響使其難以避免遭受腐蝕，在地震作用下結(jié)構(gòu)極易發(fā)生應(yīng)力集中。因此高落差埋地含腐蝕缺陷管道的抗震安全性也就備受關(guān)注。

自20世紀(jì)60年代，大量學(xué)者就缺陷管道的失效行為和失效評估進(jìn)行了研究，形成了多個評價腐蝕缺陷管道性能的標(biāo)準(zhǔn)和規(guī)范，如美國的ASME B31G［2］和API 579［3］、挪威船級社的DNV RP-F101［4］、英國的BS 7910［5］等。天然氣管道腐蝕是一個連續(xù)可變的過程，隨著服役時間的延長，天然氣管道由于內(nèi)部介質(zhì)以及土壤因素等出現(xiàn)不同程度的腐蝕，導(dǎo)致力學(xué)和抗震性能持續(xù)下降。因此，如果不考慮管道發(fā)生腐蝕這一情況，可能會導(dǎo)致破壞性甚至災(zāi)難性的后果。2007年，劉威等［6］建立了地震作用下動態(tài)線性腐蝕管道分析模型，并解釋其作用。2008年，LIU W.等［7］基于彈性地基梁法和隨機(jī)攝動法建立了地震作用下腐蝕管道響應(yīng)方法，并用該方法計算了地震作用下管道位移和應(yīng)力的均值和標(biāo)準(zhǔn)差，且驗證了其準(zhǔn)確性。2012年，QI H.B.等［8］采用有限元法分別對有、無缺陷埋地輸油管道進(jìn)行地震響應(yīng)分析，結(jié)果表明，腐蝕缺陷對管道地震響應(yīng)影響較大。2015年，管友海等［9］建立了地震作用下架空管道與內(nèi)部流體的液固耦合模型，指出無論是內(nèi)腐蝕還是外腐蝕，當(dāng)管道腐蝕深度超過管壁厚度的70％時，在地震作用下管道都會有破壞危險。2016年，譚現(xiàn)東等［10］建立了地震作用下含腐蝕缺陷的架空輸油管道有限元模型，分析了不同腐蝕深度、腐蝕位置以及場地土下的地震反應(yīng)規(guī)律。2017年，王桂萱等［11］以埋地腐蝕管線為研究對象，探究了腐蝕深度、寬度、腐蝕位置和埋深等參數(shù)對管線抗震性能的影響。2018年，LI Q.Z.等［12］利用有限元軟件建立了熱流固多場耦合力學(xué)模型，分析了流體壓力、流體溫度、腐蝕缺陷和地震波對管道力學(xué)性能的影響，結(jié)果指出，相對于腐蝕深度，流體壓力和流體溫度對腐蝕管道的力學(xué)性能有明顯影響。2020年，劉海水［13］運用有限元軟件研究了地震作用下不同溫度梯度對管道腐蝕缺陷處應(yīng)力分布的影響。2021年，WANG Y.H.等［14］研究了管道地震作用下局部腐蝕對管道安全性的影響。同年，鄭山鎖等［15］采用數(shù)值模擬軟件進(jìn)行非線性動力時程分析，揭示了不同土壤環(huán)境下多齡期埋地鋼管的地震響應(yīng)規(guī)律。

綜上所述，目前對腐蝕管道抗震研究各有側(cè)重點：一方面重點關(guān)注腐蝕在時域上的變化，評估地震作用下腐蝕管道的安全性；另一方面關(guān)注熱流固多場耦合力學(xué)，期望達(dá)到數(shù)值模擬的真實性。但都只是關(guān)注埋地直管的腐蝕，缺乏對腐蝕彎管的研究，對于高落差埋地含腐蝕管道在地震作用下的力學(xué)研究更是少之又少。為此，筆者根據(jù)山區(qū)管道腐蝕類型，選取了管道腐蝕簡化模型；基于ABAQUS軟件建立了管土相互作用模型，腐蝕缺陷建立在高應(yīng)力區(qū)域；探討了不同腐蝕深度、腐蝕寬度及腐蝕長度等因素對高落差埋地管道力學(xué)特性的影響規(guī)律；通過敏感性分析得到高落差含腐蝕缺陷埋地管道在地震作用下腐蝕參數(shù)對管道應(yīng)力的影響程度，以期為山區(qū)長輸管道的設(shè)計施工、運營安全和工程管理提供理論參考。

1 管道腐蝕類型的選取

由于天然氣管道服役時間的延長，管道腐蝕現(xiàn)象日益嚴(yán)重，即使管道企業(yè)采取加陰極保護(hù)系統(tǒng)、涂覆防腐層等多種防護(hù)措施，仍然無法避免腐蝕。根據(jù)某檢測機(jī)構(gòu)的檢測數(shù)據(jù)［16］，檢測發(fā)現(xiàn)腐蝕缺陷共計20 706處，其中內(nèi)腐蝕20 635處占比98.2%，外腐蝕71處占比1.8%。因此，本文采用內(nèi)腐蝕缺陷來模擬高落差埋地含腐蝕X80管道的地震動力響應(yīng)特性。

內(nèi)腐蝕缺陷形狀主要有蝕坑和溝槽狀腐蝕，常用橢圓形、均勻壁厚、拋物線形進(jìn)行規(guī)則化處理，處理形式如圖1所示。對于壁厚為t的天然氣管道，拋物線形腐蝕缺陷簡化模型是以實際腐蝕最大深度a和腐蝕缺陷的投影長度l來表示，簡化的面積為A=23la。等深度矩形簡化模型是以實際腐蝕深度平均值和腐蝕缺陷投影長度來表示，簡化后面積為A=0.85la，即簡化的平均深度為0.85a。

等深度腐蝕缺陷在模擬中腐蝕位置的剩余厚度基本相同，腐蝕底部是與管道同軸的圓柱面；橢圓形腐蝕為內(nèi)凹的球面，通常用來模擬缺陷長度和寬度較小的蝕坑，腐蝕中間深兩邊淺；拋物線腐蝕底面類似拋物線，對于寬度較大的缺陷不易模擬。為分析含缺陷管道在缺陷處附近應(yīng)力和應(yīng)變的變化規(guī)律，需研究不同腐蝕深度、腐蝕長度及腐蝕寬度對管道的影響，并且腐蝕長、寬、深變化范圍較大。因此，本文選取等深度腐蝕作為腐蝕缺陷的簡化模型。

2 地震作用下高落差含腐蝕管道有限元模型

2.1 有限元模型建立

由于西南地區(qū)某天然氣管道沿線地形地質(zhì)條件復(fù)雜，多數(shù)管線處于抗震設(shè)防烈度為6烈度及以上的山區(qū)，具有長距離、大口徑、高壓力及高落差的特點，所以本文基于西南地區(qū)某一典型高落差埋地X80管道實際工況進(jìn)行模型建立，并且選取有限元軟件ABAQUS進(jìn)行建模分析。

2.1.1 條件假設(shè)

在數(shù)值建模時，做出如下基本假設(shè)：

（1）管道埋地段采用同一種土體，土體不產(chǎn)生大變形，且未考慮材料的各向異性。

（2）不考慮溫度及管材隨機(jī)性對管道性能的影響，以及管道連接方式對管道力學(xué)性能的影響。

（3）假設(shè)地震載荷空間一致性，即管道兩端同步振動。

2.1.2 材料參數(shù)

（1）管道參數(shù)。管道設(shè)計運行壓力為10 MPa，鋼材采用X80級鋼管，管道外徑為1 016 mm，彈性模量為207 GPa，密度為7 850 kg/m3，泊松比為0.3，屈服應(yīng)力為555 MPa。采用隨動雙線性硬化彈塑性模型，該模型遵從von Mises屈服準(zhǔn)則，其表達(dá)式如下：

（2）土體參數(shù)。為了充分了解土壤物理性能以及為后續(xù)有限元數(shù)值模擬提供重要的土體參數(shù)，按照土工試驗方法標(biāo)準(zhǔn)［18］進(jìn)行土的物理性能試驗。試驗先對土樣進(jìn)行測定，當(dāng)土樣的基本形狀符合試驗要求后再進(jìn)行大量取土。本文主要利用電子天平、沙浴鍋、全自動氣壓固結(jié)儀、應(yīng)變控制式直剪儀等儀器測定土樣的天然密度、壓縮模量以及內(nèi)摩擦角等基本參數(shù)。

經(jīng)試驗測得土體的密度為1 880 kg/m3，內(nèi)摩擦角為14.7°，壓縮模量為6.58 MPa，黏聚力為20 kPa，為某種黏土。根據(jù)文獻(xiàn)［19］取其彈性模量為18 MPa，泊松比為0.38。由于DruckerPrager彈塑性本構(gòu)模型對土體的彈塑性特征描述較為準(zhǔn)確，尤其是可以反映土體對管道的摩擦與支撐作用，所以土體的本構(gòu)模型采用ABAQUS軟件中自帶的DP模型。

2.1.3 模型簡化

本文選取的高落差埋地管道可分為上彎管和中部直管以及下彎管3個部分。管道高程為H=30 m，敷設(shè)角度α=50°，管材為X80管道，外徑D=1 016 mm，壁厚為t=15.3 mm，埋深為h=2 m，上彎管和下彎管曲率半徑ρ=6 m。由于計算效率的需求，在模型中只能考慮有限長度的管道，所以對管道進(jìn)行了幾何模型簡化，如圖3所示。

2.1.4 管土接觸

管土相互作用采用非線性接觸模型，該理論可以通過有限元軟件ABAQUS中的“硬”接觸以及主從接觸面來有效地實現(xiàn)管土接觸這種非線性行為。通過表面與表面接觸進(jìn)行力學(xué)約束，設(shè)置剛體較大的管道外表面為主表面，包裹管道的場地土表面為從表面。接觸屬性定義法向行為和切向行為，其中法向行為主要對土壓力進(jìn)行約束，切向行為主要對管道與場地土間摩擦力進(jìn)行約束。法向行為定義為“硬”接觸，接觸后可分離；切向行為通過“罰”函數(shù)進(jìn)行定義，摩擦因數(shù)取0.5［20］。

2.1.5 網(wǎng)格劃分

基于團(tuán)隊前期研究成果，高應(yīng)力區(qū)域主要位于下彎管附近，因此，為了分析最不利工況，本文將腐蝕缺陷建立在直管與下彎管連接處頂部。含腐蝕管道網(wǎng)格劃分需要確保計算的準(zhǔn)確性并縮短計算時間，具體的實現(xiàn)方法是：在腐蝕部位密集地布置網(wǎng)格，并對非腐蝕區(qū)域進(jìn)行粗化網(wǎng)格劃分。由于不同的尺寸模型劃分的網(wǎng)格數(shù)量不同，腐蝕處網(wǎng)格沿軸向、圓周方向以及徑向的網(wǎng)格尺寸分別為30、50以及2 mm，如圖4所示。

模型響應(yīng)結(jié)果可行性的必要條件是執(zhí)行網(wǎng)格敏感性分析，以落差高度30 m、敷設(shè)角度45°、彎管曲率6 m、腐蝕深度0.3t、腐蝕寬度0.3D、腐蝕長度0.3D為例，對高落差埋地含腐蝕管道依次加載重力（g=9.8 m/s2）、內(nèi)壓（10 MPa）及地震載荷（0.2g）。通過反復(fù)改變圓周方向、縱向方向和軸向腐蝕的網(wǎng)格尺寸，計算腐蝕管道的最優(yōu)網(wǎng)格數(shù)量，結(jié)果如表1所示。由表1可知，當(dāng)總網(wǎng)格數(shù)量達(dá)到12萬后，最大等效應(yīng)力趨于穩(wěn)定，所以本文采用方案4進(jìn)行數(shù)值模擬計算。

2.1.6 邊界條件

邊界條件一般包括簡單邊界、黏性邊界及黏彈性邊界3種。一般來講，黏彈性邊界的精度要高于黏性邊界和簡單邊界［21］。因此，本文在模型底部約束豎向位移，場地四周邊界則采用目前較為成熟的黏彈性人工邊界來處理場地邊界，人工邊界法向與切向的彈簧剛度參考谷音等［22］的研究成果：

2.1.7 地震動輸入

采用動力響應(yīng)法進(jìn)行地震反應(yīng)分析時，一般采用地震加速度時程作為地震動輸入。ABAQUS中地震波一般用表格（Table）的文本格式代替時程加速度曲線讀入到軟件中，然后將時程加速度作為邊界約束條件施加到模型底部上作為地震動輸入，分析步時間設(shè)置為0.02 s。水平方向上的地震動多用于地上結(jié)構(gòu)，垂直方向上的地震動多用于地下結(jié)構(gòu)。因此，本文選擇沿場地土底部向Y軸方向施加一致地震波激勵。地震波（EI波）輸入如圖5所示。

2.2 模型驗證

相似的數(shù)值模型或試驗可以執(zhí)行有限元模型驗證。在這項工作中，有限元模型模擬了遭受地震的高落差埋地腐蝕管道，可以將其驗證分解為地震響應(yīng)和局部腐蝕管道的爆破模型。本團(tuán)隊前期研究成果已證明前者，后者基于A.C.BENJIAMIN等［23］的爆破試驗數(shù)據(jù)。利用等效應(yīng)力對應(yīng)的極限屈服強(qiáng)度作為破壞準(zhǔn)則，數(shù)值模擬的破壞壓力為20.3 MPa。用下式計算其偏差值，結(jié)果為4.2×10-2％。由此證實了所建立的模型可行。

3 腐蝕缺陷對力學(xué)響應(yīng)的影響分析

3.1 缺陷尺寸選取

含內(nèi)腐蝕管道的腐蝕尺寸示意圖如圖6所示。a、w、l分別為腐蝕缺陷的深度、腐蝕缺陷的寬度及腐蝕缺陷的長度。通過上述變化設(shè)計多工況下的算例，開展地震作用下高落差埋地含腐蝕管道地震動力響應(yīng)分析。

由圖7可以看出，應(yīng)力云圖左右兩邊大致對稱，但上下不對稱。這是高落差埋地管道在Z方向上管道尺寸和腐蝕尺寸以及載荷對稱，而在Y方向上幾何以及載荷都不對稱造成的。

含腐蝕天然氣管道等效應(yīng)力云圖可以分為2個區(qū)域：遠(yuǎn)離腐蝕區(qū)域和腐蝕附近區(qū)域。遠(yuǎn)離腐蝕區(qū)域應(yīng)力基本均勻分布，而腐蝕附近區(qū)域應(yīng)力分布情況復(fù)雜多變，出現(xiàn)最大等效應(yīng)力，此時管道除了承受內(nèi)壓引起的薄膜應(yīng)力外，還要承受地震作用力以及腐蝕處幾何尺寸變化所引起的不連續(xù)應(yīng)力?？拷g尺寸邊界處深受不連續(xù)應(yīng)力的影響，極易發(fā)生強(qiáng)度失效［24］。

為了對比不同腐蝕深度下腐蝕區(qū)域最大應(yīng)力在管道環(huán)向的分布情況，圖8給出了等效應(yīng)力最大處YZ平面等效應(yīng)力分布隨腐蝕深度變化的關(guān)系曲線。由本團(tuán)隊的前期研究成果得知，內(nèi)壓和敷設(shè)角度分別是載荷和管道易損性的重要參數(shù)。為分析其影響規(guī)律，圖9給出了不同敷設(shè)角度以及不同內(nèi)壓下埋地管道最大等效應(yīng)力隨腐蝕深度的變化曲線。

由圖8可以看出，腐蝕區(qū)域兩端由于尺寸變化引起應(yīng)力集中，導(dǎo)致沿管道環(huán)向方向兩端應(yīng)力最大。隨著腐蝕深度增加，管道腐蝕區(qū)域的等效應(yīng)力增加，而非腐蝕區(qū)域的等效應(yīng)力幾乎沒有改變。

由圖9可以看出：不同腐蝕深度下的管道最大等效應(yīng)力呈線性增長，總體上隨著腐蝕深度的增加，最大等效應(yīng)力呈現(xiàn)線性向上增長趨勢，運行壓力和敷設(shè)角度并不影響其應(yīng)力增長趨勢，但會影響管道發(fā)生腐蝕后的起始等效應(yīng)力值；在敷設(shè)角度為30°和45°時，運行壓力對管道最大等效應(yīng)力的增幅程度遠(yuǎn)大于敷設(shè)角度，在60°時則反之；當(dāng)敷設(shè)角度為30°時，在運行壓力8 MPa下的管道最大等效應(yīng)力值增幅最大，約為160 MPa，但此工況下的管道均處于彈性階段。敷設(shè)角度為45°時，最大等效應(yīng)力隨著腐蝕深度的增加幾乎呈線性增長，在運行壓力為10 MPa下，腐蝕深度為0.5t時，管道達(dá)到屈服狀態(tài)，其余工況管道均處于彈性階段。敷設(shè)角度為60°時，最大等效應(yīng)力隨著腐蝕深度的增加大致呈線性增加，直至管道達(dá)到屈服狀態(tài)后，最大等效應(yīng)力增長逐漸平緩，在運行壓力分別為10、8和6 MPa下，腐蝕深度分別對應(yīng)達(dá)到0.2t、0.3t和0.5t時管道達(dá)到屈服狀態(tài)。因此，對于在役的高落差埋地X80管道發(fā)生腐蝕后，應(yīng)及時采取腐蝕修復(fù)措施來降低管道最大等效應(yīng)力增加的趨勢。

3.3 腐蝕寬度

為探究地震作用下高落差含不同腐蝕寬度埋地管道的力學(xué)影響規(guī)律，對不同腐蝕寬度（0.1D、0.2D、0.3D、0.4D、0.5D）高落差埋地管道依次加載重力（g=9.8 m/s2）、內(nèi)壓載荷（10 MPa）和地震載荷（0.2g），計算地震作用下不同腐蝕寬度對高落差埋地管道的影響。圖10給出了EI波下2.32 s時，不同腐蝕寬度下高落差埋地管道腐蝕區(qū)域的應(yīng)力云圖。

由圖10可以看出，隨著寬度增加，最大等效應(yīng)力逐漸減小，主要原因是腐蝕缺陷寬度增加之后，有更長的腐蝕缺陷減薄區(qū)來承受因內(nèi)壓而產(chǎn)生的塑性應(yīng)變，相同壓應(yīng)力集中水平略有下降，所以會出現(xiàn)隨著寬度增加最大等效應(yīng)力逐漸減小的趨勢。

為了對比不同腐蝕寬度下腐蝕區(qū)域最大應(yīng)力在管道環(huán)向的分布情況，圖11給出了等效應(yīng)力最大處YZ平面腐蝕區(qū)域等效應(yīng)力分布隨腐蝕寬度變化的關(guān)系曲線圖。圖12給出了不同敷設(shè)角度以及不同內(nèi)壓下埋地管道最大等效應(yīng)力隨腐蝕寬度的變化曲線。

由圖11可以看出：腐蝕缺陷中部等效應(yīng)力明顯比缺陷環(huán)向方向兩端應(yīng)力小；隨著腐蝕寬度增加，管道腐蝕區(qū)域的等效應(yīng)力減小，而非腐蝕區(qū)域的等效應(yīng)力幾乎沒有改變。由此可見，腐蝕寬度的增加僅會影響腐蝕區(qū)域應(yīng)力分布。

由圖12可以看出：不同腐蝕寬度下的管道最大應(yīng)力呈現(xiàn)分層現(xiàn)象，總體上隨著腐蝕寬度的增加，最大等效應(yīng)力呈現(xiàn)緩慢向下減小趨勢；在敷設(shè)角度為30°和45°時，運行壓力對管道最大等效應(yīng)力的增幅程度遠(yuǎn)大于敷設(shè)角度，在60°時則反之。在運行壓力為10 MPa時，敷設(shè)角度在30°和45°工況下，不同腐蝕寬度的最大等效應(yīng)力變化趨勢大體一致，都在緩慢減小，而在敷設(shè)角度為60°工況下，管道已經(jīng)達(dá)到屈服狀態(tài)，且腐蝕寬度的增加對最大等效應(yīng)力影響極小。在運行壓力為8 MPa時，敷設(shè)角度在30°和45°工況下，最大等效應(yīng)力減小趨勢較緩，在60°工況下，腐蝕寬度小于等于0.2D時，管道處于屈服狀態(tài)，在腐蝕寬度大于0.2D后，管道最大等效應(yīng)力減小趨勢增大。當(dāng)運行壓力為6 MPa時，敷設(shè)角度在30°和45°以及60°工況下，整個管道都處于彈性狀態(tài)。因此，當(dāng)在役高落差埋地X80管道出現(xiàn)腐蝕后，調(diào)低內(nèi)壓運行可以大幅降低管道的最大等效應(yīng)力。對于設(shè)計階段的高落差埋地管道應(yīng)盡量采取小于等于45°的角度進(jìn)行敷設(shè)。

3.4 腐蝕長度

為探究地震作用下高落差含不同腐蝕長度埋地管道的力學(xué)影響規(guī)律，對不同腐蝕長度（0.1D、0.2D、0.3D、0.4D、0.5D）高落差埋地管道依次加載重力（g=9.8 m/s2）、內(nèi)壓載荷（10 MPa）、地震載荷（0.2g），計算地震作用下不同腐蝕長度對高落差埋地管道的影響。圖13給出了EI波下2.32 s時，不同腐蝕長度下高落差埋地管道腐蝕區(qū)域的應(yīng)力云圖。

由圖13可以看出：高落差埋地含腐蝕管道在不同腐蝕長度下，高應(yīng)力區(qū)仍處于腐蝕尺寸邊界處；隨著腐蝕長度的增加，最大等效應(yīng)力增大，并且不斷向腐蝕缺陷兩端集中。為了對比腐蝕區(qū)域不同腐蝕長度下最大應(yīng)力在管道環(huán)向的分布情況，圖14給出了等效應(yīng)力最大處YZ平面等效應(yīng)力分布隨腐蝕深度變化的關(guān)系曲線。圖15給出了不同敷設(shè)角度以及不同內(nèi)壓下埋地管道最大等效應(yīng)力隨腐蝕長度的變化曲線。

由圖14可以看出：隨著腐蝕長度的增加，最大等效應(yīng)力逐漸增大；結(jié)合圖14當(dāng)腐蝕長度為0.1D時，最大等效應(yīng)力出現(xiàn)在腐蝕區(qū)域中間位置；當(dāng)腐蝕長度大于0.1D后，最大等效應(yīng)力出現(xiàn)在環(huán)向方向兩端。

由圖15可以看出：不同腐蝕長度下的管道最大應(yīng)力呈增長現(xiàn)象，總體上隨著腐蝕長度的增加，最大等效應(yīng)力呈現(xiàn)先增加后逐漸平緩趨勢；在運行壓力為6 MPa時，此工況下管道均處于彈性狀態(tài)，但在敷設(shè)角度為60°時，管道最大等效應(yīng)力一直處于增加，且腐蝕長度大于0.3D后的最大等效應(yīng)力增幅大于0.3D之前的。這說明在腐蝕長度超過0.3D后，敷設(shè)角度對管道最大等效應(yīng)力的影響超過了腐蝕長度。在運行壓力為8 MPa下，敷設(shè)角度為60°，腐蝕長度小于0.3D時，管道最大等效應(yīng)力上升趨勢較快，當(dāng)腐蝕長度為0.3D時，管道已經(jīng)達(dá)到屈服破壞，而敷設(shè)角度為45°以及30°時，此時等效應(yīng)力的增長速率逐漸減緩。在運行壓力為10 MPa下，比較3種埋地管道敷設(shè)角度，隨著腐蝕長度的增加，最大等效應(yīng)力增大。敷設(shè)角度在30°和45°工況下，不同腐蝕長度的最大等效應(yīng)力變化率趨近一致，腐蝕長度在0.1D～0.3D時最大等效應(yīng)力增長較快，在腐蝕長度大于0.3D后，最大等效應(yīng)力增長較慢，然后趨于平緩。在運行壓力為10 MPa，埋地管道敷設(shè)角度為60°時，腐蝕深度為0.3t的工況下，無論腐蝕長度大小，高落差埋地管道都處于屈服階段。這說明敷設(shè)角度和腐蝕深度對天然氣管道抗震性能非常重要，在天然氣管道設(shè)計、運行階段，應(yīng)避免敷設(shè)角度大于45°且重點監(jiān)測其腐蝕深度。

按式（7）分別計算各參數(shù)對管道最大等效應(yīng)力的敏感度系數(shù)，得到敏感度系數(shù)與各參數(shù)變化量的關(guān)系曲線，如圖16所示，其中黑色標(biāo)簽為對應(yīng)參數(shù)的參數(shù)變化量。

由圖16可以看出，管道最大等效應(yīng)力與腐蝕深度、腐蝕長度呈正方向變化，與腐蝕寬度呈反方向變化。在相同參數(shù)變化下，計算這3項參數(shù)的敏感性系數(shù)平均絕對值并做歸一化處理。由結(jié)果可得，腐蝕深度對地震作用下高落差埋地X80管道最大等效應(yīng)力影響最大，占比0.71；其次為腐蝕長度，占比0.27；腐蝕寬度對管道應(yīng)力影響最小，占比0.02。因此，對于山區(qū)天然氣管道的日常運行應(yīng)當(dāng)加強(qiáng)腐蝕監(jiān)測管理，對于出現(xiàn)腐蝕的管段，可采用瀝青、內(nèi)涂層油管［25］、3PE防腐層與外加電流陰極保護(hù)相結(jié)合的方法進(jìn)行腐蝕防護(hù)。

5 結(jié)論及認(rèn)識

（1）在無腐蝕管壁上，等效應(yīng)力分布比較均勻，越靠近腐蝕尺寸邊界，其等效應(yīng)力突變越大，最終最大等效應(yīng)力發(fā)生在腐蝕區(qū)域沿管道環(huán)向的兩個端點處。

（2）隨著腐蝕深度的增加，高落差埋地X80管道最大等效應(yīng)力呈線性向上增長趨勢，運行壓力和敷設(shè)角度并不影響其應(yīng)力增長趨勢，但會影響管道發(fā)生腐蝕后的起始等效應(yīng)力值，在敷設(shè)角度為30°和45°時，運行壓力對管道最大等效應(yīng)力的增幅程度遠(yuǎn)大于敷設(shè)角度，在60°時則反之。因此，在天然氣管道設(shè)計和施工階段，應(yīng)避免敷設(shè)角度大于45°，對于在役管段出現(xiàn)腐蝕處應(yīng)重點監(jiān)測其腐蝕深度。

（3）隨著腐蝕寬度的增加，高落差埋地X80管道最大等效應(yīng)力逐漸減小，減小的趨勢隨著腐蝕寬度的增加逐漸增大。當(dāng)運行壓力為6 MPa時，敷設(shè)角度在30°、45°以及60°工況下，整個管道都處于彈性狀態(tài)。因此，當(dāng)在役高落差埋地X80管道出現(xiàn)腐蝕后，調(diào)低內(nèi)壓運行可以大幅降低管道的最大等效應(yīng)力。

（4）隨著腐蝕長度的增加，高落差埋地X80管道最大等效應(yīng)力增大，增大的趨勢隨著腐蝕長度的增加逐漸減小。當(dāng)敷設(shè)角度為60°時，腐蝕長度超過0.3D后，敷設(shè)角度對管道最大等效應(yīng)力的影響超過了腐蝕長度，因此在高落差管道設(shè)計與運行階段，應(yīng)盡量減小敷設(shè)角度。

（5）通過參數(shù)敏感性分析得出，腐蝕深度對地震作用下高落差埋地含腐蝕X80管道等效應(yīng)力影響最大，占比0.71，其次為腐蝕長度，占比0.27，腐蝕寬度對管道應(yīng)力影響最小，占比0.02。因此，對于山區(qū)天然氣管道的日常管理應(yīng)當(dāng)加強(qiáng)腐蝕監(jiān)測，對于出現(xiàn)腐蝕的管段，需采取必要的腐蝕防護(hù)措施。

［1］ 徐潔，黃坤，廖檸，等．基于應(yīng)力分析的高陡邊坡輸氣管道安全性分析［J］．中國安全生產(chǎn)科學(xué)技術(shù)，2015，11（12）：110-115．

XU J，HUANG K，LIAO N，et al.Safety analysis on gas pipeline at high and steep slope based on stress analysis［J］.Journal of Safety Science and Technology，2015，11（12）：110-115.

［2］ ASME.Manual for determining the remaining strength of corroded pipelines：ASME B31G—2009［S］.New York：ASME，2009.

［3］ API.Fitnessforservice：API 579-1［S］.Washington，DC：American Petroleum Institute，2021.

［4］ DNV.Corroded pipelines：DNV-RP-F101［S］.Oslo：DNV，2021.

［5］ BSI.Guidance on methods for assessing the acceptability of flaws in fusion welded structures：PD 6493：1991［S］.London：British Standards Institution，1991.

［6］ 劉威，李杰．考慮隨機(jī)腐蝕作用的埋地管線地震反應(yīng)分析［J］．土木工程學(xué)報，2007，40（2）：104-108．

LIU W，LI J.Seismic response analysis of buried pipelines with stochastic corrosions［J］.China Civil Engineering Journal，2007，40（2）：104-108.

［7］ LIU W，LI J.Stochastic seismic response of pipelines with corrosion［J］.Journal of Earthquake Engineering，2008，12（6）：914-931.

［8］ QI H B，WANG Q S，WANG L L，et al.Seismic response research on oil pipeline with corrosion defects［J］.Applied Mechanics and Materials，2012，166/169：2124-2127.

［9］ 管友海，王亞洲，賈娟娟，等．含不同腐蝕缺陷輸油管道抗震性能分析［J］．化工設(shè)備與管道，2015，52（1）：62-66．

GUAN Y H，WANG Y Z，JIA J J，et al.Analysis of seismic property of oil transportation pipeline with various corrosion defects［J］.Process Equipment & Piping，2015，52（1）：62-66.