黃 熠 ，王爾鈞 ，張 崇，孟文波，同武軍，趙維青，趙蘇文，段明付

(1.中海石油(中國(guó))有限公司 湛江分公司，廣東 湛江524057；2.中海油能源發(fā)展股份有限公司工程技術(shù)深水鉆采技術(shù)公司，廣東 湛江524057)

自從海洋油氣開(kāi)采提升到國(guó)家戰(zhàn)略，各油田公司紛紛加大在海洋的油氣開(kāi)采力度。水下采油樹(shù)是水下生產(chǎn)系統(tǒng)進(jìn)行海洋油氣開(kāi)采的關(guān)鍵裝備，在石油開(kāi)采中扮演者不可替代的作用[1-3]，其具有質(zhì)量大、結(jié)構(gòu)復(fù)雜等特點(diǎn)，多變的海況給下放安裝帶來(lái)極大的挑戰(zhàn)[4]。

鋼絲繩直接吊裝和鉆桿下放安裝是水下采油樹(shù)的主要安裝方法，目前已有很多高校和研究人員對(duì)水下采油樹(shù)的安裝進(jìn)行了相關(guān)研究。周美珍等[5]基于平臺(tái)月池安裝，并對(duì)平臺(tái)偏移、采油樹(shù)質(zhì)量、環(huán)境參數(shù)、鉆桿參數(shù)等因素展開(kāi)對(duì)水下采油樹(shù)不同安裝階段進(jìn)行分析研究；歐陽(yáng)賽賽等[6]建立了三維流固耦合計(jì)算模型，研究了淺水采油樹(shù)鉆桿下放過(guò)程的應(yīng)力情況和橫向偏移響應(yīng)；脫浩虎等[7]基于水動(dòng)力學(xué)軟件OrcaFlex建立了水下采油樹(shù)鉆桿下放的分析模型，研究了環(huán)境因素對(duì)下放過(guò)程動(dòng)態(tài)響應(yīng)的影響；張偉國(guó)等[8]基于鋼絲纜繩和鉆桿2種下放方案的特點(diǎn)，研究了這2種下放安裝方法下入500 m水深，水下采油樹(shù)運(yùn)動(dòng)響應(yīng)情況與鉆桿受力分布情況，同時(shí)研究了環(huán)境因素對(duì)2種下放安裝方法的敏感度；湯建冬等[9]基于材料力學(xué)小變形梁理論構(gòu)建采油樹(shù)鉆桿下放模型，采用有限差分法數(shù)值求解，并與OrcaFlex軟件計(jì)算結(jié)果展開(kāi)對(duì)比驗(yàn)證研究；Thomas等[10]詳細(xì)分析了水下采油樹(shù)的下放安裝工藝；R.VOSS等[11]詳細(xì)闡述了水下采油樹(shù)的安裝工藝，并記錄下放過(guò)程的相關(guān)技術(shù)問(wèn)題。王瑩瑩等[12]利用OrcaFlex水動(dòng)力學(xué)軟件研究了不同環(huán)境下1 500 m水深水下管匯安裝工藝過(guò)程中鉆桿的受力情況，計(jì)算得到鉆桿頂端部分屬于危險(xiǎn)部分；顧純巍等[13]研究了不同影響因素對(duì)水下應(yīng)急封井裝置下入過(guò)程中鋼絲纜繩最大受力和裝置運(yùn)動(dòng)響應(yīng)的影響。

由此可見(jiàn)，水下采油樹(shù)的鉆桿下放受力極為復(fù)雜，為了解決水下采油樹(shù)鉆桿下放過(guò)程中出現(xiàn)的技術(shù)問(wèn)題，以及下放安裝過(guò)程中鉆桿的相關(guān)受力問(wèn)題，筆者基于OrcaFlex水動(dòng)力學(xué)軟件進(jìn)行數(shù)值仿真，針對(duì)東方1-1氣田海域水深73 m，開(kāi)展水下采油樹(shù)下放過(guò)程影響因素敏感性分析。研究結(jié)果將對(duì)水下采油樹(shù)的下放安裝起到一定的參考作用。

1 鉆桿下放受力模型分析

水下采油樹(shù)鉆桿下放過(guò)程主要分為5個(gè)階段。第1階段為裝船固定階段，第2個(gè)階段為提升操作階段，第3個(gè)階段為穿越飛濺區(qū)階段，第4個(gè)階段為深水下放階段，第5個(gè)為與井口對(duì)接階段。不同階段所受載荷和邊界條件不同，本文主要對(duì)水下采油樹(shù)鉆桿下放過(guò)程環(huán)境影響展開(kāi)研究，因此選擇第4階段研究。目前對(duì)于鉆桿下放海洋裝備動(dòng)力學(xué)研究主要利用材料力學(xué)小變形梁理論和有限元分析方法[14]。

