海底柔性復(fù)合管隆起屈曲有限元分析

， ，

(1.中海石油(中國(guó))有限公司天津分公司 天津 300450； 2.石油管材及裝備材料服役行為與結(jié)構(gòu)安全國(guó)家重點(diǎn)實(shí)驗(yàn)室 陜西 西安 710065； 3.西安交通大學(xué)機(jī)械工程學(xué)院 陜西 西安 710049)

·試驗(yàn)研究·

海底柔性復(fù)合管隆起屈曲有限元分析

龔海潮1，魏斌2，李兵3

(1.中海石油(中國(guó))有限公司天津分公司 天津 300450； 2.石油管材及裝備材料服役行為與結(jié)構(gòu)安全國(guó)家重點(diǎn)實(shí)驗(yàn)室 陜西 西安 710065； 3.西安交通大學(xué)機(jī)械工程學(xué)院 陜西 西安 710049)

柔性復(fù)合管在海洋油氣輸送領(lǐng)域發(fā)揮著重要作用，管道隆起屈曲分析是海洋管道設(shè)計(jì)的重要技術(shù)之一。柔性復(fù)合管具有多層結(jié)構(gòu)，每層材料性質(zhì)差異較大，且層間應(yīng)力關(guān)系復(fù)雜，使得隆起屈曲分析有別于傳統(tǒng)的單層剛性管道。利用有限元分析軟件ANSYS，研究了管道內(nèi)流體溫度、壓力以及埋泥深度三種因素對(duì)多層柔性管道隆起屈曲程度的影響。結(jié)果表明：對(duì)于多層柔性管道，輸送流體溫度越高，管道隆起高度越高，輸送流體壓力升高將限制管道的隆起，埋泥深度增大將限制管道的隆起。研究工作和分析結(jié)果將為海底管道鋪設(shè)提供理論依據(jù)和設(shè)計(jì)參考。

海底柔性管；隆起屈曲；有限元分析；溫度；壓力；埋泥深度

0 引 言

柔性復(fù)合管是連接海底石油或天然氣井口與浮式平臺(tái)輸送油氣或注水的重要海洋裝備,一般是由聚合物阻隔層和各種金屬螺旋鎧甲層、非金屬耐磨層等結(jié)構(gòu)層復(fù)合而成。柔性復(fù)合管主要優(yōu)點(diǎn)有：良好的柔性、安裝速度快、安全可靠性高、可采用模塊化的設(shè)計(jì)方式、抗腐蝕性好、抗壓性能好、可回收重新利用等七方面。柔性管在海洋石油資源開發(fā)中具有良好的環(huán)境適應(yīng)能力和經(jīng)濟(jì)效益，在國(guó)內(nèi)外多個(gè)海洋油氣田的開發(fā)過程中得到了廣泛的應(yīng)用[1]。

近年來，海底管線輸送流體的工作溫度隨著油氣井開發(fā)深度的增大而不斷地提高，很多管線設(shè)計(jì)溫度普遍達(dá)到或者超過80 ℃，甚至出現(xiàn)120 ℃的高溫[2]。海底埋設(shè)管道在這種高溫高壓條件下運(yùn)行時(shí)，有時(shí)會(huì)從原來位置上突然隆起，甚至拱出埋設(shè)土層，這種現(xiàn)象稱為海底管道隆起屈曲[3]，如圖1所示。該過程類似于梁在軸向載荷達(dá)到臨界值時(shí)發(fā)生的歐拉屈曲[4]。當(dāng)管道隆起屈曲時(shí)，鋼管會(huì)產(chǎn)生較大的彎矩以及較大的塑性變形，甚至可能進(jìn)一步出現(xiàn)裂紋、疲勞、局部屈曲[5]。這些問題給管道安全運(yùn)行帶來巨大安全隱患，挪威船級(jí)社(DNV)規(guī)范RP4-110《Global buckling of submarine pipelines-structural design due to high temperature and high pressure pipelines》明確要求在海洋管道設(shè)計(jì)時(shí)對(duì)這種情況進(jìn)行詳細(xì)分析。

圖1 管道隆起示意圖

現(xiàn)有的管道設(shè)計(jì)規(guī)范雖然都詳細(xì)的描述了管道溫度應(yīng)力的計(jì)算方法，但都是以單層管為對(duì)象[6]。并且當(dāng)存在初始缺陷情況下，解析及求解非常困難，因此通常采用數(shù)值模擬分析具有初始隆缺陷的海底管道的隆起屈曲[7]。由于多層柔性管道在結(jié)構(gòu)上與傳統(tǒng)的單層剛性管道差別較大，然而傳統(tǒng)的分析單層剛性管道方法運(yùn)用在多層柔性管道中有著計(jì)算精度低，計(jì)算量大等問題。

本文針對(duì)一新型七層柔性海洋管道，采用有限元軟件ANSYS，建立了七層管道的高精度有限元模型，分析了管道溫度、壓力和埋泥深度等不同因素對(duì)其隆起屈曲變形的影響，研究工作為海底管道施工提供重要的理論依據(jù)和技術(shù)參考。

1 柔性海洋管有限元建模

海洋柔性管的內(nèi)徑為140 mm，外徑為250 mm，主要結(jié)構(gòu)為7層，端面結(jié)構(gòu)如圖2所示，有里層到外層依次為內(nèi)壓屏蔽層、壓力鎧裝層、抗拉鎧裝層、中間屏蔽層、抗外壓鎧裝層和外保護(hù)層。各層材料類型和功能見表1。

圖2 多層海洋柔性管道橫截面圖

層數(shù)名稱材料功能1內(nèi)壓屏蔽層聚偏氟乙烯提供內(nèi)部流體完整性通道2壓力鎧裝層Q195L鋼承受內(nèi)壓荷載3抗拉鎧裝層Q235B鋼4抗拉鎧裝層Q235B鋼平衡扭轉(zhuǎn),承受軸向及部分環(huán)向荷載5中間屏蔽層耐熱聚乙烯阻止內(nèi)外介質(zhì)二次滲透,保護(hù)金屬結(jié)構(gòu)層6抗外壓鎧裝層Q195鋼承受外壓荷載7外保護(hù)層耐熱聚乙烯阻止外部流體進(jìn)入,保護(hù)柔性管內(nèi)部結(jié)構(gòu)

多層柔性管道主要受軸向載荷影響，軸向載荷的大小取決于多種因素，如鋪設(shè)過程中產(chǎn)生的張力、管道運(yùn)行過程中內(nèi)壓、溫度的變化所產(chǎn)生的軸向載荷、以及管道與海底的摩擦等。溫度載荷是管道發(fā)生隆起屈曲的關(guān)鍵控制因素[8-10]。假設(shè)有初始形狀的管道鋪設(shè)在海床上，兩端無軸向位移，內(nèi)管受溫度載荷作用。對(duì)該問題建立溫度載荷有限元分析模型，內(nèi)外管采用ANSYS中Elbow290單元[11,12]。該單元是ANSYS公司最新開發(fā)的專門針對(duì)多層管道有限元建模的多層管單元。該單元能夠適應(yīng)線性大轉(zhuǎn)角情況，同時(shí)也適用于非線性大應(yīng)力等非線性情況。該單元可以定義每一層的材料性質(zhì)進(jìn)行建模分析，適用于復(fù)合材料管道分析，其復(fù)合管道計(jì)算精確度受第一剪切應(yīng)力準(zhǔn)則控制。計(jì)算中采用的復(fù)合管各層材料性質(zhì)見表2。

表2 有限元計(jì)算中應(yīng)用的柔性復(fù)合管各層材料性質(zhì)

由于管道隆起對(duì)于管道長(zhǎng)度有所限制，長(zhǎng)度過短會(huì)影響隆起部分的隆起形狀和應(yīng)力分布，因此建模長(zhǎng)度選取為200 m。初始幾何缺陷是海底管道整體屈曲分析中非常關(guān)鍵的因素，基于有限元分析軟件建立的海底管道整體屈曲分析的有限元模型必須引入幾何缺陷才能進(jìn)行整體屈曲分析[13]。管道初始形狀可用式(1)表示[14]：

(1)

其中，ω0為初始缺陷高度，L0為初始缺陷長(zhǎng)度，缺陷形狀假設(shè)關(guān)于x=0對(duì)稱，初始缺陷高度為0.3 m，初始缺陷長(zhǎng)度為10 m。

在隆起屈曲-5 m至5 m處每0.1 m劃分一個(gè)單元，其余每1 m劃分一個(gè)單元。管道沿周向劃分40等份，建模完成共劃分290個(gè)單元，建模后效果如圖3和圖4所示。

圖3 管道建模整體示意圖

圖4 管道建模截面示意圖

2 邊界條件及載荷施加

計(jì)算過程中需要考慮兩種邊界條件，一是管道接頭兩端的固定作用，二是地面的支撐作用。柔性復(fù)合管單支長(zhǎng)度一般為50～100 m，中間靠金屬接頭連接，因此計(jì)算過程中可將其簡(jiǎn)化為管道兩端的六自由度完全約束。其物理意義在于接頭會(huì)限制管道的軸向的移動(dòng)和轉(zhuǎn)動(dòng)，垂直于地面的移動(dòng)和轉(zhuǎn)動(dòng)還有平行于地面方向的移動(dòng)和轉(zhuǎn)動(dòng)，恰好完全限制了六個(gè)自由度。對(duì)于地面支撐，沒有隆起的管道直接限制垂直于地面的移動(dòng)。根據(jù)管道三種實(shí)際工況分別施加如下載荷：

1)溫度載荷：對(duì)于Elbow290單元，溫度載荷對(duì)其來說為體載荷，可直接設(shè)定每一層表面的溫度進(jìn)行計(jì)算。

2)內(nèi)壓載荷：對(duì)于Elbow290單元，壓力載荷對(duì)其為面載荷，可將壓力直接施加在某一層的表面上。處理時(shí)將壓力載荷直接施加在管道內(nèi)層表面上，壓力載荷正方向的定義是垂直于表面向外為正方向，對(duì)于內(nèi)表面來說，內(nèi)壓方向?yàn)榇怪庇趦?nèi)表面向內(nèi)，因此內(nèi)壓載荷的設(shè)定為負(fù)。

3)埋泥載荷：對(duì)于Elbow290單元，可以將埋泥深度轉(zhuǎn)化為節(jié)點(diǎn)力直接施加在變形單元的節(jié)點(diǎn)上，其線載荷可表述為式(2)。

Rclay=γHD+D2(1/2-π/8)+2SU(H+D/2)

(2)

式中：γ為土壤浮重，γ取6 900 N/m3；H為管道上部覆蓋土體高度(海床至管道最上端)；D為管道外徑；Su為管道中心到溝槽頂部的平均不排水剪切強(qiáng)度。Su取2 000 N/m2。管道輸送流體溫度為80 ℃，內(nèi)部壓力為20 MPa。

3 計(jì)算結(jié)果

3.1 熱應(yīng)力對(duì)管道隆起屈曲的影響

由于需要對(duì)七層管道每一層施加溫度載荷，因此在管道屈曲分析之前進(jìn)行傳熱分析，以確定每一層管道的每個(gè)面的溫度情況。進(jìn)行建模之后，施加邊界條件，內(nèi)層溫度為變量從80～30 ℃變化，外層溫度定為0 ℃，進(jìn)行計(jì)算得到溫度分布云圖如圖5所示。

圖5 管道橫截面溫度分布云圖

將以上熱分析結(jié)果帶入到Elbow290單元體載荷中進(jìn)行及計(jì)算，打開非線性大變形開關(guān)，將管道橫截面溫度分布數(shù)據(jù)帶入，以初始不平度高度為0.3 m，初始不平度長(zhǎng)度為10 m，內(nèi)層溫度為80 ℃為例，管道屈曲形狀如圖6和圖7所示。最大應(yīng)力以及截面應(yīng)力分布如圖8所示。

圖6 管道變形圖

圖7 管道變形局部變形云圖

圖8 管道屈曲界面應(yīng)力分布分布云圖

改變內(nèi)層溫度，經(jīng)計(jì)算得到溫度與管道隆起屈曲垂直方向位移關(guān)系如圖9所示。由圖9可以看出，內(nèi)層溫度直接影響了管道的隆起屈曲程度。隨著溫度的升高，隆起屈曲的位移不斷變大。產(chǎn)生上述結(jié)果的原因是溫度升高使材料發(fā)生了膨脹，并在軸向方向產(chǎn)生了壓應(yīng)力，又由于初始不平度的影響導(dǎo)致的歐拉桿的失穩(wěn)，最終導(dǎo)致管道發(fā)生隆起屈曲。另一方面是由于壓力鎧裝層的材料彈性模量較大，在同等變形的情況下會(huì)產(chǎn)生更大的應(yīng)力，容易造成整個(gè)結(jié)構(gòu)的塑性屈服。但是熱導(dǎo)率很低的最內(nèi)層材料聚偏氟乙烯減小了向外的熱量傳遞，限制了壓力鎧裝層的溫度升高，進(jìn)一步限制了壓力鎧裝層因?yàn)闇囟鹊淖兓斐傻呐蛎洠罱K減小了整個(gè)結(jié)構(gòu)的位移變形，因此內(nèi)層材料使用熱導(dǎo)率低的材料是十分重要的。從圖8可以看出，最大應(yīng)力發(fā)生在第六層的外壓鎧裝層。由于第六層厚度比較小，會(huì)受到整個(gè)內(nèi)層的擠壓而變形，并且因?yàn)槠涮幱谕鈱游恢茫韵鄳?yīng)的變形量和應(yīng)力都較大。

