CN110673210A - 一种地震原始数据信噪比定量分析评价方法 - Google Patents

一种地震原始数据信噪比定量分析评价方法 Download PDF

Info

Publication number
CN110673210A
CN110673210A CN201910970434.7A CN201910970434A CN110673210A CN 110673210 A CN110673210 A CN 110673210A CN 201910970434 A CN201910970434 A CN 201910970434A CN 110673210 A CN110673210 A CN 110673210A
Authority
CN
China
Prior art keywords
shot
seismic
frequency band
noise ratio
data
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
Application number
CN201910970434.7A
Other languages
English (en)
Other versions
CN110673210B (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.)
Shenzhen Chengze Information Technology Co ltd
MCC Wukan Engineering Technology Co Ltd
Original Assignee
Northeast Petroleum University
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 Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN201910970434.7A priority Critical patent/CN110673210B/zh
Publication of CN110673210A publication Critical patent/CN110673210A/zh
Application granted granted Critical
Publication of CN110673210B publication Critical patent/CN110673210B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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

本发明涉及的是一种地震原始数据信噪比定量分析评价方法，它包括:一是求取炮域地震原始数据20‑40Hz频带范围内的信息能量；二是求取炮域地震原始数据0‑10Hz频带范围内的信息能量；三是求取炮域地震原始数据150Hz以上频带范围内的信息能量；四是定量求取低频信噪比与高频信噪比，并计算相应的数学期望与离散度；五是基于上述数据开展地震原始数据信噪比定量分析评价。本发明能够在空间平面上描述地震原始数据有效信号与噪声的分布及发育情况，可以定量的对地震原始数据信噪比开展分析与评价，为后续地震数据去噪处理提供参考，为地震原始数据品质评价提供了一种定量手段，本发明对油气、矿产资源勘探具有重要应用价值。

Description

