地球物理资料数字处理第二十二十一讲.docVIP

最小熵反褶积 (Minimum Entropy Deconvolution) 地球物理技术具有高科技密集，更新速度快的特点。每一次新技术的出现及其应用都对数字处理技术产生巨大的冲击。尤其是随着现代信息处理技术的发展和计算机技术的进步，更加促进了新技术在地震资料数字处理中的应用。最小熵反褶积就是将信息论中熵的概念引入到反褶积中的技术。与最小平方反褶积不同的是，它不需要假定反射系数是白噪，子波为最小相位的条件，从这个意义上讲，它比最小平方反褶积方法先进，从另外的角度看，它的算法较为复杂，这又是这种方法的缺点。这里讨论最小熵反褶积技术，意义不在于其方法本身。目的是把反褶积技术引入更广阔的领域。 最小熵反褶积的目的是压缩地震子波，提高地震记录的分辨率。 本节要点： 熵的基本概念； 最小熵反褶积原理； 最小熵反褶积原理的步骤。 熵的基本概念(The Basic Concept of Entropy) 研究最小熵反褶积，必须首先理解熵的概念。熵，最早并不是地震资料数字处理技术中的专有名词。它是由物理学中热学的概念借用而来。在物理热力学中，研究分子运动状态，用熵来衡量分子热运动的秩序。从分子运动的观点看，分子热运动要从有序趋于混乱，这时分子的熵发生变化，熵越大，分子运动越无秩序，规律性不强，反之熵越小，分子运动越有秩序。可见，最早熵是用来衡量分子运动秩序的物理量。 在信息论中，熵用来衡量系统内信息的不确定性或不可预测性，熵越大系统内信息的不确定性就越大，反之，则不确定性就越小 最小熵反褶积原理（The Principle OF Minimum Entropy Deconvolution） 根据地震褶积模型，地震记录，最小熵反褶积问题归结为 1．寻找一线性（Linear Deconvolution Factor）反褶积算子 2．使得地震记录与反褶积算子褶积后其输出为(Wave Shape)波形简单，(Polarity)和延迟时间(Time Delay)未知的尖脉冲序列。即输出为反射系数序列，则有下式成立 （5-4-1） 从熵的角度分析，最小熵反褶积的输出准则是希望得到褶积输出是由一系列极性和延迟时间未知（由地下介质的性质决定）的尖脉冲序列组成，且应使信号的秩序最好，规律性强，亦即熵最小。因此，最小熵反褶积的关键是如何使褶积的输出达到最小。 为了提高反褶积算子的精度，实际确定反褶积算子时，不是使用一道地震记录，而是选区一组地震道组成一个道集，如图（5-11）。这时，反褶积的输出亦为一个道集，则 （5-4-2） 其中，求取反褶积因子的道集，为道序号,为记录样点采样序号(,m为反褶积因子的长度,下面的问题是如何衡量输出道集的熵达到最小。 为了确定输出波形达到熵最小，构造函数，成为最大方差范数 图5-11 确定最小熵反褶积因子输入道集 （5-4-3） （5-4-4） 可以证明 分析最大方差范数的特点，若取输出为尖脉冲时，即输出波形最为简单，亦即熵最小的情况，则，若输出为幅度相等的2个尖脉冲时，V=0.50。若输出为 10个幅度相等的尖脉冲时，V=0.10。虽然，当脉冲幅度相等时，脉冲的个数越多，则最大方差范数V的值越小，且最大方差范数不受脉冲的极性和出现时间影响。对于输出波形为幅度不等的尖脉冲时，可以参看图（5-12）了解最大方差范数的特点。可以理解，最大方差范数V可以作为衡量输出信号熵大小的准则。 下面讨论如何确定反褶积因子。根据上面建立的最小熵反褶积准则及最大方差范数的特点，所求的反褶积因子必满足使最大方差范数达到最大，也就是输出地震记录的熵最小，所以必须满足下式 （5-4-5） 图5-12简单脉冲序列最大方差值数V的特征 所以 (5-4-6) 为简化求导，令，则（5-4-6）可写为 （5-4-7） 由（5-4-2）式知则（5-4-7）又可写成下式