慢漂運動下多成分錨泊線阻尼參數研究

楊 勇，陳佳欣，黃 山

(1．上海外高橋造船有限公司，上海200137;2．上海外高橋造船海洋工程設計有限公司，上海200137;3．上海振華重工(集團)股份有限公司海上重工設計研究院，上海200125;4．哈爾濱工程大學船舶工程學院，黑龍江哈爾濱150001)

0 引言

海洋浮式平臺的定位性能與運動控制是保證油氣安全開采與生產的關鍵因素之一，與動力定位系統(tǒng)相比，傳統(tǒng)的錨泊定位系統(tǒng)環(huán)境適應性好、可靠性高，在當今世界范圍內占據主導地位。某些浮式鉆井平臺由于作業(yè)需要，其水平位移一般不能超過水深的5%［1］，但是浮式平臺由于受到二階波浪力的作用會產生大位移的慢漂運動，因此最大能夠占到系統(tǒng)總阻尼80%的錨泊阻尼［2］對于消耗系統(tǒng)能量、限制運動幅度的作用十分關鍵。Ormberg［3］對浮式平臺進行耦合分析，發(fā)現(xiàn)錨泊阻尼是系統(tǒng)總阻尼的主要成分，錨泊阻尼對預判平臺的運動響應起著重要的作用，應用全耦合分析法將錨泊阻尼效應考慮進整個系統(tǒng)的運動響應分析中是十分必要。Matsumoto［4］通過比較錨泊阻尼和其他阻尼后發(fā)現(xiàn)，在低頻運動范圍內錨泊阻尼是系統(tǒng)總阻尼的主要成分。Johanning［5］總結了錨泊阻尼的所有數值計算方法，分別是基于錨泊線靜力分析的準靜態(tài)法，基于錨泊線動力分析的時域全動態(tài)有限元法，以及模型試驗法。Kitney［6］應用模型試驗法研究了單一成分錨泊線的阻尼特性。Webster［7］對單一成分錨泊線在上部結構物產生波頻牽引運動時的錨泊阻尼進行了參數研究。喬東生［8］以單一成分錨泊線為研究對象，應用有限元法比較了上部浮體在發(fā)生慢漂運動、波頻運動以及兩者組合運動時的錨泊阻尼。本文通過對更廣泛應用于深水平臺的多成分錨泊線進行阻尼參數研究，得到各個參數對錨泊阻尼的影響，以及水動力系數對錨泊線各成分阻尼的影響，所得結論對浮式平臺的錨泊系統(tǒng)設計有很大的參考意義。

1 錨泊阻尼定義與計算原理

錨泊阻尼主要由海床與錨泊線摩擦阻尼、錨泊線內阻尼和錨泊線流體動力阻尼3部分組成。海床與錨泊線發(fā)生摩擦會產生一部分的系統(tǒng)能量消耗，錨泊線內阻尼由錨泊線內部材料決定，但是這兩部分阻尼成分較小，錨泊線流體動力阻尼由作用在錨泊線上的拖曳力引起，這一部分阻尼是錨泊阻尼的主要部分［8］，本文研究重點也是錨泊線流體動力阻尼。

系泊海洋浮式結構物的水平運動方程為:

式中:M為上部結構物質量;Madd為附加水重量;A為運動加速度;Bstr為上部結構物線性阻尼系數;Bmoor為錨泊線性阻尼系數;V為運動速度;K(X)為結構剛度矩陣;X為運動位移;F為運動力。

錨泊線頂端水平張力瞬時值Tx與錨泊線性阻尼系數Bmoor的關系為:

假設錨泊線頂端導纜孔處運動時程與浮式結構物運動時程相同，波浪作用下浮式結構物運動響應為正弦運動，設為X(t)=Aosin(ωt)，其中Ao為運動振幅，運動周期為τ。單根錨泊線在頂端結構物一個正弦振蕩周期內所消耗的系統(tǒng)能量E可通過下式進行推導:

因此，錨泊線性阻尼系數Bmoor可表示為:

錨泊線所消耗的系統(tǒng)能量E需要借助于有限元軟件計算所得結果，再結合指示圖相關理論來進行求解，具體原理如圖1所示。在軟件中定義錨泊線頂端水平位移Sx隨時間的變化曲線，即輸入錨泊線頂端慢漂運動方程，根據有限元動力計算結果，輸出錨泊線頂端水平張力瞬時值Tx隨時間的變化曲線，并繪制錨泊線頂端水平張力－水平位移曲線，該曲線所圍成的面積，即錨泊線在浮式結構物一個縱蕩振蕩周期內所耗散的系統(tǒng)能量E。

圖1 錨泊線頂端水平張力－水平位移指示圖Fig.1 Indicator diagram for horizontal tension-horizontal displacement of mooring line top end

2 計算模型與思路

以單根多成分錨泊線為研究對象，頂部段與底部段為鋼鏈，中間段為鋼索，設計水深400 m，錨泊線屬性如表1所示，錨泊線計算模型如圖2所示。

圖2 研究模型Fig.2 Model for research

表1 錨泊線屬性Tab.1 Properties of the mooring line

利用全動態(tài)非線性動力計算軟件OrcaFlex進行數值分析，該軟件被國外海洋工程界廣泛的應用于錨系與立管的水動力分析中。許多文獻證實［5，9］，該軟件所得結果與試驗所得結果十分相近，可靠性很高。錨泊線頂端牽引運動為浮式結構物主要水平運動縱蕩運動，通過運行軟件進行時域非線性動力計算，繪制指示圖得出耗散能量，并計算出錨泊線性阻尼系數Bmoor。

3 計算結果與分析

3.1 海底摩擦系數的影響

考慮到錨泊線流體動力阻尼和海床摩擦阻尼都與錨泊線的頂端運動方程相關，為了進一步理清兩者的比例關系，通過表2列舉的設計工況研究了存在海底摩擦力和沒有海底摩擦力時的錨泊阻尼。從表中結果可以看出，兩者之間的結果相差很小，海底摩擦力系數選取0.5時，所有工況的錨泊阻尼增加均不超過10%，這一結論與文獻［10］研究成果也相一致，再一次驗證了流體動力阻尼是錨泊阻尼的主要成分。在本文后續(xù)研究中，海底摩擦系數均設置為0.5。

表2 海底摩擦系數計算參數與結果Tab.2 Calculation parameters and results for variation of seabed friction coefficient

3.2 慢漂運動振幅與周期的影響

表3設計工況研究了錨泊線頂端結構物牽引運動振幅與周期對錨泊阻尼的影響。以300 s為周期不變，頂端運動振幅從30 m到50 m進行變化，以30 m為振幅不變，頂端運動周期從300 s到100 s進行變化。圖3(a)和圖3(b)分別描述了運動振幅變化對錨泊阻尼影響的指示圖和趨勢圖，圖4(a)和圖4(b)分別描述了運動周期變化對錨泊阻尼影響的指示圖和趨勢圖。

從圖3可發(fā)現(xiàn)，隨著慢漂運動振幅的增加，錨泊線所消耗的能量增加，錨泊阻尼呈現(xiàn)出近似線性上升趨勢。從圖4可發(fā)現(xiàn)，隨著慢漂運動周期的增加，錨泊線所消耗的能量減小，錨泊阻尼呈現(xiàn)出非線性下降趨勢。

表3 慢漂運動振幅與周期計算參數與結果Tab.3 Calculation parameters and results for variation of amplitude and period of slow drift motion

圖3 慢漂運動振幅的影響Fig.3 Effect of amplitude of slow drift motion

圖4 慢漂運動周期的影響Fig.4 Effect of period of slow drift motion

3.3 流速大小與分布形式的影響

選擇2組慢漂運動，海底流速Cspeed分別選取0 m/s，0.5 m/s，1.0 m/s和 1.5 m/s，研究流速變化對錨泊阻尼的影響，設計工況如表4所示。同時，針對其中一組慢漂運動研究流速分布形式Cprofile對錨泊阻尼的影響，流速分布形式沿水深H變化如圖5所示，平均流速均為1.0 m/s。

從圖6(a)可知，隨著流速的增加，錨泊阻尼呈現(xiàn)出非線性上升趨勢，無流速情況下2組慢漂運動所產生的錨泊阻尼值相差較大，但隨著流速的不斷增加，兩者之間越來越相近，流速成為決定錨泊阻尼大小的主導因素。從圖6(b)可看出，海流均速相同時，不同流速分布形式對錨泊阻尼的影響十分的小，可忽略。

表4 流速大小與分布形式計算參數與結果Tab.4 Calculation parameters and results for variation of current speed and current profile

圖5 流速分布形式Fig.5 Current profiles

圖6 流速大小與分布形式的影響Fig.6 Effect of current speed and current profile

3.4 垂向拖曳力系數的影響

鋼鏈的垂向拖曳力系數分別選取1.1，1.6，2.4和3.2，鋼索的垂向拖曳力系數分別選取1.2，1.6，2.4和3.2，出于研究目的所取數值并不完全在規(guī)范［11］推薦取值范圍之內。本文研究錨泊線是由3個部分組成，但是頂部段R4鋼鏈長度僅為10 m，驗證后對錨泊阻尼影響不大，故不在此給出具體結果。表5設計工況將對剩余兩部分垂向拖曳力系數開展參數研究，其中S2代表中間段螺旋鋼絲索，S3代表底部段R3S鋼鏈。

從圖7可以看出，隨著拖曳力系數的增加，錨泊阻尼呈現(xiàn)出近乎線性的增加趨勢，并且底部段R3S鋼鏈的拖曳力系數改變對錨泊阻尼的影響較大。導致這一現(xiàn)象產生的原因主要有2點:鋼鏈的拖曳直徑要大于鋼索的拖曳直徑，因此利用莫里森方程計算拖曳力時，鋼鏈拖曳系數的改變導致作用在其上的拖曳力改變更大;鋼索的運動響應位移與鋼鏈相比也較小，如圖8所示，劃線上部為鋼索 (除頂端有10 m鋼鏈)，下部為鋼鏈。

表5 垂向拖曳力系數計算參數與結果Tab.5 Calculation parameters and results for variation of normal drag coefficient

圖7 垂向拖曳力系數的影響Fig.7 Effect of normal drag coefficient

圖8 錨泊線法向位移比較Fig.8 Comparison of transverse displacement of the line element

3.5 垂向附加質量系數的影響

表6設計工況研究了錨泊線垂向附加質量系數改變對錨泊阻尼的影響，從表中所列結果可看出，無論是改變中間段螺旋鋼絲索還是底部段R3S鋼鏈的垂向附加質量系數，錨泊阻尼值幾乎不變，該系數對錨泊阻尼影響不大。

表6 垂向附加質量系數計算參數與結果Tab.6 Calculation parameters and results for variation of normal added mass coefficient

4 結語

通過對一根應用于實際工程的多成分懸鏈式錨泊線進行時域非線性動力分析，計算能量耗散與錨泊阻尼值，可得出以下結論:

1)海底摩擦力所引起的海床摩擦阻尼與拖曳力引起的流體動力阻尼相比很小，存在海底摩擦力時，錨泊阻尼大約增加5%至7%;

2)隨著慢漂運動振幅的增加亦或是周期的減小，錨泊線消耗能量上升，錨泊阻尼值變大;

3)隨著流速的增加，錨泊阻尼值呈現(xiàn)非線性上升趨勢，并且成為主導錨泊阻尼值大小的決定性因素，在海流均速相同的情況下，不同流速分布形式對錨泊阻尼值的影響很小;

4)錨泊線垂向拖曳力系數的改變對錨泊阻尼值的影響很大，隨著該系數的增加，錨泊阻尼值呈現(xiàn)出近似線性上升趨勢，并且底部段R3S鋼鏈的垂向拖曳力系數改變對錨泊阻尼值影響更大;

5)錨泊線垂向附加質量系數的改變對錨泊阻尼值的影響很小。

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

［2］ HUSE E，MATSUMOTO K．Mooring line damping due to first-and second-order vessel motion［C］//Offshore Technology Conference．Houston，1989:135 －148．

［3］ ORMBERG H，LARSEN K．Coupled analysis of floater motion and mooring dynamics for a turret－moored ship［J］．Applied Ocean Research，1998，20:55 －67．

［4］ MATSUMOTO K．The influence of mooring line damping on the prediction of low-frequency vessel motions at sea［C］//Offshore Technology Conference．Houston，1991:119－128．

［5］ JOHANNING L，SMITH G H，WOLFRAM J．Measurements of static and dynamic mooing line damping and their importance for floating WEC devices［J］．Ocean Engineering，2007，34:1918 －1934．

［6］ KITNEY N，BROWN D T．Experimental investigation of mooring line loading using large and small-scale models［J］．Journal of Offshore Mechanics and Arctic Engineering，2001，123:1 －9．

［7］ WEBSTER W C．Mooring-induced damping［J］．Ocean Engineering，1995，22(6):571 －591．

［8］ 喬東升，歐進萍．深水懸鏈錨泊線粘性阻尼計算［J］．海洋工程，2009，27(4):16 －22．QIAO Dong-sheng，OU Jin-ping．Calculation on viscous damping of deepwater catenary mooring line［J］．The Ocean Engineering，2009，27(4):16 －22．

［9］ BROWN D T，MAVRAKOSS．Comparative study on mooing line dynamic loading［J］．Marine Structures，1999(12):131－151．

［10］喬東升，歐進萍．深水懸鏈錨泊線阻尼計算［J］．振動與沖擊，2011，30(2):24 －31．QIAO Dong-sheng，OU Jin-ping．Damping calculation of a deepwater catenary mooring line［J］．Journal of Vibration and Shock，2011，30(2):24 －31．

［11］ DNV．Environmental conditions and environmental loads(DNV －RP－C205)［S］．Norway，2010:52－62．