帶月池深海鉆井船搖蕩運動性能分析

張曉宇胡開業(yè)周雯雯（.中國石油渤海裝備研究院 海工裝備分院 盤錦400； .哈爾濱工程大學 船舶工程學院 哈爾濱5000）

我院主持召開首屆國際噴水推進技術研討會

?

帶月池深海鉆井船搖蕩運動性能分析

張曉宇1胡開業(yè)2周雯雯2
（1.中國石油渤海裝備研究院 海工裝備分院 盤錦124010； 2.哈爾濱工程大學 船舶工程學院 哈爾濱150001）

［摘 要］海上鉆井作業(yè)過程要求鉆井船具備良好的穩(wěn)定性能，而月池內(nèi)流體柱的不規(guī)則振蕩以及鉆井船細長的船體外形都將對船體運動性能造成不利影響，從而限制其在惡劣海域的作業(yè)適應性。文章基于三維勢流理論對帶月池鉆井船的搖蕩運動進行預報，并通過船模試驗對預報結果進行有效性驗證。在此基礎上，提出鉆井船搖蕩運動衡準值，對鉆井船海上作業(yè)適應性進行研究，分析月池尺度和形狀對鉆井船運動性能的影響。

［關鍵詞］深水鉆井船；運動預報；月池

胡開業(yè)（1980-）,男，博士，講師，研究方向：水動力性能分析。

周雯雯（1989-），女，碩士，研究方向：船舶工程。

引 言

與其他鉆井平臺相比，鉆井船的機動性好、自持力強、具有自航遷移能力，更能適應深水及超深水作業(yè)。截止2013年10月，全球現(xiàn)役鉆井船總量為91艘，另外還有88艘正在建造中。在91艘鉆井船中，83艘鉆井船具備在4 000 ft（約合1 219 m）以上深水作業(yè)的能力，比例為91.2%，68艘鉆井船適應于10 000 ft（約3 048 m）以上的超深水海域作業(yè)，比例達74.7%［1］。與國外蓬勃的發(fā)展現(xiàn)狀相比，國內(nèi)鉆井船技術僅處于起步階段。由于鉆井船高附加值的特點，其相關技術壟斷嚴重，研究成果為商業(yè)秘密，增加了國內(nèi)鉆井船技術的發(fā)展難度［2］。本文從運動性能的角度對鉆井船技術進行研究，為國內(nèi)自主船型的開發(fā)提供技術支撐。

海上鉆井作業(yè)對平臺運動性能要求苛刻，過大的運動幅度不僅影響鉆井精度與作業(yè)效率，而且將縮短鉆桿、隔水管等鉆井設備的使用壽命，甚至導致皸裂、脫扣等安全事故，造成巨大的經(jīng)濟損失。然而，鉆井船具有細長的船體外型，受波浪的常態(tài)影響，對運動性能造成不良影響。另外，鉆井船中部布置有貫穿船底的月池結構，月池內(nèi)流體振蕩規(guī)律復雜、產(chǎn)生漩渦等非線性現(xiàn)象，影響流體域的壓力分布及鉆井船搖蕩運動幅度，限制了鉆井船在惡劣海域的作業(yè)適應性［3］。因此，鉆井船在波浪中的搖蕩運動預報與分析研究，對鉆井船的船型論證具有實際意義。

本文應用基于三維零航速勢流理論的AQWA海洋工程水動力分析工具，進行鉆井船搖蕩運動預報，并完成鉆井船搖蕩運動模型試驗。在此基礎上，提出鉆井船搖蕩運動衡準，結合短期預報理論，對相應海況下的作業(yè)適應性進行分析。通過計算對比無月池及不同尺度月池下的鉆井船搖蕩運動響應，分析月池對鉆井船運動性能的影響。

1 鉆井船搖蕩運動預報分析理論

隨著計算機技術的發(fā)展，三維勢流理論日臻完善，在海洋工程結構物水動力分析領域應用廣泛［4］。

1.1 三維勢流理論基礎

對于無航速浮體，流場擾動速度勢為：（1）式中：、分別為入射勢、繞射勢和輻射勢。

