CN114564778A - 一种裂缝性地层计算安全钻井密度的方法 - Google Patents

一种裂缝性地层计算安全钻井密度的方法 Download PDF

Info

Publication number
CN114564778A
CN114564778A CN202210216104.0A CN202210216104A CN114564778A CN 114564778 A CN114564778 A CN 114564778A CN 202210216104 A CN202210216104 A CN 202210216104A CN 114564778 A CN114564778 A CN 114564778A
Authority
CN
China
Prior art keywords
density
drilling fluid
model
drilling
borehole
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.)
Granted
Application number
CN202210216104.0A
Other languages
English (en)
Other versions
CN114564778B (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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202210216104.0A priority Critical patent/CN114564778B/zh
Publication of CN114564778A publication Critical patent/CN114564778A/zh
Application granted granted Critical
Publication of CN114564778B publication Critical patent/CN114564778B/zh
Priority to US18/118,229 priority patent/US11822046B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B21/00—Methods or apparatus for flushing boreholes, e.g. by use of exhaust air from motor
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/002—Survey of boreholes or wells by visual inspection
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
    • E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Geophysics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Structural Engineering (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Mechanical Engineering (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Earth Drilling (AREA)

Abstract

本发明公开了一种裂缝性地层计算安全钻井密度的方法，包括以下步骤：S1、图像处理，识别井下裂缝；S2、通过井下裂缝参数建立三维地质模型，利用井眼尺寸和长度建立钻井井眼模型；S3、赋值模型材料参数、边界条件、初始钻井液密度上下界，然后计算精度；S4、利用3DEC离散元软件，求解三维地质模型，判定井壁稳定性；S5、利用二分法确定钻井液密度上下界；S6、重复执行步骤S4~S5；S7、达到设定的精度条件后，保存并输出安全钻井液密度。本发明以成像测井图为基础，建立三维井眼模型，使用离散单元理论，计算安全钻井液密度，可在破碎带、断层、裂缝性地层的钻井液密度计算更加准确，更好的指导裂缝性地层的钻井工程钻进。

Description

一种裂缝性地层计算安全钻井密度的方法
技术领域
本发明涉及石油与天然气工程中的钻井实践、地热钻探、深地科学探索领域，尤其涉及一种裂缝性地层计算安全钻井密度的方法。
背景技术
随着石油与天然气工程的发展，钻井深度越来越深。钻井遇到的地层，往往不是完整的岩石，而是包含结构面的岩体。在钻遇破碎带、断层、裂缝时，这种问题现象尤为突出。钻井井眼需要钻井液压力的支撑，不然就会出现井壁缩颈、坍塌等问题。进而造成卡钻，埋钻杆等问题，严重情况时将报废整个井眼。
因此，钻井所需的安全钻井液压力或者安全钻井液的密度对于钻井工程的成败起着决定性的作用。一般来说，目前计算钻井液密度是利用连续性介质理论，开展井壁应力计算，通过研究不同钻井液密度对井壁稳定的影响程度，利用二分法后者牛顿迭代法计算得到所需钻井液密度。常见的计算流程如下：
1.利用Farihurst方程计算井壁应力；
2.使用Mohr-Coulomb强度准则判定井壁稳定性；
3.利用二分法或者牛顿迭代法计算钻井液密度。
在钻遇破碎带、断层、裂缝等地层时，由于基本假设和理论上的不同，无法考虑地层的几何产状等因素对井壁稳定和钻井液密度的作用。因此，上述方法不完全适用于破裂带、断层、裂缝性地层钻井工程的钻井液密度计算。
发明内容
本发明的目的在于克服现有技术的不足，提供一种裂缝性地层计算安全钻井密度的方法，该方法在破碎带、断层、裂缝性地层的钻井液密度计算更加准确，更好的指导裂缝性地层的钻井工程钻进。
本发明的目的是通过以下技术方案来实现的：
一种裂缝性地层计算安全钻井密度的方法，包括以下步骤：
S1、图像处理，识别井下裂缝；
S2、通过井下裂缝参数建立三维地质模型，利用井眼尺寸和长度建立钻井井眼模型；
S3、赋值模型材料参数、边界条件、初始钻井液密度上下界，然后计算精度；
S4、利用3DEC离散元软件，求解三维地质模型，判定井壁稳定性；
S5、利用二分法确定钻井液密度上下界；
S6、重复执行步骤S4～S5；
S7、达到设定的精度条件后，保存并输出安全钻井液密度。
进一步的，所述S1具体为先将地层成像测井图转换为灰度图，同时得到灰度颜色直方图，利用阈值方法等分割算法处理灰度图像得到二值图，得到裂缝二值图后，拾取判定裂缝的高点坐标和低点坐标，利用公式得到裂缝参数；高低和低点坐标与裂缝参数公式如下：
平面二维裂缝公式如下所示：
z＝Ax+By+C (1)
式中，z为高度，其单位为m；x，y分别是X坐标和Y坐标；A，B，C为平面控制参数；
井眼方程如公式如下所示：
Figure BDA0003534757860000021
式中，x，y分别是X坐标和Y坐标；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
通过公式和，得到井壁上的三角函数曲线形式的裂缝曲线函数，如公式所示：
Figure BDA0003534757860000022
式中，A，B，C为平面控制参数；z为高度，其单位为m；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
已知高点坐标H(x0,y0)和低点坐标L(x1,y1)，则将公式改写为如下公式：
Figure BDA0003534757860000031
式中，θ为井周角，其单位为rad；z为高度，其单位为m；x0，y0为高点坐标，x1，y1为低坐标；
因此裂缝参数中的倾角、方位、倾向和曲线长度分别如以下公式所示：
Figure BDA0003534757860000032
β＝min(x0,x1) (6)
Figure BDA0003534757860000033
Figure BDA0003534757860000034
式中，α为倾角，其单位为rad；β为走向，其单位为rad；λ为倾向，其单位为rad；l为曲线长度，其单位为m。
进一步的，所述S2将处理得到裂缝参数，构建平面裂缝模型，利用井眼尺寸，构建地层井眼模型；所述井眼模型为长方体井眼模型，所述井眼模型的井眼半径与实际的井眼半径一致；所述长方体井眼模型的长度取决于模型的长度，其宽度大于3倍井眼直径；为满足精度条件和快速求解运算，设置为5倍。
进一步的，所述S3中，在3DEC离散元软件中赋值地层材料参数，包括弹性模量、泊松比、密度、摩擦角、内聚力和抗拉强度；为裂缝赋值力学参数，包括法向刚度和剪切刚度；为模型赋值边界条件，包括水平最大主应力、水平最小主应力、垂向应力，初始应力、边界位移和钻井液密度上下界限；最后设定二分法求解精度，确定迭代次数n；迭代次数的公式如下：
Figure BDA0003534757860000041
式中，x*为真值；Xk为第k次迭代计算值；b为钻井液密度上界；a为钻井液密度下界；ε为求解精度。
进一步的，所述S4中，利用3DEC离散元软件，确定钻井液密度r，r等于钻井液密度上界+钻井液密度下界/2；设定求解步长，运算模型，直到最大不平衡力保持不变或者为零；软件求解结束，后处理判定井壁单元是否满足稳定性条件，可设定判定条件如下：
①在单元的形心位置小于阈值1时，屈服单元的数量小于阈值2；
②在单元的形心位置小于位置阈值1时，井壁单元的最大水平位移小于阈值3；
③在单元的形心位置小于位置阈值1时，屈服单元的数量小于阈值且井壁单元的最大水平位移小于阈值3。
更进一步的，所述S5中，根据判定条件，当判定条件为真，钻井液密度上界为当前模型钻井液密度，钻井液密度下界为当前模型钻井液密度下界；当判定条件为假，钻井液密度上界为当前模型钻井液密度上界，钻井液密度下界为当前模型钻井液密度；更新钻井液密度上下界限。
更进一步的，所述S6中，根据钻井液密度上下界限，确定模型钻井液密度；重新唤起模型，执行n-1次步骤S4～S5。
本发明的有益效果：
本发明提出的一种裂缝性地层计算安全钻井密度的方法以地层成像测井图为基础，根据离散单元法，可准确计算裂缝性地层井壁稳定性，确定裂缝性地层的安全钻井液密度。
附图说明
图1为本发明的流程框图；
图2为实施例1中图像处理结果图，其图2中的(a)为原始成像测井图，(b)为原始成像测井灰度图，(c)为原始成像测井二值图，(d)为图像处理后裂缝井壁三角函数形式图；
图3为实施例1中步骤S2中建立的三维地质模型，其图3中的(a)为裂缝模型，(b)最后形成的运算模型；
图4为实施例1中运算结果图，其图4中的(a)为原始运算模型，(b)显示经过7次迭代计算，(c)为破坏单元在井壁的呈现；
图5为实施例2中运算结果图，其图5中的(a)为原始运算模型，(b)为示经过7次迭代计算，(c)为破坏单元在井壁的呈现。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图说明本发明的具体实施方式。
实施例一、本实施例中，如图1至图5所示，一种裂缝性地层计算安全钻井密度的方法，包括以下步骤：
S1、图像处理，识别井下裂缝；
S2、通过井下裂缝参数建立三维地质模型，利用井眼尺寸和长度建立钻井井眼模型；
S3、赋值模型材料参数、边界条件、初始钻井液密度上下界，然后计算精度；
S4、利用3DEC离散元软件，求解三维地质模型，判定井壁稳定性；
S5、利用二分法确定钻井液密度上下界；
S6、重复执行步骤S4～S5；
S7、达到设定的精度条件后，保存并输出安全钻井液密度。
进一步的，所述S1具体为先将地层成像测井图转换为灰度图，同时得到灰度颜色直方图，利用阈值方法等分割算法处理灰度图像得到二值图，得到裂缝二值图后，拾取判定裂缝的高点坐标和低点坐标，利用公式得到裂缝参数；高低和低点坐标与裂缝参数公式如下：
平面二维裂缝公式如下所示：
z＝Ax+By+C (1)
式中，z为高度，其单位为m；x，y分别是X坐标和Y坐标；A，B，C为平面控制参数；
井眼方程如公式如下所示：
Figure BDA0003534757860000051
式中，x，y分别是X坐标和Y坐标；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
通过公式和，得到井壁上的三角函数曲线形式的裂缝曲线函数，如公式所示：
Figure BDA0003534757860000061
式中，A，B，C为平面控制参数；z为高度，其单位为m；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
已知高点坐标H(x0,y0)和低点坐标L(x1,y1)，则将公式改写为如下公式：
Figure BDA0003534757860000062
式中，θ为井周角，其单位为rad；z为高度，其单位为m；x0，y0为高点坐标，x1，y1为低坐标；
因此裂缝参数中的倾角、方位、倾向和曲线长度分别如以下公式所示：
Figure BDA0003534757860000063
β＝min(x0,x1) (6)
Figure BDA0003534757860000064
Figure BDA0003534757860000065
式中，α为倾角，其单位为rad；β为走向，其单位为rad；λ为倾向，其单位为rad；l为曲线长度，其单位为m。
图像处理伪码：
Figure BDA0003534757860000071
本实施例进一步设置为：所述S2将处理得到裂缝参数，构建平面裂缝模型，利用井眼尺寸，构建地层井眼模型；所述井眼模型为长方体井眼模型，所述井眼模型的井眼半径与实际的井眼半径一致；所述长方体井眼模型的长度取决于模型的长度，其宽度大于3倍井眼直径；为满足精度条件和快速求解运算，设置为5倍。
通过3DEC离散元软件建立三维地质模型和钻井井眼模型，运行伪码如下：
Figure BDA0003534757860000072
本实施例进一步设置为：所述S3中，在3DEC离散元软件中赋值地层材料参数，包括弹性模量、泊松比、密度、摩擦角、内聚力和抗拉强度；为裂缝赋值力学参数，包括法向刚度和剪切刚度；为模型赋值边界条件，包括水平最大主应力、水平最小主应力、垂向应力，初始应力、边界位移和钻井液密度上下界限；最后设定二分法求解精度，确定迭代次数n；迭代次数的公式如下：
Figure BDA0003534757860000073
式中，x*为真值；Xk为第k次迭代计算值；b为钻井液密度上界；a为钻井液密度下界；ε为求解精度。
运行伪码如下：
Figure BDA0003534757860000081
本实施例进一步设置为：所述S4中，利用3DEC离散元软件，确定钻井液密度r，r等于钻井液密度上界+钻井液密度下界/2；设定求解步长，运算模型，直到最大不平衡力保持不变或者为零；软件求解结束，后处理判定井壁单元是否满足稳定性条件，可设定判定条件如下：
①在单元的形心位置小于阈值1时，屈服单元的数量小于阈值2；
②在单元的形心位置小于位置阈值1时，井壁单元的最大水平位移小于阈值3；
③在单元的形心位置小于位置阈值1时，屈服单元的数量小于阈值且井壁单元的最大水平位移小于阈值3。
S4利用离散元方法，求解三维地质模型，判定井壁稳定性，3DEC运行伪码下：
Figure BDA0003534757860000082
Figure BDA0003534757860000091
本实施例进一步设置为：所述S5中，根据判定条件，当判定条件为真，钻井液密度上界为当前模型钻井液密度，钻井液密度下界为当前模型钻井液密度下界；当判定条件为假，钻井液密度上界为当前模型钻井液密度上界，钻井液密度下界为当前模型钻井液密度；更新钻井液密度上下界限，即在计算机程序中，新的钻井液密度上下界限值替换旧的钻井液密度上下界限值。
S5利用二分法确定钻井液密度上下界，运行伪码如下：
Figure BDA0003534757860000092
本实施例进一步设置为：所述S6中，根据钻井液密度上下界限，确定模型钻井液密度；重新唤起模型，执行n-1次步骤S4～S5。运行伪码如下：
Figure BDA0003534757860000093
步骤S7，达到精度条件，保存输出安全钻井液密度。
实施例二、以下以某井裂缝性地层8400m-8404m段为例，对本发明的特征和性能作进一步的详细描述。
步骤S1，图像处理，识别井下裂缝；
执行S1图像处理伪码或函数，输入井眼直径0.075，井眼长度4m，读取成像测井图片，如图2(a)所示。利用rgb2gray函数转换为灰度图片，如图2(b)所示；输入阈值100，分割图像，得到二值图像，如图2(c)所示；在图像上识别裂缝，利用鼠标选择三角函数曲线最大值点和最小值点，得到裂缝几何参数。
Figure BDA0003534757860000101
裂缝参数如表1所示：
表1
Figure BDA0003534757860000102
步骤S2，通过井下裂缝参数，建立三维地质模型和钻井井眼模型；
首先利用ploy tunnel和delete range cy命令建立三维钻井井眼模型；井眼直径0.075m，长方体尺寸为0.75×0.75×4m；在步骤S1中，我们得到了9条裂缝，因此，使用重复使用9次Jset命令，建立裂缝模型；最后网格划分采用四面体网格划分，划分尺寸为0.03m，设置gen函数参数；裂缝模型如图3(a)所示，模型如图3(b)。
3DEC运行伪码下:
Figure BDA0003534757860000111
S3，赋值模型材料参数、边界条件、初始钻井液密度上下界、计算精度；
基质岩石参数，弹性模量40GPa,，密度2600kg/m3，泊松比0.21，摩擦角40°，抗拉强度8.5MPa，内聚力30.64MPa；
裂缝参数法向刚度9GPa，剪切刚度6GPa；
初始原地应力，最大水平主应力199MPa，最小水平主应力137MPa，垂向应力206MPa；
边界应力，最大水平主应力199MPa，最小水平主应力137MPa，垂向应力206MPa；
位移边界，xz平面yvel为0，yz平面xvel为0；
初始位移和初始速度均为0；
设定初始钻井液密度上下界限为2.09g/cm3和0g/cm3，保存至当前程序文件夹地址。设定求解精度0.01。
执行以下3DEC伪码或函数：
Figure BDA0003534757860000112
S4，利用离散元方法，求解三维地质模型，判定井壁稳定性。
设定运行步长2000次，井壁条件为单元形心小于0.08m；单元状态为非弹性就统计单元失效单元数量；以下3DEC代码中threshold1就为0.08，threshold2为1。
执行以下3DEC伪码或函数：
Figure BDA0003534757860000121
S5，利用二分法确定钻井液密度上下界；
在这一步中，模型稳定性的条件为单元失效数量占比全部单元数量小于10％，认为满足模型稳定性条件。执行以下3DEC伪码或函数：
Figure BDA0003534757860000122
Figure BDA0003534757860000131
S6，重复执行步骤S4～S5；我们设定的精度为0.01，钻井液密度上下界限为2.09g/cm3和0g/cm3，根据公式(8)，得到重复计算的次数为7次，在3DEC代码执行窗口复制S4和S5代码7份。
执行以下3DEC伪码或函数：
Figure BDA0003534757860000132
步骤S7，达到精度条件，保存输出安全钻井液密度，模拟结果如图4所示，最后安全钻井液密度为1.36g/cm3。
本发明也可对无裂缝地层计算安全钻井液密度，以下以无裂缝地层为例，对本发明的特征和性能作进一步的详细描述
步骤S1，图像处理，识别井下裂缝；
由于地层没有裂缝地层，可缺省第一步操作。
S2，通过井下裂缝参数，建立三维地质模型和钻井井眼模型；
缺省井下裂缝参数，直接建立三维钻井井眼模型；井眼直径0.075m，长方体尺寸为0.75×0.75×4m；网格划分采用四面体网格划分，划分尺寸为0.03m；设置函数参数，执行以下3DEC伪码或函数：
Figure BDA0003534757860000133
S3，赋值模型材料参数、边界条件、初始钻井液密度上下界、计算精度；
基质岩石参数，弹性模量40GPa,，密度2600kg/m3，泊松比0.21，摩擦角40°，抗拉强度8.5MPa，内聚力30.64MPa；
裂缝参数法向刚度9GPa，剪切刚度6GPa；
初始原地应力，最大水平主应力199MPa，最小水平主应力137MPa，垂向应力206MPa；
边界应力，最大水平主应力199MPa，最小水平主应力137MPa，垂向应力206MPa；
位移边界，xz平面yvel为0，yz平面xvel为0；
初始位移和初始速度均为0；
设定初始钻井液密度上下界限为2.09g/cm3和0g/cm3，保存至当前程序文件夹地址。设定求解精度0.01。
执行以下3DEC伪码或函数：
Figure BDA0003534757860000141
S4，利用离散元方法，求解三维地质模型，判定井壁稳定性。
设定运行步长2000次，井壁条件为单元形心小于0.08m；单元状态为非弹性就统计单元失效单元数量；以下3DEC代码中threshold1就为0.08，threshold2为1。
执行以下3DEC伪码或函数：
Figure BDA0003534757860000142
S5，利用二分法确定钻井液密度上下界；
在这一步中，模型稳定性的条件为单元失效数量占比全部单元数量小于10％，认为满足模型稳定性条件。执行以下3DEC伪码或函数：
Figure BDA0003534757860000143
Figure BDA0003534757860000151
S6，重复执行步骤S4～S5；我们设定的精度为0.01，钻井液密度上下界限为2.09g/cm3和0g/cm3，根据公式(8)，得到重复计算的次数为7次，在3DEC代码执行窗口复制S4和S5代码7份。
执行以下3DEC伪码或函数：
Figure BDA0003534757860000152
步骤S7，达到精度条件，保存输出安全钻井液密度，模拟结果如图5所示，最后安全钻井液密度为0.74g/cm3。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护的范围由所附的权利要求书及其等效物界定。

Claims (7)

1.一种裂缝性地层计算安全钻井密度的方法，其特征在于，包括以下步骤：
S1、图像处理，识别井下裂缝；
S2、通过井下裂缝参数建立三维地质模型，利用井眼尺寸和长度建立钻井井眼模型；
S3、赋值模型材料参数、边界条件、初始钻井液密度上下界，然后计算精度；
S4、利用3DEC离散元软件，求解三维地质模型，判定井壁稳定性；
S5、利用二分法确定钻井液密度上下界；
S6、重复执行步骤S4～S5；
S7、达到设定的精度条件后，保存并输出安全钻井液密度。
2.根据权利要求1所述的一种裂缝性地层计算安全钻井密度的方法，其特征在于：所述S1具体为先将地层成像测井图转换为灰度图，同时得到灰度颜色直方图，利用阈值方法等分割算法处理灰度图像得到二值图，得到裂缝二值图后，拾取判定裂缝的高点坐标和低点坐标，利用公式得到裂缝参数；高低和低点坐标与裂缝参数公式如下：
平面二维裂缝公式如下所示：
z＝Ax+By+C (1)
式中，z为高度，其单位为m；x，y分别是X坐标和Y坐标；A，B，C为平面控制参数；
井眼方程如公式如下所示：
Figure FDA0003534757850000011
式中，x，y分别是X坐标和Y坐标；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
通过公式和，得到井壁上的三角函数曲线形式的裂缝曲线函数，如公式所示：
Figure FDA0003534757850000012
式中，A，B，C为平面控制参数；z为高度，其单位为m；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
已知高点坐标H(x0,y0)和低点坐标L(x1,y1)，则将公式改写为如下公式：
Figure FDA0003534757850000021
式中，θ为井周角，其单位为rad；z为高度，其单位为m；x0，y0为高点坐标，x1，y1为低坐标；
因此裂缝参数中的倾角、方位、倾向和曲线长度分别如以下公式所示：
Figure FDA0003534757850000022
β＝min(x0,x1) (6)
Figure FDA0003534757850000023
Figure FDA0003534757850000024
式中，α为倾角，其单位为rad；β为走向，其单位为rad；λ为倾向，其单位为rad；l为曲线长度，其单位为m。
3.根据权利要求1所述的一种裂缝性地层计算安全钻井密度的方法，其特征在于：所述S2将处理得到裂缝参数，构建平面裂缝模型，利用井眼尺寸，构建地层井眼模型；所述井眼模型为长方体井眼模型，所述井眼模型的井眼半径与实际的井眼半径一致；所述长方体井眼模型的长度取决于模型的长度，其宽度大于3倍井眼直径；为满足精度条件和快速求解运算，设置为5倍。
4.根据权利要求1所述的一种裂缝性地层计算安全钻井密度的方法，其特征在于：所述S3中，在3DEC离散元软件中赋值地层材料参数，包括弹性模量、泊松比、密度、摩擦角、内聚力和抗拉强度；为裂缝赋值力学参数，包括法向刚度和剪切刚度；为模型赋值边界条件，包括水平最大主应力、水平最小主应力、垂向应力，初始应力、边界位移和钻井液密度上下界限；最后设定二分法求解精度，确定迭代次数n；迭代次数的公式如下：
Figure FDA0003534757850000031
式中，x*为真值；Xk为第k次迭代计算值；b为钻井液密度上界；a为钻井液密度下界；ε为求解精度。
5.根据权利要求1所述的一种裂缝性地层计算安全钻井密度的方法，其特征在于：所述S4中，利用3DEC离散元软件，确定钻井液密度r，r等于钻井液密度上界+钻井液密度下界/2；设定求解步长，运算模型，直到最大不平衡力保持不变或者为零；软件求解结束，后处理判定井壁单元是否满足稳定性条件，可设定判定条件如下：
①在单元的形心位置小于阈值1时，屈服单元的数量小于阈值2；
②在单元的形心位置小于位置阈值1时，井壁单元的最大水平位移小于阈值3；
③在单元的形心位置小于位置阈值1时，屈服单元的数量小于阈值且井壁单元的最大水平位移小于阈值3。
6.根据权利要求5所述的一种裂缝性地层计算安全钻井密度的方法，其特征在于：所述S5中，根据判定条件，当判定条件为真，钻井液密度上界为当前模型钻井液密度，钻井液密度下界为当前模型钻井液密度下界；当判定条件为假，钻井液密度上界为当前模型钻井液密度上界，钻井液密度下界为当前模型钻井液密度；更新钻井液密度上下界限。
7.根据权利要求6所述的一种裂缝性地层计算安全钻井密度的方法，其特征在于：所述S6中，根据钻井液密度上下界限，确定模型钻井液密度；重新唤起模型，执行n-1次步骤S4～S5。
CN202210216104.0A 2022-03-07 2022-03-07 一种裂缝性地层计算安全钻井密度的方法 Active CN114564778B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210216104.0A CN114564778B (zh) 2022-03-07 2022-03-07 一种裂缝性地层计算安全钻井密度的方法
US18/118,229 US11822046B2 (en) 2022-03-07 2023-03-07 Method for calculating safe drilling fluid density in fractured formation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210216104.0A CN114564778B (zh) 2022-03-07 2022-03-07 一种裂缝性地层计算安全钻井密度的方法

Publications (2)

Publication Number Publication Date
CN114564778A true CN114564778A (zh) 2022-05-31
CN114564778B CN114564778B (zh) 2022-11-18

Family

ID=81718193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210216104.0A Active CN114564778B (zh) 2022-03-07 2022-03-07 一种裂缝性地层计算安全钻井密度的方法

Country Status (2)

Country Link
US (1) US11822046B2 (zh)
CN (1) CN114564778B (zh)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110214879A1 (en) * 2010-03-03 2011-09-08 Baker Hughes Incorporated Tactile pressure sensing devices and methods for using same
US20130275099A1 (en) * 2012-04-17 2013-10-17 Schlumberger Technology Corporation Determining A Limit Of Failure In A Wellbore Wall
CN104806233A (zh) * 2015-02-12 2015-07-29 中国石油大学(北京) 一种预测弱面地层坍塌压力当量密度窗口的方法
CN108071389A (zh) * 2017-12-10 2018-05-25 西南石油大学 一种页岩强度准则及页岩气井井眼坍塌压力预测方法
CN108843313A (zh) * 2018-06-23 2018-11-20 东北石油大学 页岩地层钻井安全钻井液密度窗口设计方法
CN109377101A (zh) * 2018-11-30 2019-02-22 西南石油大学 一种基于风险控制模型的井壁稳定定量评价方法
CN109555515A (zh) * 2018-10-08 2019-04-02 中国石油大学(北京) 地层坍塌压力确定方法和装置
CN110032777A (zh) * 2019-03-25 2019-07-19 西南石油大学 一种基于不确定性分析的钻井液密度安全窗口估计方法
CN209327131U (zh) * 2019-06-20 2019-08-30 西南石油大学 一种煤岩三轴破坏测试实验装置
CN110259443A (zh) * 2019-07-12 2019-09-20 西南石油大学 一种基于3dec离散元的煤层井壁稳定性预测方法
CN111797490A (zh) * 2019-04-04 2020-10-20 中国石油天然气股份有限公司 钻井液安全密度设计的方法及设备
CN113221347A (zh) * 2021-05-07 2021-08-06 中国石油大学(北京) 一种井壁稳定性钻井优化方法、装置及设备
CN113935092A (zh) * 2021-10-23 2022-01-14 西南石油大学 一种裂缝性地层井壁稳定双重孔隙热弹性耦合评价方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110214879A1 (en) * 2010-03-03 2011-09-08 Baker Hughes Incorporated Tactile pressure sensing devices and methods for using same
US20130275099A1 (en) * 2012-04-17 2013-10-17 Schlumberger Technology Corporation Determining A Limit Of Failure In A Wellbore Wall
CN104806233A (zh) * 2015-02-12 2015-07-29 中国石油大学(北京) 一种预测弱面地层坍塌压力当量密度窗口的方法
CN108071389A (zh) * 2017-12-10 2018-05-25 西南石油大学 一种页岩强度准则及页岩气井井眼坍塌压力预测方法
CN108843313A (zh) * 2018-06-23 2018-11-20 东北石油大学 页岩地层钻井安全钻井液密度窗口设计方法
CN109555515A (zh) * 2018-10-08 2019-04-02 中国石油大学(北京) 地层坍塌压力确定方法和装置
CN109377101A (zh) * 2018-11-30 2019-02-22 西南石油大学 一种基于风险控制模型的井壁稳定定量评价方法
CN110032777A (zh) * 2019-03-25 2019-07-19 西南石油大学 一种基于不确定性分析的钻井液密度安全窗口估计方法
CN111797490A (zh) * 2019-04-04 2020-10-20 中国石油天然气股份有限公司 钻井液安全密度设计的方法及设备
CN209327131U (zh) * 2019-06-20 2019-08-30 西南石油大学 一种煤岩三轴破坏测试实验装置
CN110259443A (zh) * 2019-07-12 2019-09-20 西南石油大学 一种基于3dec离散元的煤层井壁稳定性预测方法
CN113221347A (zh) * 2021-05-07 2021-08-06 中国石油大学(北京) 一种井壁稳定性钻井优化方法、装置及设备
CN113935092A (zh) * 2021-10-23 2022-01-14 西南石油大学 一种裂缝性地层井壁稳定双重孔隙热弹性耦合评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
屈平等: "三维离散元在煤层水平井井壁稳定中的应用", 《石油学报》 *
张光福等: "基于3DEC 离散元的煤层井壁稳定性", 《科学技术与工程》 *
陈平等: "基于井壁稳定分析的井眼轨迹优化方法", 《天然气工业》 *

Also Published As

Publication number Publication date
CN114564778B (zh) 2022-11-18
US20230280493A1 (en) 2023-09-07
US11822046B2 (en) 2023-11-21

Similar Documents

Publication Publication Date Title
RU2602403C1 (ru) Управление микросейсмическими данными для отслеживания трещин
Deng et al. Pillar design by combining finite element methods, neural networks and reliability: a case study of the Feng Huangshan copper mine, China
Vardakos et al. Parameter identification in numerical modeling of tunneling using the Differential Evolution Genetic Algorithm (DEGA)
AU2013402201B2 (en) In-situ wellbore, core and cuttings information system
EP2583214B1 (en) Systems and methods for wellbore optimization
EP3455654A1 (en) Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof
EP3455459A2 (en) Three-dimensional fracture abundance evaluation of subsurface formations
Zhang et al. A geometrically and locally adaptive remeshing method for finite difference modeling of mining-induced surface subsidence
CN113221347B (zh) 一种井壁稳定性钻井优化方法、装置及设备
CN114564778B (zh) 一种裂缝性地层计算安全钻井密度的方法
US10650107B2 (en) Three-dimensional subsurface formation evaluation using projection-based area operations
CN116611248A (zh) 基于巴西圆盘实验考虑的各向异性岩石真实剪切断裂计算方法
Ali et al. Numerical analysis of stress distribution during tunneling in clay stone rock
CN113356843B (zh) 针对地层的井壁稳定性分析方法、装置、介质及设备
Wei et al. Numerical study of damage to rock surrounding an underground coal roadway excavation
EP3427085A1 (en) Determining rupture envelopes of a fault system
Zhang et al. Numerical simulation of rock failure process with a 3D grain-based rock model
Rahaman et al. Lower bound limit analysis for an axisymmetric problem in rock mass using power cone programming
Florkowska The application of numerical analysis in determining the state of the rock mass around directional wells
Chowdhury et al. Mechanical borehole stability analysis using phase2
Yazdi et al. Effect of nail layout variability on soil nailed wall analysis
CN116842789B (zh) 一种考虑生产时间窗的加密井井壁稳定分析方法
Zoughy et al. Assessing the Minimum Safe Drilling Mud Weight to Avoid Damage Zone by Frictional-Strengthening Cohesion-Weakening Models
Yufin et al. Several aspects of formation of structures of computer programs, selection and algorithmization of medium models in geomechanics
Wang et al. Tunneling issues regarding the rock tunnel-shaft intersection in Taiwan

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant