CN112764107B - 地表高速覆盖地区的地震资料静校正方法及装置 - Google Patents
地表高速覆盖地区的地震资料静校正方法及装置 Download PDFInfo
- Publication number
- CN112764107B CN112764107B CN201911064857.9A CN201911064857A CN112764107B CN 112764107 B CN112764107 B CN 112764107B CN 201911064857 A CN201911064857 A CN 201911064857A CN 112764107 B CN112764107 B CN 112764107B
- Authority
- CN
- China
- Prior art keywords
- speed
- igneous rock
- coverage area
- seismic data
- correction
- 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 168
- 230000003068 static effect Effects 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 67
- 239000011435 rock Substances 0.000 claims abstract description 188
- 230000006698 induction Effects 0.000 claims description 22
- 238000003860 storage Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 abstract description 13
- 239000010410 layer Substances 0.000 description 68
- 238000004891 communication Methods 0.000 description 17
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 11
- 230000001052 transient effect Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000002344 surface layer Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012546 transfer 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
-
- 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/37—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy specially adapted for seismic systems using continuous agitation of the ground, e.g. using pulse compression of frequency swept signals for enhancement of received signals
- G01V1/375—Correlating received seismic signals with the emitted source signal
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是本申请实施例中的地表高速覆盖地区的地震资料静校正方法的步骤010至030的流程示意图；
图3是本申请实施例中的地表高速覆盖地区的地震资料静校正方法的步骤301至303的流程示意图；
图4是本申请实施例中的地表高速覆盖地区的地震资料静校正装置的结构示意图；
图5是本申请实施例中的火成岩覆盖区的近地表模型示意图；
图6是本申请实施例中的瞬变电磁勘探成果数据-地表火成岩厚度图；
图7是本申请具体应用实例中的地表高速火成岩校正过程示意图；
图8是地表高速火成岩校正后单炮初至拾取方法示意图；
图9是应用现有静校正方法与本申请具体应用实例中的地震资料静校正方法，获取的共炮检距初至和叠加剖面的对比图；
图10为本申请实施例的电子设备9600的***构成示意框图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
以某地区为例，该地区地表存在大面积的火山喷出岩，这种火山岩是某一地质时期火山活动，火山熔岩流沿断层向上溢流，覆盖在地表冷凝而成，该区域的近地表结构存在明显的速度反转现象，巨厚的高速火成岩覆盖在正常低降速度带之上。由于地表覆盖的火成岩的厚度较大，且地震波速度明显高于下伏地层，造成基于微测井和初至拾取的传统的静校正方法难以得到火成岩下伏低降速带的速度信息，把该种近地表结构等同于高速地层出露的近地表结构，将地表认为是高速顶，对于火成岩覆盖地区的炮点、检波点只是实现了高程校正，而没有对地表异常高速介质及其下伏的低降速带进行校正，难以较好地解决该地区的静校正问题，严重影响了该地区地下构造的成像质量。
基于此，本申请提出一种瞬变电磁勘探数据和地震勘探数据联合作用，获得近地表速度模型，提取静校正量，进而实现地表高速和下伏低降速带同时校正的静校正方法，能够适应上述特殊的近地表结构。
根据上述内容，本申请实施例中提供了一种地表高速覆盖地区的地震资料静校正装置，该装置可以是一服务器或客户端设备，所述客户端设备可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备和智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表和智能手环等。
在实际应用中，进行地震资料静校正的部分可以在如上述内容所述的服务器侧执行，也可以所有的操作都在所述客户端设备中完成。具体可以根据所述客户端设备的处理能力，以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备中完成，所述客户端设备还可以包括处理器。
上述的客户端设备可以具有通信模块(即通信单元)，可以与远程的服务器进行通信连接，实现与所述服务器的数据传输。所述服务器可以包括任务调度中心一侧的服务器，其他的实施场景中也可以包括中间平台的服务器，例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备，也可以包括多个服务器组成的服务器集群，或者分布式装置的服务器结构。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信，包括在本申请提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然，所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol，远程过程调用协议)、REST协议(Representational State Transfer，表述性状态转移协议)等。
具体通过下述各个实施例进行说明。
为了能够实现地表高速和下伏低降速带同时校正，进而改善地震资料成像质量，本申请提供一种执行主体是地表高速覆盖地区的地震资料静校正装置的地表高速覆盖地区的地震资料静校正方法，参见图1，该方法具体包含有：
步骤100：基于预获取的目标地表高速覆盖地区的电阻率剖面获取待测点的电阻率值，并将该电阻率值输入预获取的地表火成岩模型，以获取该待测点对应的火成岩层的实际厚度。
具体地，所述目标地表高速覆盖地区由地表向下依次覆盖有高速火成岩和低降速度带。本申请能够对高速火成岩和低降速度带对地震资料的影响同时进行校正，目标地表高速覆盖地区具体范围可根据实际情况进行设定，本申请对此不作限制。具体地，根据所述电阻率剖面，能够获取待测点水平位置坐标对应的电阻率。将该电阻率输入预获取的电阻率与火成岩层实际厚度对应关系，能够获得该电阻率对应的待测点水平位置坐标，以及该水平位置坐标对应的火成岩层的实际厚度。
步骤200：应用所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正。
具体地，根据所述地表火成岩模型和地震单炮直达波获取所述火成岩速度。应用所述火成岩层的实际厚度和预获取的火成岩速度，对地震资料进行炮、检点高速时差校正，将原来火成岩之上激发接收的数据近似转变为火成岩底部激发接收的数据，实现地表火成岩层对应的地震资料的高速校正。
步骤300：根据所述待测点的水平位置坐标，从预获取的近地表结构模型获取对应的下伏低降速层的静校正量，以实现对所述待测点的下伏低降速层对应的地震资料进行低速校正。
为了进一步提高地震资料的准确性与成像质量，在本申请一个或多个实施例中，
步骤200包含有：
步骤201：根据所述火成岩层的实际厚度和预获取的火成岩速度，获取炮、检点高速时差。
步骤202：应用所述炮、检点高速时差，对所述待测点的火成岩层对应的地震资料进行高速校正。
为了获取目标地表高速覆盖地区的电阻率剖面，进而提高地震资料静校正的，在步骤100之前还包含有：
步骤001：接收目标地表高速覆盖地区内的垂直磁感应场的归一化感应电动势。
步骤002：基于所述垂直磁感应场的归一化感应电动势获取所述目标地表高速覆盖地区的电阻率剖面。
也就是说，应用电磁勘探技术能够获取目标地表高速覆盖地区的电阻率剖面，具体地，采集所述归一化感应电动势，并将所述归一化感应电动势进行数据处理和反演，获得地下地层的电阻率剖面；根据表层火成岩和下伏地层的电性差异，进行电性层位标定和解释，得到近地表电性层位剖面，近地表电性层位剖面和地下地层的电阻率剖面结合，获取所述目标地表高速覆盖地区的电阻率剖面。
参见图2，为了获取所述地表火成岩模型，以确定待测点的火成岩层实际厚度，进而地震资料成像质量，在步骤100之前还包含有：
步骤010：获取火成岩厚度采集点火成岩层的实际厚度，其中，该火成岩厚度采集点在所述预获取的目标地表高速覆盖地区内。
具体地，所述火成岩厚度采集点包含有火成岩层和低降速带。
步骤020：从所述预获取的目标地表高速覆盖地区的电阻率剖面中获取所述火成岩厚度采集点的电阻率。
步骤030：根据所述火成岩厚度采集点的电阻率及火成岩厚度采集点火成岩层的实际厚度获取所述地表火成岩模型，其中，所述地表火成岩模型包括：所述预获取的目标地表高速覆盖地区内的电阻率与火成岩层实际厚度对应关系。
参见图3，为了生成所述近地表结构模型，进而提高地震资料在低降速带静校正的准确性，在本申请一个或多个实施例中，在步骤300之前还包含有：
步骤301：采集所述目标地表高速覆盖地区内经过高速校正后的地震资料对应的初至时间。
步骤302：剔除该初至时间中的高速直达波初至时间，以获取所述初至时间中的折射波初至时间。
步骤303：对所述折射波初至时间进行层析反演，以生成所述近地表结构模型，其中，该近地表模型的横坐标表示水平位置坐标，纵坐标表示剔除火成岩层实际厚度后的深度位置坐标，色标为速度。
参见图7，在本申请一个实施例中，所述预获取的目标地表高速覆盖地区的地震资料静校正方法，包含有：
利用瞬变电磁勘探得到的火成岩厚度数据分别提取出炮点、检波点对应位置的火成岩厚度，使用地震波在该地区地表火成岩中的传播速度对炮检点进行表层高速校正，近似实现将原来火成岩之上激发接收的数据转变为火成岩底部激发接收的数据。
从软件层面上来说，为了能够实现地表高速和下伏低降速带同时校正，进而改善地震资料成像质量，本申请提供一种用于执行所述地表高速覆盖地区的地震资料静校正方法中全部或部分内容的地表高速覆盖地区的地震资料静校正装置的实施例，所述地表高速覆盖地区的地震资料静校正装置，参见图4，具体包含有如下内容：
获取实际厚度模块10，用于基于预获取的目标地表高速覆盖地区的电阻率剖面获取待测点的电阻率值，并将该电阻率值输入预获取的地表火成岩模型，以获取该待测点对应的火成岩层的实际厚度。
高速校正模块20，用于应用所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正；
低速校正模块30，用于根据所述待测点的水平位置坐标，从预获取的近地表结构模型获取对应的下伏低降速层的静校正量，以实现对所述待测点的下伏低降速层对应的地震资料进行低速校正。
在本申请一个或多个实施例中，所述高速校正模块20包含有：
获取高速时差单元，用于根据所述火成岩层的实际厚度和预获取的火成岩速度，获取炮、检点高速时差；
高速校正单元，用于应用所述炮、检点高速时差，对所述待测点的火成岩层对应的地震资料进行高速校正。
在本申请一个或多个实施例中，所述的地表高速覆盖地区的地震资料静校正装置，还包含有：
接收二次感应电动势模块，用于接收目标地表高速覆盖地区内的垂直磁感应场的归一化感应电动势；
获取电阻率剖面模块，用于基于所述垂直磁感应场的归一化感应电动势获取所述目标地表高速覆盖地区的电阻率剖面。
在本申请一个或多个实施例中，所述的地表高速覆盖地区的地震资料静校正装置，还包含有：
采集厚度模块，用于获取火成岩厚度采集点火成岩层的实际厚度，其中，该火成岩厚度采集点在所述预获取的目标地表高速覆盖地区内；
获取电阻率模块，用于从所述预获取的目标地表高速覆盖地区的电阻率剖面中获取所述火成岩厚度采集点的电阻率；
获取对应关系模块，用于获取根据所述火成岩厚度采集点的电阻率及火成岩厚度采集点火成岩层的实际厚度获取所述地表火成岩模型，其中，所述地表火成岩模型包括：目标地表高速覆盖地区内的电阻率与火成岩层实际厚度对应关系。
在本申请一个或多个实施例中，所述的地表高速覆盖地区的地震资料静校正装置，还包含有：
采集初至模块，用于采集所述目标地表高速覆盖地区内经过高速校正后的地震资料对应的初至时间；
获取折射波初至模块，用于剔除该初至时间中的高速直达波初至时间，以获取所述初至时间中的折射波初至时间；
生成近地表结构模型模块，用于对所述折射波初至时间进行层析反演，以生成所述近地表结构模型，其中，该近地表模型的横坐标表示水平位置坐标，纵坐标表示剔除火成岩层实际厚度后的深度位置坐标，色标为速度。
为了进一步实现地表高速火成岩层和下伏低降速带同时校正，进而改善地震资料成像质量，本申请提供一种地表高速覆盖地区的地震资料静校正方法的具体应用实例，参见图5，在该具体应用实例中地表高速覆盖地区由地表向下依次覆盖有高速火成岩、低降速度带和高速层；该具体应用实例包含有如下内容：
S100：建立地表火成岩模型。
具体地，实施瞬变电磁勘探，采用人工源激发，利用接地导线或不接地回线发送直流电流，然后突然切断电源，线框内的电流将发生一个瞬间突变，这种瞬态变化的电流将在空间产生一个瞬态的磁场，感应在地下形成二次涡流场，涡流场随时间的推移不断向地下扩散，扩散的速度与地下地层的电阻率有关，不同时间扩散到不同深度；用接收线圈来观测二次电磁场信号，通过观测二次场的变化规律来研究地下介质的电性分布特征。瞬变电磁仪野外观测数据是垂直磁感应场的归一化感应电动势。通过数据处理、反演得到地下地层的电阻率剖面，根据表层火成岩和下伏地层的电性差异，进行电性层位标定、解释，得到近地表电性层位剖面，火成岩在电阻率剖面上表现为明显的高阻特征，结合实际钻井信息对电阻率反演剖面上玄武岩底界进行标定，从而可以准确查清地表火成岩的厚度和空间分布，建立地表火成岩的模型。
举例来说，火成岩覆盖区兰23井位于瞬变电磁测线上，结合该井揭示的玄武岩厚度，对电阻率反演剖面上玄武岩底界进行标定，对玄武岩的厚度变化进行预测，可以得到火成岩相对精确的厚度，进而建立该区近地表火成岩模型，本区主要为两条带状火成岩体，火成岩体宽度200～2600米、厚度在0～140米之间。与ZK06、ZK10、ZK14、兰9、兰13等5口探井钻遇火成岩资料对比，火成岩厚度的平均误差5.1％，最大误差6％。
S200：进行表层高速异常体校正。
具体地，利用瞬变电磁勘探得到的火成岩厚度模型和地震单炮直达波测量火成岩速度，对地震数据进行炮、检点高速时差校正，将原来火成岩之上激发接收的数据近似转变为火成岩底部激发接收的数据，实现地表火成岩的高速校正。
S300：对地表高速校正后的地震单炮数据进行初至拾取。
具体地，首先利用初至拾取软件自动拾取初至，然后对自动拾取到的高速直达波初至时间进行人工剔除，不参与层析反演，最后利用折射波初至时间进行层析反演，得到地表火成岩介质剥离后的近地表结构模型。参见图8，虚线示意的是初至拾取的位置8a。
S400：进行火成岩下伏的低降速层校正。
具体地，利用火成岩高速校正后数据体的折射初至时间层析反演得到高速火成岩剥离后的表层结构模型计算高速火成岩下伏低降速层的静校正量。
最终实现地表高速和下伏低降速带的同时校正，能够解决地表存在高速火成岩覆盖地区的静校正问题。
参见图6，反映本申请一个具体应用实例中瞬变电磁勘探成果数据-地表火成岩厚度的分布；坐标为大地坐标，单位为米，曲线表示为地表火成岩的等厚度曲线，单位为米；不同的曲线表示不同的厚度。
参见图9，为在同一目标地表高速覆盖地区，采用现有的模型静校正方法、层析静校正方法和本申请提供的地表高速覆盖地区的地震资料静校正方法得到共炮检距初至和叠加剖面，由图9可知，本申请提供的地表高速覆盖地区的地震资料静校正方法得到的共炮检距初至波形的平滑程度更高，叠加剖面同相轴波形的连续性更好，静校正效果更好。
从上述描述可知，本申请提供的地表高速覆盖地区的地震资料静校正方法及装置，能够实现地表高速和下伏低降速带的同时校正，解决地表存在高速火成岩覆盖地区的静校正问题，改善地震资料成像质量。
从硬件层面来说，为了能够实现地表高速和下伏低降速带的同时校正，改善地震资料成像质量，本申请提供一种用于实现所述地表高速覆盖地区的地震资料静校正方法中的全部或部分内容的电子设备的实施例所述电子设备具体包含有如下内容：
处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于实现所述地表高速覆盖地区的地震资料静校正装置以及用户终端等相关设备之间的信息传输；该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照实施例用于实现所述地表高速覆盖地区的地震资料静校正方法的实施例及用于实现所述地表高速覆盖地区的地震资料静校正装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。
图10为本申请实施例的电子设备9600的***构成的示意框图。如图10所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图10是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
在本申请一个或多个实施例中，地震资料静校正功能可以被集成到中央处理器9100中。其中，中央处理器9100可以被配置为进行如下控制：
步骤100：基于预获取的目标地表高速覆盖地区的电阻率剖面获取待测点的电阻率值，并将该电阻率值输入预获取的地表火成岩模型，以获取该待测点对应的火成岩层的实际厚度。
步骤200：应用所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正。
步骤300：根据所述待测点的水平位置坐标，从预获取的近地表结构模型获取对应的下伏低降速层的静校正量，以实现对所述待测点的下伏低降速层对应的地震资料进行低速校正。
从上述描述可知，本申请的实施例提供的电子设备，能够提高配置资管风控引擎的自动化程度，进而提高对投资交易的风险预警的准确性。
在另一个实施方式中，地表高速覆盖地区的地震资料静校正装置可以与中央处理器9100分开配置，例如可以将地表高速覆盖地区的地震资料静校正装置配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现地震资料静校正的配置功能。
如图10所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图10中所示的所有部件；此外，电子设备9600还可以包括图10中没有示出的部件，可以参考现有技术。
如图10所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。
上述描述可知，本申请的实施例提供的电子设备，能够实现地表高速和下伏低降速带的同时校正，改善地震资料成像质量。
本申请的实施例还提供能够实现上述实施例中的地表高速覆盖地区的地震资料静校正方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的地表高速覆盖地区的地震资料静校正方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤100：基于预获取的目标地表高速覆盖地区的电阻率剖面获取待测点的电阻率值，并将该电阻率值输入预获取的地表火成岩模型，以获取该待测点对应的火成岩层的实际厚度。
步骤200：应用所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正。
步骤300：根据所述待测点的水平位置坐标，从预获取的近地表结构模型获取对应的下伏低降速层的静校正量，以实现对所述待测点的下伏低降速层对应的地震资料进行低速校正。
从上述描述可知，本申请实施例提供的计算机可读存储介质，能够实现地表高速和下伏低降速带的同时校正，改善地震资料成像质量。
本申请中上述方法的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。相关之处参见方法实施例的部分说明即可。
本领域内的技术人员应明白，本申请的实施例可提供为方法、***、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本申请中应用了具体实施例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想；同时，对于本领域的一般技术人员，依据本申请的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本申请的限制。
Claims (10)
1.一种地表高速覆盖地区的地震资料静校正方法，其特征在于，包括：
基于预获取的目标地表高速覆盖地区的电阻率剖面获取待测点的电阻率值，并将该电阻率值输入预获取的地表火成岩模型，以获取该待测点对应的火成岩层的实际厚度，所述地表火成岩模型包括：目标地表高速覆盖地区内的电阻率与火成岩层实际厚度对应关系；
应用所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正；
根据所述待测点的水平位置坐标，从预获取的近地表结构模型获取对应的下伏低降速层的静校正量，以实现对所述待测点的下伏低降速层对应的地震资料进行低速校正；
在所述从预获取的近地表结构模型获取对应的下伏低降速层的静校正量之前，还包括：
采集所述目标地表高速覆盖地区内经过高速校正后的地震资料对应的初至时间；
剔除该初至时间中的高速直达波初至时间，以获取所述初至时间中的折射波初至时间；
对所述折射波初至时间进行层析反演，以生成所述近地表结构模型，其中，该近地表结构模型的横坐标表示水平位置坐标，纵坐标表示剔除火成岩层实际厚度后的深度位置坐标，色标为速度。
2.根据权利要求1所述的地表高速覆盖地区的地震资料静校正方法，其特征在于，所述根据所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正，包括：
根据所述火成岩层的实际厚度和预获取的火成岩速度，获取炮、检点高速时差；
应用所述炮、检点高速时差，对所述待测点的火成岩层对应的地震资料进行高速校正。
3.根据权利要求1所述的地表高速覆盖地区的地震资料静校正方法，其特征在于，还包括：
接收目标地表高速覆盖地区内的垂直磁感应场的归一化感应电动势；
基于所述垂直磁感应场的归一化感应电动势获取所述目标地表高速覆盖地区的电阻率剖面。
4.根据权利要求1所述的地表高速覆盖地区的地震资料静校正方法，其特征在于，还包括：
获取火成岩厚度采集点火成岩层的实际厚度，其中，该火成岩厚度采集点在所述预获取的目标地表高速覆盖地区内；
从所述预获取的目标地表高速覆盖地区的电阻率剖面中获取所述火成岩厚度采集点的电阻率；
根据所述火成岩厚度采集点的电阻率及火成岩厚度采集点火成岩层的实际厚度获取所述地表火成岩模型。
5.一种地表高速覆盖地区的地震资料静校正装置，其特征在于，包括：
获取实际厚度模块，用于基于预获取的目标地表高速覆盖地区的电阻率剖面获取待测点的电阻率值，并将该电阻率值输入预获取的地表火成岩模型，以获取该待测点对应的火成岩层的实际厚度，所述地表火成岩模型包括：目标地表高速覆盖地区内的电阻率与火成岩层实际厚度对应关系；
高速校正模块，用于应用所述火成岩层的实际厚度和预获取的火成岩速度，对所述待测点的火成岩层对应的地震资料进行高速校正；
低速校正模块，用于根据所述待测点的水平位置坐标，从预获取的近地表结构模型获取对应的下伏低降速层的静校正量，以实现对所述待测点的下伏低降速层对应的地震资料进行低速校正；
采集初至模块，用于采集所述目标地表高速覆盖地区内经过高速校正后的地震资料对应的初至时间；
获取折射波初至模块，用于剔除该初至时间中的高速直达波初至时间，以获取所述初至时间中的折射波初至时间；
生成近地表结构模型模块，用于对所述折射波初至时间进行层析反演，以生成所述近地表结构模型，其中，该近地表结构模型的横坐标表示水平位置坐标，纵坐标表示剔除火成岩层实际厚度后的深度位置坐标，色标为速度。
6.根据权利要求5所述的地表高速覆盖地区的地震资料静校正装置，其特征在于，所述高速校正模块包括：
获取高速时差单元，用于根据所述火成岩层的实际厚度和预获取的火成岩速度，获取炮、检点高速时差；
高速校正单元，用于应用所述炮、检点高速时差，对所述待测点的火成岩层对应的地震资料进行高速校正。
7.根据权利要求5所述的地表高速覆盖地区的地震资料静校正装置，其特征在于，还包括：
接收二次感应电动势模块，用于接收目标地表高速覆盖地区内的垂直磁感应场的归一化感应电动势；
获取电阻率剖面模块，用于基于所述垂直磁感应场的归一化感应电动势获取所述目标地表高速覆盖地区的电阻率剖面。
8.根据权利要求5所述的地表高速覆盖地区的地震资料静校正装置，其特征在于，还包括：
采集厚度模块，用于获取火成岩厚度采集点火成岩层的实际厚度，其中，该火成岩厚度采集点在所述预获取的目标地表高速覆盖地区内；
获取电阻率模块，用于从所述预获取的目标地表高速覆盖地区的电阻率剖面中获取所述火成岩厚度采集点的电阻率；
获取对应关系模块，用于获取根据所述火成岩厚度采集点的电阻率及火成岩厚度采集点火成岩层的实际厚度获取所述地表火成岩模型。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述的地表高速覆盖地区的地震资料静校正方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述指令被执行时实现权利要求1至4任一项所述的地表高速覆盖地区的地震资料静校正方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911064857.9A CN112764107B (zh) | 2019-11-04 | 2019-11-04 | 地表高速覆盖地区的地震资料静校正方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911064857.9A CN112764107B (zh) | 2019-11-04 | 2019-11-04 | 地表高速覆盖地区的地震资料静校正方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112764107A CN112764107A (zh) | 2021-05-07 |
| CN112764107B true CN112764107B (zh) | 2024-05-28 |
Family
ID=75692383
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911064857.9A Active CN112764107B (zh) | 2019-11-04 | 2019-11-04 | 地表高速覆盖地区的地震资料静校正方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112764107B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2007139977A (ru) * | 2007-10-29 | 2009-05-10 | Закрытое акционерное общество "КРАСНОЯРСКГЕОФИЗИКА" (RU) | Способ определения скоростного строения среды |
| CN103728666A (zh) * | 2012-10-16 | 2014-04-16 | 中国石油化工股份有限公司 | 一种确定静校正中间参考面替换速度的方法 |
| CN105093320A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 针对高速结晶盐壳覆盖区层析静校正初至拾取方法 |
| CN107817525A (zh) * | 2017-09-29 | 2018-03-20 | 中国石油天然气集团公司 | 一种近地表常速填充方法及装置 |
| CN108333627A (zh) * | 2018-01-16 | 2018-07-27 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
| CN109581503A (zh) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及*** |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1128373C (zh) * | 2001-01-21 | 2003-11-19 | 周熙襄 | 用初至波对地震勘探资料进行短波长静校正处理的方法 |
| US20050256648A1 (en) * | 2004-05-11 | 2005-11-17 | West Michael P | Velocity determination of the near-surface layers in the earth using exploration 2D or 3D seismic data |
-
2019
- 2019-11-04 CN CN201911064857.9A patent/CN112764107B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2007139977A (ru) * | 2007-10-29 | 2009-05-10 | Закрытое акционерное общество "КРАСНОЯРСКГЕОФИЗИКА" (RU) | Способ определения скоростного строения среды |
| CN103728666A (zh) * | 2012-10-16 | 2014-04-16 | 中国石油化工股份有限公司 | 一种确定静校正中间参考面替换速度的方法 |
| CN105093320A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 针对高速结晶盐壳覆盖区层析静校正初至拾取方法 |
| CN109581503A (zh) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及*** |
| CN107817525A (zh) * | 2017-09-29 | 2018-03-20 | 中国石油天然气集团公司 | 一种近地表常速填充方法及装置 |
| CN108333627A (zh) * | 2018-01-16 | 2018-07-27 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
Non-Patent Citations (6)
| Title |
|---|
| Seismic exploration pratice in igneous rock region of Fushan sag;Lei Dong;《Nonrecurring Meetings 2011: International Geophysical Conference, Shenzhen, Chnia》;第1339-1347页 * |
| 交互减去法压制多次波技术――以北部湾盆地福山凹陷为例;谷跃民;张进铎;孙维昭;赖江德;张素红;;《石油勘探与开发》;第37卷(第01期);第77-82页 * |
| 正交迭代静校正方法;李章明, 刘贤功, 王顺国, 梁秀文;《石油地球物理勘探》;第32卷(第01期);第124-129页 * |
| 火成岩下伏地层深度预测方法研究——以渤海歧北地区为例;蔡纪琰等;《中国石油学会2017年物探技术研讨会论文集》;第900-903页 * |
| 火成岩发育区静校正技术与应用;李素闪;程元方;刘建红;杨会廷;崔向飞;《石油地球物理勘探》;第52卷(第S2期);第45-49页 * |
| 表层火成岩覆盖地区的表层结构调查;丁建荣;毛凤鸣;郝天珧;《地球物理学进展》;第22卷(第03期);第872-878页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112764107A (zh) | 2021-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Scognamiglio et al. | Fast determination of moment tensors and rupture history: What has been learned from the 6 April 2009 L'Aquila earthquake sequence | |
| CN106873033B (zh) | 一种基于神经网络的多属性地震信息融合裂缝预测的方法 | |
| EA021972B1 (ru) | Множественная инверсия анизотропных параметров для трансверсально-изотропной модели земли с наклонной осью симметрии | |
| CN111164462A (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
| CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
| CN110197318A (zh) | 露天矿山品位模型更新方法 | |
| CN110389380A (zh) | 地震剖面的同相轴自动拾取方法、装置及存储介质 | |
| CN105093318B (zh) | 一种自适应波动方程波场延拓静校正方法 | |
| CN112764107B (zh) | 地表高速覆盖地区的地震资料静校正方法及装置 | |
| CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
| Manea et al. | A high‐resolution site amplification map for Wellington, New Zealand | |
| CN112255679B (zh) | 地震资料绕射深度偏移处理方法及装置 | |
| CN112305588B (zh) | 初至波自动拾取方法及装置 | |
| CN113933901B (zh) | 深度域地震数据拼接方法及装置 | |
| CN113721296A (zh) | 远震数据处理方法及装置 | |
| CN110632650B (zh) | 全波形反演数据的分时窗提取方法及装置 | |
| CN113325476B (zh) | 地震资料静校正方法及装置 | |
| CN113077074A (zh) | 基于储层预测因子的储层预测方法及装置 | |
| CN114578428A (zh) | 一种速度场高效构建方法及装置 | |
| CN112541607A (zh) | 火山岩油藏产能预测方法及装置 | |
| CN114578445A (zh) | 基于重力资料确定断裂位置的方法及装置 | |
| CN112363246B (zh) | 一种火成岩岩性识别方法及装置 | |
| CN112748458B (zh) | 用于现场精细变观的在线实时面元获取方法及装置 | |
| CN114442169B (zh) | 地震资料中近源信号压制方法及装置 | |
| CN114114417B (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 |