建模時(shí)做如下假設(shè)：

1) 鉆桿的材料為均勻各向同性的。

2) 環(huán)境參數(shù)的運(yùn)動(dòng)與平臺(tái)及鉆桿的運(yùn)動(dòng)在同1個(gè)平面內(nèi)。

3) 鉆桿受外力變形屬于小變形。

水下采油樹(shù)鉆桿下放過(guò)程受力分析模型如圖1所示。

圖1 水下采油樹(shù)鉆桿下放過(guò)程受力分析模型

由鉆桿的受力分析模型可知，鉆桿主要受到海流力、波浪力、重力、浮力、張緊力和彎矩等載荷，將受力變形后的鉆桿劃分為n個(gè)微分單元，由挪威船級(jí)社DNV RP-C205規(guī)范[15]可計(jì)算鉆桿第i個(gè)微分單元所受水動(dòng)力載荷Fi和采油樹(shù)所受水動(dòng)力Fc。

(1)

(2)

式中：i=1，3，…，n；ρ為海水密度，kg/m3；Cd為鉆桿水平拖曳力系數(shù)；CD為采油樹(shù)水平拖曳力系數(shù)；CM1為鉆桿慣性力系數(shù)；CM2為采油樹(shù)性力系數(shù)；d為鉆桿外徑，m；vi為第i個(gè)微分單元處的海流流速，m/s；μ為波浪引起的水質(zhì)點(diǎn)速度，m/s；Ap為采油樹(shù)水平拖曳力面積，m2；H為纜繩微分單元的高度，m；V為水下采油樹(shù)的排水體積，m3。

鉆桿微分單元的自重力Gi可表示為：

Gi=ωg·dx

(3)

式中：ω為線密度，kg/m。

微分段的力學(xué)微分方程可表示為：

(4)

Ti=Fxicosθ+Fyisinθ

(5)

Fxi=(Gc-Ff)+ω(n-i)dx·g

(6)

(7)

式中：EI為鉆桿的抗彎剛度，N/m2；Ti為軸向張力，N；θ為鉆桿微元與豎直方向的夾角，rad；Gc采油樹(shù)為采油樹(shù)產(chǎn)生的重力，N；Ff為水下采油樹(shù)的浮力，N。

2 水下采油樹(shù)鉆桿下放仿真分析

水下采油樹(shù)主要由采油樹(shù)下放工具、本體、閥組和框架等組成?；诎霛撌狡脚_(tái)模型參數(shù)，利用OrcaFlex軟件進(jìn)行繪制，該平臺(tái)的水動(dòng)力學(xué)參數(shù)主要由RAO算子決定，水下采油樹(shù)的本體和下放工具利用6D buoy浮體建模，框架和相關(guān)閥組利用Shape中的Elastic solid進(jìn)行建模，再將不同部分連接。水下采油樹(shù)下放安裝耦合系統(tǒng)模型如圖2所示。

圖2 水下采油樹(shù)下放安裝耦合系統(tǒng)模型

2.1 環(huán)境參數(shù)的設(shè)定

本文主要研究東方1-1氣田的風(fēng)、浪、流等環(huán)境參數(shù)對(duì)水下采油樹(shù)鉆桿下放安裝過(guò)程中的受力和運(yùn)動(dòng)響應(yīng)問(wèn)題，所以空氣、海面、海床等選取默認(rèn)。

2.1.1 風(fēng)參數(shù)

對(duì)于水下采油樹(shù)所受環(huán)境載荷，風(fēng)載荷的比例占15%。在黏性阻力系數(shù)和迎風(fēng)面積一定的情況下，風(fēng)速是影響風(fēng)載荷的關(guān)鍵因素。東方1-1氣田風(fēng)參數(shù)如表1所示。

表1 東方1-1氣田風(fēng)參數(shù)

2.1.2 浪參數(shù)

表2 東方1-1氣田波浪參數(shù)

2.1.3 海流參數(shù)

由于流速數(shù)據(jù)很少，因此OrcaFlex軟件通過(guò)在Environment中設(shè)置表層流速和海床流速，使用軟件自帶的Power Law方法形成0～73 m連續(xù)不同深度流速曲線圖，具體海流參數(shù)如表3所示。

表3 東方1-1氣田海流參數(shù)

2.2 鉆桿建模

