，

(華南理工大學(xué) 土木與交通學(xué)院，廣州 510641)

海洋平臺在遭遇惡劣海況時不能主動避航，并且在進(jìn)行海上油氣作業(yè)時要求其具有低幅運動特性。為了保證平臺作業(yè)的安全性而進(jìn)行平臺定位系統(tǒng)的選擇極為重要，其中錨泊定位系統(tǒng)由于其經(jīng)濟性和便于維護(hù)修理的優(yōu)勢，在各種動力定位系統(tǒng)中仍然處在不可或缺的地位。

對于海洋平臺的水動力性能及其錨泊系統(tǒng)的設(shè)計模擬分析，集中在比較研究Morison方程和勢流理論兩種方法，并對比實驗結(jié)果，得到這2種方法除了在低頻范圍相差較大外，在主要頻率內(nèi)符合良好的結(jié)論[1]；對半潛式平臺用勢流理論并結(jié)合Morison方程進(jìn)行粘性力修正的方法，不僅在頻域和時域內(nèi)計算結(jié)果與試驗結(jié)果符合良好，且考慮非線性波后時域結(jié)果與試驗結(jié)果也有很好的符合[2]；運用Rankine源的邊界元法對半潛式海洋平臺的波浪繞射問題進(jìn)行了時域模擬[3]；采用靜力方法對懸鏈線式系泊和張緊式系泊的纜繩張力進(jìn)行比較，發(fā)現(xiàn)兩者隨水深的變化規(guī)律大不相同[4]，利用全耦合時域數(shù)值模型來研究深水系泊鏈的阻力系數(shù)[5]；進(jìn)行了風(fēng)浪流同時作用下的 FPSO 與轉(zhuǎn)塔式系泊系統(tǒng)以及立管系統(tǒng)的耦合模型的動力分析方法及其模型試驗的研究[6]；用管線的有限元模擬方式模擬系泊纜，將管線與浮體的連接看作是節(jié)點，用頻域計算方法分析平臺主體及其錨泊系統(tǒng)在一階波浪力作用下的運動響應(yīng)[7]；建立考慮平臺主體、錨泊系統(tǒng)和立管系統(tǒng)相互作用的半潛式平臺耦合分析數(shù)值模型[8]；分析懸鏈線系泊系統(tǒng)的錨纜直徑、長度與預(yù)張力角度等因素對系泊纜的動力特性的影響[9]；研究系泊系統(tǒng)的一般設(shè)計方法并基于彈性桿理論的非線性有限元法對繃緊索系泊系統(tǒng)進(jìn)行動力分析[10]；采用集中質(zhì)量法建立系泊纜索的力學(xué)模型，得出了浮體和海流流速對纜繩受力變形影響較大，海底錨鏈長度至少要保持200 m左右的長度的結(jié)論[11]。

在前人對水動力及系泊系統(tǒng)分析研究的基礎(chǔ)上，進(jìn)行目標(biāo)半潛平臺典型工況水動力分析和深水張緊式系泊方案的設(shè)計，并在相同預(yù)張力情況下就該系泊方案與懸鏈線系泊進(jìn)行對比分析。

1 計算理論

1.1 平臺頻域運動響應(yīng)分析

在本半潛平臺水動力計算中，波浪采用Airy線性波，平臺的運動響應(yīng)以響應(yīng)函數(shù)RAO的形式來表達(dá)，平臺在頻域中的響應(yīng)譜SR(ω)可表示為

(1)

根據(jù)DNV- RP規(guī)范規(guī)定，平臺運動的最大響應(yīng)幅值可以用以下公式來進(jìn)行預(yù)報。

(2)

1.2 平臺時域非線性響應(yīng)分析

在半潛平臺的時域分析中取用了平臺與系泊系統(tǒng)的耦合運動模型，考慮風(fēng)、浪、流聯(lián)合作用下整體結(jié)構(gòu)的非線性耦合響應(yīng)。平臺及其系泊系統(tǒng)耦合運動方程可表示為

(3)

式中：M為平臺的質(zhì)量矩陣；C為平臺的阻尼矩陣；K為平臺的剛度矩陣；Fst+Fwf+Fsd+Fm分別為靜態(tài)力、波頻力、低頻慢漂力以及系泊力。靜態(tài)力包括靜態(tài)風(fēng)力、靜態(tài)流力及平均二階漂移力，其中平均漂移力與波高的平方成正比。通常平臺波浪漂移力的計算方法分為遠(yuǎn)場解法和近場解法。對比只能求解3個水平分量(縱蕩、橫蕩和艏搖)的遠(yuǎn)場解法來說，近場解法采用在物體濕表面直接進(jìn)行壓力積分的計算方法，只要節(jié)點足夠多就可以對6個分量上獲得令人滿意的結(jié)果。近場法更加適用于復(fù)雜形狀的海洋結(jié)構(gòu)物的波浪漂移力的計算。

