陳希恰，宋磊建，杜夏英，萬 軍，付世曉

(1. 中海石油(中國)天津分公司渤海石油研究院，天津 300452；2. 上海交通大學(xué)海洋工程國家重點(diǎn)實(shí)驗(yàn)室，上海 200240)

深海臍帶纜總體響應(yīng)特性研究

陳希恰1，宋磊建2，杜夏英1，萬 軍1，付世曉2

(1. 中海石油(中國)天津分公司渤海石油研究院，天津 300452；2. 上海交通大學(xué)海洋工程國家重點(diǎn)實(shí)驗(yàn)室，上海 200240)

對(duì)某深海臍帶纜的總體布局進(jìn)行了選型設(shè)計(jì)研究。臍帶纜的設(shè)計(jì)研究正逐漸成為國內(nèi)海洋工程領(lǐng)域的研究熱點(diǎn)，對(duì)于臍帶纜設(shè)計(jì)而言總體響應(yīng)分析是其關(guān)鍵課題。首先推導(dǎo)了懸鏈線和緩波形兩種布局的線型方程。進(jìn)而基于推導(dǎo)的方程，在有限元軟件OrcaFlex中建立了臍帶纜懸鏈線和緩波形兩種布局的數(shù)值模型，對(duì)比分析了兩種布局形態(tài)下的有效張力和彎曲曲率等響應(yīng)特性。研究發(fā)現(xiàn)相對(duì)于懸鏈線布局，緩波形布局線型較復(fù)雜，但能有效降低懸掛點(diǎn)處有效張力。浮力塊能有效緩解上部動(dòng)態(tài)響應(yīng)向觸地點(diǎn)處的傳播。此外，對(duì)緩波形布局中浮子段長度、浮力因子和浮子段起始位置等關(guān)鍵參數(shù)對(duì)總體響應(yīng)的影響進(jìn)行了敏感性分析。所得結(jié)論便于設(shè)計(jì)人員理解臍帶纜總體響應(yīng)特性,選擇合適的參數(shù)用于臍帶纜總體設(shè)計(jì)。

深海臍帶纜；總體響應(yīng)；懸鏈線；緩波形；有效張力；彎曲曲率

0 引 言

連接上部浮體和水下生產(chǎn)系統(tǒng)的深海臍帶纜，為水下生產(chǎn)系統(tǒng)提供電能、控制信號(hào)以及化學(xué)藥劑，是深海油氣開發(fā)中的重要裝備。隨著我國海洋油氣開發(fā)不斷向深海邁進(jìn)，臍帶纜的國產(chǎn)化以及相關(guān)基礎(chǔ)研究逐漸成為國內(nèi)海洋工程領(lǐng)域的熱點(diǎn)。

深海臍帶纜所受到的主要載荷包括自重引起的拉伸載荷和波、浪、流及船體運(yùn)動(dòng)引起的動(dòng)態(tài)彎曲載荷。合理的總體布局能夠有效降低臍帶纜受到的載荷。因此，總體布局設(shè)計(jì)是臍帶纜研究工作中的重要課題。在該領(lǐng)域，眾多學(xué)者已經(jīng)開展了大量的研究[1-5]。懸鏈線布局與緩波形布局是臍帶纜常用的兩種布局形式。研究表明，在一般水深的作業(yè)環(huán)境中，由于具有結(jié)構(gòu)簡單、施工方便、成本較低和對(duì)浮體運(yùn)動(dòng)適應(yīng)能力強(qiáng)等特點(diǎn)，懸鏈線布局是臍帶纜布局形式的首選[6-7]。當(dāng)水深增加時(shí)，隨著長度和重量的增加，臍帶纜頂部懸掛點(diǎn)處的張力會(huì)相應(yīng)增大，從而提高了對(duì)頂部張緊器的要求；當(dāng)工作環(huán)境惡劣時(shí)，頂部浮體和波浪運(yùn)動(dòng)的動(dòng)態(tài)干擾將會(huì)顯著地傳播到臍帶纜觸地點(diǎn)處進(jìn)而引發(fā)發(fā)臍帶纜的疲勞損傷。而采用在臍帶纜下半段設(shè)置浮力塊的緩波形布局形式，可以有效地降低頂部懸掛點(diǎn)處的張力，同時(shí)也可緩解上部平臺(tái)運(yùn)動(dòng)和波浪作用對(duì)觸地點(diǎn)處臍帶纜的影響。但相對(duì)于懸鏈線布局，緩波形布局更加復(fù)雜，安裝和制造費(fèi)用更高[8-11]。