根據(jù)線性波浪理論，波浪入射勢可表示為：

、（2）式中：A為波幅，m；k為波數(shù)，rad/m；h為水深，m；g為重力加速度，m/s2；β為浪向與x軸正方向的夾角，°，規(guī)定β=180°為迎浪。

入射勢、繞射勢根據(jù)Green函數(shù)法求取：（3）式中：S為物體邊界，為分布于濕表面上的源強，為物面S上的動點，為場點。

根據(jù)速度勢與Lagrange積分式（2）-式（4），即可獲得流場壓力分布。（4）

流體壓力沿結構物表面積分，得到結構物所受到的總體流體作用力，包括波浪激勵力、輻射力以及流體靜回復力，其中輻射力通常用附加質(zhì)量系數(shù)和輻射阻尼系數(shù)表征。因此，船體在頻域下的一階搖蕩方程表達式為：（5）式中：J和k為運動模態(tài) （j=1～6，k=1～6）；xj為船體運動位移；mkj為船體質(zhì)量；μkj為船體附加質(zhì)量；λkj為船體阻尼系數(shù)；ckj為船體回復系數(shù)；Fwk為船體的波浪干擾力。

1.2 鉆井船頻域運動計算模型

將鉆井船與月池內(nèi)的流體視為一個體系中2個獨立又相互影響的組成部分［5］。相對于鉆井船而言，船體運動受到月池內(nèi)流體運動的影響，如以月池內(nèi)流體垂向運動為影響因素，建立鉆井船運動方程為：（6）式中：x8為月池內(nèi)流體垂向運動位移；ak8、bk8分別為月池垂向運動引起的附加質(zhì)量影響系數(shù)、阻尼影響系數(shù)。

鉆井船以及月池內(nèi)流體體系的運動示意如下頁圖1所示。

月池內(nèi)流體的絕對運動可分解為：

圖1　體系運動示意圖

式中：x3為月池內(nèi)流體與船體一致的運動；x′為月池內(nèi)流體與船體的相對運動。

鉆井船運動預報過程中，由于月池運動x8未知且量級較小，因此視x8≈x3，式（6）可表示為：

同理，計入月池六自由度運動影響的鉆井船搖蕩運動方程為：

式中：akj、bkj分別為月池k自由度上的運動對鉆井船j自由度上的附加質(zhì)量影響系數(shù)、附加阻尼影響系數(shù)。

據(jù)A·B·Aalbers研究成果顯示，ω＜0.3 rad/s范圍內(nèi)，附加質(zhì)量影響系數(shù)、阻尼影響系數(shù)均可視為常數(shù)。因此，為簡化預報過程且提高預報精度，AQWA數(shù)值計算過程中，將結合附加質(zhì)量影響系數(shù)、阻尼影響系數(shù)、鉆井船粘性阻尼系數(shù)進行水動力系數(shù)修正［6-8］。

1.3 搖蕩運動短期預報理論

為了獲得鉆井船在隨機波中的搖蕩運動統(tǒng)計值，根據(jù)譜分析理論對鉆井船搖蕩運動幅度進行短期預報。主要過程為：

（1）搖蕩運動幅頻響應函數(shù)（RAO）。

式中：ω0為譜峰頻率；γ為譜峰升高因子；σ為譜峰形狀參數(shù)；α為無因次常數(shù)。

（4）搖蕩運動位移單幅有義值：

1.4 鉆井船搖蕩運動衡準

鉆井船具有細長的船體外形且水線面積較大，在波浪中運動幅值較大，將可能影響海上正常鉆井作業(yè)。因此，在船型設計階段，需要結合搖蕩運動衡準及運動幅值預報結果，對不同海況下的鉆井船作業(yè)適應性進行分析。然而，目前尚無權威機構規(guī)定明確的鉆井船搖蕩運動衡準選取標準，這是鉆井船作業(yè)適應性分析的技術瓶頸。

然而，鉆井船是具有船體外形的海上鉆井平臺，本文在結合常規(guī)船舶耐波性要求及海洋平臺作業(yè)特點基礎上，提出鉆井船作業(yè)過程中（安裝防噴器）的搖動運動衡準如表1所示［9］。

