基于D-P準(zhǔn)則的天然氣水合物井壁穩(wěn)定性研究

鉆井過(guò)程中，鉆井液對(duì)水合物地層產(chǎn)生熱交換作用，同時(shí)鉆井液侵入會(huì)影響水合物穩(wěn)定，兩者均引起井周應(yīng)力場(chǎng)變化，導(dǎo)致井壁坍塌等復(fù)雜情況，嚴(yán)重制約著天然氣水合物資源高效開(kāi)發(fā)。基于線(xiàn)性熱彈性多孔介質(zhì)理論建立了熱流固耦合的水合物斜井井周應(yīng)力模型，引入歐拉變換考慮三維地應(yīng)力大小和方向的隨機(jī)性，采用D-P準(zhǔn)則強(qiáng)化中間主應(yīng)力對(duì)水合物井壁穩(wěn)定的影響，并開(kāi)展了井斜角/方位角、水合物飽和度、井壁滲透性以及鉆井液溫度等因素對(duì)水合物井壁穩(wěn)定的研究。研究結(jié)果表明：井斜角對(duì)井壁的穩(wěn)定性具有更大的影響，井斜角每增加30°坍塌壓力變化0.155 MPa，方位角每增加30°坍塌壓力變化0.112 MPa。水合物從高到低飽和度的分解過(guò)程，前期地層強(qiáng)度緩慢降低，后期快速降低。水合物地層的有效孔隙度增加使得鉆井液侵入量增加，致使地層有效應(yīng)力降低，導(dǎo)致地層的變形和破裂，增加井壁失穩(wěn)風(fēng)險(xiǎn)。水合物地層對(duì)高于相平衡溫度（285 K）的鉆井液較為敏感，每升高1 K坍塌壓力當(dāng)量密度約增加0.006 9 g/cm³，相反則表現(xiàn)較為遲鈍。在進(jìn)行鉆井作業(yè)時(shí)，選用合適的高抑制性、低溫鉆井液有助于控制水合物的分解以及減少鉆井液侵入，進(jìn)而降低水合物井壁失穩(wěn)風(fēng)險(xiǎn)。研究結(jié)果可為天然氣水合物井壁穩(wěn)定性研究提供參考。

天然氣水合物；井壁穩(wěn)定；D-P準(zhǔn)則；坍塌壓力；熱流固耦合

TE52

A

005

Research on Gas Hydrate Wellbore Stability Based on D-P Criterion

Sun Xiaofeng^1，2 Zhao Yuanzhe^1，2 Li Zhipeng³ Hu Qiaobo² Tao Liang^1，2 Yao Di²

（1.SANYA Offshore Oil amp; Gas Research Institute，Northeast Petroleum University；2.Key Laboratory of Enhanced Oil and Gas Recovery of Educational Ministry，Northeast Petroleum University；3.No.4 Drilling Engineering Company，CNPC Bohai Drilling Engineering Company Limited）

During drilling operations，drilling fluid generates heat exchange with hydrate formation，and the invasion of drilling fluid affects the stability of the hydrate.Both situations cause changes in the stress field around the well，leading to complexities such as wellbore collapse，which restrict the efficient development of gas hydrate resources.In the paper，a thermo-fluid-solid coupling circumferential stress model of deviated hydrate well was built according to the theory of linear thermoelastic porous media.Then，by introducing the Euler’s transformation and considering the random magnitudes and directions of 3D in-situ stress，the D-P criterion was used to strengthen the influence of intermediate principal stress on hydrate wellbore stability，and understand the influences of the factors such as borehole inclination/azimuth，hydrate saturation，sidewall permeability and drilling fluid temperature on hydrate wellbore stability.The results show that the inclination has a greater influence on the wellbore stability.For every 30° increase in inclination，the collapse pressure changes by 0.155 MPa；for every 30° increase in azimuth，the collapse pressure changes by 0.112 MPa.In the decomposition process of hydrate from high to low saturation，the formation strength slowly decreases in the early stage and rapidly decreases in the later stage.The increase in effective porosity of hydrate formation leads to an increase in drilling fluid invasion，resulting in a decrease in effective stress of the formation，which further leas to the deformation and fracturing of the formation and an increase in the risk of wellbore instability.Hydrate formation is more sensitive to drilling fluids above the phase equilibrium temperature （285 K）.For every 1 K increase，the equivalent density of collapse pressure is increased by about 0.006 9 g/cm³；conversely，hydrate formation is rather dull to drilling fluids.At the time of drilling，the selection of appropriate high inhibitory and low-temperature drilling fluids can help control the decomposition of hydrates and reduce drilling fluid invasion，thereby reducing the risk of hydrate wellbore instability.The research results provide reference for investigating gas hydrate wellbore stability.

natural gas hydrate；wellbore stability；D-P criterion； collapse pressure；thermo-fluid-solid coupling

0 引 言

天然氣水合物是由天然氣和水合物在高壓低溫下形成的晶體狀化合物，廣泛分布于海洋深水地層中，儲(chǔ)量較大，燃燒效率高，被譽(yù)為新的、有前景的非常規(guī)替代能源^［1］。中國(guó)、日本、加拿大和美國(guó)等國(guó)家均投入巨資進(jìn)行水合物勘探和開(kāi)采方面的研究^［2］，但水合物地層強(qiáng)度低、孔隙度高以及水合物狀態(tài)容易受到溫壓影響發(fā)生相變^［3］，這導(dǎo)致井壁坍塌等復(fù)雜情況頻發(fā)，嚴(yán)重影響了鉆井效率。

在水合物井壁穩(wěn)定性研究方面，國(guó)內(nèi)外學(xué)者進(jìn)行了大量的相關(guān)研究。李慶超等^［4］分析水合物分解對(duì)井眼穩(wěn)定的影響。熊超等^［5］分析鉆井液溫度對(duì)水合物地層的影響。SUN J.X.等^［6］分析鉆井液性質(zhì)和初始儲(chǔ)層條件對(duì)井壁穩(wěn)定的影響。M.KURIHARA等^［7］分析了溫度對(duì)井壁周?chē)鷳?yīng)力的影響。M.SALEHABADI等^［8］研究了鉆井液侵入分解后的井筒穩(wěn)定性以及非均勻載荷作用下水合物井壁垮塌誘因。GUO Z.Y.等^［9］研究了鉆井液性質(zhì)對(duì)水合物影響及水合物分解對(duì)地層強(qiáng)度的影響。以上研究工作多建立于Mohr-Coulomb準(zhǔn)則的基礎(chǔ)上，而且未充分考慮中間應(yīng)力對(duì)水合物井壁穩(wěn)定性的影響。本文在前人研究的基礎(chǔ)上加入歐拉變換，充分考慮三維地應(yīng)力方向的隨機(jī)性，引入線(xiàn)性熱彈性多孔介質(zhì)理論，建立熱流固耦合的斜井井周應(yīng)力模型，采用Drucker Prager（D-P）準(zhǔn)則強(qiáng)化中間主應(yīng)力對(duì)水合物井壁的影響，進(jìn)一步分析井斜角/方位角、水合物飽和度、鉆井液溫度以及地層孔隙度等因素對(duì)水合物地層井壁的影響。研究結(jié)果可為天然氣水合物井壁穩(wěn)定性研究提供參考。

孫曉峰，等：基于D-P準(zhǔn)則的天然氣水合物井壁穩(wěn)定性研究

1 模型建立及強(qiáng)度準(zhǔn)則

1.1 天然氣水合物井周?chē)鷰r應(yīng)力模型

天然氣水合物的應(yīng)力狀態(tài)由巖體本身載荷狀態(tài)和其經(jīng)歷的地質(zhì)歷史應(yīng)力路徑?jīng)Q定，因此其大小、方向具有隨機(jī)性和區(qū)域差異性^［10］。筆者引入歐拉變換實(shí)現(xiàn)地應(yīng)力坐標(biāo)系與井眼坐標(biāo)系的轉(zhuǎn)換，來(lái)描述天然氣水合物斜井井壁圍巖的應(yīng)力狀態(tài)及井壁穩(wěn)定性。

在圖1中，假設(shè)在水合物地層中主地應(yīng)力方向與圖中標(biāo)注方向相同。圖1中：σ_H為水平最大地應(yīng)力，MPa；σ_h為水平最小地應(yīng)力，MPa；σ_v為上覆巖層壓力，MPa。根據(jù)歐拉變換，地應(yīng)力坐標(biāo)系與為井眼坐標(biāo)系轉(zhuǎn)換表達(dá)式如下：

式中：φ為井斜角，（°）；ψ為井斜方位角，（°）；［σ_xyz］為井眼坐標(biāo)系下的地應(yīng)力張量，MPa；［σ_x］為地應(yīng)力，MPa;［M_{（φ，ψ）}］為地應(yīng)力坐標(biāo)系的地應(yīng)力張量，MPa；［M^T_{（φ，ψ）}］為［M_{（φ，ψ）}］的轉(zhuǎn)置矩陣，MPa^-1。

為構(gòu)建相對(duì)簡(jiǎn)單易用的水合物井壁穩(wěn)定模型，對(duì)模型進(jìn)行如下假設(shè)：①假設(shè)水合物地層為多孔彈塑性地層，均質(zhì)、各向同性材料，同時(shí)井壁圍巖處于平面應(yīng)變狀態(tài)；②鉆井過(guò)程中產(chǎn)生的氣體或者水相對(duì)較少，研究中忽略了水合物分解產(chǎn)生氣體或者水的影響。

基于上述假設(shè)并結(jié)合歐拉變換，得到水合物井壁圍巖應(yīng)力分布方程^［11-12］：

式中：σ_r、σ_θ、σ_z分別為距離井眼r處地層所受的徑向應(yīng)力、周向應(yīng)力和垂向應(yīng)力，MPa；α為有效應(yīng)力系數(shù)；p（r）為水合物井壁周?chē)貙拥膶?shí)時(shí)地層孔隙壓力，MPa；p₀為地層初始孔隙壓力，MPa;p_i為液柱壓力，MPa；R為井眼半徑，m；r為與井軸的距離，m；θ為方位角，（°）；ν為泊松比，無(wú)量綱;τ_rθ為井壁上的剪切應(yīng)力，MPa。

在鉆井過(guò)程中，鉆井液在正壓差作用下向地層發(fā)生徑向流動(dòng)，造成巖石內(nèi)部裂縫擴(kuò)展延伸，流體滲流空間變大，產(chǎn)生滲流應(yīng)力變化?；诳紫稄椥粤W(xué)原理^［13］，其相關(guān)表達(dá)式如下：

式中：σ_r^p、σ_θ^p、σ_z^p分別為地層滲流在井壁周?chē)a(chǎn)生的徑向應(yīng)力、周向應(yīng)力和垂向應(yīng)力，MPa；p^f（r）為地層孔隙壓力，MPa；γ為水合物地層孔隙度。

在鉆井過(guò)程中，鉆井液傳熱、鉆頭切削地層巖石及井底鉆具與水合物井壁摩擦都會(huì)產(chǎn)生熱，從而導(dǎo)致水合物發(fā)生分解，進(jìn)而影響水合物地層的膠結(jié)及骨架，使得井壁發(fā)生坍塌?；跓釓椥岳碚摼軣釕?yīng)力的表達(dá)式如下^［14］：

式中：σ^T_r、σ^T_θ、σ^T_z分別為地層溫度變化在井壁周?chē)a(chǎn)生的徑向應(yīng)力、周向應(yīng)力和垂向應(yīng)力，MPa；E為彈性模量，MPa;α_T為水合物地層熱膨脹系數(shù)，K^-1；T（r）為水合物井壁周?chē)貙拥膶?shí)時(shí)地層溫度，K；T₀為水合物地層原始溫度，K。

聯(lián)立式（3）～式（5），建立熱流固耦合的水合物斜井井周應(yīng)力模型如下：

根據(jù)天然氣水合物井周應(yīng)力模型，得到井壁主應(yīng)力為：

式中：δ為滲透系數(shù)，當(dāng)井壁滲透時(shí)δ=1，井壁不滲透時(shí)δ=0；σ_i、σ_j、σ_k分別為井壁上的3個(gè)主應(yīng)力，MPa；σ_θz為周向應(yīng)力在垂向應(yīng)力上的投影，MPa。

1.2 強(qiáng)度準(zhǔn)則

