CN112925024B - 一种地震记录的方波压制方法 - Google Patents

一种地震记录的方波压制方法 Download PDF

Info

Publication number
CN112925024B
CN112925024B CN202110102152.2A CN202110102152A CN112925024B CN 112925024 B CN112925024 B CN 112925024B CN 202110102152 A CN202110102152 A CN 202110102152A CN 112925024 B CN112925024 B CN 112925024B
Authority
CN
China
Prior art keywords
square wave
sample
value
square
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.)
Active
Application number
CN202110102152.2A
Other languages
English (en)
Other versions
CN112925024A (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.)
Geophysical Survey Research Institute jianghan Oilfield Branch china Petrochemical Corp
China Petroleum and Chemical Corp
Original Assignee
Geophysical Survey Research Institute jianghan Oilfield Branch china Petrochemical Corp
China Petroleum and Chemical Corp
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 Geophysical Survey Research Institute jianghan Oilfield Branch china Petrochemical Corp, China Petroleum and Chemical Corp filed Critical Geophysical Survey Research Institute jianghan Oilfield Branch china Petrochemical Corp
Priority to CN202110102152.2A priority Critical patent/CN112925024B/zh
Publication of CN112925024A publication Critical patent/CN112925024A/zh
Application granted granted Critical
Publication of CN112925024B publication Critical patent/CN112925024B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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
    • 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

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)

Abstract

本发明公开了一种地震记录的方波压制方法，属于石油勘探开发领域，该方法包括以下步骤：在单炮地震数据体上逐道、逐样点判断方波是否存在；如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值保持不变，接着处理下一个样点。本发明利用方波能量强，且持续时间长，取值恒定的特点，采取差分结合能量的方式进行判断，可以较为准确地判断方波是否存在；在预测方波时，通过均值结合中值技术可以快速的预测方波值，由于是逐样点处理，去噪比较彻底；该技术克服了常规方法的边界效应及对有效信号的破坏，保真度较高。

Description

一种地震记录的方波压制方法
技术领域
本发明属于石油勘探开发领域，更具体地，涉及一种地震记录中方波的压制方法。
背景技术
在地震记录中，方波通常是由于仪器或检波器漏码故障引起，在地震记录上表现为一个非常大的常数干扰值叠加在有效信号上，形同废道废炮。不同道的干扰值相差较大，可能只是影响局部数据，也可能对整道数据都有影响。
常规处理方法往往有去道去炮，内切除，高通滤波，中值滤波等等。去道去炮和内切除都造成了有效信号的丧失和经济上的损失。而高通滤波因为大值的影响会损害有效信号，且会造成严重的边界效应。中值滤波往往耗时较大，当方波连续多道出现时，去噪结果往往存在较大误差，且在没有方波的地方对有效信号有较大影响。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提出了一种地震记录的方波压制方法，用于解决去噪对有效信号损害较大的问题。
为实现上述目的，按照本发明的一个方面，提供了一种地震记录的方波压制方法，包括：
(1)在单炮地震数据体上逐道、逐样点进行处理，以各样点为中心，通过预设长度时窗内的数据进行方波存在与否判断；
(2)如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值保持不变。
在一些可选的实施方案中，所述通过预设长度时窗内的数据进行方波存在与否判断，包括：
通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断。
在一些可选的实施方案中，所述通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断，包括：
设地震记录为si，其中，下标i为样点序号，设置滑动时窗长度参数为2*L+1，识别方波的门槛因子为C，令： L为自然数；
对于第i个样点，在Ci＞C时，则第i个样点存在方波，否则第i个样点不存在方波。
在一些可选的实施方案中，所述进行方波预测与压制，包括：
通过均值结合中值技术进行方波预测与压制。
在一些可选的实施方案中，所述通过均值结合中值技术进行方波预测与压制，包括：
取步骤(1)中的同一个时窗内样点数据，令
将sj，(j＝i-L，...，i+L)分为两个序列Ai，(i＝1，...，M)和Bi，(i＝1，...，N)，M+N＝2*L+1，其中，Ai为序列sj中大于的样点值，Bi为序列sj中小于/>的样点值，M，N为自然数；
若M＞N，令IM＝M-L，则在序列Ai中找到第IM小的值为样点i处方波值，记为Vi；若M≤N，令IN＝N-L，在序列Bi中找到第IN大的值为样点i处方波值，记为Vi；
由s′i＝si-Vi确定去除方波后的结果s′i。
按照本发明的另一方面，提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述方法的步骤。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
利用方波能量强，且持续时间长，取值恒定的特点，采取差分结合能量的方式进行判断，可以较为准确地判断方波是否存在；在预测方波时，通过均值结合中值技术可以快速的预测方波值，由于是逐样点处理，去噪比较彻底；本发明克服了常规方法的边界效应及对有效信号的破坏，保真度较高。
附图说明
图1是本发明实施例提供的一种地震记录的方波压制方法的流程示意图；
图2是本发明实施例提供的一种压制方波前的地震数据图；
图3是本发明实施例提供的一种压制方波后的地震数据图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示，本发明实施例提供的一种地震记录方波压制方法，主要包括以下步骤：
S1：在单炮地震数据体上逐道、逐样点进行处理，以该样点为中心，通过预设长度时窗内的数据进行方波存在与否判断；
在本发明实施例中，步骤S1可以通过以下方式实现：
通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断。
具体地，设地震记录为si，其中，i为样点序号；
设置滑动时窗长度参数为2*L+1，识别方波的门槛因子为C，令：
当Ci＞C时，存在方波，否则不存在方波，L为自然数。
S2：如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值不变，接着处理下一个样点，以此类推。
在本发明实施例中，步骤S2可以通过以下方式实现：
通过均值结合中值技术进行方波预测与压制。
具体地，取步骤S1中的同一个时窗内样点数据，令：
将sj，(j＝i-L，...，i+L)分为两个序列Ai，(i＝1，...，M)和Bi，(i＝1，...，N)，M+N＝2*L+1。其中Ai为序列sj中大于的样点值，Bi为序列sj中小于/>的样点值，M，N为自然数；
若M＞N，令IM＝M-L，则在序列Ai中找到第IM小的值为样点i处方波值，记为Vi，若M≤N，令IN＝N-L，在序列Bi中找到第IN大的值为样点i处方波值，记为Vi。令：
s′i＝si-Vi (5)
s′i即为去除方波后的结果。
图2是压制方波前的地震数据图。由于仪器漏码引起的方波能量非常强，几乎将有效信号完全淹没，如果不做去噪处理将严重影响后续子波处理、速度分析以及叠加效果。
图3是压制方波后的地震数据图。与图2相比，方波被完全去除，掩盖在方波下的地震信息也显现出来。去噪效果非常明显。
综上，本发明的一种地震记录方波压制方法主要改进包括三点：(1)利用方波能量强，且持续时间长，取值恒定的特点，采取差分结合能量的方式进行判断，可以较为准确地判断方波是否存在；(2)在预测方波时，通过均值结合中值技术可以快速的预测方波值，由于是逐样点处理，去噪比较彻底；(3)该技术克服了常规方法的边界效应及对有效信号的破坏，保真度较高。
本申请还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，其上存储有计算机程序，程序被处理器执行时实现方法实施例中的地震记录的方波压制方法。
需要指出，根据实施的需要，可将本申请中描述的各个步骤/部件拆分为更多步骤/部件，也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件，以实现本发明的目的。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (4)

1.一种地震记录的方波压制方法，其特征在于，包括：
(1)在单炮地震数据体上逐道、逐样点进行处理，以各样点为中心，通过预设长度时窗内的数据进行方波存在与否判断；
(2)如果该样点存在方波则进行方波预测与压制，如果不存在方波，则该样点值保持不变；
所述进行方波预测与压制，包括：
通过均值结合中值技术进行方波预测与压制；
所述通过均值结合中值技术进行方波预测与压制，包括：
取步骤(1)中的同一个时窗内样点数据，令
将sj,j＝i-L,…,i+L分为两个序列Ai,i＝1,…,M和Bi,i＝1,…,N,M+N＝2*L+1，其中，Ai为序列sj中大于的样点值，Bi为序列sj中小于/>的样点值，M，N，L为自然数，地震记录为s，下标为样点序号；
若M>N，令IM＝M-L，则在序列Ai中找到第IM小的值为样点i处方波值，记为Vi；若M≤N，令IN＝N-L，在序列Bi中找到第IN大的值为样点i处方波值，记为Vi；
由s′i＝si-Vi确定去除方波后的结果s′i。
2.根据权利要求1所述的方法，其特征在于，所述通过预设长度时窗内的数据进行方波存在与否判断，包括：
通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断。
3.根据权利要求2所述的方法，其特征在于，所述通过预设长度时窗内的数据，采用差分结合能量的方式进行方波存在与否判断，包括：
设地震记录为si，其中，下标i为样点序号，设置滑动时窗长度参数为2*L+1，识别方波的门槛因子为C，令： L为自然数；
对于第i个样点，在Ci>C时，则第i个样点存在方波，否则第i个样点不存在方波。
4.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至3任一项所述方法的步骤。
CN202110102152.2A 2021-01-26 2021-01-26 一种地震记录的方波压制方法 Active CN112925024B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110102152.2A CN112925024B (zh) 2021-01-26 2021-01-26 一种地震记录的方波压制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110102152.2A CN112925024B (zh) 2021-01-26 2021-01-26 一种地震记录的方波压制方法

Publications (2)

Publication Number Publication Date
CN112925024A CN112925024A (zh) 2021-06-08
CN112925024B true CN112925024B (zh) 2023-10-20

Family

ID=76166110

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110102152.2A Active CN112925024B (zh) 2021-01-26 2021-01-26 一种地震记录的方波压制方法

Country Status (1)

