馮加偉，周立臣，梁思穎

(中海石油(中國)有限公司 湛江分公司 工程建設中心，廣東 湛江 524000)

0 引 言

自升式鉆井船具有移動便利性，被廣泛應用于海洋油氣勘探開發(fā)中。為降低成本，許多海洋油氣生產(chǎn)/處理平臺不設鉆機模塊，在投產(chǎn)前需要用自升式鉆井船進行打井，在投產(chǎn)后需要用自升式鉆井船進行修井。自升式鉆井船在靠近海洋平臺打井/修井作業(yè)時，需要將樁靴插入泥中一定深度，為鉆井船提供承載力。為增加鉆井船操作穩(wěn)定性，鉆井船樁靴一般設計成大樁靴，直徑可達20 m。在鉆井船插樁過程中，樁靴會排擠大量土體，對鄰近的海洋平臺樁基產(chǎn)生影響(額外載荷)，需對此進行研究。

鉆井船樁靴貫入土體對鄰近海洋平臺樁基影響的研究通常有模型試驗研究法和仿真數(shù)值研究法。試驗研究需要一定的條件，付出較高的成本。隨著計算機硬件和有限元分析方法的不斷發(fā)展，土壤本構模型不斷完善，有限元方法等數(shù)值方法能夠有效地模擬各種復雜的樁基貫入問題[1]。國內(nèi)外學者[2-3]對土體大變形等復雜巖土工程問題的計算方法進行大量研究，研究表明采用耦合歐拉-拉格朗日(Coupled Euler-Lagrange，CEL)法計算土體大變形問題獲得的數(shù)值結果與離心模擬試驗結果一致，可運用CEL有限元數(shù)值仿真方法分析實際工程中鉆井船插樁對鄰近海洋平臺樁基的影響。因此，本文用CEL法對鉆井船插樁對海洋平臺樁基的影響進行研究。

1 研究方法

1.1 耦合歐拉-拉格朗日法

鉆井船插樁過程是一個樁靴貫入泥土引起土壤大變形的問題，若運用常規(guī)的有限元算法，則在計算過程中有限元離散網(wǎng)格會發(fā)生嚴重畸變，導致計算中斷，無法模擬整個連續(xù)插樁過程累積的擠土效應對鄰近樁基的影響。

拉格朗日法基于粒子思想，研究一個選定的質(zhì)點在空間運動過程中各物理量的變化規(guī)律；歐拉法基于網(wǎng)格思想，以流場為對象，研究各時刻各流體質(zhì)點在流場中的變化規(guī)律[4]。CEL方法是歐拉有限元方法與拉格朗日有限元方法互相耦合的有限元算法，結合歐拉網(wǎng)格中網(wǎng)格固定而材料可以在網(wǎng)格中自由運動的方式建立模型，歐拉體材料與拉格朗日體材料之間的接觸采用基于罰函數(shù)法的通用接觸算法離散。這種方法基于加強的侵入邊界方法，使拉格朗日結構融入歐拉體中的空單元，拉格朗日體和歐拉體之間的分界面會被自動計算和跟蹤，能夠很好地解決涉及大變形的高度非線性的接觸問題[5-6]。

1.2 數(shù)值模型

建立鉆井船樁靴和海洋平臺鋼樁主要結構模型，鉆井船樁靴直徑為23 m，立面設計如圖1所示。

單位：mm圖1 樁靴立面圖

海洋平臺鋼樁直徑為2.438 m，材料為DH36，樁長為138 m，入泥深度為114.5 m，樁基壁厚隨樁深h變化情況如表1所示。

表1 鋼樁壁厚隨樁深變化

樁靴與鋼樁的相對位置如圖2所示，海洋平臺A1腿的3根裙樁A1-1、A1-2、A1-3中心距樁靴中心的距離分別為24.949 m、18.792 m、18.261 m。

圖2 樁靴與鋼樁相對位置示例

土體模型尺寸為138 m(長)×138 m(寬)×150 m(高)，如圖3所示。將泥面及泥面以下45 m的土體設為歐拉體，余下的土體設為拉格朗日體，采用六面體進行網(wǎng)格劃分。樁靴整體設置成剛體，在泥面之上設置高度為15 m的空穴，為土體提供向上隆起的空間。

圖3 土體有限元模型

鉆井船就位地質(zhì)參數(shù)如表2所示，將土體取為理想彈塑性：砂土層彈性模量取20 MPa，泊松比取0.300，黏聚力取1 kPa；黏土層彈性模量E取150倍的不排水抗剪強度，泊松比取0.495。

表2 地質(zhì)條件設計參數(shù)

根據(jù)樁靴和鋼樁的材料屬性，樁靴和鋼樁的彈性模量為210 GPa，密度為7 850 kg/m3，泊松比為0.250。

1.3 邊界條件

將歐拉土體和空穴在3個方向邊界節(jié)點的水平速度設置為0 m/s，并將邊界設置為歐拉吸收邊界。將拉格朗日土體區(qū)域側面和底面的位移設置為0 m/s。對鋼樁的頂部進行位移約束，僅允許樁發(fā)生豎向位移。

在計算時，采用位移控制方式設定樁靴貫入土層的速率為0.5 m/s，樁靴貫入深度為12 m。樁的加載方式如下：將泥面處樁的截面中心與截面耦合在一起，給截面中心施加載荷。

2 計算結果與分析

2.1 樁靴貫入引起的土體響應

樁靴貫入歐拉土體引起周圍土體擠壓隆起和回淤通過歐拉體積分數(shù)(Eulerian Volume Fraction，EVF)進行顯示。插樁結束后歐拉體積分數(shù)如圖4所示。在樁靴貫入過程中，樁靴不斷擠壓兩側土體使其向外隆起，樁靴上部形成空穴，影響范圍約為2倍樁靴直徑。當貫入深度為12 m時，受群樁的影響，靠近樁側土體的回淤現(xiàn)象比遠離樁側明顯，如圖5所示。

圖4 插樁結束后歐拉體積分數(shù)

圖5 樁靴貫入土體響應

由樁靴貫入土體速度矢量圖(見圖6和圖7)可知，在鉆井船樁靴貫入過程中，樁靴下部的土由底部向兩側運動。在樁坑頂部，靠近海洋平臺鋼樁側泥面處的土體受到鄰近海洋平臺鋼樁的影響，其運動明顯大于遠離鋼樁側。

圖6 樁A1-1和樁A1-2周圍土體速度矢量圖

圖7 樁A1-3周圍土體速度矢量圖

圖8為樁靴貫入不同深度(4 m、8 m、12 m)后樁A1-3周圍土體等效塑性應變分布圖。由圖8可知，由于樁基的存在，隨著樁靴逐步貫入，樁靴兩側的塑性區(qū)由對稱分布變成不對稱分布，這說明樁基影響土體塑性區(qū)的分布，從而對樁基承載力產(chǎn)生影響。

圖8 樁A1-3周圍土體塑性應變

2.2 樁靴貫入對鋼樁的影響

圖9從左至右分別為樁A1-1、A1-2、A1-3在樁靴貫入不同深度(4 m、8 m、12 m)下的側向位移分布圖。由圖9可知，A1-3的位移較大，最大位移出現(xiàn)在樁靴貫入最大深度處。

圖9 樁基側向位移分布

圖10為樁靴貫入不同深度下的樁基應力分布圖。由圖10可知，隨著樁靴貫入深度增大，樁基應力逐漸增大，且最大應力位置向下移動。與樁A1-1和樁A1-2相比，樁A1-3的應力較大。

圖10 樁基應力分布

由計算結果可知，最危險的樁基是最靠近樁靴的樁基，即樁A1-3。對樁靴貫入不同深度(4 m、8 m、12 m)的樁身數(shù)據(jù)進行對比可知，當樁靴安裝就位，即樁靴貫入達12 m時，樁身應力最大，應力最大處為泥面以下20 m左右，可見樁靴貫入對土體存在一定的影響范圍，樁靴對其斜下方處鋼樁應力和水平位移影響最大。

2.3 樁靴貫入對樁基承載力的影響

在鉆井船插樁過程中，樁靴對鄰近海洋平臺樁基周圍土體和鋼樁本身都會產(chǎn)生一定的影響。圖11和圖12為A1-3樁基在樁靴貫入前后樁基軸向力和軸向位移變化分布圖。由圖12可知，樁靴貫入對樁身軸向力的影響有限，軸向位移增加約3 mm。

圖11 插樁前后A1-3樁基軸向力變化

圖12 插樁前后A1-3樁基軸向位移變化

鉆井船插樁結束(入泥深度12 m)后，海洋平臺樁基周圍和鉆井船樁靴下部土體出現(xiàn)明顯的塑性區(qū)(見圖13)，樁基周圍地基土塑性區(qū)從泥面處沿樁身展開，樁靴底部塑性區(qū)由兩側底邊緣展開并貫通。樁基周圍最大展開深度Zmax=12.72 m(泥面處Z=0 m)，樁靴底部塑性區(qū)最大開展深度Zmax=15.14 m。在泥面以下11.90 m，樁周出現(xiàn)最大塑性應變。

圖13 A1-3樁切面塑性區(qū)分布云圖

由塑性區(qū)分布可知，樁靴插入對于樁基周圍土體影響深度為樁靴的插入深度。樁基入泥114.50 m，采用保守的計算方法，塑性區(qū)土體不提供相應的側摩阻力，由此而產(chǎn)生的豎向承載力的降低很小。

綜合樁基樁身軸力變化、樁基豎向位移和塑性區(qū)展開的情況，得到鉆井船插樁對樁基承載力影響較小。

3 結 論

用CEL方法對鉆井船插樁過程進行數(shù)值模擬，得到鉆井船插樁過程引起的土體響應規(guī)律，對鉆井船樁靴對鄰近海洋平臺(距離小于1倍樁靴直徑)樁基的影響進行分析得到：

(1)在鉆井船樁靴貫入過程中，樁靴擠壓兩側土體使其向外隆起，樁靴上部形成空穴，樁靴下部土體向兩側運動，受鄰近平臺樁基影響，靠近樁側土體運動較遠離樁基側明顯，樁靴兩側的塑性區(qū)由對稱分布變成了不對稱分布，靠近樁側土體回淤更快。

(2)鉆井船樁靴貫入對距離樁靴最近的樁基影響最大，隨著鉆井船樁靴貫入深度增大，樁基應力逐漸增大，樁基最大應力位置逐漸下移，鉆井船樁靴斜下方處樁基應力和水平位移所受影響最大。

(3)鉆井船插樁過程對鄰近平臺樁基周圍土體影響深度為樁靴插入深度，樁靴對樁基承載力和軸向位移影響較小。