陶磊 趙虹宇 鄧嵩 朱碩 耿浩然 王相 紀(jì)夢(mèng)楠

常州大學(xué)

深水淺部地層的膠結(jié)質(zhì)量差，固井難度大，海底地震來臨時(shí)，對(duì)表層套管柱和海底管線可能會(huì)造成較大的破壞。據(jù)統(tǒng)計(jì)美國(guó)威名頓、中國(guó)大慶以及俄羅斯薩格勒布的油藏有超過1 000口由于地震引起的套管損毀井[1]。在鉆井過程中，水下井口-導(dǎo)管面臨著隔水管張力、淺層土耦合作用以及地震載荷的影響[2]。

深水水下井口-導(dǎo)管與海底土耦合振動(dòng)的研究大多以樁基力學(xué)為基礎(chǔ)，最典型的樁土模型是Winker和Novak模型[3]，但這兩個(gè)模型都忽略了土層中縱向應(yīng)力的變化。胡昌斌等建立了一個(gè)三維軸對(duì)稱模型[4]；朱幸科等研究了深水采油套管柱的縱向振動(dòng)[5]，但沒有考慮地震力的影響；ELSAYED等進(jìn)行了地震力下導(dǎo)管架平臺(tái)的分析，其沒有涉及到水下井口-導(dǎo)管的分析[6]；有部分學(xué)者做了地震活動(dòng)下海底管線的動(dòng)態(tài)規(guī)律研究，但沒考慮縱向上與海底土的耦合作用[7]。綜上，水下井口-導(dǎo)管振動(dòng)分析還處在初步階段，大多數(shù)沒有考慮到耦合作用以及地震力的影響。本文建立了一個(gè)考慮地震力的水下井口-導(dǎo)管與海底土的縱向振動(dòng)耦合模型，在模型中同時(shí)考慮了井口上段隔水管的張力、鉆井液上返引起的縱向摩擦力以及地震力，將海底淺層土假設(shè)為固液兩相的飽和黏性土，運(yùn)用Flourier變換得到了水下井口頻域內(nèi)的振動(dòng)特征，并將計(jì)算結(jié)果與ABAQUS模擬結(jié)果以及單相土的計(jì)算結(jié)果對(duì)比，準(zhǔn)確度較高。

1 模型建立

1.1 模型假設(shè)

模型假設(shè)如下：①海底土是均勻、各向同性、具有黏彈性的飽和土；②海底土壤介質(zhì)由土壤骨架和孔隙流體組成；③隔水管導(dǎo)管材料是彈性材料，上返鉆井液是黏彈性流體；④導(dǎo)管和土形成的系統(tǒng)間的縱向位移很小，系統(tǒng)接觸面變形和力是連續(xù)的；⑤地震力為海底土底部瞬時(shí)的簡(jiǎn)諧激振力；⑥導(dǎo)管的BOP頂部有一個(gè)由隔水管縱向張力引起的向上的固定力。

耦合作用簡(jiǎn)化模型如圖1所示。

圖1 水下井口-導(dǎo)管與海底土耦合作用模型Fig.1 Coupling model of underwater wellhead-pipe and seafloor soil

1.2 控制方程

（1） 海底土層縱向振動(dòng)方程。海底土層為固液兩相飽和黏性土，依據(jù)土動(dòng)力學(xué)原理，建立軸對(duì)稱柱坐標(biāo)系的動(dòng)態(tài)平衡方程如下

其中

式中：ur、uz分別為土壤骨架的徑向和垂直位移，mm；wr、wz分別為孔隙流體的徑向和垂直位移，mm； ρs、ρf分別為土壤骨架和孔隙流體的密度，g/cm3； ks、kf、kb分別為土壤骨架、孔隙流體和套管柱的體積彈性模量，MPa；λ、G為土壤的拉梅常數(shù)；c為土壤的黏性阻尼；ρ為飽和土壤的質(zhì)量密度，g/cm3；n為土壤的孔隙度。