為了深入了解深海臍帶纜在不同布局形式下的結(jié)構(gòu)響應(yīng)，本文首先推導(dǎo)了懸鏈線和緩波形兩種布局的線型方程。進(jìn)而基于推導(dǎo)的方程，在有限元軟件OrcaFlex[12]中分別建立了臍帶纜在懸鏈線和緩波形兩種布局形式下的數(shù)值模型。對(duì)比分析了兩種布局形式下臍帶纜有效張力和彎曲曲率的結(jié)構(gòu)響應(yīng)特性。此外，還對(duì)緩波形布局中浮子段長度、浮力因子和浮子段起始位置等關(guān)鍵參數(shù)對(duì)臍帶纜總體響應(yīng)的影響進(jìn)行了敏感性分析。

1 總體布局線型理論推導(dǎo)

臍帶纜空間線型是其總體響應(yīng)分析的基礎(chǔ)。基于錨鏈的懸鏈線方程，可推導(dǎo)獲得臍帶纜懸鏈線布局和緩波形布局理論公式。

1.1 懸鏈線布局

圖1為懸鏈線線型示意圖，其空間坐標(biāo)位置表達(dá)式為

圖1 懸鏈線布局線型示意圖Fig.1 Shape of the free hanging catenary

式中：TH為懸掛點(diǎn)的水平張力；φw為頂部懸掛角；h為懸掛點(diǎn)處水深；W為管線單位長度濕重。

由式(1)和式(2)可知，獲得頂部懸掛角、懸掛點(diǎn)處水深和管線單位長度濕重三個(gè)參數(shù)，即可確定臍帶纜任意一點(diǎn)的位置坐標(biāo)。

1.2 緩波形布局

緩波形布局總體線型分布如圖2。整條臍帶纜可以分為懸垂段AC、浮子段CE、下降段EF以及觸地段FP。其中懸垂段可以分為懸掛段AB和跳接段BC；浮子段可以分為上升段CD和拖曳段DE；下降段為EF。圖2中的S1、S2、S3、S4和S5表示的是各段臍帶纜的長度；x1～x5和y1～y5分別為懸掛段AB、跳接段BC、上升段CD、拖曳段DE和下降段EF的水平跨距和垂直長度。

圖2 緩波形布局線型示意圖Fig.2 Shape of the lazy wave

緩波形布局中，懸掛段AB、跳接段BC、升舉段CD和拖曳段DE以及下降段EF均可看作懸鏈布局。根據(jù)懸鏈線方程可以得到緩波形線形空間坐標(biāo)理論公式[3]：

式中：αB、αD和αF分別為點(diǎn)B、點(diǎn)D和點(diǎn)F處的曲率。

根據(jù)懸垂段AC、浮子段CE、下降段EF的長度S12、S34和S5，各段單位長度濕重W12、W34和W5，可計(jì)算出各段的x、y坐標(biāo)以及αB、αD和αF：

式中：n為浮力因子，其計(jì)算公式為

式中：Wf為CE段浮力塊提供的浮力。

臍帶纜緩波形布局設(shè)計(jì)中，一般給出的參數(shù)包括：頂部懸掛點(diǎn)處水深h，頂端懸掛角θ，管線單位長度的濕重W，懸垂段長度S12，浮子段長度S34以及浮力因子n。下降段的長度S5是未知的。因此，為了得到緩波形的線型布局，需要進(jìn)行迭代計(jì)算[3]。具體計(jì)算過程如下：

(1) 假定S5=S0。

(2) 將S5代入式(4)即可得到αB、αD和αF。

(3) 利用式(5)計(jì)算得到x1～x5和y1～y5。

(4) 根據(jù)公式H=y1+y4+y5-y2-y3得到H。

(5) 當(dāng)|H-h|/h<0.1%時(shí)，終止計(jì)算；否則，令S5=S0+ΔS，重復(fù)步驟(2)～(5)。

將各參數(shù)的計(jì)算結(jié)果代入式(3)即可獲得緩波形線型空間位置坐標(biāo)。

2 數(shù)值模型建立

臍帶纜結(jié)構(gòu)組成由內(nèi)到外層次分明，如圖3所示，外層構(gòu)件以一定的角度螺旋纏繞于內(nèi)層構(gòu)件上，外層與內(nèi)層的構(gòu)件采用非粘結(jié)方式接觸，從而導(dǎo)致臍帶纜具有多層、多接觸、可滑動(dòng)的特性。雖然臍帶纜的結(jié)構(gòu)組成非常復(fù)雜，但其總體力學(xué)行為與結(jié)構(gòu)構(gòu)造簡單的立管相類似，因此在總體響應(yīng)分析時(shí)可將臍帶纜等效為結(jié)構(gòu)組成單一的鋼制管道[13]。表1所示為某臍帶纜等效結(jié)構(gòu)參數(shù)。

圖3 深海臍帶纜結(jié)構(gòu)Fig.3 Components of deep-sea umbilical

針對(duì)表1的臍帶纜進(jìn)行懸鏈線布局和緩波形布局設(shè)計(jì)，表2給出了臍帶纜在兩種布局下的基本布局參數(shù)。

表2 懸鏈線布局和緩波形布局基本參數(shù)Table 2 Main parameters of the free hanging catenary shape and the lazy-wave shape

基于表1和表2中的參數(shù)，根據(jù)式(3)和式(5)可分別計(jì)算出臍帶纜懸鏈線布局和緩波形布局的空間位置坐標(biāo)。進(jìn)而，在有限元軟件OrcaFlex中建立臍帶纜兩種布局下的總體模型，如圖4和圖5所示。臍帶纜采用Line單元建立而緩波形中浮子段采用Line with Floats單元模擬。假設(shè)臍帶纜躺底段平鋪于海底，其與海底的作用采用線性假設(shè)理論，土壤的垂直剛度為100 kPa/m2，土壤與臍帶纜的切向摩擦系數(shù)為0.5。臍帶纜的頂端懸掛于浮體底部，底端固定在海底井口處，連接處均可自由轉(zhuǎn)動(dòng)。

圖4 懸鏈線布局有限元模型Fig.4 Finite element model of the free hanging catenary

圖5 緩波形布局有限元模型Fig.5 Finite element model of the lazy wave

3 懸鏈線和緩波形布局總體響應(yīng)對(duì)比

3.1 靜態(tài)響應(yīng)對(duì)比

利用建立的數(shù)值模型，分別計(jì)算臍帶纜在懸鏈線和緩波形兩種布局下的靜態(tài)結(jié)構(gòu)響應(yīng)特性，結(jié)果對(duì)比如圖6和圖7所示。

圖6 靜態(tài)有效張力分布對(duì)比Fig.6 Comparison of static effective tension

圖6為臍帶纜兩種布局下有效張力對(duì)比圖。由圖可見，懸鏈線布局下有效張力沿著臍帶纜長度方向不斷減?。欢鴮?duì)于緩波形布局，由于浮子段浮力作用，有效張力先減小、再增大最后又減小，導(dǎo)致有效張力在浮子段末尾點(diǎn)出現(xiàn)一個(gè)極大值。同時(shí)，兩種布局下臍帶纜有效張力的最大值均發(fā)生在頂部懸掛點(diǎn)，且懸鏈線布局有效張力最大值明顯大于緩波形布局有效張力的最大值。因此，當(dāng)水深增加到一定深度后，臍帶纜采用懸鏈線布局滿足不了張力設(shè)計(jì)要求時(shí)，可以選用緩波形布局。

圖7 靜態(tài)彎曲曲率分布對(duì)比Fig.7 Comparison of static bending curvature

圖7為臍帶纜在懸鏈線布局和緩波形布局下的彎曲曲率響應(yīng)。由圖可知，懸鏈線布局下臍帶纜彎曲曲率最大值發(fā)生在觸地點(diǎn)處，而緩波形布局下彎曲曲率的最大值發(fā)生在浮子段上升段的末端，且緩波形布局形式下的最大彎曲曲率大于懸鏈線布局形式的。這是由于浮子段處產(chǎn)生拱形線型，造成了較大的局部彎曲。因此，在臍帶纜緩波布局設(shè)計(jì)時(shí)，應(yīng)該注意浮力塊參數(shù)的選取，避免產(chǎn)生較大的局部彎曲。

3.2 動(dòng)態(tài)響應(yīng)對(duì)比

分別計(jì)算臍帶纜兩種布局在表3中環(huán)境條件下的動(dòng)態(tài)響應(yīng)，對(duì)比結(jié)果如圖8～11所示。

表3 環(huán)境參數(shù)Table 3 Environment parameters

圖8為懸鏈線臍帶纜動(dòng)態(tài)有效張力分布圖。從圖中可以看出，頂部浮體的運(yùn)動(dòng)以及波浪載荷對(duì)整體臍帶纜的有效張力分布均有影響，其中對(duì)臍帶纜頂端張力的影響最大，越往下越小。而如圖9所示，外部動(dòng)態(tài)擾動(dòng)對(duì)緩波形臍帶纜懸垂段部分的有效張力有較大的影響，對(duì)浮子段及其以下部分則幾乎沒有影響。

圖8 懸鏈線臍帶纜動(dòng)態(tài)張力分布圖Fig.8 Dynamic effective tension of the free hanging catenary

圖9 緩波形臍帶纜動(dòng)態(tài)張力分布圖Fig.9 Dynamic effective tension of the lazy wave

圖10為懸鏈線臍帶纜動(dòng)態(tài)彎曲曲率分布圖。外部擾動(dòng)對(duì)臍帶纜上彎曲曲率的影響主要集中在頂部懸掛點(diǎn)區(qū)域以及觸地點(diǎn)區(qū)域，而對(duì)懸垂段和躺地段處的彎曲曲率影響很小。彎曲曲率的極值出現(xiàn)在頂部懸掛點(diǎn)和觸底點(diǎn)上，而對(duì)于頂部懸掛點(diǎn)出現(xiàn)的彎曲曲率極值，可以通過在臍帶纜的懸掛點(diǎn)處安裝抗彎裝置(Bend stiffer)來解決。

圖10 懸鏈線臍帶纜動(dòng)態(tài)曲率分布圖Fig.10 Dynamic bending curvature of the free hanging catenary

相對(duì)于懸鏈線布局，緩波形布局動(dòng)態(tài)響應(yīng)主要集中在頂部懸掛點(diǎn)區(qū)域，其余部分包括觸地點(diǎn)處的彎曲曲率所受影響很小，如圖11所示。對(duì)于頂部懸掛點(diǎn)出現(xiàn)的彎曲曲率極值，可以通過在懸掛點(diǎn)處安裝抗彎裝置來緩解。

圖11 緩波形臍帶纜動(dòng)態(tài)曲率分布圖Fig.11 Dynamic bending curvature of the lazy wave

4 緩波形布局參數(shù)敏感性分析

緩波形布局形式較為復(fù)雜，影響其線型分布的因素很多。下面對(duì)緩波形布局臍帶纜總體響應(yīng)進(jìn)行敏感性分析，分析浮子段長度、浮力因子和浮子段起始位置對(duì)臍帶纜總體響應(yīng)的影響。

4.1 浮子段長度敏感性分析

分別計(jì)算浮子段長度為200 m、300 m、400 m和500 m時(shí)臍帶纜的總體響應(yīng)。對(duì)應(yīng)的有效張力和彎曲曲率分布如圖12和圖13所示。

圖12 不同浮子斷長度下彎曲曲率分布對(duì)比Fig.12 Comparision of bending curvature for different lengths of buoyancy catenary

圖13 不同浮子斷長度下有效張力分布對(duì)比Fig.13 Comparison of effective tension for different lengths of buoyancy catenary

如圖12所示，浮子段越長，臍帶纜上最大的彎曲曲率以及觸地點(diǎn)處的彎曲曲率也越大。對(duì)于有效張力而言，如圖13所示，隨著浮子段的增長，頂部懸掛點(diǎn)處的有效張力減小，而浮子段末尾處以及下降段上的有效張力增大。這是因?yàn)橛捎诟×Φ脑龃螅瑴p小了臍帶纜對(duì)懸掛點(diǎn)處張力的需求，然而由于浮力增大，使得立管的下降段增長，從而使得浮子段末尾和下降段上的有效張力變大。

