CN112241021B - 初至拾取方法及装置 - Google Patents
初至拾取方法及装置 Download PDFInfo
- Publication number
- CN112241021B CN112241021B CN201910639172.6A CN201910639172A CN112241021B CN 112241021 B CN112241021 B CN 112241021B CN 201910639172 A CN201910639172 A CN 201910639172A CN 112241021 B CN112241021 B CN 112241021B
- Authority
- CN
- China
- Prior art keywords
- offset
- range
- fitting
- arrival
- result
- 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 72
- 238000012545 processing Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 13
- 238000003708 edge detection Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 29
- 230000006870 function Effects 0.000 description 9
- 238000001514 detection method Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000005728 strengthening 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- 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/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/41—Arrival times, e.g. of P or S wave or first break
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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)
- Image Analysis (AREA)
Abstract
本发明提供了一种初至拾取方法及装置,该方法包括:获取目标单炮的初至时间和炮检距;确定目标单炮的分界点炮检距,并将第一炮检距至分界点炮检距间的范围作为第一范围,将分界点炮检距至第二炮检距间的范围作为第二范围;根据最小二乘法对第一范围内的目标单炮的炮检距和初至时间进行拟合,得到第一拟合结果,根据最小二乘法对第二范围内的目标单炮的炮检距和初至时间进行拟合,得到第二拟合结果;根据第一拟合结果和第二拟合结果确定目标炮检距的目标初至时间;根据目标初至时间进行分析拾取,得到初至拾取结果。本发明可以得到拾取成功率更高的拾取结果,提高了初至拾取的效率。
Description
技术领域
本发明涉及地震资料处理技术领域,尤其是涉及一种初至拾取方法及装置。
背景技术
在勘探区域的地表以及相应的地质情况越来越复杂的趋势下,初至拾取已是地震资料处理过程中必不可少的环节,初至不仅可以用于时间偏移处理中的静校正计算,同时还是反演深度偏移近地表模型的关键。因此为了做好地震资料的静校正处理工作和深度域精细速度模型建立工作,必须要对初至时间进行拾取,越精准的初至拾取越能为后续处理提供保障。
目前地震资料多是高密度三维数据,具有炮数多、排列长、道数多的特点,但现有软件普遍存在自动拾取成功率低的问题,必须要对初至进行人工的手动修改,从而严重耗时,影响项目处理周期,因此急切需要找到一种更高效的方法来解决这一问题。
发明内容
本发明提供了一种初至拾取方法及装置,可以提高拾取的准确率和成功率,进而提高拾取效率,减少初至拾取工作的整体时间。
第一方面,本发明实施例提供了一种初至拾取方法,该方法包括:确定所述目标单炮的分界点炮检距,并将第一炮检距至所述分界点炮检距间的范围作为第一范围,将所述分界点炮检距至第二炮检距间的范围作为第二范围;所述分界点炮检距介于第二炮检距和第一炮检距之间;所述第一炮检距小于所述第二炮检距;根据最小二乘法对所述第一范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第一拟合结果,根据最小二乘法对所述第二范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第二拟合结果;根据所述第一拟合结果和所述第二拟合结果确定目标炮检距的目标初至时间;根据所述目标初至时间进行分析拾取,得到初至拾取结果。
第二方面,本发明实施例还提供一种初至拾取装置,该装置包括:获取模块,用于获取目标单炮的初至时间和炮检距;分界模块,用于确定所述目标单炮的分界点炮检距,并将第一炮检距至所述分界点炮检距间的范围作为第一范围,将所述分界点炮检距至第二炮检距间的范围作为第二范围;所述分界点炮检距介于第二炮检距和第一炮检距之间;所述第一炮检距小于所述第二炮检距;拟合模块,用于根据最小二乘法对所述第一范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第一拟合结果,根据最小二乘法对所述第二范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第二拟合结果;确定模块,用于根据所述第一拟合结果和所述第二拟合结果确定目标炮检距的目标初至时间;拾取模块,用于根据所述目标初至时间进行分析拾取,得到初至拾取结果。
第三方面,本发明实施例还提供一种计算机设备,包括存储器、处理器,所述存储器中存储有可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述初至拾取方法。
第四方面,本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质,所述程序代码使所述处理器执行上述初至拾取方法。
本发明实施例带来了以下有益效果:本发明实施例提供了一种初至拾取方法及装置,获取目标单炮的初至时间和炮检距之后,该方法首先确定目标单炮的分界点炮检距,并将第一炮检距至分界点炮检距间的范围作为第一范围,将分界点炮检距至第二炮检距间的范围作为第二范围,进而,根据最小二乘法对第一范围内的目标单炮的炮检距和初至时间进行拟合,得到第一拟合结果,根据最小二乘法对第二范围内的目标单炮的炮检距和初至时间进行拟合,得到第二拟合结果,根据第一拟合结果和第二拟合结果确定目标炮检距的目标初至时间,从而对得到的目标初至时间进行分析拾取,得到初至拾取结果。本发明实施例通过最小二乘法对初至时间和炮检距进行分段拟合,可提高拟合的准确度,从而基于拟合结果得到准确度更高的目标初至时间,进而得到拾取成功率更高的拾取结果,提高了初至拾取的效率。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的初至拾取方法流程图;
图2为本发明实施例提供的初至点两边的振幅值差异示意图;
图3为本发明实施例提供的未分段最小二乘公式拟合结果;
图4为本发明实施例提供的第一种范围的近炮检距初至拟合结果示意图;
图5为本发明实施例提供的第二种范围的近炮检距初至拟合结果示意图;
图6为本发明实施例提供的第三种范围的近炮检距初至拟合结果示意图;
图7为本发明实施例提供的近炮检距范围的多个相关系数随偏移距(即炮检距)变化趋势;
图8为本发明实施例提供的第一种范围的远炮检距初至拟合结果示意图;
图9为本发明实施例提供的第二种范围的远炮检距初至拟合结果示意图;
图10为本发明实施例提供的第三种范围的远炮检距初至拟合结果示意图;
图11为本发明实施例提供的远炮检距范围的多个相关系数随偏移距(即炮检距)变化趋势;
图12为本发明实施例提供的相关系数R随炮检距变化趋势示意图;
图13为本发明实施例提供的单炮上拟合约束前后的初至自动拾取结果示意图;
图14为本发明实施例提供的全工区内拟合约束前后的初至自动拾取结果示意图;
图15为本发明实施例提供的初至拾取装置结构框图;
图16为本发明实施例提供的计算机设备结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,现有初至波自动拾取软件的基本原理多是基于地震时窗属性特征的初至拾取方法,主要通过图像处理中的“边界检测”确定出边界,初步确定初至时间;然后在此基础上,通过对相邻能量边界进行扫描,剔除异常点,进而确定初至时间。其关键点在于准确的判定边界位置。以下公式是单道能量边界检测的公式。其中,fi为第i个样点求取的边界检测因子;Sk为第k个样点的振幅值;n为计算的点数。从以上公式中可以看出,判定初至时间主要依赖初至点两边样点的振幅值差异,振幅值差异越大越容易判定,参见图2所示的初至点两边的振幅值差异示意图。然而在实际资料中,地震资料的采集会受到多种因素的影响,会导致很多情况下初至能量较弱,且噪声干扰严重,从而造成初至拾取成功率低的问题。如果沿初至起跳点定义增益时窗,增强起跳点及其之后的能量差异,则可大大提高现有软件的初至拾取效率。所以如何具有针对性的对初至部分的能量进行调整,就成为一个关键的问题。
基于此,本发明实施例提供的一种初至拾取方法及装置可以提高初至拾取的效率。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种初至拾取方法进行详细介绍。
本发明实施例提供了一种初至拾取方法,参见图1所示的一种初至拾取方法流程图,该方法包括以下步骤:
步骤S102,获取目标单炮的初至时间和炮检距。
在本发明实施例中,目标单炮是通过调查分析全工区的单炮初至情况,选取的有代表性的并能够控制全工区变化的单炮。可以预先对目标单炮进行手工精细的初至拾取,得到初至时间和炮检距。
步骤S104,确定目标单炮的分界点炮检距,并将第一炮检距至分界点炮检距间的范围作为第一范围,将分界点炮检距至第二炮检距间的范围作为第二范围。
在本发明实施例中,目标单炮的分界点炮检距是介于第二炮检距和第一炮检距之间的炮检距。第一炮检距可以是目标单炮的最小炮检距,第二炮检距可以是目标单炮的最大炮检距。通过分界点炮检距将目标单炮的炮检距分成两个范围。
步骤S106,根据最小二乘法对第一范围内的目标单炮的炮检距和初至时间进行拟合,得到第一拟合结果,根据最小二乘法对第二范围内的目标单炮的炮检距和初至时间进行拟合,得到第二拟合结果。
在本发明实施例中,炮检距与初至时间一一对应,因此,可以对第一范围内的炮检距及其对应的初至时间进行拟合,得到第一拟合结果,第一拟合结果是时间关于炮检距的函数,炮检距在第一范围内取值。对第二范围内的炮检距及其对应的初至时间进行拟合,得到第二拟合结果,第二拟合结果是时间关于炮检距的函数,炮检距在第二范围内取值。
步骤S108,根据第一拟合结果和第二拟合结果确定目标炮检距的目标初至时间。
确定目标炮检距所属于的炮检距范围区段,如果目标炮检距属于第一拟合结果所对应的炮检距范围区段,即第一范围内,则使用第一拟合结果计算该目标炮检距的目标初至时间,如果目标炮检距属于第二拟合结果所对应的炮检距范围区段,即第二范围内,则使用第二拟合结果计算该目标炮检距的目标初至时间。从而,可以对每一个给定的目标炮检距得到一个对应的目标初至时间。
步骤S110,根据目标初至时间进行分析拾取,得到初至拾取结果。
在本发明实施例中,可以使用现有方法根据目标初至时间进行分析拾取,得到初至拾取结果。
本发明实施例提供了一种初至拾取方法,获取目标单炮的初至时间和炮检距之后,该方法首先确定目标单炮的分界点炮检距,并将第一炮检距至分界点炮检距间的范围作为第一范围,将分界点炮检距至第二炮检距间的范围作为第二范围,进而,根据最小二乘法对第一范围内的目标单炮的炮检距和初至时间进行拟合,得到第一拟合结果,根据最小二乘法对第二范围内的目标单炮的炮检距和初至时间进行拟合,得到第二拟合结果,根据第一拟合结果和第二拟合结果确定目标炮检距的目标初至时间,从而对得到的目标初至时间进行分析拾取,得到初至拾取结果。本发明实施例通过最小二乘法进行分段拟合,可以提高拟合的准确度,从而基于拟合结果得到准确度更高的目标初至时间,进而得到拾取成功率更高的拾取结果,提高了初至拾取的效率。
参见图3所示的未分段最小二乘公式拟合结果,图中黑色细曲线为拟合结果,黑色细曲线周围的散点粗曲线为初至分布情况,可见,曲线两端已经与该部分的初至数据存在一定偏差。考虑到为了减小拟合的误差,确定目标单炮的分界点炮检距,可以按照如下步骤执行:
(1)对目标单炮的炮检距进行分段得到第一组炮检距范围和第二组炮检距范围;第一组炮检距范围包括多个炮检距范围,且每个炮检距范围的起点为第一炮检距;第二组炮检距范围包括多个炮检距范围,且每个炮检距范围的终点为第二炮检距。
在本发明实施例中,第一组炮检距范围和第二组炮检距范围中均包括多个炮检距范围。例如,第一组炮检距范围包括的炮检距范围可以为:0~2000m,0~2500m,0~3000m,0~3500m,0~6000m等;第二组炮检距范围包括的炮检距范围可以为:5000~8000m,4000~8000m,3000~8000m,2500~8000m,2000~8000m等,其中,炮检距8000m处为第二炮检距,炮检距0m处为第一炮检距。
需要说明的是,第一炮检距、第二炮检距、第一组炮检距范围中每个炮检距范围的终点以及第二组炮检距范围中每个炮检距范围的起点的取值均可以根据实际需求进行设置,本发明实施例对此不作具体限定。
(2)根据最小二乘法对第一组炮检距范围和第二组炮检距范围的每个炮检距范围内的初至时间和炮检距进行拟合,得到多个拟合结果。
在本发明实施例中,对每个炮检距范围分别进行拟合,对每个拟合结果均获取其相关系数,可以得到与炮检距范围一一对应的多个相关系数。例如,对0~2000m,0~2500m,0~3000m,0~3500m,0~6000m以及5000~8000m,4000~8000m,3000~8000m,2500~8000m,2000~8000m这些炮检距范围可以分别利用最小二乘法对初至时间和炮检距进行拟合,得到多个拟合结果。
(3)获取每个拟合结果的相关系数,并根据相关系数确定目标单炮的分界点炮检距。
在本发明实施例中,根据相关系数确定分界点炮检距可以包括以下步骤:确定第一范围的多个相关系数随炮检距变化第一曲线,并确定第二范围的多个相关系数随炮检距变化第二曲线;将第一曲线和第二曲线交点的炮检距作为分界点炮检距。
为使拟合误差最小,处理中将数据进行不同炮检距拟合,分界点炮检距需通过实验来确定。图4-图11展示了对不同范围炮检距的初至拟合结果,通过相关系数随炮检距变化的曲线可以看出,随炮检距范围的增大,相关系数逐渐趋近于1。
参见图4-图7所示的不同范围的近炮检距初至拟合结果示意图,图4示出了0~2000m炮检距范围的拟合曲线,图5示出了0~2500m炮检距范围的拟合曲线,图6示出了0~3000m炮检距范围的拟合曲线,图7示出了近炮检距范围的多个相关系数随偏移距(即炮检距)变化趋势。参见图8-图11所示的不同范围的远炮检距初至拟合结果示意图,图中8示出了5000~8000m炮检距范围的拟合曲线,图9示出了4000~8000m炮检距范围的拟合曲线,图10示出了3000~8000m炮检距范围的拟合曲线,图11示出了远炮检距范围的多个相关系数随偏移距(即炮检距)变化趋势。
参见图12所示的相关系数R随炮检距变化趋势示意图,将两部分拟合结果的相关系数曲线放在一起,可见在炮检距大小为3000m附近时,近、远炮检距的拟合相关系数达到一个均衡点。在这一点处,近炮检距和远炮检距范围都具有较好的拟合效果。因此选定3000m作为近、远炮检距的分界点炮检距。
本发明实施例将距离起点炮检距不同远近范围的炮检距与起点炮检距分别组成多个炮检距范围,将距离终点炮检距不同远近范围的炮检距与终点炮检距分别组成多个炮检距范围,再获取各个炮检距范围内拟合结果的相关系数,根据相关系确定分界点炮检距,从而对第一炮检距范围和第二炮检距范围分别进行拟合,相比于直接对起点到终点炮检距这个范围的炮检距与初至时间直接进行拟合准确度更高。
为了得到更精确的初至拾取结果,根据目标初至时间进行分析拾取,得到初至拾取结果可以按照如下步骤执行:
对目标初至时间进行增益处理,得到增益结果;根据边缘检测算法对增益结果进行分析拾取,得到初至拾取结果。
用目标初至时间作为约束起始点,给定增益时窗对目标初至时间进行增益处理来加强初至起跳点能量,进而加强初至起跳点两侧能量差异。利用边缘检测方法对全区范围内单炮增益结果进行自动分析拾取,得到相对更加精准的初至拾取结果,从而提高初至拾取的成功率。
参见图13所示的单炮上拟合约束前后的初至自动拾取结果示意图,该图中上半图为应用本方法之前的初至自动拾取情况,图中圈出的散点即为拾取结果,存在多个跳点以及不能正确识别初至时间的点,圆圈标识了自动拾取较差的区域,这些必须要进行手动修改。该图中下图是应用本发明实施例提供的方法后的初至自动拾取情况,可以明显的看到初至拾取成功率大大提高,大部分区域不需再进行手工修改。
参见图14所示的全工区内拟合约束前后的初至自动拾取结果示意图,图中,上半图是使用本发明实施例方法之前的初至自动拾取结果示意图,下半图是使用本发明实施例方法之后的初至自动拾取结果示意图。可见,使用本发明实施例方法的初至自动拾取结果明显优于现有方法的结果,其收敛性更好,随机跳点明显减少,同时效率也大大提高,应用本发明实施例的方法既可以保证质量,还可以提高大约5倍效率。
图13-图14是某盆地某三维地震资料使用该方法后的实际效果图。从该项目的应用效果看,使用初至拟合法约束后的自动拾取成功率明显提高,在投入同等资源情况下,总拾取时间由50天缩短至10天,提高至原来的5倍,解决了该地区普遍存在的初至拾取效率低问题,粗略计算在人工费和设备费上节约5万余元。该方法的提出为提高拾取的质量,缩短处理周期,减少人工耗时及费用作出了巨大贡献,同样适用于其他沙漠地区或者地表情况相对稳定地区的初至拾取,具有良好的推广应用前景。
本发明实施例还提供一种初至拾取装置,参见图15所示的初至拾取装置结构框图,该装置包括:获取模块71,用于获取目标单炮的初至时间和炮检距;分界模块72,用于确定目标单炮的分界点炮检距,并将第一炮检距至分界点炮检距间的范围作为第一范围,将分界点炮检距至第二炮检距间的范围作为第二范围;分界点炮检距介于第二炮检距和第一炮检距之间;第一炮检距小于第二炮检距;拟合模块73,用于根据最小二乘法对第一范围内的目标单炮的炮检距和初至时间进行拟合,得到第一拟合结果,根据最小二乘法对第二范围内的目标单炮的炮检距和初至时间进行拟合,得到第二拟合结果;确定模块74,用于根据第一拟合结果和第二拟合结果确定目标炮检距的目标初至时间;拾取模块75,用于根据目标初至时间进行分析拾取,得到初至拾取结果。
分界模块,具体用于:对目标单炮的炮检距进行分段得到第一组炮检距范围和第二组炮检距范围;第一组炮检距范围包括多个炮检距范围,且每个炮检距范围的起点为第一炮检距;第二组炮检距范围包括多个炮检距范围,且每个炮检距范围的终点为第二炮检距;根据最小二乘法对第一组炮检距范围和第二组炮检距范围的每个炮检距范围内的初至时间和炮检距进行拟合,得到多个拟合结果;获取每个拟合结果的相关系数,并根据相关系数确定目标单炮的分界点炮检距。
分界模块,具体用于:确定第一范围的多个相关系数随炮检距变化第一曲线,并确定第二范围的多个相关系数随炮检距变化第二曲线;将第一曲线和第二曲线交点的炮检距作为分界点炮检距。
拾取模块,具体用于:对目标初至时间进行增益处理,得到增益结果;根据边缘检测算法对增益结果进行分析拾取,得到初至拾取结果。
本发明实施例还提供一种计算机设备,参见图16所示的计算机设备结构示意框图,该计算机设备包括存储器81、处理器82,存储器中存储有可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述任一种方法的步骤。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的计算机设备的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述
本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质,程序代码使处理器执行上述任一种方法的步骤。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种初至拾取方法,其特征在于,包括:
获取目标单炮的初至时间和炮检距;
确定所述目标单炮的分界点炮检距,并将第一炮检距至所述分界点炮检距间的范围作为第一范围,将所述分界点炮检距至第二炮检距间的范围作为第二范围;所述分界点炮检距介于第二炮检距和第一炮检距之间;所述第一炮检距小于所述第二炮检距;
根据最小二乘法对所述第一范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第一拟合结果,根据最小二乘法对所述第二范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第二拟合结果;
根据所述第一拟合结果和所述第二拟合结果确定目标炮检距的目标初至时间;
根据所述目标初至时间进行分析拾取,得到初至拾取结果。
2.根据权利要求1所述的方法,其特征在于确定所述目标单炮的分界点炮检距,包括:
对所述目标单炮的炮检距进行分段得到第一组炮检距范围和第二组炮检距范围;所述第一组炮检距范围包括多个炮检距范围,且每个炮检距范围的起点为所述第一炮检距;所述第二组炮检距范围包括多个炮检距范围,且每个炮检距范围的终点为第二炮检距;
根据最小二乘法对所述第一组炮检距范围和所述第二组炮检距范围的每个炮检距范围内的初至时间和炮检距进行拟合,得到多个拟合结果;
获取每个所述拟合结果的相关系数,并根据所述相关系数确定目标单炮的分界点炮检距。
3.根据权利要求2所述的方法,其特征在于,根据所述相关系数确定目标单炮的分界点炮检距,包括:
确定所述第一范围的多个相关系数随炮检距变化第一曲线,并确定所述第二范围的多个相关系数随炮检距变化第二曲线;
将所述第一曲线和所述第二曲线交点的炮检距作为分界点炮检距。
4.根据权利要求1所述的方法,其特征在于,包括:根据所述目标初至时间进行分析拾取,得到初至拾取结果,包括:
对所述目标初至时间进行增益处理,得到增益结果;
根据边缘检测算法对所述增益结果进行分析拾取,得到初至拾取结果。
5.一种初至拾取装置,其特征在于,包括:
获取模块,用于获取目标单炮的初至时间和炮检距;
分界模块,用于确定所述目标单炮的分界点炮检距,并将第一炮检距至所述分界点炮检距间的范围作为第一范围,将所述分界点炮检距至第二炮检距间的范围作为第二范围;所述分界点炮检距介于第二炮检距和第一炮检距之间;所述第一炮检距小于所述第二炮检距;
拟合模块,用于根据最小二乘法对所述第一范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第一拟合结果,根据最小二乘法对所述第二范围内的所述目标单炮的炮检距和所述初至时间进行拟合,得到第二拟合结果;
确定模块,用于根据所述第一拟合结果和所述第二拟合结果确定目标炮检距的目标初至时间;
拾取模块,用于根据所述目标初至时间进行分析拾取,得到初至拾取结果。
6.根据权利要求5所述的装置,其特征在于,所述分界模块,具体用于:
对所述目标单炮的炮检距进行分段得到第一组炮检距范围和第二组炮检距范围;所述第一组炮检距范围包括多个炮检距范围,且每个炮检距范围的起点为所述第一炮检距;所述第二组炮检距范围包括多个炮检距范围,且每个炮检距范围的终点为第二炮检距;
根据最小二乘法对所述第一组炮检距范围和所述第二组炮检距范围的每个炮检距范围内的初至时间和炮检距进行拟合,得到多个拟合结果;
获取每个所述拟合结果的相关系数,并根据所述相关系数确定目标单炮的分界点炮检距。
7.根据权利要求6所述的装置,其特征在于,所述分界模块,具体用于:
确定所述第一范围的多个相关系数随炮检距变化第一曲线,并确定所述第二范围的多个相关系数随炮检距变化第二曲线;
将所述第一曲线和所述第二曲线交点的炮检距作为分界点炮检距。
8.根据权利要求5所述的装置,其特征在于,所述拾取模块,具体用于:
对所述目标初至时间进行增益处理,得到增益结果;
根据边缘检测算法对所述增益结果进行分析拾取,得到初至拾取结果。
9.一种计算机设备,包括存储器、处理器,所述存储器中存储有可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至4任一项所述的方法的步骤。
10.一种具有处理器可执行的非易失的程序代码的计算机可读介质,其特征在于,所述程序代码使所述处理器执行上述权利要求1-4任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910639172.6A CN112241021B (zh) | 2019-07-16 | 2019-07-16 | 初至拾取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910639172.6A CN112241021B (zh) | 2019-07-16 | 2019-07-16 | 初至拾取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112241021A CN112241021A (zh) | 2021-01-19 |
CN112241021B true CN112241021B (zh) | 2023-05-26 |
Family
ID=74166487
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910639172.6A Active CN112241021B (zh) | 2019-07-16 | 2019-07-16 | 初至拾取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112241021B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101178439A (zh) * | 2007-03-22 | 2008-05-14 | 孙素梅 | 精确地表调查方法 |
CN101551464A (zh) * | 2008-03-31 | 2009-10-07 | 中国石油天然气集团公司 | 一种复杂地表三维勘探初至波时间确定方法 |
CN101630015A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高初至波拾取精度和效率的方法 |
CN101930080A (zh) * | 2009-06-26 | 2010-12-29 | 中国石油集团东方地球物理勘探有限责任公司 | 曲面拟合海底电缆二次定位方法 |
CN104570122A (zh) * | 2015-01-14 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | 基于反射波的地震数据静校正方法及装置 |
CN105242311A (zh) * | 2015-10-22 | 2016-01-13 | 中国石油天然气集团公司 | 预估初至拾取的时窗的方法和装置 |
CN108897036A (zh) * | 2018-05-15 | 2018-11-27 | 中国石油天然气集团有限公司 | 一种地震数据处理方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1128373C (zh) * | 2001-01-21 | 2003-11-19 | 周熙襄 | 用初至波对地震勘探资料进行短波长静校正处理的方法 |
US8995222B2 (en) * | 2010-05-06 | 2015-03-31 | Bp Corporation North America Inc. | System and method for accurate determination of ocean bottom seismometer positioning and timing |
US10920585B2 (en) * | 2017-12-26 | 2021-02-16 | Saudi Arabian Oil Company | Determining sand-dune velocity variations |
-
2019
- 2019-07-16 CN CN201910639172.6A patent/CN112241021B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101178439A (zh) * | 2007-03-22 | 2008-05-14 | 孙素梅 | 精确地表调查方法 |
CN101551464A (zh) * | 2008-03-31 | 2009-10-07 | 中国石油天然气集团公司 | 一种复杂地表三维勘探初至波时间确定方法 |
CN101630015A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高初至波拾取精度和效率的方法 |
CN101930080A (zh) * | 2009-06-26 | 2010-12-29 | 中国石油集团东方地球物理勘探有限责任公司 | 曲面拟合海底电缆二次定位方法 |
CN104570122A (zh) * | 2015-01-14 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | 基于反射波的地震数据静校正方法及装置 |
CN105242311A (zh) * | 2015-10-22 | 2016-01-13 | 中国石油天然气集团公司 | 预估初至拾取的时窗的方法和装置 |
CN108897036A (zh) * | 2018-05-15 | 2018-11-27 | 中国石油天然气集团有限公司 | 一种地震数据处理方法及装置 |
Non-Patent Citations (5)
Title |
---|
RTM scattered artifacts removal by illumination compensation method – application in streamer data;Wang Wei et al.;《2016 SEG High Performance Computing Workshop》;全文 * |
The Use of Wavelet-Based Denoising Techniques to Enhance the First-Arrival Picking on Seismic Traces;Said Gaci;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;第52卷(第8期);全文 * |
基于波形相似性的远震初至拾取方法进展与对比研究;林凡生等;《地球物理学进展》;第32卷(第4期);全文 * |
最小二乘曲线拟合的微地震初至优化拾取方法及应用;刘腾蛟等;《石油地球物理勘探》;第53卷;全文 * |
结合CEEMDAN和主成分分析的低信噪比微地震初至信号检测;胡瑞卿等;《石油地球物理勘探》;第54卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112241021A (zh) | 2021-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN104330828B (zh) | 甜点储层的预测方法及预测装置 | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN106443774A (zh) | 一种提高起伏地表叠前深度偏移成像精度的方法 | |
CN109884700B (zh) | 多信息融合地震速度建模方法 | |
CN109187744B (zh) | 一种基于静力触探试验的剪切波速估算方法 | |
CN109239781A (zh) | 一种地震数据校正方法及装置 | |
CN103954997A (zh) | 一种地震勘探资料近地表静校正处理方法 | |
CN112241021B (zh) | 初至拾取方法及装置 | |
CN111257935B (zh) | 一种加快层析反演速度收敛的速度融合方法及处理终端 | |
CN109490961B (zh) | 起伏地表无射线追踪回折波层析成像方法 | |
CN109725354B (zh) | 各向异性速度建模方法及系统 | |
CN111399037B (zh) | 高速顶界面提取的方法和装置 | |
CN115081485B (zh) | 一种基于ai的漏磁内检测数据自动分析方法 | |
CN106842325B (zh) | 利用检波点静校正迭代法提高初至波拾取效率及准确率的方法 | |
CN112415601B (zh) | 表层品质因子q值的确定方法及装置 | |
CN115079257A (zh) | 基于融合网络的q值估计及地震衰减补偿方法 | |
CN111352157B (zh) | 一种横波静校正方法及系统 | |
CN112394413A (zh) | 三维初至波剩余静校正方法及装置 | |
CN115327630A (zh) | 声波测井校正方法、装置、设备、存储介质及程序产品 | |
CN112363222A (zh) | 叠后自适应宽带约束波阻抗反演方法及装置 | |
CN112444861B (zh) | 速度模型更新方法、计算机存储介质及计算机系统 | |
CN111624662A (zh) | 速度场校正方法、装置和设备 | |
CN117872475A (zh) | 用于时延勘探资料的静校正处理方法、装置、设备及介质 | |
CN109633746A (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 |