CN112925025B - 基于面波的地震数据静校正方法、装置及电子设备 - Google Patents
基于面波的地震数据静校正方法、装置及电子设备 Download PDFInfo
- Publication number
- CN112925025B CN112925025B CN202110114539.XA CN202110114539A CN112925025B CN 112925025 B CN112925025 B CN 112925025B CN 202110114539 A CN202110114539 A CN 202110114539A CN 112925025 B CN112925025 B CN 112925025B
- Authority
- CN
- China
- Prior art keywords
- wave
- model
- data
- static correction
- wave 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.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 129
- 230000003068 static effect Effects 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims abstract description 96
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims description 20
- 239000006185 dispersion Substances 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000004587 chromatography analysis Methods 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 16
- 238000010586 diagram Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration 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
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 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
- 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/53—Statics correction, e.g. weathering layer or transformation to a datum
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
技术领域
本发明涉及数据处理技术领域，尤其是涉及一种基于面波的地震数据静校正方法、装置及电子设备。
背景技术
相关技术中，通常采用从同一地震波数据中提取的初至波对反射波进行静校正。具体可以采用基于层状模型的反演方法或基于网格的层析法，上述方法需要采用表层速度或初始模型；表层速度一般通过直达波估算得到，初始速度模型可以为通过直达波和折射波的分支点建立的简单层状模型；然而由于采用的表层速度或初始模型不准确，上述方法对反射波的静校正结果较差。
发明内容
有鉴于此，本发明的目的在于提供一种基于面波的地震数据静校正方法、装置及电子设备，以提高对反射波的静校正应用效果。
第一方面，本发明实施例提供了一种基于面波的地震数据静校正方法，包括：获取地震波数据；根据地震波数据，确定反射波数据、初至波数据及面波数据；根据初至波数据及面波数据，确定静校正量；通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。
结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，根据初至波数及面波数据，确定静校正量的步骤，包括：基于面波数据，生成横波速度模型；基于横波速度模型及地震波数据，计算静校正量。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第二种可能的实施方式，其中，基于面波数据，生成横波速度模型的步骤，包括：对面波数据进行预处理；通过第一预设算法从经过预处理的面波数据中提取频散曲线；根据频散曲线和预先获取的半波长深度经验关系建立水平层状初始模型；基于第二预设算法及水平层状初始模型进行反演，得到横波速度模型。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，基于面波数据，生成横波速度模型的步骤，包括：对面波数据进行预处理；基于第三预设算法及经过预处理的所述面波数据生成全部频散谱；根据全部频散谱和预先获取的半波长深度经验关系建立水平层状初始模型；基于第四预设算法及水平层状初始模型进行反演，得到横波速度模型。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第四种可能的实施方式，其中，基于横波速度模型及初至波数据，计算静校正量的步骤，包括：从所述初至波数据中拾取初至时间；基于横波速度模型得到模型参数；基于初至时间及模型参数，生成近地表模型；基于近地表模型、预设的最终基准面及预设的替换速度，计算静校正量。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第五种可能的实施方式，其中，上述模型参数包括表层速度或初始速度模型；基于初至时间及模型参数，生成近地表模型的步骤，包括：基于预先获取的层状模型、表层速度及初至时间，生成近地表模型；或者，通过基于网格的层析法，对模型参数及初至时间进行处理，得到近地表模型。
结合第一方面的第四种可能的实施方式，本发明实施例提供了第一方面的第六种可能的实施方式，其中，上述地震波数据基于纵波勘探过程生成；基于横波速度模型及初至波数据，计算静校正量之前，方法还包括：基于初至波数据，生成初至波反演结果及面波反演结果；将初至波反演结果与面波反演结果的比值确定为比例因子；通过比例因子对横波速度模型进行速度转换，得到转换后的横波速度模型。
第二方面，本发明实施例还提供一种基于面波的地震数据静校正装置，包括：数据获取模块，用于获取地震波数据；数据确定模块，根据地震波数据，确定反射波数据、初至波数据及面波数据；静校正量确定模块，用于根据初至波数据及面波数据，确定静校正量；静校正模块，用于通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。
第三方面，本发明实施例还提供一种电子设备，包括处理器和存储器，存储器存储有能够被处理器执行的机器可执行指令，处理器执行机器可执行指令以实现上述方法。
第四方面，本发明实施例还提供一种机器可读存储介质，机器可读存储介质存储有机器可执行指令，机器可执行指令在被处理器调用和执行时，机器可执行指令促使处理器实现上述方法。
本发明实施例带来了以下有益效果：
本发明实施例提供了一种基于面波的地震数据静校正方法、装置及电子设备，在获取地震波数据后，首先根据地震波数据，确定反射波数据、初至波数据及面波数据；然后根据初至波数据及面波数据，确定静校正量；进而通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。该方式在对反射波进行静校正的过程中采用了面波数据，提高了对反射波的静校正的应用效果。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种原始地震单炮记录数据的处理方法的流程图；
图2为本发明实施例提供的一种基于面波的地震数据静校正方法的流程图；
图3为本发明实施例提供的另一种基于面波的地震数据静校正方法的流程图；
图4为本发明实施例提供的一种基于面波信息的原始地震单炮记录数据的处理方法的流程图；
图5为本发明实施例提供的弹炮记录数据的示例图；
图6为本发明实施例提供的面波数据的示例图；
图7为本发明实施例提供的一种地震单炮记录中面波反演的横波速度结构剖面；
图8为本发明实施例提供的现有方法的初始模型的示例图；
图9为本发明实施例提供的基于现有方法的初始模型的初至层析反演结果示例图；
图10为本发明实施例提供的基于面波的地震数据静校正方法的初始模型的示例图；
图11为本发明实施例提供的一种基于基于面波的地震数据静校正方法的初始模型的初至层析反演结果示例图；
图12为本发明实施例提供的未应用静校正的反射波的示例图；
图13为本发明实施例提供的现有方法的反射波静校正效果的示例图；
图14为本发明实施例提供的基于面波的地震数据静校正方法的反射波静校正效果示例图；
图15为本发明实施例提供的基于面波的地震数据静校正方法的速度模型和静校正结果对比示例图；
图16为本发明实施例提供的一种基于面波的地震数据静校正装置的结构示意图；
图17为本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
广泛用于油气、煤炭、金属矿、工程地质和地壳深部结构领域的反射地震勘探多次覆盖技术包括数据采集、数据处理和数据解释3个步骤，其中数据处理是将采集的原始地震单炮记录数据(也可称为地震波数据)经过去噪、反褶积、静校正、动校正和叠加、偏移等处理，获得能够反映地下地质信息的地震剖面或地震体。
在复杂的近地表环境下，如巨厚的覆盖层、低速带剧烈的横向变化等，静校正处理是反射地震资料处理中最重要的一步,静校正的处理效果会直接影响最终叠加剖面的信噪比和同相轴构造的真实性。
基于单炮初至拾取和反演的近表面模型重建静校正技术，已经成为地震反射数据处理中的行业标准工作流程。如图1所示，该流程主要包括利用野外原始数据(即上述原始地震单炮记录数据，该数据经过观测系统定义)，通过滤波等数据处理手段得到反射波、直达波、折射波和面波。在目前的静校正过程中，没有使用面波。通常基于初至波进行初至时间拾取，其中，初至波在近炮检距时为直达波。随着炮检距的增大被浅层折射波所代替；然后得到初始模型或表层速度，最后生成近地表速度模型；进而采用该模型对反射波进行静校正。
单炮数据初至静校正技术主要分为两类：一类是基于层状模型的反演方法，主要包括广义互换法(Palmer 1981)和广义线性反演(Hampson and Russell 1984)；另一类是基于网格的层析法，包括射线追踪层析静校正(Zhang and Toksoz 1998)和无射线追踪层析静校正(Osypov 1999)两种。
基于层状模型的折射静校正和基于网格的初至层析静校正是目前地震资料处理中常用的两种静校正方法。对于层状模型折射静校正方法，需要预先给出表层速度，校正效果一定程度上依赖于表层速度的准确度；对于层析静校正，则需提供初始速度-深度模型，由于模型参数与明显依赖于射线路径的旅行时之间的非线性关系，层析反演的结果对初始模型很敏感，虽然理论上对于全局性的反演方法，层析射线追踪反演结果不依赖初始模型，但是，实际中的反演次数不可能进行无限次，因此其结果对初始模型的依赖性较强。
总之，在实际中，更接近真实结构的表层速度和初始模型对初至折射静校正方法的应用效果具有重要作用。然而，在现有技术方案中，如图1所示，单炮初至拾取和反演的近表面模型重建静校正技术的表层速度一般通过直达波估算，初始速度模型也是通过直达波和折射波的分支点建立的简单层状模型，显然这些初始模型在复杂地区显得过于简单，影响静校正的应用效果。
基于此，本发明实施例提供的一种基于面波的地震数据静校正方法、装置以及电子设备，可以应用于地震数据的处理过程中。
为便于对本实施例进行理解，首先对本发明实施例所公开的一种基于面波的地震数据静校正方法进行详细介绍。
本发明实施例提供了一种基于面波的地震数据静校正方法，如图2所示，该方法包括以下步骤：
步骤S100，获取地震波数据。
该地震波数据通常是通过爆炸物埋藏在地下并引爆后，由检测装置检测到的。该地震波数据包括全波场信息，包括反射波、直达波、折射波、面波、绕射波、声波等。
步骤S102，根据地震波数据，确定反射波数据、初至波数据及面波数据。
上述反射波数据和面波数据可以通过一维频率滤波、二维速度滤波等数据处理方法得到。反射是波动在传播媒质的两种波阻抗不同的分界面中才会发生。地震波波前首先到达某个观测点，此点介质的质点开始发生振动的时刻，并记录下来的波，称为初至波，或波的初至(时间)。利用初至波可以研究风化层的速度和厚度。面波(surface wave)是地震波的一种，主要沿地表附近传播，能量最大，速度低，往往最后被记录到。面波具有能量强、频率低和明显的频散特性，也携带地下地质信息。
步骤S104，根据初至波数据及面波数据，确定静校正量。
通常情况下，可以基于面波数据，生成横波速度模型(也称为“横波速度结构”或“剪切波速度-深度模型”)，然后再基于横波速度模型及初至时间，计算静校正量。具体而言，当采用基于层状模型的折射静校正方法时，可以根据初至波数据得到初至时间，然后由初至时间和通过横波速度模型确定的表层速度生成近地表模型，再基于近地表模型确定静校正量。当采用基于网格的初至层析静校正方法时，需要用到地震波数据中的初至波数据，再结合通过横波速度模型确定的初始速度模型，生成近地表模型，再基于近地表模型确定静校正量。
步骤S106，通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。
地震勘探解释的理论假定激发点与接收点是在一个水平面上，并且地层速度是均匀的；但实际上地面通常为不平坦的，且地下速度也是不均匀的。各个激发点深度也可能不同，所以必将影响实测的时距曲线形状。因此，在地震数据处理过程中，需要将地震数据校正到一个统一基准面，一般为水平面，该矫正过程为静校正。当得到静校正量后，采用静校正量对反射波数据进行一定的数学运算，则可以实现矫正过程，得到静校正后的反射波数据。
本发明实施例提供了一种基于面波的地震数据静校正方法，在获取地震波数据后，首先根据地震波数据，确定反射波数据、初至波数据及面波数据；然后根据初至波数据及面波数据，确定静校正量；进而通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。该方法在对反射波进行静校正的过程中采用了面波数据，提高了对反射波的静校正的应用效果。
本发明实施例还提供了另一种基于面波的地震数据静校正方法，该方法在图1所示的方法基础上实现。如图3所示，该方法包括以下步骤：
步骤S200，获取地震波数据。
步骤S202，根据地震波数据，确定反射波数据、初至波数据及面波数据。
步骤S204，对面波数据进行预处理；该过程用于提纯地震波数据的面波信息，提高面波信息的信噪比，以压制反射波、折射波和声波等干扰波。
步骤S206，通过第一预设算法从经过预处理的面波数据中提取频散曲线；具体而言，该过程将时空域的面波数据转换到频率-相速度域，计算频散谱，提取频散曲线。上述第一预设算法包括相移法、F-K法(频率波数域方法)、倾斜叠加法、τ-p法等。
步骤S208，根据频散曲线和预先获取的半波长深度经验关系建立水平层状初始模型。在具体实施过程中，上述步骤S206中可以从经过预处理的面波数据中提取全部频散谱，则步骤S208中，需要根据全部频散谱和预先获取的半波长深度经验关系建立水平层状初始模型。
步骤S210，基于第二预设算法及水平层状初始模型进行反演，得到横波速度模型；上述第二预设算法(也称为“反演方法”)通常包括两类，一类是基于优化设计的最小二乘方法反演及阻尼最小二乘法反演，另一类是基于全局搜索的算法，如遗传算法、神经网络算法、蒙特卡洛算法等。
由面波的频散曲线反演获取的地下速度信息(即上述“横波速度模型”)为横波速度结构，当反射单炮数据为横波勘探时不需要转换。
步骤S212，从初至波数据中拾取初至时间。
步骤S214，基于横波速度模型得到模型参数；在实际实现过程中，如果上述地震波数据基于纵波勘探过程生成，在静校正应用时需要地下纵波速度结构，此时最关键的一步是需将面波反演的横波速度结构转化为纵波速度结构；速度转换可以通过将横波速度乘以随深度变化的比例因子来实现。该因子受常规初始速度模型的纵波的初至旅行时反演结果的约束，经过速度转换以后，获得具有横波速度变化趋势的纵波范围内的连续速度结构。具体而言，首先基于地震波数据，生成折射波反演结果及面波反演结果，然后将折射波反演结果与面波反演结果的比值确定为比例因子，最后通过比例因子对横波速度模型进行速度转换，得到转换后的横波速度模型。如果地震波数据基于横波勘探过程生成，则不需要进行结构转化。
步骤S216，基于初至时间及模型参数，生成近地表模型。上述模型参数可以为表层速度或初始速度模型；具体地，可以基于预先获取的层状模型、表层速度及初至时间，生成近地表模型；或者，通过基于网格的层析法，对模型参数及初至时间进行处理，得到近地表模型。
步骤S218，基于近地表模型、预设的最终基准面及预设的替换速度，计算静校正量；即通过给定的最终基准面和替换速度，根据最终的近地表模型计算静校正量。
步骤S220，通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。
简而言之，在得到了横波速度模型后，可根据不同的静校正方法采取不同的应用途径建立近地表模型：当采用基于层状模型的折射静校正时，由面波信息提供横向变化的表层速度；当采用基于网格初至层析静校正时，由面波信息提供横向变化的初始速度模型。最后通过给定的最终基准面和替换速度，根据最终的近地表模型计算静校正量。
上述方法充分利用了同一地下介质在地震反射单炮记录中产生的多种不同类型的弹性波，在提高近地表速度模型精度的基础上，改善静校正方法的应用效果。与现有技术相比，主要增加了对面波的利用，如图4所示。
上述方法存在以下优点：
(1)充分利用了地震波数据中的面波信息，增加了地下结构的信息源，降低了勘探成本。
传统的表层速度调查需要进行专门的小折射、微测井和面波勘探等工作，具有较高成本。面波在单炮记录(如图5所示，其中横坐标为道号，纵坐标为时间Time；从图5中可知，震源位置Source Station在2200桩号处，Trace Number为接收排列道数，可以看出震源位置在81道附近，第1道位置为2120桩号附近)中具有能量强、频率低和明显的频散特性(如图6所示，其中横坐标为相速度Phase velocity，纵坐标为频率Frequency)。现有技术方案只将面波视为干扰波，通过噪音压制将其滤除，并未利用其频散特性提取地下速度结构。图7(横坐标为距离Distance，纵坐标为深度Depth)为上述方法使用单炮记录中的面波信息获得的沿测线的地下横波速度剖面(S-wave velocity cross-section)。
(2)考虑到现有的两类静校正方法对表层速速和初始模型的依赖性，本发明通过同源数据的面波信息改善表层速度和初始模型的精度，提高静校正应用效果。
首先，静校正处理的应用效果取决于近地表速度模型的精确度，而地下实际速度模型只有一个，同一个速度模型同时产生直达波、折射波、面波等多种类型的弹性波，都具有相同的横向速度变化趋势(参见图7和图9)。其次，由于目前的初至反演方法在实际应用中对初始模型和表层速度的依赖性较强，现有传统方法只能建立简单的层状模型(图8，也称为初始速度模型，Initial velocity model)导致反演结果(图9，最终速度模型，Finalvelocity model)横向变化和细节不够丰富，不够接近实际地下结构。因此可利用单炮记录中与初至同源的面波信息构建初始模型(图10)，提高最终反演近地表速度模型的精度(图11)，改善静校正应用效果(参见图14和图15)，其改善程度可以与未应用静校正的反射波(参见图12，其中，CMP为Common Middle Point的缩写，即CMP道集)及现有方法的反射波静校正效果(参见图13)比较而得。
对应于上述方法实施例，本发明实施例还提供一种基于面波的地震数据静校正装置，如图16所示，该装置包括：
数据获取模块1000，用于获取地震波数据；
数据确定模块1002，根据地震波数据，确定反射波数据、初至波数据及面波数据；
静校正量确定模块1004，用于根据初至波数据及面波数据，确定静校正量；
静校正模块1006，用于通过静校正量对反射波数据进行静校正，得到静校正后的反射波数据。
本发明实施例提供的基于面波的地震数据静校正装置，与上述实施例提供的基于面波的地震数据静校正方法具有相同的技术特征，所以也能解决相同的技术问题，达到相同的技术效果。
本发明实施例还提供了一种电子设备，参见图17所示，该电子设备包括处理器130和存储器131，该存储器131存储有能够被处理器130执行的机器可执行指令，该处理器130执行机器可执行指令以实现上述基于面波的地震数据静校正方法。
进一步地，图17所示的电子设备还包括总线132和通信接口133，处理器130、通信接口133和存储器131通过总线132连接。
其中，存储器131可能包含高速随机存取存储器(RAM，Random Access Memory)，也可能还包括非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。通过至少一个通信接口133(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。总线132可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图17中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
处理器130可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器130中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器130可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DigitalSignal Processing，简称DSP)、专用集成电路(Application Specific IntegratedCircuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器131，处理器130读取存储器131中的信息，结合其硬件完成前述实施例的方法的步骤。
本发明实施例还提供了一种机器可读存储介质，该机器可读存储介质存储有机器可执行指令，该机器可执行指令在被处理器调用和执行时，该机器可执行指令促使处理器实现上述基于面波的地震数据静校正方法，具体实现可参见方法实施例，在此不再赘述。
本发明实施例所提供的基于面波的地震数据静校正方法、装置和电子设备的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，网关电子设备，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (6)
1.一种基于面波的地震数据静校正方法，其特征在于，包括：
获取地震波数据；所述地震波数据基于纵波勘探过程生成；
根据所述地震波数据，确定反射波数据、初至波数据及面波数据；
基于所述面波数据，生成横波速度模型；
基于所述初至波数据数据，生成初至波反演结果及面波反演结果；
将所述初至波反演结果与所述面波反演结果的比值确定为比例因子；
通过所述比例因子对所述横波速度模型进行速度转换，得到转换后的横波速度模型；
从所述初至波数据中拾取初至时间；
基于所述横波速度模型得到模型参数；所述模型参数包括表层速度或初始速度模型；
基于预先获取的层状模型、所述表层速度及所述初至时间，生成近地表模型；
或者，通过基于网格的层析法，对所述模型参数及初至时间进行处理，得到近地表模型；
基于所述近地表模型、预设的最终基准面及预设的替换速度，计算静校正量；
通过所述静校正量对所述反射波数据进行静校正，得到静校正后的反射波数据。
2.根据权利要求1所述的方法，其特征在于，基于所述面波数据，生成横波速度模型的步骤，包括：
对所述面波数据进行预处理；
通过第一预设算法从经过预处理的所述面波数据中提取频散曲线；
根据所述频散曲线和预先获取的半波长深度经验关系建立水平层状初始模型；
基于第二预设算法及所述水平层状初始模型进行反演，得到横波速度模型。
3.根据权利要求1所述的方法，其特征在于，基于所述面波数据，生成横波速度模型的步骤，包括：
对所述面波数据进行预处理；
基于第三预设算法及经过预处理的所述面波数据生成全部频散谱；
根据所述全部频散谱和预先获取的半波长深度经验关系建立水平层状初始模型；
基于第四预设算法及所述水平层状初始模型进行反演，得到横波速度模型。
4.一种基于面波的地震数据静校正装置，其特征在于，包括：
数据获取模块，用于获取地震波数据；所述地震波数据基于纵波勘探过程生成；
数据确定模块，根据所述地震波数据，确定反射波数据、初至波数据及面波数据；
静校正量确定模块，用于基于所述面波数据，生成横波速度模型；基于所述初至波数据数据，生成初至波反演结果及面波反演结果；将所述初至波反演结果与所述面波反演结果的比值确定为比例因子；通过所述比例因子对所述横波速度模型进行速度转换，得到转换后的横波速度模型；从所述初至波数据中拾取初至时间；基于所述横波速度模型得到模型参数；所述模型参数包括表层速度或初始速度模型；基于预先获取的层状模型、所述表层速度及所述初至时间，生成近地表模型；或者，通过基于网格的层析法，对所述模型参数及初至时间进行处理，得到近地表模型；基于所述近地表模型、预设的最终基准面及预设的替换速度，计算静校正量；
静校正模块，用于通过所述静校正量对所述反射波数据进行静校正，得到静校正后的反射波数据。
5.一种电子设备，其特征在于，包括处理器和存储器，所述存储器存储有能够被所述处理器执行的机器可执行指令，所述处理器执行所述机器可执行指令以实现权利要求1-3任一项所述的方法。
6.一种机器可读存储介质，其特征在于，所述机器可读存储介质存储有机器可执行指令，所述机器可执行指令在被处理器调用和执行时，所述机器可执行指令促使所述处理器实现权利要求1-3任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110114539.XA CN112925025B (zh) | 2021-01-26 | 2021-01-26 | 基于面波的地震数据静校正方法、装置及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110114539.XA CN112925025B (zh) | 2021-01-26 | 2021-01-26 | 基于面波的地震数据静校正方法、装置及电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112925025A CN112925025A (zh) | 2021-06-08 |
| CN112925025B true CN112925025B (zh) | 2023-02-10 |
Family
ID=76167349
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110114539.XA Active CN112925025B (zh) | 2021-01-26 | 2021-01-26 | 基于面波的地震数据静校正方法、装置及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112925025B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101907727A (zh) * | 2010-08-17 | 2010-12-08 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
| CN104133244A (zh) * | 2013-05-02 | 2014-11-05 | 中国石油化工股份有限公司 | 一种震源横波勘探的近地表静校正处理方法 |
| WO2017106127A1 (en) * | 2015-12-18 | 2017-06-22 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101819279B (zh) * | 2010-03-26 | 2012-05-23 | 中国石油集团西北地质研究所 | 基于变差函数拟合重构的高精度三维静校正方法 |
| CN102937721B (zh) * | 2012-11-07 | 2015-07-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用初至波走时的有限频层析成像方法 |
| CN103630937B (zh) * | 2013-12-09 | 2016-08-03 | 中国石油天然气集团公司 | 对地震记录进行横波静校正处理的方法和装置 |
| CN104730579B (zh) * | 2013-12-18 | 2018-02-13 | 中国石油化工股份有限公司 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
-
2021
- 2021-01-26 CN CN202110114539.XA patent/CN112925025B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101907727A (zh) * | 2010-08-17 | 2010-12-08 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
| CN104133244A (zh) * | 2013-05-02 | 2014-11-05 | 中国石油化工股份有限公司 | 一种震源横波勘探的近地表静校正处理方法 |
| WO2017106127A1 (en) * | 2015-12-18 | 2017-06-22 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
Non-Patent Citations (1)
| Title |
|---|
| 基于交叉梯度约束的地震初至纵波与瑞雷面波联合反演;傅磊，刘四新;《地球物理学报》;20161231;第59卷;正文第4469页第4段、第4470页第5段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112925025A (zh) | 2021-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101086535B (zh) | 表征油藏随时间演变的方法和程序 | |
| US9470811B2 (en) | Creating a high resolution velocity model using seismic tomography and impedance inversion | |
| AU2006279961B2 (en) | Method of accounting for wavelet stretch in seismic data | |
| JP5379163B2 (ja) | 地震探査データのスペクトルシェーピングインバージョン法及びマイグレーション法 | |
| CN107817526B (zh) | 叠前地震道集分段式振幅能量补偿方法及系统 | |
| US9952341B2 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| EP2113792A1 (en) | Spectral shaping inversion and migration of seismic data | |
| WO2017136133A1 (en) | Efficient seismic attribute gather generation with data synthesis and expectation method | |
| RU2570827C2 (ru) | Гибридный способ для полноволновой инверсии с использованием способа одновременных и последовательных источников | |
| CN113740901A (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
| CN108508489B (zh) | 一种基于波形微变化匹配的地震反演方法 | |
| Ding et al. | Reverse time migration (RTM) imaging of iron-oxide deposits in the Ludvika mining area, Sweden | |
| EA030770B1 (ru) | Система и способ адаптивной сейсмической оптики | |
| CN112925025B (zh) | 基于面波的地震数据静校正方法、装置及电子设备 | |
| Li et al. | Wavelet-based detection of singularities in acoustic impedances from surface seismic reflection data | |
| US6418380B1 (en) | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration | |
| US8380440B2 (en) | 3D residual binning and flatness error correction | |
| Onnis et al. | Retrieving shallow shear-wave velocity profiles from 2D seismic-reflection data with severely aliased surface waves | |
| CN113126146A (zh) | 针对峡谷地区复杂地质的探测方法、控制装置和存储介质 | |
| Wu et al. | A new method for high resolution well-control processing of post-stack seismic data | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| CN104166162B (zh) | 基于迭代三参数小波变换的缝洞发育带检测方法 | |
| CN112394393B (zh) | 一种crp道集数据体重构的方法 | |
| CN113267808B (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 |