祖 巍 馬冬輝 李彥麗 楊肖龍

（中海油能源發(fā)展裝備技術(shù)有限公司工程設(shè)計研發(fā)中心 天津300452）

引 言

隨著海洋石油工業(yè)向越來越深的水域的發(fā)展，傳統(tǒng)的淺海固定式平臺已經(jīng)不能適應(yīng)深水油田的開發(fā)。因此，發(fā)展適應(yīng)深水石油開發(fā)的移動式平臺尤為重要。其中，張力腿平臺（TLP）是深水順應(yīng)式平臺的一種典型形式。張力腿平臺在20世紀80年代初獲得蓬勃發(fā)展，可看作是帶有張力系泊系統(tǒng)的半潛式平臺。張力腿平臺不僅垂蕩運動較?。?］，而且控制方向的張力對非控制方向的運動具有牽制作用，所以漂移搖擺的幅度也比一般半潛式平臺小，具有波浪中運動性能好、抵抗惡劣環(huán)境能力強的優(yōu)點，成為海洋工業(yè)深水采油平臺的主要形式。［2］由于張力腿平臺的運動對結(jié)構(gòu)設(shè)計帶來的新問題和新挑戰(zhàn)是研究的主要難點之一，目前我國也在積極自主研發(fā)深水張力腿平臺。

對于淺水固定式平臺上的鉆井設(shè)備模塊，大多采用靜力分析，對于張力腿平臺上的鉆井設(shè)備模塊，必須考慮平臺的運動對其產(chǎn)生的影響。本文使用有限元分析軟件SEASAM分析張力腿平臺運動響應(yīng)，得到各個波浪頻率和不同浪向下的運動響應(yīng)幅值（RAO）、極值加速度等水動力參數(shù)；并在此基礎(chǔ)上，使用有限元分析軟件ANSYS對TLP平臺的上部鉆井設(shè)備模塊進行結(jié)構(gòu)準靜力分析。

1 張力腿平臺極值加速度計算

張力腿平臺水動力主要采用DNV開發(fā)的計算軟件SESAM來進行計算，其中以使用GeniE模塊和HydroD模塊為主。GeniE模塊主要用于建立TLP的幾何模型以及進行濕表面網(wǎng)格的劃分。HydroD模塊主要將GeniE建立的有限元模型作為輸入，建立TLP水動力計算模型，定義TLP質(zhì)量模型，并采用Wadam模塊對船體以及海洋結(jié)構(gòu)物的波浪載荷、運動響應(yīng)進行分析計算，［3］采用Postresp模塊對水動力計算結(jié)果進行后處理；從而得到張力腿平臺在不同工況下的極值加速度，并作為下一步對張力腿平臺上的鉆井設(shè)備模塊進行準靜力分析的輸入。

1.1 張力腿平臺模型

TLP平臺的船體結(jié)構(gòu)建成濕表面模型，采用Panel單元模擬，根據(jù)單位勢流理論計算波浪力；張力腿筋鍵和立管采用Morison單元模擬［4］，根據(jù)莫里森方程計算波浪力?？紤]正常作業(yè)工況，不可預(yù)期作業(yè)工況和可預(yù)期極端工況三種不同的工況。各工況下平臺的船體濕表面模型、張力腿筋鍵的水動力模型如圖1所示。

圖1 張力腿平臺水動力模型

1.2 張力腿平臺的極值加速度

HydroD基于三維勢流理論［5］可計算得到平臺運動響應(yīng)幅值RAO。結(jié)合TLP平臺RAO計算結(jié)果，可計算平臺的極值加速度，這就需要計算平臺的加速度傳遞函數(shù)，且將六個自由度方向的加速度轉(zhuǎn)化為X、Y、Z三個自由度方向的線性加速度，轉(zhuǎn)化公式（僅以Z方向自由度為例）見式（1）：

根據(jù)上述轉(zhuǎn)化公式，可得平臺三個方向的加速度傳遞函數(shù)。

TLP平臺的極值加速度采用譜方法計算，根據(jù)計算得到的加速度傳遞函數(shù)， 輸入波浪譜，可計算各波浪入射角度的加速度響應(yīng)譜［6］。根據(jù)加速度響應(yīng)譜模擬3個小時的短期海況，取最大值，則為加速度極值。

