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

-第四节第五节第六节 地球物理资料数字处理(第二十、二十一讲)-第 8 页第七节 最小熵反褶积(Minimum Entropy Deconvolution) 地球物理技术具有高科技密集，更新速度快的特点。每一次新技术的出现及其应用都对数字处理技术产生巨大的冲击。尤其是随着现代信息处理技术的发展和计算机技术的进步，更加促进了新技术在地震资料数字处理中的应用。最小熵反褶积就是将信息论中熵的概念引入到反褶积中的技术。与最小平方反褶积不同的是，它不需要假定反射系数是白噪，子波为最小相位的条件，从这个意义上讲，它比最小平方反褶积方法先进，从另外的角度看，它的算法较为复杂，这又是这种方法的缺点。这里讨论最小熵反褶积技术，意义不在于其方法本身。目的是把反褶积技术引入更广阔的领域。 最小熵反褶积的目的是压缩地震子波，提高地震记录的分辨率。本节要点：l 熵的基本概念；l 最小熵反褶积原理；l 最小熵反褶积原理的步骤。一、 熵的基本概念(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）又可写成下式对上式整理，又可写成下式令 将代入（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）并求解，可得一新的褶积算子 5)重复(2)(4)步骤,直至得到满意的输出为止,这时得到的可做为反褶积算子。上述的迭代过程一般46次就可以结束。由于所解方程为非线性，这时输出的反褶积算子并不能保证最大方差范数V是唯一的极大值，但此时的反褶积算子仍不失其使用意义。在迭代的过程中为确保最小熵准则的一致性，必须对每次迭代所得的反褶积算子进行归一化处理，以使每次计算的最大方差范数V是使用统一的数量标准比较。在建立最小熵反褶积准则时，我们仅考虑输出具有简单外形的条件，并没有考虑输出脉冲的极性和其时间延迟，这可能导致输出脉冲的极性与实际反射系数的极性相反，或与实际的时间延迟不一致的情况出现，这是最小熵反褶积方法的一个缺点，只要我们合理选择初始的反褶积算子，我们可以避开这个问题。一般，若选线性反褶积算子的初始值为一尖脉冲时，其间脉冲出现的位置应接近实际位置。图（5-13）是利用人工合成记录进行最小熵反褶积的算例。图（5-13）（a）是原始模型， （a） (b) (c) 图（5-13） 最小熵反褶积算例上部是输入子波，下部是反射系数序列。图(5-13)(b)上部是由所得反褶积算子计算的最小平方逆，下部是输入子波与反射数褶积所得的纪录道集，图(5-13)(c)是最小熵反褶积的输出结果，上部是所求的线性反褶积算子。值得一提的是反褶积算子的最小平方逆与模型的输入子波形非常相似，这说明，最小熵反褶积方法的效果令人满意。第八节 极大似然反褶积(Mostlikelihood Deconvolution)极大似然反褶积是建立在概率论中极大似然估计理论基础之上的一种反褶积新技术。由于反褶积准则是建立在概率基础之上，不需假设子波的最小相位特征，只需假设反射系数满足高斯-贝努力序列条件，而该条件符合实际情况，因此，极大似然反褶积方法适应性强。且所得每一个反射脉冲都比较准确地对应于反射界面，地震剖面的精度和分辨率均得到提高。由于极大似然反褶积算法比较复杂，这里只介绍基本原理。本节要点：l 极大似然反褶积原理；l 分块最优化迭代法。 一、极大似然反褶积原理 根据地震记录褶积模型，对于任意一地震道记录，可以写成下式现对于波和反射系数做如下假设 1对于任意相位的地震子波和反射系数，均可由ARMA（自由归平均）模型表示，即子波的Z变换为 （5-5-1）其中，m,n为子波参数 2假设，地层的反射系数序列满足高斯贝努力序列，且可用乘积模型表示，则 （5-5-2）其中是与地层散射序列有关的背景散射序列，是零均值的高斯白噪，它表示反射系数幅度的大小，是序列，它表示反射是否存在。，表示有反射，表示无反射，且其概率密度满足 （5-5-3） 其中，称为反射频次 对于地震记录的褶积模型，若江地震记录写到向量形式，则有 (5-5-4)其中 表示列向量，N表示记录长度。对于(5-5-4)式，可以进一步表示成下式(5-5-5) 为使表达式简洁，对所有的变量均以列向量的形式表示，则有其中，表示统计参数列向量，表示方差，表示变量的方差，表示的方差，表示的方差，为反射频次。设地震记录的似然函数为，根据概率理论，似然函数可以表示成下式(5-5-6) 其中,P表示概率密度，由条件概率的乘法定理，有(5-5-7)若设与无关，则(5-5-7)式可写为(5-5-8) 现在计算上式中的似然函数。根据概率分布性质，有其中 (5-5-10) (5-5-11) (5-5-12)上面各式中括号外的表示转置将(5-5-9)式(5-5-1)式代入(5-5-8)式，则似然函数的表达式为(5-5-13)由于似然函数与似然函数的自然对数等价，则有 （5-5-14） 图5-14 分块最优化迭代法示意图 上式表明，若求得的子波参数，反射系数参数，以及统计参数和与实际的地震记录的参数一致，则似然函数会出现极大值。可见极大似然反褶积准则很好地解决了子波压缩问题和反射脉冲出现的时间问题。进一步分析表明，极大似然反褶积可以同时获得子波参数和反射系数序列。此时，反褶积问题归结为，如何计算我们所关心的子波参数反射系数序列以及其他参数，使得似然函数达到极大值。分析(5-5-14)式知，计算中的似然函数需求解个未知参数，当地震数据很大时，直接求解几乎不可能。分块最优化迭代可以较好地解决上述问题，下面仅介绍最优化迭代的基本原理。由于采用分块最优化迭代法，至少可以保证所求出的子波参数和反射系数参数使似然函数收敛到某一个极值，也就是说，分块最优化迭代至少可以达到局部最优，不保证达到全局最优，这并不是极大似然反褶积方法本身的缺陷，但这也不是说一定不达到全局最优，只要我们合理选择初始参数，可以获得非常好的效果。 （a） （b） 图（5-15） 极大似然反褶积实例 图（5-15）（a）是输入的实际地震资料，选自我国的辽河油田，（b）是极大似然反褶积结果，可见，极大似然反褶积方法的精度是相当高的。