CN113031072A - 虚同相轴层间的多次波压制方法、装置及设备 - Google Patents
虚同相轴层间的多次波压制方法、装置及设备 Download PDFInfo
- Publication number
- CN113031072A CN113031072A CN202110266472.1A CN202110266472A CN113031072A CN 113031072 A CN113031072 A CN 113031072A CN 202110266472 A CN202110266472 A CN 202110266472A CN 113031072 A CN113031072 A CN 113031072A
- Authority
- CN
- China
- Prior art keywords
- target interface
- multiples
- interbed
- wave
- seismic data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000001629 suppression Effects 0.000 claims abstract description 49
- 239000011229 interlayer Substances 0.000 claims description 61
- 239000010410 layer Substances 0.000 claims description 47
- 238000003825 pressing Methods 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 18
- 238000010276 construction Methods 0.000 claims description 16
- 238000011410 subtraction method Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000007493 shaping process Methods 0.000 claims description 5
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000000875 corresponding effect Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- RVRCFVVLDHTFFA-UHFFFAOYSA-N heptasodium;tungsten;nonatriacontahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W] RVRCFVVLDHTFFA-UHFFFAOYSA-N 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 101100134058 Caenorhabditis elegans nth-1 gene Proteins 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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
技术领域
本发明涉及层间多次波压制技术领域，尤其涉及一种虚同相轴层间的多次波压制方法、装置及设备。
背景技术
相比于海上全程多次波，产生陆上层间多次波的强反射来源多、传播类型更多，使得多次波场特征复杂，周期规律性不好，特别是深层地震有效信号能量普遍较弱，多次波和有效波的能量、速度差异小，识别和压制难度更大。而层间多次波的存在对后续地震资料的解释、反演和应用都会带来不利影响，导致构造和油气识别的精度下降，钻探成功率降低，对油气的勘探开发带来极大的风险。
传统的虚同相轴层间多次波压制方法通过一步预测相减可有效压制某一界面产生的所有层间多次波，但是对另一界面产生的层间多次波也需要通过相应的虚同相轴重新构建，这对于实际数据处理来说计算量是非常大的，且在实际数据中很难完整的提取构成层间多次波的一次波，因此该方法有很大的局限性。
发明内容
有鉴于此，本发明提供一种虚同相轴层间的多次波压制方法、装置及设备，至少部分解决现有技术中存在的问题。
第一方面，本发明提供了一种虚同相轴层间的多次波压制方法，方法包括：
通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型；所述虚同相轴，是根据目标界面逆时一次反射波及位于目标界面以下界面所产生的一次波互相关得到；
基于初始层间多次波模型，通过预设的卷积因子对目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据；
从包含一次波和层间多次波的地震数据中，减去预测层间多次波数据，得到目标压制结果。
在一些实施例中，通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型，具体包括：
通过目标界面的一次反射波构建虚同相轴；虚同相轴，为目标界面对应逆时一次反射波与第一地震数据的褶积；第一地震数据不含目标界面及其以上界面的一次反射波；
构建目标界面的层间多次波；层间多次波，为虚同相轴与第二地震数据的褶积，第二地震数据是切除目标界面及其以上界面所产生的所有反射波得到的在一些实施例中，通初始层间多次波模型包括：
目标界面的层间多次波第n次迭代后的第三地震数据，为第n-1迭代后的第四地震数据与卷积因子的乘积。
在一些实施例中，在所述基于所述初始层间多次波模型，通过预设的卷积因子对所述目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据中，具体包括：确定卷积因子。
在一些实施例中，确定卷积因子，具体包括：
通过共聚焦点层间多次波压制法，得到目标界面的层间多次波压制结果；层间多次波压制结果，为基于自适应相减法，由上一次多次波压制结果得到的；
通过界面算子，将目标界面的层间多次波压制结果转换为关于卷积因子的关系式，求取卷积因子。
第二方面，本发明提供了一种多次波压制装置，装置包括：
第一构建模块，用于通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型；所述虚同相轴，是根据目标界面逆时一次反射波及位于所述目标界面以下界面所产生的一次波互相关得到；
迭代反演模块，用于基于所述初始层间多次波模型，通过预设的卷积因子对所述目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据；
计算模块，用于从包含所述一次波和所述层间多次波的地震数据中，减去所述预测层间多次波数据，得到目标压制结果。
第三方面，本发明提供了一种多次波压制设备，设备包括：
处理器，以及存储有计算机程序指令的存储器；处理器读取并执行计算机程序指令，以实现如上述任意一实施例的多次波压制方法。
第四方面，本发明提供了一种计算机存储介质，计算机存储介质上存储有计算机程序指令，计算机程序指令被处理器执行时实现如上述任意一实施例的多次波压制方法。
第五方面，本发明提供了一种地震数据的处理方法，方法包括：
读取野外采集数据，得到含有坐标信息的原始地震数据；
对原始地震数据进行预处理，并进行子波整形和振幅补偿，得到待压制层间多次波数据；
对待压制层间多次波数据进行处理，处理包括如上述任意一实施例的多次波压制方法。
第六方面，本发明提供了一种地震数据的处理设备，包括：处理器，以及存储有计算机程序指令的存储器；处理器读取并执行计算机程序指令，以实现上述的地震数据的处理方法。
本发明提供的虚同相轴层间的多次波压制方法、装置及设备，通过虚同相轴构建目标界面的层间多次波，并引入卷积因子构建初始层间多次波模型，将层间多次波压制转换为层间多次波的迭代反演过程，这样只需将某一目标界面的层间多次波作为初始多次波模型进行迭代以得到所有层的层间多次波，计算量较小。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的虚同相轴层间的多次波压制方法流程示意图；
图2A是图1所示的步骤S101在一个示例中的流程示意图；
图2B是构建关于界面A的虚同相轴的示意图，其中“★”表示褶积；
图2C是基于图2B构建虚同相轴后，构建关于界面A的层间多次波；
图2D是根据图2C得到的界面A的相关的层间多次波；
图2E是构建关于界面B的虚同相轴；
图2F是基于图2E构建虚同相轴后，构建关于界面B的层间多次波；
图3是图1所示的步骤S102在一个示例中的流程示意图；
图4是本发明实施例提供的虚同相轴层间的多次波压制装置结构示意图；
图5是本发明实施例提供的虚同相轴层间的多次波压制设备硬件结构示意图；
图6是本发明实施例提供的地震数据处理方法的流程示意图；
图7A是本发明验证时的速度模型图的示意图；
图7B是单炮记录的示意图；
图7C是虚同相轴层间多次波压制前后的地震记录示意图；
图8是本发明实施例提供的地震数据处理设备硬件结构示意图。
具体实施方式
下面结合附图对本发明实施例进行详细描述。
需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合；并且，基于本公开中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。
需要说明的是，下文描述在所附权利要求书的范围内的实施例的各种方面。应显而易见，本文中所描述的方面可体现于广泛多种形式中，且本文中所描述的任何特定结构及/或功能仅为说明性的。基于本公开，所属领域的技术人员应了解，本文中所描述的一个方面可与任何其它方面独立地实施，且可以各种方式组合这些方面中的两者或两者以上。举例来说，可使用本文中所阐述的任何数目个方面来实施设备及/或实践方法。另外，可使用除了本文中所阐述的方面中的一或多者之外的其它结构及/或功能性实施此设备及/或实践此方法。
多次波的产生依赖于传播介质的波阻抗差异。自由表面多次波在地表或海水表面发生下行反射，通常表现有强能量、明显的周期性等性质，而层间多次波的下行反射界面全部位于地下。当地下存在强反射界面时，反射波向上传播到该层界面会再次向下反射，由于地下岩层速度大，因此层间多次波与对应的一次波时差差距较小，在接收记录中会严重干扰目标层的有效信号，使得解释工作者难以辨别目标地质体的反射波，进而为地震属性提取、反演等增加困难。
多次波压制方法主要分为两类，一类是滤波方法，另一类是预测相减方法。滤波方法主要利用一次波与多次波的时差关系和周期特征进行多次波的识别和压制，如抛物Radon滤波和双曲Radon滤波等，该类方法计算效率高，算法容易实现，当有效波和多次波之间的动校正时差较大时，可取得满意的多次波压制效果。但对于复杂介质，如速度梯度较小(或速度反转)的情况，或构造变化剧烈的介质，应用滤波法难以有效识别有效波和多次波，往往得不到理想的多次波压制结果。预测相减方法基于波动理论，能更好的适应于复杂介质的情况。该类多次波压制方法包括逆散射级数法，该方法无需地下的先验信息，预测一次可以得到与所有界面相关的同阶层间多次波，在没有有效手段区分有效波和多次波时，是层间多次波压制较为有效的方法，但是该方法存在计算量大，预测远偏移距多次波效果差的缺点。共聚焦点技术适合复杂介质条件的多次波压制，但该方法一次只能预测得到与某一界面相关的层间多次波，在一定程度上依赖初始速度模型以求取准确的聚焦算子。基于Marchenko的层间多次波压制方法在模型数据上取得了很好的效果，也可应用于实际数据，但是对于层间多次波发育且信噪比较低的陆上地震资料多次波压制还存在计算不稳定的情况。传统技术中还有一种数据驱动构造层间多次波的方法，将地下散射点移到表面，但是存在对于实际数据无法进行精确的层间多次波压制问题，在一些技术中解决该问题可以通过引入虚源点的概念即基于虚同相轴的层间多次波压制，但这种通过一步预测相减有效压制某一界面产生的所有层间多次波，对另一界面产生的层间多次波则需要通过相应的虚同相轴重新构建，这对于实际数据处理来说计算量是非常大的，且对观测系统要求较高，不能摆脱对人工操作的依赖。
为此，本公开实施例提供了一种虚同相轴层间的多次波压制方法、装置及设备，旨在解决传统的虚同相轴层间多次波压制方法计算量大、实际数据中难以提取一次波构建层间多次波的弊端，同时，提高层间多次波预测精度，也在一定程度上解决后续多次波自适应相减方法中的非线性问题。
图1示出了本公开一个实施例提供的虚同相轴层间的多次波压制方法的流程示意图。如图1所示，该多次波压制方法包括：
S101.通过构建目标界面的虚同相轴，得到关于目标界面的初始层间多次波模型；虚同相轴，是根据目标界面逆时一次反射波及位于所述目标界面以下界面所产生的一次波互相关得到；
S102.基于初始层间多次波模型，通过预设的卷积因子对目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据；
S103.从包含一次波和层间多次波的地震数据中，减去预测层间多次波压制数据，得到目标压制结果。
本实施例中，通过步骤S101～S103，利用预设的卷积因子，通过对初始层间多次波模型进行迭代反演得到所有层的层间多次波，然后利用自适应匹配相减法，将预测的层间多次波从原始数据中减去。该方法与虚同相轴方法相比无需构建每一层的层间多次波，只需将某一层的层间多次波作为初始多次波模型进行迭代以得到所有层的层间多次波，使层间多次波预测过程简化并降低了计算量，有利于实际数据的应用。由于本发明在层间多次波预测中隐含的考虑了表面算子的空间变化，并在层间多次波预测不断迭代反演更新中实现层间算子的空间变化，因此在迭代反演过程中，层间多次波预测精度提高，未被预测或者预测不完整的层间多次波随迭代次数的增加逐渐完善，同时在一定程度上解决了后续多次波自适应相减方法中的非线性问题。
示例性的，在步骤S101中通过构建目标界面的虚同相轴，得到关于目标界面的初始层间多次波模型，如图2A所示，具体可以包括：
S201.通过目标界面的一次反射波构建虚同相轴；虚同相轴，为目标界面对应逆时一次反射波与第一地震数据的褶积；第一地震数据不含目标界面及其以上界面的一次反射波；
S202.构建目标界面的层间多次波；层间多次波，为虚同相轴与第二地震数据的褶积，第二地震数据是切除目标界面及其以上界面所产生的所有反射波得到的。
在步骤S201.中，如果P0(xS,xR)为某一层地下界面(也即目标界面)相关的一次反射波，P0'(xS,xR)为不含P0(xS,xR)的第一地震数据，也即第一地震数据中不含目标界面及其以上界面的一次反射波。那么虚同相轴即可定义为逆时一次反射波与不含该一次波的数据的褶积，虚同相轴可以表示为：
其中上标“—”代表复共轭。PV(xS,xR)为构建的虚同相轴，x表示地表界面上的任意一点，xR表示检波点水平方向坐标，xS表示炮点水平方向坐标。
虚同相轴是一种无法在地震记录中被观察到的波，但是虚同相轴的存在可以用来构建与该地下界面相关的层间多次波，其中“相关”指的是，在目标界面发生过至少一次向下散射。将虚同相轴与切除了该界面一次反射波后的第二地震数据做褶积，就可以得到与该界面相关的层间多次波。
因此步骤S202中，可以将虚同相轴与切除了对应界面一次反射波后的数据(也即第二地震数据，该数据是切除目标界面及其以上界面所产生的所有反射波得到的)做褶积，就可以得到与该界面相关的层间多次波，即式(2)：
MI(xS,xR)＝∫PV(xS,x)P′0(x,xR)dx (2)。
参考图2B，是构建关于界面A的虚同相轴，图2B中“★”表示褶积，在其他图中也表示褶积；同时参考图2C是基于虚同相轴，构建关于界面A的层间多次波，图2D是根据图2C得到的界面A的相关的层间多次波。图2E是构建关于界面B的虚同相轴，图2F是基于图2E示出的构建虚同相轴，构建关于界面B的层间多次波。
从虚同相轴的构建中可以看出，对于不同界面产生的层间多次波需要相应的一次波和不含这个一次波的数据褶积重新构建，这对构造复杂的地震数据来说计算量是非常大的且一次波的提取也非常困难，因此传统技术中仅通过虚同相轴进行层间多次波压制很难广泛的应用于实际数据。
本实施例中，构建得到目标界面的层间多次波后，可以得到初始层间多次波模型，后续基于预设的卷积因子进行迭代。示例性的，初始层间多次波模型包括：目标界面的层间多次波第n次迭代后的第三地震数据，为第n-1迭代后的第四地震数据与卷积因子的乘积。
在本示例中，卷积因子可以借助于现有的共聚焦点层间多次波压制方法中的界面算子推算得到。即本实施例中，得到初始层间多次波模型，具体可以包括：
设经过n次迭代压制层间多次波后的地震记录{P(zS0,zR0)}n′可以由上一次迭代结果
与类似于界面算子
的卷积因子T(n-1)的乘积表示，则有
也即，式(4)为基于卷积因子，构建的初始层间多次波模型。
本实施例中，通过步骤S102.基于初始层间多次波模型，通过预设的卷积因子对目标界面的层间多次波进行迭代压制反演，可以得到所有层的预测层间多次波压制数据。迭代压制反演的过程即卷积因子的确定过程，通过对目标界面的层间多次波进行迭代压制反演，得到确定的卷积因子，继而可以得到确定的目标界面层间多次波模型。
示例性的，本发明方法中步骤S102可以包括：S310.确定卷积因子。如图3所示，步骤S310具体可以包括：
S301.通过共聚焦点层间多次波压制法，得到目标界面的层间多次波压制结果；层间多次波压制结果，为基于自适应相减法，由上一次多次波压制结果得到的；
S302.通过界面子，将目标界面的层间多次波压制结果转换为关于卷积因子的关系式，求取卷积因子。
其中，根据共聚焦点层间多次波压制技术法,关于界面zk的层间多次波压制后可表示为：
由式(5)可知，界面算子
为：
将式(6)代入式(5)后，可以得到
在实际应用中，矩阵的逆可以根据最小平方的形式进行求取，则式(7)又可以表示为
式(8)中
表示关于界面zk的层间多次波，若考虑所有的层间多次波，公式(8)可以进一步表示为：
式(9)中，H表示共轭转置，则卷积因子T(n-1)为：
确定卷积因子后，通过步骤S103根据自适应相减法，从包含所述一次波和所述层间多次波的地震数据中，减去预测层间多次波压制数据，得到目标压制结果。则基于式(10)，可以得到最终的目标压制结果，该目标压制结果可以表示为：
由上可知，通过包含式(9)至式(10)的MSI(Multiple suppression inversion，多重抑制反演)算法中，不需要显式的界面算子和显式的震源矩阵，界面算子被原始地震数据
及先前两步的多次波压制结果
和
所替代。
并且本发明引入卷积因子T后，通过对初始层间多次波模型进行迭代反演可以得到任意一层(也即所有层)的层间多次波，然后利用自适应匹配相减法，将预测的层间多次波从原始数据中减去。该方法与传统虚同相轴方法相比无需构建每一层的层间多次波，只需将某一层的层间多次波作为初始多次波模型进行迭代以得到所有层的层间多次波，相当于完全的数据驱动，可以使层间多次波预测过程简化并降低了计算量，有利于实际数据的应用。由于本发明在层间多次波预测中实际隐含的考虑了界面算子的空间变化，并在层间多次波预测不断迭代反演更新中实现层间界面算子的空间变化，因此在迭代反演过程中，层间多次波预测精度提高，未被预测或者预测不完整的层间多次波随迭代次数的增加逐渐完善，同时在一定程度上解决了后续多次波自适应相减方法中的非线性问题。
故而，由于传统的虚同相轴层间多次波压制方法，对于不同界面产生的层间多次波需要相应的一次波和不含这个一次波的数据褶积重新构建，这对构造复杂的地震数据来说提取一次波是非常困难的，所以很难广泛的应用于实际数据。但本发明方法无需提取每一层的一次波来构建每一层产生的层间多次波，因此本发明的方法适应性强。并且传统的虚同相轴层间多次波压制方法需构建每一层的层间多次波，计算量是非常大，但本发明通过引入卷积因子T，将层间多次波压制转变为迭代反演的问题，只需将某一层的层间多次波作为初始多次波模型进行迭代以得到所有层的层间多次波，经过实践验证，本发明方法中一般迭代3次就可以得到较好的结果，计算量较小，具有良好的实际地震资料层间多次波压制的应用前景。
图4示出了本发明提供的一种多次波压制装置结构示意图。如图4所示的，装置包括：
第一构建模块401，用于通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型；所述虚同相轴，是根据目标界面逆时一次反射波及位于所述目标界面以下界面所产生的一次波互相关得到；
迭代反演模块402，用于基于所述初始层间多次波模型，通过预设的卷积因子对所述目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据；
计算模块403，用于从包含所述一次波和所述层间多次波的地震数据中，减去所述预测层间多次波数据，得到目标压制结果。
示例性的，第一构建模块401可以执行上述图1中示出的步骤S101，迭代反演模块402可以执行上述图1中示出的步骤S102，计算模块403可以执行上述图1中示出的步骤S103。
需要说明的是，上述方法实施例涉及的各步骤的所有相关内容均可以援引到对应功能模块的功能描述，并能达到其相应的技术效果，为简洁描述，在此不再赘述。
图5示出了本发明提供的一种多次波压制设备的硬件结构示意图。如图5所示的，设备可以包括：
处理器501以及存储有计算机程序指令的存储器502。
具体地，上述处理器501可以包括中央处理器(Central Processing Unit，CPU)，或者特定集成电路(Application Specific Integrated Circuit，ASIC)，或者可以被配置成实施本公开实施例的一个或多个集成电路。
存储器502可以包括用于数据或指令的大容量存储器。举例来说而非限制，存储器502可包括硬盘驱动器(Hard Disk Drive，HDD)、软盘驱动器、闪存、光盘、磁光盘、磁带或通用串行总线(Universal Serial Bus，USB)驱动器或者两个或更多个以上这些的组合。在一个实例中，存储器502可以包括可移除或不可移除(或固定)的介质，或者存储器502是非易失性固态存储器。存储器502可在综合网关容灾设备的内部或外部。
存储器502可以包括只读存储器(ROM)，随机存取存储器(RAM)，磁盘存储介质设备，光存储介质设备，闪存设备，电气、光学或其他物理/有形的存储器存储设备。因此，通常，存储器包括一个或多个编码有包括计算机可执行指令的软件的有形(非暂态)计算机可读存储介质(例如，存储器设备)，并且当该软件被执行(例如，由一个或多个处理器)时，其可操作来执行参考根据本公开的一方面的方法所描述的操作。
处理器501通过读取并执行存储器502中存储的计算机程序指令，以实现图1所示实施例中的方法/步骤S101至S104，并达到图1所示实例执行其方法/步骤达到的相应技术效果，为简洁描述在此不再赘述。
另外，结合上述实施例中的次波压制方法，本公开实施例可提供一种计算机存储介质来实现。该计算机存储介质上存储有计算机程序指令；该计算机程序指令被处理器执行时实现上述实施例中的任意一种次波压制方法。
图6示出了本发明提供的一种地震数据的处理方法的流程示意图，如图6所示，方法包括：
S601.读取野外采集数据，得到含有坐标信息的原始地震数据；
S602.对原始地震数据进行预处理，并进行子波整形和振幅补偿，得到待压制层间多次波数据；
S603.对待压制层间多次波数据进行处理；其中处理包括如图1示出的多次波压制方法中的步骤S101～S103。
在实际地震数据的处理过程中，通过对野外采集数据读取，道头编辑以及观测系统加载的流程运行之后，得到了含有坐标信息的原始地震数据，首先对地震数据进行了去除野值和异常能量的处理，然后进行子波整形和振幅补偿，得到需要压制层间多次波的数据，使用相应的多次波压制方法进行压制。对于模型数据无需进行野值、异常振幅及子波整形等处理，可以直接进行多次波压制。
本发明在方法验证中采用模型数据进行验证，首先建立如图7A所示的速度模型图，图中灰色圆圈中L反射点产生的层间多次波都表示与界面A相关的层间多次波，对应图7B示出的单炮记录中箭头AL所指示的同相轴，灰色圆圈K反射点产生的层间多次波都表示与界面B相关的层间多次波，对应图7B示出的单炮记录中箭头AK所指示的同相轴。借助于有限差分正演方法，对速度模型模型四周采用吸收边界得到只含有一次波和层间多次波的待测模拟地震数据，参考图7B为待测模拟地震数据的单炮记录。并利用地震波的传播规律对产生的层间多次波进行分类分析，最后利用改进的虚同相轴层间多次波压制方法对其进行层间多次波压制，得到不含层间多次波的地震记录，如图7C所示虚同相轴层间多次波压制前后，其中(a)为原始地震记录，(b)为迭代1次的地震记录，(c)为迭代3次的地震记录。
图8示出了本发明提供的一种地震数据的处理设备的硬件结构示意图。如图8所示的，设备可以包括：
处理器801以及存储有计算机程序指令的存储器802。
具体地，上述处理器801可以包括中央处理器(Central Processing Unit，CPU)，或者特定集成电路(Application Specific Integrated Circuit，ASIC)，或者可以被配置成实施本公开实施例的一个或多个集成电路。
存储器802可以包括用于数据或指令的大容量存储器。举例来说而非限制，存储器802可包括硬盘驱动器(Hard Disk Drive，HDD)、软盘驱动器、闪存、光盘、磁光盘、磁带或通用串行总线(Universal Serial Bus，USB)驱动器或者两个或更多个以上这些的组合。在一个实例中，存储器802可以包括可移除或不可移除(或固定)的介质，或者存储器802是非易失性固态存储器。存储器802可在综合网关容灾设备的内部或外部。
存储器802可以包括只读存储器(ROM)，随机存取存储器(RAM)，磁盘存储介质设备，光存储介质设备，闪存设备，电气、光学或其他物理/有形的存储器存储设备。因此，通常，存储器包括一个或多个编码有包括计算机可执行指令的软件的有形(非暂态)计算机可读存储介质(例如，存储器设备)，并且当该软件被执行(例如，由一个或多个处理器)时，其可操作来执行参考根据本公开的地震数据的处理方法所描述的操作。
处理器801通过读取并执行存储器802中存储的计算机程序指令，以实现图6所示实施例中的方法/步骤S601至S603，并达到图6所示实例执行其方法/步骤达到的相应技术效果，为简洁描述在此不再赘述。
另外，结合上述实施例中的次波压制方法，本公开实施例可提供一种计算机存储介质来实现。该计算机存储介质上存储有计算机程序指令；该计算机程序指令被处理器执行时实现上述实施例中的任意一种地震数据的处理方法。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
Claims (10)
1.一种虚同相轴层间的多次波压制方法，其特征在于，包括：
通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型；所述虚同相轴，是根据目标界面逆时一次反射波及位于所述目标界面以下界面所产生的一次波互相关得到；
基于所述初始层间多次波模型，通过预设的卷积因子对所述目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据；
从包含所述一次波和所述层间多次波的地震数据中，减去所述预测层间多次波数据，得到目标压制结果。
2.根据权利要求1所述的虚同相轴层间的多次波压制方法，其特征在于，所述通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型，具体包括：
通过目标界面的一次反射波构建所述虚同相轴；所述虚同相轴，为所述目标界面对应逆时一次反射波与第一地震数据的褶积；所述第一地震数据不含所述目标界面及其以上界面的一次反射波；
构建所述目标界面的层间多次波；所述层间多次波，为所述虚同相轴与第二地震数据的褶积，所述第二地震数据是切除所述目标界面及其以上界面所产生的所有反射波得到的。
3.根据权利要求1所述的虚同相轴层间的多次波压制方法，其特征在于，所述初始层间多次波模型包括：
所述目标界面的层间多次波第n次迭代后的第三地震数据，通过第n-1迭代后的第四地震数据与所述卷积因子的乘积得到。
4.根据权利要求3所述的虚同相轴层间的多次波压制方法，其特征在于，在所述基于所述初始层间多次波模型，通过预设的卷积因子对所述目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据中，具体包括：确定所述卷积因子。
5.根据权利要求4所述的虚同相轴层间的多次波压制方法，其特征在于，所述确定所述卷积因子，具体包括：
通过共聚焦点层间多次波压制法，得到所述目标界面的层间多次波压制结果；所述层间多次波压制结果，为基于自适应相减法，由上一次多次波压制结果得到的；
通过界面算子，将目标界面的层间多次波压制结果转换为关于所述卷积因子的关系式，求取所述卷积因子。
6.一种多次波压制装置，其特征在于，所述装置包括：
第一构建模块，用于通过构建目标界面的虚同相轴，得到关于所述目标界面的初始层间多次波模型；所述虚同相轴，是根据目标界面逆时一次反射波及位于所述目标界面以下界面所产生的一次波互相关得到；
迭代反演模块，用于基于所述初始层间多次波模型，通过预设的卷积因子对所述目标界面的层间多次波进行迭代压制反演，得到所有层的预测层间多次波数据；
计算模块，用于从包含所述一次波和所述层间多次波的地震数据中，减去所述预测层间多次波数据，得到目标压制结果。
7.一种多次波压制设备，其特征在于，所述设备包括：
处理器，以及存储有计算机程序指令的存储器；所述处理器读取并执行所述计算机程序指令，以实现如权利要求1-5任意一项所述的多次波压制方法。
8.一种计算机存储介质，其特征在于，所述计算机存储介质上存储有计算机程序指令，所述计算机程序指令被处理器执行时实现如权利要求1-5任意一项所述的多次波压制方法。
9.一种地震数据的处理方法，其特征在于，包括：
读取野外采集数据，得到含有坐标信息的原始地震数据；
对所述原始地震数据进行预处理，并进行子波整形和振幅补偿，得到待压制层间多次波数据；
对所述待压制层间多次波数据进行处理，所述处理包括如权利要求1-5任意一项所述的多次波压制方法。
10.一种地震数据的处理设备，其特征在于，包括：处理器，以及存储有计算机程序指令的存储器；所述处理器读取并执行所述计算机程序指令，以实现如权利要求9所述的地震数据的处理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110266472.1A CN113031072B (zh) | 2021-03-11 | 2021-03-11 | 虚同相轴层间的多次波压制方法、装置及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110266472.1A CN113031072B (zh) | 2021-03-11 | 2021-03-11 | 虚同相轴层间的多次波压制方法、装置及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113031072A true CN113031072A (zh) | 2021-06-25 |
| CN113031072B CN113031072B (zh) | 2024-03-15 |
Family
ID=76470197
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110266472.1A Active CN113031072B (zh) | 2021-03-11 | 2021-03-11 | 虚同相轴层间的多次波压制方法、装置及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113031072B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113687417A (zh) * | 2021-08-23 | 2021-11-23 | 北京大学 | 一种三维叠前地震数据层间多次波预测和压制方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
-
2021
- 2021-03-11 CN CN202110266472.1A patent/CN113031072B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
Non-Patent Citations (3)
| Title |
|---|
| YANGHUA WANG: "Multiple prediction through inversion: A fully data-driven concept for surface-related multiple attenuation", GEOPHYSICS, vol. 69, no. 2, pages 547 - 553 * |
| 包培楠 等: "基于迭代反演的层间多次波压制方法", 2020年中国地球科学联合学术年会论文集（十四）—专题四十：地震波传播与成像、专题四十一：高压实验矿物学、岩石学与地球化学、专题四十二：地球物理人工智能和信息技术进展, pages 1222 - 1224 * |
| 崔永福 等: "虚同相轴方法及其在陆上地震层间多次波压制中的应用", 石油地球物理勘探, vol. 54, no. 6, pages 1228 - 1236 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113687417A (zh) * | 2021-08-23 | 2021-11-23 | 北京大学 | 一种三维叠前地震数据层间多次波预测和压制方法 |
| CN113687417B (zh) * | 2021-08-23 | 2022-07-22 | 北京大学 | 一种三维叠前地震数据层间多次波预测和压制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113031072B (zh) | 2024-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772583B (zh) | 一种地震绕射波分离方法和装置 | |
| CN102998706B (zh) | 一种衰减地震数据随机噪声的方法及系统 | |
| KR101564094B1 (ko) | 지하구조 영상의 품질 향상을 위해 절대값 함수가 적용된 탄성 매질에서의 역시간 구조보정 장치 및 방법 | |
| CN103926622B (zh) | 一种基于l1范数多道匹配滤波压制多次波的方法 | |
| CN104345343B (zh) | 一种复杂海底相关的层间多次波预测方法 | |
| US9952341B2 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| CN109061727A (zh) | 一种频率域粘声介质全波形反演方法 | |
| CN107065013A (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| US9594176B1 (en) | Fast beam migration using plane-wave destructor (PWD) beam forming | |
| Datta et al. | Full-waveform inversion of salt models using shape optimization and simulated annealing | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| US20230037886A1 (en) | Method and system for determination of seismic propagation velocities using nonlinear transformations | |
| Xue et al. | Reverse-time migration using multidirectional wavefield decomposition method | |
| CN113031072A (zh) | 虚同相轴层间的多次波压制方法、装置及设备 | |
| Chen et al. | Joint data and model-driven simultaneous inversion of velocity and density | |
| CN109521470B (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
| CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
| CN105353409A (zh) | 一种用于抑制全波形反演震源编码串扰噪音的方法和系统 | |
| CN107479091A (zh) | 一种提取逆时偏移角道集的方法 | |
| CN106547021A (zh) | 基于单井褶积算法建立初始模型的方法和装置 | |
| US8699298B1 (en) | 3D multiple prediction and removal using diplets | |
| CN111914609B (zh) | 井震联合叠前地质统计学弹性参数反演方法及装置 | |
| CN113791447A (zh) | 一种反射结构导引的反射波层析反演方法 | |
| CN110858000B (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 |