趙晶瑞，范 模，段夢蘭

(1.中海油研究總院，北京 100027;2.中國石油大學(xué)(北京)，北京 102249)

深水半潛式生產(chǎn)平臺耦合動力特性對比研究

趙晶瑞1，2，范 模1，段夢蘭2

(1.中海油研究總院，北京 100027;2.中國石油大學(xué)(北京)，北京 102249)

應(yīng)用數(shù)值模擬與模型試驗(yàn)相結(jié)合的方法研究半潛式生產(chǎn)平臺系泊狀態(tài)下的耦合動力特性。建立耦合分析模型，時域內(nèi)計(jì)算求解平臺的動力響應(yīng)，選取縮尺比為1∶60，采用等效截?cái)嗄Ｐ头椒▽?shù)值模擬結(jié)果進(jìn)行驗(yàn)證。通過對比模型試驗(yàn)與數(shù)值模擬結(jié)果發(fā)現(xiàn):等效截?cái)嘞挡聪到y(tǒng)可以較好地模擬平臺的位移響應(yīng)，但在系泊張力方面卻差異較大，此外極端海況下平臺的甲板上浪問題也必須得到充分重視。

半潛式生產(chǎn)平臺;耦合動力分析;等效截?cái)?模型試驗(yàn);數(shù)值模擬

半潛式生產(chǎn)平臺是一種新型的深水順應(yīng)式生產(chǎn)平臺，由于具有較大的甲板面積和相對較低的建造難度，以及抗風(fēng)浪能力強(qiáng)等特點(diǎn)，是目前數(shù)量最多，更新?lián)Q代最快的一類深水浮式平臺，被廣泛應(yīng)用于世界各地的深水油氣開發(fā)項(xiàng)目［1］。生產(chǎn)平臺長期固定于無遮蔽的深水海域，承受極端的環(huán)境載荷，其運(yùn)動響應(yīng)、纜繩張力以及甲板上浪等技術(shù)指標(biāo)是立管系統(tǒng)安裝乃至整體油田開發(fā)方案制定的重要依據(jù)。一直以來，業(yè)界通常采用數(shù)值模擬與模型試驗(yàn)相結(jié)合的方法來獲得以上數(shù)據(jù)［2］。

在數(shù)值模擬方面，隨著水深的不斷增加，系泊纜繩與平臺主體間的耦合動力特性已達(dá)到必須考慮的程度［3］，必須將平臺主體與系泊系統(tǒng)的動力響應(yīng)作為一個整體求解，其中主體的水動力系數(shù)、一階波浪力以及二階波浪力首先在頻域中計(jì)算，然后將其轉(zhuǎn)換到時域內(nèi)，聯(lián)立求解方程得到耦合系統(tǒng)的整體運(yùn)動響應(yīng)［4-6］。而在模型試驗(yàn)方面，考慮到平臺實(shí)際作業(yè)水深不斷加大和試驗(yàn)水池水深有限的矛盾，國內(nèi)外學(xué)者提出了混合截?cái)嗍侥Ｐ驮囼?yàn)技術(shù)解決這一難題［7-10］，其基本思想［11］是采用數(shù)值模擬方法計(jì)算平臺在實(shí)際海洋環(huán)境下的運(yùn)動響應(yīng)與纜繩張力，在進(jìn)行模型試驗(yàn)之前，根據(jù)給定比例尺將纜繩截?cái)?，其截?cái)嗵幚|繩運(yùn)動與張力根據(jù)計(jì)算結(jié)果進(jìn)行數(shù)值重構(gòu)，從而保證模型實(shí)驗(yàn)結(jié)果能夠較為準(zhǔn)確地模擬實(shí)際情況。等效設(shè)計(jì)的難點(diǎn)在于:1)纜繩垂向與水平向恢復(fù)特性的準(zhǔn)確模擬，即截?cái)嗨钕挡蠢|與水面夾角、纜繩周向角等隨水平位移的變化率將大于全水深系泊纜;2)非均勻材質(zhì)與特殊裝置的使用產(chǎn)生摩擦阻尼難以量化;3)纜繩本身的質(zhì)量動力特性與粘性阻尼的準(zhǔn)確模擬。因此該領(lǐng)域的科研工作也成為海洋工程界的熱點(diǎn)問題之一。

