梁 寧， 黃維平， 周 陽， 常 爽， 孫遠(yuǎn)坤

(1.中國海洋大學(xué)， 山東 青島 266100； 2.國家海洋局第二海洋研究所， 杭州 310012)

?

基于抵消周期設(shè)計(jì)的干樹半潛式平臺水動(dòng)力性能研究

梁 寧1， 黃維平1， 周 陽2， 常 爽1， 孫遠(yuǎn)坤1

(1.中國海洋大學(xué)， 山東 青島 266100； 2.國家海洋局第二海洋研究所， 杭州 310012)

基于抵消周期設(shè)計(jì)了干樹半潛式平臺，并對其水動(dòng)力性能進(jìn)行了研究。深吃水式干樹半潛式平臺不采用垂蕩板等附屬結(jié)構(gòu)，將干式開發(fā)系統(tǒng)和半潛式平臺的優(yōu)點(diǎn)相結(jié)合。通過對結(jié)構(gòu)主尺度的優(yōu)化來降低平臺的垂蕩響應(yīng)，使其能夠支持干采油樹系統(tǒng)。采用三維勢流理論，運(yùn)用SESAM軟件對干樹半潛式平臺進(jìn)行了水動(dòng)力性能分析。運(yùn)用HydroD模塊對平臺進(jìn)行頻域分析，得到平臺附加質(zhì)量、勢流阻尼、幅值響應(yīng)函數(shù)等參數(shù)。運(yùn)用DeepC模塊對平臺進(jìn)行時(shí)域耦合分析，得到平臺加系泊系統(tǒng)后的時(shí)域響應(yīng)。

干樹半潛式平臺；水動(dòng)力性能；結(jié)構(gòu)設(shè)計(jì)

0 引言

目前,在深海油氣資源開發(fā)中，采用干采油樹的浮式平臺主要有張力腿平臺(TLP)和單柱式平臺(SPAR)[1，2]。在水深超過1 500 m后,張力腿平臺張力筋腱的設(shè)計(jì)將面臨極大的挑戰(zhàn)，其造價(jià)也將大幅度增加，進(jìn)而造成整體開發(fā)方案的經(jīng)濟(jì)性降低。并且張力腿平臺需要海上施工組裝，因此海上施工成本和風(fēng)險(xiǎn)性是張力腿平臺的另一個(gè)制約點(diǎn)。單柱式平臺能夠支撐的甲板面積有限，其超長的立柱決定了它必須要經(jīng)過拖航后主船體的扶正和上部組塊的海上組裝[3，4]。與上述兩種平臺相比，干樹半潛式平臺能夠提供更大的甲板面積、支撐更大的上部荷載，能夠?qū)崿F(xiàn)上部組塊的碼頭安裝進(jìn)而整體拖航。因此干樹半潛式平臺成為深海油氣資源開發(fā)方案中的一個(gè)新的選擇[5]。

半潛式平臺作為一種典型的深海浮式平臺已經(jīng)歷過七代共四十余年的發(fā)展歷程[6]。它具有甲板可變載荷更大、上部組塊碼頭安裝、平臺整體拖航無需海上施工等優(yōu)點(diǎn)[7]。干采油樹位于平臺井口甲板上，通過張緊器與頂張式立管相連接，具有便于控制、作業(yè)維護(hù)費(fèi)用低、生產(chǎn)效率高等特點(diǎn)。傳統(tǒng)半潛式平臺垂蕩固有周期與極端海況的波浪周期接近，導(dǎo)致平臺垂蕩運(yùn)動(dòng)較大，無法支持干采油樹。干樹半潛式平臺通過增加吃水、結(jié)構(gòu)優(yōu)化等方式改善平臺垂蕩運(yùn)動(dòng)性能，使其可以將半潛式平臺和干采油樹的優(yōu)點(diǎn)結(jié)合起來。

1 干樹半潛式平臺設(shè)計(jì)

1.1 抵消周期

通過對主尺度的選取，可以在垂向荷載的抵消周期中降低平臺的垂向運(yùn)動(dòng)響應(yīng)。當(dāng)作用在穿越水面立柱上的Froude-Krylov力和作用在全浸沒浮箱上的繞射力相等時(shí)，平臺所受的垂向荷載抵消進(jìn)而造成平臺垂向運(yùn)動(dòng)最小[8]。

垂向向上的Froude-Krylov力在波浪峰值處的大小可以近似估算為：

(1)

式中：A為立柱水線面面積；D為吃水深度(取正值)；H為波高；T為波浪周期。

垂向向下的浮箱繞射力在波浪峰值處的大小可以近似估算為：

(2)

式中：V為浮箱排水體積；z為浮箱重心高度(取負(fù)值)；CM為附加質(zhì)量系數(shù)。

對于確定的波高，作用在立柱上的垂向Froude-Krylov力隨著波浪周期的提高而增大。與它作用方向相反的浮箱的繞射力隨著波浪周期的提高呈現(xiàn)先增大后減小的趨勢[8]。通過改變干樹半潛式平臺的水線面面積與浮箱排水體積比，可以選定使抵消周期發(fā)生在某個(gè)特定周期。對于生產(chǎn)平臺來講，由于極端環(huán)境條件下立管系統(tǒng)還是與平臺主體保持連接，因此抵消周期的選擇應(yīng)當(dāng)偏向于極端環(huán)境條件下盡可能地減小立管沖程。對于該文的干樹半潛式平臺，針對南海百年一遇極端環(huán)境條件下的譜峰周期為16.1 s。由于平臺是對稱結(jié)構(gòu)，通過改變單一立柱水線面面積和浮箱排水體積比即可選定平臺抵消周期。最終選定結(jié)構(gòu)抵消周期為21 s，即圖1中兩條曲線交點(diǎn)的橫坐標(biāo)值。

圖1 選定平臺的抵消周期

圖2 干樹半潛式平臺示意圖

1.2 干樹半潛式平臺構(gòu)造

該文的干樹半潛式平臺具有深吃水的特點(diǎn)，從而減小了浮箱上的波浪激勵(lì)力[9]。平臺由上部組塊、主船體、系泊系統(tǒng)、立管系統(tǒng)組成，如圖2所示。其主船體包括四個(gè)方形立柱和四個(gè)矩形浮箱，采用16纜張緊式系泊。在立柱選取上需要考慮截面形式、水線面面積、立柱間距等因素。方形截面具有便于和浮箱連接、建造方便等優(yōu)勢，且通過對方形立柱截面倒角的方式使其能夠達(dá)到圓形立柱截面的水動(dòng)力性能。增加立柱水線面可以增加恢復(fù)力矩，進(jìn)而改善平臺穩(wěn)性。但過大的水線面面積會(huì)降低平臺垂蕩運(yùn)動(dòng)固有周期和抵消周期，不利于改善平臺的垂蕩運(yùn)動(dòng)性能。增大立柱間距會(huì)使浮箱與立柱排水比增大，這會(huì)直接提高平臺的抵消周期，但過大的立柱間距會(huì)對甲板設(shè)計(jì)造成極大的挑戰(zhàn)。浮箱主要提供平臺的壓載，降低平臺重心，增大浮箱體積會(huì)提高平臺抵消周期[10]。

1.3 干樹半潛式平臺尺寸

平臺主尺度在滿足半潛式平臺要求的前提下，基于平臺的抵消周期對干樹半潛式平臺主尺度進(jìn)行設(shè)計(jì)、調(diào)整。當(dāng)作用在浮箱上的波浪力同立柱上的波浪力相等時(shí)，平臺的垂蕩運(yùn)動(dòng)最小。在百年一遇極端環(huán)境條件下，平臺的抵消周期與波浪周期接近，此時(shí)，可以使平臺在極端環(huán)境條件下垂蕩運(yùn)動(dòng)幅值較小進(jìn)而保護(hù)頂張式立管系統(tǒng)不被破壞[8，11]。綜合考慮浮箱與立柱排水比、總排水量、上部甲板等因素，選定平臺的主要參數(shù)見表1。