圖9 隆起屈曲高度與輸送流體溫度關(guān)系圖

圖10 隆起屈曲高度與管道內(nèi)壓關(guān)系圖

3.2 內(nèi)壓載荷對(duì)管道隆起屈曲的影響

使用SFE面載荷命令將壓力加載到內(nèi)層管道中，管道內(nèi)部壓力從0 MPa至20 MPa每隔2 MPa進(jìn)行一次計(jì)算，以輸送流體壓力為20 MPa，初始不平度高度為0.3 m，初始不平度長(zhǎng)度為10 m，內(nèi)層液體溫度為80 ℃。改變輸送流體壓力，計(jì)算后得到內(nèi)部液體壓力與管道隆起屈曲的關(guān)系如圖10所示。由圖10可以看出，隆起高度隨管道內(nèi)壓的升高而降低，提高壓力能夠減小管道隆起屈曲的程度。分析原因是因?yàn)楣艿纼?nèi)部壓力越大，管道會(huì)因?yàn)閴毫Φ淖饔糜幸欢ǖ捏w積膨脹，此時(shí)層與層之間的正壓力就會(huì)增加，正壓力的增加就會(huì)間接地增加最大靜摩擦力以及滑動(dòng)摩擦力。而摩擦力的增加會(huì)阻礙管道的層間滑動(dòng)。因此增大內(nèi)部液體壓力會(huì)限制管道的隆起屈曲。

3.3 埋泥深度對(duì)管道隆起屈曲的影響

根據(jù)式(2)能夠計(jì)算出在不同埋泥深度下，管道隆起需要克服的垂直方向作用力。經(jīng)過計(jì)算埋泥深度和垂直方向作用力的關(guān)系，改變不同的埋泥深度，計(jì)算后得到埋泥深度與管道隆起高度之間的關(guān)系如圖11所示。在管道上方埋泥不僅有重力的作用，并且還有管道突破土層需要克服排水剪切壓力，因此埋泥是最能夠限制管道隆起屈曲的因素。但是考慮到生產(chǎn)成本和管道鋪設(shè)難度的問題，埋泥深度不宜過大。計(jì)算結(jié)果表明埋泥深度在0.3 m左右，管道隆起現(xiàn)象已經(jīng)被極大地抑制。

圖11 隆起屈曲高度與埋泥深度關(guān)系圖

4 結(jié) 論

采用有限元軟件ANSYS，建立了七層柔性復(fù)合管的隆起屈曲有限元模型，分析了管內(nèi)溫度、壓力和埋泥深度三種主要因素對(duì)管道隆起屈曲變形的影響，得到如下結(jié)論：

1)管內(nèi)溫度越高，管道隆起高度越高；

2)提高輸送流體壓力能夠減小管道隆起屈曲的程度，內(nèi)壓在14～20 MPa之間對(duì)管道隆起屈曲有較好抑制效果；

3)埋泥也能夠限制管道的隆起屈曲，埋泥深度超過0.3 m，管道隆起已被明顯抑制。

[1] 余志兵，孫國(guó)民，王 輝，等. 海底長(zhǎng)輸油氣管道的隆起屈曲分析[J]. 中國(guó)海洋平臺(tái)，2015，30(2):64-69.

[2] 劉新帥，孫國(guó)民，劉志剛，等. 海底埋設(shè)高溫管道隆起屈曲分析研究[J]. 中國(guó)石油與化工標(biāo)準(zhǔn)與質(zhì)量，2016，36(15):93-94.

[3] 曾霞光，段夢(mèng)蘭，車小玉. 海底埋設(shè)雙層管管道隆起屈曲分析[J]. 海洋工程，2014，32(2):72-77.