以兩類目前最流行的半潛式生產(chǎn)平臺作為研究對象，結(jié)合南海某油氣田特定的環(huán)境條件，首先采用數(shù)值模擬方法計(jì)算平臺極端載荷下的運(yùn)動響應(yīng)與系泊張力，之后在上海交通大學(xué)海洋深水試驗(yàn)池中，采用被動等效截?cái)嗄Ｐ驮囼?yàn)技術(shù)進(jìn)行平臺模型試驗(yàn)，從而對之前的數(shù)值模擬結(jié)果進(jìn)行驗(yàn)證，分析等效截?cái)嘞挡聪到y(tǒng)與全水深系泊系統(tǒng)的性能差異，最終對于兩類半潛式生產(chǎn)平臺能否與多點(diǎn)系泊系統(tǒng)相配合，在我國南海深水海區(qū)作為井口生產(chǎn)平臺的方案可行性進(jìn)行前期論證。

1 平臺主尺度與設(shè)計(jì)環(huán)境條件

兩座半潛式生產(chǎn)主體采用了最流行的深吃水設(shè)計(jì)理念，具有雙浮箱與雙橫梁船體，立柱部分分別為對心傾斜四立柱式與垂直六立柱式，其主尺度參數(shù)如表1所示。

表1 兩類半潛式平臺主尺度參數(shù)Tab.1 Main scale particulars of two semi platforms

兩座半潛式生產(chǎn)平臺均采用4×4共16根鋼質(zhì)懸鏈線系泊纜進(jìn)行定位，鋼纜直徑為180 mm，破斷張力為21 300 kN;錨鏈直徑為152 mm，破斷張力為20 157 kN，系泊系統(tǒng)設(shè)計(jì)參數(shù)如表2所示。

表2 兩類半潛式平臺系泊系統(tǒng)參數(shù)Tab.2 Main parameters of mooring system of two semi platforms

以四立柱半潛式平臺為例，圖1、2是該平臺多點(diǎn)系泊系統(tǒng)的總體布置圖與耦合計(jì)算模型圖。

圖1 四立柱平臺系泊系統(tǒng)總體布置Fig.1 The general arrangement of mooring system of 4-leg semi platforms

圖2 四立柱平臺耦合計(jì)算模型示意Fig.2 The coupled calculation model of 4-leg semi platforms

半潛式生產(chǎn)平臺的假想作業(yè)海區(qū)選取為南海某油氣田，最大作業(yè)水深為1 500 m，擬懸掛SCR剛性立管，根據(jù)國外工程經(jīng)驗(yàn)，須保證平臺在設(shè)計(jì)環(huán)境條件下的最大水平偏移小于8%的作業(yè)水深。此次試驗(yàn)中對于平臺在該海區(qū)不同重現(xiàn)期的環(huán)境條件下的運(yùn)動特性進(jìn)行了時域模擬，分別包括一年一遇，百年一遇與五百年一遇等，其多點(diǎn)擴(kuò)展式系泊系統(tǒng)的設(shè)計(jì)載荷強(qiáng)度取為百年一遇環(huán)境條件，具體參數(shù)如表3所示。

表3 海洋環(huán)境主要參數(shù)Tab.3 Main parameters of environment condition

2 數(shù)值模擬方法

采用挪威船級社(DNV)出品的SESAM軟件與美國Ultramarine公司出品的MOSES軟件數(shù)值分析半潛式平臺的系泊耦合動力特性，根據(jù)主尺度圖紙建立兩座平臺的水動力計(jì)算模型如圖3、4所示。

圖3 四立柱半潛式平臺水動力模型示意Fig.3 Hydrodynamic model of 4-leg semi platform

圖4 六立柱半潛式平臺水動力模型示意Fig.4 Hydrodynamic model of 6-leg semi platform

在依據(jù)三維勢流理論計(jì)算得到平臺船體的附加質(zhì)量、阻尼系數(shù)及波浪激勵力等數(shù)據(jù)之后，采用耦合時域方法計(jì)算平臺與系泊系統(tǒng)耦合動力響應(yīng)，半潛式平臺船體6自由度時域運(yùn)動方程可表示:

式中:M為平臺船體質(zhì)量矩陣，A(ω)為頻域內(nèi)附加質(zhì)量矩陣，C為頻域內(nèi)勢流阻尼矩陣，D2為平方項(xiàng)粘性阻尼矩陣，K(x)為非線性恢復(fù)剛度矩陣，q為外力矢量，可展開表示為

其中，qwind為風(fēng)力為一階波浪力為二階波浪力，qcur為海流力，qmoor為由系泊系統(tǒng)以及立管系統(tǒng)所產(chǎn)生的耦合作用力。

之后根據(jù)Cummins脈沖理論對于耦合系統(tǒng)進(jìn)行時域求解計(jì)算，其中平臺與系泊系統(tǒng)的運(yùn)動響應(yīng)在每個時間步中聯(lián)合求解，完全考慮了兩者的全耦合動力特性，動力平衡方程可表示為

其中，RI、RD和RS分別表示整體慣性力、整體阻尼力與相互作用力為平臺船體與系泊系統(tǒng)整體的位移、速度與加速度矢量;RE表示外部激勵力矢量，主要包括重浮力、與平臺位移相關(guān)的回復(fù)力，系泊纜繩水動力載荷與特定集中力等。

3 模型試驗(yàn)

模型試驗(yàn)在上海交通大學(xué)海洋工程國家重點(diǎn)實(shí)驗(yàn)室進(jìn)行，其海洋深水試驗(yàn)池最大水深為10 m，配有先進(jìn)的造波、造流、造風(fēng)系統(tǒng)，可以模擬各種實(shí)際海洋環(huán)境條件，通過大面積可升降假底，可根據(jù)試驗(yàn)要求改變水深。本次試驗(yàn)按照1∶60海洋工程常規(guī)比例尺設(shè)計(jì)，其系泊系統(tǒng)由上海交通大學(xué)深水工程試驗(yàn)室進(jìn)行設(shè)計(jì)，系泊等效截?cái)嘣O(shè)計(jì)時主要考慮以下因素:1)系泊系統(tǒng)水平/垂向回復(fù)力特性一致;2)保證系泊系統(tǒng)垂向回復(fù)力/位移特性一致;3)保證具有代表性的系泊纜張力特性一致。以四立柱平臺180°方向?yàn)槔?，圖5、6是該平臺全水深系泊系統(tǒng)截?cái)嗨钕挡聪到y(tǒng)的整體水平回復(fù)力與最大纜繩張力隨平臺水平偏移量的變化對比曲線。

通過圖5與圖6發(fā)現(xiàn)，當(dāng)平臺出現(xiàn)100 m左右的水平偏移時，數(shù)值模擬與模型試驗(yàn)的結(jié)果極為接近，即在小位移情況下，截?cái)嗍较挡聪到y(tǒng)基本上可以模擬全水深條件下系泊系統(tǒng)的整體剛度與靜態(tài)纜繩張力。模型試驗(yàn)包括靜水自由衰減試驗(yàn)與風(fēng)、浪、流組合下的加載試驗(yàn)。在加載試驗(yàn)中，假設(shè)各類環(huán)境載荷入射方向相同，考慮3種典型角度(180°迎浪方向，225°首斜浪方向與270°橫浪方向)，以及不同作業(yè)狀態(tài)，對兩座平臺以下數(shù)據(jù)進(jìn)行了測量和處理:1)平臺6個自由度的運(yùn)動，即縱蕩、橫蕩、垂蕩、縱搖、橫搖及首搖;2)每根纜繩導(dǎo)纜孔處纜繩張力;3)平臺主甲板下緣各處的氣隙高度值。

圖5 四立柱平臺水平回復(fù)剛度曲線Fig.5 Horizontal stiffness curves of 4-leg semi platform

圖6 四立柱平臺纜繩張力變化曲線Fig.6 The maximum tension of 4-leg semi platform

4 數(shù)值模擬與模型試驗(yàn)結(jié)果的對比與分析

以下是兩種方法得到的平臺運(yùn)動響應(yīng)與纜繩張力對比結(jié)果。以四立柱半潛式平臺為例，表4是該平臺不同方向百年一遇環(huán)境載荷作用下的各自由度運(yùn)動響應(yīng)最大值情況;表5為此時最大系泊纜繩張力情況。圖7為四立柱平臺縱蕩、縱搖與纜繩張力譜情況。

表4 四立柱平臺生存狀態(tài)運(yùn)動響應(yīng)最大值情況Tab.4 The maximum motion responses of 4-leg semi platform in survival condition

表5 四立柱平臺系泊系統(tǒng)最大纜繩張力情況Tab.5 The maximum tension of 4-leg semi platform mooring system

通過表4與表5中數(shù)據(jù)結(jié)果進(jìn)行對比并結(jié)合圖7可以發(fā)現(xiàn):

1)平臺6自由度運(yùn)動響應(yīng)方面，兩種方法得到的結(jié)果吻合較好，特別是水平向運(yùn)動響應(yīng)，最大幅值的相對誤差僅在5%以內(nèi)，其中低頻與定常頻率成分占主導(dǎo)地位，數(shù)值模擬與模型試驗(yàn)結(jié)果均小于8%的作業(yè)水深，滿足生產(chǎn)平臺懸掛SCR立管的設(shè)計(jì)要求;垂蕩響應(yīng)方面，運(yùn)動幅值在180°方向達(dá)到最大，響應(yīng)的譜峰周期與波浪譜峰周期接近，實(shí)驗(yàn)值略大于計(jì)算值;而縱橫搖響應(yīng)方面，波頻成分與低頻成分所占比重較為接近，現(xiàn)象表明，極端載荷下平臺可能出現(xiàn)低頻共振搖擺現(xiàn)象。數(shù)值模擬方法得到的計(jì)算結(jié)果可以較為準(zhǔn)確地模擬半潛式平臺惡劣海況下的運(yùn)動響應(yīng)。

2)纜繩張力方面，數(shù)值模擬與模型試驗(yàn)結(jié)果顯示，當(dāng)環(huán)境載荷呈225°入射時，系泊纜繩將出現(xiàn)最大張力，兩種方法得到的結(jié)果相差較大，其中試驗(yàn)結(jié)果比計(jì)算結(jié)果大17%左右，通過對于張力時間歷程的譜分析發(fā)現(xiàn)，兩者在除波頻與定常段符合較好，而在低頻段差異較為明顯，因此導(dǎo)致如此大的誤差，其可能原因:被動截?cái)嗍较挡聪到y(tǒng)基于靜力等效原則進(jìn)行設(shè)計(jì)，系泊纜繩的質(zhì)量動力與粘性阻尼特性與全水深系泊系統(tǒng)相差較大，此外，大比例尺模型試驗(yàn)時流體的粘性作用被放大無疑會加劇這一現(xiàn)象。

圖7 四立柱平臺生存狀態(tài)運(yùn)動響應(yīng)譜與最大系泊張力譜Fig.7 The motion and maximum tension spectra of 4-leg semi platform in survival condition

3)在最小氣隙值方面，采用頻域法計(jì)算得到的平臺最小氣隙值大于3 m，滿足規(guī)范要求，而采用模型試驗(yàn)方法得到的氣隙值小于零，表明此時平臺已經(jīng)出現(xiàn)甲板上浪。通過現(xiàn)場觀測發(fā)現(xiàn)，極端環(huán)境載荷下半潛式生產(chǎn)平臺產(chǎn)生了較大的定常傾角，瞬時波面與平臺各自由度運(yùn)動間也會出現(xiàn)較大的相位差，以及波浪沿立柱表面的局部爬升現(xiàn)象等非線性問題，都使上浪問題變得更為復(fù)雜和難以預(yù)測，因此采用數(shù)值模擬方法計(jì)算甲板上浪較為困難，須結(jié)合模型試驗(yàn)才能得出較為可信的結(jié)果(見圖8)。

圖8 百年一遇環(huán)境條件下兩座平臺的甲板上浪情況Fig.8 The green water condition of two semi platforms in a hundred years condition

表6為兩座半潛式平臺百年一遇環(huán)境載荷下最大水平位移與最大纜繩張力的對比情況。

表6 兩座平臺最大水平偏移與最大纜繩張力對比Tab.6 The comparison of horizontal offsets and mooring line tension between two platforms

通過對比發(fā)現(xiàn)，相對于六立柱平臺中，四立柱平臺水平向的最大偏移量與纜繩張力較小，由此可見，盡量減小平臺的水線面面積，能有效降低平臺所承受的水平方向流體載荷，增大纜繩的安全裕度，有助于提高深水浮式平臺極端環(huán)境下的操作安全性。

5 結(jié)語

以兩座預(yù)想安裝在我國南海，作業(yè)水深為1 500 m的半潛式生產(chǎn)平臺為研究對象，采用數(shù)值模擬與模型試驗(yàn)相結(jié)合方法對平臺的運(yùn)動響應(yīng)特性進(jìn)行對比研究，得到結(jié)論如下:

