CN114970235B - 裂缝性储层地质力学非均质-各向异性建模方法 - Google Patents
裂缝性储层地质力学非均质-各向异性建模方法 Download PDFInfo
- Publication number
- CN114970235B CN114970235B CN202210394190.4A CN202210394190A CN114970235B CN 114970235 B CN114970235 B CN 114970235B CN 202210394190 A CN202210394190 A CN 202210394190A CN 114970235 B CN114970235 B CN 114970235B
- Authority
- CN
- China
- Prior art keywords
- rock
- stress
- fracture
- fractured
- model
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 239000011435 rock Substances 0.000 claims abstract description 113
- 238000002474 experimental method Methods 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 238000013178 mathematical model Methods 0.000 claims abstract description 15
- 238000011160 research Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 13
- 206010017076 Fracture Diseases 0.000 description 79
- 208000010392 Bone Fractures Diseases 0.000 description 74
- 238000012360 testing method Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241000269319 Squalius cephalus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Classifications
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Computer Graphics (AREA)
- Geochemistry & Mineralogy (AREA)
- Software Systems (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及油气田勘探开发领域,尤其是裂缝性储层地质力学非均质‑各向异性建模方法。通过裂缝面力学实验,建立裂缝面法向、切向刚度系数与正应力间的数学模型,依据裂缝面三维激光扫描实验结果,建立含裂缝面有限元几何模型,依据裂缝密度与岩体等效杨氏模量、泊松比之间的数学模型,建立裂缝性岩体地质力学非均质模型,并在裂缝性岩体地质力学非均质模型不同方向上赋予不同的岩石力学参数,建立裂缝性储层地质力学非均质‑各向异性模型,得到裂缝性岩体应力大小和方向。本发明专利基于物理模拟与数值模拟的方法提出了一种裂缝性储层地质力学非均质‑各向异性建模方法,实现了不连续岩体地质力学精细建模。
Description
技术领域
本发明涉及油气田勘探开发领域,尤其是裂缝性储层地质力学非均质-各向异性建模方法。
背景技术
在我国的含油气盆地中,裂缝性油气藏是一种重要的油气藏类型。这类油藏的探明地质储量已经超过40亿吨,占目前探明油气资源总量的1/3以上,另据有关专家预测,在我国130多亿吨的剩余油气资源量中,大约有60%左右的油气资源量分布与储层裂缝有关。因此,裂缝性油气藏的勘探开发在我国石油工业中的地位越来越重要,如何有效地勘探、合理地动用与高效地开发这些裂缝性油气藏,对我国石油工业的持续稳定发展以及增强我国能源安全供应的保障能力,确保国家安全与经济协调、持续快速发展具有长远的战略意义。
储层裂缝研究已有近百年的历史,国外以Nelson、Warren、Narr、Pollard、VanGolf-Racht、Lorenz、Aguilera以及Смехов等,国内以曾联波、宋惠珍、袁士义、刘建中、穆龙新、戴俊生、周新桂以及侯贵廷等为代表的学者们开展了大量的工作,为储层裂缝的研究奠定了良好的理论基础。由于储层裂缝的复杂性以及人们认识水平和技术水平的限制,在很长一段时间内裂缝性岩体地质力学建模以均质建模为主,即不考虑裂缝对岩体力学参数的影响。实际上,岩石的力学性质会影响天然裂缝发育程度与成因机制,同样,裂缝发育也会影响岩石力学参数的大小与各向异性。因此裂缝性储层地质力学建模与应力场数值模拟,应该考虑其非均质性与各向异性。其中,非均质性是指模型中不同单元体力学性质大小的差异性,而各向异性是指模型中同一个单元体的不同方向岩体力学性质的差异性。
本发明专利通过裂缝面力学实验,建立裂缝面法向、切向刚度系数与正应力的数学模型,开展裂缝面三维激光扫描实验,建立含裂缝面有限元几何模型,依据不同裂缝的几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,赋予裂缝性岩体几何模型等效力学参数,建立裂缝性岩体地质力学非均质模型;在裂缝性岩体地质力学非均质模型不同方向上赋予不同的岩石力学参数,建立裂缝性储层地质力学非均质-各向异性模型,最后依据研究区应力大小、方向,确定边界应力大小与施力方向,模拟确定裂缝性岩体应力大小和方向。
发明内容
本发明旨在解决上述问题,提供了裂缝性储层地质力学非均质-各向异性建模方法,它实现了裂缝性岩体地质力学非均质与各向异性精细建模。
本发明的技术方案为:裂缝性储层地质力学非均质-各向异性建模方法,具体步骤如下:
第一步依据裂缝的形成时期,针对不同组系的裂缝,选取含裂缝面岩石样品,开展不同正应力、剪应力条件下的裂缝面力学实验,得到含裂缝面岩石样品的应力-应变关系曲线;
第二步利用不同正应力、剪应力条件下岩石应力-应变关系曲线,计算裂缝面的法向、切向刚度系数,分别建立裂缝面法向、切向刚度系数与正应力之间的数学模型;
第三步针对不同组系的裂缝面,开展裂缝面三维激光扫描实验,得到裂缝面三维几何数据体;
第四步选取不含裂缝的样品,开展岩石三轴力学实验,确定完整岩石的杨氏模量、泊松比、密度,所述的完整岩石是指裂缝不发育的岩石;
第五步利用裂缝面三维几何数据体,结合裂缝的几何样式,建立含裂缝面有限元几何模型,利用岩石三轴力学实验,开展数值模拟,建立不同裂缝的几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,同时,揭示岩石力学参数与空间方位之间的数学关系,明确不同方向岩石力学参数大小;
第六步结合裂缝的几何样式与裂缝密度,建立裂缝性岩体几何模型,依据岩石三轴力学实验结果,对裂缝性岩体几何模型网格划分,赋予不同单元体相应的力学参数,建立裂缝性岩体均质模型;
第七步依据不同裂缝几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,结合裂缝性岩体几何模型中裂缝密度分布,赋予裂缝性岩体均质模型中不同单元体对应的等效力学参数,建立裂缝性岩体地质力学非均质模型;
第八步依据岩石力学参数与空间方位之间的数学关系,在裂缝性岩体地质力学非均质模型不同方向上赋予不同的岩石力学参数,建立裂缝性储层地质力学非均质-各向异性模型;
第九步利用研究区应力大小、方向,确定边界应力大小与施力方向,模拟确定裂缝性岩体应力大小和方向。
本发明的有益效果是:通过裂缝面力学实验,建立裂缝面法向、切向刚度系数与正应力的数学模型,开展裂缝面三维激光扫描实验,建立含裂缝面有限元几何模型,依据不同裂缝的几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,赋予裂缝性岩体几何模型等效力学参数,建立裂缝性岩体地质力学非均质模型;在裂缝性岩体地质力学非均质模型不同方向上赋予不同的岩石力学参数,建立裂缝性储层地质力学非均质-各向异性模型,最后依据研究区应力大小、方向,确定边界应力大小与施力方向,模拟确定裂缝性岩体应力大小和方向。本发明专利从物理模拟与数值模拟的角度提出了一种裂缝性储层地质力学非均质-各向异性建模方法,具有较高的实用价值,可操作性强,能大幅提高裂缝性储层地质力学建模精度,揭示不连续岩体地应力分布规律,为我国目前裂缝性油气藏地质工程一体化实践提供技术支持。
附图说明
图1为裂缝性储层地质力学非均质-各向异性建模方法的流程图。
图2为裂缝面岩石力学实验样品照片。
图3为裂缝面法向应力–法向位移曲线。
图4为不同样品裂缝面切向应力–位移曲线。
图5为裂缝面三维激光扫描结果。
图6为完整岩石在60MPa围压下的应力、应变关系曲线。
图7为裂缝线密度与岩体等效力学参数之间的数学关系。
图8为不同方向上裂缝性岩体等效力学参数差异性。
图9为裂缝性岩体几何模型。
图10为裂缝发育带地质力学模型。
图11为裂缝发育带地质力学模型边界条件。
图12为裂缝性岩均质模型、非均质模型、非均质-各向异性模型。
具体实施方式
下面结合附图说明本发明的具体实施方式:
本发明专利以裂缝性致密砂岩为例为例,说明本发明具体实施过程。
第1步依据裂缝的形成时期,针对不同组系的裂缝,选取含裂缝面岩石样品(图2),开展不同正应力、剪应力条件下的裂缝面力学实验,得到含裂缝面岩石样品的应力-应变关系曲线。第2步针对不同组系的裂缝,开展不同正应力、剪应力条件下岩石应力-应变测试(图3),计算裂缝面的法向、切向刚度系数,分别建立裂缝面法向、切向刚度系数与正应力之间的数学模型(图4)。
在岩石力学参数中,主要考虑岩石的杨氏模量、泊松比对应力的影响,而裂缝的发育同样能够改变岩石的力学性质。通过力学实验模拟不同法向应力作用下的裂缝面剪应力-剪切变形实验曲线,实验为高压岩石三轴动态测试系统。通过裂缝性岩体试样进行室内闭合试验,得到裂缝法向应力-法向位移关系曲线,进行模型进行数据处理得到,法向刚度与法向应力的关系,不同样品参数测试结果如表1所示。
Kn=Aσn B (1)
公式(1)中,Kn为裂缝面法向刚度系数,σn为裂缝面法向应力,A、B为法向刚度与法向应力的拟合系数。
通过测量不同法向应力作用下的裂缝面剪应力-剪切变形试验曲线。定义剪切位移0.1mm时对应的斜率为裂缝面的剪切刚度,得到裂缝面剪切刚度与法向应力的关系,不同样品参数测试结果如表2所示。
Ks=Cσn D (2)
公式(2)中,Ks为裂缝面剪切刚度系数,σn为裂缝面法向应力,C、D为法向刚度与法向应力的拟合系数。
表1不同样品参数测试结果
A | B | R2 | |
样品1 | 167.3 | 0.32 | 0.94 |
样品2 | 183.1 | 0.28 | 0.99 |
样品3 | 163.6 | 0.33 | 0.94 |
样品4 | 172.9 | 0.29 | 0.99 |
表2不同样品参数测试结果
实验结果表明,不同样品之间裂缝面剪切刚度数值受法向应力影响,两者呈幂函数变化。
第3步针对不同组系的裂缝面,开展裂缝面三维激光扫描实验,得到裂缝面三维几何数据体。
岩石断裂表面的粗糙形态表现出良好的自相似性,即从不同的尺度观测它,均可以观测到相似的粗糙不平的特征,因此,研究很小的断裂面便可以反映宏观裂缝的力学行为。利用三维激光扫描设备获取裂缝面的三维数据,拾取裂缝面关键点的坐标。三维激光扫描实验在三维激光扫描仪测试完成,最后,采用裂缝面三维数据处理软件得到上、下裂缝面数据的三维坐标(图5)。
第4步选取不含裂缝的样品,开展岩石三轴力学实验,依据致密砂岩的埋深确定岩石力学实验围压,得到岩石应力-应变曲线(图6),确定完整岩石的杨氏模量、泊松比、密度。
第5步利用裂缝面三维几何数据体,研究区裂缝为平行裂缝,建立含裂缝面有限元几何模型,在有限元软件中建立上、下两个裂缝面的三维分布数据点,建立裂缝面几何模型。利用岩石三轴力学实验,开展数值模拟,建立不同裂缝的几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型(图7),对四个样品裂缝面密度、杨氏模量变化率、泊松比变化率进行数值统计并绘制点状图,可以看出,随着裂缝面密度的增加,四个样品杨氏模量变化率均减小,泊松比变化率均增大(图7)。由此可知,随着裂缝面密度的增加,裂缝性岩体等效杨氏模量呈线性减小,泊松比呈线性增大。同时,揭示岩石力学参数与空间方位之间的数学关系,明确不同方向岩石力学参数大小(图8)。
第6步结合裂缝的几何样式与裂缝密度,建立裂缝性岩体几何模型(图9),依据岩石三轴力学实验结果,采用SOLID45单元类型最符合砂岩储层岩石的力学特性。网格边长的选择要视具体情况而定,划分之后形成的网格不宜过多或过少,过多会使单次数值模拟计算的时间变长,而且对模拟精度的提高并无太大意义,太少则会使数值模拟结果过于粗糙,影响数值模拟效果。相对于非裂缝区域,裂缝带的宽度一般较小,因此网格需要细分。根据上述原则,对地质模型进行了网格划分,共划分出27625个节点、89523个单元(图10);最终通过裂缝性岩体几何模型网格划分,赋予不同单元体相应的力学参数,建立裂缝性岩体均质模型。
第7步依据不同裂缝几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,结合裂缝性岩体几何模型中裂缝密度分布,赋予裂缝性岩体均质模型中不同单元体对应的等效力学参数,建立裂缝性岩体地质力学非均质模型。
第8步依据岩石力学参数与空间方位之间的数学关系,在裂缝性岩体地质力学非均质模型不同方向上赋予不同的岩石力学参数,建立裂缝性储层地质力学非均质-各向异性模型。
第9步利用研究区应力大小、方向,确定边界应力大小与施力方向,模拟确定裂缝性岩体应力大小和方向。
利用测井资料解释得到的现今地应力值是单个井点的数值,受构造起伏等因素的影响,单井地应力值不能直接作为区域应力载荷施加在力学模型上,而是需要通过数值模拟进行拟合,从而得到合适的区域地应力值。按照以上思路,对裂缝发育带地质力学模型的现今地应力场进行了数值模拟。区域地应力值主要依靠已确定的关键井地应力,结合研究区的实际地质情况,利用线性叠加原理,经多次反演拟合后最终确定。由于应力分布的复杂性,使数值模拟时的边界条件难以一次性地准确设定,因此地应力场的数值模拟实际上是一个正演与反演多次反复计算、修正的过程。反演的标准是关键井地应力的数值模拟结果与实际测算结果误差最小,且要兼顾应力方位和应力大小两方面,最终确定应力,得到边界条件,如图11所示,在模型北西南东方向施加130MPa的压应力,在北东-南西方向施加170MPa的压应力。
如图12所示,均值、非均质、非均质-各向异性裂缝性岩体模型采取不同的建模方式,模拟结果显示出明显的差异性。均质裂缝性岩体建模应力结果显示出基本均匀的应力分布的特征,局部降低,应力大小范围为152-148MPa,应力方向呈北东向。非均质裂缝性岩体建模应力结果显示出裂缝带内部应力值高,仅靠裂缝带应力低,外围又升高的特征,应力大小范围为178-149MPa,应力方向呈北东向。非均质-各向异性裂缝性岩体模型应力结果显示裂缝发育带应力较低,向外围呈两个环带逐渐升高的特征,应力大小为154-141MPa,靠近裂缝发育带,应力方向呈弧线式相切与裂缝发育带,远离裂缝发育带,应力方向呈北东向。
上面以举例方式对本发明进行了说明,但本发明不限于上述具体实施例,凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。
Claims (1)
1.裂缝性储层地质力学非均质-各向异性建模方法,实现的步骤如下:
第一步 依据裂缝的形成时期,针对不同组系的裂缝,选取含裂缝面岩石样品,开展不同正应力、剪应力条件下的裂缝面力学实验,得到含裂缝面岩石样品的应力-应变关系曲线;
第二步 利用不同正应力、剪应力条件下岩石应力-应变关系曲线,计算裂缝面的法向、切向刚度系数,分别建立裂缝面法向、切向刚度系数与正应力之间的数学模型;
第三步 选取不同组系的裂缝面,开展裂缝面三维激光扫描实验,得到裂缝面三维几何数据体;
第四步 选取不含裂缝的样品,开展岩石三轴力学实验,确定完整岩石的杨氏模量、泊松比以及密度,所述的完整岩石是指裂缝不发育的岩石;
第五步 利用裂缝面三维几何数据体,结合裂缝的几何样式,建立含裂缝面有限元几何模型,利用岩石三轴力学实验,开展数值模拟,建立不同裂缝的几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,揭示岩石力学参数与空间方位之间的数学关系,明确不同方向岩石力学参数大小;
第六步 结合裂缝的几何样式与裂缝密度,建立裂缝性岩体几何模型,依据岩石三轴力学实验结果,对裂缝性岩体几何模型网格划分,赋予不同单元体相应的力学参数,建立裂缝性岩体均质模型;
第七步 依据不同裂缝几何样式下,裂缝密度与裂缝性岩体等效杨氏模量、泊松比之间的数学模型,结合裂缝性岩体几何模型中裂缝密度分布,赋予裂缝性岩体均质模型中不同单元体对应的等效力学参数,建立裂缝性岩体地质力学非均质模型;
第八步 依据岩石力学参数与空间方位之间的数学关系,在裂缝性岩体地质力学非均质模型不同方向上赋予不同的岩石力学参数,建立裂缝性储层地质力学非均质-各向异性模型;
第九步 利用研究区应力大小、方向,确定边界应力大小与施力方向,模拟确定裂缝性岩体应力大小和方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210394190.4A CN114970235B (zh) | 2022-04-15 | 2022-04-15 | 裂缝性储层地质力学非均质-各向异性建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210394190.4A CN114970235B (zh) | 2022-04-15 | 2022-04-15 | 裂缝性储层地质力学非均质-各向异性建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114970235A CN114970235A (zh) | 2022-08-30 |
CN114970235B true CN114970235B (zh) | 2024-04-02 |
Family
ID=82976707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210394190.4A Active CN114970235B (zh) | 2022-04-15 | 2022-04-15 | 裂缝性储层地质力学非均质-各向异性建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114970235B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115795828B (zh) * | 2022-11-16 | 2024-03-19 | 武汉大学 | 基于数据驱动的裂隙岩体变形计算方法及装置 |
CN117725807B (zh) * | 2024-02-18 | 2024-05-14 | 中国科学院地质与地球物理研究所 | 基于边界元计算非均质储层裂缝应力的方法、系统及设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110779795A (zh) * | 2019-11-04 | 2020-02-11 | 中国石油大学(华东) | 裂缝性储层地质力学建模网格单元大小确定方法 |
CN113919196A (zh) * | 2021-09-26 | 2022-01-11 | 中国石油大学(华东) | 储层三维应力场模拟方法、模拟系统、终端、存储介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2690991C (en) * | 2007-08-24 | 2013-12-24 | Exxonmobil Upstream Research Company | Method for multi-scale geomechanical model analysis by computer simulation |
US20170145793A1 (en) * | 2015-08-20 | 2017-05-25 | FracGeo, LLC | Method For Modeling Stimulated Reservoir Properties Resulting From Hydraulic Fracturing In Naturally Fractured Reservoirs |
CN110704888B (zh) * | 2019-07-04 | 2022-07-29 | 成都理工大学 | 一种非常规油气藏加密井体积压裂施工参数优化设计方法 |
-
2022
- 2022-04-15 CN CN202210394190.4A patent/CN114970235B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110779795A (zh) * | 2019-11-04 | 2020-02-11 | 中国石油大学(华东) | 裂缝性储层地质力学建模网格单元大小确定方法 |
CN113919196A (zh) * | 2021-09-26 | 2022-01-11 | 中国石油大学(华东) | 储层三维应力场模拟方法、模拟系统、终端、存储介质 |
Non-Patent Citations (2)
Title |
---|
页岩储层地应力研究进展;林长城;丁文龙;王兴华;谷阳;尹帅;刘敬寿;;科技通报;20180731(07);全文 * |
页岩储层构造应力场模拟与裂缝分布预测方法及应用;丁文龙;曾维特;王濡岳;久凯;王哲;孙雅雄;王兴华;;地学前缘;20160315(02);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114970235A (zh) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110779795B (zh) | 裂缝性储层地质力学建模网格单元大小确定方法 | |
CN114970235B (zh) | 裂缝性储层地质力学非均质-各向异性建模方法 | |
CN113901681B (zh) | 一种全寿命周期页岩气储层双甜点三维可压性评估方法 | |
CN110147561B (zh) | 一种含天然裂缝致密油气储层体积压裂缝网预测方法 | |
CN109583113B (zh) | 一种岩石地层压实系数和有效孔隙体积压缩系数计算方法 | |
CN105424466A (zh) | 一种评价波浪作用下砂土对埋置海底管线抗力的方法 | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
CN111006987B (zh) | 一种储层裂缝孔渗参数预测方法 | |
CN104865610A (zh) | 一种深部储层岩心的数值表征方法 | |
Fang et al. | Analysis of Changes in the Micromorphology of Sandstone Joint Surface under Dry‐Wet Cycling | |
Zhao et al. | A new method to investigate the size effect and anisotropy of mechanical properties of columnar jointed rock mass | |
CN108875115B (zh) | 一种确定岩石强度的方法 | |
Zhang et al. | Fracability evaluation method of a fractured-vuggy carbonate reservoir in the shunbei block | |
CN111680380A (zh) | 基于地质力学特征空间展布的全三维压裂设计方法 | |
CN113280951B (zh) | 一种建立峡谷区斜坡地应力场分布的方法 | |
CN106638508A (zh) | 高地应力条件下大型洞室轴线选择方法 | |
CN108442918A (zh) | 一种压裂裂缝双翼缝长定量预测方法 | |
Bessinger et al. | The role of compressive stresses in jointing on Vancouver Island, British Columbia | |
Ai et al. | Study on the Method for Measuring Stress on Rock Mass Excavation Surface under Extremely High Stress Conditions | |
Zhao et al. | Cyclic impact damage and water saturation effects on mechanical properties and Kaiser effect of red sandstone under uniaxial cyclic loading and unloading compression | |
Zhang et al. | Cenozoic fault distribution characteristics and evolution in Qikou sag of Bohai Basin, China | |
CN116050293B (zh) | 一种确定稠油泡沫调驱技术参数最优数值指示图版的方法 | |
CN112418597B (zh) | 一种适用于页岩气老区开发调整井的储层可压性评价方法 | |
Zhengzheng et al. | Experimental study on the fracture surface morphological characteristics and permeability characteristics of sandstones with different particle sizes | |
CN106321081A (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 |