CN106054251B - 一种初至波拾取方法及装置 - Google Patents
一种初至波拾取方法及装置 Download PDFInfo
- Publication number
- CN106054251B CN106054251B CN201610443272.8A CN201610443272A CN106054251B CN 106054251 B CN106054251 B CN 106054251B CN 201610443272 A CN201610443272 A CN 201610443272A CN 106054251 B CN106054251 B CN 106054251B
- Authority
- CN
- China
- Prior art keywords
- arrival
- seismic data
- autocorrelation function
- correction
- 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 55
- 238000012937 correction Methods 0.000 claims abstract description 118
- 238000012545 processing Methods 0.000 claims abstract description 27
- 238000005311 autocorrelation function Methods 0.000 claims description 81
- 230000003068 static effect Effects 0.000 claims description 32
- 238000007781 pre-processing Methods 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 16
- 239000000523 sample Substances 0.000 claims description 8
- 238000002203 pretreatment Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 208000011580 syndromic disease Diseases 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001052 transient effect Effects 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
技术领域
本申请涉及地球物理勘探数据处理技术领域,特别涉及一种初至波拾取方法及装置。
背景技术
在地震数据处理中折射波静校正、层析静校正等地震勘探技术都要用到地震初至波信息,初至波拾取的质量在很大程度上影响静校正计算的精度和质量。而静校正的应用效果直接影响地震数据分辨率的提高,如果更精确地拾取初至对于地震数据分辨率的提高非常重要。地震勘探野外采时由于地表条件复杂变化,采集到的地震数据初至波背景噪音强,初至波信息往往被强能量噪音掩盖,尤其是远偏移距的初至,这将在很大程度上降低常规初至波自动拾取初至的精度和效率,初至波人工修改初至的工作量将会大量增加,尤其是对于大数据量低信噪比初至波的拾取,需要花费巨大的人力。
发明人发现现有技术中至少存在如下问题:现有的初至波拾取技术得到的初至波信噪比低。
发明内容
本申请实施例的目的是提供一种初至波拾取方法及装置,以提高初至波的信噪比。
为解决上述技术问题,本申请实施例提供一种初至波拾取方法及装置是这样实现的:
一种初至波拾取方法,包括:
获取采集的地震数据,生成地震数据的炮集记录;
对所述地震数据进行校正处理,得到校正后的地震数据;
基于自相关算法对所述校正后的地震数据进行初至拾取前处理;
确定初始拾取时间,在所述初始拾取时间对所述初至拾取前处理得到的地震数据进行初至波拾取。
优选方案中,所述对地震数据进行校正处理,得到校正后的地震数据,包括:
从所述地震数据中获取炮点静校正量和波点校正量;
根据所述炮点静校正量和检波点校正量,对所述采集的地震数据进行数据校正;
对所述数据校正后的地震数据进行线性校正,将所述线性校正得到的地震数据作为校正后的地震数据。
优选方案中,所述对数据校正后的地震数据进行线性校正,包括:计算样点时间长度,将所述地震数据移动所述样点时间长度。
优选方案中,所述计算样点时间长度具体包括:获取所述地震数据的偏移距,确定预设初至波速度和采样率的比值,利用所述偏移距除以所述比值得到所述样点时间长度。
优选方案中,所述基于自相关算法对所述校正后的地震数据进行初至拾取前处理,具体包括:
选取所述炮集记录中第一炮记录的部分初至波;
确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;
选取所述炮集记录中所述第一炮记录的第一地震道数据;
对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;
根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;
对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果。
优选方案中,所述时间窗长度的取值包括:5个样点~50个样点。
优选方案中,所述第一自相关函数和所述第二自相关函数的自相关长度相同,所述自相关长度等于所述时间窗导长度的2倍加1。
优选方案中,所述预设规则包括:第一自相关函数与所述第三自相关函数的乘积等于第二自相关函数。
优选方案中,所述确定初始拾取时间包括:确定第一初至时间,根据所述校正处理确定校正时间,将所述第一初至时间减去所述校正时间得到初始拾取时间。
优选方案中,所述校正时间包括:静校正时间和线性校正时间的总和。
优选方案中,所述确定第一初至时间采用下述中的至少一种方法:能量比值法、最大能量法。
一种初至波拾取装置,包括:地震数据获取模块、校正模块、初至波拾取前处理模块和初至波拾取模块;其中,
所述地震数据获取模块,用于获取采集的地震数据,生成地震数据的炮集记录;
所述校正模块,用于对所述地震数据进行校正处理,得到校正后的地震数据;
所述初至波拾取前处理模块,用于基于自相关算法对所述校正后的地震数据进行初至拾取前处理;
所述初至波拾取模块,用于确定初始拾取时间,在所述初始拾取时间对所述初至拾取前处理得到的地震数据进行初至波拾取。
优选方案中,所述校正模块包括:校正量获取子模块、静校正子模块和线性校正子模块;其中,
所述校正量获取子模块,用于从所述地震数据中获取炮点静校正量和波点校正量;
所述静校正子模块,用于根据所述炮点静校正量和检波点校正量,对所述采集的地震数据进行数据校正;
所述线性校正子模块,用于对所述数据校正后的地震数据进行线性校正,将所述线性校正得到的地震数据作为校正后的地震数据。
优选方案中,所述初至波拾取前处理模块包括:第一自相关函数确定子模块、第二自相关函数确定子模块、第三自相关函数确定子模块和运算子模块;其中,
所述第一自相关函数确定子模块,用于选取所述炮集记录中第一炮记录的部分初至波,确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;
所述第二自相关函数确定子模块,用于选取所述炮集记录中所述第一炮记录的第一地震道数据,对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;
所述第三自相关函数确定子模块,用于根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;
所述运算子模块,用于对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果。
由以上本申请实施例提供的技术方案可见,本申请实施例提供的一种初至波拾取方法及装置,可以先对采集的地震数据进行校正处理和初至拾取前处理,以提高初至波信噪比和分辨率,再对处理后的地震数据进行初至波拾取,可以提高拾取的初至波的信噪比,从而可以提高初至波的自动拾取效率。进一步地,还可以提高静校正的计算精度,并减少人工修改初至波的工作量。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请初至波拾取方法一个实施例的流程图;
图2是本申请初至波拾取方法实施例中地震数据的示意图;
图3是本申请对图2所示的地震数据进行校正该处理后的结果;
图4是本申请方法实施例中对图2所示的地震数据进行初至波拾取得到的结果;
图5是本申请初至波拾取装置一个实施例的模块图;
图6是本申请装置实施例中校正模块的组成示意图;
图7是本申请装置实施例中初至波拾取前处理模块的组成示意图。
具体实施方式
本申请实施例提供一种初至波拾取方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请初至波拾取方法一个实施例的流程图。如图1所示,所述初至波拾取方法可以包括:
S101:获取采集的地震数据,生成地震数据的炮集记录。
所述地震数据可以包括:地震波的时间和振幅。一个炮产生的地震数据可以生成一个炮集记录。所述获取到的地震数据可以是多个炮产生的地震数据,可以对应有多个炮集记录。
S102:对所述地震数据进行校正处理,得到校正后的地震数据。
具体地,可以包括:从所述地震数据中获取炮点静校正量和波点校正量;根据所述炮点静校正量和检波点校正量,对所述采集的地震数据进行数据校正;对所述数据校正后的地震数据进行线性校正,将所述线性校正得到的地震数据作为校正后的地震数据。
所述炮点静校正量和所述波点校正量可以根据地震数据的道头获取。
在另一个实施方式中,所述炮点静校正量和所述波点校正量的值可以采用高程校正量的值。
在另一个实施例方式中,所述炮点静校正量和所述波点校正量的值还可以通过计算得到。具体地,可以用炮点或者检波点基准面高程除以给定的近地表速度得到时间,所述得到的时间就是所述炮点静校正量或者检波点校正量。
所述对采集的地震数据进行数据校正可以包括:用所述炮点静校正量和所述检波点校正量把所述地震数据校正到基准面上,获得校正到基准面上的地震数据。通过对所述地震数据进行数据校正,可以把整个工区的地震数据校正到同一个基准面上,进行数据校正后初至波可以变得更连续、光滑,从而可以提高初至波的拾取效率。进一步地,数据校正后所有初至波的最小时间和最大时间都在大致相同的范围内,在很大程度上可以不受高程的变化的影响,从而可以提高初至波的拾取效率。
其中,所述对数据校正后的地震数据进行线性校正,可以包括:计算样点时间长度,将所述地震数据移动所述样点时间长度。所述将地震数据移动所述样点时间长度可以包括:将所述地震数据向时间值较小的方向整体移动。
进一步地,所述计算样点时间长度具体可以包括:获取所述地震数据的偏移距,确定预设初至波速度和采样率的比值,利用所述偏移距除以所述比值得到所述样点时间长度。
S103:基于自相关算法对所述校正后的地震数据进行初至拾取前处理。
所述基于自相关算法对所述校正后的地震数据进行初至拾取前处理,具体可以包括:选取所述炮集记录中第一炮记录的部分初至波,确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;选取所述炮集记录中所述第一炮记录的第一地震道数据,对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果。
其中,所述预设规则可以包括:第一自相关函数与所述第三自相关函数的乘积等于第二自相关函数。
所述时间窗长度的取值包括:5个样点~50个样点。所述第一自相关函数和所述第二自相关函数的自相关长度相同,所述自相关长度等于所述时间窗导长度的2倍加1。
例如,可以选取所述炮集记录中第一炮记录的部分初至波,确定时间窗长度为L个样点,可以对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数,可以用C(N)表示,N为0,1,2,3……N-1,其中N为相关长度,所述自相关长度大小一般情况可以为N=2*L+1。可以选取所述炮集记录中所述第一炮记录的第一地震道数据,对所述时间窗长度(即L个样点)内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数,可以用A(N)表示,其中N=2*L+1;
利用下述公式1可以计算第三自相关函数,所述第三自相关函数可以用B(M)表示,所述M的取值可以为0,1,2,3……L+1。
利用公式(1)求得所述第三自相关函数B(M)后,可以将所述第三自相关函数和所述第一地震道数据进行褶积运算。
S104:确定初始拾取时间,在所述初始拾取时间对所述初至拾取前处理得到的地震数据进行初至波拾取。
所述确定初始拾取时间具体可以包括:确定第一初至时间,根据所述校正处理确定校正时间,将所述第一初至时间减去所述校正时间得到初始拾取时间。
所述确定第一初至时间采用下述中的至少一种方法:能量比值法、最大能量法。
所述校正时间可以包括:静校正时间和线性校正时间的总和。
下面结合一个实际的例子来说明本申请的技术方案。首先获取地震数据,图2是本申请初至波拾取方法实施例中地震数据的示意图,其显示的数据是一个三维地震记录中一个排列。图3是本申请对图2所示的地震数据进行校正该处理后的结果。可以明显发现图3中校正处理后的地震数据初至能量差别减少,噪声能量减弱,初至变得更光滑连续。选取其中的第一炮记录的部分初至波,例如,可以选定地震道范围为2180地震道至2220地震道这段初至波(图中箭头所指位置),时间窗选择为400毫秒到700毫秒,对本申请进行初至拾取前处理。对经过所述初至拾取前处理得到的地震数据进行初至波拾取。图4是本申请方法实施例中对图2所示的地震数据进行初至波拾取得到的结果。通过图4可以看出,采用本申请的方法进行拾取的初至波,由于初至波光滑连续,能量差别小,噪声干扰弱,可以提高拾取的初至波的信噪比,保证较好的初至波拾取效果。
上述实施例公开的一种初至波拾取方法,可以先对采集的地震数据进行校正处理和初至拾取前处理,以提高初至波信噪比和分辨率,再对处理后的地震数据进行初至波拾取,可以提高拾取的初至波的信噪比,从而可以提高初至波的自动拾取效率。进一步地,还可以提高静校正的计算精度,并减少人工修改初至波的工作量。
图5是本申请初至波拾取装置一个实施例的模块图。如图5所示,所述初至波拾取装置可以包括:地震数据获取模块501、校正模块502、初至波拾取前处理模块503和初至波拾取模块504;其中,
所述地震数据获取模块501,可以用于获取采集的地震数据,生成地震数据的炮集记录;
所述校正模块502,可以用于对所述地震数据进行校正处理,得到校正后的地震数据;
所述初至波拾取前处理模块503,可以用于基于自相关算法对所述校正后的地震数据进行初至拾取前处理;
所述初至波拾取模块504,可以用于确定初始拾取时间,在所述初始拾取时间对所述初至拾取前处理得到的地震数据进行初至波拾取。
图6是本申请装置实施例中校正模块的组成示意图。如图6所示,所述校正模块502可以包括:校正量获取子模块5021、静校正子模块5022和线性校正子模块5023;其中,
所述校正量获取子模块5021,可以用于从所述地震数据中获取炮点静校正量和波点校正量;
所述静校正子模块5022,可以用于根据所述炮点静校正量和检波点校正量,对所述采集的地震数据进行数据校正;
所述线性校正子模块5023,可以用于对所述数据校正后的地震数据进行线性校正,将所述线性校正得到的地震数据作为校正后的地震数据。
图7是本申请装置实施例中初至波拾取前处理模块的组成示意图。如图7所示,所述初至波拾取前处理模块503可以包括:第一自相关函数确定子模块5031、第二自相关函数确定子模块5032、第三自相关函数确定子模块5033和运算子模块5034;其中,
所述第一自相关函数确定子模块5031,可以用于选取所述炮集记录中第一炮记录的部分初至波,确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;
所述第二自相关函数确定子模块5032,可以用于选取所述炮集记录中所述第一炮记录的第一地震道数据,对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;
所述第三自相关函数确定子模块5033,可以用于根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;
所述运算子模块5034,可以用于对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果。
上述实施例公开的初至波拾取装置与本申请方法实施例相对应,可以实现本申请的方法实施例,并取得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (13)
1.一种初至波拾取方法,其特征在于,包括:
获取采集的地震数据,生成地震数据的炮集记录;
对所述地震数据进行校正处理,得到校正后的地震数据;
基于自相关算法对所述校正后的地震数据进行初至拾取前处理;所述基于自相关算法对所述校正后的地震数据进行初至拾取前处理,具体包括:选取所述炮集记录中第一炮记录的部分初至波;确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;选取所述炮集记录中所述第一炮记录的第一地震道数据;对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果;
确定初始拾取时间,在所述初始拾取时间对所述初至拾取前处理得到的地震数据进行初至波拾取。
2.如权利要求1所述的一种初至波拾取方法,其特征在于,所述对地震数据进行校正处理,得到校正后的地震数据,包括:
从所述地震数据中获取炮点静校正量和波点校正量;
根据所述炮点静校正量和检波点校正量,对所述采集的地震数据进行数据校正;
对所述数据校正后的地震数据进行线性校正,将所述线性校正得到的地震数据作为校正后的地震数据。
3.如权利要求2所述的一种初至波拾取方法,其特征在于,所述对数据校正后的地震数据进行线性校正,包括:计算样点时间长度,将所述地震数据移动所述样点时间长度。
4.如权利要求3所述的一种初至波拾取方法,其特征在于,所述计算样点时间长度具体包括:获取所述地震数据的偏移距,确定预设初至波速度和采样率的比值,利用所述偏移距除以所述比值得到所述样点时间长度。
5.如权利要求1所述的一种初至波拾取方法,其特征在于,所述时间窗长度的取值包括:5个样点~50个样点。
6.如权利要求5所述的一种初至波拾取方法,其特征在于,所述第一自相关函数和所述第二自相关函数的自相关长度相同,所述自相关长度等于所述时间窗导长度的2倍加1。
7.如权利要求1所述的一种初至波拾取方法,其特征在于,所述预设规则包括:第一自相关函数与所述第三自相关函数的乘积等于第二自相关函数。
8.如权利要求1所述的一种初至波拾取方法,其特征在于,所述确定初始拾取时间包括:确定第一初至时间,根据所述校正处理确定校正时间,将所述第一初至时间减去所述校正时间得到初始拾取时间。
9.如权利要求8所述的一种初至波拾取方法,其特征在于,所述校正时间包括:静校正时间和线性校正时间的总和。
10.如权利要求8所述的一种初至波拾取方法,其特征在于,所述确定第一初至时间采用下述中的至少一种方法:能量比值法、最大能量法。
11.一种初至波拾取装置,其特征在于,包括:地震数据获取模块、校正模块、初至波拾取前处理模块和初至波拾取模块;其中,
所述地震数据获取模块,用于获取采集的地震数据,生成地震数据的炮集记录;
所述校正模块,用于对所述地震数据进行校正处理,得到校正后的地震数据;
所述初至波拾取前处理模块,用于基于自相关算法对所述校正后的地震数据进行初至拾取前处理;所述基于自相关算法对所述校正后的地震数据进行初至拾取前处理,具体包括:选取所述炮集记录中第一炮记录的部分初至波;确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;选取所述炮集记录中所述第一炮记录的第一地震道数据;对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果;
所述初至波拾取模块,用于确定初始拾取时间,在所述初始拾取时间对所述初至拾取前处理得到的地震数据进行初至波拾取。
12.如权利要求11所述的一种初至波拾取装置,其特征在于,所述校正模块包括:校正量获取子模块、静校正子模块和线性校正子模块;其中,
所述校正量获取子模块,用于从所述地震数据中获取炮点静校正量和波点校正量;
所述静校正子模块,用于根据所述炮点静校正量和检波点校正量,对所述采集的地震数据进行数据校正;
所述线性校正子模块,用于对所述数据校正后的地震数据进行线性校正,将所述线性校正得到的地震数据作为校正后的地震数据。
13.如权利要求11所述的一种初至波拾取装置,其特征在于,所述初至波拾取前处理模块包括:第一自相关函数确定子模块、第二自相关函数确定子模块、第三自相关函数确定子模块和运算子模块;其中,
所述第一自相关函数确定子模块,用于选取所述炮集记录中第一炮记录的部分初至波,确定时间窗长度,对所述时间窗长度内的所述部分初至波对应的地震数据进行自相关运算,确定与所述部分初至波对应的第一自相关函数;
所述第二自相关函数确定子模块,用于选取所述炮集记录中所述第一炮记录的第一地震道数据,对所述时间窗长度内的所述地震道数据进行自相关运算,得到与所述地震道数据对应的第二自相关函数;
所述第三自相关函数确定子模块,用于根据预设规则,利用第一自相关函数和第二自相关函数求取第三自相关函数;
所述运算子模块,用于对所述第三自相关函数和所述第一地震道数据进行褶积运算,将运算得到的结果作为初至拾取前处理的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610443272.8A CN106054251B (zh) | 2016-06-20 | 2016-06-20 | 一种初至波拾取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610443272.8A CN106054251B (zh) | 2016-06-20 | 2016-06-20 | 一种初至波拾取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106054251A CN106054251A (zh) | 2016-10-26 |
CN106054251B true CN106054251B (zh) | 2018-08-10 |
Family
ID=57168496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610443272.8A Active CN106054251B (zh) | 2016-06-20 | 2016-06-20 | 一种初至波拾取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106054251B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106855640B (zh) * | 2017-02-03 | 2019-05-07 | 中国石油天然气集团公司 | 一种海底电缆地震数据初至计算方法及装置 |
CN109143329B (zh) * | 2017-06-27 | 2021-02-19 | 中国石油化工股份有限公司 | 初至拾取处理方法 |
CN109143364B (zh) * | 2017-06-27 | 2020-04-14 | 中国石油化工股份有限公司 | 初至时间拾取方法 |
CN109143367B (zh) * | 2017-06-28 | 2020-06-19 | 中国石油化工股份有限公司 | 一种初至波自动拾取方法 |
CN108279436B (zh) * | 2017-12-20 | 2019-09-10 | 中国石油天然气股份有限公司 | 一种拾取地震波初至的方法及装置 |
CN108387936B (zh) * | 2018-03-20 | 2020-01-10 | 中国石油集团东方地球物理勘探有限责任公司 | 一种可控震源资料初至高效拾取方法及层析静校正方法 |
CN109725345B (zh) * | 2018-11-15 | 2020-08-11 | 中国石油天然气集团有限公司 | 一种初至波正演模拟方法及装置 |
CN112305593B (zh) * | 2019-07-23 | 2023-07-25 | 中国石油天然气集团有限公司 | 异常初至波识别方法及装置 |
CN111077569B (zh) * | 2019-12-23 | 2022-05-06 | 中国石油天然气股份有限公司 | 全波形反演中分时窗提取数据的方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073067A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气集团公司 | 一种提高地震数据初至波自动拾取效率的方法 |
CN104570076A (zh) * | 2013-10-17 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于二分法的地震波初至自动拾取方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6850462B2 (en) * | 2002-02-19 | 2005-02-01 | Probe Technology Services, Inc. | Memory cement bond logging apparatus and method |
US7706211B2 (en) * | 2006-02-06 | 2010-04-27 | Shell Oil Company | Method of determining a seismic velocity profile |
-
2016
- 2016-06-20 CN CN201610443272.8A patent/CN106054251B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073067A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气集团公司 | 一种提高地震数据初至波自动拾取效率的方法 |
CN104570076A (zh) * | 2013-10-17 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于二分法的地震波初至自动拾取方法 |
Non-Patent Citations (1)
Title |
---|
复杂地表静校正方法研究及应用;熊定钰;《中国博士论文全文数据库 基础科学辑》;20130115(第1期);第53页倒数第1段-第54页倒数第2段 * |
Also Published As
Publication number | Publication date |
---|---|
CN106054251A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106054251B (zh) | 一种初至波拾取方法及装置 | |
CN108279436B (zh) | 一种拾取地震波初至的方法及装置 | |
CN107563297B (zh) | 一种波形匹配方法及装置 | |
CN107765298B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN107193045B (zh) | 一种地震数据处理方法及装置 | |
CN106842304B (zh) | 一种叠前深度偏移方法及装置 | |
CN106646611B (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN108897040B (zh) | 一种确定陆检与海底之间的耦合性的方法及装置 | |
CN106338763A (zh) | 一种确定地震剖面显示数据的方法及装置 | |
CN107329170B (zh) | 一种地震波形匹配方法及装置 | |
CN108181644A (zh) | 一种采样方法及装置 | |
CN108035708B (zh) | 一种去除地层界面反射波的方法及装置 | |
CN106353816B (zh) | 一种地震采集脚印噪音压制方法及系统 | |
CN111239832B (zh) | 3d3c-vsp成像处理方法、装置及设备 | |
CN106291691A (zh) | 一种地震偏移成像方法及装置 | |
CN107561588B (zh) | 一种地震数据噪声压制方法及装置 | |
CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
CN109143338A (zh) | 一种地震数据合并的方法及装置 | |
CN111596358B (zh) | 一种多次波的压制方法、设备以及系统 | |
CN106772596B (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN108107470B (zh) | 一种地震数据处理方法及装置 | |
CN106054253A (zh) | 一种道集同相轴校正的方法及装置 | |
CN108983289B (zh) | 一种确定地震波走时的方法及装置 | |
CN111751881A (zh) | 一种海上采集地震数据旅行时的校正方法、装置及系统 | |
CN106125137B (zh) | 一种转换波地震数据处理方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |