深水鉆井船錨泊定位分析方法研究

魏躍峰 潘方豪（中國船舶及海洋工程設計研究院 上海200011）

特邀專欄

深水鉆井船錨泊定位分析方法研究

魏躍峰 潘方豪
（中國船舶及海洋工程設計研究院 上海200011）

采用準動態方法和動態耦合方法對作業于1 500 m水深的鉆井船開展錨泊定位分析，比較兩種方法的結果，并研究兩種錨泊定位分析方法的特點。研究表明：不同設計階段采用不同的錨泊定位分析方法能夠提高設計的效率和精度。

深水鉆井船；錨泊定位；準動態；動態耦合

引　言

深水鉆井船的主要功能是進行海上油氣勘探、鉆井作業。為保護鉆井隔水管及下部的球形接頭不受損壞，鉆井船的船體運動應限制在一定的范圍內，這對深水鉆井船的定位性能提出了很高的要求［1］。鉆井船的定位方式有動力定位［2］和錨泊定位［3-4］兩種。動力定位鉆井船具有快速作業、快速撤離的優勢，但運營成本較高。作業于相對溫和海域的鉆井船，可以采用錨泊定位方式。

錨泊定位分析方法是深水鉆井船錨泊定位系統設計的關鍵。本文分別采用準動態方法和動態耦合分析方法對深水鉆井船開展錨泊定位分析，計算完整系泊和單根錨索破斷時鉆井船錨泊系統的定位能力，研究兩種不同方法的特點及在錨泊定位系統設計不同階段的作用。

1　理論方法

1.1 準動態法

準動態方法首先根據三維頻域勢流理論計算船體水動力系數、一階波浪力、二階波浪力和運動幅值響應傳遞函數（RAOs），然后求解船體時域運動。在求解船體時域運動時，每一個時間步內的低頻運動（縱蕩、橫蕩和首搖）和波頻運動（垂蕩、橫搖和縱搖）均需分別計算。低頻運動響應通過求解時域運動方程得到，船體低頻運動求解方程見式（1）［5］：式中：M為船體水平運動方向的質量矩陣；Fwave（2）為船體所受到的二階波浪漂移力；Fcurrent為海流力；Fwind為風力；Fmoor為錨泊力。準動態方法中假定錨索處于靜止狀態，錨泊力通過求解靜力學方程得到［6］。

在低頻運動位置和船體首搖角的基礎上，根據頻域計算得到運動幅值響應結果，結合波浪時歷獲得波頻運動響應時歷。

準動態方法雖考慮到船體附加質量和二階波浪漂移力，但忽略了船體阻尼、一階波浪力以及錨索自身拖曳力和慣性力等諸多影響。

1.2 動態耦合

動態耦合分析方法首先通過三維頻域勢流理論計算船體水動力系數、一階波浪力和二階波浪力；然后通過快速Fourier變換（FFT）將頻域分析得到的結果轉換為相應的時間序列；最后結合每一個時間步內算得的錨泊力，求解船體時域運動方程。時域耦合運動方程為［7］：

式中：M為船體六自由度運動方向的質量矩陣；μ為附加質量；K（）為延遲函數；Fwave（1）為船體所受到的一階波浪力；Fwave（2）為船體所受到的二階波浪漂移力；Fcurrent為海流力；Fwind為風力；Fmoor為系泊力。動態耦合方法中，錨泊力通過求解動力學方程得到［7］。

動態耦合分析不僅考慮船體水動力系數、一階波浪力、二階波浪漂移力以及錨索自身拖曳力和慣性力的影響，而且還考慮每一個時間步船體運動對錨索張力的影響。

2　深水鉆井船及錨泊定位系統

本文研究的深水鉆井船船體主尺度參數見表1。

表1　鉆井船主尺度參數

圖1　錨泊系統布置圖

鉆井船采用8根錨索（分為4組、每組2根），這2根錨索與船體中線的夾角分別為40°和60°。錨泊系統布置如圖1所示。每一根錨索由鋼絲繩和錨鏈組成。與船體相連的一段為鋼絲繩，與錨相連的一段為錨鏈。錨索參數見表2。

表2　錨索參數

3　環境條件

深水鉆井船作業水深1 500 m，作業海域環境條件見表3。在進行錨泊定位分析計算時，波浪譜采用JONSWAP譜，風采用API風譜，流采用定常流。對于多點系泊，當風、浪、流等環境載荷方向相同時，錨索受力最大，因此本文計算時也設定風、浪、流、方向相同，見圖1。

表3　設計環境條件

4　結果與分析

文中分別采用BV船級社基于準動態方法開發的ARIANE軟件以及DNV船級社基于動態耦合分析方法開發的DeepC軟件計算了一年一遇季風和一年一遇臺風情況，分別在迎浪、首斜浪和橫浪下，完整系泊和單根錨索（第四根）的破斷工況，以及鉆井船運動偏移、導纜孔處錨索張力和錨泊點處錨索張力。本文僅以首斜浪（風、浪、流同向且為157.5°）為例，給出計算結果，其他方向結論與之相同。

