深水測試管柱動力學分析

摘 要

摘要：深水測試對于及時發(fā)現(xiàn)和準確評價海洋深水油氣藏具有重要意義，而深水測試的成功率和效果則是通過深水測試管柱力學研究的結(jié)果來評價的，其中測試管柱動力學分析可用以評價測
摘要：深水測試對于及時發(fā)現(xiàn)和準確評價海洋深水油氣藏具有重要意義，而深水測試的成功率和效果則是通過深水測試管柱力學研究的結(jié)果來評價的，其中測試管柱動力學分析可用以評價測試管柱在各種海況下的安全性。為此，提出了使用幅值響應算子(RAO)來計算鉆井船縱蕩水平運動和升沉運動的位移隨時間響應的方法。把通過大鉤補償系統(tǒng)補償后的升沉運動載荷作為測試管柱頂端的力邊界條件和鉆井船的縱蕩水平運動位移作為測試管柱頂端的位移邊界條件，應用有限元方法對測試管柱進行了動力學分析。研究得出以下結(jié)論：①在波浪的作用下，鉆井船的升沉運動在3m內(nèi)，而鉆井船的縱蕩位移較大；②隨著波高的增大，測試管柱的位移和應力波動幅度增大，而應力平均值則變化不大；③波高在10m以內(nèi)，測試管柱是安全可靠的。
關鍵詞：海洋深水油氣藏；測試管柱；鉆井船；波浪；動力學；應力；安全性；效果
    旨在發(fā)現(xiàn)海洋深水油氣藏的深水測試的成功率和效果是通過深水測試管柱力學研究結(jié)果來評價的。測試管柱動力學分析是在分析測試管柱所受的環(huán)境載荷及其約束條件下，通過有限元軟件建立深水測試管柱有限元模型，進而確定在不同的環(huán)境條件下測試管柱的動力響應的分析方法。
1 鉆井船運動的模擬
1.1 對海洋波浪的模擬
   實際海面上的波浪極為復雜，海面上經(jīng)常出現(xiàn)的波浪是由風產(chǎn)生的風浪。風的方向是多變的。任意不規(guī)則的隨機波浪模擬可以簡化為將若干個不同周期、不同波高、不同相位角的簡單正弦波進行線性疊加。P-M譜代表充分成長的風浪譜，在海浪研究及有關工程問題中得到廣泛應用。國際船模試驗室會議(ITTC)確定的標準波浪譜—P-M譜可以用[1]下式表達：
 
式中Sη(ω)表示P-M單邊譜密度；叫表示圓頻率，1/s；H1/3為有效波高(即波浪總數(shù)中1/3的最高波浪的平均波高)，m。
使用描述波浪的郎格-赫金斯(Longuet-Higgins)模型，波浪可由若干個不同周期、不同波高、不同相位角的簡單余弦波進行線性疊加來模擬。長峰波浪可以表示為：
 
式中η(x，t)表示x處的水面相對于靜水面在t時刻的高度，m；ai為第i個組成波的波幅，m；ki表示t時刻在這個水面高度的波的數(shù)量，就是第i個組成波的數(shù)量；ωi表示第i個組成波的圓頻率，1/s；φi表示第i個組成波的初始相位，(°)。
對于單個波的波幅可以用下式來確定：
 
式中ai表示單個諧波的波幅，m；Sηη(ω)表示所需模擬的波浪頻譜，也稱為靶譜。頻譜范圍ωL～ωH的選取取決于所要求的精度。對于靶譜為P-M譜(本文取千分之二，即μ=0.002)，ωL與ωH則與可按下式確定[2]：
 
    由于波浪是窄帶隨機過程，頻率很小和頻率很大的組成波提供能量很小，能量主要由狹窄頻率帶內(nèi)的組成波提供。因此，隨意地增大頻譜范圍反而會使模擬的精度下降。
1.2 鉆井船運動的計算
    鉆井船運動可分解為：水平縱蕩運動和升沉運動。
    為了模擬鉆井船的水平縱蕩運動，模型必須是平均鉆井船偏移，長期慢漂運動，以及鉆井船對不規(guī)則波浪的瞬時響應三者的疊加。鉆井船水平方向運動數(shù)學模型如式(6)所示[3～6]：
 
式中S(t)表示鉆井船運動響應，為時間t的函數(shù)，m；S0表示鉆井船平均偏移，m；SL表示鉆井船慢漂運動的幅值，m；Sn為第n個組成波的波幅，m；an表示鉆井船與波浪的相位差，(°)。
鉆井船長期慢漂運動表示為：
 
式中SL表示鉆井船慢漂的單邊幅值，m；TL表示鉆井船慢漂運動的周期，s；aL表示慢漂運動于波浪之間的相位差(其值常取為0)，(°)。
    鉆井船的響應幅值算子(Response Amplitude Operator，簡稱RAO)可以用來模擬鉆井船對特定波的響應。響應幅值算子包括幅值傳遞函數(shù)與相位傳遞函數(shù)。前者確定了不同波浪周期時，鉆井船響應與對應組成波的幅值比；后者確定了不同波浪周期時，鉆井船響應與對應組成波的相位差。圖1中是鉆井船縱蕩RAO幅值比和相位差。
    在鉆井船的升沉運動過程中，隨機海浪強迫船體或平臺做受迫振動。因此鉆井船的升沉運動同樣可由余弦波疊加而得到。根據(jù)鉆井船的升沉RA0可以間接得到鉆井船的升沉運動響應。API規(guī)范給出的鉆井船升沉RA0幅值比和相位差見圖2。
 

圖3是在10m波高波浪作用下鉆井船升沉位移和縱蕩位移隨時間的變化?？紤]鉆井船的平均偏移為14m，假設慢漂運動的單邊幅值為10m，慢漂運動的周期為200s，波浪的有效波高為10m。
2 測試管柱動力學計算
以我國南海某深水井為例，其水深為1400m，測試管柱外徑為Φ114.3m，壁厚為9.897mm，大鉤升沉補償機構的系數(shù)為7.88kN/m，鉆井液密度為1.05×103kg/m3，彈性模量取2.5×105MPa，泊松比取0.3，管柱重量為248N/m。
采用有限元分析方法，將鉆井船的升沉響應時間歷程通過鉆井船的懸掛系統(tǒng)轉(zhuǎn)換為軸向力的時間歷程并將其作為測試管柱頂端的力邊界條件，將鉆井船的水平縱蕩響應位移時間歷程作為測試管柱頂端的位移邊界條件。通過有限元動力學模擬得到了波高為10m的測試管柱應力隨時間變化的曲線。圖4中反映了深水測試管柱10m波高波浪作用下離轉(zhuǎn)盤面8m和1300m處應力隨時間的變化情況。圖5是不同波高下深度為8m處測試管柱平均應力與應力波動對比圖。
 
   從圖4、5可知：隨著波高的增大，測試管柱的應力的波動增大，而且測試管柱應力波動的幅度也增大；在不同波高下測試管柱應力的平均值則基本不變。
3 結(jié)論
   1) 在波浪的作用下，鉆井船的升沉運動在3m內(nèi)，而鉆井船的縱蕩位移較大。
    2) 隨著波高的增大，測試管柱的位移和應力波動幅度增大，而應力平均值變化不大。
    3) 波高在10m以內(nèi)，測試管柱是安全的。
參考文獻
[1] 竺艷蓉.海洋工程波浪力學[M].天津：天津大學出版社，1991.
[2] ROBERT M，AGBEZUGE L K.Random wave and vessel motion effects on drilling riser dynamics[C]∥paper 2650 MS presented at the Offshore Technology Conference.3-6 May，1976，Houston，Texas：Texas：Offshore Technology Conference，1976.
[3] 俞聿修.隨機海浪及其工程應用[M].大連：大連理工大學出版社，2000.
[4] 李軍強，劉宏昭，何欽象，等.波浪力作用下海洋鉆井隔水管隨機振動研究[J].機械科技與技術，2004，23(1)：7-10.
[5] 石曉兵，郭昭學，聶榮國，等.海洋深水鉆井隔水管變形及載荷分布規(guī)律研究[J].天然氣工業(yè)，2004，24(3)：88-90.
[6] 石曉兵，陳平.三維載荷對海洋深水鉆井隔水管強度的影響分析[J].天然氣工業(yè)，2004，24(12)：86-88.
 
(本文作者：謝鑫1 付建紅1 張智1 何玉發(fā)2 1.“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室·西南石油大學；2.中海石油研究中心)