CN109307890A - 基于上下行波场分解的逆时偏移方法及系统 - Google Patents
基于上下行波场分解的逆时偏移方法及系统 Download PDFInfo
- Publication number
- CN109307890A CN109307890A CN201710633034.8A CN201710633034A CN109307890A CN 109307890 A CN109307890 A CN 109307890A CN 201710633034 A CN201710633034 A CN 201710633034A CN 109307890 A CN109307890 A CN 109307890A
- Authority
- CN
- China
- Prior art keywords
- wave field
- uplink
- wave
- reverse
- downlink
- 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.)
- Pending
Links
- 230000005012 migration Effects 0.000 title claims abstract description 58
- 238000013508 migration Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 31
- 238000003384 imaging method Methods 0.000 claims abstract description 47
- 230000009466 transformation Effects 0.000 claims abstract description 32
- 230000001902 propagating effect Effects 0.000 claims abstract description 11
- 230000001419 dependent effect Effects 0.000 claims description 7
- 238000013213 extrapolation Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 5
- 230000006854 communication Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 5
- 238000001914 filtration Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- NIICMVKXWJNYSB-UHFFFAOYSA-N 13-Hydroxy-Leontalbinine Natural products OC1CC(=O)N2CC3CCCN4CCCC(=C2C1)C34 NIICMVKXWJNYSB-UHFFFAOYSA-N 0.000 description 1
- CDDHEMJXKBELBO-UHFFFAOYSA-N Leontalbinin Natural products C1CCC2CN3C(=O)CCCC3=C3C2N1CCC3 CDDHEMJXKBELBO-UHFFFAOYSA-N 0.000 description 1
- 240000007711 Peperomia pellucida Species 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production 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
- 238000012805 post-processing Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- 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
- G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
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
本发明提出了一种基于上下行波场分解的逆时偏移方法及系统，该方法包括：构建震源与地震记录的解析波场；对解析波场进行波场的正向与反向延拓；在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；由原始波场与所述变换波场求解上下行波场；选择特定方向的波场利用互相关成像条件进行成像。基于上下行波分解的逆时偏移方法能有效压制逆时偏移中的低频噪音，提高成像剖面精度；同时具有较好的子波保幅性，这也为后续保真处理提供了技术保证。
Description
技术领域
本发明属于地震资料叠前深度偏移成像技术领域，具体地涉及基于上下行波场分解的逆时偏移方法及系统。
背景技术
随着复杂地区的勘探开发，获得正确的构造认识和油气有利圈闭对叠前偏移成像的精度要求越来越高。目前叠前深度偏移算法主要分为基于射线理论和波动方程两类，Kirchhoff偏移和单程波偏移就是两类方法的代表，具有计算稳定、计算效率高等优点，但是在面对复杂的地质构造时存在着算法上的不足，如多路径、倾角限制等问题。
逆时偏移通过求解双程波动方程，实现波场的正反向传播，能够精确地描述各类地震波的传播，且不受倾角的限制，是目前最精确的地震成像方法。由于逆时偏移在应用互相关成像条件进行成像时没有考虑波场的方向性，会造成成像结果中产生低波数高振幅的偏移噪音和偏移假象，降低了偏移剖面的成像精度。
实际工业生产中，主要解决方法是在叠后进行拉普拉斯滤波，但拉普拉斯滤波只能压制低波数、高振幅的偏移噪音，不能有效消除成像假象，同时会改变成像的振幅以及频率相位信息，不具有保幅性，不利于后续地震解释工作。
发明内容
复杂构造工区进行叠前逆时偏移成像时，针对全波场互相关成像条件中产生的低波数、高振幅的偏移噪音以及相干偏移假象问题，本发明介绍了一种基于炮检点上下行波场分解的逆时偏移算法，通过构建新的成像条件，利用特定方向波场进行互相关成像，在保证成像精度的同时能有效解决低频噪音和偏移假象问题，提高成像质量，对于推动高精度地震成像具有重要意义。
根据本发明的一个方面，提供一种基于上下行波场分解的逆时偏移方法，该方法包括：
构建震源与地震记录的解析波场；
对解析波场进行波场的正向与反向延拓；
在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；
由原始波场与所述变换波场求解上下行波场；
选择特定方向的波场利用互相关成像条件进行成像。
进一步地，由震源与地震记录原始波场及其希尔伯特变换结果，构建解析波场：
式中，其中为解析波场，u为原始波场，uh为u的希尔伯特变换结果。
进一步地，基于以下波动方程组实现对解析波场的正向与反向波场延拓：
式中，为震源或地震记录，为震源或地震记录的希尔伯特变换，为传播过程中原始波场及其希尔伯特变换波场，为地震波速度。
进一步地，对解析波场进行波场的正向与反向延拓，得到每个传播时刻的原始波场与沿时间方向做完希尔伯特变换后的变换波场。
进一步地，在波场传播的每个时刻，沿成像网格深度方向对uh进行z方向的希尔伯特变换，得到每个传播时刻做完时间与深度希尔伯特变换后的变换波场：
式中，uh为u的希尔伯特变换结果，HtHz(u(x,t))为时间方向与深度方向的希尔伯特变换波场。
进一步地，由原始波场与变换波场求解上下行波场，其表达式为：
udown(x,t)＝[u(x,t)+HtHz(u(x,t))]/2
uup(x,t)＝[u(x,t)-HtHz(u(x,t))]/2
式中，Ht与Hz分别表示在时间方向与深度方向上的希尔伯特变换，u(x,t)为原始传播波场，udown(x,t)、uup(x,t)为下行波场与上行波场。
进一步地，基于上下行波场与原始波场以及变换波场间的关系式，完成炮点与检波点的上下行波场分离。
进一步地，利用传播方向相反波场进行互相关成像，成像条件为：
式中，Su(t,m)、Sd(t,m)为震源上行与下行波场；Ru(t,m)、Rd(t,m)为检波端上行与下行波场。
根据本发明的另一方面，提供一种基于上下行波场分解的逆时偏移系统，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
构建震源与地震记录的解析波场；
对解析波场进行波场的正向与反向延拓；
在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；
由原始波场与所述变换波场求解上下行波场；
选择特定方向的波场利用互相关成像条件进行成像。
进一步地，对炮点波场与检波点波场进行上下行波分离后，只对传播方向相反的波场进行相关成像。
理论模型的波场分解逆时偏移结果表明：基于上下行波分解的逆时偏移方法能有效压制逆时偏移中的低频噪音，提高成像剖面精度；同时具有较好的子波保幅性，这也为后续保真处理提供了技术保证。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了本发明实施例的处理流程图。
图2示出了单道脉冲常规逆时偏移结果。
图3示出了单道脉冲运用本发明得到的逆时偏移结果。
图4示出了炮点波场与上下行分解波场。
图5示出了Marmousi模型常规逆时偏移结果。
图6示出了本发明得到的逆时偏移结果。
图7示出了拉普拉斯滤波后结果。
图8示出了单道对比，其中圆实线为常规逆时偏移结果，三角实线为本发明结果，实线为拉普拉斯滤波后结果。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明属于地震资料叠前深度偏移成像技术领域。由于常规逆时偏移在应用全波场互相关成像条件进行成像时没有考虑波场的方向性，会造成成像结果中产生低波数、高振幅的偏移噪音和偏移假象，降低了偏移剖面的成像精度。本发明介绍了一种基于炮、检点上下行波场分解的逆时偏移算法，通过构建新的成像条件，利用特定方向波场进行互相关成像，在保证成像精度的同时可以有效解决低频噪音和偏移假象问题，提高了成像剖面精度。
在频率-波数域，利用频率和波数乘积的符号可以将上下行波场进行很好的分解，但由于时间维度为数据的最慢维，需要将炮点和检波点所有时间片波场进行存盘后再进行傅里叶变换，因此频率-波数域的波场分解对数据的I/O和硬盘的存储要求巨大，现阶段还不能满足实际生产的需要。基于解析时间波场的外推，本文提出了一种在时空域将波场显示地进行上、下行波场分解的有效方法，并对特定方向波场进行相关成像，能有效压制偏移噪音与偏移假象。
传统波场的分解可以在f-k域通过傅里叶变换进行实现，其上下行波场由频率与波场的符号共同决定，但是该方法计算效率效率低且存储大；基于希尔伯特变换(HT)与傅里叶变换(FT)以及两者间的关系，Shen and Albertin(2015，SEG)推导了新的上下行波分解公式：
udown(x,t)＝[u(x,t)+HtHz(u(x,t))]/2
uup(x,t)＝[u(x,t)-HtHz(u(x,t))]/2 (1)
其中Ht与Hz分别表示在时间方向与深度方向上的希尔伯特变换。
在频率域HT与FT存在着以下关系：
HT(x(t))＝-jgsgnFT(x(t)) (2)
其中符号函数为：
基于公式1的波场分解公式，可以利用FT与HT间的关系式2实现上下行波的分解，但与传统方法类似，计算效率也是制约该方法在实际生产中应用的瓶颈之一。
为了解决波场分解的计算效率与存储问题，必须在时空域进行研究。本方法提出类比解析信号的定义方式，得到解析波场定义为：
其中为解析波场；u为原始波场；uh为u的希尔伯特变换结果。
在时间域希尔伯特变换可以进行高效求解，其表达式为：
解析波场的外推意味着在进行逆时偏移波场延拓时，需要同时进行两套波场的波场延拓。炮点与检波点波场传播算子类似，描述波在声波介质中传播的波动方程为：
为震源或地震记录的希尔伯特变换，通过有限差分求解上述偏微分方程组，即可以得到解析波场的实部与虚部波场。
由于波场传播算子具有线性特征，因此，对震源做HT后进行波场的传播得到的波场uh相当于对波场进行了时间方向的HT，由于不需要储存每个时刻的波场，节省了大量的存储空间。在每一个传播时刻得到的波场uh是关于空间位置的函数，然后沿深度方向做HT即沿成像网格z方向做HT，此步可利用公式(5)进行高效实现。基于上述波场传播，在时空域得到了波场时间方向与深度方向的希尔伯特变换波场即HtHz(u(x,t))，最后通过求解公式(1)可以实现炮点与检波点的上、下行波场分离。
逆时偏移中偏移噪音是由于炮点与检波点方向相同的波场进行相关成像产生的，利用本发明对炮点波场与检波点波场进行上下行波分离后，只对传播方向相反的波场进行相关成像，在保证成像精度的同时有效压制了偏移噪音，与叠后拉普拉斯滤波方法相比，本方法不改变子波的相位与振幅，具有更好的保幅性。
根据本发明的一个方面，提供一种基于上下行波场分解的逆时偏移方法，该方法包括：
构建震源与地震记录的解析波场；
对解析波场进行波场的正向与反向延拓；
在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；
由原始波场与所述变换波场求解上下行波场；
选择特定方向的波场利用互相关成像条件进行成像。
具体而言，通过构建震源与地震记录的解析波场，并对解析波场进行波场的正向与反向延拓，得到每个传播时刻的原始波场与沿时间方向做完希尔伯特变换后的变换波场。在成像网格沿深度方向对每个时刻的波场进行希尔伯特变换，得到每个传播时刻最终做完时间与深度希尔伯特变换后的变换波场。基于上下行波场与原始波场以及变换波场间的关系式(1)，完成炮点与检波点的上下行波场分离。最后选择特定方向的波场利用互相关成像条件进行成像，可以有效压制逆时偏移噪音，提高成像精度。
图1为本发明实施例的实现步骤流程示意图，具体实现步骤包括：
步骤一：由震源与地震记录原始波场及其HT结果，构建解析波场：
式中，其中为解析波场，u为原始波场，uh为u的希尔伯特变换结果。
步骤二：基于以下波动方程实现对解析波场的正向与反向波场延拓：
式中，为震源或地震记录，为震源或地震记录的希尔伯特变换，为传播过程中原始波场及其希尔伯特变换波场，为地震波速度。
对解析波场进行波场的正向与反向延拓，得到每个传播时刻的原始波场与沿时间方向做完希尔伯特变换后的变换波场。
步骤三：在波场传播的每个时刻，沿成像网格深度方向对uh进行z方向的希尔伯特变换，得到每个传播时刻做完时间与深度希尔伯特变换后的变换波场：
式中，uh为u的希尔伯特变换结果，HtHz(u(x,t))为时间方向与深度方向的希尔伯特变换波场。
步骤四：由原始波场与变换波场求解上下行波场，其表达式为：
udown(x,t)＝[u(x,t)+HtHz(u(x,t))]/2
uup(x,t)＝[u(x,t)-HtHz(u(x,t))]/2
式中，Ht与Hz分别表示在时间方向与深度方向上的希尔伯特变换，u(x,t)为原始传播波场，udown(x,t)、uup(x,t)为下行波场与上行波场。
基于上下行波场与原始波场以及变换波场间的关系式，完成炮点与检波点的上下行波场分离。
步骤五：利用传播方向相反波场进行互相关成像，成像条件为：
式中，Su(t,m)、Sd(t,m)为震源上行与下行波场；Ru(t,m)、Rd(t,m)为检波端上行与下行波场。
根据本发明的另一方面，提供一种基于上下行波场分解的逆时偏移系统，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
构建震源与地震记录的解析波场；
对解析波场进行波场的正向与反向延拓；
在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；
由原始波场与所述变换波场求解上下行波场；
选择特定方向的波场利用互相关成像条件进行成像。
进一步地，对炮点波场与检波点波场进行上下行波分离后，只对传播方向相反的波场进行相关成像。
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
通过单道脉冲信号和Marmousi理论模型，运用上下行波场分解的逆时偏移处理，验证了本方法对常规互相关逆时偏移中存在的偏移噪音和假象具有很好的压制效果，提高了逆时偏移的成像精度，相对于叠后拉普拉斯滤波处理具有更好的子波保幅性。
图2为脉冲常规互相关逆时偏移结果，可以看到由于炮点与检波点同时包含上行波与下行波，利用互相关成像条件进行成像时，炮检点传播方向相反的波能有效成像，而传播方向相同的波场则形成高振幅的低频噪音，在脉冲的逆时偏移结果上表现为传播路劲上的“兔耳朵”，干扰了成像结果。
图3为本发明基于上下行波场分解后的逆时偏移结果，其中波场路径上的低频噪音得到了很好的压制；图4为脉冲逆时偏移过程中某一时刻的原始炮点波场以及运用本发明得到的上下行波场，其上下行波场被很好的分解，验证了本方法的正确性与有效性。
图5为Marmousi模型常规互相关逆时偏移结果，偏移结果特别是浅层受低频干扰严重，这与逆时偏移中低频噪音产生机制相吻合。图6为运用本发明后的逆时偏移结果，低波数、高振幅的偏移噪音得到压制，提高了成像剖面分辨率。图7为图5得到的常规逆时偏移结果进行叠后拉普拉斯滤波后的结果，由于拉普拉斯算子中包含振幅和相位信息，因此滤波后对成像剖面的振幅和相位都会产生影响。
图8为图5-7中某一个cdp的单道对比结果，其中圆实线为常规逆时偏移结果，三角实线为本发明得到的偏移结果，实线为拉普拉斯滤波后处理结果，对比可知，本方法得到的偏移结果与常规互相关成像结果具有很好的振幅一致性，而叠后拉普拉斯处理结果振幅变化较大。
理论模型的波场分解逆时偏移结果表明：基于上下行波分解的逆时偏移方法能有效压制逆时偏移中的低频噪音，提高成像剖面精度；同时具有较好的子波保幅性，这也为后续保真处理提供了技术保证。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种基于上下行波场分解的逆时偏移方法，其特征在于，该方法包括：
构建震源与地震记录的解析波场；
对解析波场进行波场的正向与反向延拓；
在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；
由原始波场与所述变换波场求解上下行波场；
选择特定方向的波场利用互相关成像条件进行成像。
2.根据权利要求1所述的基于上下行波场分解的逆时偏移方法，其特征在于，由震源与地震记录原始波场及其希尔伯特变换结果，构建解析波场：
式中，其中为解析波场，u为原始波场，uh为u的希尔伯特变换结果。
3.根据权利要求1所述的基于上下行波场分解的逆时偏移方法，其特征在于，基于以下波动方程组实现对解析波场的正向与反向波场延拓：
式中，为震源或地震记录，为震源或地震记录的希尔伯特变换，为传播过程中原始波场及其希尔伯特变换波场，为地震波速度。
4.根据权利要求3所述的基于上下行波场分解的逆时偏移方法，其特征在于，对解析波场进行波场的正向与反向延拓，得到每个传播时刻的原始波场与沿时间方向做完希尔伯特变换后的变换波场。
5.根据权利要求1所述的基于上下行波场分解的逆时偏移方法，其特征在于，在波场传播的每个时刻，沿成像网格深度方向对uh进行z方向的希尔伯特变换，得到每个传播时刻做完时间与深度希尔伯特变换后的变换波场：
式中，uh为u的希尔伯特变换结果，HtHz(u(x,t))为时间方向与深度方向的希尔伯特变换波场。
6.根据权利要求1所述的基于上下行波场分解的逆时偏移方法，其特征在于，由原始波场与变换波场求解上下行波场，其表达式为：
udown(x,t)＝[u(x,t)+HtHz(u(x,t))]/2
uup(x,t)＝[u(x,t)-HtHz(u(x,t))]/2
式中，Ht与Hz分别表示在时间方向与深度方向上的希尔伯特变换，u(x,t)为原始传播波场，udown(x,t)、uup(x,t)为下行波场与上行波场。
7.根据权利要求6所述的基于上下行波场分解的逆时偏移方法，其特征在于，基于上下行波场与原始波场以及变换波场间的关系式，完成炮点与检波点的上下行波场分离。
8.根据权利要求1所述的基于上下行波场分解的逆时偏移方法，其特征在于，利用传播方向相反波场进行互相关成像，成像条件为：
式中，Su(t,m)、Sd(t,m)为震源上行与下行波场；Ru(t,m)、Rd(t,m)为检波端上行与下行波场。
9.一种基于上下行波场分解的逆时偏移系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
构建震源与地震记录的解析波场；
对解析波场进行波场的正向与反向延拓；
在成像网格沿深度方向对每个传播时刻的波场进行希尔伯特变换，得到变换后的变换波场；
由原始波场与所述变换波场求解上下行波场；
选择特定方向的波场利用互相关成像条件进行成像。
10.根据权利要求9所述的基于上下行波场分解的逆时偏移系统，其特征在于，对炮点波场与检波点波场进行上下行波分离后，只对传播方向相反的波场进行相关成像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710633034.8A CN109307890A (zh) | 2017-07-28 | 2017-07-28 | 基于上下行波场分解的逆时偏移方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710633034.8A CN109307890A (zh) | 2017-07-28 | 2017-07-28 | 基于上下行波场分解的逆时偏移方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109307890A true CN109307890A (zh) | 2019-02-05 |
Family
ID=65205312
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710633034.8A Pending CN109307890A (zh) | 2017-07-28 | 2017-07-28 | 基于上下行波场分解的逆时偏移方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109307890A (zh) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901223A (zh) * | 2019-03-20 | 2019-06-18 | 吉林大学 | 基于稳相区间理论的逆时偏移成像方法 |
| CN110018517A (zh) * | 2019-05-07 | 2019-07-16 | 西安石油大学 | 一种多尺度地面微地震逆时干涉定位方法 |
| CN110398776A (zh) * | 2019-08-21 | 2019-11-01 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
| CN112444863A (zh) * | 2019-08-28 | 2021-03-05 | 中国石油化工股份有限公司 | 粘声介质波场延拓和波场分解方法、存储介质及计算设备 |
| CN112578455A (zh) * | 2019-09-30 | 2021-03-30 | 中国石油化工股份有限公司 | 一种空间波数混合域地震波场多方向分解方法及系统 |
| WO2021141903A1 (en) * | 2020-01-09 | 2021-07-15 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging |
| CN113126150A (zh) * | 2018-12-28 | 2021-07-16 | 中国石油化工股份有限公司 | 基于单程波动方程的增强地震成像的方法和设备 |
| CN113419282A (zh) * | 2021-07-19 | 2021-09-21 | 中国矿业大学(北京) | 垂直地震剖面波场分离方法和系统 |
| US11275190B2 (en) | 2018-05-16 | 2022-03-15 | Saudi Arabian Oil Company | Generating diffraction images based on wave equations |
| US11313988B2 (en) | 2019-12-13 | 2022-04-26 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction imaging |
| US11467303B2 (en) | 2020-03-09 | 2022-10-11 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using a post-stack seismic diffraction imaging condition |
| CN115755175A (zh) * | 2022-11-18 | 2023-03-07 | 成都理工大学 | 基于Hilbert变换的复合波场弹性波逆时偏移成像方法 |
| US11681043B2 (en) | 2019-09-03 | 2023-06-20 | Saudi Arabian Oil Company | Diffraction imaging using pseudo dip-angle gather |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706583A (zh) * | 2009-10-16 | 2010-05-12 | 西安交通大学 | 一种多偏移距vsp成像的局域化相空间方法 |
| US20120051176A1 (en) * | 2010-08-31 | 2012-03-01 | Chevron U.S.A. Inc. | Reverse time migration back-scattering noise removal using decomposed wavefield directivity |
| US20170003410A1 (en) * | 2015-07-02 | 2017-01-05 | Korea Institute Of Geoscience And Mineral Resource | Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging |
-
2017
- 2017-07-28 CN CN201710633034.8A patent/CN109307890A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706583A (zh) * | 2009-10-16 | 2010-05-12 | 西安交通大学 | 一种多偏移距vsp成像的局域化相空间方法 |
| US20120051176A1 (en) * | 2010-08-31 | 2012-03-01 | Chevron U.S.A. Inc. | Reverse time migration back-scattering noise removal using decomposed wavefield directivity |
| US20170003410A1 (en) * | 2015-07-02 | 2017-01-05 | Korea Institute Of Geoscience And Mineral Resource | Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging |
Non-Patent Citations (4)
| Title |
|---|
| PENG SHEN AND UWE ALBERTIN: "Up-Down Separation Using Hilbert Transformed Source for Causal Imaging Condition", 《SEG NEW ORLEANS ANNUAL MEETING》 * |
| SIRUI TAN，等: "A staggered-grid finite-difference scheme optimized in the time-space domain for modeling scalar-wave propagation", 《SEG》 * |
| 胡江涛，等: "基于解析时间波场外推与波场分解的逆时偏移方法研究", 《地球物理学报》 * |
| 郭鹏，等: "VTI介质弹性波波场分解的空间域算法", 《石油地球物理勘探》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11275190B2 (en) | 2018-05-16 | 2022-03-15 | Saudi Arabian Oil Company | Generating diffraction images based on wave equations |
| CN113126150B (zh) * | 2018-12-28 | 2024-02-09 | 中国石油化工股份有限公司 | 基于单程波动方程的增强地震成像的方法和设备 |
| CN113126150A (zh) * | 2018-12-28 | 2021-07-16 | 中国石油化工股份有限公司 | 基于单程波动方程的增强地震成像的方法和设备 |
| CN109901223A (zh) * | 2019-03-20 | 2019-06-18 | 吉林大学 | 基于稳相区间理论的逆时偏移成像方法 |
| CN110018517B (zh) * | 2019-05-07 | 2021-04-02 | 西安石油大学 | 一种多尺度地面微地震逆时干涉定位方法 |
| CN110018517A (zh) * | 2019-05-07 | 2019-07-16 | 西安石油大学 | 一种多尺度地面微地震逆时干涉定位方法 |
| CN110398776B (zh) * | 2019-08-21 | 2021-03-30 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
| CN110398776A (zh) * | 2019-08-21 | 2019-11-01 | 西南石油大学 | 一种基于右手法则的波场分离方法 |
| CN112444863A (zh) * | 2019-08-28 | 2021-03-05 | 中国石油化工股份有限公司 | 粘声介质波场延拓和波场分解方法、存储介质及计算设备 |
| US11681043B2 (en) | 2019-09-03 | 2023-06-20 | Saudi Arabian Oil Company | Diffraction imaging using pseudo dip-angle gather |
| CN112578455A (zh) * | 2019-09-30 | 2021-03-30 | 中国石油化工股份有限公司 | 一种空间波数混合域地震波场多方向分解方法及系统 |
| US11313988B2 (en) | 2019-12-13 | 2022-04-26 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction imaging |
| WO2021141903A1 (en) * | 2020-01-09 | 2021-07-15 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging |
| US11402529B2 (en) | 2020-01-09 | 2022-08-02 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging |
| US11467303B2 (en) | 2020-03-09 | 2022-10-11 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using a post-stack seismic diffraction imaging condition |
| CN113419282A (zh) * | 2021-07-19 | 2021-09-21 | 中国矿业大学(北京) | 垂直地震剖面波场分离方法和系统 |
| CN115755175A (zh) * | 2022-11-18 | 2023-03-07 | 成都理工大学 | 基于Hilbert变换的复合波场弹性波逆时偏移成像方法 |
| CN115755175B (zh) * | 2022-11-18 | 2024-06-11 | 成都理工大学 | 基于Hilbert变换的复合波场弹性波逆时偏移成像方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109307890A (zh) | 基于上下行波场分解的逆时偏移方法及系统 | |
| CN105974470B (zh) | 一种多分量地震资料最小二乘逆时偏移成像方法及系统 | |
| Zhang et al. | Robust source-independent elastic full-waveform inversion in the time domain | |
| US8614930B2 (en) | System and method for seismic data modeling and migration | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN105549081A (zh) | 各向异性介质共炮域高斯束偏移成像方法 | |
| CN102890290A (zh) | 一种起伏地表条件下的叠前深度偏移方法 | |
| CN110471113B (zh) | 基于非稳态地震资料的反演动校正方法、装置及存储介质 | |
| CN103995288A (zh) | 一种高斯束叠前深度偏移方法及装置 | |
| US11573342B2 (en) | Systems and methods to enhance 3-D prestack seismic data based on non-linear beamforming in the cross-spread domain | |
| CN103308941B (zh) | 一种基于任意广角波动方程的成像方法及装置 | |
| CN103926623A (zh) | 一种压制逆时偏移低频噪音的方法 | |
| CN108845357A (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
| Araujo et al. | Symplectic scheme and the Poynting vector in reverse-time migration | |
| CN105676291A (zh) | 一种基于同相轴优化追踪的多次波匹配衰减方法 | |
| CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
| CN105301648A (zh) | 一种获取共反射面元叠加参数的方法 | |
| CN105510975A (zh) | 提高地震数据信噪比的方法及装置 | |
| CN110967734B (zh) | 基于快速傅立叶变换的虚源重构方法及系统 | |
| CN105242313A (zh) | 一种弹性波逆时偏移极性反转的校正方法及系统 | |
| CN106353798A (zh) | 多分量联合高斯束叠前逆时偏移成像方法 | |
| Wang et al. | An explicit stabilization scheme for Q-compensated reverse time migration | |
| CN107340537A (zh) | 一种p-sv转换波叠前逆时深度偏移的方法 | |
| He et al. | Source-free P-SV converted-wave reverse-time migration using first-order velocity-dilatation-rotation equations | |
| CN106950600A (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190205 |
|
| RJ01 | Rejection of invention patent application after publication |