應(yīng)用J.A. Pinkster所提出的近場法計算結(jié)構(gòu)物波浪漂移力F(2)和力矩M(2)。

?S0ρ(Φ″wt+Φ″dt)nods

(4)

(5)

式中：下標(biāo)t為對時間求偏導(dǎo)；n0為平臺內(nèi)法線方向；S0為平臺濕面積；Γ0為水線；ρ為流體密度；ζr為相對波面升高；r0為平臺表面各點的坐標(biāo)；α為平臺角位移；P0為平臺上各點的位移；Φ為速度勢；Φw為波浪入射勢；Φd為波浪繞射勢。

采用Jonswap譜進(jìn)行不規(guī)則波的模擬計算，Jonswap譜譜密度表達(dá)式為

(6)

式中：H1/3為有義波高，T01為波浪的特征周期，ω為波浪圓頻率。

2 半潛平臺參數(shù)與系泊系統(tǒng)設(shè)計

2.1 平臺主要參數(shù)

以某半潛平臺為研究對象,平臺的主要結(jié)構(gòu)包括2個浮體、4個立柱、4個橫撐，以及主甲板和井架，平臺的主要參數(shù)見表1。

表1 平臺主要參數(shù) m

2.2 系泊系統(tǒng)設(shè)計

系泊系統(tǒng)需要保證平臺能夠在百年一遇的極端海況下所受到的最大張力小于其破斷力，并且需要滿足一定的安全儲備。對于半潛式平臺來說，水深小于300 m的淺水海況下采用張緊式系泊方式的平臺僅為21%，而水深在300 m到1 500 m之間的增加到40%。當(dāng)水深增加到1 500 m以上超深水海域，就有超過85%的半潛平臺采用張緊式系泊。張緊式系泊系統(tǒng)采用高強度的錨鏈和聚酯材料，選用亞星R5級錨鏈來進(jìn)行錨泊系統(tǒng)的頂端和底端設(shè)計。頂端考慮到要求中間段聚酯纜的上端要遠(yuǎn)離導(dǎo)纜器，至少保持在水面100 m以下；而海底段要防止海底沙塵等對聚酯纜造成磨損，用錨鏈來連接海底錨固系統(tǒng)且長度控制在400 m以內(nèi)。此外，為了防止系泊纜出現(xiàn)大的磨損，將系泊纜的預(yù)張力控制在纜索破斷強度的30%左右。

考慮對工作水深為1 500 m的某半潛平臺設(shè)計一套合適的錨泊系統(tǒng)，建立平臺以及系泊系統(tǒng)的耦合模型，以水動力性能分析計算為基礎(chǔ)，在時域內(nèi)研究平臺和系泊系統(tǒng)在風(fēng)、浪、流聯(lián)合作用下的運動響應(yīng)的模擬計算，并就結(jié)果進(jìn)行規(guī)范校核。該半潛平臺系泊系統(tǒng)設(shè)計采用8點系泊，根據(jù)API-RP-2SM規(guī)范[12]，系泊系統(tǒng)均采用3段式對稱布置，懸鏈線系泊采用鋼鏈(100 m)-鋼纜(1 700 m)-鋼鏈(450 m)，張緊式系泊采用鋼鏈(170 m)-聚酯纜(1 270 m)-鋼鏈(250 m)布置。纜繩布置角度為30°，夾角為30°，懸鏈線式系泊系統(tǒng)的水平輻射距離為2 660 m，張緊式系泊系統(tǒng)的水平輻射距離為1 650 m，按照統(tǒng)一的建造標(biāo)準(zhǔn)來選擇鋼鏈和聚酯纜的具體屬性。系泊纜的詳細(xì)參數(shù)如表2。

表2 系泊纜參數(shù)

將平臺的實際模型簡化后所建立的平臺主體的數(shù)值模型如圖1，系泊纜的布置方式如圖2。

圖1 半潛平臺水動力模型

圖2 系泊纜布置情況

3 平臺運動模擬分析

平臺工作水深為1 500 m，自存狀態(tài)下吃水為19 m，頻域分析中搜索浪向間隔為15°，環(huán)境載荷取用中國南海海域百年一遇海況下的波浪載荷、風(fēng)載荷及流載荷。選取JONSWAP波浪譜來模擬不規(guī)則波浪作用，有義波高為12.9 m，譜峰周期為16.1 s；風(fēng)譜選用API風(fēng)譜，參考高度為10 m，平均風(fēng)速為40.9 m/s；流載荷為剪切流，海平面流速為1.9 m/s。