兩座半潛式生產(chǎn)平臺在該油汽田百年一遇環(huán)境條件下的最大水平偏移量小于8%的作業(yè)水深，滿足設(shè)計(jì)要求，由于四立柱平臺受到的水平環(huán)境載荷更低，纜繩的安全裕度也更大，采用兩種方法得到的結(jié)果吻合校核;當(dāng)極端環(huán)境載荷呈斜45°角入射時，纜繩出現(xiàn)最大系泊張力，是系泊系統(tǒng)的最危險狀態(tài)。截?cái)嗍较挡聪到y(tǒng)難以準(zhǔn)確模擬纜繩的質(zhì)量動力與粘性阻尼特性，導(dǎo)致最大系泊張力的試驗(yàn)結(jié)果與計(jì)算結(jié)果相差較大，此外數(shù)值模擬方法不能準(zhǔn)確評估深水平臺的甲板上浪問題，需結(jié)合模型試驗(yàn)才能得出較為可信的結(jié)果。

［1］李潤培.深海平臺技術(shù)的研究現(xiàn)狀與發(fā)展趨勢［J］.中國海洋平臺，2003，18(3):1-5.

［2］童 波，楊建民，李 欣.深水半潛平臺懸鏈線式系泊系統(tǒng)耦合動力分析［J］.中國海洋平臺，2008，23(6):1-7.

［3］石 山，冉志煌，漆春茂.海洋工程深水開發(fā)浮式系統(tǒng)的耦合計(jì)算方法［J］.中國造船，2005，46(4):47-51.

［4］余 龍，王 娟.半潛式平臺深水錨泊系統(tǒng)三維時域動力分析［J］.中國海洋平臺，2007，22(6):34-37.

［5］Chen X H，Zhang J，Johnson P.Dynamic analysis of mooring lines with inserted spring［J］.Applied Ocean Research，2001，23(5):277-284.

［6］Garrett D L.Coupled analysis of floating production systems［J］.Ocean Engineering，2005，32:802-816.

［7］Stansberg C T，Oritsland O，Kleiven G.Verideep:Reliable methods for laboratory verification of mooring and station keeping in deep water［C］//Proceedings of the Annual Offshore Technology Conference.2000:683-693.

［8］Ormberg H，Baarholm R，Stansberg C T.Time-domain coupled analysis of deepwater TLP，and verification against model tests［C］//Proceedings of the International Offshore and Polar Engineering Conference.2003:145-152.

［9］蘇一華，楊建民，肖龍飛.深海平臺試驗(yàn)中水深截?cái)嘞挡蠢|動力特性對比分析［J］.海洋工程，2009，27(1):17-21.

［10］周杭霞，張火明，高明正，等.智能優(yōu)化平臺在等效水深截?cái)嘞到y(tǒng)設(shè)計(jì)中應(yīng)用［J］.海洋工程，2011，29(1):94-99.

［11］張 威，楊建民，胡志強(qiáng).深水半潛式平臺模型試驗(yàn)與數(shù)值分析［J］，上海交通大學(xué)學(xué)報(bào)，2007，41(9):1430-1434.

Comparative study of coupled dynamic characteristics for deepwater semi-submerged production platforms

ZHAO Jing-rui1，2，F(xiàn)AN Mo1，DUAN Meng-lan2
(1.CNOOC Research Institute，Beijing 100027，China;2.China University of Petroleum，Beijing 102249，China)

The characteristics of coupled dynamic response for the moored semi-submerged production platforms are studied by the model experimental approach and the numerical simulation method in this paper.The coupled analysis models are established and the dynamic responses are obtained in time domain.An equivalent truncated model test is adopted，with the model scale 1:60，and the 6 DOF motion of the two platforms and the mooring line tensions are recorded in order to validate the numerical results.By comparing the full-depth numerical results，it is concluded that the dynamic response of the platform can be modeled quite well，but the difference of the mooring line tension is significant.In addition，enough attention should be paid to the green water problem in rough sea.

semi-submerged production platform;coupled dynamic analysis;equivalent truncated;model test;numerical simulation

P751

A

1005-9865(2012)04-0049-06

2011-11-25

中海油總公司綜合科研課題“南海深水浮式平臺水動力模型試驗(yàn)研究”資助項(xiàng)目

趙晶瑞(1983－)，男，山西大同人，博士后，從事深水浮體系泊定位及設(shè)計(jì)方法研究。E-mail:zhaojr2@cnooc.com.cn