Country Link
CN (1) CN112925024B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133248A (zh) * 2014-08-08 2014-11-05 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种高保真声波干扰压制方法
CN105717543A (zh) * 2016-01-25 2016-06-29 浪潮(北京)电子信息产业有限公司 一种噪声压制方法及系统
CN105911590A (zh) * 2016-05-09 2016-08-31 刘宏杰 一种谐波噪声的压制处理方法及装置
CN107966734A (zh) * 2017-09-22 2018-04-27 中国地质大学（北京） 多分量地震数据的矢量去噪方法
CN109425897A (zh) * 2017-08-30 2019-03-05 中国石油化工股份有限公司 消除地震资料野值干扰的方法及系统
CN109425894A (zh) * 2017-08-23 2019-03-05 中国石油化工股份有限公司 一种地震异常道检测方法及装置
CN109901224A (zh) * 2019-02-15 2019-06-18 中交铁道设计研究总院有限公司 一种地震资料低频信号保护压制噪声方法
CN112014884A (zh) * 2019-05-30 2020-12-01 中国石油天然气集团有限公司 压制近炮点强能量噪声的方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9612351B2 (en) * 2012-07-10 2017-04-04 Chevron U.S.A. Inc. System and method for estimating and attenuating noise in seismic data

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133248A (zh) * 2014-08-08 2014-11-05 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种高保真声波干扰压制方法
CN105717543A (zh) * 2016-01-25 2016-06-29 浪潮(北京)电子信息产业有限公司 一种噪声压制方法及系统
CN105911590A (zh) * 2016-05-09 2016-08-31 刘宏杰 一种谐波噪声的压制处理方法及装置
CN109425894A (zh) * 2017-08-23 2019-03-05 中国石油化工股份有限公司 一种地震异常道检测方法及装置
CN109425897A (zh) * 2017-08-30 2019-03-05 中国石油化工股份有限公司 消除地震资料野值干扰的方法及系统
CN107966734A (zh) * 2017-09-22 2018-04-27 中国地质大学（北京） 多分量地震数据的矢量去噪方法
CN109901224A (zh) * 2019-02-15 2019-06-18 中交铁道设计研究总院有限公司 一种地震资料低频信号保护压制噪声方法
CN112014884A (zh) * 2019-05-30 2020-12-01 中国石油天然气集团有限公司 压制近炮点强能量噪声的方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
地震数据高维统计滤波方法;王福 等;石油物探;第58卷(第03期);第335-345页 *
基于数学形态学的大地电磁强干扰分离及应用;李晋;中国博士论文全文数据库;全文 *
夏洪瑞 等.环形非对称中值滤波技术在可控震源资料去噪处理中的应用.石油物探.2009,第48卷(第05期),第510-514、520页. *
循环中值滤波在消除可控震源地震资料噪声中的应用;张怿平 等;江汉石油职工大学学报;第22卷(第02期);第93-96页 *
环形非对称中值滤波技术在可控震源资料去噪处理中的应用;夏洪瑞 等;石油物探;第48卷(第05期);第510-520页 *

Also Published As

Publication number Publication date
CN112925024A (zh) 2021-06-08

Similar Documents

Publication Publication Date Title
Fontani et al. Hiding traces of median filtering in digital images
CN113237432B (zh) 激光焊接熔深提取方法及设备
CN113538263A (zh) 基于改进DeblurGAN模型的去运动模糊方法、介质及设备
CN115631494A (zh) 财务数据处理方法、系统、电子设备及存储介质
Wang et al. High quality impulse noise removal via non‐uniform sampling and autoregressive modelling based super‐resolution
CN112925024B (zh) 一种地震记录的方波压制方法
CN110852443A (zh) 特征稳定性检测方法、设备及计算机可读介质
Zou et al. Anti‐Forensics of Image Contrast Enhancement Based on Generative Adversarial Network
CN110929574A (zh) 一种红外弱小目标快速检测方法
CN108562936B (zh) 一种裂缝预测方法、系统、存储介质及终端
Tibi et al. Comparative Study of the Performance of Seismic Waveform Denoising Methods Using Local and Near‐Regional Data
Yang et al. Chaotic signal denoising based on energy selection TQWT and adaptive SVD
CN116756486A (zh) 基于声光电磁多源数据融合的海上目标识别方法及装置
CN112462425B (zh) 海底节点混采数据中交叉干扰源的识别方法及装置
CN109031431A (zh) 一种针对探地雷达数据的数据处理方法及系统
CN114897723A (zh) 一种基于生成式对抗网络的图像生成与加噪方法
US20150254812A1 (en) Automated Selection of Filter Parameters for Seismic Analysis
CN109656916A (zh) 一种雷达目标数据的获取方法
CN111157966A (zh) 一种假目标干扰识别方法
CN116381790A (zh) 一种三维叠前地震记录相干噪声衰减方法
CN112560457B (zh) 基于非监督的文本去噪方法、系统、电子设备及存储介质
CN117197014B (zh) 一种可降噪的肺部医学图像融合方法、系统和电子设备
CN114200523B (zh) 一种补偿可控震源单炮频率缺失的方法、装置和存储介质
LU504136B1 (en) An intelligent spectral processing system for libs
CN111240583B (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