基于有限元法的含凹陷油氣管道應(yīng)變分析

陳健++雷錚強(qiáng)等

摘要： 為對(duì)含凹陷的油氣管道進(jìn)行更合理的完整性評(píng)價(jià)，采用有限元法模擬管道在承受外部巖石擠壓時(shí)的應(yīng)變分布，討論凹陷深度、管道壁厚和擠壓體大小等參數(shù)對(duì)應(yīng)變的影響.結(jié)果表明：在相同條件下，最大等效應(yīng)變隨凹陷深度增大而增大；在同樣深度條件下，管壁較厚處的凹陷最大等效應(yīng)變較大；擠壓體半徑越大，凹陷的輪廓越光滑，應(yīng)變集中程度越??；管材等級(jí)對(duì)等效應(yīng)變的影響很小.該方法和結(jié)論可用于改進(jìn)凹陷評(píng)價(jià)方法和風(fēng)險(xiǎn)排序.

關(guān)鍵詞： 油氣管道； 凹陷； 完整性評(píng)價(jià)； 失效； 風(fēng)險(xiǎn)排序； 應(yīng)變； 有限元法

中圖分類號(hào)： TE88文獻(xiàn)標(biāo)志碼： B

0引言管道是油氣運(yùn)輸?shù)闹饕侄?，是我?guó)能源輸送的大動(dòng)脈.截止到2013年底，我國(guó)長(zhǎng)輸油氣管道總長(zhǎng)度已達(dá)10.6萬km，預(yù)計(jì)到2015年將達(dá)15萬km.油氣管道在我國(guó)國(guó)民經(jīng)濟(jì)建設(shè)中起重要作用.但是，長(zhǎng)輸油氣管道壓力高、管徑大，一旦發(fā)生失效，會(huì)導(dǎo)致嚴(yán)重的后果.

凹陷是由于管道與其他物體的物理接觸導(dǎo)致管道橫截面發(fā)生變形，是長(zhǎng)輸油氣管道的最常見缺陷之一.凹陷可能發(fā)生在管道施工期間，由搬運(yùn)、回填過程中的碰撞或巖石障礙等原因?qū)е?；也可能發(fā)生在管道服役期間，由挖掘設(shè)備、巖石等外物的壓砸等原因?qū)е?凹陷對(duì)管道安全運(yùn)行的影響表現(xiàn)在2個(gè)方面：一方面，凹陷會(huì)引起管體局部應(yīng)力、應(yīng)變集中，尤其當(dāng)凹陷與裂紋、劃傷或焊縫缺陷相關(guān)時(shí)，容易導(dǎo)致管道發(fā)生斷裂失效；另一方面，凹陷縮減管道的有效內(nèi)徑，進(jìn)而影響管內(nèi)介質(zhì)的輸送，同時(shí)阻礙清管器和內(nèi)檢測(cè)器等設(shè)備在管中的正常運(yùn)行，給管道的管理帶來困難.因此，如何對(duì)含凹陷管道進(jìn)行合理的完整性評(píng)價(jià)和有針對(duì)性的修復(fù)成為亟待解決的問題.

2基于有限元的凹陷應(yīng)變分析

2.1有限元模型

對(duì)于實(shí)際的含凹陷管道，其應(yīng)力、應(yīng)變分布很難測(cè)量，而在有限元模擬中可以很方便地獲得應(yīng)力和應(yīng)變等參數(shù)，有助于判斷含凹陷管體的最危險(xiǎn)位置，更好地研究管體的失效行為.對(duì)于管道凹陷的有限元模擬，其難點(diǎn)在于問題的非線性，同時(shí)需要考慮接觸非線性、幾何非線性和材料非線性.

本文以Φ720 mm的X70鋼管道為主要研究對(duì)象，模擬管道在巖石擠壓下的變形，重點(diǎn)分析等效應(yīng)變的分布和最大值.對(duì)于管線鋼的非線性應(yīng)力應(yīng)變關(guān)系，使用RambergOsgood本構(gòu)模型描述.擠壓巖石選用花崗巖的材料參數(shù)，由于擠壓巖石相對(duì)于被擠壓管道來說剛度要大得多，因此可以直接使用線性應(yīng)力應(yīng)變關(guān)系.模擬中使用的主要材料參數(shù)見表1.

材料密度/（kg/m3）彈性模量/GPa屈服強(qiáng)度/MPa抗拉強(qiáng)度/MPa

管材（X70）7.8206483570

巖石（花崗巖）2.6780

為更好地模擬被擠壓管道的非線性行為，同時(shí)合理減少計(jì)算時(shí)間，使用高階殼結(jié)構(gòu)單元SHELL 281作為管體的單元類型，擠壓巖石使用高階實(shí)體結(jié)構(gòu)單元SOLID 186.擠壓巖石與管道之間的接觸采用面面接觸進(jìn)行定義.由于被擠壓管道剛度較低，將其定義為“目標(biāo)面”，使用三維目標(biāo)單元TARGE 170.擠壓巖石剛度較大，將其定義為“接觸面”，使用三維8節(jié)點(diǎn)接觸單元CONTA 174.

最終劃分的有限元網(wǎng)格見圖2.由于被擠壓管道僅在局部產(chǎn)生較大變形，所以只截取被擠壓側(cè)1/2管道進(jìn)行建模.擠壓巖石假定為球形，同樣也只截取1/2建模.

2.2模擬結(jié)果和分析

圖3給出Φ720 mm×8 mm的X70鋼管道被4 mm半徑球形巖石擠壓出10%管道外徑深度凹陷后的等效應(yīng)變分布圖，可以看到等效應(yīng)變最大值出現(xiàn)在凹陷底部最深處，達(dá)到0.339，往外迅速減小，距離凹陷較遠(yuǎn)處應(yīng)變很小.

表2給出Φ720 mm×8 mm的X70鋼管道被4 mm半徑球形巖石分別擠壓出5%，10%，15%和20%管道外徑深度凹陷后的模擬結(jié)果.在同樣條件下等效應(yīng)變最大值隨著凹陷深度增大而增大，可見，在同樣條件下，凹陷越深對(duì)管道完整性的影響越大.

表3給出壁厚分別為8，9和10 mm的Φ720 mm X70鋼管道被4 mm半徑球形巖石擠壓出5%管道外徑深度凹陷后的模擬結(jié)果，可見，在同樣條件下等效應(yīng)變最大值隨著壁厚增大而增大.這是由于管道壁厚越大剛度越大，產(chǎn)生同樣深度的凹陷實(shí)際應(yīng)變集中程度越大.因此，在同樣深度條件下，應(yīng)優(yōu)先修復(fù)壁厚較大處的凹陷.

表4給出Φ720 mm×8 mm的X70鋼管道分別被4，8和10 mm半徑球形巖石擠壓出5%管道外徑深度凹陷后的模擬結(jié)果，可見，在同樣條件下等效應(yīng)變隨著擠壓巖石的半徑增大而減小.這是由于擠壓巖石半徑越大，凹陷的輪廓越平滑，應(yīng)變集中程度也越小.因此，在同樣深度條件下，應(yīng)優(yōu)先修復(fù)軸向長(zhǎng)度或環(huán)向?qū)挾容^小的凹陷.