OrcaFlex中使用Line模型可以模擬鉆桿及繩索，它提供3種類(lèi)型，分別為Homogeneous pipe、Equivalent line和General，其中Homogeneous pipe用于模擬均質(zhì)管線，Equsivalent line主要用于模擬多層復(fù)合管結(jié)構(gòu)，General則用于除這2種之外的其他類(lèi)型管線，因此選取General作為鉆桿的類(lèi)型。鉆桿參數(shù)設(shè)置如表4所示。

表4 鉆桿參數(shù)

3 基于OrcaFlex的水下采油樹(shù)鉆桿下放安裝敏感性參數(shù)仿真研究

3.1 風(fēng)速對(duì)鉆桿受力和水下采油樹(shù)運(yùn)動(dòng)響應(yīng)的影響

由表1知，重現(xiàn)期1 a時(shí)平均風(fēng)速為28.4 m/s， 重現(xiàn)期10 a年時(shí)平均風(fēng)速為41.0m/s，重現(xiàn)期50 a平均風(fēng)速為50.3 m/s，重現(xiàn)期100 a平均風(fēng)速為54.2 m/s。風(fēng)主要是直接作用于鉆井平臺(tái)和引起海水的波動(dòng)，對(duì)采油樹(shù)的下放造成不可忽視的安全問(wèn)題，因此展開(kāi)了對(duì)風(fēng)速這個(gè)單一因素的研究。不同風(fēng)速影響下的鉆桿受力狀態(tài)和采油樹(shù)運(yùn)動(dòng)響應(yīng)計(jì)算結(jié)果如表5所示。

表5 不同風(fēng)速的影響數(shù)據(jù)

鉆桿所受到應(yīng)力和彎矩極值同樣是分布在鉆桿的頂端位置，不同風(fēng)速對(duì)鉆桿的Von Mises應(yīng)力以及彎矩(包括x、y方向彎矩)影響總體較小，其中鉆桿的上端點(diǎn)的應(yīng)力以及彎矩隨風(fēng)速增大而增大，鉆桿的下端點(diǎn)的Von Mises應(yīng)力受風(fēng)速影響較小，彎矩值為0。

不同風(fēng)速對(duì)水下采油樹(shù)在x、y方向的偏移量和傾斜角度影響較小，主要是因?yàn)橹苯幼饔迷谒虏捎蜆?shù)上的載荷只有海流載荷，風(fēng)載荷只能作用于鉆井平臺(tái)，通過(guò)與鉆桿之間的耦合作用，間接影響水下采油樹(shù)。采油樹(shù)在下放過(guò)程中最大的偏斜角度在9°左右。

3.2 海流流速對(duì)鉆桿受力和水下采油樹(shù)運(yùn)動(dòng)響應(yīng)的影響

海流是制約水下結(jié)構(gòu)物下放的重要因素，其大小變化莫測(cè)，但一般海面流速都在0.50～2.04 m/s，2.04 m/s已是百年一遇，故不再對(duì)更大的流速進(jìn)行討論。本文選擇了1.09 m/s(一年一遇)、1.46 m/s(十年一遇)、1.86 m/s(五十年一遇)和2.04 m/s(百年一遇)4種流速。海流主要是作用在鉆桿和采油樹(shù)上，這會(huì)給鉆桿造成一個(gè)很大的彎矩，影響到鉆桿的正常使用。通常，鉆桿的最大彎矩出現(xiàn)在鉆臺(tái)處，所以研究海流對(duì)鉆桿受力的影響尤為必要，其他環(huán)境條件，浪高選擇1 a一遇的情況，浪高6.3 m，周期為8 s；風(fēng)速選擇1 a一遇的平均風(fēng)速28.4 m/s。不同海流流速影響下的鉆桿受力狀態(tài)和采油樹(shù)運(yùn)動(dòng)響應(yīng)計(jì)算結(jié)果如圖3～8所示。

大數(shù)據(jù)時(shí)代給管理會(huì)計(jì)工作帶來(lái)的不止是機(jī)遇也有更大的挑戰(zhàn)，新時(shí)代的特點(diǎn)就是機(jī)會(huì)與挑戰(zhàn)并存。所以，企業(yè)一定要對(duì)此加強(qiáng)重視，要緊跟時(shí)代潮流，牢牢抓住機(jī)遇，整合分析自身的問(wèn)題，并且結(jié)合實(shí)際情況對(duì)癥下藥，充分利用大數(shù)據(jù)的優(yōu)勢(shì)為企業(yè)的發(fā)展謀求更多機(jī)遇。

圖3 不同流速下鉆桿的Von Mises應(yīng)力