表1 干樹半潛式平臺主要參數(shù)

2 計(jì)算理論

2.1 三維勢流理論

在SESAM中采用三維勢流理論計(jì)算結(jié)構(gòu)物濕表面上的波浪力。勢流理論假定流體為理想流體，即假設(shè)流體具有均勻、無粘性、無旋、不可壓縮的特性。速度勢滿足拉普拉斯方程和海底條件、自由表面條件、物體濕表面條件和輻射條件。根據(jù)拉普拉斯方程以及邊界條件可以唯一的確定出速度勢，然后由伯努利方程計(jì)算出物體濕表面上的壓力，再通過對物體表面進(jìn)行壓力積分就可求出物體所受到的波浪載荷[12，13]。

流場中的總速度勢包含入射勢φI，繞射勢φD和輻射勢φR，則總速度勢函數(shù)為：

(3)

無限水深條件的入射勢為：

(4)

式中：g為重力加速度；A為入射波振幅；ω為入射波頻率；k為入射波波數(shù)；α為入射波方向。

2.2 頻域運(yùn)動(dòng)控制方程

平臺在整體坐標(biāo)系中考慮所受的波浪激勵(lì)以及附加質(zhì)量、阻尼等作用的運(yùn)動(dòng)方程為[14]：

(5)

式中：M為質(zhì)量和慣量矩陣；A(ω)為結(jié)構(gòu)的附加質(zhì)量矩陣；B(ω)p為頻域下的勢流阻尼矩陣；Bv為頻域下的線性化后的粘性阻尼矩陣；C為靜水回復(fù)力矩陣；Ce為外部回復(fù)力矩陣；F(ω,β)為外部激勵(lì)力矩陣。

2.3 時(shí)域運(yùn)動(dòng)控制方程

海洋結(jié)構(gòu)物在實(shí)際的海洋環(huán)境中會(huì)受到波浪載荷、風(fēng)載荷和海流載荷的聯(lián)合作用，故需要較全面的將各因素納入考慮范圍?？紤]風(fēng)、浪、流，以及系泊系統(tǒng)作用的平臺運(yùn)動(dòng)時(shí)域方程可表達(dá)為：

(6)

式中：Fwa(1)(t)為一階波浪載荷；Fwa(2)(t)為二階波浪載荷；Fwi(t)為風(fēng)載荷；Fcu(t)為流載荷；Fmo(t)為系泊系統(tǒng)對平臺的作用力；mij為平臺的質(zhì)量矩陣；uij為附加質(zhì)量矩陣；Cij為恢復(fù)力系數(shù)矩陣；Lij(t-τ)為平臺系統(tǒng)的延遲函數(shù)。

平臺在非線性波浪作用下的瞬時(shí)波浪力可以寫為：

(7)

(8)

(9)

(10)

(11)

圖3 干樹半潛式平臺模型

3 頻域水動(dòng)力分析

設(shè)定干樹半潛式平臺的服役海域?yàn)橹袊虾?，工作水? 500 m。頻域分析時(shí)取規(guī)則波波浪周期為5 s～50 s，間隔1 s。波浪入射方向從0°～90°，間隔15°。利用HydroD中的Wadam求解干樹半潛式平臺的頻域響應(yīng)。對于在海洋環(huán)境中的浮體，運(yùn)動(dòng)時(shí)其所受到的阻尼包含浮體的輻射阻尼、慢漂阻尼、粘性阻尼等。在實(shí)際中難以對上述懲罰逐一進(jìn)行準(zhǔn)確計(jì)算，因此在工程應(yīng)用中采用經(jīng)驗(yàn)系數(shù)的方式對阻尼進(jìn)行估算[15]。在干樹半潛式平臺的計(jì)算中，橫搖、縱搖線性阻尼系數(shù)取臨界阻尼的3%，垂蕩線性阻尼系數(shù)取臨界阻尼的5%。

