数字处理第三章课程方案.ppt

　　地表一致性反褶积的主要目标是校正子波的振幅谱，所体，这种方法适合于地表条件变化大的地区．但是真正方法并不是着重展宽频谱，所以，分辨率并不能提高很多． 3.9 反Q滤波及谱白化 由于地层的吸收作用,地震波经地层传播后,能量被衰减损耗,频率变低,特别是深层,分辨率大大下降.因此,为恢复地震波原能量,必须做吸收补偿,即Q补偿或滤波. 1、Q因子的物理含义 波在介质中传播，产生的弹性能逐渐被介质吸收，最后转换为热能，其过程称为吸收。一个波长内，原地震波能量与传输所损耗的能量之比为Q。 将 3-124 式泰勒展开，舍弃高次项，得 2、Q值的估算 1 振幅包络法； 2 李氏经验公式： 3、反Q滤波的方法实现 设反Q滤波器为最小相位，振幅谱为 求出滤波因子，便可滤波。 * 或 MATLAB 程序： 三、预测反褶积压缩反射脉冲（自学） 四、参数选择 （1）预测步长，关键参数，效果与此有关； （2）预测因子长度，要适宜，小，波形尾部波动；大，增加运算量，相邻反射及噪声影响； （3）预白化量，增加求解的稳定性，需根据噪声水平来定。 后有三图，以助说明。 预测因子长度128ms，预测步长2ms，效果最好，其余有延续时间（非脉冲信号）。 （a）为脉冲序列，（b）为合成记录 预测步长2ms,预测因子长度128ms，效果好。小长度因子，尾部波动明显 预白化量0.1，效果较好；预白化量增大，尾部波动增大，效果变差。 3.6 子波整形反褶积 １．子波的相位与分辨率 　　子波的相位通常有三种，即最小相位，混合相位和最大相位，这些子波是单边的物理可实现信号． 海上勘探和陆上勘探爆炸震源产生的地震子波都接近最小相位，所以，经常假设和讨论最小相位．还有一种相位的子波是零相位子波，虽说是一种非因果信号，是物理不可实现的，但在数字滤波、反褶积和反演中经常用到。 为什么要用到零相位子波，基于以下两个考虑： ①陆上可控震源子波因是通过自相关处理得到的，所以是零相位的； ②零相位子波，比其他相位子波的分辨率高． 　 ２．子波与反子波的时域分布的特点 　对于褶积公式 　 　假设子波为有限长度， 　 　从理论上来说，我们可以得到无限长的子波 　利用Ｚ变换，我们可以推导出有限长子波与所求反子波的关系 ３.相位对反褶积精度的影响 　　　在前面的讨论中，为了讨论问题方便，都假设了子波为最小相位．实际上子波有不同的相位．在反褶积中，其他条件相同，相位不同，反褶积的精度是不同的． 　　对于最小平方反褶积，滤波方程为 　我们分两种相位的期望输出来讨论反褶积问题． 　１）期望输出为（１，０，０） 　这时的方程为 解得反子波为 ，反褶积输出为 与期望输出（１，０，０）的误差能 量为　　． 　从中我们也可以看到，次时的误差比最小相位是１／２１大多了． ２）期望输出为（0,1,0） 这是的滤波方程为 得到的反子波为 ，反褶积输出为 这与期望输出 0,1,0 的误差能量 为　　．冲中我们可以看到此时的反褶积比上 面的精度要高 由此我们可以得到一个重要的结论： ①子波振幅相同时，最小相位子波对期望输出为零延迟的反褶积，误差最小； ②在子波为混合相位和最大相位时，期望输出的相位应与子波的相位相匹配，有一个最佳延迟，只有这样才能得到合适的反褶积结果． ４．子波整形反褶积 　　　上面关于 -0.5,1 和期望输出关系的讨论，实际上已经牵涉到子波整形反褶积问题，即如何把不同相位的子波转变成最佳子波的反褶积． PPT未做完 3.8 地表一致性反褶积 目的：消除由于近地表变化对地震子波波形的影响。 地表一致性谱分解（Taner 1981）认为任意道的频谱是对应的炮点、接收点、共中心点和炮检距的频谱响应的乘积，对其取对数，可分解炮点和接收点的振幅响应和相位响应。振幅响应被用于振幅一致性补偿，相位响应被用于时差校正。 地表一致性谱分解的基础是褶积模型和预测原理，并扩展为多道平均，把等效子波分解成多个分量。 地表一致性模型： 其中： 求 频谱 ： 其振幅谱和相位谱分别为： 对振幅谱取对数: 设 是实际地震道的振幅谱，根据地表一致性的假设 ，一个地震道由相应分量组成，而每一分量对相关的地震记录的贡献是一样的，则实际振幅谱与模型估计振幅谱之差的最小平方和为： 使误差能量达到最小，得到关于4个分量的方程组。 用迭代法求解取每个分量，在最小相位条件下，估算它们的反算子，进行反褶积。 PPT未做完，想法是将书中复杂内容用精练语言来表达。 与单道反褶积的区别：炮点（或接收点）分量里仅含与观测点有关因素，对相关的所有道是相同的。 一．方法原理 　１．褶积模型 　对于地震道有 　式中w t 为一中和子波，其他符号意义同前． 　考虑到地表四个因素： 共炮点 共接收点 共中心点 共炮检距点 （3-118）