韓 熠，楊 輝，劉 波，楊 亮，方 堃

（海洋石油工程股份有限公司，天津 300451）

0 引言

作為海上油田的生產(chǎn)儲(chǔ)油設(shè)備，大型浮式生產(chǎn)儲(chǔ)油卸油裝置（Floating Production Storage and Offloading，F(xiàn)PSO）長(zhǎng)期連續(xù)作業(yè)于固定海域，雖然多數(shù)FPSO采用船型結(jié)構(gòu)，但由于不解脫不定期進(jìn)塢的要求使得FPSO的船體強(qiáng)度和疲勞強(qiáng)度都比一般船舶的標(biāo)準(zhǔn)更高。因此，在前期設(shè)計(jì)階段中，需對(duì)FPSO整個(gè)生命周期內(nèi)可能受到的波浪誘導(dǎo)彎矩、剪力以及運(yùn)動(dòng)參數(shù)做出可靠預(yù)報(bào)。

當(dāng)FPSO系泊在一個(gè)淺水海域，浮體的吃水與水深相對(duì)接近，此時(shí)既要保持大的裝載能力，又要具有良好的運(yùn)動(dòng)性能以保證海上作業(yè)安全，因此，問(wèn)題更加復(fù)雜化。近年來(lái)，人們對(duì)淺水大型浮式結(jié)構(gòu)的研究愈發(fā)深入，賀五洲等[1]、楊玥等[2]曾提出，對(duì)于中等以上水深，F(xiàn)PSO的響應(yīng)主要取決于風(fēng)、浪、流及其系泊系統(tǒng)，與水深關(guān)系并不大，但是對(duì)于淺水而言，水深效應(yīng)逐漸明顯。隨著水池試驗(yàn)的發(fā)展，謝永和等[3]通過(guò)試驗(yàn)數(shù)據(jù)證明當(dāng)水深變淺時(shí)，垂蕩和縱搖的響應(yīng)峰值變小，彎矩和剪力則會(huì)增大，而且水深越淺波浪誘導(dǎo)垂向彎矩與剪力的增幅越明顯。由此可見(jiàn)，淺水效應(yīng)在實(shí)際工程應(yīng)用中會(huì)對(duì)結(jié)構(gòu)優(yōu)化產(chǎn)生影響。而對(duì)于淺水油田前期開(kāi)發(fā)方案的選取，在快速篩選設(shè)計(jì)工況及選定設(shè)計(jì)載荷的過(guò)程中，有以下問(wèn)題值得關(guān)注：1）淺水條件下不同裝載情況載荷預(yù)報(bào)的差異；2）使載荷響應(yīng)達(dá)到最大的環(huán)境條件是什么；3）由于淺水效應(yīng)的存在，鋼制海船入籍規(guī)范所推薦的載荷規(guī)范值是否依然保守；4）淺水效應(yīng)對(duì)結(jié)構(gòu)設(shè)計(jì)而言會(huì)帶來(lái)哪些關(guān)注點(diǎn)。

為探討以上問(wèn)題，筆者以渤海海域大型 FPSO為例，計(jì)算其在不同壓載情況下的RAO（Response Amplitude Operator，反應(yīng)量），并基于渤海環(huán)境條件統(tǒng)計(jì)響應(yīng)極值，并與挪威船級(jí)社規(guī)范（DNV GL）推薦值進(jìn)行比較以探究相關(guān)研究結(jié)論。

1 淺水規(guī)則波船體載荷計(jì)算理論

規(guī)則波中的三維源匯理論假設(shè)船體運(yùn)動(dòng)及波浪均是微幅、線(xiàn)性的，流體視為無(wú)粘性的理想流體，則船體的運(yùn)動(dòng)方程為

式中：ω為圓頻率；β為浪向角；M為廣義質(zhì)量矩陣；A(ω)為附加質(zhì)量陣；B(ω)p為勢(shì)流阻尼陣；Bv是線(xiàn)性粘性阻尼陣；C為靜水回復(fù)力矩陣；Ce為外部回復(fù)力矩陣（錨鏈力等）；F(ω,β)為波浪激勵(lì)力矩陣[4]。

計(jì)算結(jié)構(gòu)響應(yīng)首先要計(jì)算波浪力F(ω,β)，常用勢(shì)流理論的方法，速度勢(shì)可線(xiàn)性分解為3部分，入射勢(shì)、繞射勢(shì)及物體自身運(yùn)動(dòng)產(chǎn)生的輻射勢(shì)。速度勢(shì)在流場(chǎng)中滿(mǎn)足拉普拉斯方程及4類(lèi)邊界條件（定解條件）：自由表面條件、物體濕表面條件、海底條件及輻射條件（無(wú)窮遠(yuǎn)處邊界條件），再用伯努利方程計(jì)算物體濕表面壓強(qiáng)。求解速度勢(shì)常用格林函數(shù)法[2]。船舶在淺水中受力狀況有所改變，數(shù)值解析主要體現(xiàn)在有限水深的格林函數(shù)積分表達(dá)式為奇異二重積分，這樣當(dāng)水深與船舶吃水的比值小于一定值時(shí)，船舶的水動(dòng)力系數(shù)及所受波浪擾動(dòng)力等與深水情況有所不同。一般認(rèn)為，當(dāng)水深小于4倍吃水時(shí)就會(huì)出現(xiàn)淺水效應(yīng)，而小于2倍吃水時(shí)將發(fā)生很大差異。主要表現(xiàn)有：平面漂移運(yùn)動(dòng)加大、垂直及橫向搖擺運(yùn)動(dòng)減緩等[5]。

2 波浪載荷和運(yùn)動(dòng)預(yù)報(bào)的基本理論

用無(wú)數(shù)個(gè)規(guī)則波相互疊加生成不規(guī)則波，比較符合真實(shí)海面上的波浪，此時(shí)應(yīng)采用概率和隨機(jī)理論的方法計(jì)算波浪誘導(dǎo)船體運(yùn)動(dòng)和波浪載荷。本文將采用短期統(tǒng)計(jì)的方法預(yù)報(bào)100年一遇環(huán)境條件下的極值響應(yīng)。

首先應(yīng)求出船體在規(guī)則波中的傳遞函數(shù)，然后通過(guò)實(shí)際海域的波浪譜計(jì)算船體運(yùn)動(dòng)和載荷的響應(yīng)譜[6]。預(yù)報(bào)方法分為短期預(yù)報(bào)和長(zhǎng)期預(yù)報(bào)。短期預(yù)報(bào)是在假設(shè)海況條件不變的數(shù)小時(shí)內(nèi)船體的響應(yīng)值。大量事實(shí)數(shù)據(jù)表明，船體運(yùn)動(dòng)和載荷的短期響應(yīng)服從瑞利分布。

短期海況通常用波浪譜來(lái)表示其特性，海浪譜密度函數(shù)S(ω)是平穩(wěn)隨機(jī)過(guò)程的頻率描述，反應(yīng)了不規(guī)則波內(nèi)各單元諧波的能量分布情況。短期預(yù)報(bào)的波浪幅值及波浪誘導(dǎo)船體運(yùn)動(dòng)幅值、載荷幅值、應(yīng)力幅值符合瑞利分布，其概率密度函數(shù)為

概率累計(jì)函數(shù)為

三一有義值、百一平均值等響應(yīng)值均是在一定超越概率下得到的統(tǒng)計(jì)值。

3 波浪載荷和運(yùn)動(dòng)分析預(yù)報(bào)

