关于反褶积的问题_周兴元.ppt

的。对非小相位地震子波应用脉冲反褶积，其结果是失真的。正因为如此有关混合相位反褶积的文章不时地出现。混合相位反褶积方法也比较多。下面就我所了解有些公司已经在用的混合相位反褶积方法做一介绍。 目 录 一、对地震记录的模拟 二、关于地震子波的相位特性与相位谱 三、反褶积 四、混合相位反褶积 五、纯相位滤波（或称子波相位校正） 六、时变反褶积 目 录 一、对地震记录的模拟 二、关于地震子波的相位特性与相位谱 三、反褶积 四、混合相位反褶积 五、纯相位滤波（或称子波相位校正） 六、时变反褶积 这类反褶积主要有同态反褶积，最小熵反褶积和 反褶积等。 1.同态反褶积 同态反褶积是由公式(1)子波时不变的褶积模型出发，在频率域为 (11) 取对数 反傅氏变换得到 (12) 叫复赛谱。 四、混合相位反褶积 一般认为在 的原点附近，主要是子波部分，而远离原点的部分为反射系数部分。现有程序做法是在给定时空范围内，随机地取小时窗，计算复赛谱并平均，目的是为了使 ，从而得到子波的复赛谱，经逆运算得到时域的子波和反子波。人工选定一个反子波做反褶积。现在改进为，在若干个反子波中，每一个都对该时空范围内记录做反褶积，熵最小者被选中并应用。 交互同态反褶积后叠加剖面 交互同态反褶积前叠加剖面 子波 反子波 交互同态反褶积前频谱 交互同态反褶积后频谱 中部地区一条煤田叠加剖面 经精细处理叠后同态褶积剖面 2.最小熵反褶积 最小熵反褶积是直接求取混合相位反子波 并计算输出 ，使得目标函数（方差模）： (13) 这个目标函数是最早提出最小熵反褶积的作者提出来的，后来又有多种变形的目标函数。由于目标函数为非二次型, 所以不能用最小平方法得到唯一解，而是给定 的初值，用迭代法求解。 反褶积的目标函数，是使反褶积输出结果 满足 反褶积因子 也是通过迭代法求解。 这里没有将 反褶积单独作为一个小标题，是因为它不健全。对于相位而言，我们希望剩余子波（如果有的话）是零相位的。那是否零相位的 比非零相位的 小呢？ 我用1003个随机数取长度为3个点的小时窗，逐点滑动，得到1000组长度为3个点的序列，每组计算自身褶积的 范数 和自相关的 范数 ，发现有4组 ，其他 。例如3,2,-2， ， ；而对8,-6,-1则 ， ，相互间大小关系正好与上例相反。所以说它是不健全的。 但是有的公司，如前Geco-prakla将 反褶积作为特殊处理手段来应用，所以这里还是得提一下。 目 录 一、对地震记录的模拟 二、关于地震子波的相位特性与相位谱 三、反褶积 四、混合相位反褶积 五、纯相位滤波（或称子波相位校正） 六、时变反褶积 目 录 一、对地震记录的模拟 二、关于地震子波的相位特性与相位谱 三、反褶积 四、混合相位反褶积 五、纯相位滤波（或称子波相位校正） 六、时变反褶积 纯相位滤波，主要是针对子波的相位特性与反子波的相位特性不匹配时，比如对混合相位子波做小相位反褶积，对反褶积结果再做一次纯相位滤波，目的是使剩余子波尽可能向零相位靠近。 纯相位滤波到目前为止主要有两种，一种是常数相位滤波，另一种是非常数相位滤波。纯相位滤波都要用某一个目标函数，在给定步长的情况下，经过扫描选定使目标函数达到最佳者。 五、纯相位滤波 （或称子波相位校正） 关于目标函数有两种： 一是方差模准则： 式中 为滤波的输出。若将这一目标函数，应用在一段记录上，其中强反射波起主要作用，弱反射则可能产生畸变。 二是Parsimony准则： 其中 为滤波输出， 。就目标函数而言，是想使弱反射畸变程度尽可能地小。 1.常数相位滤波 若在频率域，可表示为： (14) 其中 为常数相位角。 在给定扫描步长 时，在－90o到90o范围内扫描，取使目标函数达到最佳者的相角并应用。 好像Larner在一篇文章中曾提到，零相位子波做完小相位反褶积后，相位接近90o。 2.非常数相位滤波 如前所述，若混合相位子波的