CN114526067A - 不同强度准则协同作用下定向井井壁坍塌压力评价方法 - Google Patents
不同强度准则协同作用下定向井井壁坍塌压力评价方法 Download PDFInfo
- Publication number
- CN114526067A CN114526067A CN202210187418.2A CN202210187418A CN114526067A CN 114526067 A CN114526067 A CN 114526067A CN 202210187418 A CN202210187418 A CN 202210187418A CN 114526067 A CN114526067 A CN 114526067A
- Authority
- CN
- China
- Prior art keywords
- rock
- collapse pressure
- mpa
- different
- strength criteria
- 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.)
- Pending
Links
- 230000002195 synergetic effect Effects 0.000 title claims abstract description 37
- 238000011156 evaluation Methods 0.000 title claims abstract description 15
- 239000011435 rock Substances 0.000 claims abstract description 121
- 238000004364 calculation method Methods 0.000 claims abstract description 45
- 238000012360 testing method Methods 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 23
- 239000011148 porous material Substances 0.000 claims abstract description 11
- 230000035699 permeability Effects 0.000 claims abstract description 7
- 238000005553 drilling Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 22
- 239000000126 substance Substances 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 238000013461 design Methods 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 7
- 238000011161 development Methods 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 108700041286 delta Proteins 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- BTCSSZJGUNDROE-UHFFFAOYSA-N gamma-aminobutyric acid Chemical compound NCCCC(O)=O BTCSSZJGUNDROE-UHFFFAOYSA-N 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- 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
- E21B49/006—Measuring wall stresses in the borehole
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种不同强度准则协同作用下定向井井壁坍塌压力评价方法，包括：利用测录井、钻完井、室内测试资料，确定所钻地层的垂深、地层孔隙压力、上覆岩层压力、最大水平地应力、最小水平地应力和最大水平地应力方位角、岩石内聚力、岩石内摩擦角、岩石泊松比、岩石孔隙度、Biot系数和井壁渗透系数；建立不同岩石强度准则下的坍塌压力计算模型；计算不同岩石强度准则之间的相对重要程度；确定出不同岩石强度准的权重系数；建立不同强度准则协同作用下定向井坍塌压力计算模型。本发明由于引入了权重系数，结合了其他准则的优点，获得了不同强度准则的优势协同作用，在给出合理权重系数的情况下可以更为准确的计算出坍塌压力当量密度值。
Description
技术领域
本发明涉及不同强度准则协同作用下定向井井壁坍塌压力评价方法，属于石油勘探开发技术领域。
背景技术
随着我国经济的飞速发展，工业化进程逐渐深入，对石油等资源的消耗越来越多，非常规油气和海洋油气资源受到了越来越多的关注。近年来，由于定向井、大斜度井、水平井、大位移井技术的突破，非常规油气和海洋油气成为了油气勘探开发的热点领域。但是，采用定向井、大斜度井、水平井、大位移井等复杂结构井，在钻井的过程中不得不面临一系列的工程技术难题和挑战，其中，井壁垮塌问题是非常棘手的问题之一。
从力学角度来看，由于钻开地层形成井眼后，打破了原地应力平衡，导致井眼周围地层中应力重新分布，当重新分布的应力超过地层强度时，就可能发生井壁垮塌。定向井井壁坍塌压力的计算的一般思路：综合考虑地应力、井眼轨迹、井筒压力等因素计算井壁应力分布，然后将计算得到的应力带入到岩石强度准则进行判断，当井壁应力等于岩石强度时的临界井筒压力即为维持井壁稳定所需的坍塌压力。
一般，井壁应力主要采用基于线弹性理论的Bradely公式计算，而岩石强度准则有很多种选择，常用岩石强度准则包括Mohr-Coulomb准则、Mogi-Coulomb准则、Drucker-Prager准则、Hoek-Brown准则、修正Lade准则、修正Wiebols-Cook准则等等。在井壁坍塌压力评价过程中，之所以采用不同强度准则主要是因为岩石强度准则不能完全适用于所有类型的岩石，这导致井壁坍塌压力评价方法往往具有一定的局限性，不能准确计算出任意类型岩石或地层的坍塌压力，即坍塌压力评价误差较大，只能适用于一种或几种特殊的岩石或地层，缺少一种通用的井壁坍塌压力评价方法。
为此，本发明通过综合Mohr-Coulomb准则、Mogi-Coulomb准则、Drucker-Prager准则、修正Lade准则、修正Wiebols-Cook准则等多种不同类型的强度准则，发挥不同强度准则的优势和协同作用，从而发明出一种不同强度准则协同作用下定向井井壁坍塌压力评价方法，该方法可以适用于不同类型的岩石或地层，是一种通用的井壁坍塌压力评价方法。
发明内容
为了克服现有技术中的问题，本发明提供一种不同强度准则协同作用下定向井井壁坍塌压力评价方法。
本发明解决上述技术问题所提供的技术方案是：不同强度准则协同作用下定向井井壁坍塌压力评价方法，包括以下步骤：
S1、利用测录井、钻完井、室内测试资料，确定所钻地层的垂深、地层孔隙压力、上覆岩层压力、最大水平地应力、最小水平地应力和最大水平地应力方位角；
S2、利用测录井和室内岩石力学测试资料，确定所钻地层的岩石内聚力、岩石内摩擦角、岩石泊松比、岩石孔隙度、Biot系数和井壁渗透系数；
S3、建立不同岩石强度准则下的坍塌压力计算模型；
S4、对所钻地层的N组岩样进行不同岩石强度准则下的坍塌压力计算，并与通过岩石力学试验获得的数据进行拟合，确定符合不同岩石强度准则的岩样个数，再通过计算相应的比值，得到不同岩石强度准则之间的相对重要程度；
S5、根据不同岩石强度准则之间的相对重要程度，并采用层次分析法，建立比较系数矩阵A，对所得权重系数向量进行一致性检验，并确定出不同岩石强度准的权重系数；
S6、结合不同强度准则计算的坍塌压力和步骤S5确定的权重系数，建立不同强度准则协同作用下定向井坍塌压力计算模型，并计算得到不同强度准则协同作用下定向井井壁坍塌压力。
进一步的技术方案，所述步骤S3的具体过程为：
S31、计算所钻地层的井壁应力分量；
S32、根据井壁应力分量计算所钻地层的井壁有效主应力；
S33、根据井壁有效主应力建立不同强度准则下的坍塌压力非线性函数；
S34、利用二分法计算不同强度准则下的坍塌压力非线性函数，得到不同强度准则对应的坍塌压力计算模型。
进一步的技术方案，所述井壁应力分量的计算公式为：
其中，
式中：σrr、σθθ、σzz分别为井壁径向、切向和轴向应力分量，MPa；τθz、τrθ、τrz分别为井壁3个剪应力分量，MPa；pm为井筒压力，MPa；pp为地层孔隙压力，MPa；σv、σH、σh分别为上覆岩石压力、最大水平地应力和最小水平地应力，MPa；δ表示井壁渗透系数，无因次，δ＝0为井壁完全不渗透、δ＝1为井壁完全渗透；φ为岩石孔隙度；ψ为井斜角，°；Ω为井斜方位与最大水平地应力方位夹角，°；θ为井周角，°；v为岩石泊松比；K1为渗流效应系数；α为Biot系数；A、B、C、D、E、F、G、H、J为坐标变换系数。
进一步的技术方案，所述井壁有效主应力的计算公式为：
式中：σ1为井壁最大主应力，MPa；σ2为井壁中间主应力，MPa；σ3为井壁最小主应力，MPa；α为Biot系数；pp为地层孔隙压力，MPa。
进一步的技术方案，所述不同强度准则下的坍塌压力非线性函数包括：
井壁坍塌的非线性函数1：
式中：
为Mohr-Coulomb准则计算得到的坍塌压力，MPa；c为岩石内聚力，MPa；
为岩石内摩擦角，°；
井壁坍塌的非线性函数2：
其中，
式中：
为Mogi-Coulomb准则计算得到的坍塌压力，MPa；a、b为材料常数，取决于岩石内聚力和内摩擦角；
井壁坍塌的非线性函数3：
其中，
式中：
为Drucker-Prager准则计算得到的坍塌压力，MPa；α'、k为材料常数，取决于内聚力和内摩擦角；J1为平均有效应力；J2为应力偏量第二不变量；
井壁坍塌的非线性函数4：
其中，
式中：
为修正Lade准则计算得到的坍塌压力，MPa；S、η为材料常数，取决于内聚力和内摩擦角；
井壁坍塌的非线性函数5：
其中，
式中：
为修正Wiebols-Cook准则计算得到的坍塌压力，MPa；A'、B'、C'为材料常数，取决于内聚力和内摩擦角；UCS为岩石单轴强度，MPa；J1为平均有效应力；J2为应力偏量第二不变量。
进一步的技术方案，所述步骤S4的具体过程为：
S41、对所钻地层的N组岩样进行不同岩石强度准则下的坍塌压力计算得到计算数据，再通过室内岩石力学试验，对所钻地层的N组岩样进行不同围压下的三轴岩石力学强度测试，得到试验数据；并将计算数据与试验数据进行拟合；
S42、分别统计测试结果符合不同岩石强度准则的各个岩样个数，其中计算数据与试验数据的相对误差在3％以内为符合标准；
S43、通过计算符合不同岩石强度准则的各个岩样个数之间的比值，得到两两准则之间的相对重要程度。
进一步的技术方案，所述不同岩石强度准则包括Mohr-Coulomb准则、Mogi-Coulomb准则、Drucker-Prager准则、修正Lade准则和修正Wiebols-Cook准则。
进一步的技术方案，所述步骤S5的具体步骤为：
S51、根据不同岩石强度准则之间的相对重要程度建立不同强度准则之间的比较系数矩阵A；
S52、计算比较系数矩阵A的特征向量；
S53、对所得权重系数向量进行一致性检验，并确定出不同岩石强度准的权重系数。
进一步的技术方案，所述不同强度准则协同作用下定向井坍塌压力计算模型为：
式中，W1为Mohr-Coulomb准则的权重系数；W2为Mogi-Coulomb准则的权重系数；W3为Drucker-Prager准则的权重系数；W4为修正Lade准则的权重系数；W5为修正Wiebols-Cook准则的权重系数；
为Mohr-Coulomb准则计算得到的坍塌压力，MPa；
为Mogi-Coulomb准则计算得到的坍塌压力，MPa；
为Drucker-Prager准则计算得到的坍塌压力，MPa；
为修正Lade准则计算得到的坍塌压力，MPa；
为修正Wiebols-Cook准则计算得到的坍塌压力，MPa；pc为不同强度准则协同作用下定向井井壁坍塌压力，MPa。
本发明具有以下有益效果：本发明由于引入了权重系数，结合了其他准则的优点，获得了不同强度准则的优势协同作用，在给出合理权重系数的情况下可以更为准确的计算出坍塌压力当量密度值。
附图说明
图1为实施流程图；
图2为井斜方位0°情况下坍塌压力计算结果对比图；
图3为井斜方位15°情况下坍塌压力计算结果对比图；
图4为井斜方位30°情况下坍塌压力计算结果对比图；
图5为井斜方位45°情况下坍塌压力计算结果对比图；
图6为井斜方位60°情况下坍塌压力计算结果对比图；
图7为井斜方位75°情况下坍塌压力计算结果对比图；
图8为井斜方位90°情况下坍塌压力计算结果对比图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的不同强度准则协同作用下定向井井壁坍塌压力评价方法，包括以下步骤：
S1、利用测录井、钻完井、室内测试资料，确定所钻地层的垂深、地层孔隙压力、上覆岩层压力、最大水平地应力、最小水平地应力和最大水平地应力方位角；
S2、利用测录井和室内岩石力学测试资料，确定所钻地层的岩石内聚力、岩石内摩擦角、岩石泊松比、岩石孔隙度、Biot系数和井壁渗透系数；
S3、建立不同岩石强度准则下的坍塌压力计算模型，具体包括以下步骤：
S31、计算所钻地层的井壁应力分量；
其中，
式中：σrr、σθθ、σzz分别为井壁径向、切向和轴向应力分量，MPa；τθz、τrθ、τrz分别为井壁3个剪应力分量，MPa；pm为井筒压力，MPa；pp为地层孔隙压力，MPa；σv、σH、σh分别为上覆岩石压力、最大水平地应力和最小水平地应力，MPa；δ表示井壁渗透系数，无因次，δ＝0为井壁完全不渗透、δ＝1为井壁完全渗透；φ为岩石孔隙度；ψ为井斜角，°；Ω为井斜方位与最大水平地应力方位夹角，°；θ为井周角，°；v为岩石泊松比；K1为渗流效应系数；α为Biot系数；A、B、C、D、E、F、G、H、J为坐标变换系数；
S32、根据井壁应力分量计算所钻地层的井壁有效主应力；
式中，σ1为井壁最大主应力，MPa；σ2为井壁中间主应力，MPa；σ3为井壁最小主应力，MPa；α为Biot系数；pp为地层孔隙压力，MPa；
S33、根据井壁有效主应力建立不同强度准则下的坍塌压力非线性函数；
(1)将所钻地层的井壁有效主应力带入Mohr-Coulomb准则，得到求解井壁坍塌的非线性函数1；
式中：
为Mohr-Coulomb准则计算得到的坍塌压力，MPa；c为岩石内聚力，MPa；
为岩石内摩擦角，°；
(2)将所钻地层的井壁有效主应力带入Mogi-Coulomb准则，得到求解井壁坍塌的非线性函数2；
其中，
式中：
为Mogi-Coulomb准则计算得到的坍塌压力，MPa；a、b为材料常数，取决于岩石内聚力和内摩擦角；
(3)将所钻地层的井壁有效主应力带入Drucker-Prager准则，得到求解井壁坍塌的非线性函数3；
其中，
式中：
为Drucker-Prager准则计算得到的坍塌压力，MPa；α'、k为材料常数，取决于内聚力和内摩擦角；J1为平均有效应力；J2为应力偏量第二不变量；
(4)将所钻地层的井壁有效主应力带入修正Lade准则，得到求解井壁坍塌的非线性函数4；
其中，
式中：
为修正Lade准则计算得到的坍塌压力，MPa；S、η为材料常数，取决于内聚力和内摩擦角；
(5)将所钻地层的井壁有效主应力带入修正Wiebols-Cook准则，得到求解井壁坍塌的非线性函数5；
其中，
式中：
为修正Wiebols-Cook准则计算得到的坍塌压力，MPa；A'、B'、C'为材料常数，取决于内聚力和内摩擦角；UCS为岩石单轴强度，MPa；J1为平均有效应力；J2为应力偏量第二不变量；
S34、利用二分法计算不同强度准则下的坍塌压力非线性函数，得到不同强度准则对应的坍塌压力计算模型；
S4、对所钻地层的N组岩样进行不同岩石强度准则下的坍塌压力计算，并与通过岩石力学试验获得的数据进行拟合，确定符合不同岩石强度准则的岩样个数，再通过计算相应的比值，得到不同岩石强度准则之间的相对重要程度；
S41、通过室内岩石力学试验，对所钻地层的N组岩样进行不同围压下的三轴岩石力学强度测试，然后分别采用Mohr-Coulomb准则、Mogi-Coulomb准则、Drucker-Prager准则、修正Lade准则和修正Wiebols-Cook准则对所钻地层的N组岩样进行坍塌压力计算，将测试结果与计算结果进行拟合；
S42、统计测试结果符合Mohr-Coulomb准则的岩样个数，将其标记为n1；以同样的方法统计测试结果符合Mogi-Coulomb准则、Drucker-Prager准则、修正Lade准则和修正Wiebols-Cook准则的岩样个数，分别标记为n2、n3、n4、n5；
其中，测试结果与计算结果的相对误差在3％以内是试验数据与强度准则吻合较好的标准；
S43、通过计算符合所述岩石强度准则的岩样个数n1、n2、n3、n4和n5之间的比值
得到两两准则之间的相对重要程度aij；
aij为i准则与j准则的重要性比较结果，其值由专家根据
的值，借助比例标度表给出，且有：
表1比例标度表
| 因素i比因素j | 量化值 |
| 同等重要 | 1 |
| 稍微重要 | 3 |
| 较强重要 | 5 |
| 强烈重要 | 7 |
| 极端重要 | 9 |
| 两相邻判断的中间值 | 2，4，6，8 |
S5、根据不同岩石强度准则之间的相对重要程度，并采用层次分析法，建立比较系数矩阵A，对所得权重系数向量进行一致性检验，并确定出不同岩石强度准的权重系数；
S51、建立不同强度准则之间的比较系数矩阵A；
其中，aij为步骤S43中求得的两两准则的相对重要程；度；
S52、计算比较系数矩阵A的特征向量；
首先，将矩阵每一列进行归一化处理：
其次，将归一化后的矩阵按行求和得到列向量；
然后，对求和得到的列向量进行归一化处理，得到特征向量；
S53、对权重系数进行一致性检验；
为了验证权重系数的有效性，通过计算一致性比率(CR)即一致性指数(CI)与随机指数(RI)的比值来检验一致性；定义一致性比率为：
其中，
式中，CR为一致性比率；CI为一致性指数；RI为随机一致性指标；λmax为矩阵A的最大特征根；(AW)i为AW的第i个分量；Wi为权重系数；
所述随机一致性指标RI和比较系数矩阵阶数有关，一般矩阵阶数越大，则出现一致性随机偏离的可能性越大，其对应关系如表2；
表2随机一致性指标RI标准值
| 矩阵阶数 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| RI | 0.00 | 0.00 | 0.58 | 0.90 | 1.12 | 1.24 | 1.32 | 1.41 | 1.45 | 1.49 |
当一致性比率CR<0.1时，认为比较系数矩阵A的不一致程度在允许的范围之内，有满意的一致性，通过一致性检验，可用其归一化特征向量W作为权重系数向量，否则要重新构造比较系数矩阵A，对aij加以调整；
S6、结合不同强度准则计算的坍塌压力和步骤S5确定的权重系数，建立不同强度准则协同作用下定向井坍塌压力计算模型，并计算得到不同强度准则协同作用下定向井井壁坍塌压力；
式中：W1为Mohr-Coulomb准则的权重系数；W2为Mogi-Coulomb准则的权重系数；W3为Drucker-Prager准则的权重系数；W4为修正Lade准则的权重系数；W5为修正Wiebols-Cook准则的权重系数；
为Mohr-Coulomb准则计算得到的坍塌压力，MPa；
为Mogi-Coulomb准则计算得到的坍塌压力，MPa；
为Drucker-Prager准则计算得到的坍塌压力，MPa；
为修正Lade准则计算得到的坍塌压力，MPa；
为修正Wiebols-Cook准则计算得到的坍塌压力，MPa；pc为不同强度准则协同作用下定向井井壁坍塌压力，MPa。
实施例1
本发明的一种不同强度准则协同作用下定向井井壁坍塌压力评价方法，包括：
S1、利用测录井、钻完井、室内测试资料，确定出所钻地层的垂深为3000m、地层孔隙压力为42MPa、上覆岩层压力为90MPa、最大水平地应力为80MPa、最小水平地应力为60MPa、最大水平地应力方位角为0°；
S2、利用测录井和室内岩石力学测试资料，确定出所钻地层的岩石内聚力为25MPa、岩石内摩擦角为35°、岩石泊松比为0.25、岩石孔隙度为0.1、Biot系数为0.8、井壁渗透系数为0；
S3、根据步骤S1和S2所确定的基础参数，建立不同岩石强度准则下的坍塌压力计算模型，利用二分法计算不同强度准则下的坍塌压力非线性函数，得到不同强度准则对应的坍塌压力计算模型
将坍塌压力换算为当量密度，其计算结果如图2～图8所示；
S4、开展室内岩石力学试验，确定符合所述岩石强度准则的岩样个数n1、n2、n3、n4和n5，通过计算相应的比值，得到不同准则之间的相对重要程度；
S5、采用层次分析法，建立比较系数矩阵A，对所得权重系数向量进行一致性检验，并确定出权重系数；
首先，建立不同强度准则之间的比较系数矩阵A；
其次，将矩阵A每一列进行归一化处理，得到：
再将归一化后的矩阵按行求和得到列向量；
然后对求和得到的列向量进行归一化处理，得到特征向量：
W＝[0.22 0.48 0.04 0.06 0.20]T
为了验证权重系数的有效性，通过计算一致性比率(CR)即一致性指数(CI)与随机指数(RI)的比值来检验一致性：
验证后发现，一致性比率CR＜0.1，通过一致性检验；因此，权重系数为W1＝0.22、W2＝0.48、W3＝0.04、W4＝0.06、W5＝0.20；
S6、结合不同强度准则计算的坍塌压力和步骤S5确定的权重系数，建立不同强度准则协同作用下定向井坍塌压力计算模型，并计算得到不同强度准则协同作用下定向井井壁坍塌压力；
于是，将采用Mohr-Coulomb准则、Mogi-Coulomb准则、Drucker-Prager准则、修正Lade准则、修正Wiebols-Cook准则计算得到的结果带入上式，即可得到不同强度准则协同作用下的坍塌压力pc，将坍塌压力换算为当量密度，其计算结果如图2～图8所示。
由图2～图8可以看出：(1)随着井斜角增大，坍塌压力当量密度也随之增大。在不同的方位角下，各准则所计算的坍塌压力当量密度有所差异，其中，Mohr-Coulomb准则由于忽略了中间主应力对岩石破坏的影响，使其计算的坍塌压力明显高于其他计算模型。然后依次是新模型、Mogi-Coulomb准则、修正Lade准则和修正Wiebols-Cook准则，而Drucker-Prager准则在所有情况中对坍塌压力的计算都明显过低。新模型由于引入了权重系数，结合了其他准则的优点，获得了不同强度准则的优势协同作用，只要给出的权重系数合理，就可以计算出更为准确的坍塌压力当量密度值。
以上所述，并非对本发明作任何形式上的限制，虽然本发明已通过上述实施例揭示，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
Claims (9)
1.不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，包括以下步骤：
S1、利用测录井、钻完井、室内测试资料，确定所钻地层的垂深、地层孔隙压力、上覆岩层压力、最大水平地应力、最小水平地应力和最大水平地应力方位角；
S2、利用测录井和室内岩石力学测试资料，确定所钻地层的岩石内聚力、岩石内摩擦角、岩石泊松比、岩石孔隙度、Biot系数和井壁渗透系数；
S3、建立不同岩石强度准则下的坍塌压力计算模型；
S4、对所钻地层的N组岩样进行不同岩石强度准则下的坍塌压力计算，并与通过岩石力学试验获得的数据进行拟合，确定符合不同岩石强度准则的岩样个数，再通过计算相应的比值，得到不同岩石强度准则之间的相对重要程度；
S5、根据不同岩石强度准则之间的相对重要程度，并采用层次分析法，建立比较系数矩阵A，对所得权重系数向量进行一致性检验，并确定出不同岩石强度准的权重系数；
S6、结合不同强度准则计算的坍塌压力和步骤S5确定的权重系数，建立不同强度准则协同作用下定向井坍塌压力计算模型，并计算得到不同强度准则协同作用下定向井井壁坍塌压力。
2.根据权利要求1所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述步骤S3的具体过程为：
S31、计算所钻地层的井壁应力分量；
S32、根据井壁应力分量计算所钻地层的井壁有效主应力；
S33、根据井壁有效主应力建立不同强度准则下的坍塌压力非线性函数；
S34、利用二分法计算不同强度准则下的坍塌压力非线性函数，得到不同强度准则对应的坍塌压力计算模型。
3.根据权利要求2所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述井壁应力分量的计算公式为：
其中，
式中：σrr、σθθ、σzz分别为井壁径向、切向和轴向应力分量，MPa；τθz、τrθ、τrz分别为井壁3个剪应力分量，MPa；pm为井筒压力，MPa；pp为地层孔隙压力，MPa；σv、σH、σh分别为上覆岩石压力、最大水平地应力和最小水平地应力，MPa；δ表示井壁渗透系数，无因次，δ＝0为井壁完全不渗透、δ＝1为井壁完全渗透；φ为岩石孔隙度；ψ为井斜角，°；Ω为井斜方位与最大水平地应力方位夹角，°；θ为井周角，°；v为岩石泊松比；K1为渗流效应系数；α为Biot系数；A、B、C、D、E、F、G、H、J为坐标变换系数。
4.根据权利要求2所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述井壁有效主应力的计算公式为：
σ3＝(1-δφ)pm+δφpp-αpp
式中：σ1为井壁最大主应力，MPa；σ2为井壁中间主应力，MPa；σ3为井壁最小主应力，MPa；α为Biot系数；pp为地层孔隙压力，MPa。
5.根据权利要求2所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述不同强度准则下的坍塌压力非线性函数包括：
井壁坍塌的非线性函数1：
式中：
为Mohr-Coulomb准则计算得到的坍塌压力，MPa；c为岩石内聚力，MPa；
为岩石内摩擦角，°；
井壁坍塌的非线性函数2：
其中，
式中：
为Mogi-Coulomb准则计算得到的坍塌压力，MPa；a、b为材料常数，取决于岩石内聚力和内摩擦角；
井壁坍塌的非线性函数3：
其中，
式中：
为Drucker-Prager准则计算得到的坍塌压力，MPa；α'、k为材料常数，取决于内聚力和内摩擦角；J1为平均有效应力；J2为应力偏量第二不变量；
井壁坍塌的非线性函数4：
其中，
式中：
为修正Lade准则计算得到的坍塌压力，MPa；S、η为材料常数，取决于内聚力和内摩擦角；
井壁坍塌的非线性函数5：
其中，
式中：
为修正Wiebols-Cook准则计算得到的坍塌压力，MPa；A'、B'、C'为材料常数，取决于内聚力和内摩擦角；UCS为岩石单轴强度，MPa；J1为平均有效应力；J2为应力偏量第二不变量。
6.根据权利要求5所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述步骤S4的具体过程为：
S41、对所钻地层的N组岩样进行不同岩石强度准则下的坍塌压力计算得到计算数据，再通过室内岩石力学试验，对所钻地层的N组岩样进行不同围压下的三轴岩石力学强度测试，得到试验数据；并将计算数据与试验数据进行拟合；
S42、分别统计测试结果符合不同岩石强度准则的各个岩样个数，其中计算数据与试验数据的相对误差在3％以内为符合标准；
S43、通过计算符合不同岩石强度准则的各个岩样个数之间的比值，得到两两准则之间的相对重要程度。
7.根据权利要求5所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述不同岩石强度准则包括Mohr-Coulomb准则、Mogi-Coulomb准则、Drucker-Prager准则、修正Lade准则和修正Wiebols-Cook准则。
8.根据权利要求6所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述步骤S5的具体步骤为：
S51、根据不同岩石强度准则之间的相对重要程度建立不同强度准则之间的比较系数矩阵A；
S52、计算比较系数矩阵A的特征向量；
S53、对所得权重系数向量进行一致性检验，并确定出不同岩石强度准的权重系数。
9.根据权利要求1所述的不同强度准则协同作用下定向井井壁坍塌压力评价方法，其特征在于，所述不同强度准则协同作用下定向井坍塌压力计算模型为：
式中，W1为Mohr-Coulomb准则的权重系数；W2为Mogi-Coulomb准则的权重系数；W3为Drucker-Prager准则的权重系数；W4为修正Lade准则的权重系数；W5为修正Wiebols-Cook准则的权重系数；
为Mohr-Coulomb准则计算得到的坍塌压力，MPa；
为Mogi-Coulomb准则计算得到的坍塌压力，MPa；
为Drucker-Prager准则计算得到的坍塌压力，MPa；
为修正Lade准则计算得到的坍塌压力，MPa；
为修正Wiebols-Cook准则计算得到的坍塌压力，MPa；pc为不同强度准则协同作用下定向井井壁坍塌压力，MPa。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210187418.2A CN114526067A (zh) | 2022-02-28 | 2022-02-28 | 不同强度准则协同作用下定向井井壁坍塌压力评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210187418.2A CN114526067A (zh) | 2022-02-28 | 2022-02-28 | 不同强度准则协同作用下定向井井壁坍塌压力评价方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114526067A true CN114526067A (zh) | 2022-05-24 |
Family
ID=81624864
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210187418.2A Pending CN114526067A (zh) | 2022-02-28 | 2022-02-28 | 不同强度准则协同作用下定向井井壁坍塌压力评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114526067A (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104405362A (zh) * | 2014-10-28 | 2015-03-11 | 中国石油集团西部钻探工程有限公司 | 欠平衡钻井井底压力自动控制装置及其使用方法 |
| WO2017024174A1 (en) * | 2015-08-06 | 2017-02-09 | Baker Hughes Incorporated | Downhole fluid typing |
| CA3045009A1 (en) * | 2016-12-07 | 2018-06-14 | Safekick Americas Llc | Automated model-based drilling |
| CN109377101A (zh) * | 2018-11-30 | 2019-02-22 | 西南石油大学 | 一种基于风险控制模型的井壁稳定定量评价方法 |
| CN110821481A (zh) * | 2019-11-15 | 2020-02-21 | 西南石油大学 | 一种空气钻井井壁稳定性评价方法 |
| CN110976525A (zh) * | 2019-12-24 | 2020-04-10 | 太原科技大学 | 一种厚钢板同径异速蛇形差温协同轧制力能参数计算方法 |
| CN111927446A (zh) * | 2020-09-11 | 2020-11-13 | 西南石油大学 | 一种水化页岩地层井壁坍塌失稳预测方法 |
| CN113010835A (zh) * | 2021-03-29 | 2021-06-22 | 重庆交通大学 | 一种基于断裂力学的危岩崩塌预警方法与系统 |
| CN113239577A (zh) * | 2021-06-16 | 2021-08-10 | 西南石油大学 | 一种复杂地层钻井井壁失稳风险定量评价方法 |
| CN113431496A (zh) * | 2021-05-31 | 2021-09-24 | 中国舰船研究设计中心 | 一种钻采船协同作业故障诊断及辅助决策方法 |
| US20210404139A1 (en) * | 2020-06-29 | 2021-12-30 | Zhejiang University Of Technology | System for bridge scour multi-source monitoring, monitoring method thereof, and scour depth evaluating method thereof |
-
2022
- 2022-02-28 CN CN202210187418.2A patent/CN114526067A/zh active Pending
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104405362A (zh) * | 2014-10-28 | 2015-03-11 | 中国石油集团西部钻探工程有限公司 | 欠平衡钻井井底压力自动控制装置及其使用方法 |
| WO2017024174A1 (en) * | 2015-08-06 | 2017-02-09 | Baker Hughes Incorporated | Downhole fluid typing |
| CA3045009A1 (en) * | 2016-12-07 | 2018-06-14 | Safekick Americas Llc | Automated model-based drilling |
| CN109377101A (zh) * | 2018-11-30 | 2019-02-22 | 西南石油大学 | 一种基于风险控制模型的井壁稳定定量评价方法 |
| CN110821481A (zh) * | 2019-11-15 | 2020-02-21 | 西南石油大学 | 一种空气钻井井壁稳定性评价方法 |
| CN110976525A (zh) * | 2019-12-24 | 2020-04-10 | 太原科技大学 | 一种厚钢板同径异速蛇形差温协同轧制力能参数计算方法 |
| US20210404139A1 (en) * | 2020-06-29 | 2021-12-30 | Zhejiang University Of Technology | System for bridge scour multi-source monitoring, monitoring method thereof, and scour depth evaluating method thereof |
| CN111927446A (zh) * | 2020-09-11 | 2020-11-13 | 西南石油大学 | 一种水化页岩地层井壁坍塌失稳预测方法 |
| CN113010835A (zh) * | 2021-03-29 | 2021-06-22 | 重庆交通大学 | 一种基于断裂力学的危岩崩塌预警方法与系统 |
| CN113431496A (zh) * | 2021-05-31 | 2021-09-24 | 中国舰船研究设计中心 | 一种钻采船协同作业故障诊断及辅助决策方法 |
| CN113239577A (zh) * | 2021-06-16 | 2021-08-10 | 西南石油大学 | 一种复杂地层钻井井壁失稳风险定量评价方法 |
Non-Patent Citations (3)
| Title |
|---|
| LISONG ZHANG: "A New Analytical Model to Evaluate Uncertainty of Wellbore Collapse Pressure Based on Advantageous Synergies of Different Strength Criteria", ROCK MECHANICS AND ROCK ENGINEERING, 23 January 2019 (2019-01-23), pages 2649, XP036872399, DOI: 10.1007/s00603-018-1729-1 * |
| 程军: "深井钻井工程风险评价方法及应用", 工程建设与设计, no. 7, 15 July 2016 (2016-07-15), pages 8 - 66 * |
| 马天寿;张赟;邱艺;刘阳;陈平: "基于可靠度理论的斜井井壁失稳风险评价方法", 石油学报, vol. 42, no. 011, 30 November 2021 (2021-11-30), pages 1486 - 1498 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | In-situ stress measurements and stress distribution characteristics of coal reservoirs in major coalfields in China: Implication for coalbed methane (CBM) development | |
| CN109858147B (zh) | 一种基于可靠度理论的井壁失稳风险定量评价方法 | |
| CN106855897B (zh) | 适用于压力衰竭地层的井壁稳定的研究方法 | |
| CN105893679B (zh) | 低产水平井续流修正试井解释方法 | |
| CN110595887B (zh) | 一种浅埋硬岩大型地下洞室成洞条件量化评价方法 | |
| CN112727534B (zh) | 基于真三向应力和渗透率动态变化的瓦斯钻孔布置方法 | |
| CN114372352A (zh) | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 | |
| CN109033504A (zh) | 一种油水井套管损坏预测方法 | |
| CN113109228B (zh) | 基于煤矿地质钻孔的松散承压含水层渗透系数确定方法 | |
| CN113109162B (zh) | 一种基于热流固耦合的岩石起裂压力计算方法 | |
| CN112412434B (zh) | 一种改进的疏松砂岩地应力计算方法 | |
| CN116205163B (zh) | 一种天然气水合物藏数值模拟中的地层稳定性分析方法 | |
| CN113123783B (zh) | 一种精细评价地层体积含量的方法 | |
| CN113175322A (zh) | 地层漏失压力剖面的建立方法、计算机装置及存储介质 | |
| CN114526067A (zh) | 不同强度准则协同作用下定向井井壁坍塌压力评价方法 | |
| CN112983406B (zh) | 一种天然气水合物储层参数指标评价方法 | |
| CN116658157A (zh) | 致密砂岩气藏地层压力预测方法及系统 | |
| CN111241460A (zh) | 一种复杂致密储层孔隙度计算方法 | |
| CN112561223B (zh) | 油气井增产潜力的定量评价方法 | |
| CN110671102B (zh) | 一种气井临界出砂压差的确定方法及系统 | |
| CN113554311A (zh) | 地面定向孔注浆治理推覆体下奥灰水害工程质量评价方法 | |
| CN110792425B (zh) | 一种测量地层流体含水率的方法 | |
| CN114647913A (zh) | 一种基于Mohr-Coulomb准则的最大水平主地应力获取方法 | |
| Wei et al. | Application of geological model dynamic updating method in Guanjiapu Oilfield, Dagang | |
| CN113123782B (zh) | 一种利用线性规划评价地层组分的方法 |
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 |