本文研究對(duì)象是渤海油田15萬(wàn)噸級(jí)大型船型生產(chǎn)儲(chǔ)油系統(tǒng)，該大型FPSO總長(zhǎng)L=276 m，型寬B=51 m，型深D=23.6 m，滿(mǎn)載排水量Δ=199 959.9 t，滿(mǎn)載吃水T=15.66 m，空載排水量Δ=118 705.4 t，空載吃水T=9.45 m，方形系數(shù)Cb=0.93。

在DNV GL的SESAM軟件中進(jìn)行模擬時(shí)，使用前處理模塊GENIE建立FPSO水動(dòng)力濕表面模型和質(zhì)量模型，見(jiàn)圖1和圖2；用WADAM模塊計(jì)算船體運(yùn)動(dòng)和波浪載荷的傳遞函數(shù)[6-8]。

圖1 濕表面結(jié)構(gòu)

圖2 質(zhì)量模型示例

3.1 水動(dòng)力模型

設(shè)定該FPSO將服役于渤海海域，工作水深只有22 m，水深吃水比僅為1.41，水動(dòng)力模型全部由三角形和四邊形單元組成。全船共有7 414個(gè)單元和7 728個(gè)節(jié)點(diǎn)。

用質(zhì)量棒來(lái)模擬整船質(zhì)量分布，可以得到截面載荷的分布。FPSO的總重量可分為空船重量和裝載重量，將空船重量保持一致，裝載重量根據(jù)裝載手冊(cè)進(jìn)行調(diào)整，使得慣性半徑與重量資料數(shù)據(jù)一致，可分別得到滿(mǎn)載和壓載工況下的質(zhì)量模型。

3.2 坐標(biāo)系定義

為了定義和計(jì)算淺水超大型FPSO的運(yùn)動(dòng)和波浪誘導(dǎo)載荷響應(yīng)，本文建立了淺水大型FPSO的整體坐標(biāo)系，見(jiàn)圖3。整體坐標(biāo)系原點(diǎn)位于FPSO基線(xiàn)的中心處，x軸正方形從船尾指向船首，y軸正向指向左舷，z軸向上為正。整體坐標(biāo)系用于描述入射波的方向和入射波的速度勢(shì)。

圖3 參考坐標(biāo)

3.3 裝載工況

為得到船體運(yùn)動(dòng)和載荷的預(yù)報(bào)極值，分別考慮船體滿(mǎn)載和空載兩種裝載狀態(tài)，見(jiàn)表1。

表1 裝載工況

3.4 傳遞函數(shù)及響應(yīng)統(tǒng)計(jì)

為研究淺水大型FPSO在特定水深下波浪誘導(dǎo)載荷響應(yīng)隨浪向變化的規(guī)律，選取0°～180°共13個(gè)浪向角，以15°為間隔。波浪譜參數(shù)見(jiàn)表2。

沿船長(zhǎng)方向取17個(gè)參考截面，采用JONSWAP譜做短期預(yù)報(bào)[9-10]。

重點(diǎn)考慮距船尾船首1/4L處剪力和船中位置彎矩，得到的RAO如圖4～圖6所示。重點(diǎn)位置短期統(tǒng)計(jì)極值見(jiàn)表3，運(yùn)動(dòng)加速度極值見(jiàn)表4。

表2 波浪譜參數(shù)

3.5 規(guī)范計(jì)算法與直接計(jì)算法比較

挪威船級(jí)社對(duì)于長(zhǎng)度大于100 m的船舶提供了載荷的規(guī)范值[11]。

靜水剪力為正時(shí)，

圖4 典型垂向剪力和Y向彎矩RAO曲線(xiàn)

圖5 壓載工況下垂蕩、縱搖、橫搖RAO曲線(xiàn)

圖6 滿(mǎn)載工況下垂蕩、縱搖、橫搖RAO曲線(xiàn)

表3 重點(diǎn)位置短期統(tǒng)計(jì)極值

表4 運(yùn)動(dòng)加速度極值

靜水剪力為負(fù)時(shí)，

