CN102176056B - 基于时频分析的地震有效波多域能量补偿方法 - Google Patents

基于时频分析的地震有效波多域能量补偿方法 Download PDF

Info

Publication number
CN102176056B
CN102176056B CN201110040394XA CN201110040394A CN102176056B CN 102176056 B CN102176056 B CN 102176056B CN 201110040394X A CN201110040394X A CN 201110040394XA CN 201110040394 A CN201110040394 A CN 201110040394A CN 102176056 B CN102176056 B CN 102176056B
Authority
CN
China
Prior art keywords
frequency
geological data
energy
seismic
time
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.)
Expired - Fee Related
Application number
CN201110040394XA
Other languages
English (en)
Other versions
CN102176056A (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 Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
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 Petroleum and Chemical Corp, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201110040394XA priority Critical patent/CN102176056B/zh
Publication of CN102176056A publication Critical patent/CN102176056A/zh
Application granted granted Critical
Publication of CN102176056B publication Critical patent/CN102176056B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及地震资料处理领域,是一种基于时频分析的地震有效波多域能量补偿方法。在地震波传播过程中,引起能量衰减的因素很多,衰减关系极为复杂,并且有效波和干扰波混杂在一起,利用单一的补偿方法无法取得满意的效果。本发明从地震资料自身包含的信息出发,依据频率扫描和分频处理,分级别地去除面波;然后在不同频带内进行振幅衰减关系的拟合求取并进行有效补偿;在空间进行子波一致性处理,实现空间波形和能量的均衡。一方面,本发明是集频率、空间、时间的多域综合处理,实现了多方位多级别的补偿;另一方面,该补偿的参数是来自对地震资料进行分频扫描本身,而非来自经验关系或其他数据,保证了参数选取的正确性。

Description

基于时频分析的地震有效波多域能量补偿方法
技术领域
本发明属于地震资料处理方法领域,具体是一种基于时频分析的地震有效波多域能量补偿方法。
现有技术
地震波在传播过程中,由于受到波前几何扩散和地层非完全弹性等因素的影响,来自深层和远道的反射波能量会显著降低。现有的能量补偿方法是基于经验关系和统计数据的,即使用                                                
Figure 201110040394X100002DEST_PATH_IMAGE001
的补偿和校正方法,其中参数的选取由处理员根据经验选取。另一项技术及反Q滤波技术也存在反Q因子提取不准的问题。同时,在补偿有效波的同时也将干扰波的能量提升上来。随着勘探的深入,地震数据处理的精度要求越来越高,这些方法已经不能满足要求。
发明内容
本发明的目的为了解决现有的能量补偿方法存在的精度不高、参数选取困难、不能提高信噪比等问题,提出一种从地震资料自身出发,通过时频分析,确定面波和有效波的分布范围,进行分频噪音压制,更为精确的基于时频分析的地震有效波多域能量补偿方法。
针对一个地震数据,时频空间域能量补偿是一个连续的处理过程,该发明具体的处理步骤分为以下几步:
1)输入原始SEGY地震数据。
2)通过球面扩散补偿方法,补偿原始SEGY地震数据由于地层吸收引起的空间能量衰减。
3)在原始SEGY地震数据中选取一段信噪比较好,波形较好的记录,求取子波校正因子,进行地震子波校正。
4)通过付利叶变换将时间域地震数据变换为频率域地震数据。
5)对频率域地震数据分开两步进行处理。一方面对频率域地震数据进行时间频率的对应关系分析,得到不同频率能量的分布谱;另一方面对频率域地震数据进行频谱扫描和频谱拟合,得到分段拓频因子。两个方面结合得到频率补偿后的频率域数据。
6)通过付利叶反变换将频率域地震数据变换为时间域地震数据。
上述方案中利用分频率扫描进行数据分析并去除面波,就是对通过采集的地震数据进行时频分析,在时间频率域内分析地震记录中面波成分在时间—频率图上的分布范围,然后通过相应的曲线拟合截取,达到去除面波的目的;在振幅补偿过程中,通过对地震炮集记录能量的拟合,较为准确地拟合出地震波在传播过程中的能量变化规律,从而消除炮集中由于检波器安置,随机噪音干扰等因素引起的地震波能量差异,对地震记录达到更准确的补偿; 在时频空间域能量补偿中,通过时频分析剖面对地震记录进行具体直观的数据分频,补偿,通过中心频率分布图和中心能量分布图,可以看出该方法具有较为理想的效果。
本发明的技术效果主要表现在:
在地震波传播过程中,引起能量衰减的因素很多,衰减关系极为复杂,并且有效波和干扰波混杂在一起,利用单一的补偿方法无法取得满意的效果。本发明从地震资料自身包含的信息出发,依据频率扫描和分频处理,分级别地去除面波;然后在不同频带内进行振幅衰减关系的拟合求取并进行有效补偿;在空间进行子波一致性处理,实现空间波形和能量的均衡。一方面,本发明是集频率、空间、时间的多域综合处理,实现了多方位多级别的补偿;另一方面,该补偿的参数是来自对地震资料进行分频扫描本身,而非来自经验关系或其他数据,保证了参数选取的正确性。
附图说明
图1为本发明的工作原理流程图。
图2为胜利油田某地区的原始实际地震炮集资料。
图3为利用本发明方法对其进行处理以后的结果。
图4为原始炮集资料的统计自相关分析结果。
图5为利用本发明方法处理结果的统计自相关分析结果。
图6为利用本发明方法处理前后的频谱对比。
具体实施方式
针对一个地震数据,该发明具体的处理步骤分为以下几步:
1)输入原始SEGY地震数据。
2)通过球面扩散补偿方法,补偿原始SEGY地震数据由于地层吸收引起的空间能量衰减。
3)在原始SEGY地震数据中选取一段信噪比较好,波形较好的记录,求取子波校正因子,进行地震子波校正。
4)通过付利叶变换将时间域地震数据变换为频率域地震数据。
5)对频率域地震数据分开两步进行处理。一方面对频率域地震数据进行时间频率的对应关系分析,得到不同频率能量的分布谱;另一方面对频率域地震数据进行频谱扫描和频谱拟合,得到分段拓频因子。两个方面结合得到频率补偿后的频率域数据。
6)通过付利叶反变换将频率域地震数据变换为时间域地震数据。
通过以上具体的步骤处理,达到对地震记录的时频能量补偿并对其进行具体的分析对比。图1为胜利油田某地区的原始实际地震炮集资料。可见面波发育,深层能量较弱。图2为利用基于时频分析的地震有效波多域能量补偿方法对其进行处理以后的结果。可见处理后的记录上面波不见了,信噪比有了提高。同时,深层能量有了增强。图3为原始炮集资料的统计自相关分析结果,图4为处理结果的统计自相关分析结果,对比可见,处理后子波的一致性得到了提高。图5为处理前后的频谱对比,可见经过处理后的地震资料频带明显得到了展宽,资料质量得到了提高。