表5給出Φ720 mm×8 mm的X70和X52鋼管道分別被4 mm半徑球形巖石擠壓出5%管道外徑深度凹陷后的模擬結(jié)果.結(jié)果表明：兩者的等效應(yīng)變最大值非常接近，X70鋼管道的結(jié)果稍大于X52鋼管道的.這是由于被擠壓管道的變形主要由剛度決定，不同等級(jí)管材的彈性模量相近，只是強(qiáng)度有所區(qū)別，因此，剛度相近，變形程度也相近.X70鋼管道剛度略高，因此產(chǎn)生同樣深度的凹陷實(shí)際應(yīng)變集中程度略大.還需要指出：當(dāng)管材為X52時(shí)，求解所需要的時(shí)間大幅度增加.這是由于單元?jiǎng)偠仍降?，非線性程度越大，為滿足收斂條件進(jìn)行迭代計(jì)算的次數(shù)也越多.因此，在使用有限元進(jìn)行非線性分析時(shí)，經(jīng)常會(huì)在合理的范圍內(nèi)人為提高剛度以減少計(jì)算時(shí)間.在凹陷等管體變形模擬中也可以視情況考慮使用這種做法.

3結(jié)論

針對(duì)油氣管道凹陷評(píng)價(jià)問題，采用有限元法模擬管道在承受外部巖石擠壓時(shí)的應(yīng)變分布，討論凹陷深度、管道壁厚和擠壓體大小等參數(shù)對(duì)等效應(yīng)變最大值的影響.研究結(jié)果表明：有限元模擬可以得到工程實(shí)際中難以測(cè)量的應(yīng)力和應(yīng)變等參數(shù)，有助于凹陷失效風(fēng)險(xiǎn)的評(píng)估.本文的分析方法和結(jié)論可用于驗(yàn)證和改進(jìn)基于應(yīng)變的凹陷評(píng)價(jià)方法，也對(duì)凹陷風(fēng)險(xiǎn)排序和制定修復(fù)計(jì)劃提供指導(dǎo).

參考文獻(xiàn)：

[1]楊瓊， 帥健， 左尚志. 管道凹陷研究現(xiàn)狀[J]. 油氣儲(chǔ)運(yùn)， 2009， 28（6）： 1015.

YANG Qiong， SHUAI Jian， ZUO Shangzhi. Research actuality of pipelines with dents[J]. Oil & Gas Storage & Transportation， 2009， 28（6）： 1015.

[2]李明菲， 周利劍， 任重， 等. 在役長(zhǎng)輸油氣管道管體凹陷檢測(cè)評(píng)價(jià)研究進(jìn)展[C]//2013中國(guó)國(guó)際管道會(huì)議暨第一屆中國(guó)管道與儲(chǔ)罐腐蝕與防護(hù)學(xué)術(shù)交流會(huì)論文集. 北京： 石油工業(yè)出版社， 2013： 194200.

[3]LEIS B N， FORTE T P， ZHU X. Integrity analysis for dents in pipelines[DB/OL]. （20041231）[20141001]. http：//proceedings.asmedigitalcollection.asme.org/proceeding.aspx？articleid=1646155.

[4]CSA Z66211Oil and gas pipeline systems[S].

[5]ASME B31.8—2010Gas transmission and distribution piping systems[S].

[6]ASME B31.4—2009Pipeline transportation systems for liquid hydrocarbons and other liquids[S].

[7]CFR 49192—2013Transportation of natural and other gas by pipeline： minimum federal safety standards[S].

[8]CFR 49195—2013Transportation of hazardous liquids by pipeline[S].

[9]NORONHA Jr D B， MARTINS R R， JACOB B P， et al. Procedures for the strain based assessment of pipeline dents[J]. Int J Pressure Vessels & Piping， 2010， 87（5）： 254265.

[10]LUKASIEWICZ S A， CZYZ J A， SUN C， et al. Calculation of strains in dents based on high resolution inline caliper survey[DB/OL]. （20061231）[20141001]. http：//proceedings.asmedigitalcollection.asme.org/proceeding.aspx？articleid=1596829.

[11]GAO M， McNEALY R， KRISHNAMURTHY R. Strainbased models for dent assessment： a review[DB/OL]. （20081231）[20141001]. http：//proceedings.asmedigitalcollection.asme.org/proceeding.aspx？articleid=1640691&resultClick=3.

（編輯于杰）