基于复倒谱匹配的地震子波估计方法

【作者】 李鲲鹏； 张学工； 李衍达；

【机构】 清华大学自动化系；

【摘要】 <正>反射波地震勘探的地震记录与地震子波和反射系数序列之间的关系可用一褶积模型表示:X₁=f₁·r₁+n_t （1）这里f_t为地震子波,r_t为反射系数序列,n_t为噪音.我们感兴趣的是r_t.估计和消除地震子波的影响成为地震数据处理流程中所必需的一个步骤.为了估计地震子波,通常假设地震子波是最小相位的而反射系数序列是具有零均值的高斯独立同分布随机过程.实际上,地震子波不是最小相位的而是混合相位的,反射系数序列也不是高斯独立同分布的随机过程,而是非高斯的且具有相关性.同态滤波法是可以估计混合相位地震子波的几种方法之一,该方法假设地震子波的振幅谱是光滑的,因而其复倒谱集中分布在原点附近;而反射系数序列的振幅谱是震荡的,其复倒谱分布在远离原点处.可以用一低通时窗或一高更多还原