CN112083480B - 同步激发采集地震数据的实时监控方法及装置 - Google Patents
同步激发采集地震数据的实时监控方法及装置 Download PDFInfo
- Publication number
- CN112083480B CN112083480B CN202010855259.XA CN202010855259A CN112083480B CN 112083480 B CN112083480 B CN 112083480B CN 202010855259 A CN202010855259 A CN 202010855259A CN 112083480 B CN112083480 B CN 112083480B
- Authority
- CN
- China
- Prior art keywords
- seismic
- gun
- shot
- synchronous
- excitation
- 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
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/24—Recording seismic data
- G01V1/247—Digital recording of seismic data, e.g. in acquisition units or nodes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/22—Transmitting seismic signals to recording or processing apparatus
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
- G01V1/26—Reference-signal-transmitting devices, e.g. indicating moment of firing of shot
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种同步激发采集地震数据的实时监控方法及装置,该方法包括:根据地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻,确定各个地震道的环境噪音,监控地震数据的采集质量。提高了初至时刻拾取的准确性,真实反映环境噪音水平。
Description
技术领域
本发明涉及地球物理勘探技术领域,特别涉及一种同步激发采集地震数据的实时监控方法及装置。
背景技术
宽方位和高密度地震数据采集是陆上地震勘探发展的方向和未来,同步激发地震数据采集过程中,多个炮点同步独立激发,该方法极大的压缩了地震数据、有效提高了生产效率,是实现陆上低成本宽方位和高密度地震数据采集的有效途径。但是,同步激发地震数据采集中不同炮点同步激发产生的地震波场相互干扰、混叠,现有的地震数据采集实时质量控制方法容易将同步激发产生的地震波加入到环境噪音的计算中,使得环境噪音的计算精度较低,不能真实反映地震数据采集过程中的环境噪音,难以满足同步激发地震数据采集实时质量监控的需求。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供一种同步激发采集地震数据的实时监控方法,用于提高环境噪音的计算精度,真实反映地震数据采集过程中的环境噪音,满足同步激发采集地震数据质量监控的需求,该方法包括:
根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;
获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;
根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;其中,初至时刻与炮检距之间的关联关系根据样本炮的初至时刻和炮检距关系确定,所述样本炮用于确定初至时刻与炮检距之间的关联关系;
根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻;
根据各个地震道的最小初至时刻,确定各个地震道的环境噪音;
根据各个地震道的环境噪音,监控地震数据的采集质量。
本发明实施例提供一种同步激发采集地震数据的实时监控装置,用于提高环境噪音的计算精度,真实反映地震数据采集过程中的环境噪音,满足同步激发采集地震数据质量监控的需求,该装置包括:
同步激发时刻范围确定模块,用于根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;
同步激发炮筛选模块,用于获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;
炮检距确定模块,用于根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
初至时刻确定模块,用于根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;其中,初至时刻与炮检距之间的关联关系根据样本炮的初至时刻和炮检距关系确定,所述样本炮用于确定初至时刻与炮检距之间的关联关系;
最小初至时刻确定模块,用于根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻;
环境噪音确定模块,用于根据各个地震道的最小初至时刻,确定各个地震道的环境噪音;
监控模块,用于根据各个地震道的环境噪音,监控地震数据的采集质量。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述同步激发采集地震数据的实时监控方法。
本发明实施例还提供了一种计算机可读存储介质,计算机可读存储介质存储有执行上述同步激发采集地震数据的实时监控方法的计算机程序。
本发明实施例通过:根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻,能够得到准确的最小初至时刻;根据各个地震道的最小初至时刻,确定各个地震道的环境噪音,可以提高环境噪音的计算精度,真实反映地震数据采集过程中的环境噪音;根据各个地震道的环境噪音,监控地震数据的采集质量,实现了地震数据采集过程中基于环境噪音的数据采集质量实时监控,提高了数据采集质量监控效率和精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中同步激发采集地震数据的实时监控方法流程的示意图;
图2为本发明实施例中同步激发采集地震数据的实时监控装置结构的示意图;
图3为本发明实施例中具体实施例的流程图;
图4为某地区同步激发采集得到的地震数据和采用现有技术进行地震数据质量监控得到的初至时刻;
图5为图4的局部放大图;
图6为图4对应的地震数据和采用本发明实施例进行地震数据质量监控得到的初至时刻;
图7为图6的局部放大图;
图8为采用现有技术得到的环境噪音水平曲线;
图9为采用本发明实施例得到的环境噪音水平曲线;
图10为采用现有技术得到的环境噪音水平统计柱状图;
图11为采用本发明实施例得到的环境噪音水平统计柱状图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本领域技术技术人员知道,本发明的实施方式可以实现为一种系统、装置、方法或计算机程序产品。因此,本发明公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
为了解决现有的地震数据采集的实时质量控制方法容易将同步激发产生的地震波加入到环境噪音的计算中,导致环境噪音的计算精度较低的技术问题,本发明实施例提供一种同步激发采集地震数据的实时监控方法,用于提高环境噪音的计算精度,真实反映地震数据采集过程中的环境噪音,满足同步激发采集地震数据质量监控的需求,图1为本发明实施例中同步激发采集地震数据的实时监控方法流程的示意图,如图1所示,该方法包括:
步骤101:根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;
步骤102:获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;
步骤103:根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
步骤104:根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;其中,初至时刻与炮检距之间的关联关系根据样本炮的初至时刻和炮检距关系确定;
步骤105:根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻;
步骤106:根据各个地震道的最小初至时刻,确定各个地震道的环境噪音;
步骤107:根据各个地震道的环境噪音,监控地震数据的采集质量。
如图1所示,本发明实施例通过:根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻,能够得到准确的最小初至时刻;根据各个地震道的最小初至时刻,确定各个地震道的环境噪音,可以提高环境噪音的计算精度,真实反映地震数据采集过程中的环境噪音;根据各个地震道的环境噪音,监控地震数据的采集质量,实现了地震数据采集过程中基于环境噪音的数据采集质量实时监控,提高了数据采集质量监控效率和精度。
具体实施时,在执行步骤101之前,首先需要根据样本炮的初至时刻和炮检距关系确定初至时刻与炮检距之间的关联关系,可以在地震采集作业工区内应用常规地震采集方法得到一个或多个具有代表性的地震记录并拾取初至时刻,得到一组或多组具有代表性的初至时刻和炮检距对应关系函数作为整个工区或工区内特定区域的标准初至时刻和炮检距对应关系函数,如公式(1)所示:
{Tj(δxj)|j=1,…,m} (1)
其中,j为地震道的序列号,m为地震道的数目,δxj为地震道j对应于样本炮点的炮检距,Tj为地震道j对应于样本炮点的初至时刻。
在一个实施例中,步骤101中可以按照如下公式(2)、公式(3)确定炮的同步激发时刻的范围:
Tmax=T0+T (2)
Tmin=T0-T (3)
其中,T0为本炮的激发时刻;T为地震数据的记录时长;Tmax为炮的同步激发时刻的最大值,Tmin为炮的同步激发时刻的最小值。
具体实施时,步骤101中,可以根据地震数据采集器输出的地震记录的本炮的激发时刻T0和地震数据的记录时长T应用公式(2)和公式(3)计算该地震记录中可能包含的同步激发炮的同步激发时刻的取值范围的最大值Tmax和最小值Tmin。
在一个实施例中,步骤102中,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮,包括:
将激发时刻大于或等于炮的同步激发时刻的最小值,且小于或等于炮的同步激发时刻的最大值的炮保留;
将激发时刻小于炮的同步激发时刻的最小值,和大于炮的同步激发时刻的最大值的炮删除。
具体实施时,步骤102中,可以从传送回来的多个炮的激发时刻中筛选出激发时刻落入上述炮的同步激发时刻的取值范围内的炮,得到筛选后的多个同步激发炮,并记录对应的炮的炮点位置和激发时刻,可以按照公式(4)进行筛选:
{(Xk,Tk)|k=1,…,l};Tmin≤Tk≤Tmax (4)
其中,Xk为炮的炮点位置,Tk为炮的激发时刻,k为符合炮的同步激发时刻的取值范围的炮点序列号,l为符合炮的同步激发时刻的取值范围的炮的炮点数目。
在一个实施例中,步骤103根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,包括:按照如下方式确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距:
其中,i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;j为地震道的顺序号;m为地震道的数目;xi为i炮点的东坐标;yi为i炮点的北坐标;xj为j地震道的东坐标,yj为j地震道的北坐标;δxi,j为j地震道与i炮点之间的炮检距。
具体实施时,步骤103中,预设的炮点坐标和地震道坐标存储在SPS文件中,可以将上述步骤102中得到的筛选后的多个同步激发炮的炮点位置与SPS文件进行匹配,得到各个同步激发炮的炮点坐标和各个地震道的坐标,根据各个同步激发炮的炮点坐标和各个地震道的坐标,分别计算各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距。
在一个实施例中,步骤104中,可以按照公式(6)确定各个地震道对应于本炮和各个同步激发炮的初至时刻:
Ti,j=Ti,j(δxi,j)+Ti-T0,i=1,...,n;j=1,...,m (6)
其中,T0为本炮的激发时刻;i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;i大于0时Ti为同步激发炮的激发时刻;j为地震道的顺序号;m为地震道的数目;δxi,j为j地震道与i炮点之间的炮检距;Ti,j为j地震道对应于i炮点所在炮的初至时刻。
具体实施时,步骤104中,可以根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及上述公式(1)中初至时刻与炮检距之间的关联关系,基于上述公式(6)计算各个地震道对应于本炮和各个同步激发炮的初至时刻。
具体实施时,步骤105中,对于每个地震道,可以将各个地震道对应于本炮和各个同步激发炮的初至时刻的最小值确定为该地震道的最终初至时刻,如公式(7)所示:
Tj =min (Ti,j) (7)
i=0,1,...,n;j=1,...,m
其中,i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;j为地震道的顺序号;m为地震道的数目;Ti,j为j地震道对应于i炮点所在炮的初至时刻;Tj为j地震道对应的最终初至时刻。
在一个实施例中,步骤106中,可以按照如下公式(8)确定地震道的环境噪音:
其中,j为地震道的顺序号;m为地震道的数目;Aj为j地震道的环境噪音;l为j地震道内样点的顺序号,Nj为j地震道最小初至时刻对应的样点的顺序号;Al,j为j地震道l样点的振幅值。
具体实施时,步骤106可以根据零时刻至最小初至时刻之间的样点对应的噪音振幅值,计算各个地震道的环境噪音,步骤107可以根据各个地震道的环境噪音,判断各个地震道的环境噪音是否符合预设阈值要求,从而实现对地震数据的采集质量的监控。
应当注意,尽管在附图中以特定顺序描述了本发明方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
基于同一发明构思,本发明实施例中还提供了一种同步激发采集地震数据的实时监控装置,如下面的实施例。由于同步激发采集地震数据的实时监控装置解决问题的原理与同步激发采集地震数据的实时监控方法相似,因此装置的实施可以参见方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图2为本发明实施例中同步激发采集地震数据的实时监控装置结构的示意图,如图2所示,该装置包括:
同步激发时刻范围确定模块01,用于根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;
同步激发炮筛选模块02,用于获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;
炮检距确定模块03,用于根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
初至时刻确定模块04,用于根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;其中,初至时刻与炮检距之间的关联关系根据样本炮的初至时刻和炮检距关系确定;
最小初至时刻确定模块05,用于根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻;
环境噪音确定模块06,用于根据各个地震道的最小初至时刻,确定各个地震道的环境噪音;
监控模块07,用于根据各个地震道的环境噪音,监控地震数据的采集质量。
在一个实施例中,还包括,关联关系建立模块08,用于建立初至时刻与炮检距之间的关联关系,初至时刻与炮检距之间的关联关系如下:
{Tj(δxj)|j=1,…,m};
其中,j为地震道的顺序号,δxj为地震道j对应于样本炮的炮点的炮检距,Tj为地震道j对应于样本炮的初至时刻,m为地震道的数目。
在一个实施例中,同步激发时刻范围确定模块01具体用于:
按照如下方式确定炮的同步激发时刻的范围:
Tmax=T0+T;
Tmin=T0-T;
其中,T0为本炮的激发时刻;T为地震数据的记录时长;Tmax为炮的同步激发时刻的最大值,Tmin为炮的同步激发时刻的最小值。
在一个实施例中,同步激发炮筛选模块02具体用于:
将激发时刻大于或等于炮的同步激发时刻的最小值,且小于或等于炮的同步激发时刻的最大值的炮保留;
将激发时刻小于炮的同步激发时刻的最小值,和大于炮的同步激发时刻的最大值的炮删除。
在一个实施例中,炮检距确定模块03具体用于:
按照如下方式确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距:
其中,i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;j为地震道的顺序号;m为地震道的数目;xi为i炮点的东坐标;yi为i炮点的北坐标;xj为j地震道的东坐标,yj为j地震道的北坐标;δxi,j为j地震道与i炮点之间的炮检距。
在一个实施例中,初至时刻确定模块04具体用于:
按照如下方式确定各个地震道对应于本炮和各个同步激发炮的初至时刻:
Ti,j=Ti,j(δxi,j)+Ti-T0,i=0,...,n;j=1,...,m;
其中,T0为本炮的激发时刻;i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;i大于0时Ti为同步激发炮的激发时刻;j为地震道的顺序号;m为地震道的数目;δxi,j为j地震道与i炮点之间的炮检距;Ti,j为j地震道对应于i炮点所在炮的初至时刻。
在一个实施例中,环境噪音确定模块06具体用于:
按照如下方式确定各个地震道的环境噪音:
其中,j为地震道的顺序号;m为地震道的数目;Aj为j地震道的环境噪音;l为j地震道内样点的顺序号,Nj为j地震道最小初至时刻对应的样点的顺序号;Al,j为j地震道l样点的振幅值。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述同步激发采集地震数据的实时监控方法。
本发明实施例还提供了一种计算机可读存储介质,计算机可读存储介质存储有执行上述同步激发采集地震数据的实时监控方法的计算机程序。
下面举一个具体的例子,以便于理解本发明如何实施。
图3为本发明实施例中具体实施例的流程图,如图3所示,包括:
第一步:在地震采集作业工区内应用常规地震采集方法得到一个或多个具有代表性的地震记录并拾取初至时刻,得到一组或多组具有代表性的初至时刻和炮检距对应关系函数作为整个工区或工区内特定区域的标准初至时刻和炮检距对应关系函数,如公式(1)所示;
第二步:根据地震数据采集器输出的地震记录的本炮的激发时刻T0和地震数据的记录时长T应用公式(2)和公式(3)计算该地震记录中可能包含的同步激发炮的同步激发时刻的取值范围的最大值Tmax和最小值Tmin;
第三步:按照公式(4),从传送回来的多个炮的激发时刻中筛选出激发时刻落入上述炮的同步激发时刻的取值范围内的炮,得到筛选后的多个同步激发炮,并记录对应的炮的炮点位置和激发时刻;
第四步:将筛选后的多个同步激发炮的炮点位置与SPS文件进行匹配,得到各个同步激发炮的炮点坐标和各个地震道的坐标,根据各个同步激发炮的炮点坐标和各个地震道的坐标,应用公式(5),分别计算各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
第五步:根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及上述公式(1)中初至时刻与炮检距之间的关联关系,基于上述公式(6)计算各个地震道对应于本炮和各个同步激发炮的初至时刻;
第六步:按照公式(7),将各个地震道对应于本炮和各个同步激发炮的初至时刻的最小值确定为该地震道的最终初至时刻;
第七步:根据零时刻至最小初至时刻之间的样点对应的噪音振幅值,按照公式(8)计算各个地震道的环境噪音;
第八步:根据各个地震道的环境噪音,判断各个地震道的环境噪音是否符合预设阈值要求,实现对地震数据的采集质量的实时监控。
图4为某地区同步激发采集得到的地震数据和采用现有技术进行地震数据质量监控得到的初至时刻;图5为图4的局部放大图,由图4和图5可以发现现有技术得到的初至时刻的不准确,使得初至时刻前出现同步激发产生的地震波;图6为图4对应的地震数据和采用本发明实施例进行地震数据质量监控得到的初至时刻,图7为图6的局部放大图,由图6和图7可以发现采用本发明实施例可以得到准确的初至时刻。
图8为采用现有技术得到的环境噪音水平曲线,如图8所示部分地震道由于同步激发产生的地震波参与环境噪音计算,无法真实反映环境噪音水平(图中黑框部分),图9为采用本发明实施例得到的环境噪音水平曲线,如图9所示本发明实施例能够真实反映环境噪音水平。
图10为采用现有技术得到的环境噪音水平统计柱状图,图11为采用本发明实施例得到的环境噪音水平统计柱状图,图10和图11中,横坐标表示噪音水平大小,从左至右由小到大,通过对比图10和图11可以发现图10中由于初至时刻不准确导致同步激发产生的地震波参与环境噪音水平估算,造成环境噪音水平估算失真,出现大量的高环境噪音水平道,无法真实反映环境噪音水平,而采用本发实施例可以得到准确的初至时刻,进而能够真实反映环境噪音水平。
综上所述,本发明实施例通过:根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻,能够得到准确的最小初至时刻;根据各个地震道的最小初至时刻,确定各个地震道的环境噪音,可以提高环境噪音的计算精度,真实反映地震数据采集过程中的环境噪音;根据各个地震道的环境噪音,监控地震数据的采集质量,实现了地震数据采集过程中基于环境噪音的数据采集质量实时监控,提高了数据采集质量监控效率和精度。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种同步激发采集地震数据的实时监控方法,其特征在于,包括:
根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;
获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;
根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;其中,初至时刻与炮检距之间的关联关系根据样本炮的初至时刻和炮检距关系确定,所述样本炮用于确定初至时刻与炮检距之间的关联关系;
根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻;
根据各个地震道的最小初至时刻,确定各个地震道的环境噪音;
根据各个地震道的环境噪音,监控地震数据的采集质量。
2.如权利要求1所述的方法,其特征在于,根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围,包括:按照如下方式确定炮的同步激发时刻的范围:
Tmax=T0+T;
Tmin=T0-T;
其中,T0为本炮的激发时刻;T为地震数据的记录时长;Tmax为炮的同步激发时刻的最大值,Tmin为炮的同步激发时刻的最小值。
3.如权利要求1所述的方法,其特征在于,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮,包括:
将激发时刻大于或等于炮的同步激发时刻的最小值,且小于或等于炮的同步激发时刻的最大值的炮保留;
将激发时刻小于炮的同步激发时刻的最小值,和大于炮的同步激发时刻的最大值的炮删除。
4.如权利要求1所述的方法,其特征在于,初至时刻与炮检距之间的关联关系如下:
{Tj(δxj)|j=1,…,m};
其中,j为地震道的顺序号,δxj为地震道j对应于样本炮的炮点的炮检距,Tj为地震道j对应于样本炮的初至时刻,m为地震道的数目。
5.如权利要求1所述的方法,其特征在于,根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,包括:按照如下方式确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距:
其中,i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;j为地震道的顺序号;m为地震道的数目;xi为i炮点的东坐标;yi为i炮点的北坐标;xj为j地震道的东坐标,yj为j地震道的北坐标;δxi,j为j地震道与i炮点之间的炮检距。
6.如权利要求1所述的方法,其特征在于,根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻,包括:按照如下方式确定各个地震道对应于本炮和各个同步激发炮的初至时刻:
Ti,j=Ti,j(δxi,j)+Ti-T0,i=0,...,n;j=1,...,m;
其中,T0为本炮的激发时刻;i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;i大于0时Ti为同步激发炮的激发时刻;j为地震道的顺序号;m为地震道的数目;δxi,j为j地震道与i炮点之间的炮检距;Ti,j为j地震道对应于i炮点所在炮的初至时刻。
7.如权利要求1所述的方法,其特征在于,根据各个地震道的最小初至时刻,确定各个地震道的环境噪音,包括:按照如下方式确定各个地震道的环境噪音:
其中,j为地震道的顺序号;m为地震道的数目;Aj为j地震道的环境噪音;l为j地震道内样点的顺序号,Nj为j地震道最小初至时刻对应的样点的顺序号;Al,j为j地震道l样点的振幅值。
8.一种同步激发采集地震数据的实时监控装置,其特征在于,包括:
同步激发时刻范围确定模块,用于根据地震数据采集器输出的地震记录的本炮的激发时刻和地震数据的记录时长,确定炮的同步激发时刻的范围;
同步激发炮筛选模块,用于获得多个炮的激发时刻,基于炮的同步激发时刻的范围对多个炮进行筛选,确定多个同步激发炮;
炮检距确定模块,用于根据预设的炮点坐标和地震道坐标,分别确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距;
初至时刻确定模块,用于根据本炮和各个同步激发炮的激发时刻,以及各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距,以及初至时刻与炮检距之间的关联关系,确定各个地震道对应于本炮和各个同步激发炮的初至时刻;其中,初至时刻与炮检距之间的关联关系根据样本炮的初至时刻和炮检距关系确定,所述样本炮用于确定初至时刻与炮检距之间的关联关系;
最小初至时刻确定模块,用于根据各个地震道对应于本炮和各个同步激发炮的初至时刻,确定各个地震道的最小初至时刻;
环境噪音确定模块,用于根据各个地震道的最小初至时刻,确定各个地震道的环境噪音;
监控模块,用于根据各个地震道的环境噪音,监控地震数据的采集质量。
9.如权利要求8所述的装置,其特征在于,所述同步激发炮筛选模块具体用于:
将激发时刻大于或等于炮的同步激发时刻的最小值,且小于或等于炮的同步激发时刻的最大值的炮保留;
将激发时刻小于炮的同步激发时刻的最小值,和大于炮的同步激发时刻的最大值的炮删除。
10.如权利要求8所述的装置,其特征在于,炮检距确定模块具体用于:按照如下方式确定各个地震道与本炮的炮点和各个同步激发炮的炮点之间的炮检距:
其中,i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;j为地震道的顺序号;m为地震道的数目;xi为i炮点的东坐标;yi为i炮点的北坐标;xj为j地震道的东坐标,yj为j地震道的北坐标;δxi,j为j地震道与i炮点之间的炮检距。
11.如权利要求8所述的装置,其特征在于,初至时刻确定模块具体用于:按照如下方式确定各个地震道对应于本炮和各个同步激发炮的初至时刻:
Ti,j=Ti,j(δxi,j)+Ti-T0,i=0,...,n;j=1,...,m;
其中,T0为本炮的激发时刻;i为炮点序列号,i等于0时为本炮的炮点,i大于0时为同步激发炮的炮点;n为同步激发炮的炮点数目;i大于0时Ti为同步激发炮的激发时刻;j为地震道的顺序号;m为地震道的数目;δxi,j为j地震道与i炮点之间的炮检距;Ti,j为j地震道对应于i炮点所在炮的初至时刻。
12.如权利要求8所述的装置,其特征在于,所述环境噪音确定模块具体用于:按照如下方式确定各个地震道的环境噪音:
其中,j为地震道的顺序号;m为地震道的数目;Aj为j地震道的环境噪音;l为j地震道内样点的顺序号,Nj为j地震道最小初至时刻对应的样点的顺序号;Al,j为j地震道l样点的振幅值。
13.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7任一所述方法。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至7任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010855259.XA CN112083480B (zh) | 2020-08-24 | 2020-08-24 | 同步激发采集地震数据的实时监控方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010855259.XA CN112083480B (zh) | 2020-08-24 | 2020-08-24 | 同步激发采集地震数据的实时监控方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112083480A CN112083480A (zh) | 2020-12-15 |
CN112083480B true CN112083480B (zh) | 2023-07-25 |
Family
ID=73727961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010855259.XA Active CN112083480B (zh) | 2020-08-24 | 2020-08-24 | 同步激发采集地震数据的实时监控方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112083480B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101852867A (zh) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种矢量合成检波点二次定位方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2598024A1 (en) * | 2005-02-18 | 2006-08-24 | Bp Corporation North America Inc. | System and method for using time-distance characteristics in acquisition, processing and imaging of t-csem data |
CN106814395B (zh) * | 2016-12-14 | 2019-01-18 | 中国石油天然气股份有限公司 | 确定背景炮的方法和装置 |
CN108375789B (zh) * | 2018-02-11 | 2020-10-09 | 中国石油化工股份有限公司 | 联合采集地震数据的同步匹配方法 |
CN109100795B (zh) * | 2018-06-12 | 2020-01-07 | 中国石油天然气集团有限公司 | 一种面元的炮检点布设方法、装置及系统 |
GB201818594D0 (en) * | 2018-11-14 | 2018-12-26 | Bp Exploration Operating Co Ltd | Passive seismic imaging |
CN110907991B (zh) * | 2019-12-11 | 2021-03-16 | 重庆大学 | 基于数据场势值的震源定位方法、系统及可读存储介质 |
-
2020
- 2020-08-24 CN CN202010855259.XA patent/CN112083480B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101852867A (zh) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种矢量合成检波点二次定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112083480A (zh) | 2020-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110568485B (zh) | 一种基于神经网络的多道地震连续记录分离方法 | |
CN1516814A (zh) | 处理数据的方法 | |
CN111352152B (zh) | 地震数据观测系统快速质控方法及装置 | |
CN115047417B (zh) | 一种模拟雷达回波及干扰的方法和系统 | |
CN111596347A (zh) | 一种快速求取表层纵横波速度比的方法及装置 | |
CN112083480B (zh) | 同步激发采集地震数据的实时监控方法及装置 | |
CN110703313A (zh) | 考虑传感器感度的声发射事件震级获取方法、系统及可读存储介质 | |
CN110954941A (zh) | 一种自动化初至拾取方法及系统 | |
CN111665563B (zh) | 基于聚焦分析的叠前偏移垂向分辨率评价方法 | |
CN116466396A (zh) | 一种微地震震源非线性优化定位方法 | |
CN111965729B (zh) | 可控震源组合中心实时监控方法、系统及装置 | |
CN112241020B (zh) | 稀疏地震数据采集中确定欠采样率的方法和装置 | |
CN111856558B (zh) | 一种地震记录自动增益方法 | |
CN110243939B (zh) | 一种基于机械波定位的水坝缺陷检测系统及其工作方法 | |
CN114910966B (zh) | 分方位旅行时层析反演方法及装置 | |
CN111929721A (zh) | 地震信息采集系统及方法 | |
CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
CN109684317B (zh) | 一种测月雷达探测数据的数据编辑方法 | |
CN111965702B (zh) | 一种基于地震初至的剩余静校正量处理方法及系统 | |
CN111856557A (zh) | 深度域合成地震记录制作方法及装置 | |
CN112014874B (zh) | 一种计算地震资料品质的量化评价值的方法及存储介质 | |
CN112904422B (zh) | 一种拾取炸药震源地震数据初至时间起跳点的方法及系统 | |
CN112946744B (zh) | 一种基于动态时差规整的最小二乘偏移成像方法及系统 | |
JP7114849B2 (ja) | 干渉信号での鉱山による地震のp波を識別するための方法 | |
CN112379418B (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 |