CN112946744A - 一种基于动态时差规整的最小二乘偏移成像方法及系统 - Google Patents
一种基于动态时差规整的最小二乘偏移成像方法及系统 Download PDFInfo
- Publication number
- CN112946744A CN112946744A CN201911281824.XA CN201911281824A CN112946744A CN 112946744 A CN112946744 A CN 112946744A CN 201911281824 A CN201911281824 A CN 201911281824A CN 112946744 A CN112946744 A CN 112946744A
- Authority
- CN
- China
- Prior art keywords
- seismic
- observation data
- migration
- work area
- target work
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 80
- 238000013508 migration Methods 0.000 claims abstract description 75
- 230000005012 migration Effects 0.000 claims abstract description 61
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000012937 correction Methods 0.000 claims abstract description 29
- 230000010363 phase shift Effects 0.000 claims abstract description 20
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000010606 normalization Methods 0.000 claims abstract description 3
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 9
- 230000006835 compression Effects 0.000 claims description 7
- 238000007906 compression Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 8
- 238000004088 simulation Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004904 shortening Methods 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种基于动态时差规整的最小二乘偏移成像方法及系统,该方法包括:利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果;对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间‑空间域合成单炮地震记录;逐道对所述时间‑空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录;求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量;利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
Description
技术领域
本发明涉及石油勘探技术领域,尤其涉及一种基于动态时差规整的最小二乘偏移成像方法及系统。
背景技术
地震数据处理的重要目的之一是通过深度偏移(Depth Migration)得到高质量的深度域图像。深度域图像的质量依赖于原始地震单炮数据的质量,提高地震数据的信噪比在整个地震数据处理中起着非常重要的作用。现有技术在地震数据深度域偏移成像过程中,往往利用最小二乘偏移成像。最小二乘偏移成像的基本原理是通过比较反偏移得到的模拟数据和实际观测数据的差别来确定偏移成像结果的改变量,反复迭代此过程直到模拟数据和观测数据的差别达到最小。在现有最小二乘法成像过程中,存在因速度模型不准确,而导致迭代过程中反偏移得到的模拟数据和实际观测数据存在周波跳跃现象,因此无法保证最小二乘法偏移成像结果的精确度。
发明内容
为解决因速度模型不准确而导致迭代过程中反偏移得到的模拟数据和实际观测数据存在周波跳跃的问题,本发明实施例提供了一种基于动态时差规整的最小二乘偏移成像方法及系统,从而保证最小二乘偏移成像结果的精确度。
一方面,本发明实施例提供了一种基于动态时差规整的最小二乘偏移成像方法,所述方法包括:
利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果;
对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间-空间域合成单炮地震记录;
逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录;
求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量;
利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
在一实施例中,逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录,包括:
通过压缩或者延展的方式,将所述地震观测数据在时间序列上与时间-空间域合成单炮地震记录进行对齐,使得地震观测数据向所述时间-空间域合成单炮地震记录匹配,得到匹配后的模拟地震记录。
在一实施例中,所述修正量为对所述残差进行偏移得到的梯度场。
另一方面,本发明实施例还提供了一种基于动态时差规整的最小二乘偏移成像系统,所述系统包括:
偏移成像单元,用于利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果;
反偏移单元,用于对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间-空间域合成单炮地震记录;
动态时差规整单元,用于逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录;
修正量获取单元,用于求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量;
成像结果修正单元,用于利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
在一实施例中,所述动态时差规整单元具体用于:
通过压缩或者延展的方式,将所述地震观测数据在时间序列上与时间-空间域合成单炮地震记录进行对齐,使得地震观测数据向所述时间-空间域合成单炮地震记录匹配,得到匹配后的模拟地震记录。
在一实施例中,所述修正量为对所述残差进行偏移得到的梯度场。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述的基于动态时差规整的最小二乘偏移成像方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行所述的基于动态时差规整的最小二乘偏移成像方法的计算机程序。
和传统的最小二乘算法相比,本发明实施例提供的基于动态时差规整的最小二乘偏移成像方法及系统大大改进了深度域偏移成像的精确度。解决了现有最小二乘偏移方法依赖于精确的速度模型的问题,消除了反演过程中反偏移得到的模拟波场和实际观测波场存在周波跳跃的问题,确保最小二乘偏移能够收敛到可靠的结果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于动态时差规整的最小二乘偏移成像方法的流程示意图;
图2为本发明实施例提供的一用于偏移的单炮地震记录;
图3是对图2所示地震记录应用反偏移方法得到的单炮地震记录;
图4是对图2所示地震记录应用动态时差规整方法得到的单炮地震记录;
图5为本发明实施例提供的基于动态时差规整的最小二乘偏移成像系统的结构示意图;
图6为本发明另一实施例提供的基于动态时差规整的最小二乘偏移成像的计算机设备结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的基于动态时差规整的最小二乘偏移成像方法的流程示意图。如图1所示,该方法主要包括以下步骤:
步骤S1、利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果。
步骤S2、对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间-空间域合成单炮地震记录。
步骤S3、逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整(Dynamic Time Warping),得到匹配后的模拟地震记录。
步骤S4、求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量。
步骤S5、利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
动态时差规整最早用于处理语音方面识别分类的问题,它的原理是通过将时间序列进行延伸和缩短,来计算两个时间序列之间的相似性。最小二乘偏移法的基本原理是通过比较反偏移得到的模拟数据和实际观测数据的差别来确定偏移成像结果的改变量,反复迭代此过程直到模拟数据和观测数据的差别达到最小。本发明实施例利用动态时差规整方法使得反偏移得到的地震记录和地震观测数据误差逐渐减小,通过设立目标函数,多次迭代步骤S2~步骤S5,使目标函数最小化,从而改进目标工区偏移剖面的精确度。和传统的最小二乘算法相比,利用本发明实施例计算得到模拟的单炮地震数据信噪比得到提高,基于动态时差规整的最小二乘偏移成像大大改进了深度域偏移成像的精确度。
通常地,在目标工区沿炮线方向人工激发地震波后,由检波器记录地震波波场,即地震观测数据。在步骤S1中,可将地震波观测数据及速度模型作为输入,通过运行基于相移加插值方法的单程波叠前深度偏移模块(Oneway PSDM)即可得到初始偏移成像结果。
在一实施例中,利用步骤S3逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,即通过压缩或者延展的方式,将地震观测数据在时间序列上与时间-空间域合成单炮地震记录进行对齐,使得地震观测数据向所述时间-空间域合成单炮地震记录匹配,从而得到匹配后的模拟地震记录。
在一实施例中,步骤S4中的修正量即为对所述残差进行偏移得到的梯度场,即在步骤S5中利用梯度场对初始偏移成像结果进行修正,从而得到目标工区偏移剖面的成像结果。
为了更清楚的对本发明提出的成像方法进行说明,本发明还提供了一具体实施例。首先,选择任意一工区,沿炮线方向人工激发地震波,由检波器接收并记录地震波场,得到地震观测数据ds,r=d(xs,xr,t),其中xs和xr分别代表震源和接收点的空间位置,t代表观测时间。图2是该实施例记录到的用于偏移的原始单炮地震数据记录图。
其次,将地震观测数据ds,r=d(xs,xr,t)与速度模型V作为输入数据,通过运行基于相移加插值方法的单程波叠前深度偏移模块得到偏移成像结果。
具体地,在时间-波数域求出声波方程的解析解,即由波动方程:
由此计算下行波的相移为exp(-ikzΔz),上行波的相移为exp(ikzΔz)。
将地震观测数据与速度模型作为输入数据,应用下行波相移因子exp(-ikzΔz),根据互相关成像条件:I(x,zj)=PS(x,zj,ωl)PR(x,zj,ωl)j=1,2,…,N,得到单程波叠前深度偏移成像结果。然后,通过基于相移加插值进行波场延拓的反偏移算法得到用于偏移的时间-空间域合成单炮地震记录:图3是对图2所示地震记录应用反偏移方法得到的单炮地震记录。
再次,对反偏移得到的模拟地震记录和地震观测数据ds,r=d(xs,xr,t)逐道进行动态时差规整,即将地震观测数据在时间序列上通过压缩或者延展的方法进行对齐,使得地震观测数据ds,r=d(xs,xr,t)向反偏移得到的模拟地震记录匹配,得到一个匹配后的新的模拟地震记录图4是对图2所示地震记录应用动态时差规整方法得到的单炮地震记录。
最后,对原始单炮地震记录数据ds,r=d(xs,xr,t)和动态时差规整算法得到的匹配后的模拟地震记录数据求取残差,对地震记录数据残差再进行偏移得到偏移剖面的修正量,对初始偏移成像结果进行修改。利用动态时差规整方法使得反偏移得到的模拟波场和野外记录的地震波场误差逐渐减小,通过设立目标函数,多次迭代上述过程,使目标函数最小化,从而改进偏移剖面的精确度。
利用本发明实施例提供的高效率、高精度的动态时差规整方法,可在最小二乘偏移技术的实现过程中,解决现有最小二乘偏移方法依赖于精确的速度模型的问题,消除反演过程中反偏移模拟波场和实际观测波场存在周波跳跃的问题,确保最小二乘偏移能够收敛到可靠的结果。
基于与图1所示基于动态时差规整的最小二乘偏移成像方法相同的发明构思,本发明实施例还提供了一种基于动态时差规整的最小二乘偏移成像系统,如下面实施例所述。由于该系统的实施与图1中基于动态时差规整的最小二乘偏移成像方法相似,因此该系统的实施可以参见图1中的基于动态时差规整的最小二乘偏移成像方法的实施,重复之处不再赘述。
在另一实施例中,本发明还提供了一种基于动态时差规整的最小二乘偏移成像系统,其结构如图5所示,该系统主要包括:偏移成像单元10、反偏移单元20、动态时差规整单元30、修正量获取单元40及成像结果修正单元50。
其中,偏移成像单元10用于利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果;反偏移单元20用于对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间-空间域合成单炮地震记录;动态时差规整单元30用于逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录;修正量获取单元40用于求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量;成像结果修正单元50用于利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
在一实施例中,动态时差规整单元30具体用于:通过压缩或者延展的方式,将所述地震观测数据在时间序列上与时间-空间域合成单炮地震记录进行对齐,使得地震观测数据向所述时间-空间域合成单炮地震记录匹配,得到匹配后的模拟地震记录。
在一实施例中,所述修正量为对所述残差进行偏移得到的梯度场。
图6为本发明另一实施例提供的基于动态时差规整的最小二乘偏移成像的计算机设备结构示意图。如图6所示,该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序。所述计算机程序在被所述处理器运行时执行图1所示的步骤。
另外,本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行图1所示步骤的计算机程序。
本发明实施例使最小二乘偏移流程不依赖于精确的速度模型,解决了反演过程中的反偏移模拟数据和实际观测数据存在周波跳跃的问题;并且,研究表明,针对实际资料速度模型不准确的情况下,应用新发明动态时差规整的最小二乘偏移技术后,大大提高了深度域偏移成像的精确度。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于动态时差规整的最小二乘偏移成像方法,其特征在于,所述方法包括:
利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果;
对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间-空间域合成单炮地震记录;
逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录;
求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量;
利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
2.根据权利要求1所述的基于动态时差规整的最小二乘偏移成像方法,其特征在于,逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录,包括:
通过压缩或者延展的方式,将所述地震观测数据在时间序列上与时间-空间域合成单炮地震记录进行对齐,使得地震观测数据向所述时间-空间域合成单炮地震记录匹配,得到匹配后的模拟地震记录。
3.根据权利要求1所述的基于动态时差规整的最小二乘偏移成像方法,其特征在于,所述修正量为对所述残差进行偏移得到的梯度场。
4.一种基于动态时差规整的最小二乘偏移成像系统,其特征在于,所述系统包括:
偏移成像单元,用于利用相移加插值法对目标工区的地震观测数据与速度模型进行单程波叠前深度偏移处理,获取目标工区的初始偏移成像结果;
反偏移单元,用于对所述地震观测数据、初始偏移成像结果进行基于相移加插值法的波场延拓处理实现反偏移,得到目标工区的时间-空间域合成单炮地震记录;
动态时差规整单元,用于逐道对所述时间-空间域合成单炮地震记录及地震观测数据进行动态时差规整,得到匹配后的模拟地震记录;
修正量获取单元,用于求取所述地震观测数据及模拟地震记录的残差,并对残差进行偏移以获取目标工区偏移剖面的修正量;
成像结果修正单元,用于利用所述修正量对所述初始偏移成像结果进行修正,得到目标工区偏移剖面的成像结果。
5.根据权利要求4所述的基于动态时差规整的最小二乘偏移成像系统,其特征在于,所述动态时差规整单元具体用于:
通过压缩或者延展的方式,将所述地震观测数据在时间序列上与时间-空间域合成单炮地震记录进行对齐,使得地震观测数据向所述时间-空间域合成单炮地震记录匹配,得到匹配后的模拟地震记录。
6.根据权利要求4所述的基于动态时差规整的最小二乘偏移成像系统,其特征在于,所述修正量为对所述残差进行偏移得到的梯度场。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至3任一项权利要求所述的方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至3任一项权利要求所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911281824.XA CN112946744B (zh) | 2019-12-11 | 2019-12-11 | 一种基于动态时差规整的最小二乘偏移成像方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911281824.XA CN112946744B (zh) | 2019-12-11 | 2019-12-11 | 一种基于动态时差规整的最小二乘偏移成像方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112946744A true CN112946744A (zh) | 2021-06-11 |
CN112946744B CN112946744B (zh) | 2023-06-30 |
Family
ID=76234406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911281824.XA Active CN112946744B (zh) | 2019-12-11 | 2019-12-11 | 一种基于动态时差规整的最小二乘偏移成像方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112946744B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060062083A1 (en) * | 2004-09-23 | 2006-03-23 | Shu-Schung Lee | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration |
CN102116869A (zh) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | 高精度叠前域最小二乘偏移地震成像技术 |
EP3163328A1 (en) * | 2015-11-02 | 2017-05-03 | CGG Services SA | Seismic data least-square migration method and device |
CN107918155A (zh) * | 2016-10-10 | 2018-04-17 | 中国石油化工股份有限公司 | 反偏移模拟数据时差校正方法及系统 |
CN108241173A (zh) * | 2017-12-28 | 2018-07-03 | 中国石油大学(华东) | 一种地震资料偏移成像方法及系统 |
CN108445532A (zh) * | 2018-02-12 | 2018-08-24 | 中国石油天然气集团有限公司 | 一种深度域反偏移方法及装置 |
CN108802813A (zh) * | 2018-06-13 | 2018-11-13 | 中国石油大学(华东) | 一种多分量地震资料偏移成像方法及系统 |
CN108845355A (zh) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | 地震偏移成像方法及装置 |
-
2019
- 2019-12-11 CN CN201911281824.XA patent/CN112946744B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060062083A1 (en) * | 2004-09-23 | 2006-03-23 | Shu-Schung Lee | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration |
CN102116869A (zh) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | 高精度叠前域最小二乘偏移地震成像技术 |
EP3163328A1 (en) * | 2015-11-02 | 2017-05-03 | CGG Services SA | Seismic data least-square migration method and device |
CN107918155A (zh) * | 2016-10-10 | 2018-04-17 | 中国石油化工股份有限公司 | 反偏移模拟数据时差校正方法及系统 |
CN108241173A (zh) * | 2017-12-28 | 2018-07-03 | 中国石油大学(华东) | 一种地震资料偏移成像方法及系统 |
CN108445532A (zh) * | 2018-02-12 | 2018-08-24 | 中国石油天然气集团有限公司 | 一种深度域反偏移方法及装置 |
CN108802813A (zh) * | 2018-06-13 | 2018-11-13 | 中国石油大学(华东) | 一种多分量地震资料偏移成像方法及系统 |
CN108845355A (zh) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | 地震偏移成像方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112946744B (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9495487B2 (en) | Iterative inversion of data from simultaneous geophysical sources | |
US10557956B2 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
US20190353814A1 (en) | Direct Migration of Simultaneous-Source Survey Data | |
CN111596366B (zh) | 一种基于地震信号优化处理的波阻抗反演方法 | |
CN109917454B (zh) | 基于双基准面的真地表叠前深度偏移成像方法及装置 | |
EP2260331A1 (en) | An efficient method for inversion of geophysical data | |
US20160061974A1 (en) | Full-Wavefield Inversion Using Mirror Source-Receiver Geometry | |
CN110531410B (zh) | 一种基于直达波场的最小二乘逆时偏移梯度预条件方法 | |
CN109239781B (zh) | 一种地震数据校正方法及装置 | |
CN112698390B (zh) | 叠前地震反演方法及装置 | |
CN109696705B (zh) | 深海obs上下行波场成像数据校正方法及装置 | |
AU2015256626B2 (en) | Efficient line search methods for multi-parameter full wavefield inversion | |
CN106353798A (zh) | 多分量联合高斯束叠前逆时偏移成像方法 | |
CN108957553A (zh) | 动校正量递推修正的无拉伸畸变动校正方法及装置 | |
CN111435174B (zh) | 强反射地区地震资料振幅补偿方法及装置 | |
CN112946744A (zh) | 一种基于动态时差规整的最小二乘偏移成像方法及系统 | |
CN111856557A (zh) | 深度域合成地震记录制作方法及装置 | |
CN111208568B (zh) | 一种时间域多尺度全波形反演方法及系统 | |
CN112540408B (zh) | 一种基于深度学习的地震数据静校正处理方法及系统 | |
Zhang et al. | WARpick: A 3D warping tool for seismic gather flattening | |
CN116660981B (zh) | 基于包络的各向异性参数反演方法、装置及存储介质 | |
CN115542387A (zh) | 一种叠前偏移产出张角偏移道集的方法及系统 | |
Ren et al. | Super-trace residual statics estimation technology | |
Espinoza et al. | Short turnaround reprocessing strategy for subsalt prospects |
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 |