各工況下TLP平臺重心位置極值加速度如表1所示。鑒于平臺及其上部鉆井設(shè)備模塊等結(jié)構(gòu)是共同發(fā)生運動的，極值加速度將作為下一步準靜力分析的輸入數(shù)據(jù)，并以慣性載荷的形式施加于上部結(jié)構(gòu)。

表1 TLP平臺極值加速度m / s2

2 鉆井設(shè)備模塊準靜力分析有限元模型

2.1 有限元模型

鉆井設(shè)備模塊包括井架、鉆井甲板和下底座。鉆井甲板有一層甲板和在甲板上縱橫分布的4根主梁。4根主梁相互交叉之處，是邊長9.144 m的正方形。正方形區(qū)域內(nèi)2根主梁支撐著轉(zhuǎn)盤，另2根主梁支撐轉(zhuǎn)盤大梁。正方形區(qū)域外的懸臂梁支撐著鉆井甲板上其他主要部位。井架結(jié)構(gòu)直接坐落在正方形的四角上。鉆井甲板的大部分區(qū)域都是加筋板組成。鉆井甲板上的主梁還直接支撐西側(cè)的貓道和東側(cè)的測井平臺。下底座由上至下有四層平面，依次為東西向滑道、中層甲板、下層甲板和BOP甲板。

在ANSYS中使用Beam188單元建立主梁和其他骨材［7］，其中T型鋼屈服強度為345 MPa，工字型鋼結(jié)構(gòu)屈服強度為355 MPa。使用shell 63單元建立所有的板結(jié)構(gòu)，其屈服強度為355 MPa。鉆井設(shè)備模塊上的所有設(shè)備采用MASS 21單元進行模擬。所有單元彈性模量E= 2.06e5 MPa，泊松比μ=0.3，密度為ρ=7.85 kg/m3。鉆井設(shè)備模塊邊界處連接主甲板強構(gòu)件，此處剛度較大，發(fā)生的位移與整體鉆井設(shè)備模塊相比可以忽略不計，因此邊界條件取作對4根主立柱施加鉸支約束。鉆井設(shè)備模塊有限元模型如圖2所示。

圖2 鉆井設(shè)備模塊有限元模型

2.2 結(jié)構(gòu)加載

鉆井設(shè)備模塊外載分為重力載荷和環(huán)境載荷兩大類。

重力載荷包括靜載荷和活載荷。靜載荷指結(jié)構(gòu)質(zhì)量和設(shè)備干重，其中主要設(shè)備包括絞車、鎖緊裝置、管道閥門、拖鏈等?；钶d荷指大鉤載荷、轉(zhuǎn)盤載荷、立根載荷、液體載荷和甲板可變載荷等。鉆井設(shè)備模塊重力載荷約重2 111 t。結(jié)構(gòu)鋼材自重是通過定義結(jié)構(gòu)幾何尺寸和鋼材密度由程序自行計算，定義結(jié)構(gòu)長、寬、高以及結(jié)構(gòu)構(gòu)件尺寸和材料密度。設(shè)備質(zhì)量模擬成“MASS21”單元［8］。設(shè)備質(zhì)量的模擬基于設(shè)備的長、寬、高和設(shè)備的重心位置，立根載荷的施加與設(shè)備質(zhì)量模擬方式類似，并按照布置圖施加到結(jié)構(gòu)上的不同位置。

環(huán)境載荷包括慣性載荷和風(fēng)載荷。

慣性載荷由平臺運動引起（如表1所示），慣性載荷以加速度的形式施加于鉆井設(shè)備模塊。

風(fēng)載荷參數(shù)見表2。

表2 風(fēng)載荷參數(shù)m/s

風(fēng)載荷的計算一般基于風(fēng)壓和受風(fēng)面積。為便于分析，文中僅考慮受風(fēng)面積。風(fēng)載荷根據(jù)API 4F 8.3章節(jié)中的相關(guān)規(guī)范進行計算。［9］按照45°的增幅計算0°～ 360°范圍內(nèi)8個方向的風(fēng)載荷，因為鉆井設(shè)備模塊是對稱結(jié)構(gòu)，故其他方向的風(fēng)載荷均可通過0°和90°方向風(fēng)載荷計算求得。計算求得鉆井設(shè)備模塊受到的風(fēng)載荷結(jié)果如表3所示。風(fēng)載荷以節(jié)點集中力的形式施加至模型受風(fēng)面的所有節(jié)點。

表3 風(fēng)載荷計算結(jié)果

3 鉆井設(shè)備模塊準靜力分析

TLP平臺鉆井設(shè)備模塊鉆井甲板以上部分可根據(jù)平臺作業(yè)時鉆探井位的不同而作相應(yīng)移動，因此準靜力分析共對三種不同的模型進行，這三種模型分別對應(yīng)東井位、中井位和西井位。對每一種井位，又分為五種工況、即正常作業(yè)工況（存在大鉤載荷）、正常作業(yè)工況（存在轉(zhuǎn)盤載荷）、不可預(yù)期作業(yè)工況、可預(yù)期極端工況（考慮液體載荷）和可預(yù)期極端工況（不考慮液體載荷）。不同工況的環(huán)境載荷及鉆井設(shè)備模塊上運行的設(shè)備有所區(qū)別。對每一種工況，還需考慮風(fēng)載荷和波浪載荷（已通過水動力分析轉(zhuǎn)化為慣性載荷）從不同方向作用于鉆井設(shè)備模塊上，共考慮0°、45°、90°、135°、180°、225°、270°、315°共八種不同的入射角度。

綜上所述，對鉆井設(shè)備模塊共進行120種情形下的準靜力分析。為詳細描述分析結(jié)果，將整體結(jié)構(gòu)分為鉆井甲板、中層甲板和下層甲板；分別提取120個工況下各層甲板板的最大等效應(yīng)力，將鉆井設(shè)備模塊分為鉆井甲板、中層甲板、下層甲板和下底座等部分，分別提取120個工況下各部分的最大位移值。僅以其中一個工況為例描述鉆井設(shè)備模塊應(yīng)力和變形。

3.1 鉆井設(shè)備模塊板結(jié)構(gòu)應(yīng)力分析

鉆井設(shè)備模塊等效應(yīng)力云圖參見下頁圖3。由圖3可知，下層甲板等效應(yīng)力最大，鉆井甲板次之，中層甲板最小。對比最大應(yīng)力的發(fā)生位置可知，鉆井甲板最大等效應(yīng)力發(fā)生在放置絞車的位置。絞車質(zhì)量較大，對鉆井甲板應(yīng)力分布影響較大。中層甲板等效應(yīng)力較大位置主要集中于與立柱和斜撐相連的區(qū)域，這部分結(jié)果主要是由于存在應(yīng)力集中現(xiàn)象所造成的；其次，發(fā)生在振動篩與BOP控制單元布置的區(qū)域。除此之外，中層甲板板單元等效應(yīng)力相對較小。下層甲板的最大應(yīng)力發(fā)生在板的邊緣位置與支撐骨材相連的區(qū)域，這也是由應(yīng)力集中造成的。

考慮到應(yīng)力計算結(jié)果較大位置主要為應(yīng)力集中處，而有限元模型中對這些結(jié)構(gòu)局部細節(jié)進行了簡化，因此結(jié)構(gòu)的實際應(yīng)力應(yīng)低于計算值，可以認為結(jié)構(gòu)實際強度儲備更大一些。中層甲板可考慮適度優(yōu)化，減少用鋼。

圖3 鉆井設(shè)備模塊甲板等效應(yīng)力云圖

3.2 鉆井設(shè)備模塊梁結(jié)構(gòu)應(yīng)力分析

在對鉆井設(shè)備模塊進行準靜力分析時，按照構(gòu)件截面的不同，將梁單元的所有構(gòu)件分為24組，并對其進行編號。分析24組構(gòu)件計算結(jié)果可知，編號為T5、T5C和CO2組骨材等效應(yīng)力較大，安全裕量相對較小。

T5組骨材位于鉆井甲板上，截面形式為T型，型號為T 250×580×20×10。梁截面較小，最大等效應(yīng)力為257.22 MPa（許用應(yīng)力為276 MPa），發(fā)生位置如圖4所示。該構(gòu)件長度約9 m，由于跨距大，因而其等效應(yīng)力也相對較大。

T5C組骨材位于中層甲板上，截面形式為T型，型號為T 200×494×16×10。梁截面較小，最大應(yīng)力為221.82 MPa（許用應(yīng)力為276 MPa），發(fā)生位置如圖4所示，該位置設(shè)有振動篩。

CO2組骨材位于中層甲板和下層甲板之間，起到支撐中層甲板的作用，其截面形式為圓管，尺寸為323.9 mm×12.5 mm，最大應(yīng)力為206.83 MPa（許用應(yīng)力為232 MPa），最大應(yīng)力發(fā)生位置見圖4。

圖4 鉆井設(shè)備模塊受力較大梁結(jié)構(gòu)

位于鉆井甲板上的T5骨材，其最大等效應(yīng)力多發(fā)生于正常作業(yè)工況，這是因為在正常作業(yè)工況下，鉆井甲板上立根、可變載荷、液體載荷等較大。位于中層甲板T5C和下底座CO2骨材的最大等效應(yīng)力多發(fā)生在可預(yù)期極端工況下，其原因是環(huán)境載荷對骨材的應(yīng)力產(chǎn)生較大影響。

對于這些局部受力較大的構(gòu)件，可以適當考慮使用截面更大的型材，或者增加構(gòu)件數(shù)量、減小部分構(gòu)件的跨距，使結(jié)構(gòu)更為安全。

3.3 鉆井設(shè)備模塊變形分析

鉆井設(shè)備模塊變形云圖參見下頁圖5 -圖6。質(zhì)量載荷對甲板位移影響較大。鉆井甲板平面Z向最大位移為2.39 cm，發(fā)生在絞車下部，這是由于絞車質(zhì)量對結(jié)構(gòu)位移產(chǎn)生了較大影響。中層甲板和下層甲板最大位移發(fā)生在跨距最長的主要構(gòu)件處。中層甲板Z向最大位移為2.55 cm，下層甲板Z向最大位移為3.16 cm。由于受到風(fēng)載荷的影響，下底座X方向最大位移為2.76 cm，發(fā)生在東側(cè)和西側(cè)的橫撐上，下底座Y方向最大位移為2.51 cm，發(fā)生在南側(cè)和北側(cè)的斜撐上。由于受BOP架上設(shè)備質(zhì)量影響，下底座Z向最大位移為1.56 cm，發(fā)生在BOP架附近跨距最長的構(gòu)件上。

圖5 鉆井設(shè)備模塊變形云圖

圖6 下底座框架位移云圖

4 結(jié) 論

筆者基于某深海張力腿平臺結(jié)構(gòu)，以慣性力的形式考慮平臺運動對上部結(jié)構(gòu)強度分析產(chǎn)生的影響，并對鉆井設(shè)備模塊進行準靜力分析，既避免了動力分析的巨大計算量，也使靜力分析獲得更接近現(xiàn)實的分析結(jié)果。本文分析計算結(jié)果可為類似張力腿平臺的移動式平臺上部設(shè)備結(jié)構(gòu)的強度分析提供一定的研究思路和技術(shù)支持。

［參考文獻］

［1］AHM AD S K， AHMAD S. Active control of non-line-arly coupled TLP response under wind and wave environments［J］. Computers and Structures， 1994， 72 ：735-747.

［2］張智，董艷秋. 1990 年后世界 TLP 平臺的發(fā)展狀況［J］.中國海洋平臺， 2004（2）：5-11.

［3］渠基順， 管義鋒， 盧燕祥. 基于 SESAM的導(dǎo)管架平臺結(jié)構(gòu)動力響應(yīng)分析［J］. 江蘇船舶， 2013（1）： 1- 5.

［4］“典型深水平臺概念設(shè)計研究”課題組.張力腿平臺水動力響應(yīng)分析 ［J］.中國造船， 2005（s1）：477-489.

［5］NEWMANN J N.船舶流體動力學(xué)［M］.周樹國譯. 北京：人民交通出版社，1986.

［6］CHATTERJEE P C， DAS P K. A hydro-structural analysis program for TLP ［J］. Ocean Engineering，1997（4）：313.

［7］楊樹耕，藤明清，孟昭瑛，等. 有限元分析軟件ANSYS在海洋工程中的應(yīng)用（續(xù)1） 第二講采用ANSYS程序的導(dǎo)管架平臺結(jié)構(gòu)有限元分析［J］. 中國海洋平臺，2000（5）：42-48.

［8］楊樹耕，孟昭瑛，任貴永. 有限元分析軟件ANSYS在海洋工程中的應(yīng)用［J］. 中國海洋平臺，2000（2）：41-44.

［9］API. Specification for Drilling and Well Servicing Structure［S］. Forth Edition，2013.