CN113219525B - 偏移成像去模糊化方法及装置 - Google Patents
偏移成像去模糊化方法及装置 Download PDFInfo
- Publication number
- CN113219525B CN113219525B CN202010081486.1A CN202010081486A CN113219525B CN 113219525 B CN113219525 B CN 113219525B CN 202010081486 A CN202010081486 A CN 202010081486A CN 113219525 B CN113219525 B CN 113219525B
- Authority
- CN
- China
- Prior art keywords
- vector
- wave number
- point
- target scattering
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 209
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000013598 vector Substances 0.000 claims abstract description 255
- 230000006870 function Effects 0.000 claims description 136
- 239000006185 dispersion Substances 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 3
- 230000005012 migration Effects 0.000 abstract description 15
- 238000013508 migration Methods 0.000 abstract description 15
- 238000010586 diagram Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling 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
-
- 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
-
- 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
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
技术领域
本发明涉及石油地球物理勘探技术领域，尤其涉及偏移成像去模糊化方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
地震偏移成像技术作为勘探地震学的关键技术环节之一，在地震勘探理论研究和实际应用中具有重要的地位，在日趋复杂和精细化的地下介质成像中发挥着越来越重要的作用。传统的偏移成像方法只是Born正演算子的共轭转置，对于有限观测***所采集的带限地震数据，其只能产生模糊的构造成像结果。此外，复杂的地下介质构造和地震数据空间采样的非规则性，还会导致偏移假象和非规则的成像照明，严重影响成像振幅的可靠性。基于线性化地震反演理论的最小二乘偏移技术是解决上述问题的有效途径。该技术利用最优化方法来求解最小二乘偏移所对应的线性反问题，理论上能够消除采集照明等因素对成像结果的影响。
然而，经典的数据域最小二乘偏移成像方法每次迭代需要一次常规的正演和偏移运算，完整的反演过程往往需要几十倍于常规偏移的计算量，因此庞大的计算量严重制约了该偏移成像方法的数据处理效率。成像域最小二乘偏移成像方法则需要计算和存储庞大的Hessian矩阵，即便可以通过数据编码等技术来优化相应的计算过程，其计算和存储的开销依然令人难以接受，同样会使得偏移成像的数据处理效率较低。
点扩散函数是局部化的Hessian矩阵，其描述了地震成像***对于单位点散射率的滤波和模糊效应，反映了成像***对某个成像点的照明能力。在点扩散函数已知的情况下，可以利用点扩散函数对传统的地震成像结果进行去模糊化，从而提高成像结果的分辨率和振幅保真性。然而，在利用点扩散函数进行偏移成像去模糊化时，计算点扩散函数的过程需要进行正演和偏移，计算量同样很大，导致地震偏移成像去模糊化的数据处理效率依然很低。
因此，现有的地震偏移成像去模糊化存在数据处理效率低的问题。
发明内容
本发明实施例提供一种偏移成像去模糊化方法，用以提高地震偏移成像去模糊化的数据处理效率，该偏移成像去模糊化方法包括：
在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量；
根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数；
根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数；
利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，形成地震数据的最终成像结果；
基于下述公式确定每道地震道的波数域点扩散函数：
其中，
表示每道地震道的波数域点扩散函数，F(ω)表示从地震数据中提取的震源子波，ω表示地震波频率，
表示波数域点扩散函数的波数矢量，
表示目标散射点处的成像波数矢量。
本发明实施例还提供一种偏移成像去模糊化装置，用以提高地震偏移成像去模糊化的数据处理效率，该偏移成像去模糊化装置包括：
射线追踪模块，用于在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量；
波数域确定模块，用于根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数；
空间域确定模块，用于根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数；
去模糊化模块，用于利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，形成地震数据的最终成像结果；
基于下述公式确定每道地震道的波数域点扩散函数：
其中，
表示每道地震道的波数域点扩散函数，F(ω)表示从地震数据中提取的震源子波，ω表示地震波频率，
表示波数域点扩散函数的波数矢量，
表示目标散射点处的成像波数矢量。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述偏移成像去模糊化方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述偏移成像去模糊化方法的计算机程序。
本发明实施例中，基于射线追踪确定目标散射点处的成像波数矢量，进而根据目标散射点处的成像波数矢量及震源子波，即可得到地震数据的空间域点扩散函数，最后利用得到的地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化。本发明实施例无需进行正演和偏移，只需基于目标散射点处的成像波数矢量及震源子波，即可快速确定地震数据的空间域点扩散函数，进而利用得到的空间域点扩散函数对成像结果进行去模糊化，能够提高地震偏移成像去模糊化的效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例提供的偏移成像去模糊化方法的实现流程图；
图2为本发明实施例提供的偏移成像去模糊化方法中步骤101的实现流程图；
图3-1为本发明实施例提供的偏移成像去模糊化方法中地震成像射线对示意图；
图3为本发明实施例提供的偏移成像去模糊化方法中步骤101的另一实现流程图；
图4为本发明实施例提供的偏移成像去模糊化方法中步骤102的实现流程图；
图5为本发明实施例提供的偏移成像去模糊化方法中步骤103的实现流程图；
图6为本发明实施例提供的偏移成像去模糊化装置的功能模块图；
图7为本发明实施例提供的偏移成像去模糊化装置中射线追踪模块601的结构框图；
图8为本发明实施例提供的偏移成像去模糊化装置中射线追踪模块601的另一结构框图；
图9为本发明实施例提供的偏移成像去模糊化装置中波数域确定模块602的结构框图；
图10为本发明实施例提供的偏移成像去模糊化装置中空间域确定模块603的结构框图；
图11为本发明实施例提供的利用本发明的偏移成像去模糊化方法得到的某地震数据的波数域点扩散函数；
图12为本发明实施例提供的利用本发明的偏移成像去模糊化方法得到的某地震数据的空间域点扩散函数。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1示出了本发明实施例提供的偏移成像去模糊化方法的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
如图1所示，偏移成像去模糊化方法，其包括：
步骤101，在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量；
步骤102，根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数；
步骤103，根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数；
步骤104，利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，形成地震数据的最终成像结果。
在本发明实施例中，地震数据包括多道地震道，即多道地震数据。在进行偏移成像时，将选定的散射点作为目标散射点，在确定每道地震道目标散射点处的成像波数矢量时，在每道地震道所对应的震源及接收点处进行射线追踪，从而确定目标散射点处的成像波数矢量。
另外，可以从地震数据中提取地震子波，将提取的地震子波作为初始震源子波，进而利用正则反演优化更新初始震源子波，从而得到更新后的震源子波，将更新后的震源子波作为本发明实施例中的震源子波。本领域技术人员可以理解的是，还可以通过其它方法从地震数据中提取震源子波，本发明实施例对此不作特别的限制。
在得到目标散射点对应的成像波数矢量，以及从地震数据中提取到震源子波后，可以基于目标散射点对应的成像波数矢量，及从地震数据中提取到的震源子波，得到每道地震道的波数域点扩散函数，进而再将每道地震道的波数域点扩散函数进行累加，以此形成地震数据的波数域点扩散函数。
在确定地震数据的波数域点扩散函数后，对地震数据的波数域点扩散函数进行转换，以得到地震数据的空间域点扩散函数。其中，将地震数据的波数域点扩散函数转换为地震数据的空间域点扩散函数，为本领域技术人员熟知的现有技术，本领域技术人员可以采用现有技术中的任一方法实现，对此本发明实施例不作详细赘述。
在得到地震数据的空间域点扩散函数，即可利用地震数据的空间域点扩散函数，对地震数据的初步成像结果进行去模糊化。其中，本领域技术人员可以采用现有技术中的任一地震数据偏移成像方法，以得到地震数据的初步成像结果。利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，能够提高地震偏移成像的效果及分辨率。
在本发明实施例中，基于射线追踪确定目标散射点处的成像波数矢量，进而根据目标散射点处的成像波数矢量及震源子波，即可得到地震数据的空间域点扩散函数，最后利用得到的地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化。本发明实施例无需进行正演和偏移，只需基于目标散射点处的成像波数矢量及震源子波，即可快速确定地震数据的空间域点扩散函数，进而利用得到的空间域点扩散函数对成像结果进行去模糊化，能够提高地震偏移成像去模糊化的效率。
图2示出了本发明实施例提供的偏移成像去模糊化方法中步骤101的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了提高确定成像波数矢量的效率，如图2所示，步骤101，在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量，包括：
步骤201，在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的传播矢量；
步骤202，根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量；
步骤203，根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量。
在确定目标散射点处的成像波数矢量时，首先在每道地震道对应的震源及接收点处进行射线追踪，以确定目标散射点处的传播矢量。进而在确定目标散射点处的传播矢量后，基于目标散射点处的传播矢量及地震波频散关系，可以确定目标散射点处的波数矢量。其中，地震波频散能够引起地震波的传播能量衰减和相速度频散，从而影响地震成像分辨率，因此在地震波正演和反演中需要时分重视对能量衰减和相速度频散的问题。地震波频散指的是波数与频率的关系。在得到目标散射点处的波数矢量后，即可基于目标散射点处的波数矢量确定目标散射点处的成像波数矢量。
具体的，可以采用下述公式表示地震波频散关系：
其中，
表示目标散射点处的波数矢量，ω表示地震波频率，
表示目标散射点处的传播矢量。
在本发明实施例中，首先在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的传播矢量，然后再根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量，最后根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量，能够提高确定成像波数矢量的效率，进而提高偏移成像效率。
图3-1示出了本发明实施例提供的偏移成像去模糊化方法中地震成像射线对示意，图3示出了本发明实施例提供的偏移成像去模糊化方法中步骤101的另一实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了进一步提高确定成像波数矢量的效率，目标散射点处的传播矢量包括目标散射点处的震源传播矢量及目标散射点处的接收点传播矢量，目标散射点处的波数矢量包括目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量。如图3-1及图3所示，所示，步骤202，根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量，包括：
步骤301，根据目标散射点处的震源传播矢量及地震波频散关系，确定目标散射点处的震源波数矢量；及
步骤302，根据目标散射点处的接收点传播矢量及地震波频散关系，确定目标散射点处的接收点波数矢量；
步骤203，根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量，包括：
步骤303，根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量。
在本发明实施例中，目标散射点处的传播矢量包括目标散射点处的震源传播矢量及目标散射点处的接收点传播矢量。相应的，目标散射点处的波数矢量包括目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量。此处为描述方便，我们分别采用
及
表示目标散射点处的震源传播矢量、目标散射点处的接收点传播矢量、目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量。
具体的，在确定目标散射点处的成像波数矢量时，首先在地震道对应的震源及接收点进行射线追踪，基于射线追踪确定目标散射点处的震源传播矢量
及目标散射点处的接收点传播矢量
在得到目标散射点处的震源传播矢量
后，可以根据目标散射点处的震源传播矢量
及上述提到的地震波频散关系，确定目标散射点处的震源波数矢量
具体的，可以采用下述公式表述：
其中，
表示目标散射点处的震源波数矢量，ω表示地震波频率，
表示目标散射点处的震源传播矢量。
在得到目标散射点处的接收点传播矢量
后，可以根据目标散射点处的接收点传播矢量
及上述提到的地震波频散关系，确定目标散射点处的接收点波数矢量
具体的，可以采用下述公式表述：
其中，
表示目标散射点处的接收点波数矢量，ω表示地震波频率，
表示目标散射点处的接收点传播矢量。
在分别确定目标散射点处的震源波数矢量
及目标散射点处的接收点波数矢量
后，即可基于目标散射点处的震源波数矢量
及目标散射点处的接收点波数矢量
确定目标散射点处的成像波数矢量。此处为表述方便，我们采用
表示目标散射点处的成像波数矢量。
在本发明的一实施例中，为了进一步提高确定成像波数矢量的效率，步骤303，根据目标散射点处的震源波数矢量
及目标散射点处的接收点波数矢量
确定目标散射点处的成像波数矢量
包括：
将目标散射点处的接收点波数矢量
与目标散射点处的震源波数矢量
相减的结果作为目标散射点处的成像波数矢量
具体的，可以采用如下公式确定目标散射点处的成像波数矢量
其中，
表示目标散射点处的成像波数矢量，
表示目标散射点处的接收点波数矢量，
表示目标散射点处的震源波数矢量。
在本发明实施例中，根据目标散射点处的震源传播矢量及地震波频散关系，确定目标散射点处的震源波数矢量；根据目标散射点处的接收点传播矢量及地震波频散关系，确定目标散射点处的接收点波数矢量，最后根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量，能够进一步提高确定成像波数矢量的效率，进而提高偏移成像的效率。
图4示出了本发明实施例提供的偏移成像去模糊化方法中步骤102的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了提高确定波数域点扩散函数的效率，如图4所示，步骤102，根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数，包括：
步骤401，根据目标散射点处的成像波数矢量及从地震数据中提取的震源子波，确定每道地震道的波数域点扩散函数；
步骤402，对每道地震道的波数域点扩散函数进行累加，确定地震数据的波数域点扩散函数。
在确定地震数据的波数域点扩散函数时，鉴于地震数据包括多个地震道，首先基于目标散射点处的成像波数矢量及从地震数据中提取的震源子波，确定每道地震道的波数域点扩散函数。
具体的，可以基于下述公式确定每道地震道的波数域点扩散函数：
其中，
表示每道地震道的波数域点扩散函数，F(ω)表示从地震数据中提取的震源子波，ω表示地震波频率，
表示波数域点扩散函数的波数矢量，
表示目标散射点处的成像波数矢量。其中，波数域点扩散函数的波数矢量
为已知量，本发明实施例对此不再具体阐述。
在得到每道地震道的波数域点扩散函数后，进而对每道地震道的波数域点扩散函数进行累加，即可得到地震数据的波数域点扩散函数。
在本发明实施例中，首先根据目标散射点处的成像波数矢量及从地震数据中提取的震源子波，确定每道地震道的波数域点扩散函数，进而对每道地震道的波数域点扩散函数进行累加，即可得到地震数据的波数域点扩散函数，能够提高确定波数域点扩散函数的效率。
图5示出了本发明实施例提供的偏移成像去模糊化方法中步骤103的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了提高确定空间域点扩散函数的效率，如图5所示，步骤103，根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数，包括：
步骤501，对地震数据的波数域点扩散函数进行二维反傅里叶变换，确定地震数据的空间域点扩散函数。
在将地震数据的波数域点扩散函数转换为地震数据的空间域点扩散函数，可以对地震数据的波数域点扩散函数进行二维反傅里叶变换即可。其中，对地震数据的波数域点扩散函数进行二维反傅里叶变换，以得到地震数据的空间域点扩散函数，为本领域技术人员熟知的现有技术部分，本发明实施例对此不再详细赘述。
在本发明实施例中，对地震数据的波数域点扩散函数进行二维反傅里叶变换，确定地震数据的空间域点扩散函数，能够提高确定空间域点扩散函数的效率，进而提高偏移成像的效率。
本发明实施例中还提供了一种偏移成像去模糊化装置，如下面的实施例所述。由于这些装置解决问题的原理与偏移成像去模糊化方法相似，因此这些装置的实施可以参见方法的实施，重复之处不再赘述。
图6示出了本发明实施例提供的偏移成像去模糊化装置的功能模块，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
参考图6，所述偏移成像去模糊化装置所包含的各个模块用于执行图1对应实施例中的各个步骤，具体请参阅图1以及图1对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述偏移成像去模糊化装置包括射线追踪模块601、波数域确定模块602、空间域确定模块603及去模糊化模块604。
射线追踪模块601，用于在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量。
波数域确定模块602，用于根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数。
空间域确定模块603，用于根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数。
去模糊化模块604，用于利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，形成地震数据的最终成像结果。
在本发明实施例中，射线追踪模块601基于射线追踪确定目标散射点处的成像波数矢量，进而波数域确定模块602及空间域确定模块603根据目标散射点处的成像波数矢量及震源子波，即可得到地震数据的空间域点扩散函数，最后去模糊化模块604利用得到的地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化。本发明实施例无需进行正演和偏移，只需基于目标散射点处的成像波数矢量及震源子波，即可快速确定地震数据的空间域点扩散函数，进而利用得到的空间域点扩散函数对成像结果进行去模糊化，能够提高地震偏移成像去模糊化的效率。
图7示出了本发明实施例提供的偏移成像去模糊化装置中射线追踪模块601的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了提高确定成像波数矢量的效率，参考图7，所述射线追踪模块601所包含的各个单元用于执行图2对应实施例中的各个步骤，具体请参阅图2以及图2对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述射线追踪模块601包括传播矢量确定单元701、波数矢量确定单元702及成像波数矢量确定单元703。
传播矢量确定单元701，用于在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的传播矢量。
波数矢量确定单元702，用于根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量。
成像波数矢量确定单元703，用于根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量。
在本发明实施例中，首先传播矢量确定单元701在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的传播矢量，然后波数矢量确定单元702再根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量，最后成像波数矢量确定单元703根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量，能够提高确定成像波数矢量的效率，进而提高偏移成像效率。
图8示出了本发明实施例提供的偏移成像去模糊化装置中射线追踪模块601的另一结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了进一步提高确定成像波数矢量的效率，目标散射点处的传播矢量包括目标散射点处的震源传播矢量及目标散射点处的接收点传播矢量，目标散射点处的波数矢量包括目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量。参考图6，所述射线追踪模块601所包含的各个单元用于执行图3对应实施例中的各个步骤，具体请参阅图3以及图3对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述波数矢量确定单元702包括震源波数矢量确定子单元801及接收点波数矢量确定子单元802。
震源波数矢量确定子单元801，用于根据目标散射点处的震源传播矢量及地震波频散关系，确定目标散射点处的震源波数矢量。及
接收点波数矢量确定子单元802，用于根据目标散射点处的接收点传播矢量及地震波频散关系，确定目标散射点处的接收点波数矢量。
成像波数矢量确定单元703，包括成像波数矢量确定子单元803。
成像波数矢量确定子单元803，用于根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量。
在本发明的一实施例中，为了进一步提高确定成像波数矢量的效率，成像波数矢量确定子单元，具体用于将目标散射点处的接收点波数矢量与目标散射点处的震源波数矢量相减的结果作为目标散射点处的成像波数矢量。
在本发明实施例中，震源波数矢量确定子单元801根据目标散射点处的震源传播矢量及地震波频散关系，确定目标散射点处的震源波数矢量；接收点波数矢量确定子单元802根据目标散射点处的接收点传播矢量及地震波频散关系，确定目标散射点处的接收点波数矢量，最后成像波数矢量确定子单元803根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量，能够进一步提高确定成像波数矢量的效率，进而提高偏移成像的效率。
图9示出了本发明实施例提供的偏移成像去模糊化装置中波数域确定模块602的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了提高确定波数域点扩散函数的效率，参考图9，所述波数域确定模块602所包含的各个单元用于执行图4对应实施例中的各个步骤，具体请参阅图4以及图4对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述波数域确定模块602包括单道波数域确定单元901及累加单元902。
单道波数域确定单元901，用于根据目标散射点处的成像波数矢量及从地震数据中提取的震源子波，确定每道地震道的波数域点扩散函数。
累加单元902，用于对每道地震道的波数域点扩散函数进行累加，确定地震数据的波数域点扩散函数。
在本发明实施例中，首先单道波数域确定单元901根据目标散射点处的成像波数矢量及从地震数据中提取的震源子波，确定每道地震道的波数域点扩散函数，进而累加单元902对每道地震道的波数域点扩散函数进行累加，即可得到地震数据的波数域点扩散函数，能够提高确定波数域点扩散函数的效率。
图10示出了本发明实施例提供的偏移成像去模糊化装置中空间域确定模块603的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，为了提高确定空间域点扩散函数的效率，参考图10，所述空间域确定模块603所包含的各个单元用于执行图5对应实施例中的各个步骤，具体请参阅图5以及图5对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述空间域确定模块603包括变换单元1001。
变换单元1001，用于对地震数据的波数域点扩散函数进行二维反傅里叶变换，确定地震数据的空间域点扩散函数。
在本发明实施例中，变换单元1001对地震数据的波数域点扩散函数进行二维反傅里叶变换，确定地震数据的空间域点扩散函数，能够提高确定空间域点扩散函数的效率，进而提高偏移成像的效率。
图11示出了本发明实施例提供的利用本发明的偏移成像去模糊化方法得到的某地震数据的波数域点扩散函数示意，图12示出了本发明实施例提供的利用本发明的偏移成像去模糊化方法得到的某地震数据的空间域点扩散函数示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
如图11所示，为利用本发明提供的偏移成像去模糊化方法得到的某地震数据的波数域点扩散函数示意；如图12所示，为利用本发明提供的偏移成像去模糊化方法得到的某地震数据的空间域点扩散函数示意。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述移成像去模糊化方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述移成像去模糊化方法的计算机程序。
综上所述，本发明实施例基于射线追踪确定目标散射点处的成像波数矢量，进而根据目标散射点处的成像波数矢量及震源子波，即可得到地震数据的空间域点扩散函数，最后利用得到的地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化。本发明实施例无需进行正演和偏移，只需基于目标散射点处的成像波数矢量及震源子波，即可快速确定地震数据的空间域点扩散函数，进而利用得到的空间域点扩散函数对成像结果进行去模糊化，能够提高地震偏移成像去模糊化的效率。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种偏移成像去模糊化方法，其特征在于，包括：
在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量；
根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数；
根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数；
利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，形成地震数据的最终成像结果；
基于下述公式确定每道地震道的波数域点扩散函数：
其中，
表示每道地震道的波数域点扩散函数，F(ω)表示从地震数据中提取的震源子波，ω表示地震波频率，
表示波数域点扩散函数的波数矢量，
表示目标散射点处的成像波数矢量。
2.如权利要求1所述的偏移成像去模糊化方法，其特征在于，在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量，包括：
在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的传播矢量；
根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量；
根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量。
3.如权利要求2所述的偏移成像去模糊化方法，其特征在于，目标散射点处的传播矢量包括目标散射点处的震源传播矢量及目标散射点处的接收点传播矢量，目标散射点处的波数矢量包括目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量，包括：
根据目标散射点处的震源传播矢量及地震波频散关系，确定目标散射点处的震源波数矢量；及
根据目标散射点处的接收点传播矢量及地震波频散关系，确定目标散射点处的接收点波数矢量；
根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量，包括：
根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量。
4.如权利要求3所述的偏移成像去模糊化方法，其特征在于，根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量，包括：
将目标散射点处的接收点波数矢量与目标散射点处的震源波数矢量相减的结果作为目标散射点处的成像波数矢量。
5.如权利要求1所述的偏移成像去模糊化方法，其特征在于，根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数，包括：
对地震数据的波数域点扩散函数进行二维反傅里叶变换，确定地震数据的空间域点扩散函数。
6.一种偏移成像去模糊化装置，其特征在于，包括：
射线追踪模块，用于在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的成像波数矢量；
波数域确定模块，用于根据目标散射点对应的成像波数矢量及从地震数据中提取的震源子波，确定由每道地震道的波数域点扩散函数累加形成的地震数据的波数域点扩散函数；
空间域确定模块，用于根据地震数据的波数域点扩散函数确定地震数据的空间域点扩散函数；
去模糊化模块，用于利用地震数据的空间域点扩散函数对地震数据的初步成像结果进行去模糊化，形成地震数据的最终成像结果；
基于下述公式确定每道地震道的波数域点扩散函数：
其中，
表示每道地震道的波数域点扩散函数，F(ω)表示从地震数据中提取的震源子波，ω表示地震波频率，
表示波数域点扩散函数的波数矢量，
表示目标散射点处的成像波数矢量。
7.如权利要求6所述的偏移成像去模糊化装置，其特征在于，射线追踪模块包括：
传播矢量确定单元，用于在每道地震道对应的震源及接收点进行射线追踪，确定目标散射点处的传播矢量；
波数矢量确定单元，用于根据目标散射点处的传播矢量及地震波频散关系，确定目标散射点处的波数矢量；
成像波数矢量确定单元，用于根据目标散射点处的波数矢量，确定目标散射点处的成像波数矢量。
8.如权利要求7所述的偏移成像去模糊化装置，其特征在于，目标散射点处的传播矢量包括目标散射点处的震源传播矢量及目标散射点处的接收点传播矢量，目标散射点处的波数矢量包括目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，波数矢量确定单元包括：
震源波数矢量确定子单元，用于根据目标散射点处的震源传播矢量及地震波频散关系，确定目标散射点处的震源波数矢量；及
接收点波数矢量确定子单元，用于根据目标散射点处的接收点传播矢量及地震波频散关系，确定目标散射点处的接收点波数矢量；
成像波数矢量确定单元，包括：
成像波数矢量确定子单元，用于根据目标散射点处的震源波数矢量及目标散射点处的接收点波数矢量，确定目标散射点处的成像波数矢量。
9.如权利要求8所述的偏移成像去模糊化装置，其特征在于，成像波数矢量确定子单元，具体用于将目标散射点处的接收点波数矢量与目标散射点处的震源波数矢量相减的结果作为目标散射点处的成像波数矢量。
10.如权利要求6所述的偏移成像去模糊化装置，其特征在于，空间域确定模块包括：
变换单元，用于对地震数据的波数域点扩散函数进行二维反傅里叶变换，确定地震数据的空间域点扩散函数。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述偏移成像去模糊化方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述偏移成像去模糊化方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010081486.1A CN113219525B (zh) | 2020-02-06 | 2020-02-06 | 偏移成像去模糊化方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010081486.1A CN113219525B (zh) | 2020-02-06 | 2020-02-06 | 偏移成像去模糊化方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113219525A CN113219525A (zh) | 2021-08-06 |
| CN113219525B true CN113219525B (zh) | 2023-04-25 |
Family
ID=77085580
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010081486.1A Active CN113219525B (zh) | 2020-02-06 | 2020-02-06 | 偏移成像去模糊化方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113219525B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115166817B (zh) * | 2022-06-30 | 2023-02-17 | 哈尔滨工程大学 | 一种基于冰层模态群慢度差特征的冰声定位方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108562937A (zh) * | 2018-03-15 | 2018-09-21 | 东北石油大学 | 一种地震成像方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO322089B1 (no) * | 2003-04-09 | 2006-08-14 | Norsar V Daglig Leder | Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder |
| US7860333B2 (en) * | 2007-01-09 | 2010-12-28 | University Of Utah Research Foundation | Systems and methods for deblurring data corrupted by shift variant blurring |
| US10677949B2 (en) * | 2017-01-05 | 2020-06-09 | China Petroleum & Chemical Corporation | Anisotropy matching filtering for attenuation of seismic migration artifacts |
| CN107422379B (zh) * | 2017-07-27 | 2019-01-11 | 中国海洋石油集团有限公司 | 基于局部自适应凸化方法的多尺度地震全波形反演方法 |
| CN108445532B (zh) * | 2018-02-12 | 2019-11-08 | 中国石油天然气集团有限公司 | 一种深度域反偏移方法及装置 |
-
2020
- 2020-02-06 CN CN202010081486.1A patent/CN113219525B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108562937A (zh) * | 2018-03-15 | 2018-09-21 | 东北石油大学 | 一种地震成像方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113219525A (zh) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Seismic trace interpolation for irregularly spatial sampled data using convolutional autoencoder | |
| Bonar et al. | Denoising seismic data using the nonlocal means algorithm | |
| CN110806602B (zh) | 基于深度学习的智能化地震数据随机噪声压制方法 | |
| CN102998706B (zh) | 一种衰减地震数据随机噪声的方法及*** | |
| CN105425301A (zh) | 一种频率域三维不规则地震数据重建方法 | |
| JP2020508454A (ja) | 指向性波動場イメージングを用いた地球物理学的画像生成 | |
| Chen et al. | Seismic inversion by hybrid machine learning | |
| CN103105623A (zh) | 一种地震勘探中的数据波形处理方法 | |
| CN113219525B (zh) | 偏移成像去模糊化方法及装置 | |
| CN111458749B (zh) | 应用于被动源地震勘探的面波与体波的分离方法及*** | |
| WO2021055152A1 (en) | Noise attenuation methods applied during simultaneous source deblending and separation | |
| Ferreira et al. | Internal multiple removal in offshore Brazil seismic data using the inverse scattering series | |
| Sun et al. | Full-waveform inversion using a learned regularization | |
| CN110161565A (zh) | 一种地震数据重建方法 | |
| CN113031072B (zh) | 虚同相轴层间的多次波压制方法、装置及设备 | |
| Li et al. | Flexibility-residual BiSeNetV2 for GPR image decluttering | |
| CN115170428A (zh) | 一种声波远探测成像图的降噪方法 | |
| CN102105814B (zh) | 用于地震道分析的***和方法 | |
| CN105319594A (zh) | 一种基于最小二乘参数反演的傅里叶域地震数据重构方法 | |
| Li et al. | Simple framework for the contrastive learning of visual representations-based data-driven tight frame for seismic denoising and interpolation | |
| Aghamiry et al. | Large-scale highly-accurate extended full waveform inversion using convergent Born series | |
| CN106405504A (zh) | 联合剪切波变换和奇异值分解的探地雷达数据去噪方法 | |
| Qiu et al. | Mitigating the cycle-skipping of full-waveform inversion: An optimal transport approach with exponential encoding | |
| CN114063159A (zh) | 一种地震面波速度确定方法及装置 | |
| CN113311485A (zh) | 地震沉积特征加强滤波方法及装置 |
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 |