CN111965700A - 光纤声波传感地震数据中零波数噪声的剔除方法及系统 - Google Patents

光纤声波传感地震数据中零波数噪声的剔除方法及系统 Download PDF

Info

Publication number
CN111965700A
CN111965700A CN201910418078.8A CN201910418078A CN111965700A CN 111965700 A CN111965700 A CN 111965700A CN 201910418078 A CN201910418078 A CN 201910418078A CN 111965700 A CN111965700 A CN 111965700A
Authority
CN
China
Prior art keywords
amplitude
seismic
zero
noise
sampling point
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.)
Pending
Application number
CN201910418078.8A
Other languages
English (en)
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 CN201910418078.8A priority Critical patent/CN111965700A/zh
Publication of CN111965700A publication Critical patent/CN111965700A/zh
Pending legal-status Critical Current

Links

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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364—Seismic filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/32—Noise reduction
    • G01V2210/324—Filtering

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (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

光纤声波传感地震数据中零波数噪声的剔除方法及系统
技术领域
本申请涉及地球物理勘探方法，更具体的，涉及一种光纤声波传感地震数据中零波数噪声的剔除方法及系统。
背景技术
最新的光纤声波传感技术目前已经在地震勘探中逐步得到利用，相对于常规检波器采集来说具有高密度采样和耐高温高压等特点。光纤声波传感技术主要利用激光在光纤中传播时产生的瑞利逆向散射，当光纤周围的介质有振动时，瑞利散射也会发生相应的变化，通过连续观测和对这种光信号的解调就可以得到地层的地震波信息。
光纤声波传感比较适合在井中观测地震波，但由于其本身的特点导致由光信号解调得到的地震波信号具有较强的零波数噪声，即在某一时刻所有通道输出的数据产生了零点漂移，这严重影响地震波的接收效果，如何有效去除这种强干扰是光纤声波传感在地震勘探应用的一个重点课题。调研表明，目前一般基于波数域滤波或均值滤波等方式，可以取得一定的应用效果，但不能完全剔除这种干扰，或者会伤害到一些有效信号，现存方法在实际应用中还存在许多不足之处，如何准确高效地压制这种光纤采集时的零波数噪声的需求十分迫切。
发明内容
为了解决上述不足，本申请提供一种光纤声波传感地震数据中零波数噪声的剔除方法及系统。
本申请第一方面实施例提供一种光纤声波传感地震数据中零波数噪声的剔除方法，包括：
从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
在某些实施例中，所述从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱包括：
从地震波场数据中筛选出均方根振幅低于所有均方根振幅中值的地震道组成振幅谱。
在某些实施例中，所述对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，包括：
对每一时间样点对应的所有地震道的振幅按照幅值从小到大进行排序，去除位于第一设定位置之前以及位于第二设定位置之后的地震道；
对每一时间样点对应去除后的地震道的振幅进行平均处理，输出每一时间样点的平均振幅；其中每一时间样点的平均振幅为所述零波数噪声在每一时间样点的振幅。
在某些实施例中，所述第一设定位置位于排序的中间位置之前，所述第二设定位置位于排序的中间位置之后。
在某些实施例中，基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据，包括：
用所述地震波场数据中每一时间样点的原始振幅减去对应时间样点的零波数噪声的振幅，得到去除零波数噪声的地震波场数据。
在某些实施例中，还包括：
设置所述时间样点。
在某些实施例中，所述设置所述时间样点，包括：
设置炮点接收点的位置坐标；
对所述位置坐标进行深度域向时间域的转换，形成每一个时间样点的时间位置坐标。
本申请第二方面实施例提供一种光纤声波传感地震数据中零波数噪声的剔除系统，包括：
振幅谱获取模块，从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
均值滤波处理模块，对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
零波数噪声剔除模块，基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
在某些实施例中，所述振幅谱获取模块从地震波场数据中筛选出均方根振幅低于所有均方根振幅中值的地震道组成振幅谱。
在某些实施例中，所述均值滤波处理模块，包括：
排序单元，对每一时间样点对应的所有地震道的振幅按照幅值从小到大进行排序，去除位于第一设定位置之前以及位于第二设定位置之后的地震道；
平均处理单元，对每一时间样点对应去除后的地震道的振幅进行平均处理，输出每一时间样点的平均振幅；其中每一时间样点的平均振幅为所述零波数噪声在每一时间样点的振幅。
在某些实施例中，所述第一设定位置位于排序的中间位置之前，所述第二设定位置位于排序的中间位置之后。
在某些实施例中，零波数噪声剔除模块用所述地震波场数据中每一时间样点的原始振幅减去对应时间样点的零波数噪声的振幅，得到去除零波数噪声的地震波场数据
在某些实施例中，还包括：
时间样点设置模块，设置所述时间样点。
在某些实施例中，所述时间样点设置模块包括：
接收点位置设置单元，设置炮点接收点的位置坐标；
时深转换单元，对所述位置坐标进行深度域向时间域的转换，形成每一个时间样点的时间位置坐标。
本申请第三方面提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述所述的光纤声波传感地震数据中零波数噪声的剔除方法的步骤。
本申请第四方面实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述所述的光纤声波传感地震数据中零波数噪声的剔除方法的步骤。
本申请的有益效果如下：
本申请提供一种光纤声波传感地震数据中零波数噪声的剔除方法及系统，首先利用光纤声波传感接收仪器采集得到地震波场数据，并对该数据进行处理，标注出统计均方根振幅较小的地震道形成振幅谱，并利用均值滤波方法预测零波数干扰的每个时间样点的振幅值，通过减去法达到了压制零波数干扰并突出有效信息及提高信噪比的目的，通过本方法处理后光缆零波数噪声可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本申请实施例中一种光纤声波传感地震数据中零波数噪声的剔除方法的流程示意图之一。
图2示出本申请实施例中光纤声波传感采集仪器采集的原始井中地震单炮记录。
图3示出本申请实施例中去除零波数噪声后的光纤采集地震数据。
图4示出本申请实施例中一种光纤声波传感地震数据中零波数噪声的剔除方法的流程示意图之二。
图5示出本申请实施例中图1中步骤S200的具体流程示意图。
图6示出本申请实施例中一种光纤声波传感地震数据中零波数噪声的剔除方法的流程示意图之三。
图7示出本申请实施例中一种光纤声波传感地震数据中零波数噪声的剔除方法的流程示意图之四。
图8示出本申请实施例中图7中步骤S001的具体流程示意图。
图9示出本申请实施例中一种光纤声波传感地震数据中零波数噪声的剔除系统结构示意图之一。
图10示出本申请实施例中图9中均值滤波处理模块200的具体结构示意图。
图11示出本申请实施例中一种光纤声波传感地震数据中零波数噪声的剔除系统结构示意图之二。
图12示出本申请实施例中图1中时间样点设置模块001的具体结构示意图。
图13示出适于用来实现本申请实施例的电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前一般基于波数域滤波或均值滤波等方式，可以取得一定的应用效果，但不能完全剔除这种干扰，或者会伤害到一些有效信号，现存方法在实际应用中还存在许多不足之处，如何准确高效地压制这种光纤采集时的零波数噪声的需求十分迫切。
有鉴于此，本申请提供一种光纤声波传感地震数据中零波数噪声的剔除方法及系统，首先利用光纤声波传感接收仪器采集得到地震波场数据，并对该数据进行处理，标注出统计均方根振幅较小的地震道形成振幅谱，并利用均值滤波方法预测零波数干扰的每个时间样点的振幅值，通过减去法达到了压制零波数干扰并突出有效信息及提高信噪比的目的。
下面结合附图对本申请实施例进行详细说明。
本申请第一方面实施例提供一种光纤声波传感地震数据中零波数噪声的剔除方法，如图1所示，包括：
S100:从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
S200:对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
S300:基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
本方面提供的光纤声波传感地震数据中零波数噪声的剔除方法，首先利用光纤声波传感接收仪器采集得到地震波场数据，并对该数据进行处理，标注出统计均方根振幅较小的地震道形成振幅谱，并利用均值滤波方法预测零波数干扰的每个时间样点的振幅值，通过减去法达到了压制零波数干扰并突出有效信息及提高信噪比的目的，通过本方法处理后光缆零波数噪声可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
设定幅值一般为小于所有均方根幅值的中值，该实施例中，如图4所示，步骤S100具体包括：
从地震波场数据中筛选出均方根振幅低于所有均方根振幅中值的地震道组成振幅谱。
进一步的，均值滤波处理即需要去除误差较大的样本，保留样本误差较小的样本，如图5所示，即步骤S200具体包括：
S201:对每一时间样点对应的所有地震道的振幅按照幅值从小到大进行排序，去除位于第一设定位置之前以及位于第二设定位置之后的地震道；
S202:对每一时间样点对应去除后的地震道的振幅进行平均处理，输出每一时间样点的平均振幅；其中每一时间样点的平均振幅为所述零波数噪声在每一时间样点的振幅。
该实施例中，首先剔除幅值较小的(位于第一设定位置之前)以及幅值较大的(位于第二设定位置之后)的地震道，然后对一个时间样点内的剩下的地震道取平均，进而进行了去尾均值滤波的处理。
所述第一设定位置位于排序的中间位置之前，所述第二设定位置位于排序的中间位置之后。一般第一设定位置与第二设定位置分别处于排序的1/4和3/4上，这样即能够保证时间样点中剩余的地震道相差较小，另一方面能够有足够的地震道进行取平均，提高精确度。
如图6所示，上述的步骤S300可以通过减去法实现，即用所述地震波场数据中每一时间样点的原始振幅减去对应时间样点的零波数噪声的振幅，得到去除零波数噪声的地震波场数据。
此外，时间样点可以是预先提供的，也可以是处于本方案的方法步骤中，在一实施例中，如图7所示，本申请的方法还包括：
S001:设置所述时间样点。
在一实施例中，如图8所示，步骤S001包括：
S001a：设置炮点接收点的位置坐标；
S001b：对所述位置坐标进行深度域向时间域的转换，形成每一个时间样点的时间位置坐标。
深度域向时间域的转换即为时深转换，本申请在此不予赘述。
结合附图2和3所示，图2为光纤声波传感采集仪器采集的原始井中地震单炮记录，可以看到明显的零波数噪声。图3为去除零波数噪声后的光纤采集地震数据，可以看出零波数噪声得到明显压制。
由此可知，本方面提供的光纤声波传感地震数据中零波数噪声的剔除方法，首先利用光纤声波传感接收仪器采集得到地震波场数据，并对该数据进行处理，标注出统计均方根振幅较小的地震道形成振幅谱，并利用均值滤波方法预测零波数干扰的每个时间样点的振幅值，通过减去法达到了压制零波数干扰并突出有效信息及提高信噪比的目的，通过本方法处理后光缆零波数噪声可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
本申请第二方面实施例还提供一种光纤声波传感地震数据中零波数噪声的剔除系统，如图9所示，包括：
振幅谱获取模块100，从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
均值滤波处理模块200，对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
零波数噪声剔除模块300，基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
本方面提供的光纤声波传感地震数据中零波数噪声的剔除系统，首先利用光纤声波传感接收仪器采集得到地震波场数据，并对该数据进行处理，标注出统计均方根振幅较小的地震道形成振幅谱，并利用均值滤波方法预测零波数干扰的每个时间样点的振幅值，通过减去法达到了压制零波数干扰并突出有效信息及提高信噪比的目的，通过本方法处理后光缆零波数噪声可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
设定幅值一般为小于所有均方根幅值的中值，该实施例中，所述振幅谱获取模块从地震波场数据中筛选出均方根振幅低于所有均方根振幅中值的地震道组成振幅谱。
进一步的，均值滤波处理即需要去除误差较大的样本，保留样本误差较小的样本，如图10所示，所述均值滤波处理模块200，包括：
排序单元201，对每一时间样点对应的所有地震道的振幅按照幅值从小到大进行排序，去除位于第一设定位置之前以及位于第二设定位置之后的地震道；
平均处理单元202，对每一时间样点对应去除后的地震道的振幅进行平均处理，输出每一时间样点的平均振幅；其中每一时间样点的平均振幅为所述零波数噪声在每一时间样点的振幅。
该实施例中，首先剔除幅值较小的(位于第一设定位置之前)以及幅值较大的(位于第二设定位置之后)的地震道，然后对一个时间样点内的剩下的地震道取平均，进而进行了去尾均值滤波的处理。
所述第一设定位置位于排序的中间位置之前，所述第二设定位置位于排序的中间位置之后。一般第一设定位置与第二设定位置分别处于排序的1/4和3/4上，这样即能够保证时间样点中剩余的地震道相差较小，另一方面能够有足够的地震道进行取平均，提高精确度。
在一实施例中，零波数噪声剔除模块用所述地震波场数据中每一时间样点的原始振幅减去对应时间样点的零波数噪声的振幅，得到去除零波数噪声的地震波场数据。
此外，时间样点可以是预先提供的，也可以是由本申请中的功能模块获得，在一实施例中，如图11所示，还包括：
时间样点设置模块001，设置所述时间样点。
在一实施例中，如图12所示，所述时间样点设置模块001包括：
接收点位置设置单元001a，设置炮点接收点的位置坐标；
时深转换单元001b，对所述位置坐标进行深度域向时间域的转换，形成每一个时间样点的时间位置坐标。
深度域向时间域的转换即为时深转换，本申请在此不予赘述。
结合附图2和3所示，图2为光纤声波传感采集仪器采集的原始井中地震单炮记录，可以看到明显的零波数噪声。图3为去除零波数噪声后的光纤采集地震数据，可以看出零波数噪声得到明显压制。
由此可知，本方面提供的光纤声波传感地震数据中零波数噪声的剔除系统，首先利用光纤声波传感接收仪器采集得到地震波场数据，并对该数据进行处理，标注出统计均方根振幅较小的地震道形成振幅谱，并利用均值滤波方法预测零波数干扰的每个时间样点的振幅值，通过减去法达到了压制零波数干扰并突出有效信息及提高信噪比的目的，通过本方法处理后光缆零波数噪声可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种电子设备的具体实施方式，参见图13，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；
所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
S100:从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
S200:对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
S300:基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
从上述描述可知，本申请提供的电子设备可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
S100:从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
S200:对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
S300:基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
从上述描述可知，本申请提供的计算机可读存储介质可以得到很大程度的压制和剔除，资料品质明显提升，为后续的地震数据处理解释提供了保障。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。

Claims (16)

1.一种光纤声波传感地震数据中零波数噪声的剔除方法，其特征在于，包括：
从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
2.根据权利要求1所述的零波数噪声的剔除方法，其特征在于，所述从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱包括：
从地震波场数据中筛选出均方根振幅低于所有均方根振幅中值的地震道组成振幅谱。
3.根据权利要求1所述的零波数噪声的剔除方法，其特征在于，所述对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，包括：
对每一时间样点对应的所有地震道的振幅按照幅值从小到大进行排序，去除位于第一设定位置之前以及位于第二设定位置之后的地震道；
对每一时间样点对应去除后的地震道的振幅进行平均处理，输出每一时间样点的平均振幅；其中每一时间样点的平均振幅为所述零波数噪声在每一时间样点的振幅。
4.根据权利要求3所述的零波数噪声的剔除方法，其特征在于，所述第一设定位置位于排序的中间位置之前，所述第二设定位置位于排序的中间位置之后。
5.根据权利要求1所述的零波数噪声的剔除方法，其特征在于，基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据，包括：
用所述地震波场数据中每一时间样点的原始振幅减去对应时间样点的零波数噪声的振幅，得到去除零波数噪声的地震波场数据。
6.根据权利要求1所述的零波数噪声的剔除方法，其特征在于，还包括：
设置所述时间样点。
7.根据权利要求6所述的零波数噪声的剔除方法，其特征在于，所述设置所述时间样点，包括：
设置炮点接收点的位置坐标；
对所述位置坐标进行深度域向时间域的转换，形成每一个时间样点的时间位置坐标。
8.一种光纤声波传感地震数据中零波数噪声的剔除系统，其特征在于，包括：
振幅谱获取模块，从地震波场数据中筛选出均方根振幅低于设定幅值的地震道，组成振幅谱；其中所述地震波场数据通过光纤声波传感采集仪器采集得到；
均值滤波处理模块，对振幅谱中的每一时间样点对应的所有地震道的振幅进行均值滤波处理，得到零波数噪声在每一时间样点的振幅；
零波数噪声剔除模块，基于所述地震波场数据中每一时间样点的原始振幅以及零波数噪声在每一时间样点的振幅，得到去除零波数噪声的地震波场数据。
9.根据权利要求8所述的剔除系统，其特征在于，所述振幅谱获取模块从地震波场数据中筛选出均方根振幅低于所有均方根振幅中值的地震道组成振幅谱。
10.根据权利要求8所述的剔除系统，其特征在于，所述均值滤波处理模块，包括：
排序单元，对每一时间样点对应的所有地震道的振幅按照幅值从小到大进行排序，去除位于第一设定位置之前以及位于第二设定位置之后的地震道；
平均处理单元，对每一时间样点对应去除后的地震道的振幅进行平均处理，输出每一时间样点的平均振幅；其中每一时间样点的平均振幅为所述零波数噪声在每一时间样点的振幅。
11.根据权利要求10所述的剔除系统，其特征在于，所述第一设定位置位于排序的中间位置之前，所述第二设定位置位于排序的中间位置之后。
12.根据权利要求8所述的剔除系统，其特征在于，零波数噪声剔除模块用所述地震波场数据中每一时间样点的原始振幅减去对应时间样点的零波数噪声的振幅，得到去除零波数噪声的地震波场数据。
13.根据权利要求8所述的剔除系统，其特征在于，还包括：
时间样点设置模块，设置所述时间样点。
14.根据权利要求13所述的剔除系统，其特征在于，所述时间样点设置模块包括：
接收点位置设置单元，设置炮点接收点的位置坐标；
时深转换单元，对所述位置坐标进行深度域向时间域的转换，形成每一个时间样点的时间位置坐标。
15.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至7任一项所述的光纤声波传感地震数据中零波数噪声的剔除方法的步骤。
16.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至7任一项所述的光纤声波传感地震数据中零波数噪声的剔除方法的步骤。
CN201910418078.8A 2019-05-20 2019-05-20 光纤声波传感地震数据中零波数噪声的剔除方法及系统 Pending CN111965700A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910418078.8A CN111965700A (zh) 2019-05-20 2019-05-20 光纤声波传感地震数据中零波数噪声的剔除方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910418078.8A CN111965700A (zh) 2019-05-20 2019-05-20 光纤声波传感地震数据中零波数噪声的剔除方法及系统

Publications (1)

Publication Number Publication Date
CN111965700A true CN111965700A (zh) 2020-11-20

Family

ID=73358069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910418078.8A Pending CN111965700A (zh) 2019-05-20 2019-05-20 光纤声波传感地震数据中零波数噪声的剔除方法及系统

Country Status (1)

Country Link
CN (1) CN111965700A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115542392A (zh) * 2022-10-11 2022-12-30 中国科学技术大学 基于分布式光纤声波传感的地下断层自动检测方法和装置
WO2023123951A1 (zh) * 2021-12-29 2023-07-06 中国石油天然气集团有限公司 利用局部fk变换提高das信噪比的方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158801A (zh) * 2015-07-29 2015-12-16 中国石油天然气集团公司 光缆耦合噪声的压制处理方法和装置
WO2016053582A1 (en) * 2014-10-01 2016-04-07 Halliburton Energy Services, Inc. Trace downsampling of distributed acoustic sensor data
US20170090054A1 (en) * 2014-07-17 2017-03-30 Halliburton Energy Services, Inc. Noise removal for distributed acoustic sensing data
US20170235006A1 (en) * 2014-09-12 2017-08-17 Halliburton Energy Services, Inc. Noise removal for distributed acoustic sensing data
CN107436451A (zh) * 2017-07-26 2017-12-05 西安交通大学 一种自动计算地震数据光缆耦合噪声强弱程度的振幅谱方法
WO2018044309A1 (en) * 2016-08-31 2018-03-08 Halliburton Energy Services, Inc. Real-time processing and control of das vsp surveys
WO2018136056A1 (en) * 2017-01-18 2018-07-26 Halliburton Energy Services, Inc. Gauge length effect and gauge length conversion

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170090054A1 (en) * 2014-07-17 2017-03-30 Halliburton Energy Services, Inc. Noise removal for distributed acoustic sensing data
US20170235006A1 (en) * 2014-09-12 2017-08-17 Halliburton Energy Services, Inc. Noise removal for distributed acoustic sensing data
WO2016053582A1 (en) * 2014-10-01 2016-04-07 Halliburton Energy Services, Inc. Trace downsampling of distributed acoustic sensor data
CN105158801A (zh) * 2015-07-29 2015-12-16 中国石油天然气集团公司 光缆耦合噪声的压制处理方法和装置
WO2018044309A1 (en) * 2016-08-31 2018-03-08 Halliburton Energy Services, Inc. Real-time processing and control of das vsp surveys
WO2018136056A1 (en) * 2017-01-18 2018-07-26 Halliburton Energy Services, Inc. Gauge length effect and gauge length conversion
CN107436451A (zh) * 2017-07-26 2017-12-05 西安交通大学 一种自动计算地震数据光缆耦合噪声强弱程度的振幅谱方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JIANYOU CHEN,ET AL.: "DISTRIBUTED ACOUSTIC SENSING COUPLING NOISE REMOVAL BASED ON SPARSE OPTIMIZATION", 《INTERPRETATION》 *
MENG LI,ET AL.: "Current and Future Applications of Distributed Acoustic Sensing as a New Reservoir Geophysics Tool", 《THE OPEN PETROLEUM ENGINEERING JOURNAL》 *
SUDHISH K. BAKKU,ET AL.: "Vertical Seismic Profiling Using Distributed Acoustic Sensing in a Hydrofrac TreatmentWell", 《GEOLOGY》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023123951A1 (zh) * 2021-12-29 2023-07-06 中国石油天然气集团有限公司 利用局部fk变换提高das信噪比的方法及装置
CN115542392A (zh) * 2022-10-11 2022-12-30 中国科学技术大学 基于分布式光纤声波传感的地下断层自动检测方法和装置
CN115542392B (zh) * 2022-10-11 2024-05-28 中国科学技术大学 基于分布式光纤声波传感的地下断层自动检测方法和装置

Similar Documents

Publication Publication Date Title
CN107505135B (zh) 一种滚动轴承复合故障提取方法及系统
Douglas Bandpass filtering to reduce noise on seismograms: Is there a better way?
CA2923888C (en) Dynamic threshold methods, systems, computer readable media, and program code for filtering noise and restoring attenuated high-frequency components of acoustic signals
CN107290589B (zh) 基于短时分数阶傅里叶变换的非线性信号时频分析方法
CN111965700A (zh) 光纤声波传感地震数据中零波数噪声的剔除方法及系统
CN105242309A (zh) 一种地震数据中规则干扰的压制方法及装置
CN111679328A (zh) 一种基于变分模态分解的瞬变电磁探测信号的降噪方法
CN105158801A (zh) 光缆耦合噪声的压制处理方法和装置
Zhou et al. An improved automatic picking method for arrival time of acoustic emission signals
US20230072038A1 (en) Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data
CN117111016B (zh) 复杂电磁环境下基于信道化的实时脉内分析方法及系统
CN111965704B (zh) 井中地震数据的光缆次生震荡噪声压制方法及系统
CN108363103B (zh) 规则干扰压制方法及装置
CN112464855A (zh) 基于eemd的随钻泥浆正脉冲信号处理方法及装置
CN112764109B (zh) 偶极横波反射波的分离提取方法及装置
CN105156901B (zh) 一种基于小波分析的光纤预警系统降噪方法及装置
CN109116425B (zh) 利用反射波的频谱设计滤波器去除噪音的方法
CN112782763A (zh) 一种地震品质因子估算方法、装置、设备及存储介质
CN113971415A (zh) 基于机器学习的地震数据多域处理方法及装置
CN106361328B (zh) 一种磁共振环境下的脑电信号提取方法
CN114624336B (zh) 基于轮式隔音的冲击回波声频法的信号检测装置及方法
CN112379427B (zh) 面波噪声压制方法及装置
CN113126164B (zh) 地震数据去噪方法及装置
CN112230286B (zh) 降低隧道无线地震仪噪声的方法及系统
CN112816056B (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