【技术实现步骤摘要】
一种地震波传播能量衰减补偿方法
本专利技术涉及地震勘探数据处理
,更具体的说是涉及一种地震波传播能量衰减补偿方法。
技术介绍
在石油勘探过程中,地层对地震波的吸收造成了震源子波时变,使资料的分辨率降低,特别是深层信号,由于地层的低通滤波作用,高频信号衰减严重。因此,地表记录到的不同深度反射信号及其频带,已经不能反映地下的真实情况,影响了后续地震资料的处理和解释工作,也影响了石油勘探结果的准确性。为了提高地震资料的分辨率,就必须对地震波衰减进行吸收补偿处理。目前,地震数据处理中采用的地震波衰减补偿方法主要有四种:(1)常规反Q滤波方法,该方法是一种通过级数展开近似补偿高频成分的方法;(2)截止频率法(增益控制法)反Q滤波,该方法通过在地震波有效截止频率外固定振幅补偿函数的补偿增益限来压制高频噪音,进行振幅补偿;(3)稳定因子法反Q滤波,通过引入稳定因子(和振幅补偿增益限有关的很小的常数),使反Q滤波振幅补偿函数在频率域满足高斯分布,使振幅补偿函数稳定且具有固定增益限,即振幅增益有固定的最大值,从而实现振幅补偿;(4)自适应增益限反Q滤波,该方法针对常规反Q滤波方法不能压制高频噪音以及稳定因子法反Q滤波增益限固定不变的缺陷,采用分段函数的思路,在地震波有效频带内采用常规反Q滤波补偿公式,在地震波有效频带外则采用稳定因子法反Q滤波公式,地震波有效频带最大值即为稳定因子法反Q滤波振幅补偿函数最大值对应的角频率,即截止角频率,该角频率随地震记录各道采样点的变化而变化。但是,常规反Q滤波方法没有考虑到振幅补 ...
【技术保护点】
1.一种地震波传播能量衰减补偿方法,其特征在于,包括以下步骤:/n输入实际地震资料,进行广义S变换将实际地震资料从时域变换到频域,对实际地震资料求取时变的品质因子,根据时变增益限稳定因子法反Q滤波方程进行反Q滤波,得到反Q滤波结果实现衰减补偿;/n所述时变增益限稳定因子法反Q滤波方程的获取步骤为:/nS1:求取品质因子和振幅补偿函数增益限的关系;/nS11:计算品质因子和瞬时能量的关系式;/nS12:根据地震波能量计算公式计算不同时刻地震波振幅比和能量比的关系式;/nS13:将S11中得到的品质因子和瞬时能量的关系式与S12中得到的不同时刻地震波振幅比和能量比的关系式进行结合,得到品质因子和振幅补偿函数增益限的关系式;/nS2:计算振幅补偿函数增益限和稳定因子的关系;/nS3:通过S1中的品质因子和振幅补偿函数增益限的关系和S2中振幅补偿函数增益限和稳定因子的关系,计算出品质因子和稳定因子的关系,用时变的品质因子替代时不变的稳定因子,得到时变增益限稳定因子法反Q滤波方程,具体方程为:/n
【技术特征摘要】
1.一种地震波传播能量衰减补偿方法,其特征在于,包括以下步骤:
输入实际地震资料,进行广义S变换将实际地震资料从时域变换到频域,对实际地震资料求取时变的品质因子,根据时变增益限稳定因子法反Q滤波方程进行反Q滤波,得到反Q滤波结果实现衰减补偿;
所述时变增益限稳定因子法反Q滤波方程的获取步骤为:
S1:求取品质因子和振幅补偿函数增益限的关系;
S11:计算品质因子和瞬时能量的关系式;
S12:根据地震波能量计算公式计算不同时刻地震波振幅比和能量比的关系式;
S13:将S11中得到的品质因子和瞬时能量的关系式与S12中得到的不同时刻地震波振幅比和能量比的关系式进行结合,得到品质因子和振幅补偿函数增益限的关系式;
S2:计算振幅补偿函数增益限和稳定因子的关系;
S3:通过S1中的品质因子和振幅补偿函数增益限的关系和S2中振幅补偿函数增益限和稳定因子的关系,计算出品质因子和稳定因子的关系,用时变的品质因子替代时不变的稳定因子,得到时变增益限稳定因子法反Q滤波方程,具体方程为:
其中,φ(τ,ω)为时变增益限稳定因子反Q滤波方程,为常规反Q滤波振幅补偿函数,σ(τ)2为τ时刻的稳定因子。
2.根据权利要求1所述的一种地震波传播能量衰减补偿方法,其特征在于,S11包括以下内容:
(1)根据离散时间信号的能量表达式,对地震波信号做广义S变换,得到地震波传播过程中能量在时频域的展布特征,也就是由时频振幅谱表示的地震波能量表达式,具体公式为:
E(τ,ω)=[GSTN(τ,ω)]2-[GSTN(τ+Δτ,ω)]·[GSTN(τ-Δτ,ω)]公式(2)
其中,GST为广义S变换后的时频振幅谱;
(2)根据步骤(1)中的由时频振幅谱表示的地震波能量表达式,求取地震波在某时刻的瞬时能量;地震波瞬时能量为频率域能量的最大值,具体为:
<...
【专利技术属性】
技术研发人员:阎建国,宋鑫磊,陈琪,王博睿,张雪纯,邓儒炳,张建团,
申请(专利权)人:成都理工大学,
类型:发明
国别省市:四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。