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

1.本发明涉及石油与天然气工程中的钻井实践、地热钻探、深地科学探索领域，尤其涉及一种裂缝性地层计算安全钻井密度的方法。


背景技术：

2.随着石油与天然气工程的发展，钻井深度越来越深。钻井遇到的地层，往往不是完整的岩石，而是包含结构面的岩体。在钻遇破碎带、断层、裂缝时，这种问题现象尤为突出。钻井井眼需要钻井液压力的支撑，不然就会出现井壁缩颈、坍塌等问题。进而造成卡钻，埋钻杆等问题，严重情况时将报废整个井眼。
3.因此，钻井所需的安全钻井液压力或者安全钻井液的密度对于钻井工程的成败起着决定性的作用。一般来说，目前计算钻井液密度是利用连续性介质理论，开展井壁应力计算，通过研究不同钻井液密度对井壁稳定的影响程度，利用二分法后者牛顿迭代法计算得到所需钻井液密度。常见的计算流程如下：
4.1.利用farihurst方程计算井壁应力；
5.2.使用mohr-coulomb强度准则判定井壁稳定性；
6.3.利用二分法或者牛顿迭代法计算钻井液密度。
7.在钻遇破碎带、断层、裂缝等地层时，由于基本假设和理论上的不同，无法考虑地层的几何产状等因素对井壁稳定和钻井液密度的作用。因此，上述方法不完全适用于破裂带、断层、裂缝性地层钻井工程的钻井液密度计算。


技术实现要素：