一种地震原始数据信噪比定量分析评价方法
技术领域：
本发明涉及的是地震勘探中反射地震资料处理技术领域，具体涉及的是一种地震原始数据信噪比定量分析评价方法。
背景技术：
地震数据信噪比是地震数据处理过程中的一个重要评价指标，信噪比的高低直接影响地震资料品质的好坏。目前行业内有关地震数据信噪比的求取方法较多，但都是基于二维地震剖面数据发展而来的，无法从宏观上给出全工区地震资料有效信息及高低频噪声的分布及发育情况，目前基于三维地震原始数据的信噪比定量分析评价方法研究较少。
发明内容：
本发明的目的是提供一种地震原始数据信噪比定量分析评价方法，这种地震原始数据信噪比定量分析评价方法用于解决基于二维地震剖面数据求取的地震数据信噪比，无法从宏观上给出全工区地震资料有效信息及高低频噪声的分布及发育情况的问题。
本发明解决其技术问题所采用的技术方案是：这种地震原始数据信噪比定量分析评价方法:
步骤一、获取全工区炮域数据，求取计算参数；
步骤二、对炮域地震原始数据求取20-40Hz频带范围内的信息能量；
①对于炮点坐标为(x，y)的炮集记录，在能量计算时窗范围内求取20-40Hz频带范围内所有地震数据的均方根振幅值，并累加求和，最终和值记为E0，则炮点坐标为(x，y)处的炮集记录在20-40Hz频带范围内的信息能量是E0/Nxy；
②对工区内其它炮集重复步骤①，最终得到全工区炮集记录在20-40Hz频带范围内的信息能量平面分布图(x,y,E0/Nxy)；
步骤三、对炮域地震原始数据求取0-10Hz频带范围内的信息能量；
①对于炮点坐标为(x，y)的炮集记录，在能量计算时窗范围内求取0-10Hz频带范围内所有地震数据的均方根振幅值，并累加求和，最终和值记为E1，则炮点坐标为(x，y)处的炮集记录在0-10Hz频带范围内的信息能量是E1/Nxy；
②对工区内其它炮集重复步骤①，最终得到全工区炮集记录在0-10Hz频带范围内的信息能量平面分布图(x,y,E1/Nxy)；
步骤四、对炮域地震原始数据求取150Hz以上频带范围的信息能量；
①对于炮点坐标为(x，y)的炮集记录，在能量计算时窗范围内求取150Hz以上频带范围所有地震数据的均方根振幅值，并累加求和，最终和值记为E2，则炮点坐标为(x，y)处的炮集记录在150Hz以上频带范围的信息能量是E2/Nxy；
②对工区内其它炮集重复步骤①，最终得到全工区炮集记录在150Hz以上频带范围的信息能量平面分布图(x,y,E2/Nxy)；
步骤五、定量求取低频信噪比与高频信噪比，并计算相应的数学期望与离散度；
①低频信噪比snlow(x,y)＝(x,y,E0/Nxy)/(x,y,E1/Nxy)；
②高频信噪比snhigh(x,y,)＝(x,y,E0/Nxy)/(x,y,E2/Nxy)；
低频信噪比对应的数学期望与离散度：
①数学期望
Figure BDA0002231896220000021
②离散度
Figure BDA0002231896220000022
高频信噪比对应的数学期望与离散度：
①数学期望
②离散度
步骤六、基于上述数据开展地震原始资料信噪比定量分析评价。
上述方案中步骤一具体为：
①对于炮点坐标为(x，y)的炮集记录，记录其排列个数为Nxy；
②选定能量计算时窗，其中时窗起始时间应晚于炮集初至波出现时间，时窗终止时间设定为3s；
③统计全工区炮集个数为Nshot。
本发明具有以下有益效果：
1.本发明能够在空间平面上描述地震原始数据有效信号与噪声的分布及发育情况；
2.本发明可以定量的对地震原始数据信噪比开展分析与评价；
3.本发明为后续地震数据去噪处理提供参考，为地震原始数据品质评价提供了一种定量手段。
附图说明
图1地震原始数据20-40Hz频带范围信息能量平面图。
图2地震原始数据0-10Hz频带范围信息能量平面图。
图3地震原始数据150Hz以上频带范围信息能量平面图。
图4低频信噪比平面图。
图5高频信噪比平面图。
具体实施方式
下面对本发明做进一步的说明：
实施例1：
一种地震原始数据信噪比定量分析评价方法，针对冀东油田老爷庙区块某地震原始资料为例，具体为以下步骤：
(1)获取全工区炮域数据，炮集记录排列个数为8，全工区共有炮集32000个，选定能量计算时窗，其中时窗起始时间应晚于炮集初至波出现时间，时窗终止时间设定为3s；
(2)求取炮域地震原始数据20-40Hz频带范围内的信息能量，如图1所示；
(3)求取炮域地震原始数据0-10Hz频带范围内的信息能量，如图2所示；
(4)求取炮域地震原始数据求取150Hz以上频带范围的信息能量，如图3所示；
(5)定量求取低频信噪比与高频信噪比，并计算相应的数学期望与离散度，如图4，图5所示；
(6)根据高频信噪比、低频信噪比及其对应的数学期望及离散度的大小开展信噪比分析与评价。

Claims (2)

