CN104977602B - 一种地震数据采集施工的控制方法及装置 - Google Patents

一种地震数据采集施工的控制方法及装置 Download PDF

Info

Publication number
CN104977602B
CN104977602B CN201510345883.4A CN201510345883A CN104977602B CN 104977602 B CN104977602 B CN 104977602B CN 201510345883 A CN201510345883 A CN 201510345883A CN 104977602 B CN104977602 B CN 104977602B
Authority
CN
China
Prior art keywords
mrow
single shot
noise
shot record
comparing 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
Application number
CN201510345883.4A
Other languages
English (en)
Other versions
CN104977602A (zh
Inventor
吴永国
张立军
刘凤智
尹吴海
马明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201510345883.4A priority Critical patent/CN104977602B/zh
Publication of CN104977602A publication Critical patent/CN104977602A/zh
Application granted granted Critical
Publication of CN104977602B publication Critical patent/CN104977602B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请实施例公开了一种地震数据采集施工的控制方法及装置。该方法包括根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级;将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果;将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果;根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。利用本申请实施例所公开的方法可以减少勘探成本。

Description

一种地震数据采集施工的控制方法及装置
技术领域
本申请涉及地震资料采集施工设计方法技术领域,特别涉及一种基于风干扰的地震数据采集施工的控制方法及装置。
背景技术
外界干扰波是地震勘探中不可避免的一类外界干扰波。它也称为环境噪声,是指没有地震激发就存在于记录中的由各种外部因素引起的干扰。外界干扰波包括油田设施、公路、铁路、工业设施等干扰源所产生的固定干扰波和高低频背景的微震、仪器噪音、50Hz交流电干扰以及风干扰等随机干扰波。
长期以来外界干扰波是地震勘探的大敌。外界干扰波的存在不但会引起地震资料信噪比的降低,也会影响地震资料处理的质量。外界干扰波中的风干扰更是无法避免,尤其是中国西部勘探区域这种高原干旱型气候,其自然条件恶劣,大风天气多(平均每日刮风的时长都在10个小时)。这种刮风天气严重影响了野外施工的工作效率。
在高密度地震勘探技术实施前,对于刮风这种恶劣天气,现有技术中通常是现场检测施工区域的刮风风级,当检测到风级达到四级以上或环境干扰大于7微伏时,直接采用停止地震采集项目施工这种控制方法。但随着高密度高效采集技术的推广应用,每天勘探成本投入很高,一旦因刮风达到四级,就停止野外施工,这会增加野外勘探成本。因此,需要一种新的地震数据采集施工的控制方法。
发明内容
本申请实施例的目的是提供一种地震数据采集施工的控制方法及装置,以提高判断是否可以在勘探区域进行地震数据采集施工的准确性。
为达到上述技术目的,本申请实施例提供了一种地震数据采集施工的控制方法,其包括:
根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级;
将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果;
将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果;
根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。
优选的,所述根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级包括:
获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的第一信噪比;
根据所获取的第一信噪比,设定所述门槛风级。
优选的,所述获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的第一信噪比包括构建下述公式来计算所述第一信噪比:
其中,SNR表示第一信噪比;E为地震道总能量;En为噪音的总能量;Es为有效信号的总能量;Rij为所选时窗内第i地震道第j个样点的振幅值,N为所选时窗内地震道的总道数,M为所选时窗内单个地震道上的采样点数,i和j均为正整数。
优选的,所述根据所获取的第一信噪比,设定所述门槛风级包括:
利用所获取的覆盖次数,计算第二信噪比;
将所获取的所有第一信噪比与所述第二信噪比进行对比,根据对比结果设定所述门槛风级。
优选的,所述将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果包括:
统计所述噪音值大于所述噪音阈值的地震道的数量;
计算所统计的地震道的数量与地震道的总数量之间的第一比值,
相应的,所述根据所述第一对比结果,控制所述目的区域中的地震数据采集施工包括:
将所述第一比值与第一预设阈值进行对比,判断是否能在所述目的区域中进行地震数据采集施工。
优选的,所述噪音阈值是通过获取不同等级风干扰背景下目的区域中预先获取的背景单炮记录中的平均噪音值来设定,或者通过获取所述门槛风级所对应的预先获取的目的区域中的背景单炮记录中的平均噪音值来设定。
优选的,所述将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果包括:
计算所述实验单炮记录中初至波峰起跳前的振幅值与初至波峰处的振幅值之间的第二比值,
相应的,所述根据所述第二对比结果,控制所述目的区域中的地震数据采集施工包括:
将所述第二比值与第二预设阈值进行对比,判断是否能在所述目的区域中进行地震数据采集施工。
优选的,在根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级前,所述方法还包括:
获取不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录。
优选的,所述获取不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录包括:
根据所获取的目的区域以往的气象资料,初步判断是否能在所述目的区域中进行地震数据采集施工;
在初步判断不能在所述目的区域中进行地震数据采集施工时,获取所述背景单炮记录和所述实验单炮记录。
本申请实施例还提供了一种地震数据采集施工的控制装置,包括:
设定单元,用于根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级;
第一对比单元,用于将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果;
第二对比单元,用于将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果;
控制单元,用于根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。
与现有技术相比,本申请实施例通过将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,以及将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比来控制所述目的区域中的地震数据采集施工,这提高了判断是否可以在勘探区域进行地震数据采集施工的准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例中所提供的一种地震数据采集施工的控制方法的流程图。
图2是实验单炮记录的第二信噪比与覆盖次数之间的关系曲线。
图3是本申请实施例中所提供的一种地震数据采集施工的控制装置的模块图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面结合附图对本申请实施例所提供的一种地震数据采集施工的控制方法进行详细说明。虽然本申请提供了如下述实施例或流程图所述的方法操作步骤,但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤,操作步骤之间的执行顺序并没有限制。如图1所示,本申请实施例所提供的一种地震数据采集施工的控制方法包括如下步骤:
S110:根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级。
所述目的区域可以是整个项目施工区域,也可以是整个项目施工区域中的部分区域。
所述实验单炮记录可以包括现场试验中放炮产生的有效信号以及风干扰产生的噪音。
所述门槛风级可以是指可以在所述目的区域进行地震数据采集施工的最高风级。
在获取到不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录后,可以根据所获取的所有实验单炮记录来设定可以在所述目的区域进行地震数据采集施工的门槛风级。
根据所获取的实验单炮记录来设定所述门槛风级可以是通过获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的信噪比来进行设定。具体的,可以包括以下两个子步骤:
S111:获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的第一信噪比。
获取每一等级风干扰所对应的所述实验单炮记录中的第一信噪比可以是从试验现场的施工仪器显示的信噪比曲线中提取;也可以是通过处理现场试验所得到的地震数据来获取。
所述通过处理现场试验所得到的地震数据来获取可以是利用能量叠加法、频谱估算法以及互相关法等方法对所获取的实验单炮记录进行计算来得到第一信噪比。其中,利用能量叠加法所得到的第一信噪比表达式可以表示如下:
上式中,SNR表示第一信噪比;E为地震道总能量;En为噪音的总能量;Es为有效信号的总能量;Rij为所选时窗内第i地震道第j个样点的振幅值,N为所选时窗内地震道的总道数,M为所选时窗内单个地震道上的采样点数,i和j均为正整数。
需要说明的是,计算第一信噪比并不限于以上方法。利用上述方法来计算第一信噪比的具体过程可以参考张军华等人发表的标题为“地震资料信噪比定量计算及方法比较”的文献。
S112:根据所获取的第一信噪比,设定可以在所述目的区域进行地震数据采集施工的门槛风级。
在得到所有等级风干扰所对应的第一信噪比后,可以根据所获取的所有第一信噪比来设定可以在所述目的区域进行地震数据采集施工的门槛风级。该步骤可以包括利用所获取的覆盖次数,计算实验单炮记录第二信噪比;将所获取的所有第一信噪比与所述第二信噪比进行对比,根据对比结果设定所述门槛风级。
在一具体实现方式中,可以获取根据滑动扫描高效采集项目来设计的观测系统的覆盖次数;然后通过构建覆盖次数与实验单炮记录的第二信噪比之间的匹配关系来计算实验单炮记录第二信噪比。
由于风干扰常表现为随机干扰的特性,可以采用统计规律来压制。因此,根据统计性原理,所构建的覆盖次数与所述第二信噪比之间的匹配关系可以表示如下:
上式中,fold为覆盖次数,Shot record(S/N)为实验单炮记录的第二信噪比,Expectsection(S/N)为预先设定的期望剖面信噪比。其中,Expect section(S/N)≥2、Expectsection(S/N)≥4或Expect section(S/N)≥8。优选的,所述期望剖面信噪比等于8。
在获取观测系统的覆盖次数后,可以将所获取的覆盖次数代入上式(2)中,计算出实验单炮记录的第二信噪比。
图2示出了一实际应用中期望剖面信噪比为8时,实验单炮记录的第二信噪比与覆盖次数之间的关系曲线。
在得到第二信噪比后,可以将所获取的所有第一信噪比与所述第二信噪比进行对比,找出与所述第二信噪比相同或最接近的第一信噪比,在找出所述第一信噪比后,可以将所述第一信噪比所对应的风级设定为门槛风级。
S120:将所述实验单炮记录中对应所述门槛风级的噪音值与所设定的噪音阈值进行对比,得到第一对比结果。
在确定出门槛风级后,可以从所获取的所有实验单炮记录中查找与所述门槛风级对应的实验单炮记录,然后从所找到的实验单炮记录中提取所有地震道的噪音值,再将所提取出的每个噪音值与所设定的噪音阈值进行对比,得到第一对比结果。具体的可以包括:
S121:设定可以在所述目的区域进行地震数据采集施工的噪音阈值。
可以通过获取不同等级风干扰背景下目的区域中预先获取的背景单炮记录中的平均噪音值来设定所述噪音阈值;也可以通过获取所述门槛风级所对应的预先获取的目的区域中的背景单炮记录中的平均噪音值来设定所述噪音阈值。所述背景单炮记录可以是在试验现场制造不同等级的风干扰,然后检测所述风干扰所产生的噪音来得到的。
所述通过获取不同等级风干扰背景下目的区域中预先获取的背景单炮记录中的平均噪音值来设定所述噪音阈值可以包括以下两个子步骤:
(1)从每个等级风干扰背景下的背景单炮记录中获取所有地震道的平均噪音值。
在一实施例中,可以首先从某一等级风干扰下所获取的背景单炮记录中提取出所有地震道的噪音值;然后计算所提取出的所有噪音值的平均噪音值;再按照获取该等级风干扰所对应的平均噪音值的方法,依次获取其他等级风干扰所对应的平均噪音值。
在另一实施例中,可以首先从所获取的所有背景单炮记录中提取出所有地震道的噪音值;然后再依次计算出每个等级风干扰所对应的所有噪音值的平均噪音值。
计算每一等级风干扰所对应的平均噪音值的公式可以表示如下:
其中,Average_Noise为平均噪音值(单位可以为微伏),R(i,j)为背景单炮记录中第i地震道第j个样点的噪音值,N’为背景单炮记录中地震道的总道数,M’为背景单炮记录中单个地震道上的采样点数。
在得到平均噪音值后,可以绘制不同等级风干扰所对应的平均噪音值曲线,还可以将所绘制的平均噪音值曲线进行显示,以便于仪器现场的操作员监控噪音。
(2)根据所得到的平均噪音值,设定所述噪音阈值。
在得到所有风级(即风干扰的等级)所对应的所有平均噪音值后,从所有平均噪音值中提取与所述门槛风级所对应的平均噪音值,可以将所提取的平均噪音值作为所述噪音阈值。此外,还可以根据所提取的平均噪音值,来设定所述噪音阈值,此时,所述噪音阈值与所提取的平均噪音值可以有所不同。所设定的噪音阈值可以为7微伏。
所述通过获取所述门槛风级所对应的预先获取的目的区域中的背景单炮记录中的平均噪音值来设定所述噪音阈值可以是从所述门槛风级所对应的背景单炮记录中获取所有地震道的平均噪音值,然后根据该平均噪音值来设定所述噪音阈值。具体的执行过程可以参考上述步骤(1)-(2),在此不再赘叙。
S122:从所述实验单炮记录中获取对应所述门槛风级的所有地震道的噪音值。
从所获取的所有实验单炮记录中查找与所述门槛风级对应的实验单炮记录,然后从所找到的实验单炮记录中提取所有地震道的噪音值。
S123:将所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果。
将所获取的每个噪音值与所设定的噪音阈值进行对比,统计噪音值大于所述噪音阈值的地震道的数量,计算所统计出的地震道的数量与地震道的总数量之间的第一比值,所述第一比值可以作为第一对比结果。
所述第一比值可以用百分比的形式表示,其可以称为噪音超标值百分比。计算所述第一比值的公式可以表示如下:
上式中,Noise_Ratio为第一比值,Total_Number为实验单炮记录中地震道总道数,Ultra_Number为实验单炮记录中噪音值大于所述噪音阈值的地震道数。
S130:将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果。
在确定门槛风级后,可以从所获取的所有实验单炮记录中提取出对应所述门槛风级的实验单炮记录,然后从所提取出的实验单炮记录中获取初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量,再将所提取出的第一地震波能量与第二地震波能量进行对比,得到第二对比结果。
所述第一地震波能量包括初至波峰处的有效波能量和风干扰产生的噪音能量。所述第二地震波能量可以是风干扰产生的噪音能量。
所述从所提取出的实验单炮记录中获取初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量可以是在含有初至的时窗内,从频谱上提取初至波峰处的地震波振幅值以及初至波峰起跳前(即初至波峰所对应样点的前一样点)的地震波振幅值。
所述将所提取出的第一地震波能量与第二地震波能量进行对比,得到第二对比结果可以是计算第一地震波能量中的有效波能量与初至波峰起跳前的第二地震波能量之间的第二比值,所述第二比值可以作为第二对比结果。所述第二比值也可以用百分比的形式表示,其可以称为初至能量比。计算所述第二比值的公式可以表示如下:
上式中,Firstarrival_ratio为第二比值,Rf(j)为实验单炮记录中初至波峰处第j个样点的地震波振幅值,Rn(j)为实验单炮记录中初至波峰起跳前第j个样点的地震波振幅值,L为时窗中的样点数。
需要说明的是,步骤S120和S130之间的执行顺序并没有限制。
S140:根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。
在得到所述第一对比结果和所述第二对比结果后,可以将所得到的第一对比结果和第一预设阈值进行对比,以及将所得到的第二对比结果和第二预设阈值进行对比,根据二者的对比结果判断是否可以在所述目的区域中进行地震数据采集施工。此外,还可以控制施工的时间。
在对比出第一对比结果小于第一预设阈值,第二对比结果小于第二预设阈值时,可以判断可以在所述目的区域中进行地震数据采集施工;在对比出第一对比结果大于第一预设阈值,第二对比结果小于第二预设阈值时,也可以判断可以在所述目的区域中进行地震数据采集施工;在对比出第一对比结果大于第一预设阈值,第二对比结果大于第二预设阈值时,则可以判断不可以在所述目的区域中进行地震数据采集施工。
所述第一预设阈值和所述第二预设阈值可以是根据相关标准来设定的,例如《陆上石油地质勘探资采集技术规范-SY/T 5314-2011》标准。所述第一预设阈值可以为16.7%,所述第二预设阈值可以为50%。
在另一实施例中,在步骤S110之前,该方法还可以包括:
S100:获取不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录。
所述获取在不同等级风干扰下目的区域中的背景单炮记录和实验单炮记录可以是基于触发条件,即时或预先获取在不同等级风干扰下目的区域中的背景单炮记录和实验单炮记录;也可以是根据所获取的目的区域以往的气象资料,初步判断是否能在所述目的区域进行地震数据采集施工;在初步判断不能在所述目的区域中进行地震数据采集施工时,获取在不同等级风干扰下目的区域中的所述背景单炮记录和所述实验单炮记录。
所述触发条件可以是现场施工人员的指示,也可以是收到的自动化指令。
所述气象资料可以包括所述目的区域近几年的刮风风级和/或风速,例如近三年每月每天的刮风风级和风速。所述气象资料可以是预先或即时向气象局获取,也可以是收集预先获取的气象资料。
所述根据所获取的目的区域以往的气象资料,初步判断所述目的区域是否能进行地震数据采集施工具体的可以是根据所获取的以往气象资料,统计分析所述目的区域近几年的刮风风级和/或风速,绘制风级和/或风速曲线;然后利用所绘制的风级和/或风速曲线,并根据项目施工周期,预估刮风对所述目的区域的地震数据采集施工的影响。根据预估的结果,在项目施工周期内没有出现超过预设风级(例如四级风)的天气时,可以判断不需要考虑刮风对地震数据采集项目施工的影响,在所述目的区域可以进行地震数据采集项目施工;在项目施工周期内出现超过预设风级(例如四级风)的天气时,可以初步判断不能在所述目的区域中进行地震数据采集施工。
在一具体实现方式中,所述获取在不同等级风干扰下目的区域中的所述背景单炮记录和所述实验单炮记录可以是在项目采集初期,在项目现场进行试验,录制不同等级风干扰下的背景单炮记录和实验单炮记录。具体的可以是:首先在项目现场制造某一风级(例如第一风级)的背景,录制该风级背景下的背景单炮记录;然后可以在较短时间内或录制背景单炮记录的同时,在项目现场的预设炮点放炮,录制放炮所产生的实验单炮记录;再依次录制其他风级背景下的背景单炮记录和实验单炮记录。此外,在录制实验单炮记录时,还可以从相关监控仪器中获取放炮噪音能量。
通过上述步骤可以看出,本申请实施例通过将所述实验单炮记录中对应所述门槛风级的噪音值与所设定的噪音阈值进行对比以及将所述实验单炮记录中受到风干扰前的初至能量与受到所述门槛风级的风干扰后的初至能量进行对比来控制所述目的区域中的地震数据采集施工,而不是仅通过检测风级和风速来控制目的区域中的地震数据采集施工,这可以提高判断是否可以在勘探区域进行地震数据采集施工的准确性,进而可以缩短施工周期以及减少勘探成本。例如,现有技术中在刮风四级以上,就停止施工,而利用本申请实施例所提供的技术方案可以在刮风大于四级时,还可以继续施工,这可以缩短施工周期以及减少勘探成本。
本申请实施例还提供了一种地震数据采集施工的控制装置,如图3所示。该装置包括:
设定单元310,用于根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级;
第一对比单元320,用于将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果;
第二对比单元330,用于将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果;
控制单元340,用于根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。
在一实施例中,设定单元310可以包括(图中未示出):
获取子单元,用于获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的第一信噪比;
设定子单元,用于根据所获取的第一信噪比,设定所述门槛风级。
在一实施例中,第一对比单元320可以包括(图中未示出):
统计子单元,用于统计所述噪音值大于所述噪音阈值的地震道的数量;
计算子单元,用于计算所统计的地震道的数量与地震道的总数量之间的第一比值,
相应的,所述控制单元340包括:
第一对比判断子单元,用于将所述第一比值与第一预设阈值进行对比,判断是否能在所述目的区域中进行地震数据采集施工。
在一实施例中,第二对比单元330可以包括(图中未示出):
计算子单元,用于计算所述实验单炮记录中初至波峰起跳前的振幅值与初至波峰处的振幅值之间的第二比值,
相应的,所述述控制单元340包括:
第二对比判断子单元,用于将所述第二比值与第二预设阈值进行对比,判断是否能在所述目的区域中进行地震数据采集施工。
在一实施例中,该装置还可以包括获取单元300,其用于获取不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (11)

1.一种地震数据采集施工的控制方法,其特征在于,包括:
根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级;
将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果;
将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果;
根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。
2.根据权利要求1所述的方法,其特征在于,所述根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级包括:
获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的第一信噪比;
根据所获取的第一信噪比,设定所述门槛风级。
3.根据权利要求2所述的方法,其特征在于,所述获取所有等级风干扰所对应的所述实验单炮记录中的噪音与有效信号之间的第一信噪比包括构建下述公式来计算所述第一信噪比:
<mrow> <mi>S</mi> <mi>N</mi> <mi>R</mi> <mo>=</mo> <mfrac> <msub> <mi>E</mi> <mi>s</mi> </msub> <msub> <mi>E</mi> <mi>n</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>E</mi> <mi>s</mi> </msub> <mrow> <mi>E</mi> <mo>-</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mi>N</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>R</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msup> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
其中,SNR表示第一信噪比;E为地震道总能量;En为噪音的总能量;Es为有效信号的总能量;Rij为所选时窗内第i地震道第j个样点的振幅值,N为所选时窗内地震道的总道数,M为所选时窗内单个地震道上的采样点数,i和j均为正整数。
4.根据权利要求2所述的方法,其特征在于,所述根据所获取的第一信噪比,设定所述门槛风级包括:
利用所获取的覆盖次数,计算第二信噪比;
将所获取的所有第一信噪比与所述第二信噪比进行对比,根据对比结果设定所述门槛风级。
5.根据权利要求1所述的方法,其特征在于,所述将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果包括:
统计所述噪音值大于所述噪音阈值的地震道的数量;
计算所统计的地震道的数量与地震道的总数量之间的第一比值,
相应的,所述根据所述第一对比结果,控制所述目的区域中的地震数据采集施工包括:
将所述第一比值与第一预设阈值进行对比,判断是否能在所述目的区域中进行地震数据采集施工。
6.根据权利要求1或5所述的方法,其特征在于,所述噪音阈值是通过获取不同等级风干扰背景下目的区域中预先获取的背景单炮记录中的平均噪音值来设定,或者通过获取所述门槛风级所对应的预先获取的目的区域中的背景单炮记录中的平均噪音值来设定。
7.根据权利要求1所述的方法,其特征在于,所述将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果包括:
计算所述实验单炮记录中初至波峰前一样点的振幅值与初至波峰处样点的振幅值之间的第二比值,
相应的,所述根据所述第二对比结果,控制所述目的区域中的地震数据采集施工包括:
将所述第二比值与第二预设阈值进行对比,判断是否能在所述目的区域中进行地震数据采集施工。
8.根据权利要求1所述的方法,其特征在于,所述根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工包括:
将所述第一对比结果与第一预设阈值进行对比;
将所述第二对比结果与第二预设阈值进行对比;
在对比出所述第一对比结果小于第一预设阈值,所述第二对比结果小于第二预设阈值时,判断能在所述目的区域中进行地震数据采集施工,或
在对比出所述第一对比结果大于第一预设阈值,所述第二对比结果小于第二预设阈值时,判断能在所述目的区域中进行地震数据采集施工。
9.根据权利要求1所述的方法,其特征在于,在根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级前,所述方法还包括:
获取不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录。
10.根据权利要求9所述的方法,其特征在于,所述获取不同等级风干扰背景下目的区域中的背景单炮记录和实验单炮记录包括:
根据所获取的目的区域以往的气象资料,初步判断是否能在所述目的区域中进行地震数据采集施工;
在初步判断不能在所述目的区域中进行地震数据采集施工时,获取所述背景单炮记录和所述实验单炮记录。
11.一种地震数据采集施工的控制装置,其特征在于,包括:
设定单元,用于根据所获取的不同等级风干扰背景下目的区域中的实验单炮记录,设定所述目的区域进行地震数据采集施工的门槛风级;
第一对比单元,用于将从所述门槛风级所对应的所述实验单炮记录中所获取的噪音值与所设定的噪音阈值进行对比,得到第一对比结果;
第二对比单元,用于将从所述门槛风级所对应的所述实验单炮记录中所获取的初至波峰处的第一地震波能量以及初至波峰前的第二地震波能量进行对比,得到第二对比结果;
控制单元,用于根据所述第一对比结果和所述第二对比结果,控制所述目的区域中的地震数据采集施工。
CN201510345883.4A 2015-06-19 2015-06-19 一种地震数据采集施工的控制方法及装置 Active CN104977602B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510345883.4A CN104977602B (zh) 2015-06-19 2015-06-19 一种地震数据采集施工的控制方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510345883.4A CN104977602B (zh) 2015-06-19 2015-06-19 一种地震数据采集施工的控制方法及装置

Publications (2)

Publication Number Publication Date
CN104977602A CN104977602A (zh) 2015-10-14
CN104977602B true CN104977602B (zh) 2017-09-01

Family

ID=54274277

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510345883.4A Active CN104977602B (zh) 2015-06-19 2015-06-19 一种地震数据采集施工的控制方法及装置

Country Status (1)

Country Link
CN (1) CN104977602B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105467428A (zh) * 2015-11-17 2016-04-06 南京航空航天大学 一种基于短时能量检测和频谱特征分析的地震波预警方法
CN110673210B (zh) * 2019-10-13 2021-06-04 东北石油大学 一种地震原始数据信噪比定量分析评价方法
CN112379421A (zh) * 2020-10-30 2021-02-19 中国石油天然气集团有限公司 识别地震勘探数据中的掉排列单炮的方法和装置
CN116046079B (zh) * 2023-04-03 2023-06-30 山东省地质调查院(山东省自然资源厅矿产勘查技术指导中心) 一种基于位置的地质环境专题数据集成管理系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2873823B1 (fr) * 2004-08-02 2006-09-15 Inst Francais Du Petrole Methode pour construire un modele d'un milieu heterogene decrit par plusieurs parametres a partir de donnees exprimees dans des echelles de temps differentes
CN100485412C (zh) * 2006-03-01 2009-05-06 北京中油天强勘探技术有限公司 地震勘探野外采集过程的控制方法
CN101551463B (zh) * 2008-04-03 2011-08-03 中国石油天然气集团公司 三维观测系统噪声压制估算方法

Also Published As

Publication number Publication date
CN104977602A (zh) 2015-10-14

Similar Documents

Publication Publication Date Title
CN104977602B (zh) 一种地震数据采集施工的控制方法及装置
Dai et al. The application of back‐propagation neural network to automatic picking seismic arrivals from single‐component recordings
CN110454687A (zh) 一种基于改进vmd的管道多点泄漏定位方法
CN105527650B (zh) 一种工程尺度下微震信号及p波初至自动识别算法
CN106199703B (zh) 一种微震震源自动定位及可靠性综合评价方法
CN106382981B (zh) 一种单站次声波信号识别提取方法
CN110196448B (zh) 一种滑坡次声信号识别方法
CN104535646B (zh) 一种粮食籽粒不完善度检测方法
CN109085663A (zh) 一种致密砂岩储层层理缝识别方法
CN106646637A (zh) 一种去除核磁信号中尖峰噪声的方法
CN103336305B (zh) 一种基于灰色理论划分致密砂岩储层岩石物理相的方法
CN105116442A (zh) 岩性油气藏弱反射地震信号的重构方法
CN106772646A (zh) 一种地面核磁共振信号提取方法
CN112364296B (zh) 一种基于深度学习的p波到时自动拾取方法
CN105204065B (zh) 一种拾取初至波的方法和装置
CN103344989A (zh) 可控震源地震记录中脉冲噪声干扰的分析方法
CN104316167B (zh) 一种低信噪比下结构响应异常检测方法
CN110020604A (zh) 一种蔬菜品质检测方法及系统
WO2019205725A1 (zh) 一种基于电成像测井分形维数的火山岩岩性识别方法
CN109541031A (zh) 基于声学和振动特性的水果硬度检测方法
CN105022041B (zh) 基于经验模态分解特征库的拖引干扰子类型识别方法
CN103886868B (zh) 冲击声检测方法及检测系统
CN108254788A (zh) 一种地震波初至拾取方法及系统
CN109409216B (zh) 基于子载波动态选择的速度自适应室内人体检测方法
CN107632326A (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