將規(guī)范推薦值和直接計(jì)算得到的數(shù)值相比較，所得結(jié)果如圖7所示。

圖7 2種裝載下波浪載荷與規(guī)范值的比較

1）滿(mǎn)載工況下的波浪載荷比壓載工況下的大很多。這是因?yàn)闈M(mǎn)載時(shí)吃水增加會(huì)承受更大的波浪力，然而在壓載工況下船體運(yùn)動(dòng)加速度相對(duì)更大。

2）滿(mǎn)載工況下傳遞函數(shù)的峰值區(qū)要比壓載工況下的尖銳，滿(mǎn)載工況下傳遞函數(shù)的最大值對(duì)應(yīng)的波浪周期在12 s～14 s之間，壓載工況下傳遞函數(shù)最大值所對(duì)應(yīng)的波浪周期為14 s～16 s。由于FPSO滿(mǎn)載狀態(tài)下的總慣性矩更大，導(dǎo)致滿(mǎn)載狀態(tài)的固有周期大于壓載工況下的固有周期。

3）對(duì)比2種工況下FPSO的傳遞函數(shù)，可以發(fā)現(xiàn)較大的垂蕩、橫搖和縱搖運(yùn)動(dòng)幅值主要集中在12 s～16 s之間的波浪周期范圍內(nèi)，波長(zhǎng)范圍大約在176 m～235 m，而對(duì)于波浪載荷而言，主要幅值集中在17 s～22 s之間，波長(zhǎng)范圍為249 m～323 m。

4）從計(jì)算結(jié)果來(lái)看，大部分位置產(chǎn)生的彎矩和剪力均小于規(guī)范推薦值，說(shuō)明對(duì)于水深吃水比在1.4左右的船舶而言，規(guī)范值還是相對(duì)保守的。但是在船體尾部剪力有一個(gè)突變，甚至大于規(guī)范推薦值，這是因?yàn)橐?guī)范公式并未包含水深項(xiàng)，而FPSO的裝載量要求及水平尺度都很大，這就導(dǎo)致其結(jié)構(gòu)剛度相對(duì)剪切剛度降低，淺水效應(yīng)放大了波浪載荷值，這在實(shí)際設(shè)計(jì)時(shí)應(yīng)受到重視。

4 結(jié)論

1）前期的設(shè)計(jì)工況單一考慮滿(mǎn)載或壓載都是不合理的。對(duì)于對(duì)載荷敏感的船體結(jié)構(gòu)而言，應(yīng)選取滿(mǎn)載工況為設(shè)計(jì)工況；對(duì)于對(duì)運(yùn)動(dòng)敏感的上部組塊設(shè)計(jì)而言，壓載工況應(yīng)被選為設(shè)計(jì)工況。

2）當(dāng)淺水大型 FPSO遭遇的波長(zhǎng)與船長(zhǎng)接近時(shí)，載荷幅值最大。

3）進(jìn)行結(jié)構(gòu)設(shè)計(jì)時(shí)，對(duì)于船體中部貨油艙的艙段分析，可以以載荷規(guī)范值作為設(shè)計(jì)值；對(duì)于首尾部線(xiàn)型過(guò)渡明顯的區(qū)域，則需根據(jù)載荷預(yù)報(bào)評(píng)估來(lái)確定設(shè)計(jì)載荷值。

4）渤海地區(qū)良好的海況使得載荷統(tǒng)計(jì)值遠(yuǎn)小于規(guī)范值。這說(shuō)明對(duì)于主要承擔(dān)總縱強(qiáng)度的貨油艙區(qū)域而言，其結(jié)構(gòu)選型可以適當(dāng)趨于臨界以保證其經(jīng)濟(jì)性。而淺水效應(yīng)下尾部剪力加大的情況，可以通過(guò)適當(dāng)增加船尾橫艙壁數(shù)量的方法來(lái)解決。這在前期總體布置規(guī)劃中可以給予適當(dāng)考量。