王小東,郭海燕,高秦嶺,孟 丹

(中國(guó)海洋大學(xué)工程學(xué)院,山東青島266100)

研究簡(jiǎn)報(bào)

鋼懸鏈線立管觸地點(diǎn)區(qū)域管土相互作用的有限元分析

王小東,郭海燕,高秦嶺,孟 丹

(中國(guó)海洋大學(xué)工程學(xué)院,山東青島266100)

由于鋼懸鏈線立管具有非線性特性,而海床土體又是軟黏土,因此鋼懸鏈線立管觸地區(qū)域的管土的相互作用十分復(fù)雜。根據(jù)國(guó)外相關(guān)試驗(yàn)數(shù)據(jù),采用ANSYS中的非線性彈簧單元模擬海床土體,考慮海床土體剛度退化和土吸力對(duì)管道的作用,建立海底管道擬靜力有限元計(jì)算模型,計(jì)算分析管道與海床土體的相互作用,并探討管道觸地點(diǎn)區(qū)域關(guān)鍵點(diǎn)在頂端升沉運(yùn)動(dòng)下彎矩的變化規(guī)律,為進(jìn)一步研究SCR與海床的動(dòng)力相互作用提供參考。

鋼懸鏈線立管(SCR);土體吸力;ANSYS有限元分析

隨著海洋石油天然氣需求量的增加,深海油氣資源開發(fā)中一種新型的立管結(jié)構(gòu)鋼懸鏈線立管(SCR)正被廣泛用于各種深水開發(fā)項(xiàng)目中。

觸地點(diǎn)(TDP)是鋼懸鏈線立管最初與海床接觸部位,是懸垂段和拖地段的連接點(diǎn),也是發(fā)生疲勞破壞可能性較大的位置。當(dāng)SCR的懸垂段頂端隨上部浮體運(yùn)動(dòng)時(shí),其拖地段將與海床土體發(fā)生相互作用。在拖地段的反復(fù)作用下,海床土體將形成溝槽,而溝槽的形成對(duì)SCR的拔出平面運(yùn)動(dòng)有較大影響。在溝槽形成的過程中,管線的運(yùn)動(dòng)將受到土的阻力作用,其中,阻力的大小取決于海底土的剛度,因而海底土剛度對(duì)SCR觸地區(qū)域的疲勞壽命有較大的影響[1]。深水油氣田海底土大部分屬于軟質(zhì)黏土(例如墨西哥灣的海底泥線處強(qiáng)度為2.6 kPa,隨深度呈線性增長(zhǎng),增長(zhǎng)幅度為1.5 kPa/m),而軟質(zhì)黏土在海底產(chǎn)生的土體吸力嚴(yán)重影響海底管線的上拔運(yùn)動(dòng),這雖對(duì)管線疲勞損傷的影響較小,但對(duì)其極限應(yīng)力的影響較大[2]。海底土吸力的影響因素較多,一般認(rèn)為與管道拔出的速度、管土反復(fù)作用的次數(shù)以及海底土重塑的時(shí)間等因素均有關(guān)系[3]。

關(guān)于SCR中觸地區(qū)域的研究,目前主要是試驗(yàn)的方法和數(shù)值模擬的方法。JUN G HWAN YOU[4]討論了海床與鋼懸鏈線立管的數(shù)值模型,提出了簡(jiǎn)化的彈簧支撐模型;Charles P,Jun Zhang等[5]研究了鋼懸鏈線立管與海床的相互作用,提出了P-y曲線及等價(jià)梁-彈簧模型;Bridge C[6]做了管土的試驗(yàn)并提出了土吸力模型等。這些研究大多假定海底土是線彈性的,實(shí)際深水海床土體對(duì)管道的作用以非線性模型模擬更合理。本文通過定義合適的土體剛度值,用等效的三折線模型模擬土體并分析管土相互作用,這種簡(jiǎn)約模型與實(shí)際情況更加接近,可以較好的模擬拖地段在拔起過程中海床土體剛度變化和土吸力對(duì)其運(yùn)動(dòng)的影響。

1 觸地點(diǎn)的管土相互作用計(jì)算模型

深水海洋立管鋪設(shè)在海床上,由于懸垂段頂端隨浮體運(yùn)動(dòng)觸地點(diǎn)將與海底發(fā)生相互作用,浮體運(yùn)動(dòng)包括浮體的一階運(yùn)動(dòng)(平衡位置附近的水平振蕩和升沉運(yùn)動(dòng))、浮體的二階慢漂(Slow drift)運(yùn)動(dòng),前者引起觸地點(diǎn)的改變[7],后者將引起觸地點(diǎn)的高應(yīng)力疲勞循環(huán)[8]。浮體運(yùn)動(dòng)下觸地區(qū)域管-土動(dòng)力相互作用計(jì)算十分復(fù)雜,本文主要進(jìn)行立管在頂端升沉運(yùn)動(dòng)下觸地點(diǎn)區(qū)域管土相互作用的擬靜力有限元分析。

1.1 管土相互作用

管道與海床土體豎向相互作用包括管道對(duì)土體的加卸載過程;海床土體對(duì)管道的抵抗力和吸附作用。本文計(jì)算數(shù)據(jù)來源于2001年2H Offshore Engineering在英國(guó)Watchet港口進(jìn)行的關(guān)于SCR與海床相互作用的大尺度二維模型試驗(yàn)[6],其土體特性類似于墨西哥灣,用于分析深水海床的土體特性以及與管道的相互作用。

KatherineLaver,EdClukey[2]等人根據(jù)STRIDE[9]的試驗(yàn)數(shù)據(jù)提出海底土體在管道垂向作用下特殊的滯回曲線模型——壓縮—回彈曲線模型(見圖1),土體在初始加載作用下的骨架曲線以及卸載和反向再加載下的滯回曲線,其中Audibert[10]提出土體在管道垂向作用下的滯回曲線可以選用雙曲線模型進(jìn)行近似模擬。

圖1 管土相互作用曲線(不分離)Fig.1 The curve of pipe and soil interaction(not separate)

1.2 土吸力模型

隨著管道向上運(yùn)動(dòng)與土體失去接觸,土體抵抗力消失,在這個(gè)過程中,土體對(duì)管道產(chǎn)生土吸力,抵抗管道繼續(xù)向上運(yùn)動(dòng)?，F(xiàn)有的海床土吸力模型是根據(jù)Watchet港口的模型試驗(yàn)得出的,其隨SCR相對(duì)于海床豎向位移變化的模型見圖2,分為3個(gè)線性階段[11]:吸力產(chǎn)生擴(kuò)大期 隨著管道向上運(yùn)動(dòng)土體的抵抗力消失,吸力從0增加到最大值;吸力平穩(wěn)期 隨著管道繼續(xù)向上運(yùn)動(dòng)吸力保持不變;吸力消退期 管道繼續(xù)向上運(yùn)動(dòng),吸力由最大逐漸減小,在破壞位置變?yōu)?。

圖2 二維的土吸力模型Fig.2 Two-dimensional soil suction model

土吸力變?yōu)?后,管道運(yùn)動(dòng)不受土體的影響。目前對(duì)于土吸力形成機(jī)理的研究處于起步階段,其卸載模型還不夠完善。

1.3 SCR的管土相互作用計(jì)算模型

管道與海床土體在觸地區(qū)域的相互作用包括2部分:管道對(duì)海床土體的加卸載部分和管道離開海床時(shí)土體的吸力對(duì)其產(chǎn)生的反作用力部分。不考慮管道上部浮體升沉運(yùn)動(dòng)速度的影響,觸地區(qū)域模擬土體對(duì)管道的作用模型簡(jiǎn)化為2部分:壓縮區(qū)域土體的模型用Audibert提出的雙曲線模型;拉伸區(qū)域的模型采用土吸力模型。

本文運(yùn)用ANSYS中的非線性彈簧單元對(duì)土體進(jìn)行模擬,土體的雙曲線模型用基于P-y曲線的等效彈簧系數(shù)來近似模擬,其中壓縮部分中的加載采用等效的三折線代替,根據(jù)文克爾的假定,土體在任意點(diǎn)的反力只與該點(diǎn)的位移成正比;吸力模型部分用與土吸力模型等效的非線性彈簧單元來模擬。等效彈簧的壓縮—拉伸與位移的關(guān)系曲線見圖3。

圖3 管土相互作用力—位移曲線Fig.3 The curve of force-displacement of pipe and soil interaction

2 有限元建模

采用有限元軟件ANSYS進(jìn)行建模及力學(xué)分析。其中,彈簧單元用COMBIN39進(jìn)行模擬;管道用PIPE59單元模擬。PIPE59是ANSYS程序中專門用于模擬浸沒在水中桿件結(jié)構(gòu)的單元,可以較好地模擬海洋波浪、海流對(duì)海中桿件的作用。

觸地點(diǎn)區(qū)域管道和土體的建模難度較大,本文將管道和土體在觸地點(diǎn)處分為2部分建立初始模型。圖4為立管及其觸地點(diǎn)的有限元計(jì)算模型,其中管道建模時(shí)分為初始觸地點(diǎn)D右側(cè)的懸垂段與左側(cè)的拖地段, C,I,J為觸地區(qū)域關(guān)鍵點(diǎn);C點(diǎn)、I點(diǎn)分別距離D點(diǎn)為6 m;J點(diǎn)距離D點(diǎn)為18 m。懸垂段在接觸到海床面前不受彈簧單元的反力,而拖地段置于一系列非線性彈簧單元上;土體模型對(duì)應(yīng)管道模型分為兩部分,D點(diǎn)左側(cè)土體用非線性彈簧單元模擬直接與管道接觸,右側(cè)土體在與管道接觸前不受壓力,當(dāng)管道向下運(yùn)動(dòng)到海床面土體開始與管道接觸時(shí),用非線性彈簧模擬土體。拖地段左端O點(diǎn)鉸接,頂端水平方向約束,豎向的荷載通過位移增量法施加升沉運(yùn)動(dòng)方式實(shí)現(xiàn)。模型中管道單元總數(shù)為220,結(jié)點(diǎn)數(shù)為221;非線性彈簧單元數(shù)為180。

本文管道和土體的數(shù)據(jù)來源于Watchet港口實(shí)驗(yàn)[2],管道參數(shù)見表1、土體參數(shù)見表2。

表1 管道參數(shù)Table 1 The parameters of pipe

表2 土體參數(shù)Table 2 The parameters of soil

圖4 立管及其觸地點(diǎn)的有限元計(jì)算模型Fig.4 The finite element calculation model of pipe and the touch down point

3 模型驗(yàn)證與計(jì)算結(jié)果分析

3.1 模型驗(yàn)證

管道鋪設(shè)在剛性的土體上,此時(shí)土體在管道作用下變形很小也不產(chǎn)生土吸力;在管道頂端(自由端)用增量法施加位移,對(duì)建立的有限元模型進(jìn)行擬靜力計(jì)算。管道頂端垂向由-0.8 m向上移動(dòng)到1.5 m,同時(shí)水平方向由-0.28 m運(yùn)動(dòng)到0.5 m,初始觸底點(diǎn)D及其它觸地區(qū)域關(guān)鍵點(diǎn)在管道頂端運(yùn)動(dòng)影響下產(chǎn)生了水平的移動(dòng),但是變化很小,最大變化范圍約為0.03 m,而本文中單元大小為0.5 m,因而忽略了節(jié)點(diǎn)水平位移的變化。D點(diǎn)彎矩變化如圖5,其中實(shí)線為數(shù)值模擬結(jié)果,虛線為STRIDE III J IP的試驗(yàn)數(shù)據(jù)圖[6]。