[4] 趙天奉，段夢(mèng)蘭，潘曉東，等. 雙層海底管道跨越設(shè)計(jì)的垂向研究[J]. 中國(guó)海上油氣，2010，22(3):197-201.

[5] DENG Huifeng, SONG Chunna, DAI Wanbao. The first use pipeline plough on subsea pipeline trenching in South China Sea [C]//Proceedings of the Twentieth (20m) International Offshore and Polar Engineering Conference, 20-25 May 2010, Beijing, China.

[6] 趙天奉.高溫海底管道溫度應(yīng)力計(jì)算與屈曲模擬研究[D].大連理工大學(xué),2008.

[7] 車小玉，段夢(mèng)蘭，曾霞光，等. 海底埋設(shè)高溫管道隆起屈曲數(shù)值模擬研究[J]. 海洋工程,2013,31(5):103-111.

[8] James GA Croll. A Simplified model of upheaval thermal buckling subsea pipelines[J]. Thin-Walled Structures ,1997,29 (1-4) :59-78.

[9] SUN Jason, WANG James. The advancement of FEA in confronting the deepwater pipelines under high pressure and high temperature [C]//Offshore Technology Conference, 2-5 May 2011, Houston, USA.

[10] NFA Parsloe , P Manchec, A Dorbec. Methods for assessing design issues with HP/HT pipelines at early stages a project [C]//Offshore Technology Conference, 1-4, May 2006, Houston, USA.

[11] K.L.巴 斯. 工程分析中的有限元法[M]. 北京: 機(jī)械工業(yè)出版社,1999:48-86.

[12] 李 兵. Ansys Workbench設(shè)計(jì)、仿真與優(yōu)化[M].清華大學(xué)出版社,2008:20-80.

[13] 劉羽霄，葛 濤，李 昕，等. 初始幾何缺陷在海底管道整體屈曲數(shù)值建模中的引入方法[J]. 水資源與水工程學(xué)報(bào)，2011，22(3):32-35.

[14] KARAMPOUR H, ALBERMANI F, GROSS J. On lateral and upheaval buckling of subsea pipelines[J]. Engineering Structures,2013,5(2): 317-330.

FiniteElementAnalysisandCalculationofUpheavalBucklingofSubmarineFlexibleCompositePipe

GONGHaichao1,WEIBin2,LIBing3

(1.ChinaNationalOffshoreOil(China)Co.Ltd.,TianjinBranch,Tianjin300450,China； 2.StateKeyLaboratoryforPerformanceandStructureSafetyofPetroleumTubularGoodsandEquipmentMaterials,Xi′an,Shaanxi710077,China; 3.SchoolofMechanicalEngineering,Xi′anJiaotongUniversity,Xi′an,Shaanxi710049,China)

Submarine pipes are gradually widely used in ocean engineering. However, when the seabed has a prop for the pipeline that used in high temperature, the buckling behavior will easily occur. The safety of pipeline is heavily affected by the buckling behavior. More multilayered flexible pipelines are used because of complex situation in ocean and engineering demand. And the upheaval buckling more easily occurs in multilayered flexible pipelines compared to traditional pipes. However, because of the many numbers of layers, the difference among different layers and the complexity of stress between two different layers, the analyzing of multilayer flexible pipeline is different from traditional single inflexible pipeline and more difficult. According to the problems above all, the buckling behavior affected by different factors is analyzed by software ANSYS. Finite element model is established. Three different factors are discussed. Temperature pressure and the depth which the pipeline is buried in are analyzed in this paper. The results show three conclusion. First, the higher temperature will cause the heavier buckling behavior. Second, the more depth will limit the buckling behavior. Third, the higher pressure will also limit the buckling behavior. The result will provide theory and reference with constructing submarine pipeline.

submarine flexible pipeline; upheaval buckling; FEM analysis; temperature; pressure; depth of embedded mud

龔海潮，男，1983年生，工程師，2006年畢業(yè)于河北工業(yè)大學(xué)工程管理專業(yè)，現(xiàn)從事海上油氣田開發(fā)工程項(xiàng)目管理工作。E-mail:gonghch@cnooc.com.cn

TE832.2

A

2096-0077(2017)06-0029-05

10.19459/j.cnki.61-1500/te.2017.06.008

2017-07-24

葛明君)