CN115711120A - 一种油气井水泥环完整性及损伤演化定量评价方法 - Google Patents
一种油气井水泥环完整性及损伤演化定量评价方法 Download PDFInfo
- Publication number
- CN115711120A CN115711120A CN202211542803.0A CN202211542803A CN115711120A CN 115711120 A CN115711120 A CN 115711120A CN 202211542803 A CN202211542803 A CN 202211542803A CN 115711120 A CN115711120 A CN 115711120A
- Authority
- CN
- China
- Prior art keywords
- cement
- cement sheath
- fractal dimension
- blank
- control group
- 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
Links
- 239000004568 cement Substances 0.000 title claims abstract description 286
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 14
- 239000011148 porous material Substances 0.000 claims abstract description 54
- 239000002245 particle Substances 0.000 claims abstract description 39
- 230000009471 action Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 22
- 239000004575 stone Substances 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 20
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 20
- 229910052753 mercury Inorganic materials 0.000 claims description 20
- 238000002347 injection Methods 0.000 claims description 14
- 239000007924 injection Substances 0.000 claims description 14
- 238000001878 scanning electron micrograph Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims description 5
- 239000002002 slurry Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001143 conditioned effect Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 230000033228 biological regulation Effects 0.000 claims description 2
- 239000007789 gas Substances 0.000 abstract description 20
- 238000005553 drilling Methods 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 239000003345 natural gas Substances 0.000 abstract description 2
- 239000003209 petroleum derivative Substances 0.000 abstract description 2
- 238000011156 evaluation Methods 0.000 description 8
- 230000005465 channeling Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 239000010410 layer Substances 0.000 description 4
- 230000000877 morphologic effect Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000007789 sealing Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000020477 pH reduction Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 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
- E21B47/00—Survey of boreholes or wells
- E21B47/005—Monitoring or checking of cementation quality or level
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Quality & Reliability (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种油气井水泥环完整性及损伤演化定量评价方法，其特征在于，主要基于分形理论、图像处理技术及套管‑水泥环‑地层组合体结构特征及失效机理，采用初始状态下的空白组套管‑水泥环界面形态分形维数、水泥环本体微观孔隙形态分形维数、颗粒形态分形维数及其与水泥环界面径向胶结强度、抗拉强度和抗压强度等宏观力学性能之间的相关性来定量评价空白组水泥环完整性；采用交变载荷作用后的条件对照组套管‑水泥环界面形态分形维数、水泥环本体微观孔隙形态分形维数、颗粒形态分形维数、裂纹形态分形维数及其与水泥环宏观力学性能之间的相关性来定量评价条件对照组水泥环完整性损伤演化规律。本发明适用于石油与天然气钻采工程技术领域。
Description
技术领域
本发明涉及深井、超深井及深部页岩气水平井固井水泥环完整性检测及评价领域，具体是一种油气井水泥环完整性及损伤演化定量评价方法。
背景技术
在石油天然气的钻探过程中，固井的目的是为封隔地层、加固井眼、建立密封性能良好的井内油气水流动通道，以保障继续高效安全钻进，并保证后期试油、注气/水、压裂、酸化、生产等作业的正常进行。为实现上述固井目的，则需介于地层与套管外壁之间的水泥环长期保持良好的完整性。水泥环完整性是指在油气井服役期间，水泥环保持良好的结构完整性和功能完整性。
由于深部及非常规低孔、低渗油气藏通常需要大规模水力压裂才能获取商业油气产量。然而，在大规模多级水力压裂过程中，因温度和压力的波动、连续变化、交替加载和卸载很容易导致套管-水泥环-地层组合体完整性失效。水泥环完整性失效容易引起层间封隔失效、环空带压、套管腐蚀及服役寿命降低、管外油、气、水容窜通等诸多复杂工况。
调研发现：水泥环完整性失效主要包括界面失效和水泥环本体失效两大类，其中界面失效分为第一界面分离、第二界面分离、轴向滑移，而水泥环本体失效分为环向/轴向开裂、屈服破坏。因此，为预防并控制水泥环完整性失效，迫切需要模拟真实井下温度、压力及作业特征，准确而高效评价水泥环的完整性，进而基于评价结果，及时采取具有针对性地预防及控制措施，保证工程的质量与安全，避免类似事故及风险发生，这对现场安全、高效的生产显得尤为重要。目前，国内外学者基于理论和实验的方法，在水泥环完整性失效机理及水泥环完整性评价等方面展开了诸多的研究，主要包括套管-水泥环界面胶结强度、密封性能及完整性失效形式。上述研究主要基于气窜原理，即采用水泥环内部气体连通这种破坏行为作为判据，并采用气窜评价仪(如FMA7150)开展水泥环完整性评价；分析发现，现有评价仪器及方法均只能定性评价已形成大规模沟通微环隙、裂缝、窜槽的水泥环完整性，而不能用于定量评价无沟通微环隙、裂缝、窜槽的水泥环完整性，且因评价过程中水泥环内部气体连通这种破坏行为不能保证水泥环完整性，从而影响评价结果的准确性。相反，在保持水泥环完整性的前提下，有关定量评价模拟真实工况下水泥环完整性及损伤演化特征的方法及相关研究鲜见报道。
近年来，随着深井、超深井及非常规水平井的开发，钻完井数量不断增加，油气井具有超深、超高压、高温的特点，水泥环完整性问题越发突出。为了保证水泥环在油气井服役期间保持良好的结构完整性和功能完整性，预防并控制水泥环完整性失效，迫切需要建立一种定量评价模拟工况下水泥环完整性的新方法。
因此，本发明针对目前在模拟工况下水泥环完整性难以定量评价的技术难题，基于分形理论，提出一种油气井水泥环完整性及损伤演化定量评价方法，该方法可准确获取真实工况下水泥环界面、水泥环本体孔隙、颗粒、裂纹形态的分形维数及其与宏观力学性能的相关性，并基于此评价水泥环完整性及其在交变载荷下的损伤演化规律，可为油气井固井水泥浆体系的优化设计及水泥环完整性定量评价提供支撑。
发明内容
本发明的目的是提供一种油气井水泥环完整性及损伤演化定量评价方法，以解决实际模拟工况下水泥环完整性定量评价的技术难题，该方法简单易行，且完全适用于深井、超深井及非常规油气井水泥环在实际服役工况下的完整性定量评价。
为实现上述目的，本发明所述的一种油气井水泥环完整性及损伤演化定量评价方法，其特征在于：该方法主要基于分形理论、图像处理技术及套管-水泥环-地层组合体结构特征及失效机理，采用以下原理定量评价套管-水泥环-地层组合体中水泥环完整性及损伤演化规律：(1)采用空白组水泥环界面(包括与水泥环接触的套管外壁及与套管接触的水泥环内壁)形态分形维数、水泥环本体孔隙形态分形维数、颗粒形态分形维数及其与宏观力学性能(水泥环界面径向胶结强度、抗拉强度和抗压强度)之间的相关性来定量评价空白组水泥环完整性；(2)采用条件对照组水泥环界面形态分形维数、水泥环本体孔隙形态分形维数、颗粒形态分形维数、裂纹形态分形维数及其与宏观力学性能之间的相关性来定量评价条件对照组水泥环完整性损伤演化规律，具体采取的技术方案如下：
步骤一：制备用于开展水泥环完整性及损伤演化定量评价实验样品：(1)采用现场井身结构和水泥浆体系，模拟实际温度及压力养护制备套管-水泥环-地层组合体，定义与水泥环接触的套管外壁为目标面，与套管接触的水泥环内壁为接触面，目标面和接触面为套管-水泥环界面；(2)将制备的套管-水泥环-地层组合体分为空白组和条件对照组，空白组用于定量评价水泥环完整性，条件对照组用于定量评价交变载荷下水泥环完整性损伤演化规律，需要说明的是条件对照组为交变载荷作用后得到的套管-水泥环-地层组合体，空白组为初始状态的套管-水泥环-地层组合体；(3)分离空白组的水泥环和套管，采用空白组的水泥环和套管制备水泥石力学性能测试试样、目标面及接触面三维轮廓扫描试样、水泥石SEM扫描试样、水泥石压汞测试试样；(4)分离条件对照组的水泥环和套管，采用条件对照组的水泥环和套管制备水泥石力学试样、目标面及接触面的三维轮廓试样、水泥石SEM试样、水泥石压汞试样。
步骤二：测定有无交变载荷作用的水泥环宏观力学性能：(1)采用步骤一制备的空白组套管-水泥环-地层组合体开展界面力学性能测试，获取空白组套管-水泥环界面径向胶结强度SBR；(2)采用步骤一制备的空白组水泥石力学试样开展力学性能测试，获取空白组水泥石抗拉强度QBL和抗压强度QBC；(3)采用步骤一制备的条件对照组力学水泥石试样开展力学性能测试，获取条件对照组水泥石抗拉强度QEL和抗压强度QEC；(4)采用步骤一制备的条件对照组套管-水泥环-地层组合体开展界面力学性能测试，获取条件对照组套管-水泥环界面径向胶结强度SER；
步骤三：测定及评价无交变载荷作用的套管-水泥环界面形态的分形维数：(1)采用光学衍射仪对步骤一得到的空白组及条件对照组目标面和接触面三维轮廓试样进行三维扫描，获取目标面和接触面的三维轮廓图及不同测量尺度τTF和τCF下所有测点高度H，其中τTF为目标面的测量尺度，τCF为接触面的测量尺度；(2)采用基于结构函数的分形模型LgS(τTF)＝lgCTF+(4-2DTFθ)LgτTF将目标面的测量尺度τTF及对应的结构测度函数S(τTF)绘制在双对数坐标，其中S(τTF)为目标面的结构测度函数，S(τTF)＝[H(Z+τTF，θ)-H(Z，θ)]2，Z为目标面上测点数据沿套管轴向的坐标；θ为目标面上测点数据沿套管环向的坐标；τTF为目标面的测量尺度，即两测点数据之间的轴向间隔；H(Z+τCF，θ)为目标面三维轮廓图中测点(Z+τCF，θ)的高度；H(Z,θ)为目标面三维轮廓图中测点(Z,θ)的高度；CTF为目标面的尺度系数；DTFθ为目标面沿不同角度θ的分形维数，θ为0°、90°、180°、270°；(3)采用基于结构函数的分形模型LgS(τCF)＝lgCCF+(4-2DCFα)LgτCF将接触面的测量尺度τCF及对应的结构测度函数S(τCF)绘制在双对数坐标，其中S(τCF)为接触面的结构测度函数，S(τCF)＝[H(Y+τTF，α)-H(Y，α)]2，Y为接触面上测点数据沿水泥环轴向的坐标；α为接触面上测点数据沿水泥环环向的坐标；τCF为接触面的测量尺度，即两测点数据之间的轴向间隔；H(Y+τCF，α)为接触面三维轮廓图中测点(Y+τCF，α)的高度；H(Y,α)为接触面三维轮廓图中测点(Y,θ)的高度；CCF为接触面的尺度系数；DCFα为接触面沿不同角度α的分形维数，α为0°、90°、180°、270°；(4)通过LgS(τTF)＝lgCTF+(4-2DTFθ)LgτTF曲线斜率计算空白组和条件对照组目标面沿0°、90°、180°、270°四个方向的分形维数，并取DTF0、DTF90、DTF270、和DTF360的平均值做为接触面的分形维数，记DBTF为空白组目标面的分形维数，DETF为条件对照组目标面的分形维数；(5)通过LgS(τCF)＝lgCCF+(4-2DCFα)LgτCF曲线斜率计算空白组和条件对照组接触面沿0°、90°、180°、270°四个方向的分形维数，并取DCF0、DCF90、DCF270和DCF360的平均值做为接触面的分形维数，令DBCF为空白组接触面的分形维数，DECF为条件对照组接触面的分形维数。
步骤四：测定及评价有无交变载荷作用的水泥环本体孔隙形态的分形维数：(1)采用压汞法对步骤一得到的空白组及条件对照组水泥石压汞试样进行压汞测试，获取空白组及条件对照组水泥石的真实孔隙率φ、不同注入压力Pi下进入水泥石孔隙的总体积VPi及对应压力Pi下的孔隙直径2Ri；(2)采用孔体积分形模型Lg(|dVPi/dRi|)＝(2-DP)LgRi+CP将总体积VPi随孔隙半径增量的绝对值|dVPi/dRi|与孔隙半径绘制在双对数坐标，其中VPi为不同注入压力Pi下进入水泥环试样中孔隙的总水银体积；Ri为对应注入压力Pi下对应的孔隙半径；DP为水泥环试样中孔隙形态的分形维数；CP为水泥环试样孔隙的分形模型常数；(3)通过Lg(|dVPi/dRi|)＝(2-DP)LgRi+CP曲线斜率计算空白组和条件对照组水泥环孔隙形态的分形维数，令DBP为空白组水泥环孔隙形态的分形维数，DEP为条件对照组水泥环孔隙形态的分形维数。
步骤五：测定及评价有无交变载荷作用的水泥环本体颗粒形态的分形维数：(1)采用扫描电镜对步骤一得到的空白组及条件对照组水泥石SEM试样进行面扫描，获取空白组及条件对照组水泥石在不同放大倍率下的SEM图像；(2)采用Python+OpenCV程序对SEM图像进行二值化处理，得到不同阀值下的二值化图像，白色区域代表微观颗粒，黑色区域代表微观孔隙；(3)基于步骤四采用压汞法获取的水泥石试样真实孔隙率φ，以水泥石试样的真实孔隙率为控制对象，采用基于边缘强度的阀值分割算法，对二值化图像阀值进行适配调控，直到二值化图像的三维孔隙率等于真实孔隙率，确定最佳阀值，从而优选出与水泥石试样真实孔隙率相等的目标二值化图像；(4)采用Matlab编程计算目标二值化图像中所有白色区域的面积和周长；(5)采用面积-周长分形模型Lg(AGi)＝DG*Lg(PGi)+CG将目标二值化图像中白色区域的面积和周长绘制在双对数坐标，其中PGi为目标二值化图像中任一个白色几何图形的等效周长，AGi为目标二值化图像中等效周长为PGi的几何图形等效面积，DG为水泥石试样颗粒形态的分形维数，CG为水泥环试样颗粒的分形模型常数；(6)通过Lg(AGi)＝DG*Lg(PGi)+CG曲线斜率计算空白组和条件对照组水泥环颗粒形态的分形维数，令DBG为空白组水泥环颗粒形态的分形维数，DEG为条件对照组水泥环颗粒形态的分形维数。
步骤六：测定及评价交变载荷作用后水泥环本体裂纹形态的分形维数：(1)采用扫描电镜对步骤一得到的条件对照组水泥石SEM试样进行面扫描，获取条件对照组水泥石在不同放大倍率下的SEM图像；(2)采用与步骤五相同的方法获取水泥石裂纹的目标二值化图像；(3)采用Matlab编程计算用不同边长δFi的正方形盒子覆盖目标二值化图像的盒子总数N(δFi)；(4)采用盒子模型LgN(δFi)＝DEF*LgδFi+CF将正方形盒子边长δFi和盒子总数N(δFi)绘制在双对数坐标，其中δFi为测量目标二值化图像用正方形盒子的边长，N(δFi)为用边长为δFi的盒子去覆盖目标二值化图像所需要的盒子总数；DEF为水泥环本体裂纹形态的分形维数；CEF为水泥环本体裂纹形态的分形模型常数；(5)通过LgN(δFi)＝DEF*LgδFi+CEF曲线斜率计算条件对照组水泥环本体裂纹形态的分形维数DEF。
步骤七：基于步骤二和步骤三，建立空白组界面径向胶结强度与套管水泥环界面形态分形维数之间的函数关系FB1(SBR，DBTF)和FB2(SBR，DBCF)。
步骤八：基于步骤二和步骤四，建立空白组水泥石力学性能与套管水泥环本体孔隙形态分形维数之间的函数关系FB3(QBL，DBP)和FB4(QBC，DBP)。
步骤九：基于步骤二和步骤五，建立空白组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FB5(QBL，DBG)和FB6(QBC，DBG)。
步骤十：基于步骤二和步骤三，建立条件对照组界面径向胶结强度与套管水泥环界面形态分形维数之间的函数关系FE1(SER，DETF)和FE2(SER，DECF)。
步骤十一：基于步骤二和步骤四，建立条件对照组水泥石力学性能与套管水泥环本体孔隙形态分形维数之间的函数关系FE3(QEL，DEP)和FE4(QEC，DEP)。
步骤十二：基于步骤二和步骤五，建立条件对照组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FE5(QEL，DEG)和FE6(QEC，DEG)。
步骤十三：基于步骤二和步骤六，建立条件对照组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FE7(QEL，DEF)和FE8(QEC，DEF)。
步骤十四：采用分形维数DBTF、DBCF、DBP、DBG和函数关系FB1(SBR，DBTF)、FB2(SBR，DBCF)、FB3(QBL，DBP)、FB4(QBC，DBP)、FB5(QBL，DBG)及FB6(QBC，DBG)定量评价空白组水泥环完整性；
步骤十五：采用分形维数DETF、DECF、DEP、DEG、DEF和函数关系FE1(SER，DETF)、FE2(SER，DECF)、FE3(QEL，DEP)、FE4(QEC，DEP)FE5(QEL，DEG)、FE6(QEC，DEG)、FE7(QEL，DEF)和FE8(QEC，DEF)定量评价条件对照组水泥环完整性损伤演化规律。
本发明的优点是：
基于分析理论、图形处理技术及水泥环完整性失效机理，准确构建水泥环完整性与其微观结构形态之间的相关性，实现模拟工况下套管-水泥环-地层组合体中水泥环完整性及损伤演化规律的定量评价，且能更好的预防和控制水泥环完整性失效，同时，为油气井固井水泥浆体系、施工工艺参数的优化设计及水泥环完整性评价提供一种新方法。
附图说明
图1为空白组水泥环本体孔隙Lg|dVPi/dRi|与LgRi之间的关系曲线
图2为空白组水泥环本体颗粒Lg(AGi)与Lg(PGi)之间的关系曲线
图3为15次交变温度作用后水泥环本体孔隙Lg|dVPi/dRi|与LgRi之间的关系曲线
图4为15次交变温度作用后水泥环本体颗粒Lg(AGi)与Lg(PGi)之间的关系曲线
图5为不同次数的交变温度载荷作用后水泥环本体孔隙及颗粒形态的分形维数
具体实施方式
下面结合附图及实施例，对本发明进行详细的描述。
参见附图，本发明提出了一种油气井水泥环完整性及损伤演化定量评价方法，所述方法主要包括以下步骤：
步骤一：制备用于开展水泥环完整性及损伤演化定量评价实验样品：(1)采用我国西部油田某井井身结构和水泥浆体系，在120℃*35MPa养护条件下制备套管-水泥环-地层组合体两组(甲组和乙组)，甲组为空白组，乙组分别置于温度30℃→120℃→30℃下交变作用10次、15次、20次、25次得到条件对照组；(2)分离空白组的水泥环和套管，采用空白组的水泥环和套管制备水泥石力学性能测试试样、目标面及接触面三维轮廓扫描试样、水泥石SEM扫描试样、水泥石压汞测试试样；(4)分离条件对照组的水泥环和套管，采用条件对照组的水泥环和套管制备水泥石力学试样、目标面及接触面的三维轮廓试样、水泥石SEM试样、水泥石压汞试样。
步骤二：测定有无交变温度(30℃→120℃→30℃)作用的水泥环宏观力学性能：(1)采用步骤一制备的空白组套管-水泥环-地层组合体开展界面力学性能测试，获取空白组套管-水泥环界面径向胶结强度SBR；(2)采用步骤一制备的空白组水泥石力学试样开展力学性能测试，获取空白组水泥石抗拉强度QBL和抗压强度QBC；(3)采用步骤一制备的条件对照组力学水泥石试样开展力学性能测试，获取条件对照组水泥石抗拉强度QEL和抗压强度QEC；(4)采用步骤一制备的条件对照组套管-水泥环-地层组合体开展界面力学性能测试，获取条件对照组套管-水泥环界面径向胶结强度SER；
步骤三：测定及评价有无交变温度(30℃→120℃→30℃)作用的套管-水泥环界面形态的分形维数：(1)采用光学衍射仪对步骤一得到的空白组及条件对照组目标面和接触面三维轮廓试样进行三维扫描，获取目标面和接触面的三维轮廓图及不同测量尺度τTF和τCF下所有测点高度H，其中τTF为目标面的测量尺度，τCF为接触面的测量尺度；(2)采用基于结构函数的分形模型LgS(τTF)＝lgCTF+(4-2DTFθ)LgτTF将目标面的测量尺度τTF及对应的结构测度函数S(τTF)绘制在双对数坐标；(3)同理，采用基于结构函数的分形模型LgS(τCF)＝lgCCF+(4-2DCFα)LgτCF将接触面的测量尺度τCF及对应的结构测度函数S(τCF)绘制在双对数坐标；(4)通过LgS(τTF)＝lgCTF+(4-2DTFθ)LgτTF曲线斜率计算空白组和条件对照组目标面沿0°、90°、180°、270°四个方向的分形维数，并取DTF0、DTF90、DTF270、和DTF360的平均值做为接触面的分形维数，记DBTF为空白组目标面的分形维数，DETF为条件对照组目标面的分形维数；(5)通过LgS(τCF)＝lgCCF+(4-2DCFα)LgτCF曲线斜率计算空白组和条件对照组接触面沿0°、90°、180°、270°四个方向的分形维数，并取DCF0、DCF90、DCF270和DCF360的平均值做为接触面的分形维数，令DBCF为空白组接触面的分形维数，DECF为条件对照组接触面的分形维数。
步骤四：测定及评价有无交变温度(30℃→120℃→30℃)作用的水泥环本体孔隙形态的分形维数：(1)采用压汞法对步骤一得到的空白组及条件对照组水泥石压汞试样进行压汞测试，获取空白组及条件对照组水泥石的真实孔隙率φ、不同注入压力Pi下进入水泥石孔隙的总体积VPi及对应压力Pi下的孔隙直径2Ri；(2)采用孔体积分形模型Lg(|dVPi/dRi|)＝(2-DP)LgRi+CP将总体积VPi随孔隙半径增量的绝对值|dVPi/dRi|与孔隙半径绘制在双对数坐标；(3)通过Lg(|dVPi/dRi|)＝(2-DP)LgRi+CP曲线斜率计算空白组和条件对照组水泥环孔隙形态的分形维数，附图1显示空白组水泥环孔隙形态的分形维数DBP为3.3466，附图3显示条件对照组水泥环经15次交变温度作用后孔隙形态的分形维数DEP为3.6703。
步骤五：测定及评价有无交变温度(30℃→120℃→30℃)作用的水泥环本体颗粒形态的分形维数：(1)采用扫描电镜对步骤一得到的空白组及条件对照组水泥石SEM试样进行面扫描，获取空白组及条件对照组水泥石在不同放大倍率下的SEM图像；(2)采用Python+OpenCV程序对SEM图像进行二值化处理，得到不同阀值下的二值化图像，白色区域代表微观颗粒，黑色区域代表微观孔隙；(3)基于步骤四采用压汞法获取的水泥石试样真实孔隙率φ，以水泥石试样的真实孔隙率为控制对象，采用基于边缘强度的阀值分割算法，对二值化图像阀值进行适配调控，优选出与水泥石试样真实孔隙率相等的目标二值化图像；(4)采用Matlab编程计算目标二值化图像中所有白色区域的面积和周长；(5)采用面积-周长分形模型Lg(AGi)＝DG*Lg(PGi)+CG将目标二值化图像中白色区域的面积和周长绘制在双对数坐标；(6)通过Lg(AGi)＝DG*Lg(PGi)+CG曲线斜率计算空白组和条件对照组水泥环颗粒形态的分形维数，附图2显示空白组水泥环颗粒形态的分形维数DBG为1.5323，附图4显示条件对照组水泥环经15次交变温度作用后颗粒形态的分形维数DEG为1.3926。
步骤六：测定及评价交变温度(30℃→120℃→30℃)作用后水泥环本体裂纹形态的分形维数：(1)采用扫描电镜对步骤一得到的条件对照组水泥石SEM试样进行面扫描，获取条件对照组水泥石在不同放大倍率下的SEM图像；(2)采用与步骤五相同的方法获取水泥石裂纹的目标二值化图像；(3)采用Matlab编程计算用不同边长δFi的正方形盒子覆盖目标二值化图像的盒子总数N(δFi)；(4)采用盒子模型LgN(δFi)＝DEF*LgδFi+CF将正方形盒子边长δFi和盒子总数N(δFi)绘制在双对数坐标；(5)通过LgN(δFi)＝DEF*LgδFi+CEF曲线斜率计算条件对照组水泥环本体裂纹形态的分形维数DEF。
步骤七：基于步骤二和步骤三，建立空白组界面径向胶结强度与套管水泥环界面形态分形维数之间的函数关系FB1(SBR，DBTF)和FB2(SBR，DBCF)。
步骤八：基于步骤二和步骤四，建立空白组水泥石力学性能与套管水泥环本体孔隙形态分形维数之间的函数关系FB3(QBL，DBP)和FB4(QBC，DBP)。
步骤九：基于步骤二和步骤五，建立空白组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FB5(QBL，DBG)和FB6(QBC，DBG)。
步骤十：基于步骤二和步骤三，建立条件对照组界面径向胶结强度与套管水泥环界面形态分形维数之间的函数关系FE1(SER，DETF)和FE2(SER，DECF)。
步骤十一：基于步骤二和步骤四，建立条件对照组水泥石力学性能与套管水泥环本体孔隙形态分形维数之间的函数关系FE3(QEL，DEP)和FE4(QEC，DEP)。
步骤十二：基于步骤二和步骤五，建立条件对照组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FE5(QEL，DEG)和FE6(QEC，DEG)。
步骤十三：基于步骤二和步骤六，建立条件对照组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FE7(QEL，DEF)和FE8(QEC，DEF)。
步骤十四：采用分形维数DBTF、DBCF、DBP、DBG和函数关系FB1(SBR，DBTF)、FB2(SBR，DBCF)、FB3(QBL，DBP)、FB4(QBC，DBP)、FB5(QBL，DBG)及FB6(QBC，DBG)定量评价空白组水泥环完整性；
步骤十五：采用分形维数DETF、DECF、DEP、DEG、DEF和函数关系FE1(SER，DETF)、FE2(SER，DECF)、FE3(QEL，DEP)、FE4(QEC，DEP)FE5(QEL，DEG)、FE6(QEC，DEG)、FE7(QEL，DEF)和FE8(QEC，DEF)定量评价条件对照组水泥环完整性损伤演化规律。
附图5显示套管-水泥环界面径向胶结强度与水泥石微观孔隙及颗粒形态分形维数密切相关，即随交变温度作用次数的增加，套管-水泥环界面径向胶结强度SR逐渐下降，导致水泥环更倾向发生界面失效并形成微环隙，对应的孔隙分形维数DP增加，主要是交变温度作用导致孔隙率增加、孔隙结构劣化且分布变得更加复杂，而对应的颗粒形态分形维数DG降低，主要交变温度作用导致颗粒排列结构松散，上述两方面的微观结构变化导致水泥环宏观力学性能劣化。
Claims (1)
1.一种油气井水泥环完整性及损伤演化定量评价方法，其特征在于，所述的油气井水泥环完整性及损伤演化定量评价方法包括以下步骤：
步骤一：制备用于开展水泥环完整性及损伤演化定量评价实验样品：(1)采用现场井身结构和水泥浆体系，模拟实际温度及压力养护制备套管-水泥环-地层组合体，定义与水泥环接触的套管外壁为目标面，与套管接触的水泥环内壁为接触面；(2)将制备的套管-水泥环-地层组合体分为空白组和条件对照组，空白组用于定量评价水泥环完整性，条件对照组用于定量评价交变载荷下水泥环完整性损伤演化规律；(3)分离空白组的水泥环和套管，采用空白组的水泥环和套管制备水泥石力学性能测试试样、目标面及接触面三维轮廓扫描试样、水泥石SEM扫描试样、水泥石压汞测试试样；(4)分离条件对照组的水泥环和套管，采用条件对照组的水泥环和套管制备水泥石力学试样、目标面及接触面的三维轮廓试样、水泥石SEM试样、水泥石压汞试样；
步骤二：测定有无交变载荷作用的水泥环宏观力学性能：(1)采用步骤一制备的空白组套管-水泥环-地层组合体开展界面力学性能测试，获取空白组套管-水泥环界面径向胶结强度SBR；(2)采用步骤一制备的空白组水泥石力学试样开展力学性能测试，获取空白组水泥石抗拉强度QBL和抗压强度QBC；(3)采用步骤一制备的条件对照组力学水泥石试样开展力学性能测试，获取条件对照组水泥石抗拉强度QEL和抗压强度QEC；(4)采用步骤一制备的条件对照组套管-水泥环-地层组合体开展界面力学性能测试，获取条件对照组套管-水泥环界面径向胶结强度SER；
步骤三：测定及评价有无交变载荷作用的套管-水泥环界面形态的分形维数：(1)采用光学衍射仪对步骤一得到的空白组及条件对照组目标面和接触面三维轮廓试样进行三维扫描，获取目标面和接触面的三维轮廓图及不同测量尺度τTF和τCF下所有测点高度H，其中τTF为目标面的测量尺度，τCF为接触面的测量尺度；(2)采用基于结构函数的分形模型LgS(τTF)＝lgCTF+(4-2DTFθ)LgτTF将目标面的测量尺度τTF及对应的结构测度函数S(τTF)绘制在双对数坐标，其中S(τTF)为目标面的结构测度函数，S(τTF)＝[H(Z+τTF，θ)-H(Z，θ)]2，Z为目标面上测点数据沿套管轴向的坐标；θ为目标面上测点数据沿套管环向的坐标；τTF为目标面的测量尺度，即两测点数据之间的轴向间隔；H(Z+τCF，θ)为目标面三维轮廓图中测点(Z+τCF，θ)的高度；H(Z,θ)为目标面三维轮廓图中测点(Z,θ)的高度；CTF为目标面的尺度系数；DTFθ为目标面沿不同角度θ的分形维数；(3)采用基于结构函数的分形模型LgS(τCF)＝lgCCF+(4-2DCFα)LgτCF将接触面的测量尺度τCF及对应的结构测度函数S(τCF)绘制在双对数坐标，其中S(τCF)为接触面的结构测度函数，S(τCF)＝[H(Y+τTF，α)-H(Y，α)]2，Y为接触面上测点数据沿水泥环轴向的坐标；α为接触面上测点数据沿水泥环环向的坐标；τCF为接触面的测量尺度，即两测点数据之间的轴向间隔；H(Y+τCF，α)为接触面三维轮廓图中测点(Y+τCF，α)的高度；H(Y,α)为接触面三维轮廓图中测点(Y,θ)的高度；CCF为接触面的尺度系数；DCFα为接触面沿不同角度α的分形维数；(4)通过LgS(τTF)＝lgCTF+(4-2DTFθ)LgτTF曲线斜率计算空白组和条件对照组目标面沿0°、90°、180°、270°四个方向的分形维数，并取DTF0、DTF90、DTF270、和DTF360的平均值做为接触面的分形维数，记DBTF为空白组目标面的分形维数，DETF为条件对照组目标面的分形维数；(5)通过LgS(τCF)＝lgCCF+(4-2DCFα)LgτCF曲线斜率计算空白组和条件对照组接触面沿0°、90°、180°、270°四个方向的分形维数，并取DCF0、DCF90、DCF270和DCF360的平均值做为接触面的分形维数，令DBCF为空白组接触面的分形维数，DECF为条件对照组接触面的分形维数；
步骤四：测定及评价有无交变载荷作用的水泥环本体孔隙形态的分形维数：(1)采用压汞法对步骤一得到的空白组及条件对照组水泥石压汞试样进行压汞测试，获取空白组及条件对照组水泥石的真实孔隙率
不同注入压力Pi下进入水泥石孔隙的总体积VPi及对应压力Pi下的孔隙直径2Ri；(2)采用孔体积分形模型Lg(|dVPi/dRi|)＝(2-DP)LgRi+CP将总体积VPi随孔隙半径增量的绝对值|dVPi/dRi|与孔隙半径绘制在双对数坐标，其中VPi为不同注入压力Pi下进入水泥环试样中孔隙的总水银体积；Ri为对应注入压力Pi下对应的孔隙半径；DP为水泥环试样中孔隙形态的分形维数；CP为水泥环试样孔隙的分形模型常数；(3)通过Lg(|dVPi/dRi|)＝(2-DP)LgRi+CP曲线斜率计算空白组和条件对照组水泥环孔隙形态的分形维数，令DBP为空白组水泥环孔隙形态的分形维数，DEP为条件对照组水泥环孔隙形态的分形维数；
步骤五：测定及评价有无交变载荷作用的水泥环本体颗粒形态的分形维数：(1)采用扫描电镜对步骤一得到的空白组及条件对照组水泥石SEM试样进行面扫描，获取空白组及条件对照组水泥石在不同放大倍率下的SEM图像；(2)采用Python+OpenCV程序对SEM图像进行二值化处理，得到不同阀值下的二值化图像，白色区域代表微观颗粒，黑色区域代表微观孔隙；(3)基于步骤四采用压汞法获取的水泥石试样真实孔隙率
以水泥石试样的真实孔隙率为控制对象，采用基于边缘强度的阀值分割算法，对二值化图像阀值进行适配调控，优选出与水泥石试样真实孔隙率相等的目标二值化图像；(4)采用Matlab编程计算目标二值化图像中所有白色区域的面积和周长；(5)采用面积-周长分形模型Lg(AGi)＝DG*Lg(PGi)+CG将目标二值化图像中白色区域的面积和周长绘制在双对数坐标，其中PGi为目标二值化图像中任一个白色几何图形的等效周长，AGi为目标二值化图像中等效周长为PGi的几何图形等效面积，DG为水泥石试样颗粒形态的分形维数，CG为水泥环试样颗粒的分形模型常数；(6)通过Lg(AGi)＝DG*Lg(PGi)+CG曲线斜率计算空白组和条件对照组水泥环颗粒形态的分形维数，令DBG为空白组水泥环颗粒形态的分形维数，DEG为条件对照组水泥环颗粒形态的分形维数；
步骤六：测定及评价交变载荷作用后水泥环本体裂纹形态的分形维数：(1)采用扫描电镜对步骤一得到的条件对照组水泥石SEM试样进行面扫描，获取条件对照组水泥石在不同放大倍率下的SEM图像；(2)采用与步骤五相同的方法获取水泥石裂纹的目标二值化图像；(3)采用Matlab编程计算用不同边长δFi的正方形盒子覆盖目标二值化图像的盒子总数N(δFi)；(4)采用盒子模型LgN(δFi)＝DEF*LgδFi+CF将正方形盒子边长δFi和盒子总数N(δFi)绘制在双对数坐标，其中δFi为测量目标二值化图像用正方形盒子的边长，N(δFi)为用边长为δFi的盒子去覆盖目标二值化图像所需要的盒子总数；DEF为水泥环本体裂纹形态的分形维数；CEF为水泥环本体裂纹形态的分形模型常数；(5)通过LgN(δFi)＝DEF*LgδFi+CEF曲线斜率计算条件对照组水泥环本体裂纹形态的分形维数DEF；
步骤七：基于步骤二和步骤三，建立空白组界面径向胶结强度与套管水泥环界面形态分形维数之间的函数关系FB1(SBR，DBTF)和FB2(SBR，DBCF)；
步骤八：基于步骤二和步骤四，建立空白组水泥石力学性能与套管水泥环本体孔隙形态分形维数之间的函数关系FB3(QBL，DBP)和FB4(QBC，DBP)；
步骤九：基于步骤二和步骤五，建立空白组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FB5(QBL，DBG)和FB6(QBC，DBG)；
步骤十：基于步骤二和步骤三，建立条件对照组界面径向胶结强度与套管水泥环界面形态分形维数之间的函数关系FE1(SER，DETF)和FE2(SER，DECF)；
步骤十一：基于步骤二和步骤四，建立条件对照组水泥石力学性能与套管水泥环本体孔隙形态分形维数之间的函数关系FE3(QEL，DEP)和FE4(QEC，DEP)；
步骤十二：基于步骤二和步骤五，建立条件对照组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FE5(QEL，DEG)和FE6(QEC，DEG)；
步骤十三：基于步骤二和步骤六，建立条件对照组水泥石力学性能与套管水泥环本体颗粒形态分形维数之间的函数关系FE7(QEL，DEF)和FE8(QEC，DEF)；
步骤十四：采用分形维数DBTF、DBCF、DBP、DBG和函数关系FB1(SBR，DBTF)、FB2(SBR，DBCF)、FB3(QBL，DBP)、FB4(QBC，DBP)、FB5(QBL，DBG)及FB6(QBC，DBG)定量评价空白组水泥环完整性；
步骤十五：采用分形维数DETF、DECF、DEP、DEG、DEF和函数关系FE1(SER，DETF)、FE2(SER，DECF)、FE3(QEL，DEP)、FE4(QEC，DEP)FE5(QEL，DEG)、FE6(QEC，DEG)、FE7(QEL，DEF)和FE8(QEC，DEF)定量评价条件对照组水泥环完整性损伤演化规律。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211542803.0A CN115711120B (zh) | 2022-12-02 | 2022-12-02 | 一种油气井水泥环完整性及损伤演化定量评价方法 |
| US18/313,210 US11965412B1 (en) | 2022-12-02 | 2023-05-05 | Quantitative evaluation method for integrity and damage evolution of cement sheath in oil-gas well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211542803.0A CN115711120B (zh) | 2022-12-02 | 2022-12-02 | 一种油气井水泥环完整性及损伤演化定量评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115711120A true CN115711120A (zh) | 2023-02-24 |
| CN115711120B CN115711120B (zh) | 2024-04-12 |
Family
ID=85235569
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211542803.0A Active CN115711120B (zh) | 2022-12-02 | 2022-12-02 | 一种油气井水泥环完整性及损伤演化定量评价方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US11965412B1 (zh) |
| CN (1) | CN115711120B (zh) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014148925A1 (en) * | 2013-03-22 | 2014-09-25 | Auckland Uniservices Limted | Method and system for monitoring and/or controlling fracture connectivity |
| CN104131803A (zh) * | 2013-05-10 | 2014-11-05 | 中国石油大学(北京) | 一种评价页岩压裂形成缝网能力的实验方法 |
| US20170205388A1 (en) * | 2014-04-04 | 2017-07-20 | Schlumberger Technology Corporation | Wellbore cement simulator |
| CN106991235A (zh) * | 2017-04-01 | 2017-07-28 | 中国石油天然气集团公司 | 固井水泥环完整性分析预测方法及装置 |
| CN107991165A (zh) * | 2017-11-03 | 2018-05-04 | 中国石油大学(北京) | 一种水泥环气密封性失效判断方法 |
| CN109238868A (zh) * | 2018-11-20 | 2019-01-18 | 山东农业大学 | 一种微破损现场检测混凝土结构抗压强度的测试方法 |
| CN110705112A (zh) * | 2019-10-09 | 2020-01-17 | 江西理工大学 | 基于块度分形理论的岩石损伤演化本构模型建立方法 |
| CN111080151A (zh) * | 2019-12-23 | 2020-04-28 | 北京理工大学 | 一种基于断口形貌分形特征的裂纹扩展评价方法 |
| CN111965032A (zh) * | 2020-06-04 | 2020-11-20 | 江西理工大学 | 不同配比钽铌矿废弃尾砂胶结材料的力学损伤分析方法 |
| US20210082543A1 (en) * | 2020-09-07 | 2021-03-18 | Southwest Petroleum University | Method for establishing chart for designing mechanical properties of cement stones in large-scale fracturing oil well |
| CN113075105A (zh) * | 2021-03-25 | 2021-07-06 | 河南理工大学 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
| US20210231520A1 (en) * | 2020-08-26 | 2021-07-29 | Southwest Petroleum University | Sealing integrity evaluation device for high-temperature and high-pressure casing-cement ring-formation and method thereof |
| CN114215505A (zh) * | 2021-12-09 | 2022-03-22 | 西南石油大学 | 一种载荷谱作用下水泥环环向及径向变形量化评价装置及方法 |
| CN114896775A (zh) * | 2022-04-28 | 2022-08-12 | 西南石油大学 | 一种评价固井水泥环在压裂工况下层间封隔的方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009140738A1 (en) * | 2008-05-23 | 2009-11-26 | The Australian National University | Image data processing |
| MX350511B (es) * | 2011-07-12 | 2017-09-04 | Ingrain Inc | Método para simular flujo de fase múltiple/componente múltiple fraccional a través de medios porosos. |
| US10150904B1 (en) * | 2018-03-30 | 2018-12-11 | King Fahd University Of Petroleum And Minerals | Nano zeolite cement additive and methods of use |
| CA3044153C (en) * | 2018-07-04 | 2020-09-15 | Eavor Technologies Inc. | Method for forming high efficiency geothermal wellbores |
-
2022
- 2022-12-02 CN CN202211542803.0A patent/CN115711120B/zh active Active
-
2023
- 2023-05-05 US US18/313,210 patent/US11965412B1/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014148925A1 (en) * | 2013-03-22 | 2014-09-25 | Auckland Uniservices Limted | Method and system for monitoring and/or controlling fracture connectivity |
| CN104131803A (zh) * | 2013-05-10 | 2014-11-05 | 中国石油大学(北京) | 一种评价页岩压裂形成缝网能力的实验方法 |
| US20170205388A1 (en) * | 2014-04-04 | 2017-07-20 | Schlumberger Technology Corporation | Wellbore cement simulator |
| CN106991235A (zh) * | 2017-04-01 | 2017-07-28 | 中国石油天然气集团公司 | 固井水泥环完整性分析预测方法及装置 |
| CN107991165A (zh) * | 2017-11-03 | 2018-05-04 | 中国石油大学(北京) | 一种水泥环气密封性失效判断方法 |
| CN109238868A (zh) * | 2018-11-20 | 2019-01-18 | 山东农业大学 | 一种微破损现场检测混凝土结构抗压强度的测试方法 |
| CN110705112A (zh) * | 2019-10-09 | 2020-01-17 | 江西理工大学 | 基于块度分形理论的岩石损伤演化本构模型建立方法 |
| CN111080151A (zh) * | 2019-12-23 | 2020-04-28 | 北京理工大学 | 一种基于断口形貌分形特征的裂纹扩展评价方法 |
| CN111965032A (zh) * | 2020-06-04 | 2020-11-20 | 江西理工大学 | 不同配比钽铌矿废弃尾砂胶结材料的力学损伤分析方法 |
| US20210231520A1 (en) * | 2020-08-26 | 2021-07-29 | Southwest Petroleum University | Sealing integrity evaluation device for high-temperature and high-pressure casing-cement ring-formation and method thereof |
| US20210082543A1 (en) * | 2020-09-07 | 2021-03-18 | Southwest Petroleum University | Method for establishing chart for designing mechanical properties of cement stones in large-scale fracturing oil well |
| CN113075105A (zh) * | 2021-03-25 | 2021-07-06 | 河南理工大学 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
| CN114215505A (zh) * | 2021-12-09 | 2022-03-22 | 西南石油大学 | 一种载荷谱作用下水泥环环向及径向变形量化评价装置及方法 |
| CN114896775A (zh) * | 2022-04-28 | 2022-08-12 | 西南石油大学 | 一种评价固井水泥环在压裂工况下层间封隔的方法 |
Non-Patent Citations (11)
| Title |
|---|
| DENG, KH: "Experimental study on the integrity of casing-cement sheath in shale gas wells under pressure and temperature cycle loading", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》, 31 December 2020 (2020-12-31) * |
| ZHOU, Y: "Study of the effect of loading modes on the acoustic emission fractal and damage characteristics of cemented paste backfill", 《CONSTRUCTION AND BUILDING MATERIALS》, 29 March 2021 (2021-03-29) * |
| 冯颖韬: "基于HYMOSTRUC3D的固井水泥石孔隙结构演变及力学强度发展规律", 《 科学技术与工程》, vol. 21, no. 12, 28 April 2021 (2021-04-28) * |
| 刘志伟: "生产阶段油气井井筒完整性风险评价研究", 《中国优秀硕士论文全文库工程科技I辑》, no. 5, 15 May 2020 (2020-05-15) * |
| 吴磊: "提高固井水泥石韧性方法研究", 《中国优秀硕士论文全文库工程科技I辑》, no. 4, 15 April 2021 (2021-04-15) * |
| 周照明: "气井水泥环微间隙封堵材料粒径分形级配优选方法", 《科学技术与工程》, vol. 10, no. 8, 30 April 2010 (2010-04-30) * |
| 林元华: "强交变热载荷下页岩气井水泥环完整性测试", 《天然气工业.》, vol. 40, no. 5, 31 May 2020 (2020-05-31) * |
| 王杰: "基于分形理论超细粉煤灰微珠低密度水泥浆的制备", 《硅酸盐通报》, vol. 30, no. 2, 15 April 2011 (2011-04-15) * |
| 聂臻: "水泥环密封失效机理研究及应用", 《中国优秀硕士论文全文库工程科技I辑》, no. 5, 15 May 2015 (2015-05-15) * |
| 解静: "微结构分维与孔隙率及抗压强度的关系", 《兰州大学学报(自然科学版)》, vol. 45, no. 2, 30 April 2009 (2009-04-30) * |
| 黄进: "盐层井段固井水泥石力学性能及微观结构研究", 《长江大学学报(自科版)》, vol. 11, no. 26, 15 September 2014 (2014-09-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115711120B (zh) | 2024-04-12 |
| US11965412B1 (en) | 2024-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109632510B (zh) | 一种预测水化损伤页岩强度的方法 | |
| CN110644980B (zh) | 一种超低渗透油藏储层综合分类评价方法 | |
| CN104405366B (zh) | 一种高温高压固井水泥环力学完整性测试装置及方法 | |
| CN107451325B (zh) | 深井超深井压裂套管失效风险实时定量评估方法及装置 | |
| CN108590601B (zh) | 一种注水扩容膨胀施工参数优选的实验方法 | |
| US10928282B2 (en) | Method and device for determining elasticity of cement stone utilized in well cementing of oil-gas well | |
| CN109142192B (zh) | 可视化异形固井二界面胶结质量测试系统 | |
| CN108894777B (zh) | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 | |
| CN108732076B (zh) | 一种煤层水力压裂渗透率预测方法 | |
| CN105930932B (zh) | 基于含气指数的页岩气层标准化无阻流量的获取方法 | |
| CN108414351B (zh) | 一种套管非均匀外挤能力评价方法 | |
| CN111444461A (zh) | 高水压下围岩大变形灾害等级预测方法 | |
| CN114048686A (zh) | 一种刚性堵漏材料定量评价方法 | |
| CN109241651A (zh) | 一种通用的地层破裂压力预测方法 | |
| CN115711120A (zh) | 一种油气井水泥环完整性及损伤演化定量评价方法 | |
| CN111929424B (zh) | 一种基于尺寸效应的巨跨地下洞室硬质围岩亚分级方法 | |
| CN110095236B (zh) | 一种封隔器胶筒密封性能测试装置 | |
| CN116106139A (zh) | 一种温差条件下天然裂缝剪切强度的测量方法 | |
| CN111982779B (zh) | 一种空心圆柱岩石试样模拟有压隧道渗流变形的试验方法 | |
| CN111305794A (zh) | 一种保持煤岩高渗透率的方法 | |
| Wu et al. | Analysis of the microstructure of a failed cement sheath subjected to complex temperature and pressure conditions | |
| CN112052595A (zh) | 一种腐蚀钢制管道外压临界弹塑性屈曲压力计算方法 | |
| CN112881184A (zh) | 基于全程应力应变的砂砾岩脆性指数确定方法及系统 | |
| CN115949367B (zh) | 一种深层裂缝漏失性地层堵漏封堵层承压能力预测方法 | |
| CN115874978B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |