超深水鉆井船U型減搖水艙設(shè)計(jì)*

張 琪 張利軍 曹 凱 邢永強(qiáng) 萬(wàn)家平

(中遠(yuǎn)海運(yùn)重工有限公司技術(shù)研發(fā)中心1) 大連 116600) (啟東中遠(yuǎn)海運(yùn)海洋工程有限公司海工研發(fā)中心2) 啟東 226200)

0 引 言

超深水鉆井船在作業(yè)工況下需要具備良好的耐波性，以保證船舶在惡劣海況下的鉆井作業(yè)能力.減搖水艙作為一種十分常見(jiàn)的減搖裝置，將其應(yīng)用在鉆井船上有助于改善鉆井船在波浪中的橫搖特性.根據(jù)工作原理不同，減搖水艙主要分為主動(dòng)式和被動(dòng)式兩大類.主動(dòng)式根據(jù)橫搖情況借助泵將水從一舷打向另一舷，以減小船舶的橫搖，雖然減搖效果好，但由于使用時(shí)功率消耗過(guò)大，目前一般只用于抗傾.被動(dòng)式通過(guò)結(jié)構(gòu)設(shè)計(jì)使艙內(nèi)水振蕩的固有周期接近于船舶的橫搖固有周期，依據(jù)雙共振原理實(shí)現(xiàn)減搖，這種方式結(jié)構(gòu)簡(jiǎn)單、造價(jià)低且維修方便，因此，應(yīng)用更為廣泛[1-3].在被動(dòng)式減搖水艙的設(shè)計(jì)方法中，采用經(jīng)驗(yàn)公式評(píng)估減搖水艙的固有周期往往準(zhǔn)確度不高，容易造成減搖效果不佳，采用臺(tái)架試驗(yàn)或者船模試驗(yàn)來(lái)進(jìn)行設(shè)計(jì)能夠有效地保證結(jié)果的準(zhǔn)確度和可靠性，但時(shí)間成本和花費(fèi)較高，隨著CFD數(shù)值模擬技術(shù)的高速發(fā)展，CFD技術(shù)在減搖水艙研究上應(yīng)用價(jià)值也逐漸顯現(xiàn)[4-7].文中提出了將減搖水艙概念應(yīng)用于鉆井船設(shè)計(jì)，通過(guò)CFD方法對(duì)超深水鉆井船減搖水艙的外形尺寸進(jìn)行優(yōu)化，并將CFD計(jì)算結(jié)果與臺(tái)架試驗(yàn)結(jié)果進(jìn)行了比較分析.

1 減搖水艙布置方案及主尺度確定

結(jié)合超深水鉆井船的布置特點(diǎn)，在減搖水艙位置的選擇上主要考慮了主甲板以上靠近船艉的位置和船舯月池的前方的艙室內(nèi)(高度在主甲板與燃油艙之間)兩種布置方案.第一種方案不占用船艙內(nèi)部空間，位于靠近船艉的位置，由于距離船舯較遠(yuǎn)影響減搖效果，且水艙的位置過(guò)高導(dǎo)致船舶重心的升高，對(duì)船舶穩(wěn)性不利；第二種方案雖然會(huì)損失一部分燃油艙的位置，但沿縱向位于船舯1/3船長(zhǎng)范圍內(nèi)，減搖水艙高度方向上布置在船的重心附近，位置更為合理[8]，因此，選擇第二種布置方案作為最終布置方案.在總布置圖中標(biāo)記的減搖水艙布置位置見(jiàn)圖1.

圖1 減搖水艙布置方案示意圖

選定的減搖水艙布置區(qū)域最大的空間尺寸為36.4 m×5.6 m×4.9 m(長(zhǎng)×寬×高)，占鉆井船(鉆井工況下)排水體積的1.1%，為充分利用有限的布置空間，減搖水艙的長(zhǎng)、寬、高與空間內(nèi)的最大尺度保持一致.

2 被動(dòng)式減搖水艙設(shè)計(jì)

2.1 U型減搖水艙設(shè)計(jì)原理

U型水艙通過(guò)艙內(nèi)水左右振蕩產(chǎn)生的周期性力矩實(shí)現(xiàn)減搖.船在波浪力矩的作用下發(fā)生橫搖，當(dāng)波浪力矩頻率與船舶固有頻率接近時(shí)，船舶橫搖角最大且滯后波浪力矩大約90°.同理，船舶橫搖帶動(dòng)艙內(nèi)水振蕩產(chǎn)生減搖力矩，當(dāng)船舶橫搖頻率與水艙固有頻率接近時(shí)，水艙產(chǎn)生的減搖力矩最大且滯后船舶橫搖大約90°.波浪力矩、船舶橫搖角、減搖力矩三者之間相位關(guān)系見(jiàn)圖2.

圖2 波浪力矩、船舶橫搖及減搖力矩的相位關(guān)系

上述物理現(xiàn)象稱之為雙共振，此時(shí)水艙減搖力矩最大且與波浪力矩反向，減搖效果最佳.采用水動(dòng)力計(jì)算軟件對(duì)超深水鉆井船鉆井和儲(chǔ)油兩種典型工況的固有周期進(jìn)行計(jì)算，得到的固有周期計(jì)算結(jié)果見(jiàn)表1.

表1 超深水鉆井船工作狀態(tài)下的固有周期計(jì)算結(jié)果

由表1可知，船舶工作狀態(tài)下的固有周期范圍為16.46～18.24 s，為兼顧減搖水艙在不同工況下的減搖效果，減搖水艙的固有周期應(yīng)與船舶的固有周期相接近，因此，取上述工況下固有周期的平均值17.3 s作為減搖水艙的目標(biāo)固有周期.

2.2 基于經(jīng)驗(yàn)公式的減搖水艙設(shè)計(jì)

依據(jù)雙共振原理設(shè)計(jì)U型水艙的結(jié)構(gòu)外形，考慮到減搖水艙的總高度不少于液位高度1.7倍的設(shè)計(jì)原則[9]，水艙內(nèi)的液位應(yīng)不高于2.9 m.

圖3所示為U型水艙的典型剖面及相關(guān)尺寸，固有周期的計(jì)算為

(1)

圖3 U型水艙的典型剖面及相關(guān)尺寸

根據(jù)渡邊公式，Le=2h1+2(l+r)·A0/A1，其中：h0為水艙內(nèi)的液位高度；A0為邊水艙的斷面積；A1為連通水道的斷面積.已知H=4.9 m，B=36.4 m，h1=h0-0.5a，這里液位高度取最大，即h0=2.9 m.假定l=7 m則r=5.6 m，根據(jù)經(jīng)驗(yàn)公式計(jì)算可知當(dāng)a=2 m時(shí)，該設(shè)計(jì)方案的固有周期為T=17.1 s與目標(biāo)固有周期17.3 s較為接近.

2.3 基于數(shù)值模擬方法的減搖水艙設(shè)計(jì)

2.3.1數(shù)值模擬方法

數(shù)值計(jì)算統(tǒng)一采用RANS模型框架下的SSTk-ω模型，SSTk-ω湍流模型是融合了k-ω湍流模型與k-ε湍流模型兩種湍流模型的理論構(gòu)建的，并且考慮了湍流剪切力的輸運(yùn)效應(yīng)，可以精準(zhǔn)地計(jì)算出逆向壓力梯度引起的流動(dòng)分離位置和作用區(qū)域.

自由液面的變化采用VOF法來(lái)捕捉，假設(shè)計(jì)算區(qū)域是V，將其分成兩個(gè)區(qū)域，流體A所在的區(qū)域記為V1，而流體B所在的區(qū)域記為V2.定義這樣一個(gè)函數(shù)(表征流體質(zhì)點(diǎn)所在的區(qū)域)：

(2)

對(duì)于由兩種不相容的流體組成的流場(chǎng)，α(x,t)滿足：

(3)

式中：U為流體的速度場(chǎng).

在每個(gè)網(wǎng)格Iij上定義Cij，Cij為α(x,t)在網(wǎng)格上的積分.

(4)

該函數(shù)稱之為VOF函數(shù)，若數(shù)值在0和1之間，那么該網(wǎng)格為自由液面處的網(wǎng)格，一系列這樣的網(wǎng)格單元組成了自由液面.

2.3.2網(wǎng)格的無(wú)關(guān)性驗(yàn)證

以基于經(jīng)驗(yàn)公式得到的設(shè)計(jì)方案(l=7 m，r=5.6 m，a=2 m)為研究對(duì)象，對(duì)縮尺比為20的幾何相似模型固有周期進(jìn)行CFD數(shù)值模擬，為了研究網(wǎng)格的疏密程度對(duì)計(jì)算結(jié)果的影響，通過(guò)調(diào)整網(wǎng)格的基礎(chǔ)尺寸分別劃分了網(wǎng)格數(shù)81萬(wàn)、120萬(wàn)和162萬(wàn)三組網(wǎng)格(網(wǎng)格示意圖以網(wǎng)格數(shù)120萬(wàn)為例).對(duì)水艙施加一個(gè)外界的擾動(dòng)(沿橫搖方向)作為初始條件，待水艙恢復(fù)到平衡位置停止外界干擾，使水艙內(nèi)的水自由衰減.通過(guò)監(jiān)測(cè)水艙在自由衰減過(guò)程中所承受力矩的衰減曲線(見(jiàn)圖4)來(lái)讀取水艙的固有周期，再將模型尺度計(jì)算結(jié)果換算為實(shí)尺度計(jì)算結(jié)果.不同網(wǎng)格條件下的計(jì)算結(jié)果見(jiàn)表2.

圖4 減搖水艙自由衰減曲線

表2 不同網(wǎng)格下的固有周期計(jì)算結(jié)果

由表2可見(jiàn)，三組不同網(wǎng)格得到的計(jì)算結(jié)果非常接近，差異度均小于0.1%，說(shuō)明網(wǎng)格的增加對(duì)計(jì)算結(jié)果幾乎沒(méi)有影響，81萬(wàn)網(wǎng)格數(shù)的計(jì)算精度是足夠的.而將固有周期數(shù)值模擬計(jì)算結(jié)果與經(jīng)驗(yàn)公式結(jié)果進(jìn)行對(duì)比，數(shù)值模擬結(jié)果較經(jīng)驗(yàn)公式結(jié)果約小10%，兩者相差較大.

2.3.3減搖水艙固有周期的數(shù)值計(jì)算

CFD計(jì)算結(jié)果表明目前設(shè)計(jì)的減搖水艙固有周期較目標(biāo)固有周期偏小，因此需通過(guò)改變U型水艙的外形尺寸來(lái)增大水艙的固有周期.觀察經(jīng)驗(yàn)式(1)可知A0/A1的大小對(duì)U型水艙的固有周期影響最大，因此通過(guò)降低a的高度增大固有周期是可行的.改變a的取值重新計(jì)算減搖水艙的固有周期，表3中所列經(jīng)過(guò)兩次計(jì)算最終設(shè)計(jì)方案3(a=1.5 m)固有周期為17.19 s，與目標(biāo)值相差0.64%十分接近.減搖水艙設(shè)計(jì)方案3在計(jì)算過(guò)程中水氣分界面的示意圖見(jiàn)圖5.

表3 各設(shè)計(jì)方案的結(jié)構(gòu)尺度及固有周期比較

圖5 減搖水艙自由衰減過(guò)程中水氣分界面示意圖

3 臺(tái)架試驗(yàn)

3.1 臺(tái)架試驗(yàn)?zāi)Ｐ图霸囼?yàn)裝置