圖4 不同流速下鉆桿的總彎矩

圖5 不同流速下鉆桿的不同方向彎矩

從圖3～8可以看出，在表層流速?gòu)?.09 m/s到2.04 m/s增大的過(guò)程中，鉆桿頂端受到的彎矩大小隨著海流的增大而增大，說(shuō)明不同流速對(duì)鉆桿受到的有效張力影響較大。同時(shí)，采油樹(shù)的偏移也隨著海流流速的增加而不斷增加。

隨著海流流速的增大，對(duì)鉆桿的Von Mises應(yīng)力以及彎矩(包括x、y方向彎矩)影響較大，同時(shí)也能看出鉆桿受到的最大Von Mises應(yīng)力和彎矩均位于鉆桿頂端處，百年一遇海流流速時(shí)鉆桿頂端最大Von Mises應(yīng)力為818.38 MPa (小于鉆桿(S135)的屈服極限930 MP)，最大彎矩為315.68 kN·m，x方向最大彎矩為223.07 kN·m和y方向最大彎矩為223.37 kN·m。

圖6 不同流速下的偏移曲線

圖7 不同流速下采油樹(shù)在x、y方向上最大偏移量

圖8 不同流速下采油樹(shù)的傾斜角度

在采油樹(shù)經(jīng)過(guò)飛濺區(qū)時(shí)，隨著海流流速的增大，采油樹(shù)的x、y方向偏移量以及傾斜角也增大，且影響較大。

在百年一遇的海況下采油樹(shù)在飛濺區(qū)的x、y方向偏移量以及傾斜角有最大值，分別為9.7 m 、13.2 m和16.4°。

采油樹(shù)在接近井口時(shí)，采油樹(shù)的傾斜角變化較小，x、y方向偏移量隨著鉆桿下放長(zhǎng)度的增加而緩慢增加。

3.3 浪高對(duì)鉆桿受力和水下采油樹(shù)運(yùn)動(dòng)響應(yīng)的影響

浪高是采油樹(shù)下放過(guò)程不可忽視的因素，海浪的波動(dòng)會(huì)引起鉆井平臺(tái)的上下浮動(dòng)以及左右擺動(dòng)，這就要求平臺(tái)具有一定的升沉補(bǔ)償能力，會(huì)形成1個(gè)飛濺區(qū)域，給采油樹(shù)的底部形成較大的沖擊作用力，同時(shí)會(huì)間接地影響鉆桿的受力情況，也會(huì)引起采油樹(shù)更大的傾斜。不同浪高影響下的鉆桿受力狀態(tài)和采油樹(shù)運(yùn)動(dòng)響應(yīng)計(jì)算結(jié)果如圖9～13所示。

從圖9～10可以看出，隨著波浪高度的增大，鉆桿頂端的Von Mises應(yīng)力和彎矩(包括x、y方向彎矩)均增加，說(shuō)明浪高對(duì)鉆桿Von Mises應(yīng)力和彎矩影響較大；當(dāng)浪高為百年一遇(13.4 m)時(shí)，鉆桿上下端點(diǎn)的最大Von Mises應(yīng)力分別為482.9 MPa、56.6 MPa;鉆桿的最大彎矩為175.5 kN·m，其中x方向彎矩為129.5 kN·m，y方向彎矩為134.8 kN·m。

圖9 不同浪高下鉆桿的Von Mises應(yīng)力

圖10 不同浪高下鉆桿上端點(diǎn)的彎矩

從圖11～13可以看出，隨著波浪高度的增大，在采油樹(shù)經(jīng)過(guò)飛濺區(qū)時(shí)，采油樹(shù)的x、y方向最大偏移量以及最大傾斜角也增大，且影響較大。

圖11 不同浪高下采油樹(shù)的偏移量

圖12 不同浪高下采油樹(shù)的偏移極值

圖13 不同浪高下采油樹(shù)的傾斜角

在百年一遇的浪高情況下，采油樹(shù)的x方向偏移為8.9 m，y方向偏移為9.9 m，采油樹(shù)的傾斜角為15.2°。

在接近井口時(shí)，采油樹(shù)的傾斜角以及偏移量變化較小，偏移量(x、y方向)隨著鉆桿下放長(zhǎng)度的增加而緩慢增加。但采油樹(shù)在水中的波動(dòng)幅度隨著浪高的增加而增大。

3.4 方向角對(duì)鉆桿受力和水下采油樹(shù)運(yùn)動(dòng)響應(yīng)影響

鉆井平臺(tái)由于結(jié)構(gòu)等的原因，在不同的海況方向會(huì)產(chǎn)生不同的運(yùn)動(dòng)，從而間接的影響水下采油樹(shù)的下放，不同海況方向影響下的鉆桿受力狀態(tài)和采油樹(shù)運(yùn)動(dòng)響應(yīng)如圖14～17所示。

從圖14～15可以看出，在不同的海況方向激勵(lì)下，鉆桿上端點(diǎn)的Von Mises應(yīng)力受海況方向的影響較大。當(dāng)海況方向?yàn)?°到180 °方向上時(shí)影響相對(duì)較小，45°到135°方向上時(shí)影響相對(duì)較大，其中0°方向有最小值258.8 MPa，135°有最大值552.9 MPa；海況方向?qū)︺@桿的下端點(diǎn)的影響較小，在45°有最大值80.5 MPa。鉆桿的上下端點(diǎn)的彎矩變化趨勢(shì)相同，其總彎矩值變化不大，但依然在45°到135°方向有較大值，在0°到180°方向有較小值，可能與采油樹(shù)的迎流面積相關(guān)。

圖15 不同海況方向下鉆桿的上端點(diǎn)彎矩極值

從圖16～17可以看出，在不同的海況方向激勵(lì)下，采油樹(shù)經(jīng)過(guò)飛濺區(qū)時(shí)，海況方向?qū)A斜角的影響較大，海況方向?yàn)?80°時(shí)采油樹(shù)有最大傾斜角度值14.0°；在接近井口時(shí)，海況的方向?qū)Σ捎蜆?shù)的傾斜角影響較小。

圖16 不同海況方向下采油樹(shù)的傾斜角

圖17 不同海況方向下采油樹(shù)的偏移

采油樹(shù)的偏移量受海況方向的影響較大，x方向在0°海況方向有最大偏移值5.91 m，y方向在90°海況方向有最大偏移值5.93 m。

3.5 水下采油樹(shù)敏感性參數(shù)分析

水下采油樹(shù)下放安裝耦合系統(tǒng)，在不同環(huán)境因素的影響下，展現(xiàn)出不同運(yùn)動(dòng)規(guī)律和特性[16-17]。水下采油樹(shù)鉆桿下放過(guò)程影響因素眾多，其中多變的海況是影響其安全下放的一個(gè)重要因素，為此對(duì)構(gòu)成海況的單一因素，風(fēng)速、海流流速、浪高和方向角進(jìn)行逐一研究，基于對(duì)水下采油樹(shù)下放過(guò)程穩(wěn)定性的影響程度不同，對(duì)參數(shù)敏感性進(jìn)行打分評(píng)價(jià)。水下采油樹(shù)安裝穩(wěn)定性的參數(shù)敏感性分析如表5所示。浪高對(duì)采油樹(shù)安裝穩(wěn)定性的影響最大，其次是流速，方向角的影響較小，風(fēng)速對(duì)采油樹(shù)安裝穩(wěn)定性的影響最小，其影響大小排序?yàn)槔烁?海流流速>方向角>風(fēng)速。

表5 水下采油樹(shù)敏感性參數(shù)分析

4 結(jié)論

1) 基于海工軟件OrcaFlex建立平臺(tái)-鉆桿-水下采油樹(shù)多體下放系統(tǒng)仿真模型，同時(shí)對(duì)水下采油樹(shù)采用鉆桿下放進(jìn)行模擬仿真。研究了水下采油樹(shù)在下放過(guò)程中橫向位移以及傾斜角度的變化情況，同時(shí)也采用了鋼絲繩下放中一樣的控制變量法，模擬分析了風(fēng)速、流速、浪高以及海況方向等不同因素對(duì)水下采油樹(shù)下放安裝過(guò)程中的影響。

2) 鉆桿所受到的Von Mises應(yīng)力以及彎矩都主要集中在鉆桿上部靠近鉆桿頂端處，在鉆桿與采油樹(shù)連接處存在較小應(yīng)力，沒(méi)有彎矩集中。應(yīng)當(dāng)避免在過(guò)大海浪、海流等海況下進(jìn)行采油樹(shù)下放作業(yè)，否則會(huì)造成鉆桿的彎曲變形。

3) 在敏感性參數(shù)分析方面，風(fēng)速對(duì)下放過(guò)程中鉆桿受力、彎矩以及水下采油樹(shù)的運(yùn)動(dòng)響應(yīng)影響很?。缓Ａ髁魉?、浪高、海況的方向?qū)ζ溆绊戄^大，其影響大小排序?yàn)槔烁?海流流速>方向角>風(fēng)速，其中海流流速和浪高都是正影響。