1.一种地震原始数据信噪比定量分析评价方法，其特征在于：
步骤一、获取全工区炮域数据，求取计算参数；
步骤二、对炮域地震原始数据求取20-40Hz频带范围内的信息能量；
①对于炮点坐标为(x，y)的炮集记录，在能量计算时窗范围内求取20-40Hz频带范围内所有地震数据的均方根振幅值，并累加求和，最终和值记为E0，则炮点坐标为(x，y)处的炮集记录在20-40Hz频带范围内的信息能量是E0/Nxy；
②对工区内其它炮集重复步骤①，最终得到全工区炮集记录在20-40Hz频带范围内的信息能量平面分布图(x,y,E0/Nxy)；
步骤三、对炮域地震原始数据求取0-10Hz频带范围内的信息能量；
①对于炮点坐标为(x，y)的炮集记录，在能量计算时窗范围内求取0-10Hz频带范围内所有地震数据的均方根振幅值，并累加求和，最终和值记为E1，则炮点坐标为(x，y)处的炮集记录在0-10Hz频带范围内的信息能量是E1/Nxy；
②对工区内其它炮集重复步骤①，最终得到全工区炮集记录在0-10Hz频带范围内的信息能量平面分布图(x,y,E1/Nxy)；
步骤四、对炮域地震原始数据求取150Hz以上频带范围的信息能量；
①对于炮点坐标为(x，y)的炮集记录，在能量计算时窗范围内求取150Hz以上频带范围所有地震数据的均方根振幅值，并累加求和，最终和值记为E2，则炮点坐标为(x，y)处的炮集记录在150Hz以上频带范围的信息能量是E2/Nxy；
②对工区内其它炮集重复步骤①，最终得到全工区炮集记录在150Hz以上频带范围的信息能量平面分布图(x,y,E2/Nxy)；
步骤五、定量求取低频信噪比与高频信噪比，并计算相应的数学期望与离散度；
①低频信噪比snlow(x,y)＝(x,y,E0/Nxy)/(x,y,E1/Nxy)；
②高频信噪比snhigh(x,y,)＝(x,y,E0/Nxy)/(x,y,E2/Nxy)；
低频信噪比对应的数学期望与离散度：
①数学期望
Figure FDA0002231896210000011
②离散度
Figure FDA0002231896210000021
高频信噪比对应的数学期望与离散度：
①数学期望
Figure FDA0002231896210000022
②离散度
Figure FDA0002231896210000023
步骤六、基于上述数据开展地震原始资料信噪比定量分析评价。
2.根据权利要求1所述的地震原始数据信噪比定量分析评价方法，其特征在于：所述的步骤一具体为：
①对于炮点坐标为(x，y)的炮集记录，记录其排列个数为Nxy；
②选定能量计算时窗，其中时窗起始时间应晚于炮集初至波出现时间，时窗终止时间设定为3s；
③统计全工区炮集个数为Nshot。
CN201910970434.7A 2019-10-13 2019-10-13 一种地震原始数据信噪比定量分析评价方法 Active CN110673210B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910970434.7A CN110673210B (zh) 2019-10-13 2019-10-13 一种地震原始数据信噪比定量分析评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910970434.7A CN110673210B (zh) 2019-10-13 2019-10-13 一种地震原始数据信噪比定量分析评价方法

Publications (2)

Publication Number Publication Date
CN110673210A true CN110673210A (zh) 2020-01-10
CN110673210B CN110673210B (zh) 2021-06-04

Family

ID=69082005

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910970434.7A Active CN110673210B (zh) 2019-10-13 2019-10-13 一种地震原始数据信噪比定量分析评价方法

Country Status (1)

Country Link
CN (1) CN110673210B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308551A (zh) * 2020-03-12 2020-06-19 王美琪 距离半径内炮数均值标准评价法
CN114114418A (zh) * 2020-08-28 2022-03-01 中国石油化工股份有限公司 地震信号增益方法、装置、计算机设备和存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004264A (zh) * 2010-10-18 2011-04-06 中国石油化工股份有限公司 一种地震采集资料质量定量分析与评价方法
RU2433425C2 (ru) * 2010-01-29 2011-11-10 Юрий Николаевич Жуков Способ сейсмической разведки при поиске углеводородов и способ определения залегания продуктивных на углеводороды пластов и сейсмическая станция для его осуществления
CN104216014A (zh) * 2014-09-25 2014-12-17 郭平 一种地震信号分频处理方法
CN104977602A (zh) * 2015-06-19 2015-10-14 中国石油天然气集团公司 一种地震数据采集施工的控制方法及装置
CN105785441A (zh) * 2016-03-07 2016-07-20 郑鸿明 地震数据的信噪比分析方法
CN106125134A (zh) * 2016-07-04 2016-11-16 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于双曲时窗的地震数据信噪比计算方法
WO2019030659A1 (en) * 2017-08-10 2019-02-14 Seismic Apparition Gmbh METHOD FOR ACQUIRING AND PROCESSING SEISMIC DATA

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2433425C2 (ru) * 2010-01-29 2011-11-10 Юрий Николаевич Жуков Способ сейсмической разведки при поиске углеводородов и способ определения залегания продуктивных на углеводороды пластов и сейсмическая станция для его осуществления
CN102004264A (zh) * 2010-10-18 2011-04-06 中国石油化工股份有限公司 一种地震采集资料质量定量分析与评价方法
CN104216014A (zh) * 2014-09-25 2014-12-17 郭平 一种地震信号分频处理方法
CN104977602A (zh) * 2015-06-19 2015-10-14 中国石油天然气集团公司 一种地震数据采集施工的控制方法及装置
CN105785441A (zh) * 2016-03-07 2016-07-20 郑鸿明 地震数据的信噪比分析方法
CN106125134A (zh) * 2016-07-04 2016-11-16 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于双曲时窗的地震数据信噪比计算方法
WO2019030659A1 (en) * 2017-08-10 2019-02-14 Seismic Apparition Gmbh METHOD FOR ACQUIRING AND PROCESSING SEISMIC DATA

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张军华 等: "地震资料信噪比定量计算及方法比较", 《石油地球物理勘探》 *
霍志坚 等: "地震记录信噪比的定量估算", 《江汉石油学院学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308551A (zh) * 2020-03-12 2020-06-19 王美琪 距离半径内炮数均值标准评价法
CN114114418A (zh) * 2020-08-28 2022-03-01 中国石油化工股份有限公司 地震信号增益方法、装置、计算机设备和存储介质

