CN108919356B - 一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法 - Google Patents
一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法 Download PDFInfo
- Publication number
- CN108919356B CN108919356B CN201810459119.3A CN201810459119A CN108919356B CN 108919356 B CN108919356 B CN 108919356B CN 201810459119 A CN201810459119 A CN 201810459119A CN 108919356 B CN108919356 B CN 108919356B
- Authority
- CN
- China
- Prior art keywords
- wave field
- wave
- desert
- equation
- module
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 50
- 238000013508 migration Methods 0.000 title claims abstract description 41
- 230000005012 migration Effects 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title abstract description 10
- 239000006185 dispersion Substances 0.000 claims abstract description 24
- 239000004576 sand Substances 0.000 claims abstract description 22
- 230000001788 irregular Effects 0.000 claims abstract description 16
- 230000006641 stabilisation Effects 0.000 claims abstract description 8
- 230000009466 transformation Effects 0.000 claims description 10
- 230000014509 gene expression Effects 0.000 claims description 9
- 230000000087 stabilizing effect Effects 0.000 claims description 4
- 238000004611 spectroscopical analysis Methods 0.000 claims description 3
- 239000006227 byproduct Substances 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 claims 1
- 239000003208 petroleum Substances 0.000 abstract description 2
- 230000000644 propagated effect Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
-
- 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/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法,属于石油地球物理勘探领域,输入偏移模型和炮记录;将偏移模型转换到计算域;正向延拓得到计算域波场正传波场快照;得到能够准确校正振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓方程,得到计算域波场反传的波场快照;引入一种计算域稳定化算子压制反向延拓的高频噪音;将正向延拓和反向延拓的波场快照反变换到物理域;应用物理域零相位互相关成像条件得到成像结果并输出成像结果。本发明能够准确校正不规则沙丘地表对地震波波传播的影响,同时可以校正松散的沙丘地表对地震记录能量的衰减和相位的频散,得到准确大沙漠区地震采集数据的高精度成像结果。
Description
技术领域
本发明属于石油地球物理勘探领域,具体涉及一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法。
背景技术
大沙漠地区近地表松散的沙丘地表条件和复杂的起伏地表结构给沙漠区采集数据的成像带来了极大的困难。松散的地表对地震波能量的吸收和衰减非常强烈,因此需要在成像过程中对衰减的地震波能量进行补偿。因松散的沙丘地表导致的强衰减作用,大沙漠区采集数据的初至不明显难以拾取,且数百米起伏高差的沙垄和沙丘导致初至曲线存在非常明显的扭曲畸变,因此大沙漠区采集数据存在较为突出的静校正问题,因此需要针对大沙漠采集数据研究针对不规则衰减沙丘地表的成像技术。
发明内容
针对现有技术中存在的上述技术问题,本发明提出了一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法,设计合理,克服了现有技术的不足,具有良好的效果。
为了实现上述目的,本发明采用如下技术方案:
一种稳定的大沙漠区衰减补偿逆时偏移成像系统,包括输入模块、坐标变换模块、波场正向延拓模块、波场反向延拓模块、稳定化算子模块、坐标反变换模块、互相关成像模块和输出模块;
输入模块,被配置为用于输入物理域大沙漠区偏移速度模型、Q模型和实际观测炮记录,并建立观测系统;
坐标变换模块,被配置为用于将大沙漠区的偏移速度模型和Q模型转换到计算域;
波场正向延拓模块,被配置为用于正向延拓更新波场,得到计算域波场正传波场快照;
波场反向延拓模块,被配置为用于得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,并进行逆时延拓,得到计算域波场反传波场快照;
稳定化算子模块,被配置为用于引入一种计算域稳定化算子压制波场反向延拓的高频噪音;
坐标反变换模块,被配置为用于将正向延拓的波场快照与反传延拓的波场快照反变换到物理域;
互相关成像模块,被配置为用于应用物理域零相位互相关成像条件得到成像结果;
输出模块,被配置为用于输出大沙漠区采集数据的成像结果。
此外,本发明还提到一种稳定的大沙漠区衰减补偿逆时偏移成像方法,该方法采用如上所述的稳定的大沙漠区衰减补偿逆时偏移成像系统,包括如下步骤:
步骤1:通过输入模块,输入物理域大沙漠区偏移速度模型、Q模型和实际观测炮记录,并建立观测系统;
步骤2:通过坐标变换模块,采用式(1)将大沙漠区的偏移速度模型和Q模型转换到计算域;
其中,(x,z)表示物理域的坐标表达式;(ξ,η)表示计算域的坐标表达式;
步骤3:通过波场正向延拓模块,采用波动方程进行正向延拓更新波场,在正向延拓过程中补偿振幅衰减、校正相位频散,得到计算域波场正传波场快照;
步骤4:通过波场反向延拓模块,得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,如式(8)所示,并进行逆时延拓,得到计算域波场反传波场快照;
其中,pR为反向延拓的波场快照,dobs为输入的实际观测炮记录;
步骤5:通过稳定化算子模块,引入一种计算域稳定化算子得到波场反向延拓波动方程,来压制波场反向延拓的高频噪音;
步骤6:通过坐标反变换模块,将正向延拓的波场快照与反传延拓的波场快照反变换到物理域;
步骤7:通过互相关成像模块,应用如式(19)所示的物理域零相位互相关成像条件得到成像结果,并通过输出模块,输出大沙漠区采集数据的成像结果:
其中,x=(x,z)表示物理域坐标,T表示记录时间,I表示成像结果。
优选地,在步骤3中,波动方程如式(2)所示:
其中,f表示震源,v0为零频速度,ps为正向延拓的声压场,t为时间;τ为松弛时间,由品质因子Q计算得到:A(ξ)、B(ξ,η)和C(η)分别由式(3)、(4)、(5)求得:
上式所示的波动方程中的第二项为振幅衰减项,第三项为相位频散项;当模拟振幅衰减时,第二项的符号为正号,当进行衰减补偿时,第二项的符号变为负号;在进行相位频散校正时,第三项的符号保持不变;在沙漠区不规则衰减地表条件下,振幅衰减的影响要远大于相位频散的影响,如果只对振幅衰减进行补偿,采用如下所示的波动方程进行波场正向延拓:
其中,v表示偏移速度;
采用如式(7)所示的伪谱法波场递推公式进行波场计算:
优选地,在步骤5中,波场反向延拓波动方程如式(9)所示:
为消除一阶项,定义Λt如下:
引入变量qR,使得qR满足下式:
qR(x,t)=ΛtpR(x,t) (11);
将方程(11)代入到不规则衰减沙丘地表的反向延拓波动方程中,并略去震源项可得:
的一阶偏导数和二阶偏导数为:
将方程(13)和(14)代入到方程(12)中并进行简化可得:
从方程(15)我们可看出该方程不再包含一阶时间偏导数,形成了接近于二阶声波方程的波动方程,则衰减的波场反向延拓波动方程等价于:
为压制高频噪音,引入一个稳定化项到指数项中,则方程(16)可变为:
其中,σ为规则化参数;
将方程组(17)中的两个方程结合并进行简化,可得引入稳定化算子的波场反向延拓波动方程:
本发明所带来的有益技术效果:
本发明可以在沙漠区无法准确获得高程静校正量的情况下,准确校正不规则沙丘地表对地震波波传播的影响,同时可以校正松散的沙丘地表对地震记录能量的衰减和相位的频散,在沙漠区采集数据的成像中发挥很大的作用;本发明提出了振幅衰减和相位频散校正的正向延拓衰减补偿波动方程,同时引入一种计算域的稳定化算子,提出了稳定地振幅衰减和相位频散校正的反向延拓波动方程,所有波动方程都是基于正交网格剖分的计算域中实现的,可以准确实现大沙漠区的地震数据成像。
附图说明
图1为本发明的一种稳定的大沙漠区衰减补偿逆时偏移成像方法的流程图;
图2为本发明使用的起伏地表洼陷衰减介质模型示意图。(a)物理域的速度模型;(b)物理域的Q模型;(c)计算域速度模型;(d)计算域Q模型;
图3为本发明使用的正交网格剖分图。(a)全局网格图;(b)部分放图;
图4为实际炮记录示意图。(a)本发明输入的炮记录;(b)作为对比的无衰减的炮记录;
图5为正传波场快照图。(a)本发明;(b)无衰减补偿
图6为反传波场快照图。(a)本发明;(b)无衰减补偿
图7为成像结果示意图。(a)无衰减补偿;(b)本发明;
图8为大沙漠区实际资料的偏移模型示意图。(a)偏移速度;(b)偏移Q场。
图9为大沙漠区实际资料的成像结果示意图。(a)无衰减补偿;(b)本发明;
图10为本发明中一种稳定的大沙漠区衰减补偿逆时偏移成像系统的结构示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
实施例1:
一种稳定的大沙漠区衰减补偿逆时偏移成像系统,其结构如图10所示,包括输入模块、坐标变换模块、波场正向延拓模块、波场反向延拓模块、稳定化算子模块、坐标反变换模块、互相关成像模块和输出模块;
输入模块,被配置为用于输入物理域大沙漠区偏移速度模型、Q模型和实际观测炮记录,并建立观测系统;
坐标变换模块,被配置为用于将大沙漠区的偏移速度模型和Q模型转换到计算域;
波场正向延拓模块,被配置为用于正向延拓更新波场,得到计算域波场正传波场快照;
波场反向延拓模块,被配置为用于得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,并进行逆时延拓,得到计算域波场反传波场快照;
稳定化算子模块,被配置为用于引入一种计算域稳定化算子压制波场反向延拓的高频噪音;
坐标反变换模块,被配置为用于将正向延拓的波场快照与反传延拓的波场快照反变换到物理域;
互相关成像模块,被配置为用于应用物理域零相位互相关成像条件得到成像结果;
输出模块,被配置为用于输出大沙漠区采集数据的成像结果。
实施例2:
在上述实施例的基础上,本发明还提到一种稳定的大沙漠区衰减补偿逆时偏移成像方法,其流程如图1所示,具体包括如下步骤:
步骤1:通过输入模块,输入物理域大沙漠区偏移速度模型、Q模型和实际观测炮记录,并建立观测系统;
步骤2:通过坐标变换模块,采用式(1)将大沙漠区的偏移速度模型和Q模型转换到计算域;
其中,(x,z)表示物理域的坐标表达式;(ξ,η)表示计算域的坐标表达式;
步骤3:通过波场正向延拓模块,采用波动方程进行正向延拓更新波场,在正向延拓过程中补偿振幅衰减、校正相位频散,得到计算域波场正传波场快照;
波动方程如式(2)所示:
其中,f表示震源,v0为零频速度,ps为正向延拓的声压场,t为时间;τ为松弛时间,由品质因子Q计算得到:A(ξ)、B(ξ,η)和C(η)分别由式(3)、(4)、(5)求得:
上式所示的波动方程中的第二项为振幅衰减项,第三项为相位频散项;当模拟振幅衰减时,第二项的符号为正号,当进行衰减补偿时,第二项的符号变为负号;在进行相位频散校正时,第三项的符号保持不变;在沙漠区不规则衰减地表条件下,振幅衰减的影响要远大于相位频散的影响,如果只对振幅衰减进行补偿,采用如下所示的波动方程进行波场正向延拓:
其中,v表示偏移速度;
采用如式(7)所示的伪谱法波场递推公式进行波场计算:
步骤4:通过波场反向延拓模块,得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,如式(8)所示,并进行逆时延拓,得到计算域波场反传波场快照;
其中,pR为反向延拓的波场快照,dobs为输入的实际观测炮记录;
步骤5:通过稳定化算子模块,引入一种计算域稳定化算子得到波场反向延拓波动方程,来压制波场反向延拓的高频噪音;
波场反向延拓波动方程如式(9)所示:
为消除一阶项,定义Λt如下:
引入变量qR,使得qR满足下式:
qR(x,t)=ΛtpR(x,t) (11);
将方程(11)代入到不规则衰减沙丘地表的反向延拓波动方程中,并略去震源项可得:
的一阶偏导数和二阶偏导数为:
将方程(13)和(14)代入到方程(12)中并进行简化可得:
从方程(15)我们可看出该方程不再包含一阶时间偏导数,形成了接近于二阶声波方程的波动方程,则衰减的波场反向延拓波动方程等价于:
为压制高频噪音,引入一个稳定化项到指数项中,则方程(16)可变为:
其中,σ为规则化参数;
将方程组(17)中的两个方程结合并进行简化,可得引入稳定化算子的波场反向延拓波动方程:
步骤6:通过坐标反变换模块,将正向延拓的波场快照与反传延拓的波场快照反变换到物理域;
步骤7:通过互相关成像模块,应用如式(19)所示的物理域零相位互相关成像条件得到成像结果,并通过输出模块,输出大沙漠区采集数据的成像结果:
其中,x=(x,z)表示物理域坐标,T表示记录时间,I表示成像结果。
本发明在沙漠区无法准确获得高程静校正量的情况下,准确校正不规则沙丘地表对地震波波传播的影响,同时可以校正松散的沙丘地表对地震记录能量的衰减和相位的频散。在沙漠区采集数据的成像中发挥很大的作用。
应用实验例1
本发明一种稳定的大沙漠区衰减补偿逆时偏移成像方法,应用于一个典型的洼陷不规则地表衰减模型数据,取得了理想的计算效果。输入物理域大沙漠区偏移速度模型(如图2a所示)、Q模型(如图2b所示)和实际观测炮记录(如图4a所示),并建立观测系统;生成正交曲网格(如图3所示),将大沙漠区的偏移速度模型和Q模型转换到计算域中的偏移速度模型(如图2c所示)和Q模型(如图2d所示);正向延拓更新波场,得到计算域波场正传波场快照;得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,并进行逆时延拓,得到计算域波场反传的波场快照;引入一种计算域稳定化算子压制波场反向延拓的高频噪音;将正向延拓的波场快照与反传延拓的波场快照反变换到物理域的正向延拓波场快照(如图5a所示)与反传延拓波场快照(如图6a所示),从图6a中可以看出,经过引入计算域稳定化算子压制波场反向延拓的高频噪音后,波场稳定,无高频噪音;应用物理域零相位互相关成像条件得到成像结果,并输出大沙漠区采集数据的成像结果(如图7b所示)。作为对比,给出常规无衰减炮记录(如图4b所示),从图中箭头位置所示,实际衰减炮记录能量明显衰减。常规无衰减补偿的正传波场快照和反传波场快照分别如图5b和6b所示。得到的常规无衰减补偿的成像结果如图7a所示。通过对比可以看出,本发明得到的成像结果衰减得到了很好地补偿。
应用实验例2
本发明一种稳定的大沙漠区衰减补偿逆时偏移成像方法,应用于一个大沙漠区实际数据,取得了理想的计算效果。图8a和8b分别表示计算域的偏移速度场和偏移Q场。图9a表示本发明得到的成像结果,图9b表示常规无衰减补偿的成像结果。结果证明,本发明可以准确地对大沙漠存在不规则起伏衰减沙丘地表采集的数据进行成像。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (4)
1.一种稳定的大沙漠区衰减补偿逆时偏移成像系统,其特征在于:包括输入模块、坐标变换模块、波场正向延拓模块、波场反向延拓模块、稳定化算子模块、坐标反变换模块、互相关成像模块和输出模块;
输入模块,被配置为用于输入物理域大沙漠区偏移速度模型、Q模型和实际观测炮记录,并建立观测系统;
坐标变换模块,被配置为用于将大沙漠区的偏移速度模型和Q模型转换到计算域;
波场正向延拓模块,被配置为用于正向延拓更新波场,得到计算域波场正传波场快照;
波场反向延拓模块,被配置为用于得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,并进行逆时延拓,得到计算域波场反传波场快照;
稳定化算子模块,被配置为用于引入一种计算域稳定化算子压制波场反向延拓的高频噪音;
坐标反变换模块,被配置为用于将正向延拓的波场快照与反传延拓的波场快照反变换到物理域;
互相关成像模块,被配置为用于应用物理域零相位互相关成像条件得到成像结果;
输出模块,被配置为用于输出大沙漠区采集数据的成像结果。
2.一种稳定的大沙漠区衰减补偿逆时偏移成像方法,其特征在于:采用如权利要求1所述的稳定的大沙漠区衰减补偿逆时偏移成像系统,包括如下步骤:
步骤1:通过输入模块,输入物理域大沙漠区偏移速度模型、Q模型和实际观测炮记录,并建立观测系统;
步骤2:通过坐标变换模块,采用式(1)将大沙漠区的偏移速度模型和Q模型转换到计算域;
其中,(x,z)表示物理域的坐标表达式;(ξ,η)表示计算域的坐标表达式;
步骤3:通过波场正向延拓模块,采用波动方程进行正向延拓更新波场,在正向延拓过程中补偿振幅衰减、校正相位频散,得到计算域波场正传波场快照;
步骤4:通过波场反向延拓模块,得到能够准确校正大沙漠区采集数据振幅衰减和相位频散的不规则衰减沙丘地表的反向延拓波动方程,如式(8)所示,并进行逆时延拓,得到计算域波场反传波场快照;
其中,pR为反向延拓的波场快照,dobs为输入的实际观测炮记录;
步骤5:通过稳定化算子模块,引入一种计算域稳定化算子得到波场反向延拓波动方程,来压制波场反向延拓的高频噪音;
步骤6:通过坐标反变换模块,将正向延拓的波场快照与反传延拓的波场快照反变换到物理域;
步骤7:通过互相关成像模块,应用如式(19)所示的物理域零相位互相关成像条件得到成像结果,并通过输出模块,输出大沙漠区采集数据的成像结果:
其中,x=(x,z)表示物理域坐标,T表示记录时间,I表示成像结果。
3.根据权利要求2所述的稳定的大沙漠区衰减补偿逆时偏移成像方法,其特征在于:
在步骤3中,波动方程如式(2)所示:
其中,f表示震源,v0为零频速度,ps为正向延拓的声压场,t为时间;τ为松弛时间,由品质因子Q计算得到:A(ξ)、B(ξ,η)和C(η)分别由式(3)、(4)、(5)求得:
上式所示的波动方程中的第二项为振幅衰减项,第三项为相位频散项;当模拟振幅衰减时,第二项的符号为正号,当进行衰减补偿时,第二项的符号变为负号;在进行相位频散校正时,第三项的符号保持不变;在沙漠区不规则衰减地表条件下,振幅衰减的影响要远大于相位频散的影响,如果只对振幅衰减进行补偿,采用如下所示的波动方程进行波场正向延拓:
其中,v表示偏移速度;
采用如式(7)所示的伪谱法波场递推公式进行波场计算:
4.根据权利要求2所述的稳定的大沙漠区衰减补偿逆时偏移成像方法,其特征在于:
在步骤5中,波场反向延拓波动方程如式(9)所示:
为消除一阶项,定义Λt如下:
引入变量qR,使得qR满足下式:
qR(x,t)=ΛtpR(x,t) (11);
将方程(11)代入到不规则衰减沙丘地表的反向延拓波动方程中,并略去震源项可得:
的一阶偏导数和二阶偏导数为:
将方程(13)和(14)代入到方程(12)中并进行简化可得:
从方程(15)我们可看出该方程不再包含一阶时间偏导数,形成了接近于二阶声波方程的波动方程,则衰减的波场反向延拓的波动方程等价于:
为压制高频噪音,引入一个稳定化项到指数项中,则方程(16)可变为:
其中,σ为规则化参数;
将方程组(17)中的两个方程结合并进行简化,可得引入稳定化算子的波场反向延拓波动方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810459119.3A CN108919356B (zh) | 2018-05-15 | 2018-05-15 | 一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810459119.3A CN108919356B (zh) | 2018-05-15 | 2018-05-15 | 一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919356A CN108919356A (zh) | 2018-11-30 |
CN108919356B true CN108919356B (zh) | 2019-11-05 |
Family
ID=64403869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810459119.3A Active CN108919356B (zh) | 2018-05-15 | 2018-05-15 | 一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919356B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109901223B (zh) * | 2019-03-20 | 2019-11-29 | 吉林大学 | 基于稳相区间理论的逆时偏移成像方法 |
CN110376646B (zh) * | 2019-07-19 | 2020-08-11 | 中国石油大学(华东) | 一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法 |
CN110488354B (zh) * | 2019-07-19 | 2020-08-11 | 中国石油大学(华东) | 一种q补偿的起伏地表棱柱波与一次波联合最小二乘逆时偏移成像方法 |
CN112327358B (zh) * | 2020-10-29 | 2022-05-24 | 中国石油大学(华东) | 一种粘滞性介质中声波地震数据正演模拟方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9829592B2 (en) * | 2014-12-16 | 2017-11-28 | Pgs Geophysical As | Seismic imaging with visco-acoustic reverse-time migration using pseudo-analytical method |
CN105403919B (zh) * | 2015-11-11 | 2018-02-02 | 中国石油天然气集团公司 | 一种逆时偏移成像方法及装置 |
CN106291687A (zh) * | 2016-07-21 | 2017-01-04 | 中国地质科学院地质研究所 | 各向异性多波高斯束叠前深度偏移成像方法 |
CN107367761B (zh) * | 2017-07-19 | 2019-04-26 | 中国石油化工股份有限公司 | 一种大沙漠区深层碳酸盐岩储层广角反射波成像方法 |
-
2018
- 2018-05-15 CN CN201810459119.3A patent/CN108919356B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108919356A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108919356B (zh) | 一种稳定的大沙漠区衰减补偿逆时偏移成像系统及方法 | |
AU2009282330B2 (en) | Estimation of soil properties using waveforms of seismic surface waves | |
CN107505654B (zh) | 基于地震记录积分的全波形反演方法 | |
US8560242B2 (en) | Pseudo-analytical method for the solution of wave equations | |
CN105652321B (zh) | 一种粘声各向异性最小二乘逆时偏移成像方法 | |
CN104122585B (zh) | 基于弹性波场矢量分解与低秩分解的地震正演模拟方法 | |
CN110133723B (zh) | 一种黏声一次波与分阶多次波联合成像方法和系统 | |
CN110579795B (zh) | 基于被动源地震波形及其逆时成像的联合速度反演方法 | |
CN112327358B (zh) | 一种粘滞性介质中声波地震数据正演模拟方法 | |
MX2011003850A (es) | Estimado de señal de dominio de imagen a interferencia. | |
CN111665556B (zh) | 地层声波传播速度模型构建方法 | |
CN105182414B (zh) | 一种基于波动方程正演去除直达波的方法 | |
Zhang et al. | 2D and 3D Q-compensated image-domain least-squares reverse time migration through the hybrid point spread functions and the hybrid deblurring filter | |
CN111025388B (zh) | 一种多波联合的叠前波形反演方法 | |
Gao et al. | Acquisition and processing pitfall with clipped traces in surface-wave analysis | |
CN108680957B (zh) | 基于加权的局部互相关时频域相位反演方法 | |
Bai et al. | Model-based ground-roll attenuation with updating quality factors | |
CN116882214B (zh) | 基于dfl粘弹性方程的瑞雷波数值模拟方法及系统 | |
Li et al. | Two-way wave equation depth migration using one-way propagator extrapolation | |
Colombo et al. | Tackling complex near surface conditions with surface-consistent transmission FWI | |
Zhou et al. | Seismic attenuation compensation using U-Net++ with a pruning scheme discussion | |
Fathalian et al. | Numerical simulation of seismic wave propagation in attenuative transversely isotropic media | |
Xie et al. | Full waveform inversion based on hybrid gradient | |
Bai et al. | Waveform inversion of synthetic reflection data for VTI attenuation parameters | |
Qu et al. | Viscoacoustic Wave Numerical Simulation Method in an Auxiliary System for Irregular Surface Topography based on a Hybrid Spatial Partial Derivative Scheme |
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 |