表1　鉆井船搖蕩運動衡準值

1.5 月池特性分析方法

由于月池結構貫穿船底，池內(nèi)流體柱與外部海水相連，在波浪作用下，流體柱將產(chǎn)生“活塞”和“晃蕩”兩種不規(guī)則的振蕩，影響流場壓力分布。目前，關于月池的研究均停留在以模型試驗和理論研究為主的自身特性研究。Faltinsen、Fukuda、Molin等人給出的月池流體振蕩固有頻率經(jīng)驗計算模型，Aalbers、Knott 和 Flower等人通過模型試驗對月池的共振頻率下的附加阻尼或附加質(zhì)量進行研究。English、Park J J等人以減小航行阻力提高遷移速度為目標，對月池尺度進行優(yōu)化設計［10］。

Molin通過推導月池內(nèi)流體做活塞運動和晃蕩運動時的速度勢得到其固有頻率計算公式如式（13）和式（14）所示。

式中：d為月池吃水，m；l為月池長度，m；b 為月池寬度，m；g為重力加速度，m/s2。式中：h為船體吃水，m；S為月池剖面面積，m2；λ=π/l，r=b/l，tanθ0=r-1。

由Molin經(jīng)驗公式估算目標船型月池固有頻率：活塞運動固有頻率=3.948）；晃蕩運動固有頻率ωsloshing=1.147 rad/s

為分析月池的存在及其尺度對鉆井船搖蕩運動的影響，對不同尺度下的鉆井船搖蕩運動響應進行計算分析，月池尺度如表2所示［11］。

表2　月池尺度m

2 實例計算與分析

本文以1艘科研鉆井船型為例，對該鉆井船的運動響應進行預報并開展船模搖蕩運動試驗，在此基礎上，分析指定海況下的搖蕩運動適應性以及月池尺度對鉆井船搖蕩運動的影響。

2.1 數(shù)值計算模型

實例鉆井船主要參數(shù)見表3。

表3　鉆井船主要參數(shù)

數(shù)值計算幾何模型如圖2所示。

圖2　鉆井船幾何模型

2.2 搖蕩運動數(shù)值預報與船模試驗

基于三維零航速勢流理論，應用水動力分析工具AQWA，獲得鉆井船在規(guī)則波中的搖蕩運動幅頻響應。

下頁圖3表示不同浪向下的垂蕩幅頻響應，縱坐標為單位波幅的垂蕩位移。

下頁圖4和圖5分別表示不同浪向下的橫搖、縱搖幅頻響應，縱坐標表示單位波傾角的傾斜弧度。

圖3　垂蕩運動幅頻響應

圖4　橫搖運動幅頻響應

圖5　縱搖運動幅頻響應

為了驗證數(shù)值計算結果的有效性，開展該鉆井船作業(yè)裝載工況（作業(yè)吃水見表3）與典型浪向下的搖蕩運動模型試驗，浪向角β=180°時為正迎浪。圖6為船模試驗過程截圖。

鉆井船搖蕩運動數(shù)值計算結果與模型試驗結果如圖7 -圖9所示。數(shù)值計算值與模型試驗值吻合度較好，驗證了數(shù)值計算結果的有效性及數(shù)值計算方法的可靠性。

圖6　船模試驗截圖

圖7　垂蕩運動響應（β=90°）

圖8　橫搖運動響應（β=165°）

圖9　縱搖運動響應（β=180°）

2.3 鉆井船作業(yè)適應性分析

根據(jù)譜分析理論，對不規(guī)則波中的搖蕩運動進行短期預報，并結合表1中的搖蕩運動衡準值，分析該鉆井船在相應海況下的作業(yè)適應性。

圖10 -圖12分別表示不同海況下的垂蕩、橫搖、縱搖位移以及對應的衡準限界線玫瑰圖。

圖10　垂蕩運動位移玫瑰圖

圖11　橫搖運動位移玫瑰圖

圖12　縱搖運動位移玫瑰圖