L.VERNIK等^［15］研究發(fā)現(xiàn)，Mohr Coulomb（M-C）準(zhǔn)則預(yù)測(cè)的井壁破壞區(qū)域與實(shí)際不相符，建議使用考慮中間主應(yīng)力影響的破壞準(zhǔn)則。筆者采用Drucker Prager（D-P）準(zhǔn)則^［16］強(qiáng)化中間應(yīng)力對(duì)水合物圍巖強(qiáng)度的影響。D-P準(zhǔn)則不僅考慮巖石的剪切破壞，而且考慮了中間主應(yīng)力以及靜水壓力的作用，同時(shí)具有較好的收斂性，廣泛應(yīng)用在數(shù)值模擬分析軟件中。D-P準(zhǔn)則的通用表達(dá)式如下^［17］：

其中：

式中：I₁為應(yīng)力第一不變量，MPa；J₂為偏應(yīng)力張量的第二應(yīng)力不變量，MPa²；R、K_f為關(guān)于內(nèi)摩擦角和內(nèi)聚力的常數(shù)；C為內(nèi)聚力，MPa。

2 水合物儲(chǔ)層物性參數(shù)

G.P.HOLLAND等^［18］針對(duì)水合物的微觀(guān)結(jié)構(gòu)研究發(fā)現(xiàn)，成藏的水合物沉積物具有很高的滲透性，主要是因?yàn)樘烊粴庖运衔锏男螒B(tài)存在于地層孔隙中，故地層的實(shí)際滲透率和有效孔隙度受到水合物的影響，其相關(guān)表達(dá)式如下^［19］：

式中：k為水合物沉積物有效滲透率，μm²；k₀為水合物沉積物初始滲透率，μm²；γ₀為水合物地層初始孔隙度，%;S_h為水合物飽和度，%。

鉆進(jìn)水合物地層之前，地層處于穩(wěn)定狀態(tài)，同時(shí)各點(diǎn)壓力等于原始地層孔隙壓力。隨著不斷鉆進(jìn)，水合物地層被打開(kāi)，水合物狀態(tài)發(fā)生改變，從而導(dǎo)致地層孔隙壓力發(fā)生變化，其相關(guān)表達(dá)式如下^［20］：

式中：r_∞為距離井眼無(wú)限遠(yuǎn)處，m。

在鉆井過(guò)程中，鉆井液溫度一般高于水合物地層溫度，與地層形成溫度差，從而產(chǎn)生熱量交換，水合物地層溫度發(fā)生改變。為簡(jiǎn)化模型計(jì)算，不考慮摩擦熱同時(shí)假設(shè)套管導(dǎo)熱系數(shù)無(wú)限大，相關(guān)表達(dá)式如下：

式中：T_i為鉆井液溫度，K；λ為地層相對(duì)導(dǎo)熱系數(shù)。

水合物地層彈性模量隨水合物飽和度呈線(xiàn)性變化，其相關(guān)表達(dá)式如下^［21］：

E=E₀+aS_h（20）

式中：E為水合物地層彈性模量，MPa；E₀為初始水合物地層彈性模量，MPa;a為材料參數(shù)，取值0.001 91 MPa。

在水合物開(kāi)采過(guò)程中，鉆井液對(duì)地層的侵入打破了水合物儲(chǔ)層的平衡，水合物發(fā)生分解，導(dǎo)致地層膠結(jié)性變差，內(nèi)聚力改變。其相關(guān)表達(dá)式如下^［22-23］：

式中：C₀為初始水合物儲(chǔ)層內(nèi)聚力，MPa；?為內(nèi)摩擦角，（°）；b為材料參數(shù)，取值1.8。

WANG L.等^［24］研究發(fā)現(xiàn)，含水合物沉積物的強(qiáng)度和彈性模量隨著水合物飽和度的增加而增加，但泊松比影響較小。K.MIYAZAKI等^［25］通過(guò)三軸試驗(yàn)發(fā)現(xiàn)，含水合物巖石的內(nèi)摩擦角不受水合物飽和度的影響。因此在計(jì)算過(guò)程假設(shè)內(nèi)摩擦角以及泊松比不變。

3 水合物井壁穩(wěn)定性分析

通過(guò)文獻(xiàn)調(diào)研^{［13，22］}得到了地層力學(xué)參數(shù)與水合物力學(xué)參數(shù)，該數(shù)值接近墨西哥灣的地層條件，具體數(shù)值見(jiàn)表1。

3.1 井斜角/方位角

采用D-P準(zhǔn)則結(jié)合表1數(shù)據(jù)，得水合物井周坍塌壓力云圖及井壁上最大主應(yīng)力與最小主應(yīng)力的應(yīng)力差值隨方位角的變化規(guī)律，如圖2和圖3所示。坍塌壓力當(dāng)量云圖隨著方位角呈周期性變化，故取1/2周期進(jìn)行分析，即方位角范圍在0°～180°內(nèi)。

由圖2和圖3得，在相同方位角鉆進(jìn)時(shí)，水合物井壁處應(yīng)力差隨著井斜角的增加不斷增加。水合物井壁坍塌壓力不斷增加，故在鉆進(jìn)過(guò)程中，保持水平井段水合物井壁穩(wěn)定，即可保證直井段和斜井段不發(fā)生坍塌。在相同方位下井斜角每變化30°，坍塌壓力當(dāng)量密度變化0.01 g/cm³，即坍塌壓力變化0.155 MPa。在相同井斜段，水合物坍塌壓力當(dāng)量密度云圖隨著方位角的增加呈現(xiàn)“紅-黃-紅”的變化規(guī)律。90°方位角井壁處應(yīng)力差值最低，即朝著水平最小地應(yīng)力的方向鉆進(jìn)較為安全；180° 方位角應(yīng)力差值最高，即朝著水平最大地應(yīng)力的方向鉆進(jìn)會(huì)增加井壁失穩(wěn)風(fēng)險(xiǎn)。當(dāng)井斜角相同時(shí)，方位角每變化30°，坍塌壓力當(dāng)量密度變化0.006 g/cm³，即坍塌壓力變化0.112 MPa。研究發(fā)現(xiàn)，相對(duì)于方位角，井斜角對(duì)水合物井壁穩(wěn)定具有更大的影響。因此，在鉆井過(guò)程中，為了降低井壁坍塌的風(fēng)險(xiǎn)，建議選擇朝著最小水平主應(yīng)力的方向進(jìn)行鉆進(jìn)，并且保持井斜角在一定范圍內(nèi)變化。

3.2 水合物飽和度

結(jié)合方位角與井斜角研究結(jié)果，采用90°方位角的直井段和水平段的數(shù)據(jù)，得到不同水合物飽和度下水合物坍塌壓力變化規(guī)律，如圖4所示。

由圖4可知，隨著水合物飽和度的不斷增加，水合物井壁的坍塌壓力不斷減小。主要是水合物在地層中起到支撐和膠結(jié)的作用^［26］，隨著水合物飽和度的增加，地層逐漸從松散變得牢固，進(jìn)而地層強(qiáng)度不斷增加。水合物飽和度從80%降低到60%，坍塌壓力當(dāng)量密度約增加了0.01 g/cm³，即坍塌壓力增加0.076 MPa。但是水合物飽和度從40%降低到20%，坍塌壓力當(dāng)量密度約增加了0.03 g/cm³，即坍塌壓力增加0.228 MPa。該趨勢(shì)表明，在水合物從高飽和度分解到低飽和的過(guò)程中，前期地層強(qiáng)度緩慢降低，后期地層強(qiáng)度快速降低，因此鉆井過(guò)程中應(yīng)時(shí)刻注意水合物分解情況，減少對(duì)水合物地層的擾動(dòng)，以降低井壁坍塌風(fēng)險(xiǎn)。

3.3 井壁滲透性

當(dāng)水合物地層井眼鉆開(kāi)后，鉆井液首先與井壁發(fā)生復(fù)雜的物理化學(xué)反應(yīng)，使井壁逐漸松散直至貫通，使得井周的孔隙壓力發(fā)生改變，同時(shí)在井壁上產(chǎn)生附加的滲流應(yīng)力，進(jìn)而改變井周應(yīng)力場(chǎng)。研究水合物井壁的滲透性與穩(wěn)定性之間的關(guān)系，可以有效降低水合物井壁失穩(wěn)的潛在風(fēng)險(xiǎn)。水合物井壁滲透性對(duì)水合物坍塌壓力的影響規(guī)律如圖5所示。