圖5 剛性土體D點(diǎn)處的彎矩變化圖Fig.5 D point bending variation of rigidity soil

由圖可知,隨著頂端向上運(yùn)動(dòng),數(shù)值模擬計(jì)算的D點(diǎn)彎矩變化與試驗(yàn)結(jié)果基本相同。

3.2 土吸力的影響與分析

按照本文建立的模型模擬管土相互作用,土體部分為有土吸力和無土吸力2種情況,用增量法對(duì)建立的有限元模型進(jìn)行計(jì)算,管道頂端(自由端)垂向由-0.8 m向上移動(dòng)到1.5 m處,同時(shí)水平方向由-0.28 m運(yùn)動(dòng)到0.5 m,節(jié)點(diǎn)的水平位移變化忽略不計(jì)。D點(diǎn)處的彎矩變化如圖6(有土吸力)和圖7(無土吸力),其中實(shí)線為數(shù)值模擬結(jié)果,虛線是STRIDE III J IP的試驗(yàn)數(shù)據(jù);圖8為初始觸地區(qū)域關(guān)鍵點(diǎn)J在土體包含土吸力和無土吸力作用的數(shù)值模擬結(jié)果的比較,其中實(shí)線為有土吸力作用,虛線為無土吸力作用。

圖6 D點(diǎn)處的彎矩變化圖Fig.6 D point bending variation

圖7 D點(diǎn)處的彎矩變化圖Fig.7 D point bending variation

由圖6和7可知,初始觸地點(diǎn)D處彎矩變化的數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果基本吻合,驗(yàn)證了用非線性彈簧模擬海床土體的可行性。從圖6中代表包含土吸力作用的實(shí)線可知,D點(diǎn)在頂端運(yùn)動(dòng)到-0.4 m前彎矩值基本保持不變,當(dāng)頂端運(yùn)動(dòng)到0.5 m處時(shí)彎矩值逐漸減小到-6.3 kNm,此時(shí)的彎矩值約為沒有土吸力作用時(shí)的2倍;隨著頂端繼續(xù)向上運(yùn)動(dòng)到1.2 m處時(shí)彎矩值逐漸增加到-3.6 kNm,之后基本保持不變。圖7中代表不包含土體吸力的實(shí)線可知,當(dāng)管道頂端為-0.8 m時(shí)D點(diǎn)的彎矩為5 kNm,隨著頂端向上運(yùn)動(dòng)D點(diǎn)彎矩逐漸減小;當(dāng)頂端運(yùn)動(dòng)到1.0 m處時(shí)D點(diǎn)彎矩減小到-3 kNm,之后基本保持不變。

圖8 J點(diǎn)處的彎矩變化圖Fig.8 J point bending variation

由圖8可知,隨著管道從-0.8 m向上移動(dòng)到1.5 m處,J點(diǎn)處的彎矩變化在土體有土吸力作用和無土吸力作用的數(shù)值模擬結(jié)果基本一致,其彎矩值在管道頂端移動(dòng)到0.8 m前基本保持不變,然后逐漸增大,可見遠(yuǎn)離觸地點(diǎn)處的土吸力對(duì)管道彎矩影響不大。

4 結(jié)語

由于鋼懸鏈線立管具有非線性特性,海床土體又是軟黏土,管土相互作用對(duì)觸地區(qū)域關(guān)鍵點(diǎn)的疲勞影響較大。本文用非線性彈簧模擬海床土體,建立管道與土體相互作用模型,用ANSYS有限元計(jì)算分析管道觸地區(qū)域各點(diǎn)的彎矩變化,并與國(guó)外試驗(yàn)進(jìn)行分析比較,探討了土體的非線性特性以及海床土吸力對(duì)管道彎矩的影響,驗(yàn)證了用非線性模型模擬土體的可行性。

由于海床土體大多是軟黏土,在研究管土相互作用時(shí),應(yīng)考慮土體的非線性特性以及土吸力的作用。本文用等價(jià)的三折線模型模擬海床土體的非線性特性,可以較好的反映海床土體在管道作用下剛度的變化以及對(duì)管道運(yùn)動(dòng)的影響。在對(duì)比無土吸力和存在土吸力的土體對(duì)管道彎矩的影響可知,土吸力對(duì)管道運(yùn)動(dòng)以及彎矩的變化具有明顯的作用,尤其對(duì)觸地區(qū)域的影響比較明顯。

雖然數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果比較吻合,但是本文并未考慮土體的摩擦系數(shù)對(duì)管道彎矩的影響;又由于土體吸力的大小與拖地段的拔出速度以及管土相互作用的次數(shù)、土的重塑時(shí)間都有關(guān),因此需要對(duì)土體的摩擦系數(shù)以及土體吸力的機(jī)理作進(jìn)一步的研究。

[1] 黃維平,李華軍.深水開發(fā)的新型立管系統(tǒng)——鋼懸鏈線立管(SCR)[J].中國(guó)海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2006,36(5):775-780.

[2] Bridge C,LaverK,Clukey E,et al.Steel catenary riser touchdown point vertical interaction models[C].∥Proc of the 2004 Offshore Technology Conference.Houston USA:OTC, 2004:16628.

[3] Theti R,Moros T.Soil interaction effects on simply catenary riser response[J].Pipe&Pipeline International,2001,46(3):15-24. [4] Jung Hwan You.Numerical model for steel catenary riser on seafloor support[D].Texas:Texas A&M University,2005.

[5] Charles P,Giovanna Biscontin,Jun Zhang.Seafloor interaction with steel catenary risers[D].Texas:Texas A&M University, 2006.

[6] Bridge C,Willis N.Steel Catenary Risers results and conclusions from large scale simulations of seabed interaction[C].14th Annual Conference Deep Offshore Technology,2002.

[7] Daniel Karunakaran,Kjell M Lund,Nils T Nordsve.Steel catenary riser configurations for North Sea Field developments[C]. //Proceedings Offshore Technology Conference.USA:OTC, 1999,2:331-338.

[8] Ghiath Manour.The impact of the second order vessel motion on the fatigue life of steel catenary risers[C].//Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering.Canada:ASME,2004,1:1177-1180.

[9] 2H Offshore Engineering Ltd-“STRIDE JIP-Effects Of Riser/Seabed Interaction On SCRs”[R].1500-RPT-008,Rev 02,2002.

[10] 2H Offshore Engineering Ltd-"Soil Stiffness Study"[R].1662-RPT-001,Rev 01,2004.

[11] 2H Offshore Engineering Ltd-"STRIDEJ IP-Pullout Resistance of a Pipe in a Clay Soil"[R].1500-RPT-006,Rev 02,2002.

Abstract: Because steel catenary riser(SCR)has nonlinear characteristics and the seabed soil is soft, SCR/seafloor interaction in the touch down point is very complex.This paper which is based on the related experimental data abroad uses nonlinear spring element of ANSYS to simulate seabed soil,and considers not only the effect of the degradation of stiffness of seabed soil but also the suction effect to the pipe and establishes submarine pipeline’s quasi-static finite element model.By using the model calculating the pipe interaction with the seabed soil,discussing the key points of the touch down region variation of moment under the heave motion and providing quite valuable reference for further dynamic study on the SCR/soil interaction.

Key words: steel catenary riser(SCR);soil suction;finite element analysis(ANSYS)

責(zé)任編輯 陳呈超

Finite Element Analysis of Steel Catenary Riser/seafloor Interaction in the Touch Down Zone

WANG Xiao-Dong,GUO Hai-Yan,GAO Qin-Ling,MENG Dan
(College of Engineering,Ocean University of China,Qingdao 266100,China)

TV312

A

1672-5174(2010)09Ⅱ-197-05

國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2006AA09Z356);國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(50739004)資助

2010-01-20;

2010-06-07

王小東(1983-),男,碩士生。E-mail:wangxiaodong111aa@163.com