8.本发明的目的在于克服现有技术的不足，提供一种裂缝性地层计算安全钻井密度的方法，该方法在破碎带、断层、裂缝性地层的钻井液密度计算更加准确，更好的指导裂缝性地层的钻井工程钻进。
9.本发明的目的是通过以下技术方案来实现的：
10.一种裂缝性地层计算安全钻井密度的方法，包括以下步骤：
11.s1、图像处理，识别井下裂缝；
12.s2、通过井下裂缝参数建立三维地质模型，利用井眼尺寸和长度建立钻井井眼模型；
13.s3、赋值模型材料参数、边界条件、初始钻井液密度上下界，然后计算精度；
14.s4、利用3dec离散元软件，求解三维地质模型，判定井壁稳定性；
15.s5、利用二分法确定钻井液密度上下界；
16.s6、重复执行步骤s4～s5；
17.s7、达到设定的精度条件后，保存并输出安全钻井液密度。
18.进一步的，所述s1具体为先将地层成像测井图转换为灰度图，同时得到灰度颜色直方图，利用阈值方法等分割算法处理灰度图像得到二值图，得到裂缝二值图后，拾取判定
裂缝的高点坐标和低点坐标，利用公式得到裂缝参数；高低和低点坐标与裂缝参数公式如下：
19.平面二维裂缝公式如下所示：
20.z＝ax+by+c
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
21.式中，z为高度，其单位为m；x，y分别是x坐标和y坐标；a，b，c为平面控制参数；
22.井眼方程如公式如下所示：
[0023][0024]
式中，x，y分别是x坐标和y坐标；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
[0025]
通过公式和，得到井壁上的三角函数曲线形式的裂缝曲线函数，如公式所示：
[0026][0027]
式中，a，b，c为平面控制参数；z为高度，其单位为m；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
[0028]
已知高点坐标h(x0,y0)和低点坐标l(x1,y1)，则将公式改写为如下公式：
[0029][0030]
式中，θ为井周角，其单位为rad；z为高度，其单位为m；x0，y0为高点坐标，x1，y1为低坐标；
[0031]
因此裂缝参数中的倾角、方位、倾向和曲线长度分别如以下公式所示：
[0032][0033]
β＝min(x0,x1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0034][0035][0036]
式中，α为倾角，其单位为rad；β为走向，其单位为rad；λ为倾向，其单位为rad；l为曲线长度，其单位为m。
[0037]
进一步的，所述s2将处理得到裂缝参数，构建平面裂缝模型，利用井眼尺寸，构建地层井眼模型；所述井眼模型为长方体井眼模型，所述井眼模型的井眼半径与实际的井眼半径一致；所述长方体井眼模型的长度取决于模型的长度，其宽度大于3倍井眼直径；为满
足精度条件和快速求解运算，设置为5倍。
[0038]
进一步的，所述s3中，在3dec离散元软件中赋值地层材料参数，包括弹性模量、泊松比、密度、摩擦角、内聚力和抗拉强度；为裂缝赋值力学参数，包括法向刚度和剪切刚度；为模型赋值边界条件，包括水平最大主应力、水平最小主应力、垂向应力，初始应力、边界位移和钻井液密度上下界限；最后设定二分法求解精度，确定迭代次数n；迭代次数的公式如下：
[0039][0040]
式中，x*为真值；xk为第k次迭代计算值；b为钻井液密度上界；a为钻井液密度下界；ε为求解精度。
[0041]
进一步的，所述s4中，利用3dec离散元软件，确定钻井液密度r，r等于钻井液密度上界+钻井液密度下界/2；设定求解步长，运算模型，直到最大不平衡力保持不变或者为零；软件求解结束，后处理判定井壁单元是否满足稳定性条件，可设定判定条件如下：
[0042]
①
在单元的形心位置小于阈值1时，屈服单元的数量小于阈值2；
[0043]
②
在单元的形心位置小于位置阈值1时，井壁单元的最大水平位移小于阈值3；
[0044]
③
在单元的形心位置小于位置阈值1时，屈服单元的数量小于阈值且井壁单元的最大水平位移小于阈值3。
[0045]
更进一步的，所述s5中，根据判定条件，当判定条件为真，钻井液密度上界为当前模型钻井液密度，钻井液密度下界为当前模型钻井液密度下界；当判定条件为假，钻井液密度上界为当前模型钻井液密度上界，钻井液密度下界为当前模型钻井液密度；更新钻井液密度上下界限。
[0046]
更进一步的，所述s6中，根据钻井液密度上下界限，确定模型钻井液密度；重新唤起模型，执行n-1次步骤s4～s5。
[0047]
本发明的有益效果：
[0048]
本发明提出的一种裂缝性地层计算安全钻井密度的方法以地层成像测井图为基础，根据离散单元法，可准确计算裂缝性地层井壁稳定性，确定裂缝性地层的安全钻井液密度。
附图说明
[0049]
图1为本发明的流程框图；
[0050]
图2为实施例1中图像处理结果图，其图2中的(a)为原始成像测井图，(b)为原始成像测井灰度图，(c)为原始成像测井二值图，(d)为图像处理后裂缝井壁三角函数形式图；
[0051]
图3为实施例1中步骤s2中建立的三维地质模型，其图3中的(a)为裂缝模型，(b)最后形成的运算模型；
[0052]
图4为实施例1中运算结果图，其图4中的(a)为原始运算模型，(b)显示经过7次迭代计算，(c)为破坏单元在井壁的呈现；
[0053]
图5为实施例2中运算结果图，其图5中的(a)为原始运算模型，(b)为示经过7次迭代计算，(c)为破坏单元在井壁的呈现。
具体实施方式
[0054]
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图说明本发明的具体实施方式。
[0055]
实施例一、本实施例中，如图1至图5所示，一种裂缝性地层计算安全钻井密度的方法，包括以下步骤：
[0056]
s1、图像处理，识别井下裂缝；
[0057]
s2、通过井下裂缝参数建立三维地质模型，利用井眼尺寸和长度建立钻井井眼模型；
[0058]
s3、赋值模型材料参数、边界条件、初始钻井液密度上下界，然后计算精度；
[0059]
s4、利用3dec离散元软件，求解三维地质模型，判定井壁稳定性；
[0060]
s5、利用二分法确定钻井液密度上下界；
[0061]
s6、重复执行步骤s4～s5；
[0062]
s7、达到设定的精度条件后，保存并输出安全钻井液密度。
[0063]
进一步的，所述s1具体为先将地层成像测井图转换为灰度图，同时得到灰度颜色直方图，利用阈值方法等分割算法处理灰度图像得到二值图，得到裂缝二值图后，拾取判定裂缝的高点坐标和低点坐标，利用公式得到裂缝参数；高低和低点坐标与裂缝参数公式如下：
[0064]
平面二维裂缝公式如下所示：
[0065]
z＝ax+by+c
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0066]
式中，z为高度，其单位为m；x，y分别是x坐标和y坐标；a，b，c为平面控制参数；
[0067]
井眼方程如公式如下所示：
[0068][0069]
式中，x，y分别是x坐标和y坐标；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
[0070]
通过公式和，得到井壁上的三角函数曲线形式的裂缝曲线函数，如公式所示：
[0071][0072]
式中，a，b，c为平面控制参数；z为高度，其单位为m；r为井眼半径，其单位为m；θ为井周角，其单位为rad；
[0073]
已知高点坐标h(x0,y0)和低点坐标l(x1,y1)，则将公式改写为如下公式：
[0074][0075]
式中，θ为井周角，其单位为rad；z为高度，其单位为m；x0，y0为高点坐标，x1，y1为低坐标；
[0076]
因此裂缝参数中的倾角、方位、倾向和曲线长度分别如以下公式所示：
[0077][0078]
β＝min(x0,x1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0079][0080][0081]
式中，α为倾角，其单位为rad；β为走向，其单位为rad；λ为倾向，其单位为rad；l为曲线长度，其单位为m。
[0082]
图像处理伪码：
[0083][0084]
本实施例进一步设置为：所述s2将处理得到裂缝参数，构建平面裂缝模型，利用井眼尺寸，构建地层井眼模型；所述井眼模型为长方体井眼模型，所述井眼模型的井眼半径与实际的井眼半径一致；所述长方体井眼模型的长度取决于模型的长度，其宽度大于3倍井眼直径；为满足精度条件和快速求解运算，设置为5倍。
[0085]
通过3dec离散元软件建立三维地质模型和钻井井眼模型，运行伪码如下：
[0086]
[0087]
本实施例进一步设置为：所述s3中，在3dec离散元软件中赋值地层材料参数，包括弹性模量、泊松比、密度、摩擦角、内聚力和抗拉强度；为裂缝赋值力学参数，包括法向刚度和剪切刚度；为模型赋值边界条件，包括水平最大主应力、水平最小主应力、垂向应力，初始应力、边界位移和钻井液密度上下界限；最后设定二分法求解精度，确定迭代次数n；迭代次数的公式如下：
[0088][0089]
式中，x*为真值；xk为第k次迭代计算值；b为钻井液密度上界；a为钻井液密度下界；ε为求解精度。
[0090]
运行伪码如下：
[0091][0092]
本实施例进一步设置为：所述s4中，利用3dec离散元软件，确定钻井液密度r，r等于钻井液密度上界+钻井液密度下界/2；设定求解步长，运算模型，直到最大不平衡力保持不变或者为零；软件求解结束，后处理判定井壁单元是否满足稳定性条件，可设定判定条件如下：
[0093]
①
在单元的形心位置小于阈值1时，屈服单元的数量小于阈值2；
[0094]
②
在单元的形心位置小于位置阈值1时，井壁单元的最大水平位移小于阈值3；
[0095]
③
在单元的形心位置小于位置阈值1时，屈服单元的数量小于阈值且井壁单元的最大水平位移小于阈值3。
[0096]
s4利用离散元方法，求解三维地质模型，判定井壁稳定性，3dec运行伪码下：
[0097][0098][0099]
本实施例进一步设置为：所述s5中，根据判定条件，当判定条件为真，钻井液密度上界为当前模型钻井液密度，钻井液密度下界为当前模型钻井液密度下界；当判定条件为假，钻井液密度上界为当前模型钻井液密度上界，钻井液密度下界为当前模型钻井液密度；更新钻井液密度上下界限，即在计算机程序中，新的钻井液密度上下界限值替换旧的钻井液密度上下界限值。
[0100]
s5利用二分法确定钻井液密度上下界，运行伪码如下：
[0101][0102]
本实施例进一步设置为：所述s6中，根据钻井液密度上下界限，确定模型钻井液密度；重新唤起模型，执行n-1次步骤s4～s5。运行伪码如下：
[0103][0104]
步骤s7，达到精度条件，保存输出安全钻井液密度。
[0105]
实施例二、以下以某井裂缝性地层8400m-8404m段为例，对本发明的特征和性能作进一步的详细描述。
[0106]
步骤s1，图像处理，识别井下裂缝；
[0107]
执行s1图像处理伪码或函数，输入井眼直径0.075，井眼长度4m，读取成像测井图片，如图2(a)所示。利用rgb2gray函数转换为灰度图片，如图2(b)所示；输入阈值100，分割图像，得到二值图像，如图2(c)所示；在图像上识别裂缝，利用鼠标选择三角函数曲线最大值点和最小值点，得到裂缝几何参数。
[0108][0109]
裂缝参数如表1所示：
[0110]
表1
[0111][0112]
步骤s2，通过井下裂缝参数，建立三维地质模型和钻井井眼模型；
[0113]
首先利用ploy tunnel和delete range cy命令建立三维钻井井眼模型；井眼直径0.075m，长方体尺寸为0.75
×
0.75
×
4m；在步骤s1中，我们得到了9条裂缝，因此，使用重复
使用9次jset命令，建立裂缝模型；最后网格划分采用四面体网格划分，划分尺寸为0.03m，设置gen函数参数；裂缝模型如图3(a)所示，模型如图3(b)。
[0114]
3dec运行伪码下:
[0115][0116]
s3，赋值模型材料参数、边界条件、初始钻井液密度上下界、计算精度；
[0117]
基质岩石参数，弹性模量40gpa,，密度2600kg/m3，泊松比0.21，摩擦角40
°
，抗拉强度8.5mpa，内聚力30.64mpa；
[0118]
裂缝参数法向刚度9gpa，剪切刚度6gpa；
[0119]
初始原地应力，最大水平主应力199mpa，最小水平主应力137mpa，垂向应力206mpa；
[0120]
边界应力，最大水平主应力199mpa，最小水平主应力137mpa，垂向应力206mpa；
[0121]
位移边界，xz平面yvel为0，yz平面xvel为0；
[0122]
初始位移和初始速度均为0；
[0123]
设定初始钻井液密度上下界限为2.09g/cm3和0g/cm3，保存至当前程序文件夹地址。设定求解精度0.01。
[0124]
执行以下3dec伪码或函数：
[0125][0126]
s4，利用离散元方法，求解三维地质模型，判定井壁稳定性。
[0127]
设定运行步长2000次，井壁条件为单元形心小于0.08m；单元状态为非弹性就统计单元失效单元数量；以下3dec代码中threshold1就为0.08，threshold2为1。
[0128]
执行以下3dec伪码或函数：
[0129][0130]
s5，利用二分法确定钻井液密度上下界；
[0131]
在这一步中，模型稳定性的条件为单元失效数量占比全部单元数量小于10％，认为满足模型稳定性条件。执行以下3dec伪码或函数：
[0132][0133][0134]
s6，重复执行步骤s4～s5；我们设定的精度为0.01，钻井液密度上下界限为2.09g/cm3和0g/cm3，根据公式(8)，得到重复计算的次数为7次，在3dec代码执行窗口复制s4和s5代码7份。
[0135]
执行以下3dec伪码或函数：
[0136][0137]
步骤s7，达到精度条件，保存输出安全钻井液密度，模拟结果如图4所示，最后安全钻井液密度为1.36g/cm3。
[0138]
本发明也可对无裂缝地层计算安全钻井液密度，以下以无裂缝地层为例，对本发明的特征和性能作进一步的详细描述
[0139]
步骤s1，图像处理，识别井下裂缝；
[0140]
由于地层没有裂缝地层，可缺省第一步操作。
[0141]
s2，通过井下裂缝参数，建立三维地质模型和钻井井眼模型；
[0142]
缺省井下裂缝参数，直接建立三维钻井井眼模型；井眼直径0.075m，长方体尺寸为0.75
×
0.75
×
4m；网格划分采用四面体网格划分，划分尺寸为0.03m；设置函数参数，执行以下3dec伪码或函数：
[0143][0144]
s3，赋值模型材料参数、边界条件、初始钻井液密度上下界、计算精度；
[0145]
基质岩石参数，弹性模量40gpa,，密度2600kg/m3，泊松比0.21，摩擦角40
°
，抗拉强度8.5mpa，内聚力30.64mpa；
[0146]
裂缝参数法向刚度9gpa，剪切刚度6gpa；
[0147]
初始原地应力，最大水平主应力199mpa，最小水平主应力137mpa，垂向应力206mpa；
[0148]
边界应力，最大水平主应力199mpa，最小水平主应力137mpa，垂向应力206mpa；
[0149]
位移边界，xz平面yvel为0，yz平面xvel为0；
[0150]
初始位移和初始速度均为0；
[0151]
设定初始钻井液密度上下界限为2.09g/cm3和0g/cm3，保存至当前程序文件夹地址。设定求解精度0.01。
[0152]
执行以下3dec伪码或函数：
[0153][0154]
s4，利用离散元方法，求解三维地质模型，判定井壁稳定性。
[0155]
设定运行步长2000次，井壁条件为单元形心小于0.08m；单元状态为非弹性就统计单元失效单元数量；以下3dec代码中threshold1就为0.08，threshold2为1。
[0156]
执行以下3dec伪码或函数：
[0157][0158]
s5，利用二分法确定钻井液密度上下界；
[0159]
在这一步中，模型稳定性的条件为单元失效数量占比全部单元数量小于10％，认为满足模型稳定性条件。执行以下3dec伪码或函数：
[0160]
[0161][0162]
s6，重复执行步骤s4～s5；我们设定的精度为0.01，钻井液密度上下界限为2.09g/cm3和0g/cm3，根据公式(8)，得到重复计算的次数为7次，在3dec代码执行窗口复制s4和s5代码7份。
[0163]
执行以下3dec伪码或函数：
[0164][0165]
步骤s7，达到精度条件，保存输出安全钻井液密度，模拟结果如图5所示，最后安全钻井液密度为0.74g/cm3。
[0166]
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护的范围由所附的权利要求书及其等效物界定。