CN104806233B - 一种预测弱面地层坍塌压力当量密度窗口的方法 - Google Patents

一种预测弱面地层坍塌压力当量密度窗口的方法 Download PDF

Info

Publication number
CN104806233B
CN104806233B CN201510077224.7A CN201510077224A CN104806233B CN 104806233 B CN104806233 B CN 104806233B CN 201510077224 A CN201510077224 A CN 201510077224A CN 104806233 B CN104806233 B CN 104806233B
Authority
CN
China
Prior art keywords
rock
stratum
well
weak
equivalent density
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510077224.7A
Other languages
English (en)
Other versions
CN104806233A (zh
Inventor
卢运虎
金衍
陈勉
刘铭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201510077224.7A priority Critical patent/CN104806233B/zh
Publication of CN104806233A publication Critical patent/CN104806233A/zh
Application granted granted Critical
Publication of CN104806233B publication Critical patent/CN104806233B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Abstract

本发明公开了一种预测弱面地层坍塌压力当量密度窗口的方法，按照先后顺序包括以下步骤：根据弱面地层特点，将其分为致密段和裂缝段；利用岩石三轴压缩试验测试岩石本体强度、岩石弱面强度；利用流固耦合模型反演致密型地层井周围岩的应力分布和孔隙压力分布，结合摩尔库伦破坏准则确定致密型地层坍塌压力当量密度的下限；利用岩石弱面破坏准则分析弱面地层井周围岩的破坏状态，确定裂缝型地层坍塌压力当量密度窗口的下限和上限；建立弱面地层坍塌压力当量密度窗口的计算模型，将各参数值代入模型，即可确定弱面地层坍塌压力当量密度窗口。本发明可在钻井前预测到弱面地层坍塌压力当量密度窗口，以便合理选择钻井液密度，有效阻止井壁失稳的发生。

Description

一种预测弱面地层坍塌压力当量密度窗口的方法
技术领域
本发明属于油气钻探工程技术领域，涉及预测地层坍塌压力窗口的方法，尤其涉及一种预测弱面地层坍塌压力当量密度窗口的方法。
背景技术
我国大部分油气资源集中在裂缝型储层，此类地层钻井是制约我国石油钻探的关键技术难题。由于地层受构造挤压变形等作用，形成多套裂缝体系，且区域构造应力差较大，存在高压水层，因此导致钻井过程中井壁失稳现象频发。当钻井液密度较低时，无法支撑井壁围岩，从而导致井壁失稳；当钻井液密度较高时，钻井液侵入裂缝内部，降低壁面正应力和缝内粘聚力，同样导致井壁失稳。这表明弱面地层坍塌压力自身存在当量密度窗口，若钻井液密度选择不当，将给钻井作业带来不同程度的人力、物力损失，因此在钻井前事先预测钻井液密度是非常必要的，它可以有效阻止井壁失稳，防止井下复杂情况的发生。
目前，人们采用多种方法预测钻井前地层坍塌压力，仅中国石油大学就申请了多篇相关专利，比如，公开号为1966934A的发明专利公开了一种随钻预测钻头底下地层坍塌压力的方法，包括：提取待钻井和与该待钻井相邻的已钻井旁若干道地震记录，加权处理得到该待钻井和已钻井的地震记录；对已钻井进行声波时差和密度测井，得到已钻井不同层段的测井数据；利用测井数据和地震记录，建立预测地层声波速度、波阻抗的分层模型；预测待钻井钻头底下地层的测井曲线；结合井壁稳定力学模型预测该待钻井钻头底下地层的坍塌压力。公开号为CN1588127A的发明专利公开了一种利用地震层速度钻前预测坍塌压力的方法。无论是上述两篇专利的技术方案，还是其他现有技术方案，虽然都取得了一定的技术效果，但是仅考虑了钻井液密度较低时导致井壁失稳的状态，而忽略了钻井液密度较高时也会导致井壁失稳的问题，因此不能全面预测地层坍塌压力窗口，从而使现场操作受到限制。
发明内容
为解决现有技术中存在的问题，本发明提供一种预测弱面地层坍塌压力当量密度窗口的方法，其按照先后顺序包括以下步骤：
步骤一：根据弱面地层的特点，将其分为致密段和裂缝段；
步骤二：利用岩石三轴压缩试验测试岩石弹性力学参数、岩石本体强度、岩石弱面强度，利用声发射试验测试地层地应力分布，利用测井参数反演得到地层孔隙压力分布，利用气测渗透率仪测试地层渗透率，利用测井成像观测地层裂缝产状；
步骤三：利用流固耦合力学理论反演致密型地层井周围岩的应力分布和孔隙压力分布，并结合摩尔库伦破坏准则确定致密型地层坍塌压力当量密度的下限P致密型,下限；
步骤四：利用岩石弱面破坏准则分析弱面地层井周围岩的破坏状态，并确定裂缝型地层坍塌压力当量密度窗口的下限P裂缝型,下限和上限P裂缝型,上限；
步骤五：根据致密型地层和裂缝型地层的特点，建立弱面地层坍塌压力当量密度窗口的计算模型为PC＝{max(P致密型,下限,P裂缝型,下限),P裂缝型,上限}，并将测试得到的各参数值代入计算模型中，即可确定弱面地层坍塌压力当量密度窗口。
本发明的预测弱面地层坍塌压力当量密度窗口的方法，通过测井资料和室内试验确定模型计算所需要的岩石力学参数，综合致密型地层孔隙弹性模型和裂缝型地层弱面破坏模型，在给定某一允许的井径扩大率下，得到弱面地层坍塌压力当量密度窗口，并绘制井壁围岩失稳程度预测图示，以便在钻井设计和现场施工时为确定安全钻井液密度提供科学依据，可有效阻止井壁失稳，防止井下复杂情况的发生。
优选的是，所述步骤一中，致密段地层与裂缝段地层呈互层形式存在于弱面地层中。弱面地层具有强度各向异性的特点。
在上述任一方案中优选的是，所述步骤二中，岩石弹性力学参数包括弹性模量和泊松比。
在上述任一方案中优选的是，所述步骤二中，岩石本体强度包括岩石本体粘聚力和岩石本体内摩擦角。
在上述任一方案中优选的是，所述步骤二中，岩石弱面强度包括岩石弱面粘聚力和岩石弱面内摩擦角。
在上述任一方案中优选的是，所述步骤二中，地层地应力包括上覆岩层压力、水平最大主应力和水平最小主应力。
在上述任一方案中优选的是，所述步骤二中，地层裂缝产状包括裂缝倾向和裂缝倾角。
在上述任一方案中优选的是，所述步骤三中，利用流固耦合力学理论的基本方程和拉普拉斯变换原理解析井周围岩的孔隙压力分布和应力分布。
打开井壁瞬间，井周围岩应力重新分布，孔隙形状和体积发生变化，地层内部流体与固体之间发生强烈的相互作用，因此，在弱面地层的致密段采用流固耦合力学模型预测坍塌压力当量密度。基于毕奥特(Biot)流固耦合理论，可得到非均匀地应力场下的井周围岩应力分布和孔隙压力分布，结合摩尔库伦破坏准则，即可判断某一钻井液密度下井壁是否失稳。
在上述任一方案中优选的是，所述流固耦合力学理论的基本方程包括本构方程和其他方程。本构方程中，流体压力与致密段岩石基质变形是双向耦合在一起的。
在上述任一方案中优选的是，所述流固耦合力学理论的本构方程为
其中，p——井周围岩孔隙压力，MPa；
Mij——刚度矩阵系数；
σ——井周围岩应力，MPa；
ε——井周围岩应变；
α——平行于层里面的毕奥特系数；
α′——垂直于层里面的毕奥特系数。
其中，E——平行于层里面的弹性模量，GPa；
E′——垂直于层里面的弹性模量，GPa；
v——平行于层里面的泊松比；
v'——垂直于层里面的泊松比；
Ks——致密段岩石基质的体积模量，GPa。
在上述任一方案中优选的是，所述井周围岩孔隙压力的本构方程为
p＝M[ζ-α(εxx+εyy)-α'εzz]
其中，M——毕奥模量，GPa；
ζ——流体体积变化。
在上述任一方案中优选的是，所述流固耦合力学理论的其他方程包括流体运动方程、平衡方程、几何方程、质量守恒方程和协调方程。
在上述任一方案中优选的是，所述井周围岩应力分布的计算公式为
splap＝p0+S0[(cf/2Gκ)C1K2(ξ1r)+A1C2(rw 2/r2)]cos2(θ)
其中，Kn(x)——第n类修正贝塞尔函数，n为阶数；
s——复频域下的时间因子；
ζ——流体体积变化；
θ——井周角，°；
rw——井眼半径，cm；
r——井眼中心至地层内部某一点的距离，cm；
cf——流体扩散系数，m2/s；
α——平行于层里面的毕奥特系数；
α′——垂直于层里面的毕奥特系数；
v——平行于层里面的泊松比；
v'——垂直于层里面的泊松比；
G——平行于层里面的剪切模量，GPa；
G′——垂直于层里面的剪切模量，GPa；
κ——地层渗透率，Darcy(达西)；
pw——钻井液液柱压力，MPa；
Sv——上覆岩层压力，MPa；
SH——水平最大主应力，MPa；
Sh——水平最小主应力，MPa。
在上述任一方案中优选的是，所述井周围岩应力分布的计算公式中，
P0＝(Sh+SH)/2
A1＝αM/(M11+α2M)
A2＝M11+M12+2α2M/(M11+α2M)
B1＝(M11/2Gα)K2(ξ1rw)
B2＝[1/ξ1rw]K1(ξ1rw)+[6/(ξ1rw)2]K2(ξ1rw)
B2＝2{[1/ξ1rw]K1(ξ1rw)+[3/(ξ1rw)2]K2(ξ1rw)}
C1＝4/[2A1(B3-B2)-A2B1]
C2＝-4B1/[2A1(B3-B2)-A2B1]
C3＝[2A1(B2+B3)+3A2B1]/{3[2A1(B3-B2)-A2B1]}
其中，M——毕奥模量，GPa；
Mij——刚度矩阵系数。
通过流固耦合力学理论的基本方程相互代换，并结合拉普拉斯变换原理，即可得到拉普拉斯域下井周围岩应力分布的解析值。上述结果是在复频域下的解析值，后续还应通过拉普拉斯逆变换得到时域下井周围岩应力分布的解析值。σrr,σθθ,σzz,σrθ,σrz,σθz为拉普拉斯逆变换后的井周围岩应力分布。
在上述任一方案中优选的是，所述步骤三中，基于摩尔库伦破坏准则，建立井壁稳定力学模型为K≥0，给定某一允许的井径扩大率，即可确定致密型地层坍塌压力当量密度的下限P致密型,下限，其中，K——致密型地层坍塌压力指数，MPa。
在上述任一方案中优选的是，所述致密型地层坍塌压力指数为
其中，σ1——井壁上的最大主应力，MPa；
σ3——井壁上的最小主应力，MPa；
φ——岩石本体内摩擦角，°；
C——岩石本体粘聚力，MPa。
上述结果为时域下井周围岩应力分布的解析值。
在上述任一方案中优选的是，所述步骤四中，裂缝型地层坍塌压力当量密度窗口模型的建立，按照先后顺序包括以下步骤：
(1)将井周围岩应力分布从极坐标系下变换至井筒直角坐标系下；
(2)将井周围岩应力分布从井筒直角坐标系下变换至大地坐标系下；
(3)将井周围岩应力分布从大地坐标系下变换至弱面坐标系下；
(4)基于岩石弱面破坏准则，建立井壁稳定力学模型为N≥0，给定某一允许的井径扩大率，即可确定裂缝型地层坍塌压力当量密度窗口的下限P裂缝型,下限和上限P裂缝型,上限，其中，N——裂缝型地层坍塌压力指数，MPa。N<0，代表弱面地层发生剪切滑移破坏。
弱面地层内部发育不规则的天然裂缝，其力学特性表现为一组强度较低的弱面。地层内部存在两种破坏形式：一种沿着岩石本体破坏，另一种沿着裂缝剪切破坏。基于弹性力学坐标变换技术，将井周围岩应力分布转换至弱面坐标系下，通过判断弱面剪切力与正应力引起的摩擦力之间的大小关系，判断某一钻井液密度下弱面地层是否发生剪切破坏而导致井壁失稳。
在上述任一方案中优选的是，所述步骤(1)中，在井筒直角坐标系下，井周围岩应力分布为
其中，
在上述任一方案中优选的是，所述步骤(2)中，在大地坐标系下，井周围岩应力分布为
σGCS＝ET×σCCS×E
其中，
在上述任一方案中优选的是，所述步骤(3)中，在弱面坐标系下，井周围岩应力分布为
其中，
在上述任一方案中优选的是，所述步骤(4)中，裂缝型地层坍塌压力指数为
其中，φ——岩石弱面内摩擦角，°；
Sw——岩石弱面粘聚力，MPa；
θ——井周角，°；
pw——钻井液液柱压力，MPa；
αb——井筒方位角，°；
βb——井筒斜角，°；
——岩石弱面倾向，°；
βw——岩石弱面倾角，°。
本发明的预测弱面地层坍塌压力当量密度窗口的方法，简单易懂、操作便捷、成本低廉，可以在钻井前预测到弱面地层坍塌压力当量密度的窗口，以便在钻井设计和现场施工时为确定安全钻井液密度提供科学依据，进而有效阻止井壁失稳，防止井下复杂情况的发生。
附图说明
图1为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的一优选实施例流程图；
图2为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图1所示实施例中裂缝型致密砂岩地层测井成像图；
图3为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图1所示实施例中打开井壁瞬间最小主应力处的孔隙压力分布曲线；
图4为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图1所示实施例中打开井壁瞬间最小主应力处的坍塌压力指数分布曲线；
图5为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图1所示实施例中不同坐标系之间的变换关系图；
图6为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图1所示实施例中不同钻井液密度下裂缝型地层井周围岩坍塌程度对比图；
图7为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图1所示实施例中弱面地层坍塌压力当量密度窗口模型的分析图；
图8为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的另一优选实施例裂缝型致密页岩地层中打开井壁瞬间最小主应力处的坍塌压力指数分布曲线；
图9为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图8所示实施例中不同钻井液密度下裂缝型地层井周围岩坍塌程度对比图；
图10为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的另一优选实施例裂缝型致密碳酸盐地层中打开井壁瞬间最小主应力处的坍塌压力指数分布曲线；
图11为按照本发明预测弱面地层坍塌压力当量密度窗口的方法的图10所示实施例中不同钻井液密度下裂缝型地层井周围岩坍塌程度对比图。
具体实施方式
为了更进一步了解本发明的发明内容，下面将结合具体实施例详细阐述本发明。
实施例一：
本实施例结合某区块裂缝型致密砂岩地层的实际情况，预测弱面地层坍塌压力当量密度窗口。该地层厚度为100-400m，为一套红色碎屑岩，主要出露于某坳陷北部。该地层分为三岩性段，自上而下依次为一段、二段和三段。一段厚度为20-70m，以含砾砂岩、砂岩为主，裂缝不发育，为致密型砂岩段；二段厚度为40-260m，以细砂岩夹薄层泥岩为主，有相对质纯的泥岩薄层出现；三段厚度为45-100m，主要为褐红色砾岩、砂砾岩夹含砾砂岩、粉砂岩。二段和三段裂缝发育，为裂缝型砂岩段，裂缝宽度为0.1-2.0mm。该地层的砂岩基质总体属于低孔低渗-特低孔特低渗储层，渗透率在100纳米达西至100毫米达西之间。通常情况使用一套钻井液密度钻穿储层，因此钻井液密度应同时满足裂缝型和致密型砂岩地层的井壁稳定要求。
如图1所示，一种预测弱面地层坍塌压力当量密度窗口的方法，其按照先后顺序包括以下步骤：
步骤一：根据弱面地层的特点，将其分为致密段和裂缝段；
步骤二：利用岩石三轴压缩试验测试岩石弹性力学参数、岩石本体强度、岩石弱面强度，利用声发射试验测试地层地应力分布，利用测井参数反演得到地层孔隙压力分布，利用气测渗透率仪测试地层渗透率，利用测井成像观测地层裂缝产状；
步骤三：利用流固耦合力学理论反演致密型地层井周围岩的应力分布和孔隙压力分布，并结合摩尔库伦破坏准则确定致密型地层坍塌压力当量密度的下限P致密型,下限；
步骤四：利用岩石弱面破坏准则分析弱面地层井周围岩的破坏状态，并确定裂缝型地层坍塌压力当量密度窗口的下限P裂缝型,下限和上限P裂缝型,上限；
步骤五：根据致密型地层和裂缝型地层的特点，建立弱面地层坍塌压力当量密度窗口的计算模型为PC＝{max(P致密型,下限,P裂缝型,下限),P裂缝型,上限}，并将测试得到的各参数值代入计算模型中，即可确定弱面地层坍塌压力当量密度窗口。
步骤一中，致密段地层与裂缝段地层呈互层形式存在于弱面地层中。
步骤二中，岩石弹性力学参数包括弹性模量和泊松比；岩石本体强度包括岩石本体粘聚力和岩石本体内摩擦角；岩石弱面强度包括岩石弱面粘聚力和岩石弱面内摩擦角；地层地应力包括上覆岩层压力、水平最大主应力和水平最小主应力；地层裂缝产状包括裂缝倾向和裂缝倾角。本实施例的裂缝型致密砂岩地层测井成像如图2所示。砂岩弱面地层各参数值如表1.1所示。
表1.1砂岩弱面地层各参数值
步骤三中，利用流固耦合力学理论的基本方程和拉普拉斯变换原理解析井周围岩的孔隙压力分布和应力分布。流固耦合力学理论的基本方程包括本构方程和其他方程。流固耦合力学理论的其他方程包括流体运动方程、平衡方程、几何方程、质量守恒方程和协调方程。本构方程中，流体压力与致密段岩石基质变形是双向耦合在一起的。
流固耦合力学理论的本构方程为：
其中，p——井周围岩孔隙压力，MPa；
Mij——刚度矩阵系数；
σ——井周围岩应力，MPa；
ε——井周围岩应变；
α——平行于层里面的毕奥特系数；
α′——垂直于层里面的毕奥特系数。
其中，E——平行于层里面的弹性模量，GPa；
E′——垂直于层里面的弹性模量，GPa；
v——平行于层里面的泊松比；
v'——垂直于层里面的泊松比；
Ks——致密段岩石基质的体积模量，GPa。
井周围岩孔隙压力的本构方程为：
p＝M[ζ-α(εxx+εyy)-α'εzz]
其中，M——毕奥模量，GPa；
ζ——流体体积变化。
井周围岩应力分布的计算公式为：
splap＝p0+S0[(cf/2Gκ)C1K2(ξ1r)+A1C2(rw 2/r2)]cos2(θ)
其中，Kn(x)——第n类修正贝塞尔函数，n为阶数；
s——复频域下的时间因子；
ζ——流体体积变化；
θ——井周角，°；
rw——井眼半径，cm；
r——井眼中心至地层内部某一点的距离，cm；
cf——流体扩散系数，m2/s；
α——平行于层里面的毕奥特系数；
α′——垂直于层里面的毕奥特系数；
v——平行于层里面的泊松比；
v'——垂直于层里面的泊松比；
G——平行于层里面的剪切模量，GPa；
G′——垂直于层里面的剪切模量，GPa；
κ——地层渗透率，Darcy(达西)；
pw——钻井液液柱压力，MPa；
Sv——上覆岩层压力，MPa；
SH——水平最大主应力，MPa；
Sh——水平最小主应力，MPa。
井周围岩应力分布的计算公式中：
P0＝(Sh+SH)/2
A1＝αM/(M11+α2M)
A2＝M11+M12+2α2M/(M11+α2M)
B1＝(M11/2Gα)K2(ξ1rw)
B2＝[1/ξ1rw]K1(ξ1rw)+[6/(ξ1rw)2]K2(ξ1rw)
B2＝2{[1/ξ1rw]K1(ξ1rw)+[3/(ξ1rw)2]K2(ξ1rw)}
C1＝4/[2A1(B3-B2)-A2B1]
C2＝-4B1/[2A1(B3-B2)-A2B1]
C3＝[2A1(B2+B3)+3A2B1]/{3[2A1(B3-B2)-A2B1]}
其中，M——毕奥模量，GPa；
Mij——刚度矩阵系数。
通过流固耦合力学理论的基本方程相互代换，并结合拉普拉斯变换原理，即可得到拉普拉斯域下井周围岩应力分布的解析值。上述结果是在复频域下的解析值，后续还应通过拉普拉斯逆变换得到时域下井周围岩应力分布的解析值。σrr,σθθ,σzz,σrθ,σrz,σθz为拉普拉斯逆变换后的井周围岩应力分布。
基于毕奥特(Biot)流固耦合理论，可得到非均匀地应力场下井周围岩的应力分布和孔隙压力分布，结合摩尔库伦破坏准则，即可判断某一钻井液密度下井壁是否失稳，进而确定致密型地层坍塌压力当量密度的下限。
打开井壁后不同时间时，最小主应力处的孔隙压力分布情况如图3所示。不渗透井壁(泥饼渗透率很低)与渗透型井壁均会在围岩内产生异常高的孔隙压力。采用传统模型预测时，不考虑流固耦合作用，因此孔隙压力分布为原始地层压力。
基于摩尔库伦破坏准则，建立井壁稳定力学模型为K≥0，给定某一允许的井径扩大率，即可确定致密型地层坍塌压力当量密度的下限P致密型,下限，其中，K——致密型地层坍塌压力指数，MPa。
致密型地层坍塌压力指数为：
其中，σ1——井壁上的最大主应力，MPa；
σ3——井壁上的最小主应力，MPa；
φ——岩石本体内摩擦角，°；
C——岩石本体粘聚力，MPa。
上述结果为时域下井周围岩应力分布的解析值。
不同钻井液密度下，打开井壁瞬间坍塌压力指数分布情况如图4所示。当钻井液密度提高时，曲线逐渐上移，有利于井壁稳定。若给定一允许井径扩大率，则可确定致密型地层坍塌压力当量密度的下限。
步骤四中，裂缝型地层坍塌压力当量密度窗口模型的建立，按照先后顺序包括以下步骤：
(1)将井周围岩应力分布从极坐标系下变换至井筒直角坐标系下。
在井筒直角坐标系下，井周围岩应力分布为：
其中，
(2)将井周围岩应力分布从井筒直角坐标系下变换至大地坐标系下。
在大地坐标系下，井周围岩应力分布为：
σGCS＝ET×σCCS×E
其中，
(3)将井周围岩应力分布从大地坐标系下变换至弱面坐标系下。
在弱面坐标系下，井周围岩应力分布为：
其中，
(4)基于岩石弱面破坏准则，建立井壁稳定力学模型为N≥0，给定某一允许的井径扩大率，即可确定裂缝型地层坍塌压力当量密度窗口的下限P裂缝型,下限和上限P裂缝型,上限，其中，N——裂缝型地层坍塌压力指数，MPa。N<0，代表弱面地层发生剪切滑移破坏。
裂缝型地层坍塌压力指数为：
其中，φ——岩石弱面内摩擦角，°；
Sw——岩石弱面粘聚力，MPa；
θ——井周角，°；
pw——钻井液液柱压力，MPa；
αb——井筒方位角，°；
βb——井筒斜角，°；
——岩石弱面倾向，°；
βw——岩石弱面倾角，°。
基于弹性力学坐标变换技术，将井周围岩应力分布转换至弱面坐标系下，如图5所示。通过判断弱面剪切力与正应力引起的摩擦力之间的大小关系，判断某一钻井液密度下弱面地层是否发生剪切破坏而导致井壁失稳。
图6为不同钻井液密度下井周围岩坍塌程度对比图，该图表示井周围岩屈服区域随钻井液密度增加而发生变化。当钻井液密度较低时，坍塌区域在东偏南10度至东偏南90度和西偏北10度至西偏北90度的范围内。随钻井液密度的增大，之前的屈服区域变小，但是在正北和正南的位置，新的屈服区域诞生并不断扩大。结果显示钻井液密度过高或过低均导致井壁失稳，但是屈服区域分布于井周不同区域或位置。当给定某一允许井径扩大率时，比如20％，即可确定裂缝型地层坍塌压力当量密度窗口的下限和上限，P裂缝型,下限＝1.5g/cm3，P裂缝型,上限＝1.9g/cm3。1.7g/cm3为最优选的钻井液密度。
将表1.1中各参数值代入弱面地层坍塌压力当量密度窗口模型PC＝{max(P致密型,下限,P裂缝型,下限),P裂缝型,上限}中，经过计算，并由图6可知，若给定允许井径扩大率为20％，则致密型地层坍塌压力当量密度的下限为1.8g/cm3，裂缝型地层坍塌压力当量密度窗口的下限为1.5g/cm3、上限为1.9g/cm3。因此，本实施例的弱面地层坍塌压力当量密度窗口为1.8-1.9g/cm3，现场选用的钻井液密度在此窗口内，可确保井径扩大率小于等于20％。
如图7所示，利用本实施例的预测方法分析钻进施工过程中的井壁失稳现象：某井6770m至6856m为裂缝型致密地层，发育一组天然裂缝。初期钻进时，钻井液密度为1.97g/cm3，后期增大钻井液密度至2.1g/cm3后，井壁反而垮塌更严重。测井资料包括四臂测井资料，其结果分析如下：其中C13臂走向为正东、C24臂走向为正北。6770-6833m时，钻井液密度为1.97g/cm3，C24臂实测井径与模拟结果符合良好，而C13臂实测井径大于模拟结果，这是因为6833m后采用密度为2.1g/cm3的钻井液钻进，增加的钻井液密度使6770-6833m在C13臂方向的井径扩大；6833-6856m时，钻井液密度为2.1g/cm3，测井资料所得C13臂方向的井径大于C24臂方向，且C24臂方向的井径相比于6770-6833m时的井径减小，与模拟结果相符。利用本实施例的预测方法分析可知，所得结果与现场实际井壁垮塌情况符合良好，该方法具有较高的准确性和可靠性。
实施例二：
本实施例结合某区块裂缝型致密页岩地层的实际情况，预测弱面地层坍塌压力当量密度窗口。其预测方法、原理、有益效果等与实施例一相同，不同的是：页岩弱面地层各参数值如表2.1所示。
表2.1页岩弱面地层各参数值
将表2.1中各参数值代入模型进行计算，计算结果如图8和图9所示：若给定允许井径扩大率为13％，则致密型页岩段坍塌压力当量密度的下限为1.1g/cm3；裂缝型页岩段坍塌压力当量密度的下限为0.95g/cm3、上限为1.35g/cm3。因此，本地层坍塌压力当量密度窗口为1.1-1.35g/cm3，现场实际选用钻井液密度在此窗口内时，可确保井径扩大率小于等于13％。
实施例三：
本实施例结合某区块裂缝型致密碳酸盐地层的实际情况，预测弱面地层坍塌压力当量密度窗口。其预测方法、原理、有益效果等与实施例一相同，不同的是：页岩弱面地层各参数值如表3.1所示。
表3.1碳酸盐弱面地层各参数值
将表3.1中各参数值代入模型进行计算，计算结果如图10和图11所示：若给定允许井径扩大率为12％，则致密型碳酸盐地层坍塌压力当量密度的下限为1.45g/cm3；裂缝型碳酸盐地层坍塌压力当量密度窗口的下限为1.5g/cm3、上限为1.6g/cm3。因此，本地层坍塌压力当量密度窗口为1.5-1.6g/cm3，现场实际选用的钻井液密度在此窗口内时，可确保井径扩大率小于等于12％。
本领域技术人员不难理解，本发明的预测弱面地层坍塌压力当量密度窗口的方法包括上述本发明说明书的发明内容和具体实施方式部分以及附图所示出的各部分的任意组合，限于篇幅并为使说明书简明而没有将这些组合构成的各方案一一描述。凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (20)

1.一种预测弱面地层坍塌压力当量密度窗口的方法，其按照先后顺序包括以下步骤：
步骤一：根据弱面地层的特点，将其分为致密段和裂缝段；
步骤二：利用岩石三轴压缩试验测试岩石弹性力学参数、岩石本体强度、岩石弱面强度，利用声发射试验测试地层地应力分布，利用测井参数反演得到地层孔隙压力分布，利用气测渗透率仪测试地层渗透率，利用测井成像观测地层裂缝产状；
步骤三：利用流固耦合力学理论反演致密型地层井周围岩的应力分布和孔隙压力分布，并结合摩尔库伦破坏准则确定致密型地层坍塌压力当量密度的下限P致密型,下限；
步骤四：利用岩石弱面破坏准则分析弱面地层井周围岩的破坏状态，并确定裂缝型地层坍塌压力当量密度窗口的下限P裂缝型,下限和上限P裂缝型,上限；
步骤五：根据致密型地层和裂缝型地层的特点，建立弱面地层坍塌压力当量密度窗口的计算模型为PC＝{max(P致密型,下限,P裂缝型,下限),P裂缝型,上限}，并将测试得到的各参数值代入计算模型中，即可确定弱面地层坍塌压力当量密度窗口。
2.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤一中，致密段地层与裂缝段地层呈互层形式存在于弱面地层中。
3.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤二中，岩石弹性力学参数包括弹性模量和泊松比。
4.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤二中，岩石本体强度包括岩石本体粘聚力和岩石本体内摩擦角。
5.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤二中，岩石弱面强度包括岩石弱面粘聚力和岩石弱面内摩擦角。
6.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤二中，地层地应力包括上覆岩层压力、水平最大主应力和水平最小主应力。
7.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤二中，地层裂缝产状包括裂缝倾向和裂缝倾角。
8.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤三中，利用流固耦合力学理论的基本方程和拉普拉斯变换原理解析井周围岩的孔隙压力分布和应力分布。
9.如权利要求8所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述流固耦合力学理论的基本方程包括本构方程和其他方程；所述其他方程包括流体运动方程、平衡方程、几何方程、质量守恒方程和协调方程。
10.如权利要求9所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述流固耦合力学理论的本构方程为
<mrow> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>z</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>y</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>z</mi> <mi>x</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>M</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>12</mn> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>13</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>M</mi> <mn>12</mn> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>13</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>M</mi> <mn>13</mn> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>13</mn> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>33</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>M</mi> <mn>44</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>M</mi> <mn>55</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>M</mi> <mn>55</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>z</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>y</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>z</mi> <mi>x</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>&amp;alpha;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&amp;alpha;</mi> </mtd> </mtr> <mtr> <mtd> <msup> <mi>&amp;alpha;</mi> <mo>&amp;prime;</mo> </msup> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>p</mi> </mrow>
其中，p——井周围岩孔隙压力，MPa；
Mij——刚度矩阵系数；
σ——井周围岩应力，MPa；
ε——井周围岩应变；
α——平行于层里面的毕奥特系数；
α′——垂直于层里面的毕奥特系数；
γ——切应变；
τ——切应力，MPa。
11.如权利要求10所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述井周围岩孔隙压力的本构方程为
p＝M[ζ-α(εxx+εyy)-α'εzz]
其中，M——毕奥模量，GPa；
ζ——流体体积变化。
12.如权利要求8所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述井周围岩应力分布的计算公式为
splap＝p0+S0[(cf/2Gκ)C1K2(ξ1r)+A1C2(rw 2/r2)]cos2(θ)
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>s&amp;sigma;</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <mi>cos</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>w</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <mo>&amp;lsqb;</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>{</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mo>+</mo> <mo>&amp;lsqb;</mo> <mn>6</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> <msub> <mi>K</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>-</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <msub> <mi>C</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mn>3</mn> <msub> <mi>C</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>4</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>4</mn> </msup> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mi>cos</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>s&amp;sigma;</mi> <mrow> <mi>&amp;theta;</mi> <mi>&amp;theta;</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <mi>cos</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>w</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <mo>&amp;lsqb;</mo> <mo>-</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>{</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mo>+</mo> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>+</mo> <mn>6</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> <msub> <mi>K</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <mn>3</mn> <msub> <mi>C</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>4</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>4</mn> </msup> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mi>cos</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <msubsup> <mi>s&amp;sigma;</mi> <mrow> <mi>z</mi> <mi>z</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <msub> <mi>S</mi> <mi>v</mi> </msub> <mo>+</mo> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>h</mi> </msub> <mo>+</mo> <msub> <mi>S</mi> <mi>H</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msup> <mi>&amp;alpha;</mi> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <mn>2</mn> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> <msub> <mi>p</mi> <mn>0</mn> </msub> <mo>+</mo> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msubsup> <mi>s&amp;sigma;</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>s&amp;sigma;</mi> <mrow> <mi>&amp;theta;</mi> <mi>&amp;theta;</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <msup> <mi>&amp;alpha;</mi> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <mn>2</mn> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> <msup> <mi>sp</mi> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msup> </mrow>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>s&amp;sigma;</mi> <mrow> <mi>r</mi> <mi>&amp;theta;</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <mo>&amp;lsqb;</mo> <mn>2</mn> <msub> <mi>A</mi> <mn>1</mn> </msub> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>{</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mo>+</mo> <mo>&amp;lsqb;</mo> <mn>3</mn> <mo>/</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> <msub> <mi>K</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <mi>r</mi> <mo>)</mo> </mrow> <mo>}</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> <msub> <mi>C</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mn>3</mn> <msub> <mi>C</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>r</mi> <mi>w</mi> </msub> <mn>4</mn> </msup> <mo>/</mo> <msup> <mi>r</mi> <mn>4</mn> </msup> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mi>sin</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>r</mi> <mi>z</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>&amp;theta;</mi> <mi>z</mi> </mrow> <mrow> <mi>l</mi> <mi>a</mi> <mi>p</mi> </mrow> </msubsup> <mo>=</mo> <mn>0</mn> </mrow>
其中，Kn(x)——第n类修正贝塞尔函数，n为阶数；
s——复频域下的时间因子；
ζ——流体体积变化；
θ——井周角，°；
rw——井眼半径，cm；
r——井眼中心至地层内部某一点的距离，cm；
cf——流体扩散系数，m2/s；
α——平行于层里面的毕奥特系数；
α′——垂直于层里面的毕奥特系数；
v'——垂直于层里面的泊松比；
G——平行于层里面的剪切模量，GPa；
κ——地层渗透率，Darcy；
pw——钻井液液柱压力，MPa；
Sv——上覆岩层压力，MPa；
SH——水平最大主应力，MPa；
Sh——水平最小主应力，MPa；
P0、S0、C、A——公式中的换算参数；
plap——拉普拉斯域下的井壁围岩孔隙压力，MPa；
σlap——拉普拉斯域下的井壁围岩应力，MPa。
13.如权利要求12所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述井周围岩应力分布的计算公式中，
P0＝(Sh+SH)/2
<mrow> <msub> <mi>S</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>0.5</mn> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>H</mi> </msub> <mo>-</mo> <msub> <mi>S</mi> <mi>h</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow>
<mrow> <msub> <mi>&amp;zeta;</mi> <mn>1</mn> </msub> <mo>=</mo> <msqrt> <mrow> <mi>s</mi> <mo>/</mo> <msub> <mi>c</mi> <mi>f</mi> </msub> </mrow> </msqrt> </mrow>
A1＝αM/(M11+α2M)
A2＝M11+M12+2α2M/(M11+α2M)
B1＝(M11/2Gα)K2(ξ1rw)
B2＝[1/ξ1rw]K1(ξ1rw)+[6/(ξ1rw)2]K2(ξ1rw)
B3＝2{[1/ξ1rw]K1(ξ1rw)+[3/(ξ1rw)2]K2(ξ1rw)}
C1＝4/[2A1(B3-B2)-A2B1]
C2＝-4B1/[2A1(B3-B2)-A2B1]
C3＝[2A1(B2+B3)+3A2B1]/{3[2A1(B3-B2)-A2B1]}
其中，M——毕奥模量，GPa；
Mij——刚度矩阵系数。
14.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤三中，基于摩尔库伦破坏准则，建立井壁稳定力学模型为K≥0，给定某一允许的井径扩大率，即可确定致密型地层坍塌压力当量密度的下限P致密型,下限，其中，K——致密型地层坍塌压力指数，MPa。
15.如权利要求14所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述致密型地层坍塌压力指数为
其中，σ1——井壁上的最大主应力，MPa；
σ3——井壁上的最小主应力，MPa；
φ——岩石本体内摩擦角，°；
C——岩石本体粘聚力，MPa；
p——地层孔隙压力，MPa。
16.如权利要求1所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤四中，裂缝型地层坍塌压力当量密度窗口模型的建立，按照先后顺序包括以下步骤：
步骤一：将井周围岩应力分布从极坐标系下变换至井筒直角坐标系下；
步骤二：将井周围岩应力分布从井筒直角坐标系下变换至大地坐标系下；
步骤三：将井周围岩应力分布从大地坐标系下变换至弱面坐标系下；
步骤四：基于岩石弱面破坏准则，建立井壁稳定力学模型为N≥0，给定某一允许的井径扩大率，即可确定裂缝型地层坍塌压力当量密度窗口的下限P裂缝型,下限和上限P裂缝型,上限，其中，N——裂缝型地层坍塌压力指数，MPa。
17.如权利要求16所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤一中，在井筒直角坐标系下，井周围岩应力分布为
<mrow> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>C</mi> <mi>C</mi> <mi>S</mi> </mrow> </msub> <mo>=</mo> <msup> <mi>C</mi> <mi>T</mi> </msup> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>r</mi> <mi>&amp;theta;</mi> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>r</mi> <mi>&amp;theta;</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>&amp;theta;</mi> <mi>&amp;theta;</mi> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>z</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>C</mi> <mo>=</mo> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msub> <mi>&amp;sigma;</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;sigma;</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>y</mi> <mi>z</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>z</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;tau;</mi> <mrow> <mi>y</mi> <mi>z</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>&amp;sigma;</mi> <mi>z</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中，
θ——井周角，°；
σ——正应力，MPa；
τ——切应力，MPa。
18.如权利要求17所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤二中，在大地坐标系下，井周围岩应力分布为
σGCS＝ET×σCCS×E
其中，
αb——方位角，°；
βb——井斜角，°。
19.如权利要求18所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤三中，在弱面坐标系下，井周围岩应力分布为
<mrow> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>B</mi> <mi>P</mi> <mi>C</mi> <mi>S</mi> </mrow> </msub> <mo>=</mo> <mi>W</mi> <mo>&amp;times;</mo> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>G</mi> <mi>C</mi> <mi>S</mi> </mrow> </msub> <mo>&amp;times;</mo> <msup> <mi>W</mi> <mi>T</mi> </msup> <mo>=</mo> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> <mi>w</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> <mi>w</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>x</mi> <mi>z</mi> </mrow> <mi>w</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>y</mi> <mi>x</mi> </mrow> <mi>w</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> <mi>w</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>y</mi> <mi>z</mi> </mrow> <mi>w</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>z</mi> <mi>x</mi> </mrow> <mi>w</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>z</mi> <mi>y</mi> </mrow> <mi>w</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>z</mi> <mi>z</mi> </mrow> <mi>w</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中，
σw——弱面上的正应力，MPa；
τw——弱面上的切应力，MPa；
αw——弱面倾向减去
βw——弱面倾角，°。
20.如权利要求19所述的预测弱面地层坍塌压力当量密度窗口的方法，其特征在于：所述步骤四中，裂缝型地层坍塌压力指数为
<mrow> <mi>N</mi> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>w</mi> </msub> <mo>+</mo> <mi>t</mi> <mi>g</mi> <mo>(</mo> <mi>&amp;phi;</mi> <mo>)</mo> <mo>(</mo> <mrow> <msubsup> <mi>&amp;sigma;</mi> <mrow> <mi>z</mi> <mi>z</mi> </mrow> <mi>w</mi> </msubsup> <mo>-</mo> <msub> <mi>p</mi> <mi>w</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>z</mi> <mi>x</mi> </mrow> <mi>w</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>&amp;tau;</mi> <mrow> <mi>z</mi> <mi>y</mi> </mrow> <mi>w</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
其中，φ——岩石弱面内摩擦角，°；
Sw——岩石弱面粘聚力，MPa；
pw——井底压力，MPa。
CN201510077224.7A 2015-02-12 2015-02-12 一种预测弱面地层坍塌压力当量密度窗口的方法 Active CN104806233B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510077224.7A CN104806233B (zh) 2015-02-12 2015-02-12 一种预测弱面地层坍塌压力当量密度窗口的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510077224.7A CN104806233B (zh) 2015-02-12 2015-02-12 一种预测弱面地层坍塌压力当量密度窗口的方法

Publications (2)

Publication Number Publication Date
CN104806233A CN104806233A (zh) 2015-07-29
CN104806233B true CN104806233B (zh) 2018-02-23

Family

ID=53691385

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510077224.7A Active CN104806233B (zh) 2015-02-12 2015-02-12 一种预测弱面地层坍塌压力当量密度窗口的方法

Country Status (1)

Country Link
CN (1) CN104806233B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105298478A (zh) * 2015-09-08 2016-02-03 中国石油大学(北京) 一种用于断裂构造的地层孔隙压力确定方法
CN107120107B (zh) * 2016-02-24 2020-11-13 中国石油化工股份有限公司 海底钻井的钻井液选择方法和其在钻井深度计算中的用途
CN105735980A (zh) * 2016-03-17 2016-07-06 成都创源油气技术开发有限公司 一种气体钻井井壁稳定性评价方法
CN108254262B (zh) * 2016-12-28 2021-03-30 中国石油天然气股份有限公司 岩石层理裂缝剪切参数预测方法及装置
EP3759314A1 (en) * 2018-02-27 2021-01-06 Saudi Arabian Oil Company Determining a mudweight of drilling fluids for drilling through naturally fractured formations
CN108843313B (zh) * 2018-06-23 2019-07-12 东北石油大学 页岩地层钻井安全钻井液密度窗口设计方法
CN110792418B (zh) * 2018-08-03 2022-03-01 中国石油天然气股份有限公司 井筒工作液配方优化方法及装置
CN109555515A (zh) * 2018-10-08 2019-04-02 中国石油大学(北京) 地层坍塌压力确定方法和装置
CN111155962B (zh) * 2019-01-15 2022-02-08 王瀚艺 用于提高套管井筒完整性的系统和方法
CN110145286B (zh) * 2019-05-22 2021-10-01 西南石油大学 一种低渗透油藏或气藏完井工程的设计方法
CN110826142B (zh) * 2019-12-06 2020-06-05 西南石油大学 一种裂缝性地层封堵承压能力的预测方法
CN111322064A (zh) * 2020-03-13 2020-06-23 西南石油大学 一种将井壁崩落法应用于页岩储层反演地应力的校正方法
CN113343336B (zh) * 2021-05-31 2022-03-11 西南石油大学 一种井壁坍塌渐进破坏过程的数值模拟方法
CN113935092B (zh) * 2021-10-23 2022-08-23 西南石油大学 一种裂缝性地层井壁稳定双重孔隙热弹性耦合评价方法
CN114564778B (zh) * 2022-03-07 2022-11-18 西南石油大学 一种裂缝性地层计算安全钻井密度的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5472050A (en) * 1994-09-13 1995-12-05 Union Oil Company Of California Use of sequential fracturing and controlled release of pressure to enhance production of oil from low permeability formations
CN1588128A (zh) * 2004-10-22 2005-03-02 石油大学(北京) 利用地震记录预测井眼待钻井段坍塌压力和破裂压力的方法
CN1588127A (zh) * 2004-10-22 2005-03-02 石油大学(北京) 利用地震层速度钻前预测坍塌压力与破裂压力的方法
CN1966934A (zh) * 2005-11-16 2007-05-23 中国石油大学(北京) 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法
CN102182453A (zh) * 2011-03-03 2011-09-14 中国石油集团钻井工程技术研究院 井壁坍塌分析方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5472050A (en) * 1994-09-13 1995-12-05 Union Oil Company Of California Use of sequential fracturing and controlled release of pressure to enhance production of oil from low permeability formations
CN1588128A (zh) * 2004-10-22 2005-03-02 石油大学(北京) 利用地震记录预测井眼待钻井段坍塌压力和破裂压力的方法
CN1588127A (zh) * 2004-10-22 2005-03-02 石油大学(北京) 利用地震层速度钻前预测坍塌压力与破裂压力的方法
CN1966934A (zh) * 2005-11-16 2007-05-23 中国石油大学(北京) 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法
CN102182453A (zh) * 2011-03-03 2011-09-14 中国石油集团钻井工程技术研究院 井壁坍塌分析方法

Also Published As

Publication number Publication date
CN104806233A (zh) 2015-07-29

Similar Documents

Publication Publication Date Title
CN104806233B (zh) 一种预测弱面地层坍塌压力当量密度窗口的方法
Bell Practical methods for estimating in situ stresses for borehole stability applications in sedimentary basins
CA2377467C (en) Cavity stability prediction method for wellbores
Hou et al. Hydraulic fracture initiation theory for a horizontal well in a coal seam
CN109356567B (zh) 深水浅部地层井壁稳定性预测方法
Guo et al. Evaluation of fracability and screening of perforation interval for tight sandstone gas reservoir in western Sichuan Basin
CN110580401B (zh) 判断定向井、水平井分段多簇压裂井筒暂堵次数的方法
Archer et al. A log based analysis to estimate mechanical properties and in-situ stresses in a shale gas well in North Perth Basin
CN103278614A (zh) 一种动静态岩石力学参数校正方法及装置
CN103257081A (zh) 一种油气藏岩体力学地下原位模型恢复的方法及装置
van-der-Zee et al. 3D geomechanical modeling of complex salt structures
WO2022193655A1 (zh) 一种基于井壁力学失稳垮塌的地应力大小评价方法
Yu et al. On how asymmetric stimulated rock volume in shales may impact casing integrity
CN108678672A (zh) 深水浅部复杂岩性地层导管高效置入及井口稳定预测方法
Hou et al. Prediction of wellbore stability in conglomerate formation using discrete element method
CN114547906A (zh) 一种深部含软弱结构面地层的井壁稳定测井解释方法
Fekete et al. Wellbore stability management in weak bedding planes and angle of attack in well planing
WO2013138584A1 (en) Screening potential geomechanical risks during waterflooding
CN108708715A (zh) 快速预测页岩气层地层破裂压力梯度的测录井方法
Moos et al. Impact of rock properties on the relationship between wellbore breakout width and depth
Al-Malikee et al. Indirect prediction of rock elasticity and compressibility strength using well log data at selected sites within Rumaila Oilfield, Southern Iraq
CN112412434B (zh) 一种改进的疏松砂岩地应力计算方法
Zhang et al. 3D geomechanical evaluation of jetting operations for producing hydrate chimneys in the sea of Japan
Al Shehhi et al. In-Situ Tectonic Stress Measurement from Wireline Straddle Packer Microfrac Testing: A Case Study in a Carbonate Reservoir (UAE)
CN109033698B (zh) 一种用于层状地层水平井破裂压力计算的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant