CN106291691A - 一种地震偏移成像方法及装置 - Google Patents
一种地震偏移成像方法及装置 Download PDFInfo
- Publication number
- CN106291691A CN106291691A CN201610701979.4A CN201610701979A CN106291691A CN 106291691 A CN106291691 A CN 106291691A CN 201610701979 A CN201610701979 A CN 201610701979A CN 106291691 A CN106291691 A CN 106291691A
- Authority
- CN
- China
- Prior art keywords
- data
- seismic
- hourage
- compensate
- compensation
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 94
- 230000005012 migration Effects 0.000 title claims abstract description 71
- 238000013508 migration Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 claims abstract description 54
- 238000001228 spectrum Methods 0.000 description 14
- 230000006872 improvement Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000007704 transition 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
-
- 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
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
技术领域
本申请涉及石油地球物理勘探地震数据处理技术领域,特别涉及一种地震偏移成像方法及装置。
背景技术
地震偏移成像方法是根据在地面上以一定方式进行弹性波激发,并在地面的一定范围(孔径)内记录来自地下弹性分界面的反射波,并利用所述反射波进行成像,进而研究地下地质岩层结构及其物性特征的一种方法。
现有的地震偏移方法通常包括:假设地下介质为声波或者弹性波,获取来自地下的弹性分界面的反射波。利用所述获取到的反射波来进行成像。进一步地,由于地下的地层介质并非完全弹性介质,并且由于地层中存在流体,地震波在实际地下介质传播过程中可能会发生衰减,导致波形振幅能量减弱、相位发生变化,可能会导致地震数据成像的结果不准确,例如:反射轴偏离实际位置、相对振幅关系不准确等,同时成像结果的分辨率也较低。因此,通常还可以利用Q补偿处理方法,例如反Q滤波方法、谱平衡方法等,对获取到的反射波进行处理,使得用于成像的地震数据可以精确反映出地下情况。
发明人发现现有技术中至少存在如下问题:现有的地震偏移方法中,对获取的地震波采用Q补偿处理方法进行处理,由于现有的传统补偿方法都是在时间域进行一维补偿,因此补偿精度较低,从而导致偏移成像的成像结果不准确,效果较差。
发明内容
本申请实施例的目的是提供一种地震偏移成像方法及装置,以提高偏移成像结果的准确性。
为解决上述技术问题,本申请实施例提供一种地震偏移成像方法及装置是这样实现的:
一种地震偏移成像方法,包括:
获取地震数据,根据所述地震数据确定成像速度场数据和品质因子场数据;
根据所述品质因子场数据,计算射线路径补偿旅行时间;
获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据;
根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿;
根据所述补偿后的叠前地震道数据,生成成像结果。
优选方案中,所述根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿,包括:
确定所述走时表和补偿旅行时间表中与所述叠前地震道数据对应的炮点、检波点的走时和补偿旅行时间;
利用所述与所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿。
优选方案中,所述根据补偿后的叠前地震道数据,生成成像结果,包括:将所述补偿后的能量叠加到相应的成像点道集,根据所述成像点道集中的数据,生成成像结果。
优选方案中,在对所述叠前地震道数据进行补偿前,所述方法还包括:对所述叠前地震道数据进行预处理操作;
相应地,所述根据所述走时表和补偿旅行时间表,对叠前地震道数据进行补偿,包括:根据所述走时表和补偿旅行时间表,对经过所述预处理操作后的叠前地震道数据进行补偿。
优选方案中,所述预处理操作包括:傅里叶变换操作。
优选方案中,所述根据品质因子场数据,计算射线路径补偿旅行时间,采用下述公式实现:
其中,T*表示补偿旅行时间,v0表示声波速度,Q表示品质因子,x表示成像点的位置,s表示射线路径。
优选方案中,所述利用与所述地震道数据对应的炮点、检波点的走时和补偿旅行时对所述预处理后的地震数据进行补偿,包括:利用下述公式计算补偿后的成像点的能量:
其中,
σ2=exp[-(0.23G+1.63)]
其中,xr表示检波点位置;xs表示炮点位置;x表示成像点位置;i表示虚数单位;ω表示角频率;ω0表示参考频率,取值为所述地震数据的主频;us(xr,xs,ω)表示变换到频率域的地震道数据;u's(xr,xs,x)表示经过相位、振幅补偿的成像点位置的能量;T*(xr,xs,x)表示炮点和检波点到成像点的补偿旅行时之和;F(xr,xs,x,ω)用于表示补偿滤波器;W(xr,xs,x,ω)表示振幅补偿因子,G表示增益限制。
一种地震偏移成像装置,包括:地震数据获取模块、射线路径补偿旅行时间计算模块、关联数据获取模块、补偿模块和成像模块;其中,
所述地震数据获取模块,用于获取地震数据,根据所述地震数据确定成像速度场数据和品质因子场数据;
所述射线路径补偿旅行时间计算模块,用于根据所述品质因子场数据,计算射线路径补偿旅行时间;
所述关联数据获取模块,用于获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据;
所述补偿模块,用于根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿;
所述成像模块,用于根据所述补偿后的叠前地震道数据,生成成像结果。
优选方案中,所述补偿模块包括:数据确定子模块和补偿计算子模块;其中,
数据确定子模块,用于确定所述走时表和补偿旅行时间表中与所述叠前地震道数据对应的炮点、检波点的走时和补偿旅行时间;
补偿计算子模块,用于利用所述数据确定子模块确定的与所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿。
优选方案中,所述地震偏移成像装置,还包括:预处理模块,用于对所述叠前地震道数据进行预处理操作;
相应地,所述补偿模块,用于根据所述走时表和补偿旅行时间表,对所述预处理模块得到的经过所述预处理操作后的叠前地震道数据进行补偿。
由以上本申请实施例提供的技术方案可见,本申请实施例提供的一种地震偏移成像方法及装置,通过获取的与地震数据对应的走时表和补偿旅行时间表,对所述地震数据对应的叠前地震道数据进行补偿,补偿过程中,可以控制高频成分振幅补偿的稳定性以及补偿频带的宽度,从而可以有效、稳定地补偿地层吸收衰减效应引起的在偏移成像结果上能量及频率成分的衰减。因此本申请实施例提供一种地震偏移成像方法可以提高成像结果的分辨率,从而提高了偏移成像结果的准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请地震偏移成像方法一个实施例的流程图;
图2是本申请方法实施例中不同增益限制参数对应的振幅谱补偿系数示意图;
图3示出了利用现有技术得到的预设数据的地震偏移成像结果;
图4示出了利用本申请方法实施例得到的图3的预设数据的地震偏移成像结果;
图5是图3所示的地震偏移成像结果与图4所示的偏移成像结果的波数谱对比示意图;
图6示出了利用现有技术得到的实际地震数据的地震偏移成像结果;
图7示出了利用本申请方法实施例得到的图6的实际地震数据的地震偏移成像结果;
图8是图6所示的地震偏移成像结果与图7所示的偏移成像结果的波数谱对比示意图;
图9是本申请地震偏移成像装置一个实施例的模块图;
图10是本申请装置实施例中补偿模块的一种组成结构图。
具体实施方式
本申请实施例提供一种地震偏移成像方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请地震偏移成像方法一个实施例的流程图。如图1所示,所示地震偏移方法可以包括:
S101:获取地震数据,根据所述地震数据确定成像速度场数据和品质因子场数据。
可以获取地震数据。所述地震数据可以包括:地震时间信息、位置信息以及振幅信息。
根据所述地震数据,可以确定成像速度场数据和品质因子场数据。所述确定成像速度场数据和品质因子场数据,可以采用本领域中的任意可实现方法,本申请对此并不作出限定。
S102:根据所述品质因子场数据,计算射线路径补偿旅行时间。
根据所述品质因子场数据,可以计算射线路径补偿旅行时间。具体地,可以通过下述公式(1)计算得到:
其中,T*可以表示补偿旅行时,v0可以表示声波速度,Q可以表示品质因子,x可以表示成像点的位置,s可以表示射线路径。所述计算得到的补偿旅行时间可以体现地震波的振幅衰减和相位变化。
S103:获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据。
可以获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据。
S104:根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿。
根据所述走时表和补偿旅行时间表,可以对所述叠前地震道数据进行补偿。具体地,可以包括:确定所述走时表和补偿旅行时间表中与所述叠前地震道数据对应的炮点、检波点的走时和补偿旅行时间,可以利用所述与所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿。
在另一个实施方式中,所述方法还可以包括:对所述叠前地震道数据进行预处理操作。那么,相应地,根据所述走时表和补偿旅行时间表,对叠前地震道数据进行补偿,可以包括:根据所述走时表和补偿旅行时间表,对经过所述预处理操作后的叠前地震道数据进行补偿。所述预处理操作包括:傅里叶变换操作。通过所述傅里叶变换操作,可以将所述叠前地震道数据从时域变换到频域,从而在频域中对所述叠前地震道数据进行补偿。
利用所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿,具体可以采用下述公式(2)实现:
其中,
σ2=exp[-(0.23G+1.63)];
其中,xr可以表示检波点位置。xs可以表示炮点位置。x可以表示成像点位置。i可以表示虚数单位。ω可以表示角频率。ω0可以表示参考频率,其取值可以为所述地震数据的主频。us(xr,xs,ω)可以表示变换到频率域的地震道数据。u's(xr,xs,x)表示经过相位、振幅补偿的成像点位置的能量。T*(xr,xs,x)可以表示炮点和检波点到成像点的补偿旅行时之和。F(xr,xs,x,ω)可以用于表示补偿滤波器。W(xr,xs,x,ω)可以表示振幅补偿因子。G可以表示增益限制参数。所述G可以用于控制高频成分振幅补偿系数的稳定性。图2为本申请方法实施例中不同增益限制参数对应的振幅谱补偿系数示意图,图2中示出了增益限制参数分别为10、20、30时的振幅谱补偿系数的取值。
S105:根据所述补偿后的叠前地震道数据,生成成像结果。
可以根据所述补偿后的叠前地震道数据,生成成像结果。具体地,可以包括:将所述补偿后的能量叠加到相应的成像点道集,根据所述成像点道集中的数据,生成成像结果。
上述实施例提供的地震偏移成像方法,通过获取的与地震数据对应的走时表和补偿旅行时间表,对所述地震数据对应的叠前地震道数据进行补偿,补偿过程中,可以控制高频成分振幅补偿的稳定性以及补偿频带的宽度,从而可以有效、稳定地补偿地层吸收衰减效应引起的在偏移成像结果上能量及频率成分的衰减。因此本申请实施例提供一种地震偏移成像方法可以提高成像结果的分辨率,从而提高了偏移成像结果的准确性。
图3示出了利用现有技术得到的预设数据的地震偏移成像结果。图4示出了利用本申请方法实施例得到的图3的预设数据的地震偏移成像结果。图5是图3所示的地震偏移成像结果与图4所示的偏移成像结果的波数谱对比示意图,图5中QPSDM的曲线表示图4所示的地震偏移成像结果的波数谱,常规PSDM的曲线表示图3所示的偏移成像结果的波数谱。从上述预设数据的例子的现有技术偏移成像结果与本申请方法偏移成像结果及其波数谱对比可以看出,本申请实施例提供的方法相对现有的常规偏移方法,可以有效、稳定地补偿地层吸收衰减效应引起的在偏移成像结果上能量及频率成分的衰减。
图6示出了利用现有技术得到的实际地震数据的地震偏移成像结果。图7示出了利用本申请方法实施例得到的图6的实际地震数据的地震偏移成像结果。图8是图6所示的地震偏移成像结果与图7所示的偏移成像结果的波数谱对比示意图,图8中Q偏移的曲线表示图7所示的偏移成像结果的波数谱,常规偏移的曲线表示图6所示的地震偏移成像结果的波数谱。从上述实际地震数据的现有技术偏移成像结果与本申请方法偏移成像结果及其波数谱对比可以看出,本申请实施例提供的方法的处理效果相比于现有的常规偏移方法的处理效果,提高了成像结果的分辨率。
图9是本申请地震偏移成像装置一个实施例的模块图。如图9所示,所述地震偏移成像装置可以包括:地震数据获取模块901、射线路径补偿旅行时间计算模块902、关联数据获取模块903、补偿模块904和成像模块905。其中,
所述地震数据获取模块901,可以用于获取地震数据,根据所述地震数据确定成像速度场数据和品质因子场数据。
所述射线路径补偿旅行时间计算模块902,可以用于根据所述品质因子场数据,计算射线路径补偿旅行时间。
所述关联数据获取模块903,可以用于获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据。
所述补偿模块904,可以用于根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿。
所述成像模块905,可以用于根据所述补偿后的叠前地震道数据生成成像结果。
图10是本申请装置实施例中补偿模块的一种组成结构图。如图10所示,所述补偿模块904可以包括:数据确定子模块9041和补偿计算子模块9042。其中,
所述数据确定子模块9041,可以用于确定所述走时表和补偿旅行时间表中与所述叠前地震道数据对应的炮点、检波点的走时和补偿旅行时间。
所述补偿计算子模块9042,可以用于利用所述数据确定子模块9041确定的与所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿。
在另一个实施方式中,所述地震偏移成像装置还可以包括:预处理模块,可以用于对所述叠前地震道数据进行预处理操作。那么相应地,所述补偿模块904,可以用于根据所述走时表和补偿旅行时间表,对所述预处理模块得到的经过所述预处理操作后的叠前地震道数据进行补偿。
上述实施例公开的地震偏移成像装置与本申请的地震偏移成像方法实施例相对应,可以实现本申请方法实施例并获得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (10)
1.一种地震偏移成像方法,其特征在于,包括:
获取地震数据,根据所述地震数据确定成像速度场数据和品质因子场数据;
根据所述品质因子场数据,计算射线路径补偿旅行时间;
获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据;
根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿;
根据所述补偿后的叠前地震道数据,生成成像结果。
2.根据权利要求1所述的一种地震偏移成像方法,其特征在于,所述根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿,包括:
确定所述走时表和补偿旅行时间表中与所述叠前地震道数据对应的炮点、检波点的走时和补偿旅行时间;
利用所述与所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿。
3.根据权利要求1所述的一种地震偏移成像方法,其特征在于,所述根据补偿后的叠前地震道数据,生成成像结果,包括:将所述补偿后的能量叠加到相应的成像点道集,根据所述成像点道集中的数据,生成成像结果。
4.根据权利要求1所述的一种地震偏移成像方法,其特征在于,在对所述叠前地震道数据进行补偿前,所述方法还包括:对所述叠前地震道数据进行预处理操作;
相应地,所述根据所述走时表和补偿旅行时间表,对叠前地震道数据进行补偿,包括:根据所述走时表和补偿旅行时间表,对经过所述预处理操作后的叠前地震道数据进行补偿。
5.根据权利要求1所述的一种地震偏移成像方法,其特征在于,所述预处理操作包括:傅里叶变换操作。
6.根据权利要求1所述的一种地震偏移成像方法,其特征在于,所述根据品质因子场数据,计算射线路径补偿旅行时间,采用下述公式实现:
其中,T*表示补偿旅行时间,v0表示声波速度,Q表示品质因子,x表示成像点的位置,s表示射线路径。
7.根据权利要求2所述的一种地震偏移成像方法,其特征在于,所述利用与所述地震道数据对应的炮点、检波点的走时和补偿旅行时对所述预处理后的地震数据进行补偿,包括:利用下述公式计算补偿后的成像点的能量:
其中,
σ2=exp[-(0.23G+1.63)]
其中,xr表示检波点位置;xs表示炮点位置;x表示成像点位置;i表示虚数单位;ω表示角频率;ω0表示参考频率,取值为所述地震数据的主频;us(xr,xs,ω)表示变换到频率域的地震道数据;u's(xr,xs,x)表示经过相位、振幅补偿的成像点位置的能量;T*(xr,xs,x)表示炮点和检波点到成像点的补偿旅行时之和;F(xr,xs,x,ω)用于表示补偿滤波器;W(xr,xs,x,ω)表示振幅补偿因子,G表示增益限制。
8.一种地震偏移成像装置,其特征在于,包括:地震数据获取模块、射线路径补偿旅行时间计算模块、关联数据获取模块、补偿模块和成像模块;其中,
所述地震数据获取模块,用于获取地震数据,根据所述地震数据确定成像速度场数据和品质因子场数据;
所述射线路径补偿旅行时间计算模块,用于根据所述品质因子场数据,计算射线路径补偿旅行时间;
所述关联数据获取模块,用于获取所述地震数据对应的走时表、补偿旅行时间表和叠前地震道数据;
所述补偿模块,用于根据所述走时表和补偿旅行时间表,对所述叠前地震道数据进行补偿;
所述成像模块,用于根据所述补偿后的叠前地震道数据,生成成像结果。
9.根据权利要求8所述的一种地震偏移成像装置,其特征在于,所述补偿模块包括:数据确定子模块和补偿计算子模块;其中,
数据确定子模块,用于确定所述走时表和补偿旅行时间表中与所述叠前地震道数据对应的炮点、检波点的走时和补偿旅行时间;
补偿计算子模块,用于利用所述数据确定子模块确定的与所述地震道数据对应的炮点、检波点的走时和补偿旅行时间对所述叠前地震数据进行补偿。
10.根据权利要求8所述的一种地震偏移成像装置,其特征在于,还包括:预处理模块,用于对所述叠前地震道数据进行预处理操作;
相应地,所述补偿模块,用于根据所述走时表和补偿旅行时间表,对所述预处理模块得到的经过所述预处理操作后的叠前地震道数据进行补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610701979.4A CN106291691A (zh) | 2016-08-22 | 2016-08-22 | 一种地震偏移成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610701979.4A CN106291691A (zh) | 2016-08-22 | 2016-08-22 | 一种地震偏移成像方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106291691A true CN106291691A (zh) | 2017-01-04 |
Family
ID=57662077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610701979.4A Pending CN106291691A (zh) | 2016-08-22 | 2016-08-22 | 一种地震偏移成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291691A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181649A (zh) * | 2017-11-23 | 2018-06-19 | 中国石油天然气集团公司 | 一种介质分界面振幅补偿方法及装置 |
CN109085648A (zh) * | 2018-07-16 | 2018-12-25 | 中国科学院地质与地球物理研究所 | 叠前深度偏移方法和装置 |
CN110031898A (zh) * | 2019-05-14 | 2019-07-19 | 中国科学院地质与地球物理研究所 | 数据优化方法及积分法叠前深度偏移方法 |
CN113009569A (zh) * | 2019-12-20 | 2021-06-22 | 中国石油天然气集团有限公司 | 一种地震偏移成像方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030208321A1 (en) * | 2001-07-31 | 2003-11-06 | Ruben Martinez | Relative true amplitude migration |
CN101661112A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油天然气集团公司 | 一种准确的地震构造成像方法 |
CN102109612A (zh) * | 2009-12-23 | 2011-06-29 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
CN103675898A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种快速精确估算地层吸收所致时差的方法 |
-
2016
- 2016-08-22 CN CN201610701979.4A patent/CN106291691A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030208321A1 (en) * | 2001-07-31 | 2003-11-06 | Ruben Martinez | Relative true amplitude migration |
CN101661112A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油天然气集团公司 | 一种准确的地震构造成像方法 |
CN102109612A (zh) * | 2009-12-23 | 2011-06-29 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
CN103675898A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种快速精确估算地层吸收所致时差的方法 |
Non-Patent Citations (4)
Title |
---|
MAUD CAVALCA 等: "Ray‐based tomography for Q estimation and Q compensation in complex media", 《SEG TECHNICAL PROGRAM EXPANDED ABSTRACTS 2011》 * |
YANGHUA WANG: "Inverse Q-filter for seismic resolution enhancement", 《GEOPHYSICS》 * |
孙赞东 等: "《三维三分量VSP方法原理及应用》", 31 December 2011, 石油工业出版社 * |
杨仁虎 等: "叠前逆时偏移影响因素分析", 《地球物理学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181649A (zh) * | 2017-11-23 | 2018-06-19 | 中国石油天然气集团公司 | 一种介质分界面振幅补偿方法及装置 |
CN108181649B (zh) * | 2017-11-23 | 2019-09-10 | 中国石油天然气集团公司 | 一种介质分界面振幅补偿方法及装置 |
CN109085648A (zh) * | 2018-07-16 | 2018-12-25 | 中国科学院地质与地球物理研究所 | 叠前深度偏移方法和装置 |
CN109085648B (zh) * | 2018-07-16 | 2019-09-27 | 中国科学院地质与地球物理研究所 | 叠前深度偏移方法和装置 |
CN110031898A (zh) * | 2019-05-14 | 2019-07-19 | 中国科学院地质与地球物理研究所 | 数据优化方法及积分法叠前深度偏移方法 |
CN110031898B (zh) * | 2019-05-14 | 2019-12-13 | 中国科学院地质与地球物理研究所 | 数据优化方法及积分法叠前深度偏移方法 |
CN113009569A (zh) * | 2019-12-20 | 2021-06-22 | 中国石油天然气集团有限公司 | 一种地震偏移成像方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842304B (zh) | 一种叠前深度偏移方法及装置 | |
Decker et al. | Diffraction imaging and time-migration velocity analysis using oriented velocity continuation | |
CN106291691A (zh) | 一种地震偏移成像方法及装置 | |
CN106646611B (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN107193045B (zh) | 一种地震数据处理方法及装置 | |
Sun et al. | Automatic wave-equation migration velocity analysis by focusing subsurface virtual sources | |
CN109061734A (zh) | 一种地震数据的处理方法及装置 | |
CN106054251A (zh) | 一种初至波拾取方法及装置 | |
CN108828668A (zh) | 一种叠前时间偏移数据处理方法及装置 | |
Ma et al. | Wavefield extrapolation in pseudodepth domain | |
CN106353818A (zh) | 一种地震数据时差校正的方法及系统 | |
CN107179546B (zh) | 一种确定速度场的方法及装置 | |
Lameloise et al. | Extension of migration velocity analysis to transmitted wavefields | |
Farshad et al. | From constant-to variable-density inverse extended Born modeling | |
Gadylshin et al. | Machine learning-based numerical dispersion mitigation in seismic modelling | |
Wang et al. | An explicit stabilization scheme for Q-compensated reverse time migration | |
CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
CN106353816A (zh) | 一种地震采集脚印噪音压制方法及系统 | |
Khoshnavaz et al. | Oriented prestack time migration using local slopes and predictive painting in the common-source domain for planar reflectors | |
CN105353406A (zh) | 一种生成角道集的方法和装置 | |
Faucher et al. | A priori estimates of attraction basins for velocity model reconstruction by time-harmonic full-waveform inversion and data-space reflectivity formulation | |
Mittet | Small-scale medium variations with high-order finite-difference and pseudospectral schemes | |
CN108693557A (zh) | 一种地震数据处理方法及装置 | |
Weibull et al. | Reverse-time demigration using the extended-imaging condition | |
CN106772599A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170104 |