4.2 浮力因子敏感性分析

浮力因子的大小表征了單位長度浮力塊提供的浮力大小。為了研究浮力因子對(duì)臍帶纜整體線形性能的影響，保持其他參數(shù)不變，計(jì)算當(dāng)浮力因子為2、3、4時(shí)臍帶纜的靜態(tài)響應(yīng)，得到臍帶纜的彎曲曲率和有效張力分布如圖14和圖15所示。

圖14 不同浮力因子下彎曲曲率分布對(duì)比Fig.14 Comparison of bending curvature for different float factors

圖15 不同浮力因子下有效張力分布對(duì)比Fig.15 Comparison of effective tension for different float factors

如圖14所示，浮力因子越大，臍帶纜最大彎曲曲率以及觸地點(diǎn)處的彎曲曲率也越大。對(duì)于有效張力而言，如圖15所示，浮力因子的大小對(duì)緩波形臍帶纜有效張力的分布有很大的影響。當(dāng)浮力因子增大時(shí)，頂部懸掛點(diǎn)處的有效張力減小，而浮子段末端的有效張力增大。這符合實(shí)際情況，因?yàn)楦×σ蜃釉酱螅∽佣翁峁┑母×σ苍酱?，從而減小對(duì)頂部懸掛點(diǎn)處張力的需求。同時(shí)，浮力因子越大，緩波形臍帶纜下降段的長度也越大，使得下降段對(duì)浮子段末尾處的拉力增大，從而使得浮子段末尾處的有效張力變大。

4.3 浮子段起始位置敏感性分析

浮子段起點(diǎn)距頂部懸掛點(diǎn)的距離，即浮子段的起始位置，對(duì)緩波形臍帶纜的初始形狀和靜態(tài)響應(yīng)有很大的影響。為了研究浮子段起始位置的參數(shù)敏感性，在保持其他參數(shù)不變的情況下，計(jì)算浮子段起始位置分別為700 m、800 m、900 m、1 000 m和1 500 m時(shí)緩波形臍帶纜的總體響應(yīng)，得到有效張力以及彎曲曲率，如圖16和圖17所示。

圖16 不同浮子段起始位置下彎曲曲率分布對(duì)比Fig.16 Comparison of bending curvature for different starting locations of the buoyancy catenary

圖17 不同浮子段起始位置下有效張力分布對(duì)比Fig.17 Comparison of effective tension for different starting locations of the buoyancy catenary

圖16所示為不同浮子段起始位置下臍帶纜的初始線型。由圖可以看出，浮子段起始位置值越大，臍帶纜頂部懸掛點(diǎn)距底部觸地點(diǎn)之間的跨距越大，臍帶纜的波形越緩。同時(shí)，隨著浮子段起始位置值的增大，最大彎曲曲率逐漸減小。對(duì)于有效張力，由圖17可知，隨著浮子段起始位置值變大，懸掛段長度變大，從而導(dǎo)致懸掛點(diǎn)處的有效張力變大，同時(shí)由于浮子段起始位置值變大，下降段的長度變小，從而導(dǎo)致浮子段末端處的有效張力減小。

5 結(jié) 語

本文基于數(shù)值軟件OrcaFlex計(jì)算并對(duì)比了臍帶纜懸鏈線和緩波形兩種布局形態(tài)下的總體響應(yīng)結(jié)果，討論了浮子段長度、浮力因子和浮力塊起始位置對(duì)緩波形臍帶纜總體響應(yīng)的影響。得到如下結(jié)論：

(1) 相對(duì)于懸鏈線布局形態(tài)，緩波形布局下臍帶纜線型更為復(fù)雜，然而由于浮子段的作用，可以有效降低懸掛點(diǎn)處的有效張力，故而緩波形適合應(yīng)用于作業(yè)水深較深、懸掛點(diǎn)處張力較大情況，且浮力塊可以有效降低上部動(dòng)態(tài)響應(yīng)向底部觸地點(diǎn)處的傳播。

(2) 浮子段長度和浮力因子對(duì)緩波形臍帶纜總體線形的影響是一樣的，增大浮子段的長度或增大浮力因子都會(huì)使臍帶纜的波形變陡，最大彎曲曲率變大，但會(huì)減小臍帶纜的最大有效張力，并使得懸垂段處的有效張力減小，浮子段末尾處及下降段上的有效張力增大；同時(shí)增大浮子段起始位置會(huì)使臍帶纜的波形變緩，減小臍帶纜的最大彎曲曲率，卻增大其最大有效張力。

[1] 黃維平,白興蘭,李華軍.國外深水鋼懸鏈線立管研究發(fā)展現(xiàn)狀[J].中國海洋大學(xué)學(xué)報(bào),2009,39(2):290.

[2] 謝彬,段夢蘭,秦太驗(yàn),等.海洋深水立管的疲勞斷裂與可靠性評(píng)估研究進(jìn)展[J].石油學(xué)報(bào),2004,25(3):95.

[3] 宋磊建，付世曉，陳希恰，等.深海臍帶纜總體響應(yīng)特性比較研究[J].振動(dòng)與沖擊，2014,33(1):119.

[4] da Silveira L M Y, Tanaka R L, Novaes J P Z. Rayleigh damping effects on global analysis of umbilicals[C]. OMAE, 2011: 49584.

[5] Qiu W, Cao Q, Librino F, et al. Global design and analysis of umbilical in offshore applications[C]. OMAE, 2009: 79841.

[6] Howells H. Advances in steel catenary riser design[C]. DEEPTEC′95, 1995.

[7] Hatton S. Update on the design of steel catenary riser systems[R]. 1999.

[8] Yue B, Walters D, Yu W. Lazy wave SCR on turret moored FPSO[C]. ISOPE, 2011: 16.

[9] 孫麗萍，周佳，王佳琦.深水柔性立管的緩波形布置及參數(shù)敏感性分析[J].中國海洋平臺(tái)，2011，26(3)：37.

[10] Yang H, Li H. Sensitivity analysis of fatigue life prediction for deepwater steel lazy wave catenary risers [J]. Science China Technological Sciences, 2011, 54(7): 1881.

[11] Li S, Nguyen C. Dynamic response of deepwater lazy-wave catenary riser[C]. DOT International, 2010.

[12] Orcina. OrcaFlex manual: version 9.0a[M]. Cumbria: Orcina, 2006.

[13] International Organization for Standardization. ISO 13628-5. Petroleum and natural gas industries—design and operation of subsea production systems. Part 5: subsea umbilicals[S]. 2002.

GlobalResponseAnalysisfortheDeep-SeaUmbilical

CHEN Xi-qia1, SONG Lei-jian2, DU Xia-ying1, WAN Jun1, FU Shi-xiao2

(1.BohaiOilfieldResearchInstitute,TianjinBranch,CNOOCChinaLimited,Tianjin300452,China;2.StateKeyLaboratoryofOceanEngineering,ShanghaiJiaoTongUniversity,Shanghai200240,China)

We present the global analysis for the deep-sea umbilical. The design and analysis of the umbilical has become the research focus in the offshore engineering field. Global response analysis is the key issue in the research of umbilical. We present the theoretical models of the free hanging catenary shape and the lazy-wave shape at first. Based on the theoretical models, the umbilical models of the free hanging catenary shape and lazy-wave shape are then developed in OrcaFlex. The distribution characteristics of effective tension and bending curvature are calculated and compared. The results show that the shape of the lazy-wave is more complex than that of the free hanging catenary. However, the lazy-wave shape can reduce the effective tension at the top hanging points and dynamic response at the touchdown point, which makes it more applicable in the deep sea. Meanwhile, the sensitivity analysis of the generalized configuration parameters, including the length of the buoyancy catenary, the buoyancy factor and the starting location of the buoyancy catenary, is performed and the effect of the configuration parameters on the static characteristics of the umbilical is studied. The research results may provide a reference for the configuration design of the deep-sea umbilical.

deep-sea umbilical; global response; free hanging catenary; lazy-wave shape; effective tension; bending curvature

TE52

A

2095-7297(2015)02-0111-07

2015-01-22

陳希恰(1988—)，男，碩士，主要從事海洋工程結(jié)構(gòu)物強(qiáng)度與疲勞分析方面的研究。