CN102841380B - 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 - Google Patents
一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 Download PDFInfo
- Publication number
- CN102841380B CN102841380B CN201110170561.2A CN201110170561A CN102841380B CN 102841380 B CN102841380 B CN 102841380B CN 201110170561 A CN201110170561 A CN 201110170561A CN 102841380 B CN102841380 B CN 102841380B
- Authority
- CN
- China
- Prior art keywords
- noise
- distributed
- seismologic record
- omega
- denoising
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法,属于石油地球物理勘探和地质工程调查领域。所述方法首先将分布有强能量相干噪音的地震记录与分布有弱能量相干噪音的地震记录分开,然后对分布有强能量相干噪音的地震记录加强去噪,对分布有弱能量相干噪音的地震记录降低去噪程度甚至不去噪,使陡倾角有效反射波得到保护。本发明解决了山地复杂构造地震资料叠前线性相干噪音压制对有效反射波造成损失的问题,实现了超大数据量地震数据自动检测分类及叠前空变去噪,有效保护了低频有效反射信息,为复杂构造的成像提供了基础较好的资料。
Description
技术领域
本发明属于石油地球物理勘探和地质工程调查领域,具体涉及一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法,是专门针对当前起伏地表超大数据量地震资料炮域相干噪音的。
背景技术
当前在高程变化剧烈的起伏地表及复杂构造区采集的地震资料越来越多,地表的剧烈变化产生强能量的相干噪音,相干噪音在地震记录上表现为陡倾角反射,地下构造复杂时有效波同相轴倾角大,而地表起伏变化剧烈时,相干噪音分布在空间上极不均匀,去除相干噪音容易损失有效波。
目前三维地震资料相干噪音去除技术主要有fk滤波、fx域相干噪音衰减、速度变换滤波法等多种技术,而目前针对三维地震资料叠前相干噪音压制的方法还没有解决空变问题。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法,解决山地复杂构造地震资料叠前线性相干噪音压制对有效反射波造成损失的问题,保护低频有效反射信息。
本发明是通过以下技术方案实现的:
一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法,所述方法首先将分布有强能量相干噪音的地震记录与分布有弱能量相干噪音的地震记录分开,然后对分布有强能量相干噪音的地震记录加强去噪,对分布有弱能量相干噪音的地震记录降低去噪程度甚至不去噪,使陡倾角有效反射波得到保护。
(1)地震记录分类
分析地震资料中的地震信号主频的分布规律,设定合适的门槛值,将地震信号主频低于门槛值的定义为分布有强相干噪音的地震记录,将地震信号主频高于门槛值的定义为分布有弱相干噪音的地震记录,从而将两者分开;
(2)相干噪音特征分析
对分布有强相干噪音的地震记录进行分析,分析强相干噪音的速度、频率分布范围以及与有效反射波的特征差异,进而确定强相干噪音的运动学特征;
(3)相干噪音压制
根据步骤(2)获得的强相干噪音的特征,对步骤(1)识别出来的分布有强相干噪音的地震记录采用FX域或均值加权噪音分离技术,形成噪音的模型,然后从原始地震记录中减掉噪音,得到去噪后的地震记录。
所述步骤(1)中获得地震信号主频的方法是将地震信号经过傅立叶变换到频率-振幅域,再对其求导,令其等于0,即振幅的最大值为此地震信号的主频,具体步骤如下:
假设地震信号为x(t),为连续时间信号,对其做傅立叶变换得到
其中,w为角频率w=2πf,f为频率;
X(ω)通常是复数,利用复数性质可得:
X(ω)=A(ω)eiφ(ω)
其中,A(ω)和φ(ω)分别为振幅谱和相位谱,它们是通过用下面的公式计算获得的:
其中,Xr(ω)和Xi(ω)为X(ω)的付里叶变换的实部和虚部,分别为:
Xr(ω)=A(ω)cosφ(ω),Xi(ω)=A(ω)sinφ(ω)
对X(ω)=A(ω)eiφ(ω)求导,令解此方程则得到的就是该信号主频。
与现有技术相比,本发明的有益效果是:本发明解决了山地复杂构造地震资料叠前线性相干噪音压制对有效反射波造成损失的问题,由于相干噪音压制同时可能使有效的低频有效信号能量损失,本发明将具有强线性相干噪音的地震记录加强去噪,质量好的记录不做去噪处理,这样实现了超大数据量地震数据自动检测分类及叠前空变去噪,有效保护了低频有效反射信息,为复杂构造的成像提供了基础较好的资料。
附图说明
图1是不同主频的地震记录图。
图2是本发明方法的步骤框图。
图3是本发明实施例中的去除强相干噪音的效果图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法,所述方法首先将分布有强能量相干噪音的地震记录与分布有弱能量相干噪音的地震记录分开,然后对分布有强能量相干噪音的地震记录加强去噪,对分布有弱能量相干噪音的地震记录降低去噪程度甚至不去噪,使陡倾角有效反射波得到保护。
如图2所示,所述方法包括以下步骤:
(1)地震记录分类
地震波传播时地表剧烈变化容易产生强能量低频相干噪音,对地震记录进行主频分析发现,具有强低频干扰的地震记录主频很低,而低频相干噪音弱的地震记录主频较高,分析地震资料中的主频分布规律,设定合适的门槛值,将主频低于门槛值的记录与其它记录分开,图1是不同主频的地震记录。
其中,获得地震信号主频的方法是将地震信号经过傅立叶变换到频率-振幅域,再对其求导,令其等于0,即振幅的最大值为此地震信号的主频,具体步骤如下:
假设地震信号为x(t),为连续时间信号,对其做傅立叶变换得到
其中,w为角频率w=2πf,f为频率;
X(ω)通常是复数,利用复数性质可得:
X(ω)=A(ω)eiφ(ω)
其中,A(ω)和φ(ω)分别为振幅谱和相位谱,它们是通过用下面的公式计算获得的:
其中,Xr(ω)和Xi(ω)为X(ω)的付里叶变换的实部和虚部,分别为:
Xr(ω)=A(ω)cosφ(ω),Xi(ω)=A(ω)sinφ(ω)
对X(ω)=A(ω)eiφ(ω)求导,令解此方程则得到的就是该信号主频。
(2)相干噪音特征分析
对强相干噪音记录进行分析,分析相干噪音的速度、频率分布范围以及与有效反射波特征差异,确定相干噪音的运动学特征。
(3)相干噪音压制
根据相干噪音的分布规律,将识别出来的强相干噪音记录通过FX域或均值加权噪音分离技术,形成噪音的模型,从原始地震记录中减掉噪音,得到去噪后的记录,实现提高信噪比的目标,同时使有效反射波能量不受损失。
本发明主要用于地表起伏剧烈复杂构造地区线性噪音衰减。利用本项发明技术在某地区山地地震资料处理中应用取得较好的效果,图3是该实施例中利用本发明方法去除强相干噪音的效果图。图3上半部分黑圈划定部分相干噪音能量较强,图3下半部分划定部分黑圈划定部分相干噪音得到了去除,经过该方法去除相干噪音后地震剖面上强线相干噪音得到压制,有效波没有受到损失。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (2)
1.一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法,其特征在于:所述方法首先将分布有强能量相干噪音的地震记录与分布有弱能量相干噪音的地震记录分开,然后对分布有强能量相干噪音的地震记录加强去噪,对分布有弱能量相干噪音的地震记录降低去噪程度甚至不去噪,使陡倾角有效反射波得到保护;
所述方法包括以下步骤:
(1)地震记录分类
分析地震资料中的地震信号主频的分布规律,设定合适的门槛值,将地震信号主频低于门槛值的定义为分布有强相干噪音的地震记录,将地震信号主频高于门槛值的定义为分布有弱相干噪音的地震记录,从而将两者分开;
(2)相干噪音特征分析
对分布有强相干噪音的地震记录进行分析,分析强相干噪音的速度、频率分布范围以及与有效反射波的特征差异,进而确定强相干噪音的运动学特征;
(3)相干噪音压制
根据步骤(2)获得的强相干噪音的特征,对步骤(1)识别出来的分布有强相干噪音的地震记录采用FX域或均值加权噪音分离技术,形成噪音的模型,然后从原始地震记录中减掉噪音,得到去噪后的地震记录。
2.根据权利要求1所述的地表起伏复杂构造地区地震资料分类相干噪音衰减方法,其特征在于:所述步骤(1)中获得地震信号主频的方法是将地震信号经过傅立叶变换到频率-振幅域,再对其求导,令其等于0,即振幅的最大值为此地震信号的主频,具体步骤如下:
假设地震信号为x(t),为连续时间信号,对其做傅立叶变换得到
其中,w为角频率w=2πf,f为频率;
X(ω)通常是复数,利用复数性质可得:
X(ω)=A(ω)eiφ(ω)
其中,A(ω)和φ(ω)分别为振幅谱和相位谱,它们是通过用下面的公式计算获得的:
其中,Xr(ω)和Xi(ω)为X(ω)的付里叶变换的实部和虚部,分别为:
Xr(ω)=A(ω)cosφ(ω),Xi(ω)=A(ω)sinφ(ω)
对X(ω)=A(ω)eiφ(ω)求导,令解此方程则得到的就是该信号主频。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110170561.2A CN102841380B (zh) | 2011-06-22 | 2011-06-22 | 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110170561.2A CN102841380B (zh) | 2011-06-22 | 2011-06-22 | 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102841380A CN102841380A (zh) | 2012-12-26 |
CN102841380B true CN102841380B (zh) | 2015-12-16 |
Family
ID=47368913
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110170561.2A Active CN102841380B (zh) | 2011-06-22 | 2011-06-22 | 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102841380B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375184B (zh) * | 2013-08-15 | 2017-04-05 | 中国石油天然气集团公司 | 一种高效的地震数据随机噪声衰减方法 |
CN104597501B (zh) * | 2013-11-01 | 2017-07-07 | 中国石油天然气集团公司 | 炮集记录处理方法及炮集记录处理装置 |
CN105093280B (zh) * | 2014-05-15 | 2018-04-06 | 中国石油化工股份有限公司 | 表层模型对地震数据影响的低频与高频成分的分解方法 |
CN105527649B (zh) * | 2015-11-30 | 2018-06-22 | 中国科学院地质与地球物理研究所 | 一种多域多次分离的高效采集多震源混合数据分离方法 |
CN107526107A (zh) * | 2016-06-21 | 2017-12-29 | 中国石油化工股份有限公司 | 强振幅相干噪音衰减方法及装置 |
CN111175826B (zh) * | 2020-01-15 | 2022-03-08 | 大庆油田有限责任公司 | 一种地震数据相干噪音压制质控方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
-
2011
- 2011-06-22 CN CN201110170561.2A patent/CN102841380B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
Non-Patent Citations (3)
Title |
---|
地震资料分频处理及应用;高爱荣等;《内蒙古石油化工》;20101231(第14期);11-12 * |
均值加权消除相干干扰;夏洪瑞等;《石油物探》;20000630;第39卷(第2期);77-82 * |
基于小波分频叠前相干噪声压制方法;覃天;《地球物理学进展》;20090831;第24卷(第4期);1426-1430 * |
Also Published As
Publication number | Publication date |
---|---|
CN102841380A (zh) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102841380B (zh) | 一种地表起伏复杂构造地区地震资料分类相干噪音衰减方法 | |
US11880011B2 (en) | Surface wave prediction and removal from seismic data | |
CN1129801C (zh) | 地震数据噪声的过滤方法 | |
Le Feuvre et al. | Passive multi-channel analysis of surface waves with cross-correlations and beamforming. Application to a sea dike | |
Cornou et al. | Contribution of dense array analysis to the identification and quantification of basin-edge-induced waves, Part II: Application to Grenoble basin (French Alps) | |
Strobbia et al. | Refraction microtremors: Data analysis and diagnostics of key hypotheses | |
Pamuk et al. | Soil characterization of Tınaztepe region (İzmir/Turkey) using surface wave methods and Nakamura (HVSR) technique | |
CN104614769B (zh) | 一种压制地震面波的聚束滤波方法 | |
CN103926623B (zh) | 一种压制逆时偏移低频噪音的方法 | |
CN103364832A (zh) | 一种基于自适应最优核时频分布的地震衰减定性估计方法 | |
Behm et al. | Love waves from local traffic noise interferometry | |
CN103091714A (zh) | 一种自适应面波衰减方法 | |
Zhang et al. | A case study of site effects on seismic ground motions at Xishan Park ridge in Zigong, Sichuan, China | |
CN107783191B (zh) | 多维空间时空时频峰值滤波消减地震勘探随机噪声的方法 | |
Ikeda et al. | Surface wave attenuation in the shallow subsurface from multichannel–multishot seismic data: a new approach for detecting fractures and lithological discontinuities | |
Dean et al. | The coherency of ambient seismic noise recorded during land surveys and the resulting implications for the effectiveness of geophone arrays | |
Morton et al. | Passive multichannel analysis of surface waves using 1D and 2D receiver arrays | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
Zhang et al. | Retrieval of shallow S-wave profiles from seismic reflection surveying and traffic-induced noise | |
Boaga et al. | Shear wave structural models of Venice Plain, Italy, from time cross correlation of seismic noise | |
CN115373023A (zh) | 一种基于地震反射和车辆噪声的联合探测方法 | |
Sivaram et al. | Stability assessment and quantitative evaluation of H/V spectral ratios for site response studies in Kumaon Himalaya, India using ambient noise recorded by a broadband seismograph network | |
Zhang et al. | Extraction of surface-wave phase velocities from ambient noise in the presence of local noise sources based on matched-field processing | |
CN105092343A (zh) | 一种基于叠前道集去除薄层调谐效应的方法 | |
CN103364834A (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 |