Claims (1)

1.一种基于时频分析的地震有效波多域能量补偿方法,其特征在于:
1)输入原始SEGY地震数据;
2)通过球面扩散补偿方法,补偿原始SEGY地震数据由于地层吸收引起的空间能量衰减;
3)在原始SEGY地震数据中选取一段信噪比较好,波形较好的记录,求取子波校正因子,进行地震子波校正;
4)通过付利叶变换将时间域地震数据变换为频率域地震数据;
5)对频率域地震数据分开两步进行处理,一方面对频率域地震数据进行时间频率的对应关系分析,得到不同频率能量的分布谱;另一方面对频率域地震数据进行频谱扫描和频谱拟合,得到分段拓频因子;两个方面结合得到频率补偿后的频率域地震数据;
6)通过付利叶反变换将频率域地震数据变换为时间域地震数据。
CN201110040394XA 2011-02-18 2011-02-18 基于时频分析的地震有效波多域能量补偿方法 Expired - Fee Related CN102176056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110040394XA CN102176056B (zh) 2011-02-18 2011-02-18 基于时频分析的地震有效波多域能量补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110040394XA CN102176056B (zh) 2011-02-18 2011-02-18 基于时频分析的地震有效波多域能量补偿方法

Publications (2)

Publication Number Publication Date
CN102176056A CN102176056A (zh) 2011-09-07
CN102176056B true CN102176056B (zh) 2012-11-28

Family

ID=44519258

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110040394XA Expired - Fee Related CN102176056B (zh) 2011-02-18 2011-02-18 基于时频分析的地震有效波多域能量补偿方法

Country Status (1)

Country Link
CN (1) CN102176056B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105954799A (zh) * 2016-04-27 2016-09-21 中国石油天然气股份有限公司 一种基于加权叠加的时频域地震资料处理方法

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998700B (zh) * 2011-09-08 2015-09-23 中国石油天然气集团公司 一种垂直地震剖面资料的地表一致性能量补偿方法
CN103630932B (zh) * 2013-09-23 2016-03-30 中国海洋石油总公司 一种地震数据分形保幅方法
CN105093282B (zh) * 2014-05-16 2018-01-05 中国石油化工股份有限公司 基于频率约束的能量置换面波压制方法
CN107219551B (zh) * 2016-03-22 2019-08-27 中国石油化工股份有限公司 拓宽地震数据频带的方法及装置
CN106772574B (zh) * 2016-11-28 2019-02-15 中国石油天然气集团有限公司 一种保持低频信息的谱延拓方法
CN109238454B (zh) * 2018-11-14 2020-12-08 国网河北省电力有限公司电力科学研究院 电力变压器绕组变形累积效应的检测方法及装置
CN111766631B (zh) * 2019-04-01 2022-12-02 中国石油天然气股份有限公司 一种面波压制方法及装置
CN110261899B (zh) * 2019-04-26 2021-03-23 中国石油化工股份有限公司 地震数据z字形干扰波去除方法
CN112034510B (zh) * 2019-06-04 2023-04-07 中国石油天然气股份有限公司 低频补偿方法及低频补偿装置
CN110737018B (zh) * 2019-07-09 2021-09-07 中国石油化工股份有限公司 Vsp地震资料各向异性建模方法
CN112698396B (zh) * 2019-10-23 2024-04-30 中国石油天然气集团有限公司 一种压制极强环境噪音的数据处理方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000016127A1 (en) * 1998-09-11 2000-03-23 Pgs Data Processing Inc. An improved method of dual wavefield reinforcement
CN101609161A (zh) * 2009-07-17 2009-12-23 中国石化集团胜利石油管理局 基于地震层序体理论多尺度资料联合频带拓展方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000016127A1 (en) * 1998-09-11 2000-03-23 Pgs Data Processing Inc. An improved method of dual wavefield reinforcement
CN101609161A (zh) * 2009-07-17 2009-12-23 中国石化集团胜利石油管理局 基于地震层序体理论多尺度资料联合频带拓展方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
白桦等.基于时频分析的地层吸收补偿.《石油地球物理勘探》.1999,第34卷(第6期),第642-646页. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105954799A (zh) * 2016-04-27 2016-09-21 中国石油天然气股份有限公司 一种基于加权叠加的时频域地震资料处理方法
CN105954799B (zh) * 2016-04-27 2018-04-03 中国石油天然气股份有限公司 一种基于加权叠加的时频域地震资料处理方法

Also Published As

Publication number Publication date
CN102176056A (zh) 2011-09-07

Similar Documents

Publication Publication Date Title
CN102176056B (zh) 基于时频分析的地震有效波多域能量补偿方法
CN103376464B (zh) 一种地层品质因子反演方法
CN101334483B (zh) 一种在地震数据处理中衰减瑞雷波散射噪声的方法
CN102109612B (zh) 一种地震波吸收衰减补偿方法
CN109669212B (zh) 地震数据处理方法、地层品质因子估算方法与装置
CN102323617B (zh) 一种复杂地表的二维地震资料连片处理方法
CN104950326B (zh) 基于目的层频谱的可控震源非线性扫描信号的设计方法
CN102798894B (zh) 一种基于稀疏反演的滑动扫描谐波压制方法
CN109738951B (zh) 一种基于地震同相轴子波谱的时变反褶积方法
CN102262243B (zh) 一种滤波法可控震源地震数据谐波干扰压制方法
CN102681014A (zh) 基于多项式拟合的规则线性干扰压制方法
CN104614769B (zh) 一种压制地震面波的聚束滤波方法
CN101852866B (zh) 一种叠后地震数据滤波方法
CN105445801A (zh) 一种消除二维地震资料随机噪音的处理方法
CN103913770A (zh) 基于vsp资料对地震数据进行处理的方法
CN103675901A (zh) 一种时频域可控震源近地表吸收补偿方法
CN105510975B (zh) 提高地震数据信噪比的方法及装置
CN103576197A (zh) 一种转换波角道集抽取方法
CN106950600B (zh) 一种近地表散射面波的去除方法
CN102053260B (zh) 获得地震纵波的方位速度的方法及处理地震数据的方法
CN107238864A (zh) 基于优势频带约束的地震波能量时频补偿方法
CN106371140B (zh) 一种提高中深层地震资料分辨率的方法
Bing et al. Extracting near-borehole P and S reflections from array sonic logging data
CN104502977A (zh) 一种井控保幅高分辨率地震资料处理方法
CN106054245B (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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121128

Termination date: 20200218

CF01 Termination of patent right due to non-payment of annual fee