3.1 運動響應(yīng)幅值算子RAO計算

計算發(fā)現(xiàn)半潛平臺艏搖運動的響應(yīng)幅值的量級要比其他方向的響應(yīng)小得多，所以本文主要考慮縱蕩、橫蕩、縱搖、橫搖及垂蕩5個自由度的運動響應(yīng)。平臺在不同浪向下，縱蕩、橫蕩、縱搖、橫搖以及垂蕩的RAO的響應(yīng)見圖3。由于半潛平臺的對稱性，這里只給出0°～90°浪向下的響應(yīng)。

圖3 RAO的響應(yīng)曲線

從計算結(jié)果可知，頻率響應(yīng)函數(shù)由于平臺的對稱性，縱蕩和橫蕩、縱搖和橫搖的響應(yīng)趨勢是相似的，平臺在低頻部分的響應(yīng)較大。在波浪順浪(0°)或者迎浪(180°)入射時，平臺的縱蕩、縱搖以及垂蕩為主要的運動模態(tài)；在90°橫浪入射時，平臺的橫蕩、橫搖和垂蕩為主要的運動模態(tài)；由橫蕩和縱蕩響應(yīng)曲線可以看出，平臺的震蕩運動在二階低頻波浪慢漂力和風(fēng)力的作用下的低頻運動響應(yīng)出現(xiàn)顯著波動，其低頻運動響應(yīng)遠(yuǎn)大于波頻運動響應(yīng)，此時其固有周期和常規(guī)波浪周期相接近，就會發(fā)生共振，并且在頻域計算中未考慮系泊系統(tǒng)對平臺運動的回復(fù)力的影響，這對平臺的橫蕩和縱蕩的計算結(jié)果有一定的影響；平臺縱搖和橫搖在0.53 rad/s附近出現(xiàn)響應(yīng)極值，平臺出現(xiàn)較大的波頻運動，這需要在實際設(shè)計中引起重視；垂蕩運動在0.27 rad/s附近出現(xiàn)RAO極值，由于平臺垂蕩固有周期為23 s左右，而常規(guī)波浪的周期一般為6～14 s，故平臺的垂蕩響應(yīng)的固有周期在常規(guī)波浪周期范圍之外，可以說該平臺具備了較好的垂蕩性能，能避免發(fā)生共振。

3.2 錨泊系統(tǒng)響應(yīng)分析

為使半潛平臺能夠安全地進(jìn)行鉆井作業(yè)，須保證平臺的水平活動小于作業(yè)水深的3%。根據(jù)API-RP-2SM規(guī)范[12]規(guī)定，采用耦合動力計算時，對于完好無損系泊系統(tǒng)，在計算半潛平臺遭遇百年一遇的海況條件的動力分析時，系泊系統(tǒng)最大張力的安全因子取1.67，聚酯纜最大張力的安全因子取1.82；對于一根系泊纜斷裂的情況，系泊系統(tǒng)最大張力的安全因子取1.25，聚酯纜最大張力的安全因子取1.43，其中安全因子定義為系泊纜破斷強度與最大張力的比值。

為了證實張緊式系泊系統(tǒng)的優(yōu)越性，在保持預(yù)張力不變的情況下分別設(shè)計了懸鏈線式系泊系統(tǒng)和張緊式系泊系統(tǒng)。針對風(fēng)浪流同向，橫浪90°入射的海況，模擬中國南海百年一遇海況，研究設(shè)計系泊系統(tǒng)與平臺耦合在該海況下實時動力響應(yīng)，包括平臺的位移響應(yīng)、受力響應(yīng)、錨鏈的張力響應(yīng)等(見圖4)。下面選取耦合計算結(jié)果的時歷響應(yīng)分析風(fēng)浪流同向90°入射時，平臺橫蕩、橫搖、垂蕩為運動的主要形態(tài)。

圖4 半潛平臺系泊纜最大張力

設(shè)計工況下2種系泊形式運動響應(yīng)與最大系泊張力計算結(jié)果見表3。