圖2　作業工況完整系泊時結果比較

圖3　作業工況單根錨索破斷時結果比較

圖2 — 圖5給出了采用準動態方法和動態耦合方法開展的深水鉆井船錨泊定位分析結果比較。圖中導纜孔處張力和錨點張力均為第三根錨索受力。從圖中可以看出，準動態方法計算出的船體運動偏移較動態耦合分析方法偏大，導纜孔處和錨點處的錨索張力偏小。這主要是由于耦合的影響：準動態方法每一時間步長內計算出的錨索張力作用到船體運動上，船體運動未引起錨索張力變化；動態耦合方法在每一時間步長內考慮兩者之間的相互影響，反復迭代，直至結果收斂。上述四種工況下的計算結果如下：一年一遇季風海況，完整系泊時，準動態方法和動態耦合分析方法得到的船體運動偏移相差最大（為15.54%）；一年一遇臺風海況，單根錨索破斷時，船體運動偏移相差最小（為2.51%）；一年一遇臺風海況，完整系泊時，準動態方法和動態耦合分析方法得到的導纜孔處錨索張力相差最大（為23.46%）；一年一遇季風海況，單根錨索破斷時，導纜孔處錨索張力相差最小（為2.41%）；一年一遇臺風海況，完整系泊時，準動態方法和動態耦合分析方法得到的錨點處錨索張力相差最大（為34.16%）；一年一遇季風海況，單根錨索破斷時，錨點處錨索張力相差最小（為6.61%）。

圖4　生存工況完整系泊時結果比較

圖5　生存工況單根錨索破斷時結果比較

準動態方法在計算時忽略了船體阻尼系數和錨索動力等的影響，而動態耦合分析方法均考慮了這些因素的影響，因此動態耦合分析方法計算出的結果更接近真實結果。在工程應用中，考慮到準動態方法和動態耦合方法計算結果的差別，BV船級社規范NR493［8］規定，準動態方法錨索張力安全系數不小于1.75，動態方法安全系數不小于1.67。由于動態耦合分析方法計算耗時較長，在進行錨泊系統設計的初始階段，通常采用準動態方法進行計算，得出對錨泊系統影響最大的惡劣海況，再對惡劣海況采用動態耦合分析的方法計算。

5　結　論

本文采用準動態和動態耦合兩種方法對作業于1 500 m水深的鉆井船開展錨泊定位分析，分別計算了一年一遇的季風和一年一遇的臺風情況下，完整錨泊和單根錨索破斷時的鉆井船船體運動偏移、導纜孔處錨索張力和錨點處的錨索張力。通過對兩種方法的計算結果進行比較，研究得出各自的特點：準動態方法同動態耦合分析方法相比，具有一定的精度。在深水鉆井船錨泊系統設計的初始階段可采用其對各種工況開展錨泊定位分析，而后再對惡劣海況采用動態耦合分析方法計算，從而高效準確地設計出合理的系泊系統。

［1］ API. Recommended Practice 2SK for Design and Analysis of Stationkeeping Systems for Floating Structures［S］.2005.

［2］ 何進輝，張海彬.深水鉆井船動力定位能力分析方法研究［J］.海洋工程.2014（1）： 25-31.

［3］ 麻翠榮.深水鉆井船船型開發及總體性能分析［D］.哈爾濱：哈爾濱工程大學，2010.

［4］ 曲媛.深海鉆井船的總體方案分析研究［D］.哈爾濱：哈爾濱工程大學，2010.

［5］ BV. Ariane 7 Theoretical Manual ［S］. 2007.

［6］ 黃祥鹿，陸鑫森. 海洋工程流體力學及結構動力響應［M］. 上海：上海交通大學出版社，1992.

［7］ Norwegin Marine Technology Research Institule. Marintek Report：SIMO-Theory Manual Version 3.4［R］. 2004.

［8］ BV. Classification of Mooring Systems for Performanent Offshore Units［S］. 2012.

On mooring position analysis method for deep-water drilling ship

WEI Yue-feng PAN Fang-hao
(Marine Design & Research Institute of China, Shanghai 200011, China)

The quasi-dynamic method and the dynamic coupled method are adopted to perform mooring position analysis of a deep-water drilling ship, which operates in the water depth of 1 500 m. The characteristics of the two methods are studied by comparison between the two methods. It shows that the efficiency and the precision of the mooring system design can be improved by the application of the different mooring position methods in the relevant design stage.

deep-water drilling ship; mooring position; quasi-dynamic; dynamic coupled

U653.2

A

1001-9855（2016）04-0001-04

10.19423 / j.cnki.31-1561 / u.2016.04.001

2016-03-16；

2016-04-29

魏躍峰（1981-），男，博士，工程師，研究方向：船舶及海洋工程水動力性能及系泊定位。

潘方豪（1978-），男，研究員，研究方向：船舶及海洋工程舾裝研究設計。