地震数据处理第三章：反褶积.ppt

第三章 反褶积 3.1、反褶积及褶积模型 3.2、反滤波 3.3、最佳维纳滤波与最小平方反褶积 3.4、脉冲反褶积 3.5、预测反褶积 3.6、子波整形反褶积 3.7、同态反褶积 3.8、地表一致性反褶积 3.6 子波整形反褶积 １．子波的相位与分辨率 　　子波的相位通常有三种，即最小相位，混合相位和最大相位，这些子波是单边的物理可实现信号． 海上勘探和陆上勘探爆炸震源产生的地震子波都接近最小相位，所以，经常假设和讨论最小相位．还有一种相位的子波是零相位子波，虽说是一种非因果信号，是物理不可实现的，但在数字滤波、反褶积和反演中经常用到。 为什么要用到零相位子波，基于以下两个考虑： ①陆上可控震源子波因是通过自相关处理得到的，所以是零相位的； ②零相位子波，比其他相位子波的分辨率高． 　 ２．子波与反子波的时域分布的特点 　对于褶积公式 　 　假设子波为有限长度， 　 　从理论上来说，我们可以得到无限长的子波 　利用Ｚ变换，我们可以推导出有限长子波与所求反子波的关系 ３.相位对反褶积精度的影响 　　　在前面的讨论中，为了讨论问题方便，都假设了子波为最小相位．实际上子波有不同的相位．在反褶积中，其他条件相同，相位不同，反褶积的精度是不同的． 　　对于最小平方反褶积，滤波方程为 　我们分两种相位的期望输出来讨论反褶积问题． 　１）期望输出为（0，０，1） 　这时的方程为 　解得反子波为　　　　　　　　，反褶积输出为 　　　　　与期望输出（0，0，1）的误差能 量为　　． 　从中我们也可以看到，此时的误差比最小相位是１／２１大多了． ２）期望输出为（0,1,0） 这是的滤波方程为 得到的反子波为　　　　　　　，反褶积输出为 　　　　　 这与期望输出(0,1,0)的误差能量 为　　．式中我们可以看到此时的反褶积比上 面的精度要高　　　　　 对（3-43）式描述的是频谱形式, 对其改造为 （3-43’） 为了解决带限问题,在地震信号的功率谱P(ω)中,从低频到高频统一加一白噪声。 2.预白化处理 （3-44） 当反射系数为白噪声时，子波的自相关=地震记录的自相关。 从以上方程可以看出： 第三章 反褶积 3.1、反褶积及褶积模型 3.2、反滤波 3.3、最佳维纳滤波与最小平方反褶积 3.4、脉冲反褶积 3.5、预测反褶积 3.6、子波整形反褶积 3.7、同态反褶积 3.8、地表一致性反褶积 3.4 脉冲反褶积 一、脉冲反褶积原理 当期望输出为 时， 除 外，其余全为零。 方程（矩阵形式）： （3-80） （3-79） 脉冲反褶积（对应3-80式）要求输入子波为最小相位的，实际中往往是非最小相位的，因此要作某种优化处理。 最小平方反滤波，对子波相位已无要求，可以任意。 反滤波因子在理论上应为无穷项，实际只能取有限项，并且其主要部分 反滤波因子的形状和位置要根据地震子波的形状来决定。 详情见书p78 最小延迟，m0=0，m 最大延迟，m0=-n-m 混合延迟，m0=-m1 二、参数选择 （效果好的标准是什么？子波是否被压缩） （1）反滤波因子长度m，通过试验来选择； （2）相关时窗长度m+n，常为2m~记录长度； （3）稳定常数，一般取 的百分数。 干扰小，取0.005~0.01 ； 干扰大，取0.02~0.05 。 下图是延迟脉冲反褶积效果分析图 （误差能量） 延迟时间 输入为非最小相位子波，期望输出 第三章 反褶积 3.1、反褶积及褶积模型 3.2、反滤波 3.3、最佳维纳滤波与最小平方反褶积 3.4、脉冲反褶积 3.5、预测反褶积 3.6、子波整形反褶积 3.7、同态反褶积 3.8、地表一致性反褶积 3.5 预测反褶积 一、预测滤波原理 预测问题是已知某个物理量的过去值和现在值,通过对已知信息加工处理来获得未来某个时刻的预测值。其数学描述如下: 设 为现在值，　　　　　 为过去值, 如果定义预测步长为 ，用现在值和过去值来预测将来 时刻的预测值 ，即 （3-83） 使预测值与实际未来值的误差—预测误差 最小，按最小平方原理，得如下方程： （3-86） 矩阵形式为： （3-