由圖5可知，水合物井壁的坍塌壓力與地層的有效孔隙度呈正相關(guān)。隨著地層有效孔隙度的增加，水合物的坍塌壓力當(dāng)量密度也會(huì)增加。水合物井壁的有效孔隙度增加，意味著存在著更多的孔隙和裂縫空間，這會(huì)促使鉆井液更容易進(jìn)入地層中。鉆井液的侵入不僅降低地層的強(qiáng)度和穩(wěn)定性，還可能與地層中的水合物反應(yīng)引起地層結(jié)構(gòu)的改變。同時(shí)有效孔隙度增高會(huì)導(dǎo)致井周附近的地層孔隙壓力升高，使得有效應(yīng)力減小，導(dǎo)致地層的變形和破裂。建議選用具有高度抑制作用的鉆井液^［27］，這樣不僅能有效降低其對(duì)巖石強(qiáng)度的影響，還能保護(hù)儲(chǔ)層并提高水合物井壁的穩(wěn)定性。

3.4 鉆井液溫度

在鉆進(jìn)以及循環(huán)過(guò)程中，鉆井液在壓力的作用下侵入地層，同時(shí)鉆井液與地層存在溫度差，會(huì)與地層產(chǎn)生熱交換，使得水合物地層井壁溫度場(chǎng)發(fā)生變化，導(dǎo)致水合物狀態(tài)不穩(wěn)定，進(jìn)而影響水合物井周的應(yīng)力場(chǎng)分布，應(yīng)力集中發(fā)生區(qū)域，地層發(fā)生塑性變形。不同鉆井液溫度下水合物坍塌壓力變化規(guī)律如圖6所示。

由圖6可知：當(dāng)鉆井液溫度逐漸從275 K升高到285 K（相平衡溫度）時(shí)，水合物井壁坍塌壓力當(dāng)量密度緩慢增加，每升高1 K坍塌壓力當(dāng)量密度約增加0.001 5 g/cm³，相當(dāng)于坍塌壓力增加了0.022 MPa，其增加量可忽略不計(jì);當(dāng)鉆井液溫度逐漸高于285 K時(shí)，水合物井壁坍塌壓力當(dāng)量密度快速增加，每升高1 K坍塌壓力當(dāng)量密度約增加0.006 9 g/cm³，相當(dāng)于坍塌壓力增加了0.105 MPa，其增加量為低溫（低于285 K）的4倍。其原因是當(dāng)鉆井液溫度高于地層溫度時(shí)，地層中固體水合物受到溫度影響向流體轉(zhuǎn)化^［28］，使得井壁穩(wěn)定性急劇下降。在壓力作用下進(jìn)入地層的鉆井液會(huì)逐漸增加，干擾的水合物區(qū)域會(huì)逐漸增加，井壁失穩(wěn)的區(qū)域也隨著增加。故在鉆進(jìn)過(guò)程中，根據(jù)地層溫度選擇合適的鉆井液能有效降低井壁失穩(wěn)的風(fēng)險(xiǎn)，減小井壁失穩(wěn)區(qū)域。

4 結(jié) 論

（1）基于線(xiàn)性熱彈性多孔介質(zhì)理論以及歐拉變換，建立了水合物斜井井壁的熱流固耦合模型，分析了井斜/方位角、水合物飽和度、溫度以及井壁滲透性對(duì)水合物井壁穩(wěn)定性的影響。

（2）研究發(fā)現(xiàn)，相對(duì)于方位角，井斜角對(duì)水合物井壁穩(wěn)定具有更大的影響，井斜角每變化30°，坍塌壓力變化0.155 MPa，方位角每增加30°，坍塌壓力變化0.112 MPa。水合物分解前期對(duì)地層強(qiáng)度的影響可忽略不計(jì)，但隨著水合物分解的不斷進(jìn)行，地層強(qiáng)度快速降低，當(dāng)水合物分解至40%以下，水合物每減少20%，坍塌壓力當(dāng)量密度約增加了0.03 g/cm³，即坍塌壓力增加0.228 MPa。水合物地層的有效孔隙度增加使得鉆井液侵入量增加，致使地層有效應(yīng)力降低，導(dǎo)致地層的變形和破裂，進(jìn)而增加井壁失穩(wěn)風(fēng)險(xiǎn)。當(dāng)鉆井液溫度逐漸高于相平衡溫度（285 K）時(shí)，水合物井壁坍塌壓力當(dāng)量密度快速增加，每升高1 K坍塌壓力當(dāng)量密度約增加0.006 9 g/cm³，相當(dāng)于坍塌壓力增加了0.105 MPa，其增加量為低溫（低于285 K）的4倍。

（3）在水合物地層鉆井過(guò)程中，應(yīng)選用合適的低溫、高抑制性鉆井液，防止地層中固體水合物受到溫度影響向流體轉(zhuǎn)化以及減少對(duì)鉆井液侵入對(duì)井周應(yīng)力場(chǎng)的影響，降低水合物井壁坍塌的風(fēng)險(xiǎn)。

［1］ XU C G，LI X S，YAN K F，et al.Research progress in hydrate-based technologies and processes in China：a review［J］.Chinese Journal of Chemical Engineering，2019，27（9）：1998-2013.

［2］ LIU L P，SUN Z L，ZHANG L，et al.Progress in global gas hydrate development and production as a new energy resource［J］.Acta Geologica Sinica-English Edition，2019，93（3）：731-755.

［3］ GUO X W，XU L，WANG B，et al.Optimized gas and water production from water-saturated hydrate-bearing sediment through step-wise depressurization combined with thermal stimulation［J］.Applied Energy，2020，276：115438.

［4］ 李慶超，王威，王義頃，等．水合物儲(chǔ)層鉆井水合物分解及井眼坍塌分析［J］．石油機(jī)械，2021，49（6）：43-50．

LI Q C，WANG W，WANG Y Q，et al.et al hydrate dissociation and wellbore collapse analysis in hydrate reservoir drilling［J］.China Petroleum Machinery，2021，49（6）：43-50.

［5］ 熊超，劉力，徐小龍，等．隔熱套管抑制水合物地層分解規(guī)律研究［J］．石油機(jī)械，2021，49（3）：65-71．

XIONG C，LIU L，XU X L，et al.Research on insulation casing inhibiting the decomposition of hydrate formation［J］.China Petroleum Machinery，2021，49（3）：65-71.

［6］ SUN J X，NING F L，LEI H W，et al.Wellbore stability analysis during drilling through marine gas hydrate-bearing sediments in Shenhu area：a case study［J］.Journal of Petroleum Science and Engineering，2018，170：345-367.

［7］ KURIHARA M，SATO A，F(xiàn)UNATSU K，et al.Analysis of production data for 2007/2008 mallik gas hydrate production tests in Canada［C］∥International Oil and Gas Conference and Exhibition in China.Beijing，China：SPE.2010：SPE 132155-MS.

［8］ SALEHABADI M，JIN M，YANG J H，et al.Finite element modeling of casing in Gas-Hydrate-Bearing sediments［J］.SPE Drilling amp; Completion，2009，24（4）：545-552.

［9］ GUO Z Y，WANG H N，JIANG M J.Elastoplastic analytical investigation of wellbore stability for drilling in methane hydrate-bearing sediments［J］.Journal of Natural Gas Science and Engineering，2020，79：103344.

［10］ 丁立欽，王志喬，呂建國(guó)，等．基于圍巖本體Mogi-Coulomb強(qiáng)度準(zhǔn)則的層理性巖層斜井井壁穩(wěn)定模型［J］．巖石力學(xué)與工程學(xué)報(bào)，2017，36（3）：622-632．

DING L Q，WANG Z Q，LYU J G，et al.A model for inclined borehole stability in bedding rocks based on Mogi-Coulomb criterion of rock matrix［J］.Chinese Journal of Rock Mechanics and Engineering，2017，36（3）：622-632.

［11］ MA T S，CHEN P，YANG C H，et al.Wellbore stability analysis and well path optimization based on the breakout width model and Mogi-Coulomb criterion［J］.Journal of Petroleum Science and Engineering，2015，135：678-701.

［12］ AL-AJMI A M，ZIMMERMAN R W.A new well path optimization model for increased mechanical borehole stability［J］.Journal of Petroleum Science and Engineering，2009，69（1/2）：53-62.

［13］ CHENG W，NING F L，SUN J X，et al.A porothermoelastic wellbore stability model for riserless drilling through gas hydrate-bearing sediments in the Shenhu area of the South China Sea［J］.Journal of Natural Gas Science and Engineering，2019，72：103036.

［14］ GHASSEMI A，TAO Q，DIEK A.Influence of coupled chemo-poro-thermoelastic processes on pore pressure and stress distributions around a wellbore in swelling shale［J］.Journal of Petroleum Science and Engineering，2009，67（1/2）：57-64.

［15］ VERNIK L，ZOBACK M D.Estimation of maximum horizontal principal stress magnitude from stress-induced well bore breakouts in the Cajon pass scientific research borehole［J］.Journal of Geophysical Research，1992，97（B4）：5109-5119.

［16］ SP?TH M，HERRMANN C，PRAJAPATI N，et al.Multiphase-field modelling of crack propagation in geological materials and porous media with drucker-prager plasticity［J］.Computational Geosciences，2021，25（1）：325-343.

［17］ CHEN Y，SUN T，ZHANG Y D，et al.The suitable strength criterion to determine the collapse of hydrate reservoirs with different saturation［C］∥The 29th International Ocean and Polar Engineering Conference.Honolulu，Hawaii，USA：ISOPE，2019：36.

［18］ HOLLAND G P，JENKINS J E，CREAGER M S，et al.Solid-state NMR investigation of major and minor ampullate spider silk in the native and hydrated states［J］.Biomacromolecules，2008，9（2）：651-657.

［19］ GAO Q，CHENG Y，HAN S，et al.Numerical modeling of hydraulic fracture propagation behaviors influenced by pre-existing injection and production wells［J］.Journal of Petroleum Science and Engineering，2019，172：976-987.

［20］ WANG H N，CHEN X P，JIANG M J，et al.Analytical investigation of wellbore stability during drilling in marine methane hydrate-bearing sediments［J］.Journal of Natural Gas Science and Engineering，2019，68：102885.

［21］ YONEDA J，MASUI A，KONNO Y，et al.Mechanical properties of hydrate-bearing turbidite reservoir in the first gas production test site of the eastern nankai trough［J］.Marine and Petroleum Geology，2015，66，Part 2：471-486.

［22］ FREIJ-AYOUB R，TAN C E，CLENNELL B，et al.A wellbore stability model for hydrate bearing sediments［J］.Journal of Petroleum Science and Engineering，2007，57（1/2）：209-220.

［23］ TAN T，ZHANG H，LI J，et al.Analysis on collapse pressure and fracture pressure of a borehole in natural gas hydrate formation［C］∥53rd U.S.Rock Mechanics/Geomechanics Symposium.New York：ARMA，2019：ARMA-2019.

［24］ WANG L，LI Y，SHEN S，et al.Mechanical behaviours of gas-hydrate-bearing clayey sediments of the South China Sea［J］.Environmental Geotechnics，2022，9（4）：210-222.

［25］ MIYAZAKI K，TENMA N，AOKI K，et al.Effects of confining pressure on mechanical properties of artificial methane-hydrate-bearing sediment in triaxial compression test［J］.International Journal of Offshore and Polar Engineering，2011，21（2）：ISOPE-11-21-2-148.

［26］ HYODO M，YONEDA J，YOSHIMOTO N，et al.Mechanical and dissociation properties of methane hydrate-bearing sand in deep seabed［J］.Soils and Foundations，2013，53（2）：299-314.

［27］ 明瑞卿，張時(shí)中，王越之．塔中西部二疊系井壁失穩(wěn)原因分析及對(duì)策［J］．特種油氣藏，2015，22（3）：138-140．

MING R Q，ZHANG S Z，WANG Y Z.Borehole lnstability in permian system in the west of central" Tarim basin-reasons and solutions［J］.Special Oil amp; Gas Reservoirs，2015，22（3）：138-140.

［28］ LIAO Y Q，WANG Z Y，CHAO M Z，et al.Coupled wellbore-reservoir heat and mass transfer model for horizontal drilling through hydrate reservoir and application in wellbore stability analysis［J］.Journal of Natural Gas Science and Engineering，2021，95：104216.第一

孫曉峰，教授，博士生導(dǎo)師，生于1980年，2014年畢業(yè)于東北石油大學(xué)天然氣工程專(zhuān)業(yè)，獲博士學(xué)位，現(xiàn)從事油氣井流體力學(xué)與井眼清潔的研究工作。地址：（163318）黑龍江省大慶市。電話(huà)：（0459）6503521。email：suneye@126.com。

2024-02-25

劉 鋒