3.1 平臺模型的建立

根據(jù)干樹半潛式平臺的主尺度，運(yùn)用SESAM軟件中的GeniE模塊建立干樹半潛式平臺的濕表面模型，并將其導(dǎo)入到HydroD模塊建立水動(dòng)力頻域計(jì)算所用的模型，如圖3所示。

3.2 附加質(zhì)量和勢流阻尼

平臺的附加質(zhì)量和勢流阻尼都是矩陣形式，取六個(gè)主對角線上的數(shù)值研究其變化規(guī)律。圖4、圖5分別給出了附加質(zhì)量和勢流阻尼隨波浪周期的變化規(guī)律?？傮w來看，結(jié)構(gòu)轉(zhuǎn)動(dòng)的附加質(zhì)量和勢流阻尼都遠(yuǎn)大于平動(dòng)的附加質(zhì)量和勢流阻尼。縱蕩和橫蕩附加質(zhì)量隨著波浪周期先急劇增大后減小，再增大后減小，最終趨于穩(wěn)定。另外四個(gè)自由度的附加質(zhì)量先減小后增大，而后減小并最終趨于穩(wěn)定?？v蕩和橫蕩的勢流阻尼隨著波浪周期增加先有一個(gè)大波峰，再有一個(gè)小波峰最終趨于零。垂蕩、橫搖、縱搖的勢流阻尼先有一個(gè)小波峰再有一個(gè)大波峰，最終進(jìn)而趨于零，艏搖在短周期內(nèi)有一個(gè)波峰后迅速減小為零。

圖4 平臺附加質(zhì)量

3.3 幅值響應(yīng)函數(shù)

平臺在6個(gè)自由度的運(yùn)動(dòng)可分為平面內(nèi)運(yùn)動(dòng)(縱蕩、橫蕩、艏搖)和平面外運(yùn)動(dòng)(橫搖、縱搖、垂蕩)。圖6～圖11分別給出了干樹半潛式平臺在6個(gè)自由度的運(yùn)動(dòng)響應(yīng)幅值算子(RAO)。結(jié)果表明，平臺垂蕩運(yùn)動(dòng)抵消周期為21 s，固有周期為24 s。平臺抵消周期在百年一遇極端環(huán)境條件譜峰高能量區(qū)域末端，在抵消周期處平臺垂蕩運(yùn)動(dòng)響應(yīng)最小。結(jié)構(gòu)垂蕩運(yùn)動(dòng)固有周期避開主要波浪周期；平臺橫搖、縱搖固有周期在32 s左右，其峰值也成功避開了主要波浪周期。由于頻域分析時(shí)沒有添加系泊系統(tǒng)，因此縱蕩、橫蕩兩個(gè)方向上沒有回復(fù)力，平臺運(yùn)動(dòng)隨波浪周期增大而變大。結(jié)構(gòu)艏搖在波浪周期8s時(shí)有峰值，隨著波浪周期增大艏搖運(yùn)動(dòng)逐漸減小。

圖5 平臺勢流阻尼

圖6 平臺縱蕩RAO

圖7 平臺橫蕩RAO

圖8 平臺垂蕩RAO

圖9 平臺橫搖RAO

圖10 平臺縱搖RAO

圖11 平臺艏搖RAO

4 時(shí)域水動(dòng)力分析

4.1 系泊系統(tǒng)特性和環(huán)境條件

干樹半潛式平臺采用張緊式系泊，該系泊系統(tǒng)由16根系泊纜組成，分為4組對稱布置。第一象限內(nèi)系泊纜與X軸的夾角分別為30°、40°、50°、60°，系泊纜與海底夾角為40°。每根系泊纜由錨鏈、尼龍纜、錨鏈三段組合而成，系泊纜一端與海底錨基相連另一端通過導(dǎo)纜器與張緊器相連。系泊纜參數(shù)見表2，根據(jù)南海海域海況，選取百年一遇環(huán)境條件，見表3。時(shí)域分析中保守考慮，設(shè)定風(fēng)、浪、流是同向的。