當(dāng)發(fā)生地震時(shí)[8]

式中：β為土壤的阻尼率；ω為地震的角頻率。

（2）水下井口-導(dǎo)管的縱向振動(dòng)方程。水下井口-導(dǎo)管系統(tǒng)在縱向上受到隔水管張緊力、土壤耦合摩擦力以及導(dǎo)管內(nèi)鉆井液上返流動(dòng)摩擦力的作用，振動(dòng)方程為

式中：Ec為套管柱截面的楊氏模量，GPa；uc為井口-導(dǎo)管系統(tǒng)的縱向位移，m；τ1為土壤對(duì)導(dǎo)管的垂向阻力，N；τ2為導(dǎo)管內(nèi)上返流體的向上摩擦力，N； f為隔水管的張緊力，假設(shè)為恒定力，N；r1、r2分別為導(dǎo)管的內(nèi)徑和外徑，mm。

1.3 邊界條件

（1）海底土邊界條件。海底土層頂部的位移為0，即

海底土層底部位移為地震力的振幅，即

海底土層位移在水平方向上無窮大

海底土層的初始狀態(tài)為

（2）水下井口-導(dǎo)管邊界條件。系統(tǒng)頂部和底部的位移為

水下井口-導(dǎo)管的初始狀態(tài)為

（3）水下井口-導(dǎo)管與海底土耦合條件。海底土壤與導(dǎo)管之間位移連續(xù)，方程為

海底土壤與導(dǎo)管之間作用力相反，方程為

2 模型求解

當(dāng)海底發(fā)生地震時(shí)，可以將海底地震簡(jiǎn)化為縱向上瞬時(shí)的簡(jiǎn)諧力，振幅為eiωt，此時(shí)海土層的振幅也為 eiωt。

2.1 土層振動(dòng)問題求解

首先引入勢(shì)函數(shù)，令

對(duì)式（1）進(jìn)行Laplace變換得

可以得到方程的一般解為

其中

代入海底土層的邊界條件，當(dāng)r趨向于無窮大時(shí)，應(yīng)力和應(yīng)變?yōu)?，可以得到下式

C1、C2、C3、C4、C5、C6是待定系數(shù)，它們之間的關(guān)系可以表示為

待定系數(shù)可以通過邊界條件求得，代入海底土層的位移條件、導(dǎo)管與海底土界面條件可求得

其中

2.2 井口-導(dǎo)管振動(dòng)問題求解

從前述可得到

用Herschel和Bulkley模型描述鉆井液剪切力[9]，因此τ2可表示為

將方程（18）和（19）帶入方程（7），根據(jù)Laplace變換，可得

方程（20）的解可表示為

其中，cosh(hnz)在區(qū)間[0 ， θ]內(nèi)表示為

可求得C1n，因此方程（21）的解為

3 算例分析

3.1 計(jì)算數(shù)據(jù)

用南海某井?dāng)?shù)據(jù)以及海底土壤性質(zhì)對(duì)深水導(dǎo)管與海底土縱向耦合振動(dòng)進(jìn)行計(jì)算分析。通常導(dǎo)管的入泥深度不能過大也不能太小，不然會(huì)導(dǎo)致承載力不足或耗時(shí)過多，一般為30～100 m，取導(dǎo)管下入深度為60 m，其余具體計(jì)算數(shù)據(jù)如表1所示[10]。

3.2 計(jì)算結(jié)果

將海底土假設(shè)為單相，與本文的計(jì)算結(jié)果對(duì)比可知，井口剛度和阻抗振幅與單相土計(jì)算結(jié)果基本吻合，剛度隨著頻率的增加而減小，阻抗隨著振動(dòng)頻率增加而增加；當(dāng)振動(dòng)頻率較低時(shí)，兩相海底土的振動(dòng)阻抗小于單相土（圖2）。同時(shí)運(yùn)用Abaqus軟件，將本文計(jì)算結(jié)果與Abaqus模擬結(jié)果進(jìn)行對(duì)比（圖3），兩者基本吻合。

表1 深水井及海底淺層土性質(zhì)數(shù)據(jù)Tab.1 Performance data of deepwater wells and seafloor shallow soil

3.3 影響因素分析

研究了剪切模量、孔隙度、導(dǎo)管長(zhǎng)度、鉆井液上返速度以及鉆井船縱向的振幅等參數(shù)對(duì)水下井口-導(dǎo)管振動(dòng)阻抗的影響，結(jié)果如圖4～圖8所示。

圖4所示為土壤剪切模量分別為20和10 MPa時(shí)井口-導(dǎo)管的剛度和阻抗振蕩幅度，從圖中可以看出，隨著剪切模量的增大，剛度和阻尼振蕩幅度都減小。圖5所示為土壤孔隙度對(duì)剛度和阻抗的影響，當(dāng)孔隙度較大時(shí)，剛度和阻尼振蕩幅度增大；剪切模量和孔隙度對(duì)井口-導(dǎo)管縱向振動(dòng)的影響較小。圖6所示為導(dǎo)管長(zhǎng)度對(duì)井口-導(dǎo)管振動(dòng)的影響，剛度和阻尼的變化幅度隨著導(dǎo)管長(zhǎng)度增加而減小，而且較為明顯。圖7所示為鉆井液上返速度對(duì)井口-導(dǎo)管的影響，變化幅度隨著鉆井液上返速度增加而增大，影響較小。從圖8可以看出，鉆井平臺(tái)或鉆井船縱向波動(dòng)幅度對(duì)井口-導(dǎo)管剛度和阻尼影響較大，且隨著縱向振幅的增加而增大。

圖2 兩相土與單相土對(duì)比Fig.2 Comparison of two phase soil and single phase soil

圖3 本文結(jié)果與Abaqus模擬結(jié)果對(duì)比Fig.3 Comparison of the results of this paper and the solutions of Abaqus

圖4 土壤剪切模量對(duì)振動(dòng)阻抗的影響Fig.4 Effect of the shear modulus of soil on vibration impedance

圖5 土壤孔隙度對(duì)振動(dòng)阻抗的影響Fig.5 Effect of the porosity of soil on vibaration impedance

圖6 導(dǎo)管長(zhǎng)度對(duì)振動(dòng)阻抗的影響Fig.6 Effect of the length of conductor on vibration impedance

圖7 鉆井液上返速度對(duì)振動(dòng)阻抗的影響Fig.7 Effect of the mud rising rate on vibration impedance

圖8 鉆井平臺(tái)或鉆井船振幅對(duì)振動(dòng)阻抗的影響Fig.8 Effect of the amplitude of drilling plate or vessel on vibration impedance

4 結(jié)論

（1）導(dǎo)管下入深度越深，井口-導(dǎo)管的軸向力產(chǎn)生的縱向振動(dòng)降低的幅度越大。

（2）地震頻率越大，水下井口-導(dǎo)管的縱向振動(dòng)幅度增加越大。

（3）鉆井液上返速度對(duì)井口-導(dǎo)管的縱向振動(dòng)影響較小。

（4）鉆井船或者平臺(tái)的縱向振動(dòng)對(duì)井口-導(dǎo)管的縱向振動(dòng)有影響。

（5）海底土的剪切模量和孔隙度都將影響水下井口-導(dǎo)管的縱向動(dòng)態(tài)性能。

（6）由于深水海底土膠結(jié)差，存在低溫高壓等特殊現(xiàn)象，造成對(duì)鉆井設(shè)備特性的要求有所不同；同時(shí)地震力也影響鉆井管柱的動(dòng)態(tài)性能。

本文的研究對(duì)于深水、超深井鉆完井作業(yè)具有重要的意義，但也存在著一些不足，如海流力、隔水管軸向力的變化，溫度、壓力的變化等對(duì)管柱動(dòng)態(tài)性能的影響都需要進(jìn)一步研究。