參照減搖水艙設(shè)計(jì)方案3的結(jié)構(gòu)尺寸制作有機(jī)玻璃試驗(yàn)?zāi)Ｐ?，模型的縮尺比為20，加工好的水艙模型見(jiàn)圖6.試驗(yàn)在六自由度晃蕩平臺(tái)系統(tǒng)中進(jìn)行，該試驗(yàn)裝置采用全電動(dòng)動(dòng)感可調(diào)墩平臺(tái)，又名六自由度運(yùn)動(dòng)仿真平臺(tái)(見(jiàn)圖7)，它由六自由度運(yùn)動(dòng)平臺(tái)以及平臺(tái)控制系統(tǒng)組成，六自由度晃蕩平臺(tái)系統(tǒng)能夠?qū)崿F(xiàn)六自由度激勵(lì).

圖6 減搖水艙試驗(yàn)?zāi)Ｐ?/p>

圖7 在六自由度晃蕩平臺(tái)系統(tǒng)

3.2 臺(tái)架試驗(yàn)結(jié)果

試驗(yàn)內(nèi)容包括減搖水艙自由衰減試驗(yàn)和臺(tái)架強(qiáng)迫振蕩試驗(yàn)，其中自由衰減試驗(yàn)共分為兩個(gè)步驟：①對(duì)水艙進(jìn)行激勵(lì)，使艙內(nèi)液體充分運(yùn)動(dòng)起來(lái)；②中斷激勵(lì)，讓液體自由震蕩，直至停止.對(duì)激勵(lì)過(guò)程進(jìn)行分析，計(jì)算出衰減過(guò)程中液體振蕩的平均周期.

圖8為測(cè)得的減搖水艙設(shè)計(jì)方案3的自由衰減曲線.讀取水艙所受力矩衰減曲線的平均周期，再將水艙模型尺度的固有周期換算為實(shí)尺度的固有周期，換算結(jié)果為17.33 s.試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果非常接近，僅比數(shù)值模擬結(jié)果大0.81%，而采用經(jīng)驗(yàn)公式計(jì)算設(shè)計(jì)方案3的固有周期為19.67 s，試驗(yàn)結(jié)果較其小11.9%，說(shuō)明使用經(jīng)驗(yàn)公式得到的U型水艙固有周期并不準(zhǔn)確.

圖8 減搖水艙自由衰減曲線

減搖水艙的強(qiáng)迫振蕩試驗(yàn)，以鉆井工況和儲(chǔ)油工況1作為典型工況進(jìn)行試驗(yàn)，通過(guò)設(shè)定激勵(lì)頻率及激勵(lì)幅值，測(cè)得單位激勵(lì)幅值下隨遭遇周期變化的減搖水艙減搖效果曲線見(jiàn)圖9.在共振周期附近鉆井工況和儲(chǔ)油工況1下的減搖效果分別能達(dá)到21%和23.5%.

圖9 隨遭遇周期變化的減搖效果曲線

4 結(jié) 論

1) 根據(jù)渡邊公式計(jì)算的U型水艙固有周期較實(shí)際結(jié)果明顯偏大，因此通過(guò)經(jīng)驗(yàn)公式得到的減搖水艙設(shè)計(jì)方案是不準(zhǔn)確的，但可用于確定初始設(shè)計(jì)方案，有利于提高設(shè)計(jì)效率.

2) 通過(guò)CFD數(shù)值模擬方法求解的U型減搖水艙固有周期試驗(yàn)值相差很小不到1%，說(shuō)明CFD數(shù)值模擬方法能夠準(zhǔn)確地預(yù)報(bào)減搖水艙的固有周期.

3) 由于受到目標(biāo)鉆井船布置空間的限制，減搖水艙的裝水量偏少，建議在空間足夠的情況下，保持減搖水艙的橫剖面不變，適當(dāng)增加縱向長(zhǎng)度可以增加減搖水艙的裝水量，進(jìn)而有效地提高減搖水艙的減搖效果.