表2 系泊纜屬性

表3 百年一遇環(huán)境條件

4.2 平臺時(shí)域響應(yīng)

圖12～圖15依次給出了干樹半潛式平臺在時(shí)域分析中縱蕩、垂蕩、橫搖和艏搖方向上的響應(yīng)?？梢钥闯觯脚_縱蕩最大值為27.2 m，垂蕩最大幅值為3.7 m，橫搖運(yùn)動(dòng)最大角度為3.3°，艏搖運(yùn)動(dòng)最大角度為6.8°。干樹半潛式平臺的核心控制項(xiàng)為垂蕩運(yùn)動(dòng)，若垂蕩運(yùn)動(dòng)幅值過大，則平臺無法支持頂張式立管，進(jìn)而不能采用干采油樹系統(tǒng)。計(jì)算得到的最大垂蕩響應(yīng)幅值為3.7 m，在頂張式立管張緊器的行程范圍內(nèi)，因此可以使用頂張式立管。

圖12 縱蕩響應(yīng)

圖13 垂蕩響應(yīng)

圖14 橫搖響應(yīng)

圖15 艏搖響應(yīng)

圖16給出了干樹半潛平臺垂蕩運(yùn)動(dòng)的頻譜分析結(jié)果?？梢钥闯?，平臺的垂蕩運(yùn)動(dòng)在南海百年一遇環(huán)境條件下有兩個(gè)峰值，第一個(gè)頻率較低的峰值在0.26 rad/s，換算周期為24 s，即為平臺垂蕩運(yùn)動(dòng)的固有周期，與頻域分析結(jié)果一致。第二個(gè)頻率較高的峰值在0.39 rad/s，換算周期為16 s，即為南海百年一遇環(huán)境條件下波浪的譜峰周期。在百年一遇環(huán)境條件下，波浪譜能量范圍內(nèi)造成了較大的平臺垂蕩運(yùn)動(dòng)，并且激發(fā)了平臺在固有周期的運(yùn)動(dòng)。圖17給出了平臺系泊纜中受到最大張力的12號系泊纜的張力時(shí)程曲線，可以得到在百年一遇極端環(huán)境條件下，平臺系泊纜(尼龍纜部分)受到的最大張力為8 120 kN，安全系數(shù)為1.97。系泊纜的安全系數(shù)滿足API-RP-2SM對尼龍纜完整狀態(tài)下最小安全系數(shù)為1.84的要求[16]。

圖17 12號纜繩張力結(jié)果

圖16 平臺垂蕩頻譜結(jié)果

5 結(jié)語

針對抵消周期和設(shè)計(jì)要求進(jìn)行干樹半潛式平臺的主尺度選型。在選定合理的吃水、立柱、浮箱等參數(shù)后平臺能夠滿足干采油樹生產(chǎn)時(shí)對垂蕩性能的要求。在頻域分析中只針對平臺主體進(jìn)行水動(dòng)力分析，得到了平臺六個(gè)自由度上的附加質(zhì)量、勢流阻尼和RAO，初步驗(yàn)證干樹半潛式平臺的運(yùn)動(dòng)性能。針對平臺進(jìn)行了系泊系統(tǒng)設(shè)計(jì)，在時(shí)域分析中計(jì)算平臺和系泊纜的耦合運(yùn)動(dòng)響應(yīng)。從時(shí)域結(jié)果中可以驗(yàn)證平臺具有良好的水動(dòng)力性能，其核心控制項(xiàng)垂蕩運(yùn)動(dòng)滿足頂張式立管的要求，系泊系統(tǒng)滿足規(guī)范要求。干樹半潛式平臺的提出及分析計(jì)算，對我國開發(fā)深水油氣資源提供了一定的參考意義。

