CN111965701B - 近地表结构反演方法及系统 - Google Patents
近地表结构反演方法及系统 Download PDFInfo
- Publication number
- CN111965701B CN111965701B CN201910418338.1A CN201910418338A CN111965701B CN 111965701 B CN111965701 B CN 111965701B CN 201910418338 A CN201910418338 A CN 201910418338A CN 111965701 B CN111965701 B CN 111965701B
- Authority
- CN
- China
- Prior art keywords
- wave
- surface structure
- cross
- correlation
- field 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000001514 detection method Methods 0.000 claims abstract description 44
- 238000005314 correlation function Methods 0.000 claims description 45
- 238000012937 correction Methods 0.000 claims description 35
- 230000003068 static effect Effects 0.000 claims description 35
- 238000010586 diagram Methods 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 abstract description 30
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 238000011835 investigation Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 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
- 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
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示出本申请实施例中一种近地表结构反演方法的流程示意图之一。
图2a示出本申请实施例中地面地震数据单炮Z分量记录。
图2b示出本申请实施例中地面地震数据单炮X分量记录。
图2c示出本申请实施例中地面地震数据Z分量与X分量初至互相关共检波点叠加记录。
图2d示出本申请实施例中三维地面地震数据Z分量与X分量初至互相关时差反演得到的横波静校正平面分布图。
图2e示出本申请实施例中横波静校正前的转换波共检波点叠加记录。
图2f示出本申请实施例中横波静校正后的转换波共检波点叠加记录。
图2g示出本申请实施例中横波静校正后的转换波共转换点叠加记录
图3示出本申请实施例中图1中步骤S200的具体流程示意图。
图4示出本申请实施例中图1的步骤S400的具体流程示意图。
图5示出本申请实施例中一种近地表结构反演方法的流程示意图之二。
图6示出本申请实施例中一种近地表结构反演系统结构示意图之一。
图7示出本申请实施例中图6中截取模块200的具体结构示意图。
图8示出本申请实施例中图6中反演模块400的具体结构示意图。
图9示出本申请实施例中一种近地表结构反演系统结构示意图之二。
图10示出适于用来实现本申请实施例的电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本申请主要基于多波地震三分量初至互相关叠加求取近地表结构的发明构思，首先获取三分量的地面地震波场数据，然后对设定时窗数据进行互相关，通过互相关共检波点叠加求取近地表结构。
下面结合附图对本申请实施例进行详细说明。
本申请第一方面实施例提供一种近地表结构反演方法，如图1所示，包括：
S100:获取每个检波点在每次激发地震波时接收的地面地震波场数据；其中所述地面地震波场数据包括平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数；
S200:从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，并组成一对；
S300:对每一对设定时窗数据执行互相关操作，得到对应每一对设定时窗数据的互相关函数；
S400:基于所有所述互相关函数进行近地表结构反演，得到近地表结构反演结果。
本方面提供的近地表结构反演方法，首先获取三分量的地面地震波场数据，然后对设定时窗数据进行互相关，通过互相关共检波点叠加求取近地表结构，在实际生产应用中，可以得到精细的多波工区的近地表结构，对于提高多波数据的处理精度具有重要意义。
下面具体对上述步骤进行详细说明。
在步骤S100中，利用三分量检波器地面接收，在地面利用爆炸震源或人工可控震源激发，采集得到多分量的地面地震波场数据。
多分量的地面地震波场数据中包括X轴方向、Y轴方向以及Z轴方向的分量。
在一实施例中，X轴方向对应平行于炮点连线方向，Z轴对应垂直于近地表面方向，Y轴方向对应垂直于X轴和Y轴方向。
本实施例中，需要采用平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数据。即上一实施例中的X轴和Z轴方向。
一实施例中，设定时窗为初至时窗，时窗长度与风化层的纵横波速度及旅行时有关，一般选择300-1000毫秒。
在以初至时窗为设定时窗的实施例中，如图3所示，步骤S200具体包括：
S201:获取每个检波点在每次激发地震波时接收的初至波的到达时间点；
S202:以每个到达时间点为时窗起点，从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取设定时长的波场数据组成一对设定时窗数据。
进一步的，对截取的X分量和Z分量的初至时窗数据进行互相关得到互相关函数F1后，需要对对应每一个检波点进行分选(归类)，该实施例中，如图4所示，步骤S400具体包括：
S401:将所有互相关函数中对应同一个检波点的作为一组，并将处于同一组的所有互相关函数叠加形成叠加函数；
S402:根据每个叠加函数最大值对应的互相关时差计算风化层厚度；
S403:根据所述风化层厚度计算得到三维工区的横波静校正值分布图；
S404:根据所述横波静校正值分布图，进行近地表结构反演解释，得到所述近地表结构反演结果。
该实施例中，首先对截取的Z分量与X分量初至时窗数据进行互相关，得到互相关函数F1，然后对所有炮的所有检波点的互相关函数按检波点进行分选，并按检波点叠加得到叠加后的检波点互相关函数F2，然后对共检波点叠加后的互相关函数拾取互相关函数最大值对应的互相关时差dt。对得到的每个检波点互相关时差dt利用相关的对应关系公式求取风化层厚度H0。
本申请中，共检波点互相关函数叠加主要是为了提高互相关的信噪比，对于信噪比较差的互相关函数在叠加前应剔除。
一实施例中，上述对应关系公式为：
H0＝dt×(Vp0-Vs0)
这里Vp0为风化层的纵波平均速度，Vs0为风化层的横波平均速度，二者均可通过小折射或微测井调查获得，也可通过多分量单炮记录的近道折射反演得到。
即对应地，步骤S402中，风化层厚度根据风化层的纵波平均速度、风化层的横波平均速度以及所述互相关时差计算得到，所述风化层的纵波平均速度、风化层的横波平均速度可通过测井数据得到。
利用得到的风化层厚度就可以通过剥层填充法得到后续处理需要的横波静校正值，利用静校正值应用于多波地震数据进行后续的处理解释，提供可用于地质研究的地球物理成果数据。
具体的，利用风化层厚度就可以通过剥层并利用高速折射层横波速度填充法得到了后续处理需要的横波静校正值。三维工区的静校正值分布图如图2e和图2f和图2g。利用静校正值应用于多波地震数据进行后续的处理(图2e和图2f为横波静校正前后共检波点叠加对比，图2g为横波静校正后的转换波共转换点叠加记录)解释，提供可用于地质研究的地球物理成果数据(图2f为转换波最终成像剖面)。
一优选实施例中，为了减少后续的频谱分析的处理量以及误差，经过炮点接收点坐标设置、去除随机干扰等处理，可以减少后续的数据处理量，减少处理负担，可以对未经截取的数据进行预处理，即在截取设定时窗数据之前，如图5所示，还包括：
S001：对所述地面地震波场数据进行预处理。
具体实施例中，预处理的步骤包括如下至少之一：
S010:去除所述地面地震波场数据中的随机干扰；
S020:设置所述地面地震波场数据中每个检波点的坐标。
下面结合一具体场景对本申请上述实施例进行详细说明。
为了验证本申请的效果，以一中国中部某油田实际采集多波多分量地震数据进行处理，首先利用三分量检波器接收，地面爆炸源激发，采集得到多分量地震数据，并对该数据进行预处理，利用本发明的方法进行初至时窗互相关，并通过共检波点叠加求取互相关时差反演近地表风化层厚度，进而得到横波静校正参数用于后续的多分量地震数据处理。具体如下步骤：
1)利用三分量检波器地面接收，在地面利用爆炸震源激发，采集得到多分量地震波场数据。
2)对步骤1)中所得到的采集数据进行预处理，经过炮点接收点坐标设置、去除随机干扰、拾取纵波初至等三步处理，得到预处理后数据(如图2a和图2b)。
3)利用步骤2)中所得到的预处理后数据根据拾取的纵波初至时间向下截取Z分量与X分量500毫秒初至时窗数据。
4)步骤2)截取的Z分量与X分量初至时窗数据进行互相关，得到互相关函数F1。
5)对所有炮的所有检波点的互相关函数按检波点进行分选，并按检波点叠加得到叠加后的检波点互相关函数F2(如图2c)。
这里的共检波点互相关函数叠加主要是为了提高互相关的信噪比，对于信噪比较差的互相关函数在叠加前已剔除。
6)对步骤5)的共检波点叠加后的互相关函数拾取互相关函数最大值对应的互相关时差dt。
7)对步骤6)得到的每个检波点互相关时差dt利用下式求取风化层厚度H0：
H0＝dt×(Vp0-Vs0)
这里Vp0为风化层的纵波平均速度，Vs0为风化层的横波平均速度，二者均可通过微测井调查获得。
8)利用步骤7)得到的风化层厚度就可以通过剥层并利用高速折射层横波速度填充法得到了后续处理需要的横波静校正值。三维工区的静校正值分布图如图2d。
9)利用步骤8)中所得到的静校正值应用于多波地震数据进行后续的处理(图2e和图2f为横波静校正前后共检波点叠加对比，图2g为横波静校正后的转换波共转换点叠加记录)解释，提供可用于地质研究的地球物理成果数据(图2f为转换波最终成像剖面)。
可以知晓，本场景中对三分量地面地震数据通过初至互相关和共检波点叠加进行近地表结构调查的方法，通过本方法的处理后得到近地表风化层厚度等关键参数，可以计算横波静校正，为后续的多波地震数据处理解释提供了保障，也可用于其他与近地表结构有关的应用。
基于与本申请第一方面的实施例相同的发明构思，本申请第二方面提供一种近地表结构反演系统，如图6所示，包括：
获取模块100，获取每个检波点在每次激发地震波时接收的地面地震波场数据；其中所述地面地震波场数据包括平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数据；
截取模块200，从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，并组成一对；
互相关模块300，对每一对设定时窗数据执行互相关操作，得到对应每一对设定时窗数据的互相关函数；
反演模块400，基于所有所述互相关函数进行近地表结构反演，得到近地表结构反演结果。
本方面提供的近地表结构反演系统，首先获取三分量的地面地震波场数据，然后对设定时窗数据进行互相关，通过互相关共检波点叠加求取近地表结构，在实际生产应用中，可以得到精细的多波工区的近地表结构，对于提高多波数据的处理精度具有重要意义。
一实施例中，基于与上述实施例相同的理由，如图7所示，所述截取模块200，包括：
初至时间点获取单元201，获取每个检波点在每次激发地震波时接收的初至波的到达时间点；
设定时窗数据截取单元202，以每个到达时间点为时窗起点，从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取设定时长的波场数据组成一对设定时窗数据。
一实施例中，基于与上述实施例相同的理由，如图8所示，所述反演模块400，包括：
叠加函数单元401，将所有互相关函数中对应同一个检波点的作为一组，并将处于同一组的所有互相关函数叠加形成叠加函数；
风化层厚度计算单元402，根据每个叠加函数最大值对应的互相关时差计算风化层厚度；
横波静校正值分布图生成单元403，根据所述风化层厚度生成三维工区的横波静校正值分布图；
反演解释单元404，根据所述横波静校正值分布图，进行近地表结构反演解释，得到所述近地表结构反演结果。
一实施例中，所述风化层厚度计算单元根据风化层的纵波平均速度、风化层的横波平均速度以及所述互相关时差计算所述风化层厚度；其中所述风化层的纵波平均速度、风化层的横波平均速度可通过测井数据得到。具体的，其对应关系公式为：
H0＝dt×(Vp0-Vs0)
这里Vp0为风化层的纵波平均速度，Vs0为风化层的横波平均速度，二者均可通过小折射或微测井调查获得，也可通过多分量单炮记录的近道折射反演得到。
一实施例中，如图9所示，本方面的系统还包括：预处理模块001，对所述地面地震波场数据进行预处理。
进一步的，所述预处理模块包括如下至少一个：
随机干扰去除单元，去除所述地面地震波场数据中的随机干扰；
接收点坐标设置单元，设置所述地面地震波场数据中每个检波点的坐标。
同理，结合本申请第一方面中的具体场景可以明确，对三分量地面地震数据通过初至互相关和共检波点叠加进行近地表结构调查的方法，通过该方法的处理后得到近地表风化层厚度等关键参数，可以计算横波静校正，为后续的多波地震数据处理解释提供了保障，也可用于其他与近地表结构有关的应用。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种电子设备的具体实施方式，参见图10，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；
所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
S100:获取每个检波点在每次激发地震波时接收的地面地震波场数据；其中所述地面地震波场数据包括平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数；
S200:从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，并组成一对；
S300:对每一对设定时窗数据执行互相关操作，得到对应每一对设定时窗数据的互相关函数；
S400:基于所有所述互相关函数进行近地表结构反演，得到近地表结构反演结果。
从上述描述可知，本申请提供的电子设备在实际生产应用中，可以得到精细的多波工区的近地表结构，对于提高多波数据的处理精度具有重要意义。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
S100:获取每个检波点在每次激发地震波时接收的地面地震波场数据；其中所述地面地震波场数据包括平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数；
S200:从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，并组成一对；
S300:对每一对设定时窗数据执行互相关操作，得到对应每一对设定时窗数据的互相关函数；
S400:基于所有所述互相关函数进行近地表结构反演，得到近地表结构反演结果。
从上述描述可知，本申请提供的计算机可读存储介质在实际生产应用中，可以得到精细的多波工区的近地表结构，对于提高多波数据的处理精度具有重要意义。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (12)
1.一种近地表结构反演方法，其特征在于，包括：
获取每个检波点在每次激发地震波时接收的地面地震波场数据；其中所述地面地震波场数据包括平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数据；
从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，并组成一对；
对每一对设定时窗数据执行互相关操作，得到对应每一对设定时窗数据的互相关函数；
基于所有所述互相关函数进行近地表结构反演，得到近地表结构反演结果；
所述基于所述互相关函数进行近地表结构反演，得到近地表结构反演结果，包括：
将所有互相关函数中对应同一个检波点的作为一组，并将处于同一组的所有互相关函数叠加形成叠加函数；
根据每个叠加函数最大值对应的互相关时差计算风化层厚度；
根据所述风化层厚度计算得到三维工区的横波静校正值分布图；
根据所述横波静校正值分布图，进行近地表结构反演解释，得到所述近地表结构反演结果。
2.根据权利要求1所述的近地表结构反演方法，其特征在于，所述从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，包括：
获取每个检波点在每次激发地震波时接收的初至波的到达时间点；
以每个到达时间点为时窗起点，从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取设定时长的波场数据组成一对设定时窗数据。
3.根据权利要求1所述的近地表结构反演方法，其特征在于，所述根据每个叠加函数最大值对应的互相关时差计算风化层厚度，包括：
根据风化层的纵波平均速度、风化层的横波平均速度以及所述互相关时差计算所述风化层厚度；
其中所述风化层的纵波平均速度、风化层的横波平均速度可通过测井数据得到。
4.根据权利要求1所述的近地表结构反演方法，其特征在于，在截取设定时窗数据之前，还包括：
对所述地面地震波场数据进行预处理。
5.根据权利要求4所述的近地表结构反演方法，其特征在于，所述对所述地面地震波场数据进行预处理，包括如下至少之一：
去除所述地面地震波场数据中的随机干扰；
设置所述地面地震波场数据中每个检波点的坐标。
6.一种近地表结构反演系统，其特征在于，包括：
获取模块，获取每个检波点在每次激发地震波时接收的地面地震波场数据；其中所述地面地震波场数据包括平行于炮点连线方向的第一分量波场数据和垂直于近地表面的第二分量波场数据；
截取模块，从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取一设定时窗数据，并组成一对；
互相关模块，对每一对设定时窗数据执行互相关操作，得到对应每一对设定时窗数据的互相关函数；
反演模块，基于所有所述互相关函数进行近地表结构反演，得到近地表结构反演结果；
所述反演模块，包括：
叠加函数单元，将所有互相关函数中对应同一个检波点的作为一组，并将处于同一组的所有互相关函数叠加形成叠加函数；
风化层厚度计算单元，根据每个叠加函数最大值对应的互相关时差计算风化层厚度；
横波静校正值分布图生成单元，根据所述风化层厚度生成三维工区的横波静校正值分布图；
反演解释单元，根据所述横波静校正值分布图，进行近地表结构反演解释，得到所述近地表结构反演结果。
7.根据权利要求6所述的近地表结构反演系统，其特征在于，所述截取模块，包括：
初至时间点获取单元，获取每个检波点在每次激发地震波时接收的初至波的到达时间点；
设定时窗数据截取单元，以每个到达时间点为时窗起点，从每个检波点在每次激发地震波时接收的两个分量波场数据中分别截取设定时长的波场数据组成一对设定时窗数据。
8.根据权利要求6所述的近地表结构反演系统，其特征在于，所述风化层厚度计算单元根据风化层的纵波平均速度、风化层的横波平均速度以及所述互相关时差计算所述风化层厚度；
其中所述风化层的纵波平均速度、风化层的横波平均速度可通过测井数据得到。
9.根据权利要求6所述的近地表结构反演系统，其特征在于，还包括：
预处理模块，对所述地面地震波场数据进行预处理。
10.根据权利要求9所述的近地表结构反演系统，其特征在于，所述预处理模块包括如下至少一个：
随机干扰去除单元，去除所述地面地震波场数据中的随机干扰；
接收点坐标设置单元，设置所述地面地震波场数据中每个检波点的坐标。
11.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至5任一项所述的近地表结构反演方法的步骤。
12.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至5任一项所述的近地表结构反演方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910418338.1A CN111965701B (zh) | 2019-05-20 | 2019-05-20 | 近地表结构反演方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910418338.1A CN111965701B (zh) | 2019-05-20 | 2019-05-20 | 近地表结构反演方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111965701A CN111965701A (zh) | 2020-11-20 |
| CN111965701B true CN111965701B (zh) | 2023-09-26 |
Family
ID=73358262
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910418338.1A Active CN111965701B (zh) | 2019-05-20 | 2019-05-20 | 近地表结构反演方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111965701B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1379460A (en) * | 1971-07-29 | 1975-01-02 | Continental Oil Co | Method of obtaining static corrections for seismic data |
| GB1380886A (en) * | 1973-04-13 | 1975-01-15 | Chevron Res | Static corrections for seismic traces by cross-correlation method |
| CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
| CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
| CN107678063A (zh) * | 2017-09-25 | 2018-02-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
| CN109407151A (zh) * | 2018-12-18 | 2019-03-01 | 吉林大学 | 基于波场局部相关时移的时间域全波形反演方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9348049B2 (en) * | 2012-01-05 | 2016-05-24 | Cgg Services Sa | Simultaneous joint estimation of the P-P and P-S residual statics |
| CN109923440B (zh) * | 2017-10-12 | 2021-03-05 | 南方科技大学 | 面波勘探方法及终端设备 |
-
2019
- 2019-05-20 CN CN201910418338.1A patent/CN111965701B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1379460A (en) * | 1971-07-29 | 1975-01-02 | Continental Oil Co | Method of obtaining static corrections for seismic data |
| GB1380886A (en) * | 1973-04-13 | 1975-01-15 | Chevron Res | Static corrections for seismic traces by cross-correlation method |
| CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
| CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
| CN107678063A (zh) * | 2017-09-25 | 2018-02-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
| CN109407151A (zh) * | 2018-12-18 | 2019-03-01 | 吉林大学 | 基于波场局部相关时移的时间域全波形反演方法 |
Non-Patent Citations (3)
| Title |
|---|
| 《A simple method for resolving large converted‐wave (P-SV) statics》;Peter W. Cary等;《geophysics》;19930301;第58卷(第3期);第429-433页 * |
| 《利用多分量地震数据反演近地表横波速度》;李晓明等;《石油地球物理勘探》;20120815;第47卷(第4期);第532-536页 * |
| 《多分量地震处理和解释技术》;黄中玉等;《油气地球物理技术新进展——第74届SEG年会论文概要》;20040630;第77-89页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111965701A (zh) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111221037B (zh) | 解耦弹性逆时偏移成像方法和装置 | |
| CN107894613B (zh) | 弹性波矢量成像方法、装置、存储介质及设备 | |
| CN102053261A (zh) | 一种地震数据处理方法 | |
| CN111025386B (zh) | 一种无分离假象的纵横波分离方法 | |
| CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
| CN106772593B (zh) | 绕射波的成像方法及装置 | |
| CN104730574A (zh) | 构建近地表结构模型的方法 | |
| CN101021567A (zh) | 地震资料处理的炮道集波动方程叠前深度偏移并行计算方法 | |
| CN102053262A (zh) | 获得地震转换波的方位速度的方法及处理地震数据的方法 | |
| EP4102260A1 (en) | Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data | |
| CN103941289B (zh) | 一种绕射点成像的方法和设备 | |
| CN113740906A (zh) | 一种水下垂直缆地震波干涉成像方法及装置 | |
| CN111965701B (zh) | 近地表结构反演方法及系统 | |
| CN110579799B (zh) | 一种等旅行时间间隔的地震采集观测方法及系统 | |
| CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
| CN111123361B (zh) | 垂直地震剖面地震数据规则化重建方法及装置、存储介质 | |
| CN111077577B (zh) | 井地联合储层描述方法及装置 | |
| CN109143398B (zh) | 一种自动网格层析深度域速度的建模方法 | |
| US10761228B2 (en) | Method to calculate acquisition illumination | |
| CN113721296A (zh) | 远震数据处理方法及装置 | |
| CN111965704B (zh) | 井中地震数据的光缆次生震荡噪声压制方法及系统 | |
| CN110579798B (zh) | 一种等反射角度间隔的地震采集观测方法及系统 | |
| CN112285782B (zh) | 近地表地震波吸收衰减调查方法及装置 | |
| CN111965703B (zh) | 井中地震数据中单炮记录的下行波反褶积处理方法及系统 | |
| JP7470614B2 (ja) | 切羽前方探査システムにおける切羽前方三次元速度構造推定方法 |
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 |