CN102116870A - 弹性波高斯束叠前深度偏移技术 - Google Patents
弹性波高斯束叠前深度偏移技术 Download PDFInfo
- Publication number
- CN102116870A CN102116870A CN2011100366315A CN201110036631A CN102116870A CN 102116870 A CN102116870 A CN 102116870A CN 2011100366315 A CN2011100366315 A CN 2011100366315A CN 201110036631 A CN201110036631 A CN 201110036631A CN 102116870 A CN102116870 A CN 102116870A
- Authority
- CN
- China
- Prior art keywords
- omega
- integral
- wave
- gaussian beam
- imaging
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供弹性波高斯束叠前深度偏移技术,步骤为1)根据选择的高斯束初始宽度,通过高斯窗将单炮记录划分为一系列重叠的区域,每个区域的中心为一个beam center;2)将多分量地震记录利用局部倾斜叠加分解为不同波型的局部平面波分量;3)分别从震源和beam center处试射弹性动力学高斯束,计算出每条高斯束有效范围内的走时、振幅信息;4)根据推导PP波和PS波的成像公式,利用高斯束的走时、振幅、极性信息以及地表的局部倾角信息计算对应的某个beam center的成像值;5)重复3)、4)步的运算,得到单炮的成像结果,最终的偏移结果是所有单炮成像结果的叠加;该方法可以对多次波至进行成像,具有较高的精度,并且还兼顾计算效率和灵活性,是一种具有较高应用价值的多分量地震成像方法。
Description
所属技术领域
本发明属于多分量地震资料处理领域,是一种适用于多波多分量地震数据的深度域成像方法。
背景技术
随着物探技术装备水平的提高,多波(多分量)地震勘探日益受到工业界的重视,并已逐渐进入工业化生产。多分量地震技术在非均质性储层的含油气预测、油气藏的精细描述和动态监测中已显示出独特的优势和巨大的应用潜力。然而,目前已有的适用于多分量地震数据的成像方法,还不够完善,往往存在成像精度以及计算效率之间的矛盾。
在传统的多分量地震数据处理过程中,往往首先通过波场分离方法将地震数据分为纵(P)波和横(S)波,然后利用传统的声波偏移成像方法,对分离后的P波和S波分别进行成像。然而,上述处理方法没有充分利用弹性波场的矢量特性;此外,在波场分离的过程中往往不能将不同波型的能量完全分离,残余的非本型波能量,往往导致成像结果中的大量噪声串扰,进而严重影响成像效果。
基于矢量波场延拓的弹性波地震成像方法,是解决上述问题的有效途径。由于常规单程波成像方法无法正确描述耦合的弹性波传播过程,现有的弹性波成像方法都是在时间域进行的,其大致可以分为两类:一类是弹性波逆时偏移,其基于弹性波波动方程,将多分量地震数据作为边界条件进行逆时外推,然后利用激发时成像条件或者互相关成像条件求取成像值。Sun和McMechan(1986)首先将逆时偏移应用到弹性波偏移中,此后,Chang和McMechan(1994),Botelho和Stoffa(1991),Yan和Sava(2008)以及杜启振和秦童(2009)等人将其进行推广并应用到三维以及各向异性介质的偏移成像中。另一类是弹性波Kirchhoff偏移,其通过计算本型波以及转换波的走时、振幅以及极性信息,随后沿着所选择波型的走时轨迹对波场能量进行叠加并求取成像值。Kuo和Dai(1984)基于均匀各向同性介质中位移格林函数的远场近似,首先提出了该方法。此后,Sena和Toksoz(1993),Hokstad(2000)以及Druzhinin(2003)等人对其进行了发展。上述两类方法,具有各自的优势及不足,其中弹性波逆时偏移虽然具有极高的成像精度,但是其计算量非常巨大,而且其成像结果往往伴有大量的低频干扰以及转换波成像时难以解决的极性反转问题。弹性波Kirchhoff偏移,虽然具有灵活高效的特点,但是其难以对多次波至进行成像,致使其对地下复杂地质构造的成像精度不高。
发明内容
本发明的目的在于提供一种弹性波高斯束叠前深度偏移技术,是一种兼顾精度与效率的弹性波地震成像方法,以多分量地震数据为输入,利用弹性波动力学高斯束的振幅、走时以及极性信息来进行矢量波场的延拓,并且利用波场的矢量投影来提取纵横波的成像值。
本发明所采用技术方案是:弹性波高斯束叠前深度偏移技术,实现步骤为:
(1)根据选择的高斯束初始宽度,通过高斯窗将单炮记录划分为一系列重叠的区域,每个区域的中心为一个beam center;
(2)将多分量地震记录利用局部倾斜叠加分解为不同波型的局部平面波分量;
(3)分别从震源和beam center处试射弹性动力学高斯束,计算出每条高斯束有效范围内的走时、振幅信息;
(4)根据推导PP波和PS波的成像公式,利用高斯束的走时、振幅、极性信息以及地表的局部倾角信息计算对应的某个beam center的成像值;
(5)重复(3)、(4)步的运算,得到单炮的成像结果,最终的偏移结果是所有单炮成像结果的叠加。
其中步骤(3)中的计算方法是:
在S所描述的观测面上,假设由震源xs激发,接收点xr接收到的两分量弹性波地震记录为ui(xr;ω),则反向延拓的弹性波位移场um(x;xr;ω)可以通过Kirchhoff-Helmholtz积分来表示:
其中*代表复共轭;Glm(x;xr;ω)为位移格林张量,其代表xr处l方向单位体力所造成的x处位移的m方向的分量;ti(xr)为xr处应力;∑im(x;xr)为应力格林张量;
若假设S为自由地表,则根据自由应力边界条件t(x;ω)=0,从而式(1)变为:
并取Glm(x;xr;ω)的高频近似解,代入式(2),得
(3)
权值具有以下形式
高斯束的初始波前为平面的,将多分量地震记录分解为不同波型不同初始方向的局部平面波,然后利用相对应高斯束进行波场延拓,其具体实现过程如下:
对地震记录加入一系列高斯窗,将其分为一系列重叠的局部区域,高斯窗函数具有如下性质
其中,x=L为高斯窗的中线,也是波束中心L的水平坐标;ΔL为其水平间隔,将式(5)代入式(3)得
其中,为对不同波型多分量地震记录的加窗局部倾斜叠加
权值此时为
(11)
接下来根据P波以及S波的传播特点,结合式(8),式(9)所表示的反向延拓的接收波场,定义如下成像公式
其中,Ipp(x)为P-P单炮成像值,Ips(x)为P-S单炮成像值,将所有炮记录进行偏移并叠加,便得到最终的弹性波成像结果。
本发明的有益效果是:该方法不但可以对多次波至进行成像,具有较高的精度,并且还兼顾计算效率和灵活性,是一种具有较高应用价值的多分量地震成像方法。
附图说明
图1是具体实施方式的平层模型偏移试算效果图。图中(a)水平分量单炮记录;(b)垂直分量单炮记录;(c)垂直分量记录P-P波单炮成像结果;(d)两分量记录的P-P单炮成像结果;(e)水平分量记录P-S波单炮成像结果;(f)经过极性校正后的两分量记录的P-S成像结果;(g)P-P波成像叠加结果;(h)P-S波成像叠加结果。
图2是P-S转换波在反射界面处的偏振图。
图3是实施例的复杂模型偏移试算结果图。其中(a)P波速度模型;(b)S波速度模型;(c)X分量单炮记录;(d)Z分量单炮记录;(e)P-P波成像叠加结果;(f)P-S波成像叠加结果。
具体实施方式
弹性波高斯束叠前深度偏移技术,实现步骤为:
(1)根据选择的高斯束初始宽度,通过高斯窗将单炮记录划分为一系列重叠的区域,每个区域的中心为一个beam center;
(2)将多分量地震记录利用局部倾斜叠加分解为不同波型的局部平面波分量;
(3)分别从震源和beam center处试射弹性动力学高斯束,计算出每条高斯束有效范围内的走时、振幅信息;
(4)根据推导PP波和PS波的成像公式,利用高斯束的走时、振幅、极性信息以及地表的局部倾角信息计算对应的某个beam center的成像值;
(5)重复(3)、(4)步的运算,得到单炮的成像结果,最终的偏移结果是所有单炮成像结果的叠加。
其中步骤(3)中的计算方法是:
在S所描述的观测面上,假设由震源xs激发,接收点xr接收到的两分量弹性波地震记录为ui(xr;ω),则反向延拓的弹性波位移场um(x;xr;ω)可以通过Kirchhoff-Helmholtz积分来表示:
其中*代表复共轭;Glm(x;xr;ω)为位移格林张量,其代表xr处l方向单位体力所造成的x处位移的m方向的分量;ti(xr)为xr处应力;∑im(x;xr)为应力格林张量;
若假设S为自由地表,则根据自由应力边界条件t(x;ω)=0,从而式(1)变为:
并取Glm(x;xr;ω)的高频近似解,代入式(2),得
(3)
高斯束的初始波前为平面的,将多分量地震记录分解为不同波型不同初始方向的局部平面波,然后利用相对应高斯束进行波场延拓,其具体实现过程如下:
对地震记录加入一系列高斯窗,将其分为一系列重叠的局部区域,高斯窗函数具有如下性质
其中,x=L为高斯窗的中线,也是波束中心L的水平坐标;ΔL为其水平间隔,将式(5)代入式(3)得
引入相移校正因子,将通过L处出射的高斯束来表示
(11)
接下来根据P波以及S波的传播特点,结合式(8),式(9)所表示的反向延拓的接收波场,定义如下成像公式
其中,Ipp(x)为P-P单炮成像值,Ips(x)为P-S单炮成像值,将所有炮记录进行偏移并叠加,便得到最终的弹性波成像结果。
首先通过一个三层水平层状介质模型来测试本发明的正确性。模型网格为301×401,纵横向采样间隔都为10m,两个水平反射层分别位于1400m以及2700m处。随着深度的增大,P波速度分别为2500m/s,2900m/s,3200m/s,S波速度分别为1500m/s,1800m/s,2600m/s,密度设定为常数。模拟弹性波记录主频为30Hz,共有51炮,每炮151道,道间距为20m,时间采样点数为1500,采样间隔为2ms。图1a,图1b分别为第26炮单炮记录的X,Z分量。图1c为仅对Z分量记录的P-P单炮成像结果,可以明显看到Z分量记录中的S波成分所导致的串扰噪声(箭头所指处),图1d为采用式(15)、(23)对两分量记录加权后得到的P-P单炮成像结果,此时串扰噪声得到了有效的压制。图1e为X分量记录的P-S单炮成像结果,图1f为采用式(15),(23)对两分量记录加权后的P-S单炮成像结果,可以看到,相比于图1e,图1f中的串扰噪声得到了有效的压制(箭头所指处),且成像结果中的极性反转现象得到了有效的校正(椭圆标识处)。图1g,1h分别为对所有炮偏移叠加后的P-P,P-S成像结果,可以两者均清晰的反映了界面的真实形态,且由于S波传播速度较慢,使得P-S成像结果具有相对更高的分辨率。
在对转换波进行成像时,往往出现转换波偏振方向不同而导致的成像剖面中的极性反转现象,严重影响叠加成像的效果。针对上述问题,本发明提出一种校正方法。图2显示了P-S转换波的传播过程,可以看到由于入射P波具有不同符号的入射角(α1>0,α2<0,以逆时针为正),使得O1点产生的P-S转换波具有正水平方向(以向右为正)的位移分量,O2产生的P-S转换波具有负水平方向的位移分量,最终使R1,R2点接收到的X分量地震记录具有相反的极性。也就是说X分量记录中转换S波的极性同P波入射角的正负有关,因此可以根据P波入射角α的正负直接对成像结果进行校正,为此引入符号函数,得到式(14)所示的P-S转换波成像公式。
利用一复杂构造模型来测试本发明对复杂地质模型的成像效果。速度模型如图3a,3b所示,其网格大小为1000×550,纵横向采样间隔都为10m,5m,同样将密度设定为常数。模型共分为7层,其中包含断层、洼陷等复杂构造。模拟弹性波记录主频为25Hz,共有200炮,每炮361道,道间距为10m,时间采样点数为1200,采样间隔为2.1ms。图3c,3d分别为第1炮单炮记录的X,Z分量,可以看到由于地下构造复杂导致记录中的波场信息也非常复杂。图3e,3f分别为最终的P-P,P-S成像结果,可以看到两者均良好的反映了速度模型的构造情况,模型中的断层、洼陷以及起伏的底层均得到了较好的成像,证明了本文方法的有效性。
Claims (3)
1.弹性波高斯束叠前深度偏移技术,其特征在于实现步骤为:
(1)根据选择的高斯束初始宽度,通过高斯窗将单炮记录划分为一系列重叠的区域,每个区域的中心为一个beam center;
(2)将多分量地震记录利用局部倾斜叠加分解为不同波型的局部平面波分量;
(3)分别从震源和beam center处试射弹性动力学高斯束,计算出每条高斯束有效范围内的走时、振幅信息;
(4)根据推导PP波和PS波的成像公式,利用高斯束的走时、振幅、极性信息以及地表的局部倾角信息计算对应的某个beam center的成像值;
(5)重复(3)、(4)步的运算,得到单炮的成像结果,最终的偏移结果是所有单炮成像结果的叠加。
2.根据权利要求1所述的弹性波高斯束叠前深度偏移技术,其特征在于步骤(3)中的计算方法是:
在S所描述的观测面上,假设由震源xs激发,接收点xr接收到的两分量弹性波地震记录为ui(xr;ω),则反向延拓的弹性波位移场um(x;xr;ω)可以通过Kirchhoff-Helmholtz积分来表示:
其中*代表复共轭;Glm(x;xr;ω)为位移格林张量,其代表xr处l方向单位体力所造成的x处位移的m方向的分量;ti(xr)为xr处应力;∑im(x;xr)为应力格林张量;
若假设S为自由地表,则根据自由应力边界条件t(x;ω)=0,从而式(1)变为:
并取Glm(x;xr;ω)的高频近似解,代入式(2),得
(3)
高斯束的初始波前为平面的,将多分量地震记录分解为不同波型不同初始方向的局部平面波,然后利用相对应高斯束进行波场延拓,其具体实现过程如下:
对地震记录加入一系列高斯窗,将其分为一系列重叠的局部区域,高斯窗函数具有如下性质
其中,x=L为高斯窗的中线,也是波束中心L的水平坐标;ΔL为其水平间隔,将式(5)代入式(3)得
(11)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100366315A CN102116870A (zh) | 2011-02-12 | 2011-02-12 | 弹性波高斯束叠前深度偏移技术 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100366315A CN102116870A (zh) | 2011-02-12 | 2011-02-12 | 弹性波高斯束叠前深度偏移技术 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102116870A true CN102116870A (zh) | 2011-07-06 |
Family
ID=44215724
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100366315A Pending CN102116870A (zh) | 2011-02-12 | 2011-02-12 | 弹性波高斯束叠前深度偏移技术 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102116870A (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102914796A (zh) * | 2012-08-21 | 2013-02-06 | 北京多分量地震技术研究院 | 一种基于高斯束的获取纵横波偏移速度的控制方法 |
CN103995288A (zh) * | 2014-05-13 | 2014-08-20 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN104280766A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油化工股份有限公司 | 一种利用局部数据同相轴斜率的直接偏移方法 |
CN105093281A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 一种反演框架下的地震多波建模方法 |
CN105487118A (zh) * | 2014-09-28 | 2016-04-13 | 中国石油化工股份有限公司 | 一种处理叠前宽线地震数据的方法 |
CN105549081A (zh) * | 2016-01-29 | 2016-05-04 | 中国石油大学(华东) | 各向异性介质共炮域高斯束偏移成像方法 |
CN105891885A (zh) * | 2014-10-20 | 2016-08-24 | 杨继东 | 针对于地下特定地质体的叠前保幅聚焦成像技术 |
CN106353798A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 多分量联合高斯束叠前逆时偏移成像方法 |
CN106597540A (zh) * | 2016-12-30 | 2017-04-26 | 中国科学院地质与地球物理研究所 | 高斯束偏移成像方法及装置 |
CN106842300A (zh) * | 2016-12-21 | 2017-06-13 | 中国石油大学(华东) | 一种高效率多分量地震资料真振幅偏移成像方法 |
CN106896403A (zh) * | 2017-05-05 | 2017-06-27 | 中国石油集团东方地球物理勘探有限责任公司 | 弹性高斯束偏移成像方法和*** |
CN106918838A (zh) * | 2017-03-06 | 2017-07-04 | 中国科学院地质与地球物理研究所 | 起伏地表条件下高斯束偏移成像方法及装置 |
CN108363101A (zh) * | 2018-02-02 | 2018-08-03 | 西安石油大学 | 一种斜井井间地震高斯束叠前深度偏移成像方法 |
CN109655882A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 基于高斯束波场模拟的地震波正演方法和装置 |
CN110058300A (zh) * | 2018-10-30 | 2019-07-26 | 南方科技大学 | 一次波重建方法、装置、终端设备及存储介质 |
CN111999770A (zh) * | 2020-09-03 | 2020-11-27 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及*** |
CN112415590A (zh) * | 2020-10-30 | 2021-02-26 | 中国石油天然气集团有限公司 | 局部平面波域弹性波成像方法及装置 |
CN113687415A (zh) * | 2021-08-13 | 2021-11-23 | 中国石油大学(华东) | 一种基于主频近似的声波介质时空域高斯束偏移成像方法 |
-
2011
- 2011-02-12 CN CN2011100366315A patent/CN102116870A/zh active Pending
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102914796A (zh) * | 2012-08-21 | 2013-02-06 | 北京多分量地震技术研究院 | 一种基于高斯束的获取纵横波偏移速度的控制方法 |
CN104280766A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油化工股份有限公司 | 一种利用局部数据同相轴斜率的直接偏移方法 |
CN104280766B (zh) * | 2013-07-12 | 2017-03-08 | 中国石油化工股份有限公司 | 一种利用局部数据同相轴斜率的直接偏移方法 |
CN103995288B (zh) * | 2014-05-13 | 2017-02-01 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN103995288A (zh) * | 2014-05-13 | 2014-08-20 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN105093281A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 一种反演框架下的地震多波建模方法 |
CN105487118A (zh) * | 2014-09-28 | 2016-04-13 | 中国石油化工股份有限公司 | 一种处理叠前宽线地震数据的方法 |
CN105487118B (zh) * | 2014-09-28 | 2018-03-06 | 中国石油化工股份有限公司 | 一种处理叠前宽线地震数据的方法 |
CN105891885A (zh) * | 2014-10-20 | 2016-08-24 | 杨继东 | 针对于地下特定地质体的叠前保幅聚焦成像技术 |
CN106353798A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 多分量联合高斯束叠前逆时偏移成像方法 |
CN105549081A (zh) * | 2016-01-29 | 2016-05-04 | 中国石油大学(华东) | 各向异性介质共炮域高斯束偏移成像方法 |
CN105549081B (zh) * | 2016-01-29 | 2018-09-11 | 中国石油大学(华东) | 各向异性介质共炮域高斯束偏移成像方法 |
CN106842300B (zh) * | 2016-12-21 | 2018-10-30 | 中国石油大学(华东) | 一种高效率多分量地震资料真振幅偏移成像方法 |
CN106842300A (zh) * | 2016-12-21 | 2017-06-13 | 中国石油大学(华东) | 一种高效率多分量地震资料真振幅偏移成像方法 |
CN106597540B (zh) * | 2016-12-30 | 2017-10-31 | 中国科学院地质与地球物理研究所 | 高斯束偏移成像方法及装置 |
CN106597540A (zh) * | 2016-12-30 | 2017-04-26 | 中国科学院地质与地球物理研究所 | 高斯束偏移成像方法及装置 |
US10261206B2 (en) | 2016-12-30 | 2019-04-16 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Gaussian beam migration imaging method and device |
CN106918838B (zh) * | 2017-03-06 | 2018-04-03 | 中国科学院地质与地球物理研究所 | 起伏地表条件下高斯束偏移成像方法及装置 |
CN106918838A (zh) * | 2017-03-06 | 2017-07-04 | 中国科学院地质与地球物理研究所 | 起伏地表条件下高斯束偏移成像方法及装置 |
CN106896403B (zh) * | 2017-05-05 | 2019-09-13 | 中国石油天然气集团有限公司 | 弹性高斯束偏移成像方法和*** |
CN106896403A (zh) * | 2017-05-05 | 2017-06-27 | 中国石油集团东方地球物理勘探有限责任公司 | 弹性高斯束偏移成像方法和*** |
CN109655882A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 基于高斯束波场模拟的地震波正演方法和装置 |
CN108363101B (zh) * | 2018-02-02 | 2019-07-02 | 西安石油大学 | 一种斜井井间地震高斯束叠前深度偏移成像方法 |
CN108363101A (zh) * | 2018-02-02 | 2018-08-03 | 西安石油大学 | 一种斜井井间地震高斯束叠前深度偏移成像方法 |
CN110058300A (zh) * | 2018-10-30 | 2019-07-26 | 南方科技大学 | 一次波重建方法、装置、终端设备及存储介质 |
CN111999770A (zh) * | 2020-09-03 | 2020-11-27 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及*** |
CN111999770B (zh) * | 2020-09-03 | 2024-01-16 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及*** |
CN112415590A (zh) * | 2020-10-30 | 2021-02-26 | 中国石油天然气集团有限公司 | 局部平面波域弹性波成像方法及装置 |
CN113687415A (zh) * | 2021-08-13 | 2021-11-23 | 中国石油大学(华东) | 一种基于主频近似的声波介质时空域高斯束偏移成像方法 |
CN113687415B (zh) * | 2021-08-13 | 2024-02-20 | 中国石油大学(华东) | 一种基于主频近似的声波介质时空域高斯束偏移成像方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102116870A (zh) | 弹性波高斯束叠前深度偏移技术 | |
CN101630014B (zh) | 一种利用垂直地震剖面数据对各向异性介质成像的方法 | |
CN105974470B (zh) | 一种多分量地震资料最小二乘逆时偏移成像方法及*** | |
CN102053261B (zh) | 一种地震数据处理方法 | |
CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN103995288B (zh) | 一种高斯束叠前深度偏移方法及装置 | |
CN102866421B (zh) | 识别小断距断点的散射波叠前成像方法 | |
CN102176053B (zh) | 提升波动方程叠前深度偏移成像效果的方法 | |
CN109307890A (zh) | 基于上下行波场分解的逆时偏移方法及*** | |
CN109856679B (zh) | 一种各向异性介质弹性波高斯束偏移成像方法及*** | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN105182408A (zh) | 一种合成地震记录的制作方法和装置 | |
CN106291687A (zh) | 各向异性多波高斯束叠前深度偏移成像方法 | |
CN102305941A (zh) | 由叠前时间偏移直接扫描确定地层叠加品质因子方法 | |
Glinsky et al. | Numerical simulation of topographical and geological site effects. Applications to canonical topographies and Rognes hill, South East France | |
CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
CN104730579A (zh) | 一种基于表层横波速度反演的纵横波联合静校正方法 | |
CN104765064A (zh) | 一种微地震干涉成像的方法 | |
CN102053267A (zh) | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 | |
Todorovska et al. | Amplitudes, polarity and time of peaks of strong ground motion during the 1994 Northridge, California, earthquake | |
CN110568487A (zh) | 基于天然地震波形的活跃断层结构成像方法 | |
CN107807393A (zh) | 基于地震干涉法的单台站集初至波增强方法 | |
CN105242313A (zh) | 一种弹性波逆时偏移极性反转的校正方法及*** | |
CN105093279A (zh) | 针对山前带的三维地震初至波菲涅尔体层析反演方法 | |
CN106019374A (zh) | 基于反射槽波频散相似度的断层成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110706 |