地球物理资料数字处理(第二十二十一讲).doc

《地球物理资料数字处理(第二十二十一讲).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第四节最小熵反褶积(MinimumEntropyDeconvolution)地球物理技术具有高科技密集，更新速度快的特点。每一次新技术的出现及其应用都对数字处理技术产生巨大的冲击。尤其是随着现代信息处理技术的发展和计算机技术的进步，更加促进了新技术在地震资料数字处理中的应用。最小熵反褶积就是将信息论中熵的概念引入到反褶积中的技术。与最小平方反褶积不同的是，它不需要假定反射系数是白噪，子波为最小相位的条件，从这个意义上讲，它比最小平方反褶积方法先进，从另外的角度看，它的算法较为复杂，这又是这种方法的缺点。这里讨论最小熵反褶积技术，意义不在于其方法本身。目的是把反褶积技术引入更广阔的领域

2、。最小熵反褶积的目的是压缩地震子波，提高地震记录的分辨率。本节要点：l熵的基本概念；l最小熵反褶积原理；l最小熵反褶积原理的步骤。一、熵的基本概念(TheBasicConceptofEntropy)研究最小熵反褶积，必须首先理解熵的概念。熵，最早并不是地震资料数字处理技术中的专有名词。它是由物理学中热学的概念借用而来。在物理热力学中，研究分子运动状态，用熵来衡量分子热运动的秩序。从分子运动的观点看，分子热运动要从有序趋于混乱，这时分子的熵发生变化，熵越大，分子运动越无秩序，规律性不强，反之熵越小，分子运动越有秩序。可见，最早熵是用来衡量分子运动秩序的物理量。在信息论中，熵用来衡量系统

3、内信息的不确定性或不可预测性，熵越大系统内信息的不确定性就越大，反之，则不确定性就越小二、最小熵反褶积原理（ThePrincipleOFMinimumEntropyDeconvolution）根据地震褶积模型，地震记录，最小熵反褶积问题归结为1．寻找一线性（LinearDeconvolutionFactor）反褶积算子2．使得地震记录与反褶积算子褶积后其输出为(WaveShape)波形简单，(Polarity)和延迟时间(TimeDelay)未知的尖脉冲序列。即输出为反射系数序列，则有下式成立（5-4-1）从熵的角度分析，最小熵反褶积的输出准则是希望得到褶积输出是由一系列极性和延迟时

4、间未知（由地下介质的性质决定）的尖脉冲序列组成，且应使信号的秩序最好，规律性强，亦即熵最小。因此，最小熵反褶积的关键是如何使褶积的输出达到最小。为了提高反褶积算子的精度，实际确定反褶积算子时，不是使用一道地震记录，而是选区一组地震道组成一个道集，如图（5-11）。这时，反褶积的输出亦为一个道集，则（5-4-2）其中，求取反褶积因子的道集，为道序号,为记录样点采样序号(,m为反褶积因子的长度,下面的问题是如何衡量输出道集的熵达到最小。为了确定输出波形达到熵最小，构造函数，成为最大方差范数图5-11确定最小熵反褶积因子输入道集（5-4-3）（5-4-4）可以证明分析最大方差范数的特点，若

5、取输出为尖脉冲时，即输出波形最为简单，亦即熵最小的情况，则，若输出为幅度相等的2个尖脉冲时，V=0.50。若输出为10个幅度相等的尖脉冲时，V=0.10。虽然，当脉冲幅度相等时，脉冲的个数越多，则最大方差范数V的值越小，且最大方差范数不受脉冲的极性和出现时间影响。对于输出波形为幅度不等的尖脉冲时，可以参看图（5-12）了解最大方差范数的特点。可以理解，最大方差范数V可以作为衡量输出信号熵大小的准则。下面讨论如何确定反褶积因子。根据上面建立的最小熵反褶积准则及最大方差范数的特点，所求的反褶积因子必满足使最大方差范数达到最大，也就是输出地震记录的熵最小，所以必须满足下式（5-4-5）图5

6、-12简单脉冲序列最大方差值数V的特征所以(5-4-6)为简化求导，令，则（5-4-6）可写为（5-4-7）由（5-4-2）式知则（5-4-7）又可写成下式对上式整理，又可写成下式令将代入（5-4-7）得（5-4-8）写成矩阵形式（5-4-9）分析（5-4-7）式知，该方程组并不是一个线性方程组，而是一个复杂的非线性高次方程，将（5-4-7）式写（5-4-8）式制式起到使其表达式简洁的作用，因此不能用求解线性方程组的方法求解。可以采用迭代法求解。迭代求解过程为1）给定一初始线性褶积算子2）根据计算值。3）根据计算和的值4）将计算所得的和的值代入矩阵（5-4-9）并求解，可得一新的褶积

7、算子5)重复(2)~(4)步骤,直至得到满意的输出为止,这时得到的可做为反褶积算子。上述的迭代过程一般4~6次就可以结束。由于所解方程为非线性，这时输出的反褶积算子并不能保证最大方差范数V是唯一的极大值，但此时的反褶积算子仍不失其使用意义。在迭代的过程中为确保最小熵准则的一致性，必须对每次迭代所得的反褶积算子进行归一化处理，以使每次计算的最大方差范数V是使用统一的数量标准比较。在建立最小熵反褶积准则时，我们仅考虑输出具有简单外形的条件，并没有考虑输出脉冲的