由表3可見，在保持2種系泊方式預(yù)張力都為2 400 kN的情況下，懸鏈線式系泊系統(tǒng)的安全系數(shù)要比張緊式系泊高得多，完全符合前述的API規(guī)范要求，但是懸鏈線式系泊系統(tǒng)橫蕩運動的位移最高接近75.1 m，遠(yuǎn)遠(yuǎn)超過了設(shè)計前要求的3%水深。不僅如此，懸鏈線系泊系統(tǒng)的水平輻射范圍與水深的比值高達(dá)1.7左右，而對于該比值張緊式系泊系統(tǒng)僅為1.1，張緊式系泊聚酯纜的安全系數(shù)為1.94，明顯大于API規(guī)范要求的1.82，符合要求。綜上所述，對于深水鉆井平臺來說，采用張緊式系泊比起懸鏈線式系泊不僅可有效減小系泊系統(tǒng)的占地面積，還極大程度地減小了系泊系統(tǒng)在水下和其他結(jié)構(gòu)物相互影響的危險性，最終采用張緊式為本項目的系泊方案。圖3a)、圖3b)為張緊式系泊在環(huán)境力90°入射時的平臺主要運動模態(tài)位移以及最大系泊纜張力值。

表3 設(shè)計工況下兩種系泊形式響應(yīng)結(jié)果比較

計算表明，在一定范圍內(nèi)改變系泊鏈的長度、預(yù)張力或剛度等參數(shù)，對平臺的橫蕩位移影響較大，而對垂蕩和橫搖運動影響甚微。

4 結(jié)論

1)該平臺的水動力響應(yīng)曲線的規(guī)律與實際相符。90°浪向下，橫蕩、橫搖為平臺的主要運動模態(tài)；0°和180°浪向下，縱蕩、縱搖為主要的運動模態(tài)；而平臺垂蕩運動的RAO響應(yīng)則受浪向的影響不大。

2)深水環(huán)境中，相同預(yù)張力情況下，張緊式與懸鏈線系泊方式比較，雖然張緊式系泊系統(tǒng)的平臺其系泊纜索的張力響應(yīng)的安全系數(shù)大于懸鏈線式，但其運動響應(yīng)幅值要比懸鏈線式小得多。所以，深水環(huán)境中在符合API-RP-SM規(guī)范的條件下，選取張緊式系泊平臺的定位能力具備更大優(yōu)勢。

3)張緊式系泊由于系泊纜中間段的聚酯纜變形的非線性太強，導(dǎo)致其受力分析存在較大的困難，故對懸鏈線和張緊式兩種系泊方式的更為詳細(xì)全面的受力比較分析，還有待今后的深入研究。

[1] So Y lemez M, Atlar M.A. Comparative Study of Two Practical Methods for Estimating the Hydrodynamic Loads and Motions of a Semi-Submersible[J] .Journal of Offshore Mechanics and Arctic Engineering, 2000, 122:57- 65.

[2] CLAUSS G F, SCHMITTNER C E, FREAK S K. Wave impact on semisubmersibles time- domain analysis of motions and forces [A]. Proceedings of theThirteenth(2003)International Offshore and Polar Engineering Conference，2003, JSC- 371.

[3] 徐剛，段文洋.規(guī)則波中半潛式平臺水動力的時域數(shù)值模擬[J].江蘇科技大學(xué)學(xué)報(自然科學(xué)版)，2012,26(1)：1- 7.

[4] 陳新權(quán)，譚家華.深海半潛式平臺系泊系統(tǒng)設(shè)計研究[J].中國海洋平臺,2011,26(2):43- 48.

[5] WICHERS J E W,VOOGT H J,ROELOFS H W, e5a8. Deep star- CTR 4401 benchmark model test[R].Technical Report No.16417- 1- OB,MARIN, Netherlands，2001.

[6] ARCANDRA Tahar,KIM M H.Hull/mooring/riser coupled dynamic analysis and sensitivity study of a tanker- based FPSO[J].Applied Ocean Research,2003,24:367- 382.

[7] Le Cunff,Ryu ,S.Duggal,A.Ricbourg,C.Heurtier,J.M.Derivation of CALM Buoy Coupled Motion RAOs in Frequency Domain and Experimental Validation[C].ISOPE 2007- JSC- 594.

[8] ROVERI F E.Coupled motion analysis of a semi- submersible platform incamposbasin[C].23rd 74International Conference on Offshore Mechanics and Arctic Engineering . Vancouver, British Columbia, Canada,OMAE 2004- 51120.

[9] 牛淑貞，潘新民，楊洪平，等．新一代天氣雷達(dá)災(zāi)害性天氣警報和臨近預(yù)報系統(tǒng)應(yīng)用[J].氣象，2009，35(5):105- 111.

[10] 由際昆，王言英.深水半潛平臺繃緊索系泊系統(tǒng)設(shè)計研究[J]中國海洋平臺,2009,24(1):24- 30.

[11] 趙晶瑞，范模，段夢蘭.深水系泊纜繩動力響應(yīng)研究[J].船舶工程,2013,35(4):103- 107.

[12] American Petroleum Institute. Design and analysis of station- keeping systems for floating structures[S].API recommended practiced 2SK,2005.