CN105242313B - 一种弹性波逆时偏移极性反转的校正方法及系统 - Google Patents
一种弹性波逆时偏移极性反转的校正方法及系统 Download PDFInfo
- Publication number
- CN105242313B CN105242313B CN201510561256.4A CN201510561256A CN105242313B CN 105242313 B CN105242313 B CN 105242313B CN 201510561256 A CN201510561256 A CN 201510561256A CN 105242313 B CN105242313 B CN 105242313B
- Authority
- CN
- China
- Prior art keywords
- mrow
- wave
- msub
- mover
- vector
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种弹性波逆时偏移极性反转的校正方法及系统。所述校正方法包括:步骤a:利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;步骤b:根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;步骤c:对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;步骤d:利用震源波场和检波点波场的纵波分量和横波分量做互相关分别得到PS偏移剖面和SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面。本发明易于实现,计算量小,不需要先验信息,分辨率较高。
Description
技术领域
本发明属于地震波信息处理技术领域,尤其涉及一种弹性波逆时偏移极性反转的校正方法及系统。
背景技术
极性反转是弹性波逆时偏移中经常遇到的问题之一。由于反射转换波(PS波和SP波)的反射系数是关于入射角的奇函数,因此在法向入射点的两侧,反射转换波的极性会产生反转。在进行逆时偏移成像时,将多炮的PS剖面或是SP剖面进行叠加后,极性反转会造成干涉相消,从而破坏偏移同相轴,使得偏移剖面难以解释。
极性反转的校正目前已有多种方法。由于反射转换波的极性与入射角度有关,因此一个最直接的方法便是在角度域进行极性反转的校正(Rosales和Rickett,2001;Rosales等,2008;Lu等,2010)。另外,Balch和Erdemir(1994)通过求取反射PS波波形来得到P波入射角,从而校正PS波的极性反转。Du等(2012)利用坡印廷矢量求取出校正因子,并设计滤波方法来提高校正的精度。Duan(2015)提出了一种转换波成像条件,可以自动校正极性反转。
综上所述,现有极性反转的校正方法至少存在以下技术问题:
1、角度域极性校正的方法实现起来非常复杂,并且计算量特别庞大;
2、Balch和Erdemir(1994)的方法实现起来比较复杂,利用了射线近似,很难应用到复杂介质的弹性波逆时偏移中;
3、Du等(2012)利用坡印廷矢量的方法虽然比较简单有效,但是由于坡印廷矢量计算传播方向的精度有限并且只能估算主要方向,因此极性校正后的偏移剖面分辨率受到了一定的限制;
4、Duan等(2015)的方法需要预先给出地下界面的法线方向,这对于复杂介质比较难以实现。
发明内容
本发明提供了一种弹性波逆时偏移极性反转的校正方法及系统,旨在解决现有的极性反转校正方法存在的计算量大、难以应用到复杂介质的弹性波逆时偏移中、极性校正后的偏移剖面分辨率受到限制以及对于复杂介质难以实现等技术问题。
本发明是这样实现的,一种弹性波逆时偏移极性反转的校正方法,包括:
步骤a:利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;
步骤b:根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;
步骤c:对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;
步骤d:利用震源波场和检波点波场的纵波分量和横波分量做互相关,分别得到PS偏移剖面和SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面。
本发明实施例采取的技术方案还包括:在所述步骤b中,所述根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子还包括:根据纵横波分离的弹性波方程分别得到检波点波场和震源波场的纵波偏振矢量和横波偏振矢量;弹性波方程表示为:
其中,U表示位移波场矢量,α表示纵波的传播速度,β表示横波的传播速度;此方程可以分解成以下两个方程:
其中,UP=(uPx,uPz)表示的是矢量纵波波场,US=(uSx,uSz)表示的矢量横波波场,U=UP+US;方程(6)和方程(7)的左边分别等于纵波和横波的加速度场,分别用和表示,和便是纵、横波的偏振矢量。
本发明实施例采取的技术方案还包括:在所述步骤b中,所述根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子还包括:根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量。
本发明实施例采取的技术方案还包括:所述根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量具体为:令表示纵波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量,对于震源波场得:
此时,便通过纵波的偏振矢量获得了震源波场中的纵波传播方向矢量;
对于检波点波场得:
此时,便通过纵波的偏振矢量获得了检波点波场中的纵波传播方向。
令表示横波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量,对于检波点波场得
此时,便通过横波的偏振矢量获得了检波点波场中的横波传播方向矢量。
对于震源波场得:
此时,便通过横波的偏振矢量获得了震源波场中的横波传播方向矢量。
本发明实施例采取的技术方案还包括:在所述步骤b中,所述根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子还包括:利用检波点波场的横波传播方向矢量和震源波场的纵波传播方向矢量做叉乘获得PS偏移剖面的极性反转校正因子;并利用检波点波场的纵波传播方向矢量和震源波场的横波传播方向矢量做叉乘获得SP偏移剖面的极性反转校正因子。
本发明实施例采取的技术方案还包括:所述PS偏移剖面的极性反转校正因子为:
所述SP偏移剖面的极性反转校正因子为:
其中,下角标中的Receiver表示检波点波场,Source表示震源波场。
本发明实施例采取的技术方案还包括:在所述步骤b中,所述根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子还包括:取一个二维空间窗,在窗口内对所述PS偏移剖面的极性反转校正因子和SP偏移剖面的极性反转校正因子进行滤波。
本发明实施例采取的另一技术方案为:一种弹性波逆时偏移极性反转的校正系统,包括波场重建模块、校正因子求取模块、波场分离模块和极性反转校正模块;
所述波场重建模块用于利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;
校正因子求取模块用于根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;
波场分离模块用于对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;
极性反转校正模块用于利用震源波场和检波点波场的纵波分量和横波分量做互相关,分别得到PS偏移剖面和SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面。
本发明实施例采取的技术方案还包括:所述校正因子求取模块包括偏振矢量计算单元、传播方向计算单元、校正因子计算单元和校正因子滤波单元;
所述偏振矢量计算单元用于根据纵横波分离的弹性波方程分别得到检波点波场和震源波场的纵波偏振矢量和横波偏振矢量;
所述传播方向计算单元用于根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,并根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量;
所述校正因子计算单元用于利用检波点波场的横波传播方向矢量和震源波场的纵波传播方向矢量做叉乘获得PS偏移剖面的极性反转校正因子;并利用检波点波场的纵波传播方向矢量和震源波场的横波传播方向矢量做叉乘获得SP偏移剖面的极性反转校正因子;所述校正因子滤波单元用于通过二维空间窗对PS偏移剖面和SP偏移剖面的极性反转校正因子进行滤波。
本发明实施例采取的技术方案还包括:还包括剖面叠加模块,所述剖面叠加模块用于将多炮偏移剖面进行叠加,得到最终的经过极性反转校正的偏移剖面。
本发明实施例的弹性波逆时偏移极性反转的校正方法及系统利用纵横波的偏振矢量求取极性反转校正因子,并对极性反转校正因子进行滤波处理,将滤波后的极性反转校正因子带入成像条件中,实现单炮逆时偏移剖面的极性反转校正;并将多炮偏移剖面进行叠加,得到最终的经过极性反转校正的偏移剖面。本发明易于实现,计算量小,不需要先验信息,可以应用到各种复杂介质的弹性波逆时偏移中,并且分辨率较高,能够将极性反转很好的校正过来。
附图说明
图1是本发明实施例的弹性波逆时偏移极性反转的校正方法的流程图;
图2是本发明实施例的求取极性反转校正因子的方法的流程图;
图3是入射纵波和反射横波偏振矢量与传播方向矢量示意图;
图4是图3中的DP在笛卡尔坐标系中的显示图;
图5是地堑模型图;
图6是图5所示地堑模型的单炮逆时偏移剖面;其中,图6(a)是未经过极性反转校正的PS剖面,图6(b)是利用本发明的方法进行极性反转后的PS剖面,图6(c)是基于坡印廷矢量进行极性反转校正后的PS剖面;
图7是Marmousi2模型的纵波速度分布图;
图8是图7所示的Marmousi2模型的逆时偏移叠加剖面图;其中,图8(a)是未经过极性反转校正的PS剖面,图8(b)是利用本发明的方法进行极性反转后的PS剖面,图8(c)是基于坡印廷矢量进行极性反转校正后的PS剖面;
图9是图8中偏移剖面的局部对比图;其中,图9(a)是Marmousi2模型的纵波速度分布图,图9(b)是利用本发明的方法进行极性反转后的PS剖面,图9(c)是基于坡印廷矢量进行极性反转校正后的PS剖面;
图10是本发明实施例的弹性波逆时偏移极性反转的校正系统的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
请参阅图1,是本发明实施例的弹性波逆时偏移极性反转的校正方法的流程图。本发明实施例的弹性波逆时偏移极性反转的校正方法包括以下步骤:
步骤100:利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;
步骤200:根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;
在步骤200中,弹性波逆时偏移剖面中极性的反转是由于反射转换波引起的,即PS偏移剖面的极性反转是由反射转换PS波引起的,SP偏移剖面的极性反转是由反射转换SP波引起的。因此,极性反转的关键步骤便是如何获得反射转换波的极性反转校正因子。如果用sign(PS)和sign(PS)分别表示求得的反射转换PS波的极性反转校正因子以及反射转换SP波的极性反转校正因子,那么在利用成像条件时,在每个时间步上将PS剖面和SP剖面分别乘以sign(PS)和sign(PS)便可将极性反转校正回去。本发明以互相关成像条件为例,用公式表示为:
IPS=∫SP(x,t)RS(x,t)·sign(PS)dt, (1)
ISP=∫SS(x,t)RP(x,t)·sign(SP)dt, (2)
其中,IPS和ISP分别表示PS偏移剖面和SP偏移剖面,SP(x,t)和SS(x,t)分别表示震源场的P波分量和S波分量,RP(x,t)和RS(x,t)分别表示检波点波场的P波分量和S波分量,t表示时间。本发明将方程(1)和方程(2)称为带有校正因子的成像条件,本发明实施例中只列出了互相关成像条件,对于其他成像条件,例如带有照明补偿的成像条件依然适用;
在现有技术中,Du等(2012)在利用坡印廷矢量求取反射转换波的极性反转校正因子时,得到过如下关系式:
sign(PS)=sign of FE×FS, (3)
sign(SP)=sign of FR×FS, (4)
其中,FR表示检波点波场的传播方向矢量,FS表示震源波场的传播方向矢量,signof FR×FS表示取FR×FS的符号。那么,现在问题的关键变为如何求取FR和FS。本发明利用纵、横波的偏振方向来求取这两个参数。为了清楚说明步骤200,请一并参阅图2,是本发明实施例的求取极性反转校正因子的方法的流程图。本发明实施例的求取极性反转校正因子的方法包括以下步骤:
步骤201:根据纵横波分离的弹性波方程分别得到检波点波场和震源波场的纵波偏振矢量和横波偏振矢量;
在步骤201中,本发明实施例中使用的是位移形式的弹性波方程,对于其他形式的方程也同样适用,例如速度-应力方程。在各向同性介质中,弹性波方程可以表示为:
其中,U表示位移波场矢量,α表示纵波的传播速度,β表示横波的传播速度。此方程可以分解成以下两个方程:
其中,UP=(uPx,uPz)表示的是矢量纵波波场,US=(uSx,uSz)表示的矢量横波波场,U=UP+US。方程(6)和方程(7)的左边分别等于纵波和横波的加速度场,分别用和表示,和便是纵、横波的偏振矢量。
步骤202:根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,并根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量;
在步骤202中,在各向同性介质中,纵波的偏振方向有两个,一个与传播方向相同,另一个与传播方向相反。如果能将这两个方向区分出来,保留与传播方向相同的偏振方向,校正与传播方向相反的偏振方向,那么纵波的传播方向便可以求得。具体如图3和图4所示,图3是入射纵波和反射横波偏振矢量与传播方向矢量示意图;其中,DP表示纵波的偏振矢量,DS表示横波的偏振矢量,D′S表示与横波传播方向平行的一个矢量。纵波是下行波,DP中实箭头的方向是需要保留的,虚箭头的方向是需要校正的。图4是图3中的DP在笛卡尔坐标系中的显示图示;其中,表示与纵波传播方向一致的偏振方向,表示与纵波传播方向相反的偏振方向。对于下行波来说,如果z轴的方向是垂直向下的,那么传播方向矢量的垂直分量是负值,因此可以用来区分和如果则认为此时的偏振方向是否则便认为是是需要保留的方向,是需要校正的方向。从图4可以看到和的关系为:
这说明将的两个分量都乘以负号,便将校正到了方向上。
令表示纵波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量。根据上面的分析,对于震源波场可得:
此时,便通过纵波的偏振矢量获得了震源波场中的纵波传播方向矢量。
对于反射横波来说,是其偏振方向并且垂直于其传播方向(图3中DS),而的方向平行于其传播方向(图3中D′S)。与入射纵波类似,反射横波的传播方向,便可由求得。令表示横波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量。对于上行波来说,传播方向矢量的垂直分量是负值,因此
入射横波的传播方向矢量可以用以下两式求得:
对于检波点波场,反射纵波的传播方向矢量为:
此时,便通过纵波的偏振矢量获得了检波点波场中的纵波传播方向。
步骤203:利用检波点波场的横波传播方向矢量和震源波场的纵波传播方向矢量做叉乘获得PS偏移剖面的极性反转校正因子;并利用检波点波场的纵波传播方向矢量和震源波场的横波传播方向矢量做叉乘获得SP偏移剖面的极性反转校正因子;
在步骤203中,根据方程(3),反射转换PS波的极性反转校正因子可以用下式求得:
反射转换SP波的极性反转校正因子可以用下式求得:
其中,下角标中的Receiver表示检波点波场,Source表示震源波场。
步骤204:取一个二维空间窗,在窗口内对PS偏移剖面的极性反转校正因子和SP偏移剖面的极性反转校正因子进行滤波;
在步骤204中,理论上纵波和横波的偏振方向是线性的,但是数值计算误差会使其变成非线性,如椭圆偏振;另外不同方向波场的叠加也会使偏振方向变得复杂;为了提高精度,需要对PS偏移剖面的极性反转校正因子sign(PS)和SP偏移剖面的极性反转校正因子sign(SP)进行滤波处理。
步骤300:对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;
步骤400:利用震源波场的纵波分量和检波点波场的横波分量做互相关得到PS偏移剖面,利用震源波场的横波分量和检波点波场的纵波分量做互相关得到SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面,实现单炮逆时偏移剖面的极性反转校正;
步骤500:将多炮偏移剖面进行叠加,得到最终的经过极性反转校正的偏移剖面。
为了进一步说明本发明的有效性和准确性,具体请一并参阅图5、图6、图7和图8,图5是地堑模型图;其中,α表示纵波传播速度,β表示横波传播速度。在模型(2km,0.1km)处设置一个爆炸震源,震源子波为雷克子波,主频为25赫兹。图6是图5所示地堑模型的单炮逆时偏移剖面:图6(a)是未经过极性反转校正的PS剖面,图6(b)是利用本发明的方法进行极性反转后的PS剖面,图6(c)是基于坡印廷矢量进行极性反转校正后的PS剖面。由于爆炸震源只能激发出纵波,因此只给出PS偏移剖面。从图6(a)中可以看出,在三个法向入射点的两侧,偏移同相轴都出现了极性反转。从图6(b)和图6(c)中可看出,本发明和基于坡印廷矢量进行极性反转校正的方法都能很好的将偏移剖面中出现的极性反转校正过来。
图7是一个Marmousi2模型的纵波速度分布图,横波的速度为纵波速度的震源同样为爆炸震源,主频为20Hz,子波为雷克子波。图8是图7所示的Marmousi2模型的逆时偏移叠加剖面图,图8(a)是未经过极性反转校正的PS剖面,图8(b)是利用本发明的方法进行极性反转后的PS剖面,图8(c)是基于坡印廷矢量进行极性反转校正后的PS剖面。由于爆炸震源只能激发出纵波,所以也只给出PS偏移剖面。从图8(a)中可以看出,由于极性反转的影响,使得最终偏移剖面上的同相轴变得杂乱,很多界面无法清晰成像。而图8(b)和图8(c)中的偏移剖面,同相轴光滑清晰,并且断层,褶皱等构造都清晰成像了。但对比图8(b)和图8(a)可以发现,利用本发明进行极性校正后的剖面分辨率更高些。
图9是图8中偏移剖面的局部对比图,图9(a)是Marmousi2模型的纵波速度分布图,图9(b)是利用本发明的方法进行极性反转后的PS剖面,图9(c)是基于坡印廷矢量进行极性反转校正后的PS剖面。可以发现黑色箭头所标注的界面,在图9(c)中没有显示,而在图9(b)中却很好的显示了出来。这说明本发明比基于基于坡印廷矢量进行极性反转校正的方法有更高的分辨率。
请参阅图10,是本发明实施例的弹性波逆时偏移极性反转的校正系统的结构示意图。本发明实施例的弹性波逆时偏移极性反转的校正系统包括波场重建模块、校正因子求取模块、波场分离模块、极性反转校正模块和剖面叠加模块;其中,
波场重建模块用于利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;
校正因子求取模块用于根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;其中,弹性波逆时偏移剖面中极性的反转是由于反射转换波引起的,即PS偏移剖面的极性反转是由反射转换PS波引起的,SP偏移剖面的极性反转是由反射转换SP波引起的。因此,极性反转的关键步骤便是如何获得反射转换波的极性反转校正因子。如果用sign(PS)和sign(SP)分别表示求得的反射转换PS波的极性反转校正因子以及反射转换SP波的极性反转校正因子,那么在利用成像条件时,在每个时间步上将PS剖面和SP剖面分别乘以sign(PS)和sign(SP)便可将极性反转校正回去。用公式表示为(如果使用的是互相关成像条件):
IPS=∫SP(x,t)RS(x,t)·sign(PS)dt, (1)
ISP=∫SS(x,t)RP(x,t)·sign(SP)dt, (2)
其中,IPS和ISP分别表示PS偏移剖面和SP偏移剖面,SP(x,t)和SS(x,t)分别表示震源场的P波分量和S波分量,RP(x,t)和RS(x,t)分别表示检波点波场的P波分量和S波分量,t表示时间。本发明将方程(1)和方程(2)称为带有校正因子的成像条件;
在现有技术中,Du等(2012)在利用坡印廷矢量求取反射转换波的极性反转校正因子时,得到过如下关系式:
sign(PS)=sign of FR×FS, (3)
sign(SP)=sign of FR×FS, (4)
其中,FR表示检波点波场的传播方向矢量,FS表示震源波场的传播方向矢量,signof FR×FS表示取FR×FS的符号。那么,现在问题的关键变为如何求取FR和FS。本发明利用纵、横波的偏振方向来求取这两个参数。
具体地,校正因子求取模块包括偏振矢量计算单元、传播方向计算单元、校正因子计算单元和校正因子滤波单元;其中,
偏振矢量计算单元用于根据纵横波分离的弹性波方程分别得到检波点波场和震源波场的纵波偏振矢量和横波偏振矢量;其中,在各向同性介质中,弹性波方程可以表示为:
其中,U表示位移波场矢量,α表示纵波的传播速度,β表示横波的传播速度。此方程可以分解成以下两个方程:
其中,UP=(uPx,uPz)表示的是矢量纵波波场,US=(uSx,uSx)表示的矢量横波波场,U=UP+US。方程(6)和方程(7)的左边分别等于纵波和横波的加速度场,分别用和表示,和便是纵、横波的偏振矢量。
传播方向计算单元用于根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,并根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量;其中,在各向同性介质中,纵波的偏振方向有两个,一个与传播方向相同,另一个与传播方向相反。如果能将这两个方向区分出来,保留与传播方向相同的偏振方向,校正与传播方向相反的偏振方向,那么纵波的传播方向便可以求得。具体如图3和图4所示,图3是入射纵波和反射横波偏振矢量与传播方向矢量示意图;其中,DP表示纵波的偏振矢量,DS表示横波的偏振矢量,D′S表示与横波传播方向平行的一个矢量。纵波是下行波,DP中实箭头的方向是需要保留的,虚箭头的方向是需要校正的。图4是图3中的DP在笛卡尔坐标系中的显示图示;其中,表示与纵波传播方向一致的偏振方向,表示与纵波传播方向相反的偏振方向。对于下行波来说,如果z轴的方向是垂直向下的,那么传播方向矢量的垂直分量是负值,因此可以用来区分和如果则认为此时的偏振方向是否则便认为是是需要保留的方向,是需要校正的方向。从图4可以看到和的关系为:
这说明将的两个分量都乘以负号,便将校正到了方向上。
令表示纵波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量。根据上面的分析,对于震源波场可得:
此时,便通过纵波的偏振矢量获得了震源波场中的纵波传播方向矢量。
对于反射横波来说,是其偏振方向并且垂直于其传播方向(图3中DS),而的方向平行于其传播方向(图3中D′S)。与入射纵波类似,反射横波的传播方向,便可由求得。令表示横波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量。对于上行波来说,传播方向矢量的垂直分量是负值,因此
入射横波的传播方向矢量可以用以下两式求得:
对于检波点波场得:,反射纵波的传播方向矢量为:
此时,便通过纵波的偏振矢量获得了检波点波场中的纵波传播方向。
校正因子计算单元用于利用检波点波场的横波传播方向矢量和震源波场的纵波传播方向矢量做叉乘获得PS偏移剖面的极性反转校正因子;并利用检波点波场的纵波传播方向矢量和震源波场的横波传播方向矢量做叉乘获得SP偏移剖面的极性反转校正因子;其中,根据方程(3),反射转换PS波的极性反转校正因子可以用下式求得:
反射转换SP波的极性反转校正因子可以用下式求得:
其中,下角标中的Receiver表示检波点波场,Source表示震源波场。
校正因子滤波单元用于通过二维空间窗对PS偏移剖面和SP偏移剖面的极性反转校正因子进行滤波;其中,理论上纵波和横波的偏振方向是线性的,但是数值计算误差会使其变成非线性,如椭圆偏振;另外不同方向波场的叠加也会使偏振方向变得复杂;为了提高精度,需要对PS偏移剖面的极性反转校正因子sign(PS)和SP偏移剖面的极性反转校正因子sign(SP)进行滤波处理。
波场分离模块用于对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;
极性反转校正模块用于利用震源波场的纵波分量和检波点波场的横波分量做互相关得到PS偏移剖面,利用震源波场的横波分量和检波点波场的纵波分量做互相关得到SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面;
剖面叠加模块用于将多炮偏移剖面进行叠加,得到最终的经过极性反转校正的偏移剖面。
本发明实施例的弹性波逆时偏移极性反转的校正方法及系统利用纵横波的偏振矢量求取极性反转校正因子,并对极性反转校正因子进行滤波处理,将滤波后的极性反转校正因子带入成像条件中,实现单炮逆时偏移剖面的极性反转校正;并将多炮偏移剖面进行叠加,得到最终的经过极性反转校正的偏移剖面。本发明易于实现,计算量小,不需要先验信息,可以应用到各种复杂介质的弹性波逆时偏移中,并且分辨率较高,能够将极性反转很好的校正过来。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种弹性波逆时偏移极性反转的校正方法,包括:
步骤a:利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;
步骤b:根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;
步骤c:对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;
步骤d:利用震源波场的纵波分量和检波点波场的横波分量做互相关得到PS偏移剖面,利用震源波场的横波分量和检波点波场的纵波分量做互相关得到SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面;
其中:所述步骤b包括:根据纵横波分离的弹性波方程分别得到检波点波场和震源波场的纵波偏振矢量和横波偏振矢量;根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,并根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量;利用检波点波场的横波传播方向矢量和震源波场的纵波传播方向矢量做叉乘获得PS偏移剖面的极性反转校正因子;并利用检波点波场的纵波传播方向矢量和震源波场的横波传播方向矢量做叉乘获得SP偏移剖面的极性反转校正因子;取一个二维空间窗,在窗口内对PS偏移剖面的极性反转校正因子和SP偏移剖面的极性反转校正因子进行滤波。
2.根据权利要求1所述的弹性波逆时偏移极性反转的校正方法,其特征在于,所述弹性波方程表示为:
<mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>U</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>&alpha;</mi>
<mn>2</mn>
</msup>
<mo>&dtri;</mo>
<mrow>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&beta;</mi>
<mn>2</mn>
</msup>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,U表示位移波场矢量,α表示纵波的传播速度,β表示横波的传播速度;此方程可以分解成以下两个方程:
<mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>U</mi>
<mi>P</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>&alpha;</mi>
<mn>2</mn>
</msup>
<mo>&dtri;</mo>
<mrow>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>U</mi>
<mi>S</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msup>
<mi>&beta;</mi>
<mn>2</mn>
</msup>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mo>&dtri;</mo>
<mo>&times;</mo>
<mi>U</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,UP=(uPx,uPz)表示的是矢量纵波波场,uPx表示纵波波场的水平分量,uPz表示纵波波场的垂直分量,US=(uSx,uSz)表示的矢量横波波场,uSx表示横波波场的水平分量,uSz表示横波波场的垂直分量,U=UP+US;方程(6)和方程(7)的左边分别等于纵波和横波的加速度波场,分别用 和表示,表示纵波加速度波场的水平分量,表示纵波加速度波场的垂直分量,表示横波加速度波场的水平分量,表示横波加速度波场的垂直分量,和便是纵、横波的偏振矢量。
3.根据权利要求2所述的弹性波逆时偏移极性反转的校正方法,其特征在于,所述根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量具体为:令表示纵波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量,对于震源波场得:
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Px</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Px</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Pz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Px</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Px</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Pz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
此时,便通过纵波的偏振矢量获得了震源波场中的纵波传播方向矢量;
对于检波点波场得:
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Px</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Px</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Pz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Px</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Px</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Pz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Pz</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Pz</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
此时,便通过纵波的偏振矢量获得了检波点波场中的纵波传播方向;
令表示横波的传播方向矢量,和分别表示该方向矢量的水平分量和垂直分量,对于检波点波场得
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sx</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sz</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Sx</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sx</mi>
</msub>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sx</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Sz</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sx</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Sx</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
此时,便通过横波的偏振矢量获得了检波点波场中的横波传播方向矢量;
对于震源波场得:
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sx</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sz</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Sx</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sx</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sx</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sz</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mi>Sz</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>-</mo>
<mover>
<mi>u</mi>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
</mrow>
<mi>Sx</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>if</mi>
<msub>
<mover>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>.</mo>
<mo>.</mo>
</mrow>
</mover>
<mi>Sx</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
此时,便通过横波的偏振矢量获得了震源波场中的横波传播方向矢量。
4.根据权利要求3所述的弹性波逆时偏移极性反转的校正方法,其特征在于,所述PS偏移剖面的极性反转校正因子为:
<mrow>
<msup>
<mi>sign</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mi>S</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mi> </mi>
<mi>o</mi>
<mi>f</mi>
<mi> </mi>
<msubsup>
<mi>U</mi>
<mi>S</mi>
<mo>*</mo>
</msubsup>
<msub>
<mo>|</mo>
<mrow>
<mi>Re</mi>
<mi>c</mi>
<mi>e</mi>
<mi>i</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msubsup>
<mi>U</mi>
<mi>P</mi>
<mo>*</mo>
</msubsup>
<msub>
<mo>|</mo>
<mrow>
<mi>S</mi>
<mi>o</mi>
<mi>u</mi>
<mi>r</mi>
<mi>c</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
所述SP偏移剖面的极性反转校正因子为:
<mrow>
<msup>
<mi>sign</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mi>P</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mi> </mi>
<mi>o</mi>
<mi>f</mi>
<mi> </mi>
<msubsup>
<mi>U</mi>
<mi>P</mi>
<mo>*</mo>
</msubsup>
<msub>
<mo>|</mo>
<mrow>
<mi>Re</mi>
<mi>c</mi>
<mi>e</mi>
<mi>i</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msubsup>
<mi>U</mi>
<mi>S</mi>
<mo>*</mo>
</msubsup>
<msub>
<mo>|</mo>
<mrow>
<mi>S</mi>
<mi>o</mi>
<mi>u</mi>
<mi>r</mi>
<mi>c</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
其中,表示纵波波场的传播方向矢量,表示横波波场的传播方向矢量,下角标中的Receiver表示检波点波场,Source表示震源波场。
5.一种弹性波逆时偏移极性反转的校正系统,其特征在于,包括波场重建模块、校正因子求取模块、波场分离模块和极性反转校正模块;
所述波场重建模块用于利用各向同性介质中的弹性波方程进行波场重建,得到检波点波场和震源波场;
校正因子求取模块用于根据检波点波场和震源波场的偏振矢量分别求取PS偏移剖面和SP偏移剖面对应的极性反转校正因子;
波场分离模块用于对震源波场和检波点波场进行纵横波的分离,分别得出震源波场和检波点波场的纵波分量和横波分量;
极性反转校正模块用于利用震源波场和检波点波场的纵波分量和横波分量做互相关,分别得到PS偏移剖面和SP偏移剖面,并将PS偏移剖面和SP偏移剖面分别乘以对应的极性反转校正因子,得到极性反转校正后的PS偏移剖面和SP偏移剖面;
其中:所述校正因子求取模块包括偏振矢量计算单元、传播方向计算单元、校正因子计算单元和校正因子滤波单元:
所述偏振矢量计算单元用于根据纵横波分离的弹性波方程分别得到检波点波场和震源波场的纵波偏振矢量和横波偏振矢量;
所述传播方向计算单元用于根据检波点波场和震源波场的纵波偏振矢量分别获得检波点波场和震源波场的纵波传播方向矢量,并根据检波点波场和震源波场的横波偏振矢量分别获得检波点波场和震源波场的横波传播方向矢量;
所述校正因子计算单元用于利用检波点波场的横波传播方向矢量和震源波场的纵波传播方向矢量做叉乘获得PS偏移剖面的极性反转校正因子;并利用检波点波场的纵波传播方向矢量和震源波场的横波传播方向矢量做叉乘获得SP偏移剖面的极性反转校正因子;
所述校正因子滤波单元用于通过二维空间窗对PS偏移剖面和SP偏移剖面的极性反转校正因子进行滤波。
6.根据权利要求5所述的弹性波逆时偏移极性反转的校正系统,其特征在于,还包括剖面叠加模块,所述剖面叠加模块用于将多炮偏移剖面进行叠加,得到最终的经过极性反转校正的偏移剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510561256.4A CN105242313B (zh) | 2015-09-06 | 2015-09-06 | 一种弹性波逆时偏移极性反转的校正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510561256.4A CN105242313B (zh) | 2015-09-06 | 2015-09-06 | 一种弹性波逆时偏移极性反转的校正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105242313A CN105242313A (zh) | 2016-01-13 |
CN105242313B true CN105242313B (zh) | 2017-11-07 |
Family
ID=55040013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510561256.4A Active CN105242313B (zh) | 2015-09-06 | 2015-09-06 | 一种弹性波逆时偏移极性反转的校正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105242313B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107340537A (zh) * | 2016-05-03 | 2017-11-10 | 中国石油化工股份有限公司 | 一种p-sv转换波叠前逆时深度偏移的方法 |
CN108181653B (zh) * | 2018-01-16 | 2019-11-19 | 东北石油大学 | 针对vti介质逆时偏移方法、设备及介质 |
CN110389377B (zh) * | 2018-04-23 | 2020-11-27 | 中国海洋大学 | 基于波形互相关系数相乘的微震偏移成像定位方法 |
CN110941011B (zh) * | 2018-09-25 | 2021-12-24 | 中国石油化工股份有限公司 | 直达p波振幅提取方法及装置 |
CN111221037B (zh) * | 2020-01-21 | 2021-07-23 | 中国石油大学(华东) | 解耦弹性逆时偏移成像方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173988A (zh) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | 一种确定地下油气藏构造的方法 |
CN102156296A (zh) * | 2011-04-19 | 2011-08-17 | 中国石油大学(华东) | 地震多分量联合弹性逆时偏移成像方法 |
CN103091710A (zh) * | 2013-01-15 | 2013-05-08 | 中国石油天然气股份有限公司 | 一种逆时偏移成像方法及装置 |
CN104122585A (zh) * | 2014-08-08 | 2014-10-29 | 中国石油大学(华东) | 基于弹性波场矢量分解与低秩分解的地震正演模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6819628B2 (en) * | 2003-04-07 | 2004-11-16 | Paradigm Geophysical (Luxembourg) S.A.R.L. | Wave migration by a krylov space expansion of the square root exponent operator, for use in seismic imaging |
WO2014172011A2 (en) * | 2013-04-16 | 2014-10-23 | Exxonmobil Upstream Research Company | Seismic velocity model updating and imaging with elastic wave imaging |
-
2015
- 2015-09-06 CN CN201510561256.4A patent/CN105242313B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173988A (zh) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | 一种确定地下油气藏构造的方法 |
CN102156296A (zh) * | 2011-04-19 | 2011-08-17 | 中国石油大学(华东) | 地震多分量联合弹性逆时偏移成像方法 |
CN103091710A (zh) * | 2013-01-15 | 2013-05-08 | 中国石油天然气股份有限公司 | 一种逆时偏移成像方法及装置 |
CN104122585A (zh) * | 2014-08-08 | 2014-10-29 | 中国石油大学(华东) | 基于弹性波场矢量分解与低秩分解的地震正演模拟方法 |
Non-Patent Citations (4)
Title |
---|
polarity reversal correction for elastic reverse time migration;Qizhen Du 等;《Geophysics》;20120331;第77卷(第2期);第31-41页 * |
polarity reversal correction for multicomponent joint elastic reverse time migration;Qizhen Du 等;《73rd EAGE conference and exhibition incorporating SPE EUROPEC》;20111231 * |
双复杂条件下弹性波非倾斜叠加精确束偏移方法研究;黄建平 等;《石油物探》;20150131;第54卷(第1期);第56-63页 * |
弹性波高斯束叠前深度偏移;徐少波 等;《石油地球物理勘探》;20140430;第49卷(第2期);第259-265页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105242313A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105242313B (zh) | 一种弹性波逆时偏移极性反转的校正方法及系统 | |
US20200271802A1 (en) | Determining a component of a wave field | |
Day | FINITE-ELEMENT ANALYSIS OF SEISMIC SCATTERING PROBLEMS. | |
CN111164462B (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
US20170299745A1 (en) | Prestack egs migration method for seismic wave multi-component data | |
CN103149586B (zh) | 一种倾斜层状粘弹性介质中波场正演模拟方法 | |
CN108108331B (zh) | 一种基于拟空间域弹性波方程的有限差分计算方法 | |
CN109307890A (zh) | 基于上下行波场分解的逆时偏移方法及系统 | |
CN108051855B (zh) | 一种基于拟空间域声波方程的有限差分计算方法 | |
CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN105629299A (zh) | 角度域叠前深度偏移的走时、角度表获取方法及成像方法 | |
CN105425298B (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
Zhang et al. | Shot-and angle-domain wave-equation traveltime inversion of reflection data: Synthetic and field data examples | |
Wang et al. | Accurately stable Q-compensated reverse-time migration scheme for heterogeneous viscoelastic media | |
CN106353798A (zh) | 多分量联合高斯束叠前逆时偏移成像方法 | |
CN107340537A (zh) | 一种p-sv转换波叠前逆时深度偏移的方法 | |
CN110488354A (zh) | 一种q补偿的起伏地表棱柱波与一次波联合最小二乘逆时偏移成像方法 | |
CN104898162B (zh) | 地质勘探中的裂缝检测方法 | |
CN111352153B (zh) | 一种基于瞬时相位互相关加权的微地震干涉定位方法 | |
CN110376646B (zh) | 一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法 | |
CN108363097A (zh) | 一种地震资料偏移成像方法 | |
Stolk et al. | Artifacts in Kirchhoff common image gathers | |
CN107728196A (zh) | 获取角度域共成像道集的方法及系统 | |
CN106154284A (zh) | 一种对大气湍流进行补偿的单光臂关联成像方法 | |
Tang et al. | A fast RTM implementation in TTI media |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |