CN109459498B - 一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 - Google Patents
一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 Download PDFInfo
- Publication number
- CN109459498B CN109459498B CN201811075122.1A CN201811075122A CN109459498B CN 109459498 B CN109459498 B CN 109459498B CN 201811075122 A CN201811075122 A CN 201811075122A CN 109459498 B CN109459498 B CN 109459498B
- Authority
- CN
- China
- Prior art keywords
- follows
- calculation formula
- crack
- rock
- matrix
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 82
- 239000011435 rock Substances 0.000 title claims abstract description 82
- 238000007598 dipping method Methods 0.000 title abstract 3
- 239000011159 matrix material Substances 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 20
- 238000002474 experimental method Methods 0.000 claims abstract description 6
- 239000012530 fluid Substances 0.000 claims description 10
- 230000010287 polarization Effects 0.000 claims description 8
- 229920006395 saturated elastomer Polymers 0.000 claims description 7
- 238000002591 computed tomography Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 abstract description 4
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000005755 formation reaction Methods 0.000 description 8
- 230000007547 defect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种含倾斜裂缝横向各向同性岩石的声波速度计算方法,包括以下步骤:S1:测量不含裂缝横向各向同性岩石的实验参数;S2:计算不含裂缝横向各向同性岩石的弹性刚度矩阵系数;S3:得到裂缝参数;S4:计算得到法向裂缝的法向柔度ZN、切向柔度ZT、柔度矩阵Z;S5:根据步骤S4得到的裂缝的法向柔度ZN、切向柔度ZT和柔度矩阵Z,计算得到含裂缝的横向各向同性岩石的弹性刚度矩阵Ceff;S6:根据步骤S5得到的弹性刚度矩阵Ceff的系数,计算声波速度。本发明的方法能够更加有效的预测含倾斜裂缝岩石的弹性性质,可以为地震勘探及测井识别、评价裂缝储层提供支持。
Description
技术领域
本发明涉及勘探地球物理领域,尤其涉及岩石物理中的一种含倾斜裂缝横向各向同性岩石的等效弹性性质及声波速度估算方法。
背景技术
地震勘探是评价裂缝性地层的常用方法。实际应用中,需要应用合适的岩石物理模型来建立地震测量得到的各向异性声速与地层裂缝特征参数的联系以评价地层中的裂缝特征。前人建立了很多岩石物理模型来研究含裂缝岩石的弹性性质。Hudson基于散射理论,针对含硬币状裂缝的各向同性岩石推导出了等效弹性模量的表达式。Cheng在Eshelby有关椭球形裂缝应变的研究基础上,提出了相应的模型。该模型可以用于计算含平行椭球形裂缝的各向同性岩石的等效弹性模量。上述模型都假设含裂缝岩石的背景介质为各向同性,实际含裂缝地层中,背景介质多由薄层构成而呈现明显的横向各向同性。Wang针对多种岩石声速各项异性的测量结果,当岩石不含裂缝时,各向异性比例大于10%。所以在岩石物理建模时,考察各向异性背景的影响能够更好的评价含裂缝地层以及更好的对比理论预测结果与实验测量结果。至今为止,少数岩石物理模型在建模时考察了TI背景介质的影响,Schoenberg和Helbig应用线性滑动理论推导出了含垂直裂缝的层状岩石的等效弹性模量。Sevostianov等基于Hill’s张量的推导方法得到了TI背景介质中裂缝的柔度矩阵。Xu等研究了TI背景介质中裂缝引起的散射现象并给出了含裂缝岩石的等效弹性参数。
然而,上述考虑背景各向异性的模型仅仅考察了平行或垂直于各向同性面的裂缝。实际裂缝性地层中,由于地应力的影响,裂缝通常不平行于各向同性平面,而是以一定角度与地层各向同性面斜交。而且Sevostianov和Xu等的模型并没有显性的给出裂缝柔度矩阵的表达式或表达形式极为复杂,这大大提高的了应用理论模型地震评价裂缝性地层的难度。
中国专利201410437668.2公开了一种各向同性地层弹性系数的测井计算方法,该方法包括:对岩心进行不同方向的取心,获得岩心样品;对岩心样品进行测量,获取岩心样品的纵波速度和横波速度,利用纵波速度和横波速度计算得到弹性系数C11,C13,C33,C44,C66;利用取心井的阵列声波测井数据和密度测井数据获取弹性系数C33,C44,C66;利用纵波速度和横波速度计算得到的弹性系数C11,C13,C33,C44,C66建立测井可测量的弹性系数C33,C44,C66和不可测量的弹性系数C11,C13之间的组合关系表达式;将通过取心井的阵列声波测井数据和密度测井数据获取的弹性系数C33,C44,C66代入组合关系表达式,得到测井不可测量的弹性系数C1。但是,该专利并没有给出含倾斜裂缝的横向各向同性岩石的弹性系数的计算方法。
中国专利201510674490.8公开了一种含垂直裂缝的页岩的波场传播特征的正演模拟方法,正演模拟方法包括:建立本构坐标系下的正交各项异性介质弹性矩阵;基于所构建的本构坐标系下的正交各项异性介质弹性矩阵,构建观测系统下的正交各向异性介质弹性矩阵;基于所构建的观测系统下的正交各向异性介质弹性矩阵,获得正交各向异性介质二阶波动方程;基于所述正交各向异性介质二阶波动方程,获得所述波场传播特征。但是该专利给出为垂直裂缝模型,并没有给出含倾斜裂缝的模型。
上述公开的现有技术中尽管公开了横向各向同性岩石中裂缝对岩石弹性性质的影响,但裂缝方向为垂直方向,并不适用于裂缝随机分布的岩石。因此,本次针对现有含裂缝岩石岩石物理模型的不足,推导出了一种考虑倾斜裂缝和各向异性背景介质影响的岩石物理模型,应用该模型可以快捷简便的计算含倾斜裂缝横向各向声波速度,并且可以获得等效弹性参数,应用更加方便。
发明内容
本发明的目的是针对现有含裂缝岩石物理模型的不足,提供一种含倾斜裂缝横向各向同性岩石的声波速度估算方法。该方法中,将背景介质设置为横向各向同性介质,裂缝形状设置为硬币状裂缝(天然裂缝大多纵横比较小),裂缝以一定角度与背景介质各向同性面斜交,更加符合真实含裂缝地层的情况,而且弹性参数计算公式形式更加简单,应用更加方便。
为了是实现上述目的,本发明提供了一种含倾斜裂缝横向各向同性岩石的声波速度的计算方法,包括以下步骤:
S1:测量不含裂缝横向各向同性岩石的实验参数,实验参数包括密度和各向异性声速,所述的声速包括入射方向于各向同性面平行的纵波速度Vp(90°)、入射方向于各向同性面呈45度夹角的纵波速度Vp(45°)、入射方向垂直于各向同性面的纵波速度Vp(0°)、偏振方向于各向同性面平行的横波速度VSH和偏振方向于各向同性面垂直的横波速度VSV。
S2:根据步骤S1测量的实验参数,计算不含裂缝横向各向同性岩石的弹性刚度矩阵系数;弹性刚度矩阵系数包括C11、C13、C12、C33、C44、C66和C1;
所述C11的计算公式为:
所述C12的计算公式为:
所述C33的计算公式为:
所述C44的计算公式为:
所述C13的计算公式为:
所述C66的计算公式为:
所述C1的计算公式为:
S3:对含裂缝的横向各向同性岩石做CT扫描实验,得到裂缝参数;所述的裂缝参数为裂缝密度e、裂缝纵横比α和裂缝与各向同性面的夹角θ;所述的裂缝为倾斜各向同性面硬币状裂缝;
S4:根据步骤S2中得到的刚度矩阵弹性常数和步骤S3中得到的裂缝参数计算得到裂缝的法向柔度ZN、切向柔度ZT、柔度矩阵Z;
法向柔度ZN的计算公式为:
切向柔度ZT的计算公式为:
柔度矩阵Z的计算公式为:
B3的计算公式为:
B4的计算公式为:
S5:根据步骤S4得到的法向柔度ZN、切向柔度ZT和柔度矩阵Z,计算得到含裂缝的横向各向同性岩石的弹性刚度矩阵Ceff;Ceff的计算方法为:先计算含裂缝的横向各向同性岩石的弹性柔度矩阵Meff,对矩阵Meff取逆矩阵,即得到所述的弹性刚度矩阵Ceff;
其中,Meff=M+Z,
所述M的计算公式为:
所述Ceff的计算公式为:Ceff=Meff-1。
S6:根据步骤S5得到的弹性刚度矩阵Ceff的系数,计算声波速度;
声波速度的计算公式为:
其中,C44 1、C22 1、C33 1、C55 1和C66 1是步骤S5中计算得到的含硬币裂缝岩石的弹性刚度矩阵Ceff中的弹性刚度矩阵系数;
ρ是岩石密度;
ψ是地震波相对于X3轴的入射角;
M计算公式为:
对于干燥硬币状裂缝的岩石,可以应用Brown–Korringa公式获得饱和流体状态的裂缝柔度矩阵:
其中,
Kf为流体体积模量;
为裂缝孔隙度,计算公式为:
Ks为背景介质的广义体积模量,计算公式为:
e为裂缝密度,计算公式为:
α为裂缝纵横比,a为X1-X2平面上的裂缝半径,n为裂缝数量,V为参考量。
本发明的有益效果为:
目前,针对含裂缝岩石的岩石物理模型仅仅考察了平行或垂直于各向同性面的裂缝。实际裂缝性地层中,由于地应力的影响,裂缝通常不平行于各向同性平面,而是以一定角度与各向同性面斜交。另外,现有模型并没有显性给出岩石弹性参数的计算方式,应用起来较为困难。本发明针对现有针对含裂缝介质的岩石物理模型的弊端,立足于含裂缝岩石的真实特征,基于横向各向同性的背景介质以及与各向同性面斜交且纵横比较小的硬币状裂缝推导出了含倾斜硬币状裂缝的横向各向同性岩石的弹性性质及声速计算模型。对比结果表明,本发明计算得到的等效弹性参数与其精确解吻合较好。本发明的方法能够更加有效的预测含倾斜裂缝岩石的弹性性质,可以为地震勘探及测井识别、评价裂缝储层提供支持。
附图说明
图1为含倾斜裂缝横向各向同性岩石的弹性刚度系数随裂缝倾斜角度的变化图。
图2为含倾斜裂缝横向各向同性的饱和流体岩石的弹性刚度系数随裂缝倾斜角度的变化图。
具体实施方式
一种含倾斜裂缝横向各向同性岩石的声波速度计算方法,包括如下步骤:
S1:测量不含裂缝横向各向同性岩石的实验参数,实验参数包括密度和各向异性声速,所述的声速包括入射方向于各向同性面平行的纵波速度Vp(90°)、入射方向于各向同性面呈45度夹角的纵波速度Vp(45°)、入射方向垂直于各向同性面的纵波速度Vp(0°)、偏振方向于各向同性面平行的横波速度VSH和偏振方向于各向同性面垂直的横波速度VSV。
S2:根据步骤S1测量的实验参数,计算不含裂缝横向各向同性岩石的弹性刚度矩阵系数;弹性刚度矩阵系数包括C11、C13、C12、C33、C44、C66和C1;
所述C11的计算公式为:
所述C12的计算公式为:
所述C33的计算公式为:
所述C44的计算公式为:
所述C13的计算公式为:
所述C66的计算公式为:
所述C1的计算公式为:
S3:对含裂缝的横向各向同性岩石做CT扫描实验,得到裂缝参数;所述的裂缝参数为裂缝密度e、裂缝纵横比α和裂缝与各向同性面的夹角θ;所述的裂缝为倾斜各向同性面硬币状裂缝;
S4:根据步骤S2中得到的刚度矩阵弹性常数和步骤S3中得到的裂缝参数计算得到裂缝的法向柔度ZN、切向柔度ZT、柔度矩阵Z;
法向柔度ZN的计算公式为:
切向柔度ZT的计算公式为:
柔度矩阵Z的计算公式为:
B3的计算公式为:
B4的计算公式为:
S5:根据步骤S4得到的裂缝的柔度ZN、切向柔度ZT和柔度矩阵Z,计算得到含裂缝的横向各向同性岩石的弹性刚度矩阵Ceff;Ceff的计算方法为:先计算含裂缝的横向各向同性岩石的弹性柔度矩阵Meff,对矩阵Meff取逆矩阵,即得到所述的弹性刚度矩阵Ceff;
所述Meff的计算公式为:Meff=M+Z,
所述M的计算公式为:
所述Ceff的计算公式为:Ceff=Meff-1。
S6:根据步骤S5得到的弹性刚度矩阵Ceff的系数,计算声波速度;
声波速度的计算公式为:
其中,
C44 1、C22 1、C33 1、C55 1和C66 1是步骤S5中计算得到的硬币裂缝岩石的弹性刚度矩阵Ceff中的弹性刚度矩阵系数;
ρ是岩石密度;
ψ是地震波相对于X3轴的入射角;
M计算公式为:
同时,本发明还提供了一种上述声波速度计算方法在测量含倾斜裂缝横向各向同性的饱和流体岩石的声波速度中的应用。
一种含倾斜裂缝横向各向同性的饱和流体岩石的声波速度计算方法,包括如下步骤:
S1:测量不含裂缝横向各向同性岩石的实验参数,实验参数包括密度和各向异性声速,所述的声速包括入射方向于各向同性面平行的纵波速度Vp(90°)、入射方向于各向同性面呈45度夹角的纵波速度Vp(45°)、入射方向垂直于各向同性面的纵波速度Vp(0°)、偏振方向于各向同性面平行的横波速度VSH和偏振方向于各向同性面垂直的横波速度VSV。
S2:根据步骤S1测量的实验参数,计算不含裂缝横向各向同性岩石的弹性刚度矩阵系数;弹性刚度矩阵系数包括C11、C13、C12、C33、C44、C66和C1;
所述C11的计算公式为:
所述C12的计算公式为:
所述C33的计算公式为:
所述C44的计算公式为:
所述C13的计算公式为:
所述C66的计算公式为:
所述C1的计算公式为:
S3:对含裂缝的横向各向同性岩石做CT扫描实验,得到裂缝参数;所述的裂缝参数为裂缝密度e、裂缝纵横比α和裂缝与各向同性面的夹角θ;所述的裂缝为倾斜各向同性面硬币状裂缝;
S4:根据步骤S2中得到的弹性刚度矩阵系数和步骤S3中得到的裂缝参数计算得到裂缝的法向柔度ZN、切向柔度ZT、柔度矩阵Z;
法向柔度ZN的计算公式为:
切向柔度ZT的计算公式为:
柔度矩阵Z的计算公式为:
其中,ZN替换为的计算公式为:
其中,
Kf为流体体积模量;
为裂缝孔隙度,计算公式为:
Ks为背景介质的广义体积模量,计算公式为:
e为裂缝密度,计算公式为:
α为裂缝纵横比,a为X1-X2平面上的裂缝半径,n为裂缝数量,V为参考量;
B3的计算公式为:
B4的计算公式为:
S5:根据步骤S4得到的裂缝的法向柔度ZN、切向柔度ZT和柔度矩阵Z,计算得到含裂缝的横向各向同性岩石的弹性刚度矩阵Ceff;Ceff的计算方法为:先计算含裂缝的横向各向同性岩石的弹性柔度矩阵Meff,对矩阵Meff取逆矩阵,即得到所述的弹性刚度矩阵Ceff;
所述Meff的计算公式为:Meff=M+Z,
所述M的计算公式为:
所述Ceff的计算公式为:Ceff=Meff-1。
S6:根据步骤S5得到的弹性刚度矩阵Ceff的系数,计算声波速度。
声波速度的计算公式为:
其中,
C44 1、C22 1、C33 1、C55 1和C66 1是步骤S5中计算得到的硬币裂缝岩石的弹性刚度矩阵Ceff中的弹性刚度矩阵系数;
ρ是岩石密度;
ψ是地震波相对于X3轴的入射角;
M计算公式为:
本次实施例中所用TI背景介质弹性刚度矩阵系数为C11=47.31GPa,C33=33.89GPa,C12=7.83GPa,C13=5.29GPa,C44=17.15GPa;各向同性背景介质弹性刚度矩阵系数为C11=47.31GPa和C12=7.83Gpa;密度2.504g/cm3;裂缝密度和纵横比分别为0.05和0.01;饱和流体体积模量为2.5GPa。
图1为含倾斜裂缝横向各向同性岩石的弹性刚度系数随裂缝倾斜角度的变化图。图2为含饱和流体时含倾斜裂缝横向各向同性岩石的弹性刚度系数随裂缝倾斜角度的变化图。结合附图1和2可知,其中,实心点,为实验测量得到的裂缝垂直各向同性面时的纵波速度;虚线和实线为计算出的弹性刚度矩阵系数,对比发现本发明计算的弹性刚度系数与其精确解基本一致,说明本次提出的模型能够准确计算含倾斜裂缝岩石的等效弹性系数。本发明针对现有针对含裂缝介质的岩石物理模型的弊端,立足于含裂缝岩石的真实特征,基于横向各向同性的背景介质以及与各向同性面斜交且纵横比较小的硬币状裂缝推导出了含倾斜硬币状裂缝的横向各向同性岩石的弹性性质及声速计算模型。对比结果表明,本发明计算得到的等效弹性参数与其精确解吻合较好。本发明的方法能够更加有效的预测含倾斜裂缝岩石的弹性性质,可以为地震勘探及测井识别、评价裂缝储层提供支持。
惟以上所述者,仅为本发明之较佳实施例而已,当不能以此限定本发明实施之范围,即大凡依本发明权利要求及发明说明书所记载的内容所作出简单的等效变化与修饰,皆仍属本发明权利要求所涵盖范围之内。此外,摘要部分和标题仅是用来辅助专利文件搜寻之用,并非用来限制本发明之权利范围。
Claims (8)
1.一种含倾斜裂缝横向各向同性岩石的声波速度计算方法,其特征在于,包括以下步骤:
S1:测量不含裂缝横向各向同性岩石的实验参数,所述的实验参数为密度参数和声速参数;
S2:根据步骤S1测量的实验参数,计算不含裂缝横向各向同性岩石的弹性刚度矩阵系数;
S3:对含裂缝的横向各向同性岩石做CT扫描实验,得到裂缝参数,所述的裂缝参数为裂缝密度e和裂缝纵横比α;
所述的裂缝为倾斜各向同性面硬币状裂缝;
S4:根据步骤S2中得到的弹性刚度矩阵系数和步骤S3中得到的裂缝参数,计算得到裂缝的法向柔度ZN、切向柔度ZT、柔度矩阵Z;
S5:根据步骤S4得到的裂缝的法向柔度ZN、切向柔度ZT和柔度矩阵Z,计算得到含裂缝的横向各向同性岩石的弹性刚度矩阵Ceff;
S6:根据步骤S5得到的弹性刚度矩阵Ceff的系数,计算声波速度。
2.根据权利要求1所述的声波速度计算方法,其特征在于,步骤S1中所述的声速参数为各向异性声速,所述的声速包括Vp(90°)、Vp(45°)、Vp(0°)、VSH和VSV,其中,
Vp(90°)为入射方向与各向同性面平行的纵波速度;
Vp(45°)为入射方向与各向同性面呈45度夹角的纵波速度;
Vp(0°)为入射方向垂直于各向同性面的纵波速度;
VSH为偏振方向与各向同性面平行的横波速度;
VSV为偏振方向与各向同性面垂直的横波速度。
3.根据权利要求2所述的声波速度计算方法,其特征在于,步骤S2中所述的弹性刚度矩阵系数包括C11、C13、C12、C33、C44、C66和C1;
所述C11的计算公式为:
所述C12的计算公式为:
所述C33的计算公式为:
所述C44的计算公式为:
所述C13的计算公式为:
所述C66的计算公式为:
所述C1的计算公式为:
4.根据权利要求3所述的声波速度计算方法,其特征在于,步骤S4中所述的:
法向柔度ZN的计算公式为:
切向柔度ZT的计算公式为:
柔度矩阵Z的计算公式为:
B3的计算公式为:
B4的计算公式为:
5.根据权利要求4所述的声波速度计算方法,其特征在于,步骤S5中所述的弹性刚度矩阵Ceff的计算方法为:对矩阵Meff取逆矩阵,即得到所述的弹性刚度矩阵Ceff,
其中,所述的Meff为含裂缝的横向各向同性岩石的弹性柔度矩阵;
Meff的计算公式为:Meff=M+Z,
所述M的计算公式为:
所述Ceff的计算公式为:
6.根据权利要求5所述的声波速度计算方法,其特征在于,步骤S6中所述的声波速度的计算公式为:
其中,C44 1、C22 1、C33 1、C55 1和C66 1是步骤S5中计算得到的硬币裂缝岩石的弹性刚度矩阵Ceff中的弹性刚度矩阵系数;
ρ是岩石密度;
ψ是地震波相对于X3轴的入射角;
M的计算公式为:
7.一种根据权利要求6所述的声波速度计算方法在测量含倾斜裂缝横向各向同性的饱和流体岩石的声波速度中的应用,其特征在于,将步骤S4中的ZN替换为
所述的的计算公式为:
其中,Kf是流体体积模量,
Ks是含裂缝的横向各向同性岩石的广义体积模量,计算公式为:
8.一种根据权利要求7所述的的声波速度计算方法在测量含裂缝岩石的等效弹性参数和声波速度中的应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811075122.1A CN109459498B (zh) | 2018-09-14 | 2018-09-14 | 一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811075122.1A CN109459498B (zh) | 2018-09-14 | 2018-09-14 | 一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109459498A CN109459498A (zh) | 2019-03-12 |
CN109459498B true CN109459498B (zh) | 2019-10-01 |
Family
ID=65606689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811075122.1A Expired - Fee Related CN109459498B (zh) | 2018-09-14 | 2018-09-14 | 一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109459498B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110515126B (zh) * | 2019-09-12 | 2021-01-01 | 中国石油大学(华东) | 一种含随机分布裂缝横向各向同性岩石的声速计算方法 |
CN111911146B (zh) * | 2020-08-17 | 2023-04-25 | 中国石油大学(华东) | 一种基于裂缝交汇图版的储层压裂效果评价方法 |
CN113504565B (zh) * | 2021-07-09 | 2023-11-17 | 中国石油大学(北京) | 裂缝型储层的地震波场模拟方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103576195A (zh) * | 2013-10-28 | 2014-02-12 | 西北大学 | 一种随压力变化的裂隙介质横波速度预测方法 |
CN104267435A (zh) * | 2014-08-29 | 2015-01-07 | 中国石油天然气股份有限公司 | 一种横观各向同性地层弹性系数的测井计算方法及装置 |
CN106597531A (zh) * | 2015-10-16 | 2017-04-26 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
CN107203005A (zh) * | 2016-03-18 | 2017-09-26 | 中国石油化工股份有限公司 | 一种定量化计算裂缝描述参数的方法 |
CN108399270A (zh) * | 2017-02-08 | 2018-08-14 | 中国石油化工股份有限公司 | 一种确定页岩地层中各向异性泥质比例的方法 |
-
2018
- 2018-09-14 CN CN201811075122.1A patent/CN109459498B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103576195A (zh) * | 2013-10-28 | 2014-02-12 | 西北大学 | 一种随压力变化的裂隙介质横波速度预测方法 |
CN104267435A (zh) * | 2014-08-29 | 2015-01-07 | 中国石油天然气股份有限公司 | 一种横观各向同性地层弹性系数的测井计算方法及装置 |
CN106597531A (zh) * | 2015-10-16 | 2017-04-26 | 中国石油化工股份有限公司 | 含垂直裂缝的页岩的波场传播特征的正演模拟方法 |
CN107203005A (zh) * | 2016-03-18 | 2017-09-26 | 中国石油化工股份有限公司 | 一种定量化计算裂缝描述参数的方法 |
CN108399270A (zh) * | 2017-02-08 | 2018-08-14 | 中国石油化工股份有限公司 | 一种确定页岩地层中各向异性泥质比例的方法 |
Non-Patent Citations (1)
Title |
---|
基于波速测量的龙马溪组页岩的各向异性研究;袁和义 等;《地下空间与工程学报》;20171031;第13卷(第5期);1200-1205 * |
Also Published As
Publication number | Publication date |
---|---|
CN109459498A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106053231B (zh) | 用于真三轴条件下页岩各向异性测试装置及其测试方法 | |
Yee et al. | Elastic and large-strain nonlinear seismic site response from analysis of vertical array recordings | |
RU2602409C2 (ru) | Инверсия анизотропии многоскважинной системы | |
CN109459498B (zh) | 一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 | |
Liu | Effects of fracture aperture and roughness on hydraulic and mechanical properties of rocks: implication of seismic characterization of fractured reservoirs | |
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
CA2878470C (en) | Anisotropy parameter estimation | |
Bakku et al. | Fracture compliance estimation using borehole tube waves | |
Choi et al. | The effect of surface roughness and mixed-mode loading on the stiffness ratio κ x/κ z for fractures | |
CN110515126B (zh) | 一种含随机分布裂缝横向各向同性岩石的声速计算方法 | |
CN106019375B (zh) | 一种页岩气地层层理地球物理评价方法 | |
Matsushima et al. | Ultrasonic measurements of attenuation and velocity of compressional and shear waves in partially frozen unconsolidated sediment and synthetic porous rock | |
CA2273783C (en) | Method of imaging the permeability and fluid content structure within sediment | |
CN113075728B (zh) | 一种建立致密砂岩多尺度三维岩石物理图板的方法 | |
CN109459497B (zh) | 一种含平行裂缝横向各向同性岩石的声波速度计算方法 | |
Hamzehpour et al. | Acoustic wave propagation in heterogeneous two-dimensional fractured porous media | |
He et al. | Calculations of rock matrix modulus based on a linear regression relation | |
Nadri et al. | Estimation of stress-dependent anisotropy from P-wave measurements on a spherical sample | |
CN109709610A (zh) | 一种岩石裂缝探测方法及系统 | |
Dande et al. | The effect of proppant-filled fractures on the elastic properties of 3D-printed rock models | |
Dethof et al. | Numerical investigation of near and far field ultrasound transducer radiation characteristics in concrete | |
CN102809761A (zh) | 一种获取煤岩的速度特征的控制方法 | |
Delaney et al. | Dynamic moduli and attenuation: Rhyolite and carbonate examples | |
Gupta et al. | Statistical AVO intercept-gradient analysis of direct S-waves: A methodology for quantitative fracture characterization | |
Martínez et al. | Investigating anisotropy in rocks by using pulse transmission method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191001 Termination date: 20200914 |
|
CF01 | Termination of patent right due to non-payment of annual fee |