圖13 -圖15分別表示不同海況下的垂蕩、橫搖、縱搖速度以及對應的衡準限界線玫瑰圖。

圖13　垂蕩運動速度玫瑰圖

圖14　橫搖運動速度玫瑰圖

圖15　縱搖運動速度玫瑰圖

由圖可知：若該鉆井船采用錨泊定位系統(tǒng)，迎浪方向隨機分布，則只有在有義波高h1/3=1.25 m的低海況下才能滿足衡準值要求，具備良好的作業(yè)適應性。若鉆井船采用動力定位系統(tǒng)，船首可控制在迎浪β=180°至首斜浪β=150°的范圍內(nèi)，則在有義波高上h1/3=5.0 m的高海況下，也能基本滿足衡準要求，保持較好的作業(yè)適應性。由于該鉆井船橫搖固有周期T=21 s，波浪周期范圍大多在15 s以下，因此，橫搖頻率響應與波浪能量譜存在峰值偏移，鉆井船橫搖運動緩慢且搖蕩幅度小，不影響正常鉆井作業(yè)。

2.4 月池對鉆井船搖蕩運動影響分析

建立不同尺度月池的鉆井船幾何模型，基于三維勢流理論得到不同尺度月池下的鉆井船垂蕩、橫搖、縱搖運動響應（β=135°）如圖16 -圖18所示。

圖16　垂蕩運動響應

圖17　橫搖運動響應

圖18　縱搖運動響應

由圖16可知，在月池的“活塞”運動固有頻率附近，鉆井船垂向運動與月池內(nèi)流體的活塞運動發(fā)生共振，鉆井船垂蕩運動曲線出現(xiàn)強烈振蕩，且矩形月池和圓形月池的振蕩現(xiàn)象最為劇烈。

由圖17可知，月池的存在減小了鉆井船的橫搖固有周期，導致橫搖更加劇烈，且明顯增大橫搖峰值，矩形月池的增大幅度最為顯著。

由圖18可知，月池的存在對縱搖響應無影響。

基于以上分析可知，月池結構對鉆井船垂蕩和橫搖兩個方向上的運動造成影響，垂蕩方向上易于發(fā)生共振現(xiàn)象，橫搖方向上表現(xiàn)為增大橫搖幅度?；谧V分析原理，對不同尺度下鉆井船橫搖運動進行短期預報如圖19所示。

圖19　橫搖運動短期預報值（h1/3=9.0 m）

3 結 論

通過以上計算與分析，得出如下結論：

（1）基于AQWA的鉆井船搖蕩運動數(shù)值預報方法較為可靠，鑒于月池存在引起型值和水動力系數(shù)的修正，獲得的鉆井船在波浪中的運動響應與模型試驗結果相吻合，從而驗證了預報方法的可靠性。

（2）保持船首處于迎浪方向有利于鉆井船在惡劣海況下也能保持良好作業(yè)適應性，且鉆井船橫搖方向運動緩慢、幅度又小，對正常作業(yè)并無不良影響。

（3）月池對鉆井船垂蕩和橫搖兩個方向的運動造成不良影響。由于月池內(nèi)流體柱的“活塞”振蕩，鉆井船垂蕩方向上易于發(fā)生共振現(xiàn)象，橫搖方向上表現(xiàn)為橫搖幅度增大，并且矩形月池的增大效果最為明顯，其次為方形月池和圓形月池。

［參考文獻］

［ 1 ］ Anon. Floating into Ultra-Deepwater［EB/OL］.［2014-03-20］. http://www.hellenicshippingnews.com/07be9lal-7e63-4b44-9807-b4e4b26a8573/.

［ 2 ］ 孫采微.深水鉆井船水動力性能研究［D］.上海：上海交通大學，2013：2-7.

［ 3 ］ Molin B. On the piston and sloshing modes in moonpools ［J］. Journal of Fluid Mechanics，2001（1）： 27-50.

［ 4 ］劉應中，繆國平.船舶在波浪中的運動理論［M］.上海：上海交通大學出版社，1987.

［ 5 ］Aalbers A B. The water motion in moonpool［J］.Ocean Engineering， 1984（6）：557-579.

［ 6 ］ Cueva M， Malta E B， Nishimoto K， et al. Estimation of Damping Coefficients of Moonpools for Monocolumn

信息動態(tài)

Type Units［C］//Proc.of OMAE’05， 24th International Conference on Offshore Mechanics and Arctic Engineering. Greece： OMAE， 2005.

［ 7 ］ Maisondieu， Jou?tte D C， Kimmoun C， et al. Numerical and Experimental Evaluation of the Viscous Damping in a Wellbay［C］ //Proc.of OMAE’04， 23rd International Conference on Offshore Mechanics and Arctic Engineering. Canada： OMAE，2004.

［ 8 ］ Fukda K. Behavior of water in vertical well with bottom opening of ship and its effect on ship motions ［J］. Journal of the Society of Naval Architects of Japan， 1977， 141：107-122.

［ 9 ］ Journée J M J，Massie W W. Offshore Hydromechanics ［M］. Delft University of Technology. Delft， 2001.

［10］ Park J J， Kim M S. Development of Design Technologies for Optimum Moonpool Shapes of Drillship ［C］// Proc. of OMAE’07， 26th International Conference on Offshore Mechanics and Arctic Engineering. USA： OMAE， 2005.

［11］ Fung D P K. Added Mass and Damping of Circular Moonpools［C］// Proceedings of the 6th International Offshore and Polar Engineering Conference. USA：ISOPE， 1996.

我院主持召開首屆國際噴水推進技術研討會

日前，中國船舶海洋及工程研究設計院（MARIC,以下簡稱“我院”）聯(lián)合噴水推進技術重點實驗室（JETLAB）、英國皇家造船師學會（RINA）和上海市船舶與海洋工程學會（SSNAOE）舉辦了首屆國際噴水推進技術研討會（IWWP2015），來自芬蘭、伊朗、俄羅斯、瑞典、英國和中國的20家單位的48位代表參加了會議。

會上，我院副院長兼JETLAB主任陳剛致歡迎辭；RINA首席執(zhí)行官T. Blakeley先生和SSNAOE副理事長繆國平教授分別主持了專題研討會；瑞典Chalmers科技大學A. Eslamdoost博士作主題報告，詳細介紹了噴水推進裝置與船體相互干擾水動力的最新研究進展。其他10位國內(nèi)外代表分別介紹了噴水推進裝置設計、分析、控制和應用等方面的最新成果。與會代表討論熱烈，取得了令人滿意的效果。我院基礎研究部和船舶特種推進事業(yè)部通過近半年的精心組織，順利完成了征文、評審、出版、邀請專家、接待代表和安排會議等工作，獲得了與會代表的好評。會議聯(lián)合主辦單位決定每兩年舉辦一次該項會議，使之逐步成為具有國際影響力的系列學術會議。

Motion performance of drillship with moonpool

ZHANG Xiao-yu1HU Kai-ye2ZHOU Wen-wen2
(1. Ocean Engineering Equipment Research Institute of CNPC Bohai Equipment Research Institute, Panjin 124010, China;
2. College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China)

Abstract:Great stability is demanded for a dillship during the offshore drilling operation. However, the irregular oscillation of the fluid cylinder of the moonpool and the elongated hull shape have adverse effect on the motion performance of the drillship, which may limit the drillship operation adaptability in rough sea states. The oscillation of a drillship with a moonpool is predicted by the three-dimension potential theory, and the results are validated by the ship model tests in order to put forward the criterion value of the drillship oscillation. On this basis, it analyzes the drillship operation adaptability in deep sea and the influence of the scale and shape of the moonpool on the motion performance of the drillship.

Keywords:drillship; motion prediction; moonpool

［作者簡介］張曉宇（1972-）,男，碩士,高級工程師，研究方向：海工裝備。

［收稿日期］2015-01-30；［修回日期］2015-11-18

［基金項目］中國石油天然氣集團公司攻關課題（2011B-1064）

［中圖分類號］U661.1

［文獻標志碼］A

［文章編號］1001-9855（2016）01-0029-08