数字处理－3.ppt

1、第三章反褶积, 3.1 反褶积的概念,左图是 反射系数序列与实际地震记录频谱比较，反射系数频谱明显特征是频谱很宽，称为白谱。而地震记录的频带范围是有限的（带通）。,图3-2 反射系数与地震记录剖面的比较 上为 反射系数剖面，下为 地震剖面, 3.1 反褶积的概念, 3.1 反褶积的概念, 3.1 反褶积的概念,图3-3 一道地震记录, 3.1 反褶积的概念, 3.1 反褶积的概念, 3.1 反褶积的概念,3.2 最佳维纳滤波,3.2 最佳维纳滤波,3.2 最佳维纳滤波,3.2 最佳维纳滤波,3.2 最佳维纳滤波,3.2 最佳维纳滤波,3.2 最佳维纳滤波,3.3 最小平方反褶积,3.3 最小平

2、方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,3.3 最小平方反褶积,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理,二有干扰时地震记录的最小平方反褶积原理

3、,三 最小平方反褶积中的预白化处理,三 最小平方反褶积中的预白化处理,三 最小平方反褶积中的预白化处理,三 最小平方反褶积中的预白化处理,三 最小平方反褶积中的预白化处理, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积,作业:证明3-4-24式, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积, 3.

4、4 预测反褶积,图3-13 海上一级鸣震的产生示意图, 3.4 预测反褶积, 3.4 预测反褶积,鸣震主频为37Hz左右，周期接近25ms 。, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积,图3-15 海上二级鸣震的产生, 3.4 预测反褶积, 3.4 预测反褶积, 3.4 预测反褶积,褶积因子的长度不变，取为128ms，白噪系数为0%，采样间隔为2ms，对不同预测长度进行试验。图中(a)为反射系数序列，(b)为合成地震记录，上面是不同步长的反褶积输出。比较输出的自相关和振幅谱，可以看出，当预测步长为2ms时，反褶积结果最好。, 3.4 预测反褶积,图(a)为输入道，滤波算

5、子长度为160ms，白噪系数为0.1，(b)(e)为反褶积输出,预测步长分别为(b) 12ms, (c) 32ms, (d) 64ms, (e) 128ms,可以看出(b)的效果最好,(e)基本上不起作用。, 3.4 预测反褶积, 3.4 预测反褶积,预测步长为2ms，白噪为零。预测长度太短（20ms和44ms时）尾部抖动现象十分明显，长度越大，反褶积效果越好，但达到一定长度时（图中128ms），反褶积效果的变化就十分平缓了。, 3.4 预测反褶积,预测步长为4ms，白噪为0.1%。预测因子长度:(b) 40ms, (c) 80ms, (d) 160ms, (e) 240ms , 图中到160

6、ms时，地震子波得到压缩，交混回响已被较好压制，在1.1, 1.35, 1.85, 2.15s处的主要反射没有重叠了。, 3.4 预测反褶积,图3-20 用理论模型试验白噪量,预测步长2ms，因子长度128ms。不加白噪和白噪系数很小（0.1%）时，反褶积输出基本不变，当白噪加到很大时，尾部抖动现象就较重。, 3.4 预测反褶积,图3-21是实际地震剖面的实验 反滤波因子长度160ms，预测步长为4ms(脉冲反褶积)。白噪系数:(b)1 ,(c)4 ,(d) 16, (e) 32 可以看出，噪系数较小时（ 1% ），反褶积效果较好。, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波

7、整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积,零相位子波,最小相位子波, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积,图3-23 子波与反子波时域分布关系, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积, 3.5 子波提取与子波整形反褶积,(a)仪器的脉冲响应w(t)（地震子波） (b)最佳整形反褶积算子，实线对应零相位子波，虚线对应脉冲 整形 (c)为实际的最佳输出，实线和虚线对应关系相同； (d)为不同延迟的最小平方反褶积的误差能量，零相位子波整形处理误差能量较脉冲整形小得多，而且零相位子波，在延迟为60 ms时实现了最佳反褶积。, 3.5 子波提取与子波整形反褶积,图3-26 子波整形反褶积前后资料的比较, 3.5 子波提取与子波整形反褶积,对两种不同震源资料进行整形反褶积, 3.6 同态反褶积, 3.6 同态反褶积, 3.6 同态反褶积, 3.6 同态反褶积, 3.6 同态反褶积, 3.6 同态反褶积, 3.6 同态反褶