袁洪濤，曾 驥，莫 建，康 莊，王 超，王鈺涵

(1.上海外高橋造船有限公司，上海 200137)

(2.哈爾濱工程大學(xué)船舶工程學(xué)院，黑龍江哈爾濱150001)

隨著全球變暖，極端氣候環(huán)境現(xiàn)象越來越頻繁，某些極端的風(fēng)浪流環(huán)境已經(jīng)超出了海上設(shè)施的極限能力.在科學(xué)理解極端工況，建立極端工況分析模型以及在確保工作人員安全的前提下，正確評估船體的實際抗風(fēng)浪流能力，減少船體結(jié)構(gòu)的損壞，防止海洋污染是目前急需解決的問題.

FPSO在生命周期內(nèi)都固定在目標海域，即使在百年一遇的極端海洋環(huán)境，也不會解脫進行避航作業(yè)，所以在各種海洋環(huán)境條件對FPSO的運動性能和系泊系統(tǒng)總體性能的預(yù)報也勢在必行.目前，有關(guān)極端海況條件下海上浮體運動響應(yīng)的文獻并不是很多，而且大多集中在百年一遇的工況.文獻［1］中分析了百年一遇極端工況下FPSO系泊系統(tǒng)靜態(tài)和動態(tài)響應(yīng).文獻［2］中運用先靜力在位分析，再利用海洋環(huán)境回歸參數(shù)進行倒塌分析得到導(dǎo)管架平臺在極限環(huán)境條件下的極限承載力.文獻［3］中在時域分析時考慮系泊系統(tǒng)與FPSO的耦合作用，預(yù)報南海500年一遇季風(fēng)和500年一遇臺風(fēng)極端海洋環(huán)境條件下的運動響應(yīng)和系泊纜受力，并在試驗水池進行FPSO在500年一遇臺風(fēng)和500年一遇季風(fēng)條件下的模型試驗，而且試驗結(jié)果與數(shù)值結(jié)果吻合良好.文獻［4］中對工作在墨西哥灣海域轉(zhuǎn)塔系泊方式FPSO在百年一遇極端海況極端條件下的運動響應(yīng)進行研究，分析了FPSO在百年一遇極端海洋環(huán)境中的運動性能以及系泊線受力.文獻［5］中對塔系泊方式FPSO在線性和非線性環(huán)境條件的運動響應(yīng)分別進行了數(shù)值和試驗研究.文獻［6］中研究了百年一遇的颶風(fēng)環(huán)境中內(nèi)轉(zhuǎn)塔系泊的FPSO的運動響應(yīng)，并且進行時域分析計算，得到了FPSO在颶風(fēng)環(huán)境運動時歷曲線.文獻［7］中對單點系泊系統(tǒng)FPSO在瞬時風(fēng)速可達50 m/s的暴風(fēng)海況條件下的運動響應(yīng)分析計算.文獻［8］中采用時域耦合分析方法對單點系泊系統(tǒng)FPSO在914.4 m水池、1 828.8 m水池以及3 048 m水深的運動響應(yīng)做分析研究，分析過程中考慮立管和系泊纜的耦合作用，并對不同阻尼系數(shù)對FPSO運動響應(yīng)影響進行比較分析，研究結(jié)果表明建立全耦合FPSO數(shù)值模型是有必要的;文獻［9］中定量分析了系泊系統(tǒng)不同參數(shù)對FPSO系統(tǒng)抵抗環(huán)境載荷能力的影響;文獻［10］中采用理論結(jié)合試驗的方法研究了試驗過程中，不同的系泊纜參數(shù)對實驗準確性的潛在影響;文獻［11］中探究了布鏈形式、預(yù)張力對系泊系統(tǒng)整體的影響，提出了適合于環(huán)境載荷的大型集裝箱船碼頭系泊形式.文中主要分析了500年一遇波浪主控環(huán)境條件下Nominal波峰Tp，95%Tp和 97% 有義波高 Hs，90%Tp和 92%Hs，85%Tp和87%Hs，500年一遇海況風(fēng)主控環(huán)境條件，500年一遇海況流主控環(huán)境條件，為FPSO的基本設(shè)計和詳細設(shè)計提供依據(jù)，用數(shù)值模擬的方法預(yù)報極端海洋環(huán)境下的運動響應(yīng)和系泊系統(tǒng)總體性能，獲得的結(jié)果對于今后海洋工程領(lǐng)域的研究具有參考意義.

1 錨泊系統(tǒng)的耦合分析

1.1 FPSO輻射繞射分析

在進行時域耦合分析之前，必須運用三維勢流理論對平臺進行輻射繞射分析.通過運動方程式(1)，求解在不同的波浪入射角和頻率的規(guī)則波下平臺的運動響應(yīng)，即RAO.

式中:Ms為平臺質(zhì)量矩陣;Ma為平臺附加質(zhì)量矩陣;C為平臺阻尼矩陣;Ks為平臺總剛度矩陣;F為平臺受到的波浪力矩陣;ω為入射波頻率.

1.2 系泊線受力分析

AQWA軟件采用基于懸鏈線法有限元法分析系泊線靜態(tài)受力，具體的懸鏈線方程式見(2～5)，動態(tài)受力采用集中質(zhì)量法.應(yīng)用離散思想對系泊纜進行分段，段與段之間通過有質(zhì)量的節(jié)點進行連接，段是沒有質(zhì)量的且被看作是剛體或者彈性體，同時假設(shè)系泊纜的重力、浮力、慣性力以及流體拖曳力等都集中作用于節(jié)點上，每段系泊纜的質(zhì)量、附連水質(zhì)量和流體拖曳力平均分配到每段兩端的節(jié)點上.

式中:H為系泊線頂部水平張力;l為單位長度;V為系泊線頂部垂向張力;T為總張力;L為系泊線未變形的原長;AE為單位長度系泊線剛度;X為系泊半徑;w為單位長度濕重;Z為系泊線的垂直高度.

1.3 時域耦合分析

浮體的水動力分析采用間接時域法、水動力計算軟件AQWA基于三維勢流理論，分別對浮體的幅值響應(yīng)函數(shù)、一階波浪力、波浪漂移力、附加質(zhì)量和阻尼系數(shù)等水動力參數(shù)進行數(shù)值求解.得到水動力參數(shù)后應(yīng)用傅里葉變換在全耦合分析軟件OrcaFlex中進行時域分析.浮體運動方程可表示為:

式中:Ms為浮體質(zhì)量;Ma為浮體附加質(zhì)量;K為靜水剛度;x(t)為位移矢量;h(t)為加速度卷積積分函數(shù);F(t)為總的外部力，包含波浪激勵力、系泊力、漂移力、阻尼力和定常力等.

最后利用全耦合的方法分析浮體及系泊纜的耦合運動，其運動方程為:

式中:M為系統(tǒng)慣性力矢量;C為系統(tǒng)阻尼力矢量;K為系統(tǒng)基于單元剛度的每一瞬時應(yīng)力;F包括重力、浮力、系泊力、海洋環(huán)境力和一些特定的外力等.

2 FPSO系泊系統(tǒng)以及環(huán)境參數(shù)

2.1 FPSO系統(tǒng)參數(shù)

表1為15萬噸FPSO具體參數(shù).圖1，2分別為FPSO有限元模型和系泊線布置方式，系泊線采用三段式結(jié)構(gòu)，系泊線具體參數(shù)見表2.

表1 FPSO主尺度參數(shù)及屬性Table 1 Principle dimension and properties of FPSO

圖1 FPSO有限元模型Fig.1 Finite element model of FPSO

圖2 系泊布置分布Fig.2 Mooring line arrangement

表2 15萬噸FPSO系泊系統(tǒng)參數(shù)Table 2 Mooring system parameters of 150000 DWT FPSO

2.2 環(huán)境參數(shù)

驗證15萬噸單點系泊的FPSO在500年一遇極端海況條件下FPSO的運動響應(yīng)以及單根系泊線張力響應(yīng).具體的環(huán)境工況如表3～6所示.表中:Hs為有義波高，Tp為波峰，θwα為波浪角度，Vwi為風(fēng)速，θwi為風(fēng)向，Vcu為流速，θcu為流向.

表3 500年一遇海況波浪主控環(huán)境條件Table 3 Environmental conditions under 500-year return period induced by wave

表4 500年一遇海況風(fēng)主控環(huán)境條件Table 4 Environmental conditions under 500-yearreturn period induced by wind

表5 500年一遇海況流主控環(huán)境條件Table 5 Environmental conditions under 500-yearreturn period induced by current

表6 500年一遇海況波浪主控環(huán)境條件—95%Tp，97%HsTable 6 Environmental conditions under 500-year return period induced by wave(95%Tp＆97%Hs)

3 計算結(jié)果及分析

時域計算采用數(shù)值軟件OrcaFlex，表7為500年一遇極端海況條件下FPSO的運動響應(yīng)以及單根系泊線張力F響應(yīng)，從圖3～5中可以看出，500年一遇的海況環(huán)境下FPSO壓載狀態(tài)下運動響應(yīng)幅值明顯大于滿載狀態(tài)下的運動響應(yīng)幅值，部分工況壓載狀態(tài)比滿載狀態(tài)下縱蕩位移D超出100 m;環(huán)境條件為500年一遇海況風(fēng)主控環(huán)境條件時，F(xiàn)PSO滿載載況時比500年一遇海況流主控環(huán)境條件時要縱蕩橫蕩值大，但是FPSO壓載載況時要小.這是由于壓載載況時，F(xiàn)PSO受風(fēng)面積大，風(fēng)對該FPSO影響比流的影響要大.由表8可以看出，不同Hs和Tp對FPSO動態(tài)響應(yīng)影響也不一樣，無論滿載和壓載載況，低的Hs和Tp都可以降低FPSO的縱蕩和橫蕩值.當風(fēng)浪流成夾角時，系泊線受力最大.

表7 極端海況條件下FPSO的運動響應(yīng)以及單根系泊線張力Table 7 FPSO Motion response and mooring line tension under extreme sea conditions

圖3 波浪主控環(huán)境結(jié)果對比Fig.3 Comparison results under wave dominate condition

圖4 風(fēng)主控環(huán)境結(jié)果對比Fig.4 Comparison results under wind dominate condition

圖5 流主控環(huán)境結(jié)果對比Fig.5 Comparison results under current dominate condition

表8 500年一遇海況波浪主控環(huán)境條件不同Hs和Tp下FPSO的運動響應(yīng)以及單根系泊線張力Table 8 FPSO motion response and mooring line tension under wave dominate condition

4 結(jié)論

1)應(yīng)用大型水動力軟件AQWA和OCRAFLEX對15萬噸深水單點系泊系統(tǒng)FPSO在500年一遇海況風(fēng)浪流各主控環(huán)境下進行了頻域和時域分析，經(jīng)過頻域分析得到了RAO和QTF等一系列水動力參數(shù).FPSO在極端的海洋環(huán)境中產(chǎn)生的大幅慢漂運動對FPSO的安全有著至關(guān)重要的影響，通過動態(tài)耦合分析，獲得了FPSO在極端海況下的運動響應(yīng)和系泊系統(tǒng)的總體性能.

2)FPSO壓載狀態(tài)下運動響應(yīng)幅值大于滿載狀態(tài)下的運動響應(yīng)幅值，不同的主控環(huán)境對FPSO的運動響應(yīng)影響也不一樣，波浪主控環(huán)境對FPSO運動的影響最大.風(fēng)和流的影響與FPSO船型具體的風(fēng)流系數(shù)關(guān)系密切.Tp和Hs對FPSO的運動響應(yīng)幅值也有較大的影響且風(fēng)浪流成夾角作用時對FPSO運動響應(yīng)的影響更加明顯.