[1] 董艷秋. 深海采油平臺波浪載荷及響應(yīng)[M]．天津：天津大學(xué)出版社，2005．

[2] Zhang D G, Chen Y J and Zhang T Y. Floating production platforms and their applications in the development of oil and gas fields in the south china sea[J]. Journal of Marine Science and Application, 2014, 13: 67-75.

[3] 羅勇, 高巍, 申輝, 等. 新型深水干樹半潛平臺選型及總體性能分析研究[J]. 中國造船(增刊), 2013(A02)：30-40.

[4] Williams N, Leverette S, Bian S, et al. FourstarTM Dry-Tree Semisubmersible Development[C]. OMAE Conference, 2010.

[5] 王桂林, 段夢蘭, 王瑩瑩，等. 干式井口半潛式生產(chǎn)平臺技術(shù)與應(yīng)用[J]. 中國海洋平臺, 2010, 25(5):11-14.

[6] 劉海霞. 深海半潛式鉆井平臺的發(fā)展[J]. 船舶, 2007, 6(3) : 6-10．

[7] 姜哲, 謝彬, 謝文會(huì). 新型深水半潛式生產(chǎn)平臺發(fā)展綜述[J]. 海洋工程, 2011, 29(3): 132-138.

[8] N.D.P. Barltrop. Floating Structures: a guide for design and analysis[M]. Houston: Oilfield Publications, 1998.

[9] Phillip A, Hugh B, Jun C, et al. A Deep Draft Semisubmersible with a Retractible Heave Plate[C]. Proceedings Offshore Technology Conference, Houston, 2002.

[10] John M, Arcandra T, Chan K Y. Hydrodynamics of Dry Tree Semisubmersibles[C]. ISOPE，Lisbon Portugal, 2007.

[11] API-RP-2SK. Recommended practice for design and analysis of station keeping systems for floating structures[S].2005.

[12] DNV.Wave Analaysis by Diffraction and Morison Theory[M]. 2008.

[13] 吳瀾, 吳寶山, 匡曉峰, 等. 基于AQWA的半潛式平臺水動(dòng)力特性研究[J]. 中國海洋平臺, 2014, 29(5): 29-33.

[14] DNV.SESAM User Manual HydroD Version 4.0[M]. 2008.

[15] 桂龍, 唐友剛, 秦堯, 等. 淺水FPSO限位作業(yè)系泊分析[J]. 海洋工程, 2014, 32(3): 28-35.

[16] API-RP-2SM.Recommended Practice for Design, Manufacture, Installation and Maintenance of Synthetic Fiber Ropes for Offshore Mooring[S]. 2001.

Dry-tree Semi-submersible Platform′s Hydrodynamic Performance Research Based on the Design of Cancellation Period

LIANG Ning1, HUANG Wei-ping1, ZHOU Yang2, CHANG Shuang1, SUN Yuan-kun1

(1. Ocean University of China， Shandong Qingdao 266100, China；2.The Second Institute of Oceanography, SOA, Hangzhou 310012, China.)

Based on the cancellation period, design the dry-tree semi-submersible platform and do research on its hydrodynamic performance. Deep draft dry-tree semi-submersible platform don't use accessory structure such as heave plate, it combining the advantages of dry-tree development system with semi-submersible platform. By optimizing the structure's principal dimensions to reduce the heave response so it can adopt the dry-tree system. Based on the three-dimensional potential flow theory, we used the theoretical analysis with numerical calculation to research the dynamic characteristics of the dry-tree semi-submersible by SESAM. The platform's added mass, potential damping, response amplitude operator and other parameters are calculated in frequency domain analysis by HydroD. Time domain coupling analysis was carried out on the platform by DeepC to get platform with mooring system's time domain response.

dry-tree semi-submersible platform; hydrodynamic performance; structural design

2015-12-14

國家自然科學(xué)基金(51239008,51179179)。

梁 寧(1990-)，男，碩士研究生。

1001-4500(2016)05-0084-09

P751

A