CN113253345A - 相位差异属性获取方法及系统 - Google Patents
相位差异属性获取方法及系统 Download PDFInfo
- Publication number
- CN113253345A CN113253345A CN202010084606.3A CN202010084606A CN113253345A CN 113253345 A CN113253345 A CN 113253345A CN 202010084606 A CN202010084606 A CN 202010084606A CN 113253345 A CN113253345 A CN 113253345A
- Authority
- CN
- China
- Prior art keywords
- phase difference
- data volume
- data
- phase
- window
- 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
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000004364 calculation method Methods 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000012935 Averaging Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 11
- 238000004891 communication Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000001427 coherent effect Effects 0.000 description 6
- 239000000872 buffer Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 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
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种相位差异属性获取方法及系统，所述方法包含：获取叠后地震数据体；对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性。
Description
技术领域
本发明涉及石油天然气地震勘探领域，尤指一种相位差异属性获取方法及系统。
背景技术
瞬时相位属性是由复地震道技术衍生出的“三瞬属性”之一，与振幅相比，相位属性不受反射系数大小的限制，即相位的形态与反射强弱无关，只与局部构造变化有关。因此，在发育构造的区域，相位往往会伴有突变，且微弱构造同样可以利用相位属性有效体现。
现有技术中，1.根据叠后数据进行的相干运算计算复杂，效率较低，构造边界模糊，分辨率较差；2.对于不同走向的构造表征能力有限，只能体现单一方向(与计算方向垂直)的构造；3.对于倾斜界面或起伏波动较剧烈的界面容易产生构造假象。
发明内容
本发明目的在于一种相位差异属性获取方法及系统，利用由相邻道瞬时相位并结合构造导向思想求取的相位差异属性精细表征和识别地下地质构造。
为达上述目的，本发明所提供的相位差异属性获取方法具体包含：获取叠后地震数据体；对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性。
在上述相位差异属性获取方法中，优选的，所述预设方向包含正北方向、北东南西向、东西方向、北西南东方向四个方向。
在上述相位差异属性获取方法中，优选的，对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差包含：通过相邻道互相关法对所述瞬时相位数据体进行倾角校正处理；针对处理后的所述瞬时相位数据体中每一道数据在预设方向上所在道及前后相邻两道上分别开取预设长度的时窗，并将中心道与相邻两道时窗内相位进行作差。
在上述相位差异属性获取方法中，优选的，根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性包含：对所述窗内相位差异结果求取均值置于时窗中点，获得时窗中点沿着这一方向的横向相位差异结果；比较四个方向的横向相位差异结果，将所述横向相位差异结果中最大值作为对应道数据的横向相位差异；根据所述瞬时相位数据体中各道数据的横向相位差异获得所述叠后地震数据体的相位差异属性。
在上述相位差异属性获取方法中，优选的，将中心道与相邻两道时窗内相位进行作差包含：
通过以下公式进行作差：
ω21(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω1(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt；
ω23(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω3(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt；
在上式中，ω1、ω2、ω3为相邻三条地震道对应的瞬时相位，ω2为中心道；ω21、ω23为中心道窗内相位与相邻两道窗内相位作差的结果，Δt为半窗长，时窗长度为2Δt+1。
在上述相位差异属性获取方法中，优选的，对所述窗内相位差异结果求取均值置于时窗中点包含：
通过以下公式对所述窗内相位差异结果求取均值：
ω'(t0)＝(abs(ω21(t0+Δt,t0-Δt))+abs(ω23(t0+Δt,t0-Δt)))/(2Δt+1)；
在上式中，abs为取模运算，ω'为求得均值，ω21、ω23为长度2Δt+1的向量。
在上述相位差异属性获取方法中，优选的，对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体包含：
通过以下公式计算获得瞬时相位数据体：
s(t)＝x(t)+ih(t)＝|s(t)|eiθ(t)；
上式中，t表示时间，h(t)为实信号x(t)的Hilbert变换地震道或正交地震道，
为瞬时振幅或包络，θ(t)为瞬时相位。
本发明还提供一种相位差异属性获取系统，所述系统包含采集模块、转换模块、计算模块和分析模块；所述采集模块用于获取叠后地震数据体；所述转换模块用于对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；所述计算模块用于对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；所述分析模块用于根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性。
在上述相位差异属性获取系统中，优选的，所述计算模块还包含：通过相邻道互相关法对所述瞬时相位数据体进行倾角校正处理；针对处理后的所述瞬时相位数据体中每一道数据在预设方向上所在道及前后相邻两道上分别开取预设长度的时窗，并将中心道与相邻两道时窗内相位进行作差。
在上述相位差异属性获取系统中，优选的，所述分析模块还包含：对所述窗内相位差异结果求取均值置于时窗中点，获得时窗中点沿着这一方向的横向相位差异结果；比较四个方向的横向相位差异结果，将所述横向相位差异结果中最大值作为对应道数据的横向相位差异；根据所述瞬时相位数据体中各道数据的横向相位差异获得所述叠后地震数据体的相位差异属性。
本发明还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明提出的考虑地质体构造走向和倾向的高效相位差异属性求取方法创新地以多方向的瞬时相位差异作为局部构造判别依据，与相干等传统构造识别方法相比，其自动化计算效率极高，算法更简单，同时可以达到相同的计算效果，为地震资料的高效解释提供了新的思路，为高分辨地震反演提供了地下地质结构约束信息。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1为本发明一实施例所提供的相位差异属性获取方法的流程示意图；
图2为本发明一实施例所提供的作差流程示意图；
图3为本发明一实施例所提供的相位差异属性计算流程示意图；
图4为本发明一实施例所提供的相位差异属性获取方法的逻辑原理示意图；
图5为本发明一实施例所提供的相位差异属性获取系统的结构示意图；
图6A至图6G为本发明一实施例所提供的数值模型1的验证结果示意图；
图7A至图7E为本发明一实施例所提供的数值模型2的验证结果示意图；
图8A至图8D为本发明一实施例所提供的实际数据中河道所在层位的验证结果示意图；
图9为本发明一实施例所提供的电子设备的结构示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。
另外，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
请参考图1所示，本发明所提供的相位差异属性获取方法具体包含：S101获取叠后地震数据体；S102对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；S103对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；S104根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性；其中，所述预设方向包含正北方向、北东南西向、东西方向、北西南东方向四个方向。在实际工作中，本发明所提供的上述实施例主要可采用以下流程实现：第一步在三维地震数据体对每个地震道沿四个方位局部开窗，将窗内相位与相邻道窗内的相位对应作差取模，求得窗内的局部横向相位差异。第二步，取各方向相位差异的最大值，即利用地质构造导向思想有效挖掘不同走向的构造。该发明在相位域实现构造识别属性求取，较传统方法效率突出，而且相位较振幅更加敏感，有助于确定弱反射构造和细微构造，可为高分辨地震反演提供高精度的地下地质结构约束信息；各步骤的具体实现方式将在后续一一举例，在此就不再详细说明。
在上述实施例中，对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体可包含：
s(t)＝x(t)+ih(t)＝|s(t)|eiθ(t)；
上式中，t表示时间，h(t)为实信号x(t)的Hilbert变换地震道或正交地震道，
为瞬时振幅或包络，θ(t)为瞬时相位。
基于上述求取相邻道相位差异的实现原理，本发明通过结合构造导向和倾角校正的思想对其进一步予以改进，以实现不同走向的构造识别，同时消除倾角的影响。对地震数据体沿着多方向进行相位差异属性求取，并将各个方向进行综合得到最终的相位差异属性。具体的，请参考图2所示，在本发明一实施例中上述步骤S103中对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差可包含：S201通过相邻道互相关法对所述瞬时相位数据体进行倾角校正处理；S202针对处理后的所述瞬时相位数据体中每一道数据在预设方向上所在道及前后相邻两道上分别开取预设长度的时窗，并将中心道与相邻两道时窗内相位进行作差。再请参考图3所示，在本发明一实施例中上述步骤S104中根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性可包含：S301对所述窗内相位差异结果求取均值置于时窗中点，获得时窗中点沿着这一方向的横向相位差异结果；S302比较四个方向的横向相位差异结果，将所述横向相位差异结果中最大值作为对应道数据的横向相位差异；S303根据所述瞬时相位数据体中各道数据的横向相位差异获得所述叠后地震数据体的相位差异属性。
在上述实施例中，将中心道与相邻两道时窗内相位进行作差主要是通过以下公式进行作差：
ω21(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω1(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt；
ω23(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω3(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt；
在上式中，ω1、ω2、ω3为相邻三条地震道对应的瞬时相位，ω2为中心道；ω21、ω23为中心道窗内相位与相邻两道窗内相位作差的结果，Δt为半窗长，时窗长度为2Δt+1。
在上述实施例中，对所述窗内相位差异结果求取均值置于时窗中点，主要是通过以下公式对所述窗内相位差异结果求取均值：
ω'(t0)＝(abs(ω21(t0+Δt,t0-Δt))+abs(ω23(t0+Δt,t0-Δt)))/(2Δt+1)；
在上式中，abs为取模运算，ω'为求得均值，ω21、ω23为长度2Δt+1的向量。
整体上，本发明所提供的上述实施例主要流程可参考图4所示，主要包含首先对地震数据体进行希尔伯特变换，求得瞬时相位体，再分别沿着正北方向、北东南西向、东西方向、北西南东四个方向进行以下操作：将每一道数据局部开窗，与沿上述方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，然后将窗内相位差异结果求取均值置于时窗中点，作为该点沿着这一方向的横向相位差异属性结果。最后，取四个方向对应的横向相位差异的最大值作为该点最终的横向相位差异，即体现了地质构造导向的思想：垂直构造走向进行计算方可有效体现构造，选取各方向相位最大值可以有效识别走向为任意方向的构造。此方法具有极高的自动化计算效率和相对稳定的计算效果，应用前景良好。
为更清楚的解释上述各实施例的具体应用方式，以下将以实际计算流程为例，对上述各实施例做详细说明，本领域相关技术人员当可知，该实例仅为便于理解本发明所提供的上述实施例的具体应用方式，并不对其做进一步限定。
首先，瞬时相位主要是根据式(1)的复地震道来进行计算：
s(t)＝x(t)+ih(t)＝|s(t)|eiθ(t) (1)
其中，t表示时间，h(t)为实信号x(t)的Hilbert变换地震道或正交地震道，
为瞬时振幅或包络，θ(t)即为瞬时相位。
由上式可知，瞬时相位实际上是复地震道的相角，它可以利用反三角函数求得：
ω(t)＝arctan(imag(s)/real(s)) (2)
显然，相位属性独立于振幅属性，它不受反射系数的影响，强、弱反射的相位形态应是一致的；因此，本发明从相位属性入手进行微弱反射构造的精细识别与解释工作。在按照上述方法求得瞬时相位数据体后，即可按照下面的步骤求取横向局部相位差异：
1、首先，针对每一点，在该点所在道及前后相邻两道上分别开取长度为2Δt+1的时窗，将中心道与相邻两道窗内相位对应作差：
ω21(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω1(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt (3)
ω23(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω3(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt (4)
其中，ω1、ω2、ω3为相邻三条地震道对应的瞬时相位，ω2为中心道；ω21、ω23为中心道窗内相位与相邻两道窗内相位作差的结果，Δt为半窗长。
2、然后，将求得的ω21、ω23取模值相加后，再取平均值后置于中心道时窗中点：
ω'(t0)＝(abs(ω21(t0+Δt,t0-Δt))+abs(ω23(t0+Δt,t0-Δt)))/(2Δt+1) (5)
其中，abs表示取模运算，最后求得的ω'是一个单一值，而ω21、ω23按前面所述，是长度为2Δt+1的向量。
3、最后，将每一点按上述方法进行计算，最终即可得到与原始数据体大小相同的相位差异属性体。
该方法直接采用相邻道的局部相位差异作为构造识别的依据，在相对连续的区域，由于相位差异小，这一属性应趋于0值；而在发育构造的区域，局部相位变化剧烈，这一属性应有相对较大的非0值。显然，在整个求取过程中，仅对相位属性进行了简单的加、减、取模和取平均运算，因此，这一方法较常规相干方法有着突出的计算效率。但是，由于相位属性对噪声较为敏感，这一方法要求数据具有较高的信噪比。
进一步的，由于地下构造往往具有不同的走向，例如河道和延伸较广的断层等，在不同部位往往发育方向各异。上述方法沿着某一方向(如Crossline方向)求取相邻道的相位差异，则当构造走向与这一求取方向垂直时有良好的求取效果，但是当二者相互平行时，由于求取的方向与边界方向一致，相位的差异则不够明显。也即上述方法只能体现走向沿着某一确定方向的构造，无法体现走向平行于求取方向的构造。
为了解决这一问题，本发明提供在求取过程中引入“构造导向”的思想；这种思想可以概括为，沿着多方向进行求取，然后将各方向结果进行综合得到最终结果。具体做法如下：
1、首先，我们将沿着单一方向进行相位差异求取改为沿着正北方向(0°)、北东南西向(45°)、东西方向(90°)、北西南东(135°)四个方向进行。
2、本发明将沿着四个方向求得的相位结果分别记为ω0、ω45、ω90、ω135。这是四个与原始数据等大的三维数据体。然后，取每一点四个方向的最大值作为最终的相位差异结果：
ω(t0)＝max(ω0(t0),ω45(t0),ω90(t0),ω135(t0)) (5)
这样即可保证在最终的结果上体现不同走向的构造(构造在对应于走向方向的相位结果上表现一个较大值)。
3、此外，倾角的影响也会造成局部相位发生变化，从而造成假象，干扰对构造的精细表征；因此，本发明利用互相关法对相邻地震道进行倾角扫描，该方法的原理如下：将某一道相对相邻道时移一系列距离，再对每一个时移量后的两道做互相关运算，最终找到相关值最大的时移量，这一时移即可视为局部倾角。
4、然后，基于前面的计算过程中将相邻窗也时移与局部倾角相等的量，即实现了倾角校正。
请参考图5所示，本发明还提供一种相位差异属性获取系统，所述系统包含采集模块、转换模块、计算模块和分析模块；所述采集模块用于获取叠后地震数据体；所述转换模块用于对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；所述计算模块用于对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；所述分析模块用于根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性。
在上述实施例中，所述计算模块还可包含：通过相邻道互相关法对所述瞬时相位数据体进行倾角校正处理；针对处理后的所述瞬时相位数据体中每一道数据在预设方向上所在道及前后相邻两道上分别开取预设长度的时窗，并将中心道与相邻两道时窗内相位进行作差；所述分析模块还可包含：对所述窗内相位差异结果求取均值置于时窗中点，获得时窗中点沿着这一方向的横向相位差异结果；比较四个方向的横向相位差异结果，将所述横向相位差异结果中最大值作为对应道数据的横向相位差异；根据所述瞬时相位数据体中各道数据的横向相位差异获得所述叠后地震数据体的相位差异属性。
本发明还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明提出的考虑地质体构造走向和倾向的高效相位差异属性求取方法创新地以多方向的瞬时相位差异作为局部构造判别依据，与相干等传统构造识别方法相比，其自动化计算效率极高，算法更简单，同时可以达到相同的计算效果，为地震资料的高效解释提供了新的思路，为高分辨地震反演提供了地下地质结构约束信息。
为验证本发明所提供的相位差异属性获取方法及系统，以下分别使用了两个数值模型、一个实际数据对相位差异属性获取方法及系统进行测试。
本发明设计的数值模型1以波速为2000m/s的围岩和波速为3000m/s的目标砂岩层组成，纵向总时间深度为151ms，Inline、Crossline方向采样点均为201。砂体厚12ms，位于70ms～82ms范围内，主要包括一个河道和两个冲积扇状构造。图6A为该数值模型在无噪条件下的水平切片，切片位置为砂体层位的中心，其中白色部分为围岩，黑色部分即为砂体。图6B为该数值模型沿着Crossline100方向的剖面。图6C至图6F分别为本发明求得的单方向相位差异属性，方向分别为0°、45°、90°和135°。可以看到，单一方向的结果只能体现某一方向的构造边界。以0°方向为例，走向沿着东西方向的部分可以清楚地体现，然而与这一方向垂直的部分则不甚清楚。图6G为按照上述方法得到的综合四个方向得到的最终相位差异属性，可以看到构造边界的所有细节都得到了清楚的展示。
数值模型2是1995年Aminzadeh等在SEG/EAGE年会上提出的，由30Hz的雷克子波和反射系数褶积而成。该数值模型大小为800km2×186ms，时间采样间隔为1ms，空间采样间隔为50m×50m。它是一个包含了河道和断层的逆冲推覆体模型，该模型模拟复杂的地质结构，顶界面是一个不整合面，不整合发育的12个连续层的沉积层序。它的三维立体图如图7A所示。图7B为该数值模型的水平切片，切片位置为100ms，大致与层位相吻合。图7C即为常规相干方法得到的结果；图7D即为与图7B对应的没有进行倾角校正的相位差异结果，可以看到由于倾角影响，有一定的串层现象(注意：为了方便与相干结果的比较，本发明将此后的相位差异结果取了相反数，因此在图中均为负值，较小值对应构造)；图7E为与图7B对应的带有倾角校正的相位差异结果，可以看到倾角校正有效消除了图7D中的背景假象，计算效果与相干结果接近。而这一方法的计算时间大约仅为常规相干的三分之一，由此可见此发明的高效性。
选用的实际数据包括一个清楚的河道构造。图8A为该数据中河道所在层位的沿层振幅切片，可以比较清楚地看到河道的大致位置和基本形态；图8B即为常规相干方法得到的计算结果沿层切片；图8C即为没有进行倾角校正的相位差异结果的沿层切片；图8D为带有倾角校正的相位差异结果的沿层切片。通过比较图8C、图8D可以看出倾角消除的效果，同时对于这一实际数据，本发明的计算速度比常规相干快约50％。
如图9所示，该电子设备600还可以包括：通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是，电子设备600也并不是必须要包括图9中所示的所有部件；此外，电子设备600还可以包括图9中没有示出的部件，可以参考现有技术。
如图9所示，中央处理器100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器100接收输入并控制电子设备600的各个部件的操作。
其中，存储器140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序，以实现信息存储或处理等。
输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142，该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。
存储器140还可以包括数据存储部143，该数据存储部143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132，以经由扬声器131提供音频输出，并接收来自麦克风132的音频输入，从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器130还耦合到中央处理器100，从而使得可以通过麦克风132能够在本机上录音，且使得可以通过扬声器131来播放本机上存储的声音。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种相位差异属性获取方法，其特征在于，所述方法包含：
获取叠后地震数据体；
对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；
对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；
根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性。
2.根据权利要求1所述的相位差异属性获取方法，其特征在于，所述预设方向包含正北方向、北东南西向、东西方向、北西南东方向四个方向。
3.根据权利要求2所述的相位差异属性获取方法，其特征在于，对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差包含：
通过相邻道互相关法对所述瞬时相位数据体进行倾角校正处理；
针对处理后的所述瞬时相位数据体中每一道数据在预设方向上所在道及前后相邻两道上分别开取预设长度的时窗，并将中心道与相邻两道时窗内相位进行作差。
4.根据权利要求3所述的相位差异属性获取方法，其特征在于，根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性包含：
对所述窗内相位差异结果求取均值置于时窗中点，获得时窗中点沿着这一方向的横向相位差异结果；
比较四个方向的横向相位差异结果，将所述横向相位差异结果中最大值作为对应道数据的横向相位差异；
根据所述瞬时相位数据体中各道数据的横向相位差异获得所述叠后地震数据体的相位差异属性。
5.根据权利要求4所述的相位差异属性获取方法，其特征在于，将中心道与相邻两道时窗内相位进行作差包含：
通过以下公式进行作差：
ω21(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω1(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt；
ω23(t0+Δt,t0-Δt)＝ω2(t0+Δt,t0-Δt)-ω3(t0+Δt,t0-Δt),t0＝Δt+1～t-Δt；
在上式中，ω1、ω2、ω3为相邻三条地震道对应的瞬时相位，ω2为中心道；ω21、ω23为中心道窗内相位与相邻两道窗内相位作差的结果，Δt为半窗长，时窗长度为2Δt+1。
6.根据权利要求5所述的相位差异属性获取方法，其特征在于，对所述窗内相位差异结果求取均值置于时窗中点包含：
通过以下公式对所述窗内相位差异结果求取均值：
ω'(t0)＝(abs(ω21(t0+Δt,t0-Δt))+abs(ω23(t0+Δt,t0-Δt)))/(2Δt+1)；
在上式中，abs为取模运算，ω'为求得均值，ω21、ω23为长度2Δt+1的向量。
7.根据权利要求1所述的相位差异属性获取方法，其特征在于，对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体包含：
通过以下公式计算获得瞬时相位数据体：
s(t)＝x(t)+ih(t)＝|s(t)|eiθ(t)；
上式中，t表示时间，h(t)为实信号x(t)的Hilbert变换地震道或正交地震道，
为瞬时振幅或包络，θ(t)为瞬时相位。
8.一种相位差异属性获取系统，其特征在于，所述系统包含采集模块、转换模块、计算模块和分析模块；
所述采集模块用于获取叠后地震数据体；
所述转换模块用于对所述叠后地震数据体进行希尔伯特变换获得瞬时相位数据体；
所述计算模块用于对所述瞬时相位数据体中每一道数据局部开窗，并沿预设方向的相邻两道在消除倾角影响后的等宽时窗内相位对应作差，获得所述瞬时相位数据体中各道数据的窗内相位差异结果；
所述分析模块用于根据所述瞬时相位数据体中各道数据的窗内相位差异结果获得所述叠后地震数据体的相位差异属性。
9.根据权利要求8所述的相位差异属性获取系统，其特征在于，所述计算模块还包含：
通过相邻道互相关法对所述瞬时相位数据体进行倾角校正处理；
针对处理后的所述瞬时相位数据体中每一道数据在预设方向上所在道及前后相邻两道上分别开取预设长度的时窗，并将中心道与相邻两道时窗内相位进行作差。
10.根据权利要求9所述的相位差异属性获取系统，其特征在于，所述分析模块还包含：
对所述窗内相位差异结果求取均值置于时窗中点，获得时窗中点沿着这一方向的横向相位差异结果；
比较四个方向的横向相位差异结果，将所述横向相位差异结果中最大值作为对应道数据的横向相位差异；
根据所述瞬时相位数据体中各道数据的横向相位差异获得所述叠后地震数据体的相位差异属性。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7任一所述方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至7任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010084606.3A CN113253345B (zh) | 2020-02-10 | 2020-02-10 | 相位差异属性获取方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010084606.3A CN113253345B (zh) | 2020-02-10 | 2020-02-10 | 相位差异属性获取方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113253345A true CN113253345A (zh) | 2021-08-13 |
| CN113253345B CN113253345B (zh) | 2024-06-25 |
Family
ID=77219375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010084606.3A Active CN113253345B (zh) | 2020-02-10 | 2020-02-10 | 相位差异属性获取方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113253345B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104597485A (zh) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | 一种微小断层检测方法及断层检测装置 |
| CN107144877A (zh) * | 2017-04-12 | 2017-09-08 | 中国石油天然气股份有限公司 | 异常构造识别方法及装置 |
| CN108398721A (zh) * | 2018-04-13 | 2018-08-14 | 中国石油大学(北京) | 基于地质导向的相移属性的识别方法及装置 |
| CN110095813A (zh) * | 2019-05-17 | 2019-08-06 | 中国石油大学(北京) | 用于确定复地震道的瞬时相位差的方法及装置 |
| CN110261905A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于倾角控制的复值相干微断层识别方法 |
-
2020
- 2020-02-10 CN CN202010084606.3A patent/CN113253345B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104597485A (zh) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | 一种微小断层检测方法及断层检测装置 |
| CN107144877A (zh) * | 2017-04-12 | 2017-09-08 | 中国石油天然气股份有限公司 | 异常构造识别方法及装置 |
| CN108398721A (zh) * | 2018-04-13 | 2018-08-14 | 中国石油大学(北京) | 基于地质导向的相移属性的识别方法及装置 |
| CN110095813A (zh) * | 2019-05-17 | 2019-08-06 | 中国石油大学(北京) | 用于确定复地震道的瞬时相位差的方法及装置 |
| CN110261905A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于倾角控制的复值相干微断层识别方法 |
Non-Patent Citations (1)
| Title |
|---|
| SANYI YUAN ET AL: "Geosteering Phase Attributes: A New Detector for the Discontinuities of Seismic Images", IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, vol. 16, no. 1, XP011695174, DOI: 10.1109/LGRS.2018.2866419 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113253345B (zh) | 2024-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102893182B (zh) | 用于建模3d地质结构的系统和方法 | |
| AU2012220584B2 (en) | Sensitivity kernel-based migration velocity analysis in 3D anisotropic media | |
| CN108398721B (zh) | 基于地质导向的相移属性的识别方法及装置 | |
| EP2160633A1 (en) | Creating an absorption parameter model | |
| CN109946741B (zh) | 一种TTI介质中纯qP波最小二乘逆时偏移成像方法 | |
| CN112415591A (zh) | 绕射波的成像方法、装置、电子设备及存储介质 | |
| Sherwood et al. | 3D beam prestack depth migration with examples from around the world | |
| CN110703331A (zh) | 一种基于常q粘滞声波方程的衰减补偿逆时偏移实现方法 | |
| CN106125139B (zh) | 一种三维地震数据处理方法及系统 | |
| An et al. | Characteristics of low‐frequency horizontal noise of ocean‐bottom seismic data | |
| Xiao et al. | An offshore Gabon full-waveform inversion case study | |
| CN113253345B (zh) | 相位差异属性获取方法及系统 | |
| CN107765306B (zh) | 一种vsp初始速度建模方法及装置 | |
| CN115480296A (zh) | 天然气水合物稳定带底界预测方法及装置 | |
| Shan et al. | Imaging overturned waves by plane-wave migration in tilted coordinates | |
| Waheed | A fast-marching eikonal solver for tilted transversely isotropic media | |
| Shan et al. | Plane-wave migration in tilted coordinates | |
| CN110261896B (zh) | 稳定的各向异性ti介质正演模拟方法 | |
| CN113933901B (zh) | 深度域地震数据拼接方法及装置 | |
| Li et al. | Seismic full-waveform inversion based on superresolution for high-precision prediction of reservoir parameters | |
| CN112711064A (zh) | 初始层速度建模方法及装置 | |
| CN114114417B (zh) | 地震成像方法及装置 | |
| US12050294B1 (en) | Method and apparatus for building envelope guided low frequency model for seismic exploration of subsurface formations | |
| Zhao et al. | Velocity building strategy for imaging deep low relief structure in Tarim basin | |
| CN107193047A (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 |