3)風(fēng)與流對錨泊線的影響相比于普通環(huán)境有較大程度的增強.波浪入射角的方向?qū)ο挡淳€的張力影響較大，滿載狀態(tài)下的系泊線最大張力又比壓載狀態(tài)下的系泊線張力大.

References)

［1］ 王宏偉，羅勇，孫麗萍.內(nèi)轉(zhuǎn)塔式FPSO截斷技術(shù)及其總體響應(yīng)［J］.天津大學(xué)學(xué)報，2010，43(11):989-990.Wang Hongwei，Luo Yong，Sun Liping. Truncation method of internal turret FPSO and its global response［J］.Journal of Tianjing University，2010，43(11):989-990.(in Chinese)

［2］ 熊卓，夏凡，黃誠.導(dǎo)管架平臺極限環(huán)境下承載力分析［J］.中國水運，2013(12):186-187.

［3］ 鄭翔.FPSO系統(tǒng)水動力混合模型試驗研究［D］.黑龍江哈爾濱:哈爾濱工程大學(xué)，2012.

［4］ Baar J J M，Heyl C N，Rodenbusch G.Extreme respon-ses of turret moored tankers［C］∥Proc Offshore Technology Conference.Houston，Texas:［s.n.］，2000.

［5］ Ward E G，Irani M B，Johnson R P.The behavior of a tanker-based FPSO in Hurricane waves，winds，and currents［C］∥11th International Offshore and Polar Eng.Conference.［S.l.］:ISOPE，2001:650 -653.

［6］ Arcandra T，Nurtjahyo P，Kim M H.Hull/mooring/riser coupled analysis of a turret moored FPSO 6000 ft:comparison between polyester and buoys-steel mooring lines［C］∥Proc.11th Offshore Symposium.The Texas Section of the Society of Naval Architects and Marine Engineers.USA:SNAME，2002:1 -8.

［7］ Duggal A，Heyl C，Amir H，et al.Response of FPSO system to squalls［C］∥Proceedings of the ASME 2011 30th International Conference on Offshore Mechanics and Arctic Engineering OMAE 2011.Rotterdam，The Netherlands，USA:ASME，2011:19 -24.

［8］ Wichers J E W，Devlin P V.Effect of coupling of mooring lines and risers on the design values for a turret moored FPSO in deep water of the Gulf of Mexico［C］∥Proceedings of the Eleventh International Offshore and Polar Engineering Conference Stavange.Norway:［s.n.］，2001:7 -22.

［9］ 王志東，劉曉健，陳劍文.FPSO系泊系統(tǒng)參數(shù)的動力效應(yīng)影響分析［J］.江蘇科技大學(xué)學(xué)報:自然科學(xué)版，2013，27(5):409 -418.Wang Zhidong，Liu Xiaojian，Chen Jianwen.Analysis of the influence of mooring system parameters on dynamic response of FPSO system［J］.Journal of Jiangsu University of Science and Technology:Natural Science Edition，2013，27(5):409 -418.(in Chinese)

［10］ 喬薛峰，肖龍飛，楊立軍.截斷系泊系統(tǒng)不同恢復(fù)力矩特性對于深水半潛式平臺搖擺運動的影響［J］.上海交通大學(xué)學(xué)報，2014(4):450-456.Qiao Xuefeng，Xiao Longfei，Yang Lijun.Effect of moment characteristics of truncated mooring system on roll and pitch motions of semi-submersible platform［J］.Journal of Shanghai Jiaotong University，2014(4):450-456.(in Chinese)

［11］ 張鼎，黃維，俞赟.大型集裝箱船碼頭系泊及抗臺風(fēng)計算分析［J］.船舶與海洋工程，2014(2):16-19.