1.一种时变地震衰减模型的构建方法,其特征在于依据物理原理而提出的根据wang’s子波变化特征估算地震衰减系数的全新方法,并通过约束反演衰减模型不同波长系数的方式,构建稳定的时变地震衰减模型。
2.根据权利要求1所述时变地震衰减模型的构建方法,其特征在于实现步骤是:(1)生成地震数据时频谱,估算各时窗的功率谱基本统计特性;(2)估算wang’s广义子波参数;(3)估算各时窗的地震衰减系数;(4)设置傅立叶级数拟合时变地震衰减模型;(5)约束反演衰减模型傅立叶级数的系数,构建时变的地震衰减模型。
3.根据权利要求2所述时变地震衰减模型的构建方法,其特征在于步骤(1)估算地震数据的时频谱;采用伽柏变换估算地震数据的时频谱,并针对实际地震数据存在的信噪比问题,对时频谱进行加权平均以改善时频谱的信噪比,改善依此得到的地震衰减模型的可靠性;估算各时窗地震数据的功率谱基本统计特性,即均值频率和频率方差。
4.根据权利要求2所述时变地震衰减模型的构建方法,其特征在于步骤(2)估算wang’s广义子波参数,即固有频率和分数阶值,分为两步迭代实现;首先根据每个时窗的频谱依次估算子波的固有频率和分数阶值;然后对全部固有频率值采用中值滤波获得全局化的子波固有频率,并根据全局化的子波固有频率重新推算对应于每个时窗的子波分数阶值。
5.根据权利要求2所述时变地震衰减模型的构建方法,其特征在于步骤(3)估算各时窗的地震衰减系数,本发明提出根据wang’s子波变化特征估算地震衰减系数的全新表达式如下:
式中α是地震信号的衰减系数,f0是子波的固有频率,τ是地震信号在介质中两点之间的传播时间,(ua,ub)是两时间位置的地震子波分数阶值。
6.根据权利要求5所述估算各时窗的地震衰减系数,其特征在于步骤(3)依据物理原理提出的地震衰减系数全新表达式也可近似表示为:
因此,估算地震衰减系数的表达式具有以下两个特征:首先,地震衰减系数跟(ua-ub)/τ成近似线性关系;其次,两时间位置的地震子波分数阶值(ua,ub)需要满足下列条件:
7.根据权利要求2所述时变地震衰减模型的构建方法,其特征在于步骤(4)设置傅立叶级数拟合时变地震衰减模型,采用傅立叶级数拟合各时窗的地震衰减系数,构建初始的时变衰减模型α(τ):
式中∆k是波数增量,{c0,ci,si}是傅立叶级数的系数;本发明采用不同波长及其系数表征时变衰减模型,将时变衰减模型α(τ)转换成2m+1个傅立叶级数的系数,z={c0,ci,si,i=(1,m)}。
8.根据权利要求2所述时变地震衰减模型的构建方法,其特征在于步骤(5)约束反演衰减模型傅立叶级数的系数,构建时变的地震衰减模型;根据各时窗估算的衰减系数可得到两组平均衰减系数,一组是根据两个子波直接估算得到的平均衰减系数ᾶ(τ),另一组是根据层间衰减系数用rms计算得到的平均值ᾶrms(τ);依照上述两组平均衰减系数构建反演问题,通过反演衰减模型不同波长系数的方法,将反演问题转从一个欠定问题转化成超定问题,从而最终获得稳定的时变地震衰减模型。
9.根据权利要求8所述通过反演构建时变地震衰减模型的方法,其特征在于提出反演傅立叶级数系数的方法,并提出反演衰减模型傅立叶级数系数的如下约束反演方法:
式中(ts,te)分别是纳入反演的平均衰减系数最浅时间和最深时间,z0是权利要求7中给出的层间衰减系数的初始估算模型,用作反演计算的约束条件,μ是其约束系数;根据上述反演结果z最终构建稳定的地震衰减模型z→α(τ)。