CN106556861B - 一种基于全方位地震资料的方位avo反演方法 - Google Patents
一种基于全方位地震资料的方位avo反演方法 Download PDFInfo
- Publication number
- CN106556861B CN106556861B CN201510618215.4A CN201510618215A CN106556861B CN 106556861 B CN106556861 B CN 106556861B CN 201510618215 A CN201510618215 A CN 201510618215A CN 106556861 B CN106556861 B CN 106556861B
- Authority
- CN
- China
- Prior art keywords
- azimuth
- objective function
- angles
- gradient
- omnibearing
- 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 24
- 239000011159 matrix material Substances 0.000 claims abstract description 32
- 230000035945 sensitivity Effects 0.000 claims abstract description 25
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 3
- 238000011161 development Methods 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 14
- 206010017076 Fracture Diseases 0.000 claims description 12
- 208000010392 Bone Fractures Diseases 0.000 claims description 10
- 238000002939 conjugate gradient method Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 3
- 239000002775 capsule Substances 0.000 abstract 1
- 238000001514 detection method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于全方位地震资料的方位AVO反演方法，包括：基于方位角度道集形式的全方位地震资料，计算不同方位角及入射角下的反射系数；将所述反射系数与实测的方位角度道集建立目标函数；推导目标函数的敏感度矩阵；基于目标函数和敏感度矩阵，利用最优化算法求解方位AVO反演。实现了从全方位角度道集反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向，为全方位地震资料的应用及预测裂缝分布提供重要资料。
Description
技术领域
本发明涉及石油地球物理勘探技术领域，具体地涉及基于全方位地震资料的方位AVO反演。
背景技术
裂缝分布和发育信息对石油和天然气的勘探及开发有重要指导作用，特别是碳酸盐岩储层的裂缝、孔隙和孔洞是油气的主要富集部位和运移通道，是碳酸盐岩油气成藏的必要条件之一。为了有效的识别裂缝，地球物理工作者提出了各种裂缝检测方法，比如基于叠后地震资料的应力场分析与不连续性检测、基于叠前地震资料的AVO与各向异性分析。在叠后资料应用方面，基于应力场分析的裂缝预测适用于构造引起的裂缝，基于不连续检测的各种裂缝预测方法，如相干分析、广义希尔伯特变换和基于图像边缘检测等，能够有效的圈定不规律的缝洞发育带与均匀介质的地震波场之间的界限，但受噪音影响较大。随着宽方位(或者全方位)地震勘探的应用，如何从全方位地震资料中挖掘更加丰富的地下介质信息成为一个重要的研究热点。
由于裂缝性储层常存在明显的HTI各向异性的影响，主要表现为振幅、速度、反射波形及相位随测线方位变化而变化，因此在宽方位叠前资料应用中可以利用纵波信息随方位的变化来预测裂缝。比如利用方位旅行时、方位衰减、方位速度(VVAZ，速度随方位角变化)和方位振幅变化(AVAZ，振幅随方位角变化)等来预测裂缝。在利用纵波方位各向异性特征预测裂缝方面，我国学者对此开展了较多的应用研究，如程冰洁等基于宽方位转换波地震资料利用转换波AVAZ裂缝检测方法对川西深层裂缝性气藏进行了裂缝检测；孙炜对用于纵波方位各向异性裂缝预测的地震属性优选进行了研究，并利用成像测井及岩心资料上的裂缝信息验证了预测结果的可靠性和准确性。目前的方法是分别对不同方位角下的地震数据进行属性提取或者地震反演，然后对不同方位角度下的属性进行各向异性椭圆拟合，由椭圆的长、短轴的比值和长短轴的方向来表征裂缝发育强度和方向，从而预测储层裂缝的分布特征。其中用于椭圆拟合的地震属性可以是AVO属性，也可是波阻抗等。这类基于椭圆拟合的方法受到椭圆拟合算法精度的影响很大，并且抗噪性低。
发明内容
本发明直接从全方位地震资料出发，以全方位地震资料为基础，综合利用各个方位角下的地震资料，利用P波方位各向异性分析原理进行方位AVO(AVAZ)反演，从全方位地震资料中反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向，为裂缝型储层预测提供重要分析技术手段。反演得到的各向异性梯度及裂缝方向可直接表征裂缝发育强度和方向。
本发明提供了一种基于全方位地震资料的方位AVO反演方法，包括：基于方位角度道集形式的全方位地震资料，计算不同方位角及入射角下的反射系数；将所述反射系数与实测的方位角度道集建立目标函数；推导目标函数的敏感度矩阵；基于目标函数和敏感度矩阵，利用最优化算法求解方位AVO反演。
进一步地，由P波方位各向异性分析原理中的Corrigan-Ruger近似方程推导出方位各向异性线性近似方程(1)，计算不同方位角及入射角下的反射系数：
所述公式(1)中，θ是入射角、φ是方位角、表示方位角为φ入射角为θ时的反射系数、A、B、C、D是Corrigan-Ruger近似方程中的临时变量，四个临时变量与纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向的关系为：
所述公式(2)中，Aiso是纵波阻抗反射系数、Biso是各项同性梯度、Bani是各向异性梯度、φsym是HTI介质对称轴方向。
进一步地，上述公式(1)为线性的，对于不同方位角和入射角的情况简写为矩阵形式：
d＝Gm (3)
其中，
所述公式(4)中，m是由临时变量构成的向量，G是由方位角和入射角构成的正演矩阵，d是不同方位角和入射角下的振幅值，R(φi,θj)为方位角为φi入射角为θj时的振幅值，M为方位角个数，N为入射角个数。
进一步地，所述将所述反射系数与实测的方位角度道集建立目标函数的步骤，利用公式(3)正演计算出不同方位角和入射角下的反射系数，并与相应的方位角和入射角下的观测数据求取误差来构建反演问题的目标函数：
式中，ds为不同方位角下的观测数据，即由实测全方位地震资料经过处理得到的方位角度道集，R(φi,θj)是方位角为φi，入射角为θj时的振幅值。
进一步地，所述推导目标函数的敏感度矩阵的步骤，包括对目标函数J(m)求取偏导数：
所述公式(6)中，为目标函数J(m)对临时变量A、B、C、D的偏导数，由此组成目标函数的敏感度矩阵即
进一步地，利用最优化算法中的共轭梯度法求解方位AVO反演。利用最优化算法中的共轭梯度法迭代求解，直到满足终止条件为止。
优选地，所述终止条件是目标函数误差达到预定值或目标函数敏感度矩阵范数达到预定值。
优选地，目标函数误差的预定值为0.1，目标函数敏感度矩阵范数的预定值为0.1。
本发明利用P波方位各向异性分析原理进行AVAZ反演，从全方位地震资料中反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向，反演得到的各向异性梯度及裂缝方向可直接表征裂缝发育强度和方向，为裂缝型储层预测提供重要分析技术手段，提高了裂缝型储层预测的效率和准确性。
本发明的方法是基于P波方位各向异性分析原理中的Corrigan-Ruger近似方程，综合的利用全方位地震资料，从全方位地震资料中反演得到反映裂缝发育强度和方向的数据，当某一方位较其他方位数据存在较大误差时算法也能够得到较为稳定的结果，抗噪能力强。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
图1是本发明方法的流程图。
图2是实施例1中的方位角度道集，作为输入数据给定。
图3a-d是实施例1中方位AVO(AVAZ)反演结果，依次为P波反射系数、各向同性梯度、各向异性梯度、方向。
图4是实施例1中方位AVO(AVAZ)反演结果的沿层切片，Aiso为P波反射系数、Biso为各向同性梯度、Bani为各向异性梯度。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明公开一种基于全方位地震资料的方位AVO反演方法，属于石油地球物理勘探技术领域。首先利用Corrigan-Ruger近似方程计算不同方位角及入射角下的反射系数，继而与方位角度道集建立目标函数，然后推导了目标函数的敏感度矩阵，最后利用最优化算法求解方位AVO反演(AVAZ)反演问题。实现了从全方位角度道集反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向，为全方位地震资料的应用及预测裂缝分布提供重要资料。
具体地，如图1所示，本发明提供了一种基于全方位地震资料的方位AVO反演方法，包括：基于方位角度道集形式的全方位地震资料，计算不同方位角及入射角下的反射系数；将所述反射系数与实测的方位角度道集建立目标函数；推导目标函数的敏感度矩阵；基于目标函数和敏感度矩阵，利用最优化算法求解方位AVO反演。
首先，基于方位角度道集形式的全方位地震资料，计算不同方位角及入射角下的反射系数。具体地，可由P波方位各向异性分析原理中的Corrigan-Ruger近似方程推导出方位各向异性线性近似方程，方位各向异性线性近似方程具体为公式(1)所示：
所述公式(1)中，方位各向异性线性近似方程能够计算出不同方位角及入射角下的反射系数，θ是入射角、φ是方位角、表示方位角为φ入射角为θ时的反射系数、A、B、C、D是Corrigan-Ruger近似方程中的临时变量，这四个临时变量与纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向的关系为：
所述公式(2)中，Aiso是纵波阻抗反射系数、Biso是各项同性梯度、Bani是各向异性梯度、φsym是HTI介质对称轴方向。
优选地，上述公式(1)为线性的，对于不同方位角和入射角的情况可简写为矩阵形式：
d＝Gm (3)
其中，
所述公式(4)中，m是由临时变量构成的向量，G是由方位角和入射角构成的正演矩阵，d是不同方位角和入射角下的振幅值，R(φi,θj)为方位角为φi入射角为θj时的振幅值，方位角个数为M，入射角个数为N。
接下来，将所述反射系数与实测的方位角度道集建立目标函数。具体地，利用公式(3)正演计算出不同方位角和入射角下的反射系数，并与相应的方位角和入射角下的观测数据求取误差来构建反演问题的目标函数，具体由公式(5)：
式中，ds为不同方位角下的观测数据，即由实测全方位地震资料经过处理得到的方位角度道集，m是由临时变量A、B、C、D构成的向量，即m＝(A,B,C,D)，G是由方位角和入射角构成的正演矩阵，具体形式见上部步骤中公式(4)。R(φi,θj)是方位角为φi，入射角为θj时的振幅值。
接下来，推导目标函数的敏感度矩阵。对反演问题的目标函数J(m)求取偏导数，详见公式(6)：
所述公式(6)中，为目标函数J(m)对临时变量A、B、C、D的偏导数，由此组成目标函数的敏感度矩阵即
最后，基于目标函数和敏感度矩阵，利用最优化算法求解方位AVO反演。可利用最优化算法中的共轭梯度法求解。具体地，基于目标函数和敏感度矩阵，利用最优化算法中的共轭梯度法迭代求解，直到满足终止条件为止，终止条件是目标函数误差达到较小的值(例如0.1)或目标函数敏感度矩阵范数达到较小的值(例如0.1)。
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
下面以某地区的实际地震资料为例，利用本发明方法从全方位地震资料中反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向，进而说明本发明的效果。
图2是实施例1中的方位角度道集，作为输入数据给定。图中黑色方框内为一个方位角度道集，由若干角度道集组成，每个角度道集所对应的方位角不同。图中方框的方位角道集由5个方位角(分别为0、36、72、108、144度)，每个方位角对应一组角度道集(入射角度道集，图中入射角为0度到45度)。
图3a-d是实施例1中方位AVO(AVAZ)反演结果，以图2中方位角度道集为输入数据，进行方位AVO(AVAZ)反演，得到图3a中P波反射系数、图3b中各向同性梯度、图3c中各向异性梯度、图3d中方向。
图4是实施例1中方位AVO(AVAZ)反演结果的沿层切片，Aiso为P波反射系数、Biso为各向同性梯度、Bani为各向异性梯度。
实施例1
本实施例是本发明用于从全方位地震资料中反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向的实施实例，一种基于全方位地震资料的方位AVO反演方法的步骤如下：
本发明中以方位角度道集的形式给定全方位地震资料，给定临时变量(A、B、C、D)的初始模型，然后利用方位各向异性线性近似方程计算不同方位角及入射角下的反射系数；
将所述反射系数与实测的方位角度道集建立目标函数；
推导目标函数的敏感度矩阵；
基于目标函数和敏感度矩阵，利用最优化算法中的共轭梯度法迭代求解上述步骤，直到满足终止条件为止，终止条件是目标函数误差达到较小的值(如0.1)或目标函数敏感度矩阵范数达到较小的值(如0.1)。
本发明直接从全方位地震资料出发，以全方位地震资料为基础，综合利用各个方位角下的地震资料，利用P波方位各向异性分析原理进行AVAZ反演，从全方位地震资料中反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向，为裂缝型储层预测提供重要分析技术手段。反演得到的各向异性梯度及裂缝方向可直接表征裂缝发育强度和方向。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (6)
1.一种基于全方位地震资料的方位AVO反演方法，包括：
基于方位角度道集形式的全方位地震资料，计算不同方位角及入射角下的反射系数；
将所述反射系数与实测的方位角度道集建立目标函数；
推导目标函数的敏感度矩阵；
基于目标函数和敏感度矩阵，利用最优化算法求解方位AVO反演；
其中由P波方位各向异性分析原理中的Corrigan-Ruger近似方程推导出方位各向异性线性近似方程(1)，计算不同方位角及入射角下的反射系数：
所述公式(1)中，θ是入射角、φ是方位角、表示方位角为φ入射角为θ时的反射系数、A、B、C、D是Corrigan-Ruger近似方程中的临时变量，四个临时变量与纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向的关系为：
所述公式(2)中，Aiso是纵波阻抗反射系数、Biso是各项同性梯度、Bani是各向异性梯度、φsym是HTI介质对称轴方向；
其中上述公式(1)为线性的，对于不同方位角和入射角的情况简写为矩阵形式：
d＝Gm (3)
其中，
所述公式(4)中，m是由临时变量构成的向量，G是由方位角和入射角构成的正演矩阵，d是不同方位角和入射角下的振幅值，R(φi,θj)为方位角为φi入射角为θj时的振幅值，M为方位角个数，N为入射角个数；其中所述将所述反射系数与实测的方位角度道集建立目标函数的步骤，利用公式(3)正演计算出不同方位角和入射角下的反射系数，并与相应的方位角和入射角下的观测数据求取误差来构建反演问题的目标函数：
式中，ds为不同方位角下的观测数据，即由实测全方位地震资料经过处理得到的方位角度道集，R(φi,θj)是方位角为φi，入射角为θj时的振幅值，dsij为ds中第i行第j列的数据，i为ds的行数，j为ds的列数；
其中所述推导目标函数的敏感度矩阵的步骤，包括对目标函数J(m)求取偏导数：
所述公式(6)中，为目标函数J(m)对临时变量A、B、C、D的偏导数，由此组成目标函数的敏感度矩阵即
2.根据权利要求1所述的方位AVO反演方法，其中利用最优化算法中的共轭梯度法求解方位AVO反演。
3.根据权利要求2所述的方位AVO反演方法，其中利用最优化算法中的共轭梯度法迭代求解，直到满足终止条件为止。
4.根据权利要求3所述的方位AVO反演方法，其中所述终止条件是目标函数误差达到预定值或目标函数敏感度矩阵范数达到预定值。
5.根据权利要求4所述的方位AVO反演方法，其中目标函数误差的预定值为0.1，目标函数敏感度矩阵范数的预定值为0.1。
6.根据权利要求1所述的方位AVO反演方法，其中反演得到的各向异性梯度及裂缝方向表征裂缝发育强度和方向。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510618215.4A CN106556861B (zh) | 2015-09-24 | 2015-09-24 | 一种基于全方位地震资料的方位avo反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510618215.4A CN106556861B (zh) | 2015-09-24 | 2015-09-24 | 一种基于全方位地震资料的方位avo反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106556861A CN106556861A (zh) | 2017-04-05 |
| CN106556861B true CN106556861B (zh) | 2018-11-13 |
Family
ID=58413820
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510618215.4A Active CN106556861B (zh) | 2015-09-24 | 2015-09-24 | 一种基于全方位地震资料的方位avo反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106556861B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107092029B (zh) * | 2017-04-26 | 2019-01-18 | 中国石油大学(北京) | 一种地震反演方法和装置 |
| CN109655905A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 基于全方位角数据叠前avaz反演裂缝预测方法及系统 |
| CN107918151B (zh) * | 2017-10-19 | 2019-10-11 | 中国石油天然气集团公司 | 基于敏感方位地震数据的断裂识别方法及装置 |
| CN109143352B (zh) * | 2018-07-04 | 2019-10-01 | 中国石油大学(华东) | 一种各向异性介质地震反射特征方程建立方法 |
| CN109471171B (zh) * | 2018-09-21 | 2020-06-09 | 中国石油天然气集团有限公司 | 一种混叠地震数据分离的方法、装置及系统 |
| CN111352154B (zh) * | 2018-12-20 | 2021-12-28 | 中国石油化工股份有限公司 | 一种基于宽方位地震的储层预测方法 |
| CN113075727B (zh) * | 2021-02-20 | 2023-09-26 | 中国石油天然气股份有限公司 | 基于振幅方位导数的储层裂缝分布特征确定方法及系统 |
| CN113933898B (zh) * | 2021-09-23 | 2022-07-19 | 中国地质大学（武汉） | 裂缝储层各向异性特征识别方法、装置、设备及存储介质 |
| CN114002738B (zh) * | 2021-11-02 | 2023-08-15 | 中海石油(中国)有限公司 | 高角度断裂约束的叠前裂缝预测方法、装置及存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102004263A (zh) * | 2010-09-17 | 2011-04-06 | 中国石油大学(华东) | 一种基于角度部分叠加地震数据的流体识别的方法 |
| CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
| CN102736103A (zh) * | 2011-12-23 | 2012-10-17 | 电子科技大学 | 基于角度梯度弹性阻抗的储层预测方法 |
| CN104005760A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 基于方位各向异性弹性阻抗的裂缝检测方法 |
| EP2846175A2 (en) * | 2013-09-06 | 2015-03-11 | Services Petroliers Schlumberger | Seismic survey analysis |
-
2015
- 2015-09-24 CN CN201510618215.4A patent/CN106556861B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102004263A (zh) * | 2010-09-17 | 2011-04-06 | 中国石油大学(华东) | 一种基于角度部分叠加地震数据的流体识别的方法 |
| CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
| CN102736103A (zh) * | 2011-12-23 | 2012-10-17 | 电子科技大学 | 基于角度梯度弹性阻抗的储层预测方法 |
| EP2846175A2 (en) * | 2013-09-06 | 2015-03-11 | Services Petroliers Schlumberger | Seismic survey analysis |
| CN104005760A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 基于方位各向异性弹性阻抗的裂缝检测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 裂缝介质中多方位AVO特征分析技术及应用举例;莫午零 等;《天然气地球科学》;20071231;第18卷(第6期);第813-818页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106556861A (zh) | 2017-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106556861B (zh) | 一种基于全方位地震资料的方位avo反演方法 | |
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| KR102021276B1 (ko) | 진폭 보존을 갖는 fwi 모델 도메인 각도 스택들 | |
| US9869783B2 (en) | Structure tensor constrained tomographic velocity analysis | |
| AU2014227491B2 (en) | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields | |
| US7952960B2 (en) | Seismic imaging with natural Green's functions derived from VSP data | |
| US20150301208A1 (en) | Seismic data processing | |
| Haldorsen et al. | Locating microseismic sources using migration-based deconvolution | |
| US20150168574A1 (en) | Seismic trace attribute | |
| CN103487831B (zh) | Avo地震正演计算方法 | |
| US10215869B2 (en) | System and method of estimating anisotropy properties of geological formations using a self-adjoint pseudoacoustic wave propagator | |
| US20140200816A1 (en) | Seismic data processing | |
| US9658354B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
| Thiel et al. | Comparison of acoustic and elastic full‐waveform inversion of 2D towed‐streamer data in the presence of salt | |
| US9921324B2 (en) | Systems and methods employing upward beam propagation for target-oriented seismic imaging | |
| Wang et al. | Full‐waveform inversion of high‐frequency teleseismic body waves based on multiple plane‐wave incidence: Methods and practical applications | |
| CN109143352B (zh) | 一种各向异性介质地震反射特征方程建立方法 | |
| Li et al. | Gaussian beam imaging of fractures near the wellbore using sonic logging tools after removing dispersive borehole waves | |
| Karalliyadda et al. | S‐wave splitting in the offshore S outh I sland, N ew Z ealand: Insights into plate‐boundary deformation | |
| Nakayama et al. | Joint Processing of Multi-Well and Multi-Fiber Walkaway DAS-VSP Datasets from a CCUS Pilot Field Onshore Japan | |
| CN108802822B (zh) | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 | |
| Kimura et al. | Stress Field Estimation From S‐Wave Anisotropy Observed in Multi‐Azimuth Seismic Survey With Cabled Seafloor Seismometers Above the Nankai Trough Megathrust Zone, Japan | |
| US10739483B2 (en) | Absolute strength and absolute sensitivity in seismic data | |
| EP3087417B1 (en) | Methods and systems of joint inversion determination of surface-consistent amplitude corrections | |
| GB2519856A (en) | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
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 |