CN109239777B - 一种利用联合反演方法检测构造煤发育的方法 - Google Patents
一种利用联合反演方法检测构造煤发育的方法 Download PDFInfo
- Publication number
- CN109239777B CN109239777B CN201810958361.5A CN201810958361A CN109239777B CN 109239777 B CN109239777 B CN 109239777B CN 201810958361 A CN201810958361 A CN 201810958361A CN 109239777 B CN109239777 B CN 109239777B
- Authority
- CN
- China
- Prior art keywords
- wave
- ava
- model
- gather
- inversion
- 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
- 239000003245 coal Substances 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000011161 development Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 31
- 238000003384 imaging method Methods 0.000 claims description 12
- 230000002194 synthesizing effect Effects 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 abstract description 6
- 238000005755 formation reaction Methods 0.000 description 22
- 230000015572 biosynthetic process Effects 0.000 description 21
- 230000006870 function Effects 0.000 description 4
- 238000013499 data model Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003034 coal gas Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用联合反演方法检测构造煤发育的方法，应用最小二乘法理论，把PP和PS波的AVO道集联合反演成相应的纵横波速度比值，另外利用P‑S波速度和密度测井产生的先验模型在联合反演中得到应用，以提高反演精度。经过模型试验结果表明，本发明方法具有较高的抗噪声能力和良好的构造煤识别能力。
Description
技术领域
本发明涉及地震勘探技术领域，具体涉及一种利用多分量地震数据模型的联合反演方法检测构造煤发育的方法。
背景技术
煤作为一种断裂、易碎的介质，在构造应力作用下，煤很容易破碎变形。在世界各地的煤田中，构造变形煤(TDC)是导致煤矿灾难的因素之一。TDC被认为与煤炭和天然气的突出有关。通常情况下，由于较大的吸附表面，TDC比未变形的煤更容易吸收更多的气体。在中国，很大一部分煤矿是瓦斯高突出的煤矿。在隧道施工前确定TDC的开发和分布将有助于避免煤和瓦斯突出的灾难。地震方法经常用于工程安全评估。在煤田地震勘探中，从地震资料中提取与TDC发育和分布有关的信息已成为一个热点和难点的研究课题。
许多学者根据地震勘探数据来探测TDC进行了大量的研究。对于一些浅煤层，当煤区变形时，可以通过对地震叠加剖面的解释来找到相应的指示。这些预测包括相位异常、振幅变化和同相轴不连续。然而，在煤田地震学中，构造变形被认为是煤层的固有属性，这是很难直接识别出来的。因此，一些学者试图通过地震属性分析来确定煤层中发育异常的地质构造。Ge等人用一种层间地震法来找出煤质中的异常现象。王等人应用极限学习机和主成分分析来预测TDC厚度。其他学者试图通过将地震数据转换成地层岩性的属性数据来确定煤块内发育的TDC。李等人只使用P波叠前反演，推导出了描述TDC的弹性阻抗、声阻抗和孔隙度等评价参数。
当目标煤层随着深度变得越来越薄时，在地震勘探中只使用P波会导致对煤层属性的解释有很大的不确定性。为了弥补传统的P波勘探的不足，煤田多分量地震勘探已经进行了试验。多组分地震勘探是可行的，可为储层预测提供可靠的P-S波速度比(α/β值)。Stewart等人指出，多组分地震的目的是记录和利用反射的P和S波(PP和PS波)，以获得与岩石属性相关的更多信息。在多分量地震方法中，联合PP和PS波的振幅偏移距变化(AVO)反演是一个岩石性质特征的高精度分析方法。与传统的PP波AVO反演相比，在反演中加入PS波将减少多解问题。Veire和landro应用最小二乘理论联合PP和PS反演。Kurt提出了一种通用的线性算法，用于联合PP和PS的AVO反演。Du和Yan应用联合PP和PS的AVO反演来检测流体。Lu等人在联合PP和PS的AVO反演中使用了泰勒展开方法。在煤田勘探领域，Lu等人论证了煤系地层预测的成功案例，其中联合PP和PS的AVO反演方法发挥了重要作用。
然而，可靠的地震反演依赖于超前的岩石物理知识。近年来，对煤的测井曲线和岩石物理特性进行了深入的研究。李等人通过对测井资料的分析，对储层渗透率进行了评价。徐等人系统地研究了TDCs的测井曲线特征。对煤岩物理试验的研究主要集中在微观结构，力学性能和弹性性质。然而，对于地震反演来说，与地震波速度有关的经验关系是最重要的。王等人的岩石物理试验表明，在相似的地层条件下，TDCs的声速和体积密度低于未变形的煤层。然而，对于在不同位置收集的煤炭样品，不管P或S波的速度是很难找到确定的阈值来区分TDC和未变形的煤。P-S波速度和密度之间的关系很差，不变形和变形的煤之间的密度差并不明显。Morcote等人关于动态弹性特性的测试数据表明，与单P或S波速度相比，煤的纵横波速度比(α/β)值受到压力的影响较小。因此，对于TDC检测来说，α/β值是一个有效且稳定的因素。从王等人的测试中可以发现，大部分变形的煤样的α/β值小于1.7。从陈等人的实验数据，对TDC的超声速度进行了分析，得出了相同的经验关系。该文提出了一种基于模型的多分量地震数据联合反演的TDC检测方法。基于模型的反演可以将高频岩石物理模型与地震反演结果联系起来。Mallick应用基于模型的反演将高频测井信息添加到AVO反演结果中。Spikes和Dvorkin使用基于模型的反演来推导岩性和物理参数。然而，关于基于模型的联合PP和PS的AVO反演理论的研究却很少。
发明内容
针对现有技术的不足，本发明旨在提供一种利用多分量地震数据模型的联合反演方法检测构造煤发育的方法，应用最小二乘法理论，把PP和PS波的AVO道集联合反演成相应的α/β值，另外利用P-S波速度和密度测井产生的先验模型在联合反演中得到应用，以提高反演精度。
为了实现上述目的，本发明采用如下技术方案：
一种利用联合反演方法检测构造煤发育的方法，包括如下步骤：
S1、从采集的三分量地震数据中分离出PP波与PS波炮集数据；
S2、对PP波与PS波炮集数据进行叠前预处理；
S3、对叠前预处理之后的PP波与PS波炮集数据进行速度分析，并抽成PP波共成像点道集和PS波共成像点道集；
S4、将PP波共成像点道集与PS波共成像点道集分别转换为PP波AVA道集和PS波AVA道集；
S5、用P波测井曲线与PP波井旁道子波进行合成得到合成的PP波AVA道集，并用合成的PP波AVA道集标定P波测井曲线与步骤S4中得到的PP波AVA道集；
S6、用P波与S波测井曲线与PS波井旁道子波进行合成得到合成的PS波AVA道集，并用来标定P波与S波测井曲线与步骤S4中得到的PS波AVA道集；
S7、根据标定后的PP波与PS波反射时间关系，将步骤S4中得到的PS波AVA道集压缩至PP波零偏移距反射时间上；
S8、将测井数据沿着地震层位外推，插值形成P波速度、S波速度以及密度的初始模型；
S9、分别估计步骤S4中得到的PP波AVA道集与步骤S7中压缩至PP波零偏移距反射时间上的PS波AVA道集的地震子波；
S10、建立联合反演的雅各比矩阵J，并建立目标函数；
S11、高斯-牛顿反演并更新初始模型；若更新后的模型能够满意则输出模型；若更新后的模型不满意，则将该模型作为初始模型，并回到步骤S10；
S12、将更新后输出的模型作为反演模型，计算纵横波速度比模型；
S13、基于纵横波速度比模型预测构造煤发育与分布。
本发明的有益效果在于：经过模型试验结果表明，本发明方法具有较高的抗噪声能力和良好的TDCs识别能力。
附图说明
图1为当P波入射到地下弹性界面时，反射波和透射波的示意图。
图2为PP和PS子波的时间样本图。
图3为本发明实施例提供的用基于模型的PP和PS道集联合反演TDC预测的实施流程图；
图4为含煤地层的理论模型示意图；
图5为合成AVA道集示意图，其中(a)为纯PP波合成AVA道集，(b)为纯PS波合成AVA道集，(c)为加入了10％随机噪声的PP波合成AVA道集，(d)为加入
图6为用纯合成的PP波和PS波AVA道集得出的基于模型的反演结果示意图，其中(a)表示P波速度(α)，(b)表示S波速度(β)，(c)表示密度(ρ)，(d)表示P、S波速度比(α/β)。
图7为加入了10％随机噪声的合成PP波和PS波道集得出的基于模型的反演结果示意图，其中(a)代表P波速度(α)，(b)代表S波速度(β)，(c)代表密度(ρ)，(d)代表P、S波速度比(α/β)。
具体实施方式
以下将结合附图对本发明作进一步的描述，需要说明的是，本实施例以本技术方案的前提，给出了详细的实施方式和具体的操作过程，但本发明的保护范围并不限于本实施例。
以下先对基于模型的联合PP和PS的AVO反演理论作进一步的描述。
1、基于模型的联合反演
在多分量地震勘探中，如图1所示，P波源常用于引发地震波，并在地面上布置多分量检波器，以记录反射波。通常，多分量检波器有三个分量：一个垂直分量(Z)和两个水平分量(X和Y)，反射的PP和PS波是由多分量检波器在地面上同时记录的。由于P波的偏振方向与传播方向平行，所以大部分的PP波的能量都是由Z分量记录的。然而，由于S波的偏振方向与传播方向正交，所以大部分的PS波能量被X和Y分量接收。
对于地下弹性界面(图1)，使用α1，β1，ρ1来表示纵波速度、横波速度和上地层的密度，α2，β2，ρ2表示纵波速度、横波速度、下地层密度，PP和PS反射系数(RPP和RPS)由Zoeppritz方程分别给出：
其中i1表示入射(或反射)P波的角度，j1表示反射的S波角度，i2和j2分别表示透射的P波和S波角度，另外：
其中p为射线参数[41]；而其他参数a，b，c，d，F，H，D没有具体的物理意义。它们只是用来简化RPP和RPS的表达式。在反演时间窗内的所有样本中都使用了P波和S波速度以及密度参数，模型参数矩阵E可以表示为：
E＝(α β ρ)T, (4)
其中：
α＝(α1 α2 L αm), (5)
β＝(β1 β2 L βm), (6)
ρ＝(ρ1 ρ2 L ρm). (7)
矩阵α，β，ρ是由纵波速度、横波速度和密度在时间样品1，2，…m的时间窗内联合反演。基于方程(1)(3)，考虑每个时间样本作为一个界面，分别计算PP和PS波反射系数矩阵RPP和RPS。然后，合成PP波和PS波的振幅随入射角变化(AVA)集合可以推出：
SPP＝WPPRPP, (8)
SPS＝WPSRPS, (9)
其中：
其中,
分别是PP和PS子波事件样品的振幅(图2)。n是子波时间样本的最大数量，应该不到m。在实际应用程序中，PP和PS子波可以通过在频率域中构造振幅和相位谱提取出来。
在提出的反演方法中，第一步是直接计算模型参数矩阵E，根据阻尼最小二乘理论，高斯牛顿公式化简如下：
E(k)＝E(k-1)+ΔE(k), (13)
H＝JTJ, (15)
k表示迭代次数，E(k)是第k次迭代之后最新的模型参数矩阵。ΔE(k)是第k次模型的修改矩阵：
ΔE(k)＝(Δα(k) Δβ(k) Δρ(k))T. (16)
J和H分别是Jacobian和Hessian矩阵[45，46]。在联合AVO反演中，如果PP和PS波AVA道集的入射角(在等式(1)和(2)中的i1)是θ1，θ2，…θx，则Jacobian矩阵J形式如下：
在等式(12)中，λ是一个阻尼参数，I是一个单位矩阵。这两个参数用在稀疏最小二乘法问题中以确保总是存在一个唯一解。一开始给出先验模型矩阵E(0)，则第k次迭代之后的矩阵E(k)和ΔE(k)就可以导出来。通过控制迭代次数，建立了联合AVO反演的目标函数如下：
其中
和
分别是输入的PP和PS波AVA道集；
和
分别是第k次迭代之后合成的PP和PS波AVA道集。权重因子ω大小从0到1。如果ω大于0.5，地震数据更多的涉及到约束，否则，先验模型对约束的控制更大。当Q(k)达到合适的最小值或者几乎不变时，迭代就可以停止了，与此相对应得E(k)则是输出的反演模型矩阵。
基于上述基于模型的联合PP和PS的AVO反演理论，本实施例提供了一种利用多分量地震数据模型的联合反演方法检测构造煤发育的方法,包括如下步骤(如图3所示)：
S1、从采集的三分量地震数据中分离出PP波与PS波炮集数据；
S2、对PP波与PS波炮集数据进行叠前预处理；
S3、对叠前预处理之后的PP波与PS波炮集数据进行速度分析，并抽成PP波共成像点道集和PS波共成像点道集；
S4、将PP波共成像点道集与PS波共成像点道集分别转换为PP波AVA道集和PS波AVA道集；
S5、用P波测井曲线与PP波井旁道子波进行合成得到合成的PP波AVA道集，并用合成的PP波AVA道集标定P波测井曲线与步骤S4中得到的PP波AVA道集；
S6、用P波与S波测井曲线与PS波井旁道子波进行合成得到合成的PS波AVA道集，并用来标定P波与S波测井曲线与步骤S4中得到的PS波AVA道集；
S7、根据标定后的PP波与PS波反射时间关系，将步骤S4中得到的PS波AVA道集压缩至PP波零偏移距反射时间上；
S8、将测井数据沿着地震层位外推，插值形成P波速度、S波速度以及密度的初始模型；
S9、分别估计步骤S4中得到的PP波AVA道集与步骤S7中压缩至PP波零偏移距反射时间上的PS波AVA道集的地震子波；
S10、建立联合反演的雅各比矩阵J，并建立目标函数；
S11、高斯-牛顿反演并更新初始模型；若更新后的模型能够满意则输出模型；若更新后的模型不满意，则将该模型作为初始模型，并回到步骤S10；
S12、将更新后输出的模型作为反演模型，计算纵横波速度比模型；
S13、基于纵横波速度比模型预测构造煤发育与分布。
以下通过模型试验对本实施例方法的性能作进一步的说明。
本试验中使用了一个合成数据集来证明本实施例所提出的TDC检测方法的有效性。含煤地层的理论模型如图4所示，其中3-5层为煤层。从表1中可以看到层5是TDC煤层，其弹性参数远远低于未变形的煤层(4和6层)。5层的α/β值小于周围的岩石。在PP反射时域内，分别由地震子波(主频率30Hz)矩阵和PP和PS反射系数矩阵的乘积合成纯PP波AVA道集(图5a)和PS波AVA道集(图5b)。对应的入射角度为5°到30°，间隔为5°。然后用10％的随机噪声(图5c，d)分别加入纯AVA道集。
表1：含煤地层的理论模型参数
如图6所示，虽然煤层很薄，初始模型与实际模型相差很大，但是基于模型的联合反演在无噪声条件下也能得到准确的结果。变形煤在任何反演模型上都很容易识别。将10％的随机噪声添加到合成的AVA道集中，如图7所示，在反演模型中有很多异常。然而，反演出的α/β模型(图7d)仍比其他反演模型稳定(图7c)，变形煤层仍是清晰的。
本实施例提出了一种基于模型的多分量地震数据联合反演的TDC检测方法。α/β值是用来区分TDCs和未变形的煤。在反演中同时使用了声波和测井资料。因此，测井数据揭示的TDC发育的先验知识可以添加到反演结果中。模型试验结果表明，基于模型的联合反演能有效地识别出在噪声条件下的煤层内的TDC。
专业术语：
TDC构造变形煤
P-wave压缩波/纵波
S-wave剪切波/横波
PP-wave转换纵波
PS-wave转换横波
α地层中纵波速度
α1上部地层纵波速度
α2下部地层纵波速度
β地层中横波速度
β1上部地层横波速度
β2下部地层横波速度
α/β纵横波速度比
ρ地层密度
ρ1上部地层的密度
ρ2下部地层的密度
i1纵波入射角
j1横波反射角
i2纵波透射角
j2横波透射角
p射线参数
RPP转换纵波的反射系数
RPS转换横波的反射系数
RPP转换纵波的反射系数矩阵
RPS转换横波的反射系数矩阵
α地层的纵波速度矩阵
β地层的横波速度矩阵
ρ地层的密度矩阵
E地层的模型参数矩阵
ΔE地层的模型矫正矩阵
E(0)地层的先验模型矩阵
E(k)地层第k次迭代之后的矫正模型矩阵
WPP转换纵波子波矩阵
WPS转换横波子波矩阵
J雅克比矩阵
H海塞矩阵
Q(k)第k次迭代的目标函数
ω权重系数
对于本领域的技术人员来说，可以根据以上的技术方案和构思，作出各种相应的改变和变形，而所有的这些改变和变形都应该包括在本发明权利要求的保护范围之内。
Claims (1)
1.一种利用联合反演方法检测构造煤发育的方法，其特征在于，包括如下步骤：
S1、从采集的三分量地震数据中分离出PP波与PS波炮集数据；
S2、对PP波与PS波炮集数据进行叠前预处理；
S3、对叠前预处理之后的PP波与PS波炮集数据进行速度分析，并抽成PP波共成像点道集和PS波共成像点道集；
S4、将PP波共成像点道集与PS波共成像点道集分别转换为PP波AVA道集和PS波AVA道集；
S5、用P波测井曲线与PP波井旁道子波进行合成得到合成的PP波AVA道集，并用合成的PP波AVA道集标定P波测井曲线与步骤S4中得到的PP波AVA道集；
S6、用P波与S波测井曲线与PS波井旁道子波进行合成得到合成的PS波AVA道集，并用来标定P波与S波测井曲线与步骤S4中得到的PS波AVA道集；
S7、根据标定后的PP波与PS波反射时间关系，将步骤S4中得到的PS波AVA道集压缩至PP波零偏移距反射时间上；
S8、将测井数据沿着地震层位外推，插值形成P波速度、S波速度以及密度的初始模型；
S9、分别估计步骤S4中得到的PP波AVA道集与步骤S7中压缩至PP波零偏移距反射时间上的PS波AVA道集的地震子波；
S10、建立联合反演的雅各比矩阵J，并建立目标函数；其中，J形式如下：
其中，矩阵α，β，ρ分别是由纵波速度、横波速度和密度在时间样本1，2，…m的时间窗内联合反演；PP和PS波AVA道集的入射角是θ1，θ2，…θx；
SPP＝WPPRPP，
SPS＝WPSRPS，
其中：
WPP和WPS分别是PP和PS子波矩阵，它们的形式是：
其中，
分别是PP或PS子波事件样本的振幅，n是子波时间样本的最大数量；
联合AVO反演的目标函数如下：
其中，
和
分别是输入的PP和PS波AVA道集；
和
分别是第k次迭代之后合成的PP和PS波AVA道集；权重因子ω大小从0到1；E(k)是第k次迭代之后最新的模型参数矩阵；E(0)是先验模型矩阵；
S11、高斯-牛顿反演并更新初始模型；若更新后的模型能够满意则输出模型；若更新后的模型不满意，则将该模型作为初始模型，并回到步骤S10；
S12、将更新后输出的模型作为反演模型，计算纵横波速度比模型；
S13、基于纵横波速度比模型预测构造煤发育与分布。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810958361.5A CN109239777B (zh) | 2018-08-22 | 2018-08-22 | 一种利用联合反演方法检测构造煤发育的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810958361.5A CN109239777B (zh) | 2018-08-22 | 2018-08-22 | 一种利用联合反演方法检测构造煤发育的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109239777A CN109239777A (zh) | 2019-01-18 |
| CN109239777B true CN109239777B (zh) | 2020-05-29 |
Family
ID=65068231
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810958361.5A Active CN109239777B (zh) | 2018-08-22 | 2018-08-22 | 一种利用联合反演方法检测构造煤发育的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109239777B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117055115B (zh) * | 2023-10-11 | 2023-12-26 | 煤炭科学研究总院有限公司 | 探测煤岩体结构异常区的方法、装置、设备和介质 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2006262697B2 (en) * | 2005-06-24 | 2010-12-23 | Exxonmobil Upstream Research Company | Method for obtaining porosity and shale volume from seismic data |
| CN102359924B (zh) * | 2011-09-19 | 2013-03-06 | 中国地质大学（北京） | 一种基于多波地震数据的煤岩强度的检测方法 |
| CA2957818C (en) * | 2014-08-19 | 2023-02-21 | Cgg Services Sa | Joint inversion of compressional and shear seismic data in native time domains |
| CN104965224B (zh) * | 2015-06-03 | 2017-10-27 | 北京多分量地震技术研究院 | 用平均入射角道集进行pp波与ps波联合avo反演方法 |
| CN106886045A (zh) * | 2015-12-15 | 2017-06-23 | 师素珍 | 一种基于叠前同步反演进行构造煤分布预测的方法 |
-
2018
- 2018-08-22 CN CN201810958361.5A patent/CN109239777B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109239777A (zh) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101329405B (zh) | 一种简单的多参数地震反演方法 | |
| Piatti et al. | Constrained 1D joint inversion of seismic surface waves and P‐refraction traveltimes | |
| WO2017035104A1 (en) | Velocity model seismic static correction | |
| US20140336940A1 (en) | Estimation of q-factor in time domain | |
| CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
| CN105938203A (zh) | 一种储层特性的检测方法及装置 | |
| Lynn et al. | Relationship of P-wave seismic attributes, azimuthal anisotropy, and commercial gas pay in 3-DP-wave multiazimuth data, Rulison Field, Piceance Basin, Colorado | |
| Diaz-Acosta et al. | Investigation of fractured carbonate reservoirs by applying shear-wave splitting concept | |
| Faraji et al. | Integration of core data, well logs and seismic attributes for identification of the low reservoir quality units with unswept gas in the carbonate rocks of the world’s largest gas field | |
| Mari | Well seismic surveying and acoustic logging | |
| CN104422955A (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
| CN109239777B (zh) | 一种利用联合反演方法检测构造煤发育的方法 | |
| Herwanger et al. | Anisotropic velocity changes in seismic time-lapse data | |
| Bates et al. | The Seismic Evaluation of a Naturally Fractured Tight-gas Sand Reservoir in the Wind River Basin, Wyoming | |
| Zong | Elastic properties of rock salt: Lab measurements, numerical and physical modeling, and Gulf of Mexico well logs and VSP | |
| Huang et al. | The in-situ stress prediction of a fractured shale reservoir based on AVAZ inversion: A case study from Southwest China | |
| Li et al. | Seismic attenuation and velocities of P-and S-waves in the German KTB area | |
| Suroso et al. | Hydrocarbon identification by evaluating anisotropy parameters estimated from crosswell seismic data | |
| Masri | Amplitude Versus Offset analysis–A possible useful tool for geothermal exploration | |
| Xu | Oil reservoir assessment using multicomponent seismic data | |
| Boonyasatphan | Reservoir characterization for unconventional resource potential, Pitsanulok basin, onshore Thailand | |
| Ampomah et al. | Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques | |
| Mur et al. | Tracking Subsurface Supercritical CO 2 Using Advanced Reflection Seismic and Well Log‐Based Workflows Incorporating Fluid Density and Pore Pressure Effects: Relevance to Reservoir Monitoring and CO 2 EOR | |
| Gordon et al. | Characterizing multiple formations for gas production and CO2 sequestration | |
| Zong et al. | Examining the significance of the converted-modes beneath an ultra-high-velocity caprock: Physical modeling and a Gulf of Mexico salt proximity VSP survey |
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 |