Also Published As

Publication number Publication date
CN110673210B (zh) 2021-06-04

Similar Documents

Publication Publication Date Title
CN107817527B (zh) 基于块稀疏压缩感知的沙漠地震勘探随机噪声压制方法
CN108845352B (zh) 基于vmd近似熵与多层感知机的沙漠地震信号去噪方法
CN109669212B (zh) 地震数据处理方法、地层品质因子估算方法与装置
CN110058305A (zh) 一种基于卷积神经网络的das地震数据降噪方法
CN104932017B (zh) 一种异常振幅压制方法及装置
CN110673210B (zh) 一种地震原始数据信噪比定量分析评价方法
CN106873031B (zh) 一种三维地震观测系统垂向分辨率定量分析评价方法
CN110632662A (zh) 一种利用DCNN-Inception网络自动识别微地震信号的算法
CN110823356A (zh) 基于梅尔频谱的分布式光纤入侵检测方法
CN107436451A (zh) 一种自动计算地震数据光缆耦合噪声强弱程度的振幅谱方法
CN113221746A (zh) 基于改进小波阈值函数的微震信号去噪方法
Li et al. Distributed acoustic sensing vertical seismic profile data denoising based on multistage denoising network
CN112558159B (zh) 一种声波测井初至拾取方法
CN103592684A (zh) 一种保持空间属性信息的海量地震数据压缩方法及装置
CN108562936A (zh) 一种裂缝预测方法、系统、存储介质及终端
CN113960660B (zh) 一种基于正演模拟的动校正畸变区域自动识别及切除方法
CN110007342B (zh) 一种用于低信噪比地震信号的时频域直接拾取初至方法及系统
CN115932959A (zh) 一种基于UNet的面波频散质量快速评估方法
CN106054245B (zh) 地震资料高分辨率处理定量评价方法
CN115456040A (zh) 一种基于卷积神经网络的p波捡拾算法
CN113219526A (zh) 一种地震波初至自动拾取与筛选方法
CN114563771A (zh) 基于聚类分析的双阈值激光雷达云层检测算法
CN112526611A (zh) 表层地震波品质因子的提取方法及装置
CN114063167B (zh) 一种用于压制外源干扰的地震资料处理方法
CN113721297B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221219

Address after: 430080 No. 17 metallurgical Avenue, Qingshan District, Hubei, Wuhan

Patentee after: MCC Wukan Engineering Technology Co.,Ltd.

Address before: Room 2011-2015, Building C1, 1983 Creative Town, No. 15, Nanxin Road, Nanling Village Community, Nanwan Street, Longgang District, Shenzhen, Guangdong 518100

Patentee before: Shenzhen Chengze Information Technology Co.,Ltd.

Effective date of registration: 20221219

Address after: Room 2011-2015, Building C1, 1983 Creative Town, No. 15, Nanxin Road, Nanling Village Community, Nanwan Street, Longgang District, Shenzhen, Guangdong 518100

Patentee after: Shenzhen Chengze Information Technology Co.,Ltd.

Address before: 163319 No. 99 Xuefu Street, Daqing Hi-tech Development Zone, Heilongjiang Province

Patentee before: NORTHEAST PETROLEUM University