CN104122587A - 一种基于炮检距域的异常初至识别方法以及系统 - Google Patents
一种基于炮检距域的异常初至识别方法以及系统 Download PDFInfo
- Publication number
- CN104122587A CN104122587A CN201410305273.7A CN201410305273A CN104122587A CN 104122587 A CN104122587 A CN 104122587A CN 201410305273 A CN201410305273 A CN 201410305273A CN 104122587 A CN104122587 A CN 104122587A
- Authority
- CN
- China
- Prior art keywords
- arrival
- big gun
- time
- examine
- big
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于炮检距域的异常初至识别方法及系统,所述方法包括:采集地震数据;从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点;对抽取出的多个炮检对的初至应用基础静校正量,得到所述炮检对的第一时间;获取预先设定的阈值;根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。本发明利用在一个中心点平面中,应用基础静校正量的初至时间在共炮检距道集中是一个光滑平面的原理,避免了在共激发点道集中对初至进行多项式拟合的不足,适应二维勘探和三维勘探。
Description
技术领域
本发明关于地球物理勘探技术领域,特别是关于地震勘探中的静校正处理技术,具体的讲是一种基于炮检距域的异常初至识别方法及系统。
背景技术
地震勘探中的静校正处理技术主要是用来解决地震波在近地表介质旅行时间变化影响反射波叠加成像的问题。地震发生后,地震观测点最先接收到的波称初至波。由于各种地震波的传播速度不同,传播到观测点的时间也就有先后。利用初至波的折射静校正方法或层析静校正方法可以反演出相对精确的近地表模型,解决静校正问题。而初至的拾取精度,直接影响折射或层析反演的结果。地震波波前到达某个观测点,在观测点上,检波器检测到质点振动的时刻称为波的初至时间,简称初至。
为了提高初至拾取的精度,首要的任务是剔除异常的初至。现有技术在剔除异常初至的方法中,大都采取了在共激发点道集中对初至进行曲线拟合的方法,如果拟合值与初至相差较大,则认为是异常初至。
上述方法在拟合曲线的过程中,由于异常初至的存在,导致拟合曲线形态不符合初至波传播规律,以至于删除了正常的初至,有时还保留了异常初至。为了避免相同炮检距的初至来自不同的层或折射界面倾斜的情况,在计算过程中要求对初至进行区分左右支,对每支的初至进行多项式拟合,而多项式的幂由于异常初至的存在,在实际操作中较难选择。
因此,如何提出一种全新的初至拾取方案,其能够剔除异常的初至,提高初至拾取的精度,进而利用初至波的折射静校正方法或层析静校正方法反演出相对精确的近地表模型,解决静校正问题是本领域亟待解决的技术难题。
发明内容
为了克服现有技术中的初至拾取方案在剔除异常初至过程中存在的缺陷,本发明提供了一种基于炮检距域的异常初至识别方法及系统,是在基础静校正量的基础上删除异常初至时间的方案,利用在一个中心点平面中,应用基础静校正量的初至时间在共炮检距道集中是一个光滑平面的原理,解决了现有技术中的初至拾取方案在共激发点道集中对初至进行多项式拟合的不足,适应二维勘探和三维勘探。
本发明的目的之一是,提供一种基于炮检距域的异常初至识别方法,包括:采集地震数据;从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点;对抽取出的多个炮检对的初至应用基础静校正量,得到所述炮检对的第一时间;获取预先设定的阈值;根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。
本发明的目的之一是,提供了一种基于炮检距域的异常初至识别系统,包括:地震数据采集装置,用于采集地震数据;炮检对抽取装置,用于从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点;第一时间确定装置,用于对抽取出的多个炮检对的初至应用基础静校正量,得到所述的炮检对的第一时间;阈值获取装置,用于获取预先设定的阈值;异常初至识别装置,用于根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。
本发明的有益效果在于,提供了一种基于炮检距域的异常初至识别方法及系统,是地球物理勘探中在炮检距道集中删除异常初至的方案,利用相同炮检距的初至、在中心点分布的平面中,初至应该是一个比较光滑平面的原理,解决了现有技术中的初至拾取方案中在共激发点道集中对初至进行多项式拟合的不足,适应二维勘探和三维勘探,其能够剔除异常的初至,提高初至拾取的精度,进而实现了利用初至波的折射静校正方法或层析静校正方法反演出相对精确的近地表模型,解决静校正问题。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于炮检距域的异常初至识别方法的实施方式一的流程图;
图2为图1中的步骤S102的具体流程图;
图3为图1中的步骤S103的具体流程图;
图4为图1中的步骤S105的具体流程图;
图5为图4中的步骤S405的具体流程图;
图6为图4中的步骤S406的具体流程图;
图7为本发明实施例提供的一种基于炮检距域的异常初至识别方法的实施方式二的流程图;
图8为本发明实施例提供的一种基于炮检距域的异常初至识别系统的实施方式一的结构框图;
图9为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的炮检对抽取装置200的具体结构框图;
图10为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的第一时间确定装置300的具体结构框图;
图11为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的异常初至识别装置500的具体结构框图;
图12为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的第二时间确定模块505的具体结构框图;
图13为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的异常初至识别模块506的具体结构框图;
图14为本发明实施例提供的一种基于炮检距域的异常初至识别系统的实施方式二的结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明是地球物理勘探中利用统计的方法在炮检距道集中删除异常初至的方法。利用在一个中心点平面中,应用基础静校正量的初至时间在共炮检距道集中应该是一个光滑平面的原理。本发明避免了在共激发点道集中对初至进行多项式拟合的不足,适应二维勘探和三维勘探。
图1为本发明提出的一种基于炮检距域的异常初至识别方法的实施方式一的具体流程图,由图1可知,在实施方式一中,所述的方法包括:
S101:采集地震数据。
在具体的实施例中,可利用现有仪器采集得到地震数据。
S102:从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点,图2为步骤S102的具体流程图。
S103:对抽取出的多个炮检对的初至应用基础静校正量,得到所述炮检对的第一时间。图3为步骤S103的具体流程图。
S104:获取预先设定的阈值,阈值可根据不同的使用场景预先设定。
S105:根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。图4为步骤S105的具体流程图。
图2为步骤S102的具体流程图,由图2可知,该步骤具体包括:
S201:获取预先设定的炮检距,具体的实施例中,诸如在预先设定炮检距D=200m。
S202:获取预先设定的浮动值,具体的实施例中,诸如在预先设定浮动值d=10m。
S203:根据所述的炮检距以及浮动值确定炮检距范围,具体的实施例中,预先设定浮动值d=10,炮检距D=200,则炮检距范围为[190m,210m]。
S204:确定所述地震数据中的所有炮检对的炮检距,所述的炮检对包括初至、中点坐标、激发点以及接收点,所述的炮检距为激发点与接收点之间的距离;
S205:从所述的地震数据中抽取出炮检距在所述炮检距范围内的炮检对;
S206:抽取出的炮检对的中点坐标以及初至形成共炮检距道集。具体的实施例中,炮检对的中点坐标诸如为(x,y),形成的共炮检距道集为G。
图3为步骤S103的具体流程图,由图3可知,该步骤具体包括:
S301:获取抽取出的炮检对的激发点的基础静校正量;
S302:获取抽取出的炮检对的接收点的基础静校正量;
S303:确定所述激发点的基础静校正量与所述接收点的基础静校正量的和;
S304:将所述的和应用到所述的炮检对的初至上,得到所述的炮检对的第一时间。在具体的实施方式中,在共炮检距道集G中,将激发点的基础静校正量与所述接收点的基础静校正量的和加上所述初至即形成第一时间t1。
图4为步骤S105的具体流程图,由图4可知,步骤S105具体包括:
S401:从抽取出的多个炮检对中选出一个炮检对,作为选定炮检对;
S402:获取预先设定的半径,在具体的实施方式中,半径可以预先设定,诸如半径R可取炮检距道集G的最大炮检距的一半。
S403:将所述的选定炮检对的中点坐标作为圆心;
S404:根据所述的半径、所述的圆心做圆;
S405:分别确定抽取出的多个炮检对中除所述选定炮检对之外的其他炮检对的中点坐标在所述圆中的第二时间,图5为步骤S405的具体流程图。
S406:根据所述的第二时间以及所述的阈值从所述的多个初至中识别出异常初至。图6为步骤S406的具体流程图。
图5为图4中的步骤S405的具体流程图,由图5可知,该步骤具体包括:
S501:确定抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离;
S502:判断所述的距离是否小于所述圆的半径;
S503:当判断为是时,当前炮检对的第一时间即为第二时间;
S504:否则,当前炮检对的第二时间为零。
即,在炮检距道集G中,当抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离小于半径R时,炮检对的第二时间t2即为第一时间t1;当抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离大于等于半径R时,炮检对的第二时间即为0。
图6为图4中的步骤S406的具体流程图,由图6可知,该步骤具体包括:
S601:对抽取出的多个炮检对中除第二时间为零的炮检对之外的炮检对的第二时间求平均值,得到平均时间,所述的平均时间即为选定炮检对的平均时间,平均时间可通过下述公式求取:
其中,为选定炮检对的平均时间,t2是指在炮检距道集G中,抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离小于半径R时,其他炮检对的第二时间,也即炮检对中点坐标到当前炮检对中点坐标的距离小于半径R的时间t1,N为抽取出的多个炮检对中除第二时间为零的炮检对之外的炮检对的总个数,i为序号。
S602:获取所述选定炮检对的初至;
S603:确定所述初至与所述平均时间的差值的绝对值;
S604:当所述差值的绝对值超出所述的阈值时,所述选定炮检对的初至即为异常初至。
图7为本发明实施例提供的一种基于炮检距域的异常初至识别方法的实施方式二的流程图,由图7可知,在实施方式二中,该方法具体包括:
S701:采集地震数据。
在具体的实施例中,可利用现有仪器采集得到地震数据。
S702:从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点。
S703:对抽取出的多个炮检对的初至应用基础静校正量,得到所述炮检对的第一时间。
S704:获取预先设定的阈值,阈值可根据不同的使用场景预先设定。
S705:根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。
S706:将所述的异常初至从所述的初至中剔除。
也即,如果选定炮检对的初至与的差值大于给定的门槛值,则认为是异常初至,可以剔除。如此,从所述的地震数据中抽取出多个炮检对对应的多个初至剔除异常初至之后,即可得到正常初至。
如上所述,即为本发明提供的一种基于炮检距域的异常初至识别方法,利用的相同炮检距的初至、在中心点分布的平面中,初至应该是一个比较光滑平面的原理,实现了在共炮检距域中剔除异常初至,避免了在共激发点道集中对初至进行多项式拟合的不足。该方法可以同时适应二维勘探和三维勘探。
图8为本发明实施例提供的一种基于炮检距域的异常初至识别系统的实施方式一的结构框图,由图8可知,在实施方式一中,所述的系统包括:
地震数据采集装置100,用于采集地震数据。
在具体的实施例中,可利用现有仪器采集得到地震数据。
炮检对抽取装置200,用于从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点,图9为炮检对抽取装置200的具体结构框图。
第一时间确定装置300,用于对抽取出的多个炮检对的初至应用基础静校正量,得到所述炮检对的第一时间。图10为第一时间确定装置300的具体结构框图。
阈值获取装置400,用于获取预先设定的阈值,阈值可根据不同的使用场景预先设定。
异常初至识别装置500,用于根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。图11为异常初至识别装置500具体结构框图。
图9为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的炮检对抽取装置200的具体结构框图,由图9可知,该炮检对抽取装置200具体包括:
炮检距获取模块201,用于获取预先设定的炮检距,具体的实施例中,诸如在预先设定炮检距D=200m。
浮动值获取模块202,用于获取预先设定的浮动值,具体的实施例中,诸如在预先设定浮动值d=10m。
炮检距范围确定模块203,用于根据所述的炮检距以及浮动值确定炮检距范围,具体的实施例中,预先设定浮动值d=10,炮检距D=200,则炮检距范围为[190m,210m]。
炮检距确定模块204,用于确定所述地震数据中的所有炮检对的炮检距,所述的炮检对包括初至、中点坐标、激发点以及接收点,所述的炮检距为激发点与接收点之间的距离;
炮检对抽取模块205,用于从所述的地震数据中抽取出炮检距在所述炮检距范围内的炮检对;
共炮检距道集形成模块206,用于抽取出的炮检对的中点坐标以及初至形成共炮检距道集。具体的实施例中,炮检对的中点坐标诸如为(x,y),形成的共炮检距道集为G。
图10为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的第一时间确定装置300的具体结构框图,由图10可知,该第一时间确定装置300具体包括:
第一静校正量获取模块301,用于获取抽取出的炮检对的激发点的基础静校正量;
第二静校正量获取模块302,用于获取抽取出的炮检对的接收点的基础静校正量;
静校正量和确定模块303,用于确定所述激发点的基础静校正量与所述接收点的基础静校正量的和;
第一时间确定模块304,用于将所述的和应用到所述的炮检对的初至上,得到所述的炮检对的第一时间。在具体的实施方式中,在共炮检距道集G中,将激发点的基础静校正量与所述接收点的基础静校正量的和加上所述初至即形成第一时间t1。
图11为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的异常初至识别装置500的具体结构框图,由图11可知,异常初至识别装置500具体包括:
炮检对选定模块501,用于从抽取出的多个炮检对中选出一个炮检对,作为选定炮检对;
半径获取模块502,用于获取预先设定的半径,在具体的实施方式中,半径可以预先设定,诸如半径R可取炮检距道集G的最大炮检距的一半。
圆心选定模块503,用于将所述的选定炮检对的中点坐标作为圆心;
圆确定模块504,用于根据所述的半径、所述的圆心做圆;
第二时间确定模块505,用于分别确定抽取出的多个炮检对中除所述选定炮检对之外的其他炮检对的中点坐标在所述圆中的第二时间,图12为第二时间确定模块505的具体结构框图。
第二时间确定模块506,用于根据所述的第二时间以及所述的阈值从所述的多个初至中识别出异常初至。图13为异常初至识别模块506的具体结构框图。
图12为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的第二时间确定模块505的具体结构框图,由图12可知,该第二时间确定模块具体包括:
距离确定单元5051,用于确定抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离;
判断单元5052,用于判断所述的距离是否小于所述圆的半径;
第一确定单元5053,用于当所述的判断单元判断为是时,当前炮检对的第一时间即为第二时间;
第二确定单元5054,用于当所述的判断单元判断为否时,当前炮检对的第二时间为零。
即,在炮检距道集G中,当抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离小于半径R时,炮检对的第二时间t2即为第一时间t1;当抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离大于等于半径R时,炮检对的第二时间即为0。
图13为本发明实施例提供的一种基于炮检距域的异常初至识别系统中的异常初至识别模块506的具体结构框图,由图13可知,该异常初至识别模块506具体包括:
平均时间确定单元5061,用于对抽取出的多个炮检对中除第二时间为零的炮检对之外的炮检对的第二时间求平均值,得到平均时间,所述的平均时间即为选定炮检对的平均时间,平均时间可通过下述公式求取:
其中,为选定炮检对的平均时间,t2是指在炮检距道集G中,抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离小于半径R时,其他炮检对的第二时间,也即炮检对中点坐标到当前炮检对中点坐标的距离小于半径R的时间t1,N为抽取出的多个炮检对中除第二时间为零的炮检对之外的炮检对的总个数,i为序号。
第一时间获取单元5062,用于获取所述选定炮检对的初至;
绝对值确定单元5063,用于确定所述初至与所述平均时间的差值的绝对值;
异常初至判定单元5064,用于当所述差值的绝对值超出所述的阈值时,所述选定炮检对的初至即为异常初至。
图14为本发明实施例提供的一种基于炮检距域的异常初至识别系统的实施方式二的结构框图,由图14可知,在实施方式二中,该系统还包括:
初至剔除装置600,用于将所述的异常初至从所述的初至中剔除。
也即,如果选定炮检对的初至与的差值大于给定的门槛值,则认为是异常初至,可以剔除。如此,从所述的地震数据中抽取出多个炮检对对应的多个初至剔除异常初至之后,即可得到正常初至。
如上所述,即为本发明提供的一种基于炮检距域的异常初至识别系统,利用的相同炮检距的初至、在中心点分布的平面中,初至应该是一个比较光滑平面的原理,实现了在共炮检距域中剔除异常初至,避免了在共激发点道集中对初至进行多项式拟合的不足。该方案可以同时适应二维勘探和三维勘探。
下面结合具体的实施例,详细介绍本发明的技术方案。
将在共炮检距域中剔除异常初至的方法应用于某地震勘探实际应用中,道距为10米、炮检距200m到2000m,则:
1)采集地震数据,抽取共炮检距道集:给出一个炮检距(激发点到接收点之间的距离)D=200和浮动范围d=10,抽取所有炮检距在范围[190,210]中的初至和相应炮检对(一个激发点和一个接收点)的中点坐标(x,y),形成一个共炮检距道集G。
2)应用基础静校正量:可以是高程静校正量,在道集G中,把已知激发点和接收点基础静校正量的和,应用到初至上,形成时间t1。
3)删除异常初至:在道集G中,以一个炮检对的中点坐标(x,y)为圆心,以给出的平滑半径R(如1000m)为半径形成一个圆,计算平均时间:
其中,t1i是指在道集G中,其它炮检对中点坐标到当前炮检对中点坐标的距离小于给定半径R的t1,N是相应的炮检对的总个数。
如果给定炮检对的初至与t2的差值的绝对值大于用户给定的门槛值,则认为是异常初至,可以剔除。
综上所述,本发明提出的一种基于炮检距域的异常初至识别方法及系统,是地球物理勘探中在炮检距道集中删除异常初至的方案,利用相同炮检距的初至、在中心点分布的平面中,初至应该是一个比较光滑平面的原理,解决了现有技术中的初至拾取方案中在共激发点道集中对初至进行多项式拟合的不足,适应二维勘探和三维勘探,其能够剔除异常的初至,提高初至拾取的精度,进而实现了利用初至波的折射静校正方法或层析静校正方法反演出相对精确的近地表模型,解决静校正问题。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一般计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
本领域技术人员还可以了解到本发明实施例列出的各种功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (14)
1.一种基于炮检距域的异常初至识别方法,其特征是,所述的方法具体包括:
采集地震数据;
从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点;
对抽取出的多个炮检对的初至应用基础静校正量,得到所述炮检对的第一时间;
获取预先设定的阈值;
根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。
2.根据权利要求1所述的方法,其特征是,从所述的地震数据中抽取多个炮检对具体包括:
获取预先设定的炮检距;
获取预先设定的浮动值;
根据所述的炮检距以及浮动值确定炮检距范围;
确定所述地震数据中的所有炮检对的炮检距,所述的炮检对包括初至、中点坐标、激发点以及接收点,所述的炮检距为激发点与接收点之间的距离;
从所述的地震数据中抽取出炮检距在所述炮检距范围内的炮检对;
抽取出的炮检对的中点坐标以及初至形成共炮检距道集。
3.根据权利要求2所述的方法,其特征是,对抽取出的多个炮检对的初至应用基础静校正量,得到所述的炮检对第一时间具体包括:
获取抽取出的炮检对的激发点的基础静校正量;
获取抽取出的炮检对的接收点的基础静校正量;
确定所述激发点的基础静校正量与所述接收点的基础静校正量的和;
将所述的和应用到所述的炮检对的初至上,得到所述的炮检对的第一时间。
4.根据权利要求3所述的方法,其特征是,根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的多个初至中识别出异常初至具体包括:
从抽取出的多个炮检对中选出一个炮检对,作为选定炮检对;
获取预先设定的半径;
将所述的选定炮检对的中点坐标作为圆心;
根据所述的半径、所述的圆心做圆;
分别确定抽取出的多个炮检对中除所述选定炮检对之外的其他炮检对的中点坐标在所述圆中的第二时间;
根据所述的第二时间以及所述的阈值从所述的多个初至中识别出异常初至。
5.根据权利要求4所述的方法,其特征是,分别确定抽取出的多个炮检对中除所述选定炮检对之外的其他炮检对的中点坐标在所述圆中的第二时间具体包括:
确定抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离;
判断所述的距离是否小于所述圆的半径;
当判断为是时,当前炮检对的第一时间即为第二时间;
否则,当前炮检对的第二时间为零。
6.根据权利要求5所述的方法,其特征是,根据所述的第二时间以及所述的阈值从所述的初至中识别出异常初至具体包括:
对抽取出的多个炮检对中除第二时间为零的炮检对之外的炮检对的第二时间求平均值,得到平均时间,所述的平均时间即为选定炮检对的平均时间;
获取所述选定炮检对的初至;
确定所述初至与所述平均时间的差值的绝对值;
当所述差值的绝对值超出所述的阈值时,所述选定炮检对的初至即为异常初至。
7.根据权利要求1至6任意一项所述的方法,其特征是,在识别出异常初至之后,所述的方法还包括:
将所述的异常初至从所述的初至中剔除。
8.一种基于炮检距域的异常初至识别系统,其特征是,所述的系统具体包括:
地震数据采集装置,用于采集地震数据;
炮检对抽取装置,用于从所述的地震数据中抽取出多个炮检对,所述的炮检对包括初至、中点坐标、激发点以及接收点;
第一时间确定装置,用于对抽取出的多个炮检对的初至应用基础静校正量,得到所述的炮检对的第一时间;
阈值获取装置,用于获取预先设定的阈值;
异常初至识别装置,用于根据所述的第一时间、抽取出的多个炮检对以及所述的阈值从抽取出的多个炮检对的初至中识别出异常初至。
9.根据权利要求8所述的系统,其特征是,所述的炮检对抽取装置具体包括:
炮检距获取模块,用于获取预先设定的炮检距;
浮动值获取模块,用于获取预先设定的浮动值;
炮检距范围确定模块,用于根据所述的炮检距以及浮动值确定炮检距范围;
炮检距确定模块,用于确定所述地震数据中的所有炮检对的炮检距,所述的炮检对包括初至、中点坐标、激发点以及接收点,所述的炮检距为激发点与接收点之间的距离;
炮检对抽取模块,用于从所述的地震数据中抽取出炮检距在所述炮检距范围内的炮检对;
共炮检距道集形成模块,用于抽取出的炮检对的中点坐标以及初至形成共炮检距道集。
10.根据权利要求9所述的系统,其特征是,所述的第一时间确定装置具体包括:
第一静校正量获取模块,用于获取抽取出的炮检对的激发点的基础静校正量;
第二静校正量获取模块,用于获取抽取出的炮检对的接收点的基础静校正量;
静校正量和确定模块,用于确定所述激发点的基础静校正量与所述接收点的基础静校正量的和;
第一时间确定模块,用于将所述的和应用到所述的炮检对的初至上,得到所述的炮检对的第一时间。
11.根据权利要求10所述的系统,其特征是,所述的异常初至识别装置具体包括:
炮检对选定模块,用于从抽取出的多个炮检对中选出一个炮检对,作为选定炮检对;
半径获取模块,用于获取预先设定的半径;
圆心选定模块,用于将所述的选定炮检对的中点坐标作为圆心;
圆确定模块,用于根据所述的半径、所述的圆心做圆;
第二时间确定模块,用于分别确定抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标在所述圆中的第二时间;
异常初至识别模块,用于根据所述的第二时间以及所述的阈值从所述的多个初至中识别出异常初至。
12.根据权利要求11所述的系统,其特征是,所述的第二时间确定模块具体包括:
距离确定单元,用于确定抽取出的多个炮检对中除所述选定炮检对之外的炮检对的中点坐标到所述选定炮检对的中点坐标的距离;
判断单元,用于判断所述的距离是否小于所述圆的半径;
第一确定单元,用于当所述的判断单元判断为是时,当前炮检对的第一时间即为第二时间;
第二确定单元,用于当所述的判断单元判断为否时,当前炮检对的第二时间为零。
13.根据权利要求12所述的系统,其特征是,所述的异常初至识别模块具体包括:
平均时间确定单元,用于对抽取出的多个炮检对中除第二时间为零的炮检对之外的炮检对的第二时间求平均值,得到平均时间,所述的平均时间即为选定炮检对的平均时间;
第一时间获取单元,用于获取所述选定炮检对的第一时间;
绝对值确定单元,用于确定所述第一时间与所述平均时间的差值的绝对值;
异常初至判定单元,用于当所述差值的绝对值超出所述的阈值时,所述选定炮检对的初至即为异常初至。
14.根据权利要求8至13任意一项所述的系统,其特征是,所述的系统还包括:
初至剔除装置,用于将所述的异常初至从所述的初至中剔除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410305273.7A CN104122587B (zh) | 2014-06-30 | 2014-06-30 | 一种基于炮检距域的异常初至识别方法以及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410305273.7A CN104122587B (zh) | 2014-06-30 | 2014-06-30 | 一种基于炮检距域的异常初至识别方法以及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104122587A true CN104122587A (zh) | 2014-10-29 |
CN104122587B CN104122587B (zh) | 2016-10-19 |
Family
ID=51768074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410305273.7A Active CN104122587B (zh) | 2014-06-30 | 2014-06-30 | 一种基于炮检距域的异常初至识别方法以及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104122587B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375186A (zh) * | 2014-12-08 | 2015-02-25 | 中国石油天然气集团公司 | 一种去除异常初至的方法及系统 |
CN106842302A (zh) * | 2016-12-26 | 2017-06-13 | 中国石油天然气股份有限公司 | 一种批量编辑初至的方法及装置 |
CN104459779B (zh) * | 2014-12-09 | 2017-07-28 | 中国石油天然气集团公司 | 一种异常地震道的自动识别的方法和装置 |
CN107797147A (zh) * | 2017-10-25 | 2018-03-13 | 中国石油化工股份有限公司 | 一种地震初至异常值快速剔除方法 |
CN109143329A (zh) * | 2017-06-27 | 2019-01-04 | 中国石油化工股份有限公司 | 初至拾取处理方法 |
CN109738946A (zh) * | 2018-11-21 | 2019-05-10 | 中国石油天然气集团有限公司 | 两步法剔除异常地震初至的方法及系统 |
CN112666606A (zh) * | 2019-10-15 | 2021-04-16 | 中国石油天然气集团有限公司 | 确定近地表方位各向异性速度的方法及装置 |
CN114114391A (zh) * | 2020-09-01 | 2022-03-01 | 中国石油天然气股份有限公司 | 地震数据负面样本识别方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020075759A1 (en) * | 2001-01-21 | 2002-06-20 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN102478666A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种提高可控震源地震记录初至拾取精度的方法 |
CN103823241A (zh) * | 2012-11-16 | 2014-05-28 | 中国石油天然气集团公司 | 一种在炮检距域中求取折射剩余静校正量的方法 |
CN103837895A (zh) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 拟合初至波获取短波长静校正量的方法 |
-
2014
- 2014-06-30 CN CN201410305273.7A patent/CN104122587B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020075759A1 (en) * | 2001-01-21 | 2002-06-20 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN102478666A (zh) * | 2010-11-23 | 2012-05-30 | 中国石油天然气集团公司 | 一种提高可控震源地震记录初至拾取精度的方法 |
CN103823241A (zh) * | 2012-11-16 | 2014-05-28 | 中国石油天然气集团公司 | 一种在炮检距域中求取折射剩余静校正量的方法 |
CN103837895A (zh) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 拟合初至波获取短波长静校正量的方法 |
Non-Patent Citations (1)
Title |
---|
徐钰 等: ""浅层折射波勘探中初至自动拾取新算法"", 《石油地球物理勘探》, vol. 47, no. 2, 30 April 2012 (2012-04-30) * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375186A (zh) * | 2014-12-08 | 2015-02-25 | 中国石油天然气集团公司 | 一种去除异常初至的方法及系统 |
CN104459779B (zh) * | 2014-12-09 | 2017-07-28 | 中国石油天然气集团公司 | 一种异常地震道的自动识别的方法和装置 |
CN106842302A (zh) * | 2016-12-26 | 2017-06-13 | 中国石油天然气股份有限公司 | 一种批量编辑初至的方法及装置 |
CN106842302B (zh) * | 2016-12-26 | 2019-06-11 | 中国石油天然气股份有限公司 | 一种批量编辑初至的方法及装置 |
CN109143329A (zh) * | 2017-06-27 | 2019-01-04 | 中国石油化工股份有限公司 | 初至拾取处理方法 |
CN107797147A (zh) * | 2017-10-25 | 2018-03-13 | 中国石油化工股份有限公司 | 一种地震初至异常值快速剔除方法 |
CN109738946A (zh) * | 2018-11-21 | 2019-05-10 | 中国石油天然气集团有限公司 | 两步法剔除异常地震初至的方法及系统 |
CN112666606A (zh) * | 2019-10-15 | 2021-04-16 | 中国石油天然气集团有限公司 | 确定近地表方位各向异性速度的方法及装置 |
CN112666606B (zh) * | 2019-10-15 | 2024-06-25 | 中国石油天然气集团有限公司 | 确定近地表方位各向异性速度的方法及装置 |
CN114114391A (zh) * | 2020-09-01 | 2022-03-01 | 中国石油天然气股份有限公司 | 地震数据负面样本识别方法及装置 |
CN114114391B (zh) * | 2020-09-01 | 2023-08-22 | 中国石油天然气股份有限公司 | 地震数据负面样本识别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104122587B (zh) | 2016-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104122587A (zh) | 一种基于炮检距域的异常初至识别方法以及系统 | |
CN105607121B (zh) | 一种煤炭陷落柱识别方法和装置 | |
CN104656141B (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN101980054B (zh) | 一种在高密度地震静校正处理中建立近地表速度模型的方法 | |
CN105022031B (zh) | 一种区域岩体微震震源的分层速度定位方法 | |
CN105785439B (zh) | 小尺度非均匀地质体空间分布位置的预测方法和装置 | |
CN104181599A (zh) | 一种基于近地表层的折射波静校正处理方法以及系统 | |
US20090067287A1 (en) | Complex analysis of kinematics for non-hyperbolic moveout corrections | |
CN105093299A (zh) | 一种基于炮检距向量片技术优化观测系统的方法及装置 | |
CN105954797B (zh) | 地震数据的断裂识别方法和装置 | |
CN104483705A (zh) | 一种三维剩余静校正方法 | |
CN104199103B (zh) | 转换波静校正方法及装置 | |
CN104123786A (zh) | 纸币卡钞判断系统及方法 | |
CN106772593A (zh) | 绕射波的成像方法及装置 | |
CN105427451A (zh) | 一种检测纸币分钞过程中倾斜的方法及系统 | |
CN105866839B (zh) | 基于共姿态道集的静校正方法及装置 | |
CN104977615A (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN106199705B (zh) | 地震勘探采集数据的检测方法 | |
CN102053275B (zh) | 一种用于单点地震室内组合的相对静校正量计算方法 | |
CN104216016A (zh) | 一种纵波约束扫描的转换波动校正方法以及系统 | |
CN104155703A (zh) | 评价三维观测系统的方法和设备 | |
CN103217708A (zh) | 获得海上多缆拖缆地震采集全单缆单源地震剖面的方法 | |
Wang et al. | BSPASS: A Beam Search‐Based Phase Association and Source Scanning Earthquake Location Method | |
CN104865601B (zh) | 确定地震勘探中采集的地震记录的类型的方法 | |
CN103837890B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |