地震资料处理中的反褶积处理PPT.ppt

地震资料处理中的反褶积处理;地震资料处理中的反褶积处理;反褶积概述;预备知识;信号;频谱的计算;信号的离散化;奈魁斯特频率; 频率折叠示意图;Z变换;信号的相位特征;相关分析;相关函数的定义;相关与褶积的关系;相关函数的频谱;两种特殊信号;反信号;预测反褶积的基本原理和计算方法;脉冲反褶积;1、脉冲反褶积的假设条件;2、脉冲反褶积的原理; 在反射系数函数是白噪声的前提下，有： rbb(t)=rxx(t)； 在地震子波b(t)为最小相位物理可实现信号时，有：当t0时，a(t)=0。 于是上面的方程变成为： 再将方程两端同除以b(0),则有：; 该方程可以求解，所得的解与反子波算子a(t)只差一个常数{b(0)}倍。;3、脉冲反褶积的计算;预测反褶积的原理和计算方法;什么叫预测;预测的条件;预测滤波; 令s=j+k,上式变为： 再令 得到： 上式为一褶积表达式，它说明：x’(t+τ)是c(s)对x(t)的过去和现在值的滤波结果，称它为x(t+τ)的预测值，c(s)称为预测滤波因子。实际值与预测值的差 e(t+τ)=x(t+τ)-x’ (t+τ) 称为预测误差。 τ 叫做预测间隙、预测步长或预测距离。 ;预测反褶积;预测反褶积的计算;选取c(s)，使Q达到最大。为此令 或 令 于是有： ;将以上方程写成矩阵形式就是： 以上方程的系数矩阵和左端的向量均由 x(t) 的自相关函数构成。该方程叫做预测方程，求解此方程，即得到最小平方意义下的预测滤波因子c(s)，用c(s) 对x(t) 滤波，若输出x’(t+τ),就是预测滤波，若输出e(t+τ),就是预测反滤波或预测反褶积。 因为 所以反预测积算子为 ： ;预测反褶积的几个主要参数;（1）算子长度;（2）自相关长度;（3）白化因子; 在实际应用中，并非仅仅是为了使方程有唯一确定解。因为数学上的解有时并不适合实际物理问题的要求。以脉冲反褶积为例： 在脉冲反褶积中，反褶积算子a(t)是地震子波b(t)的反信号： b(t)*a(t)=δ(t) 在频率域就是： 或 显然要使上式成立，对任何频率f,必须有B(f) ≠0，此外， B(f) 也不太接近0，否则会使A(f)的值在这一频率点上过大。因此预先白噪化在频率域就是将B(f) 加上一个小 数w，使其不那么接近0，这时有： W就是白噪因子。 ; 在反褶积处理中，大多都需要预先白噪化这一步骤，但它不是反褶积理论推导中的必然步骤，而是根据我们的需要人为地加上去的。因此白噪因子只能是一个很小的数，即只能对理论公式做少许修改。过大的白噪因子可能把理论公式改得面目全非。 在实际处理中，白噪因子参数应根据资料的具体情况由试验确定。;（4）预测距离;Omega系统的主要反褶积模块;预测反褶积处理模块 ;预测反褶积谱分析;预测反褶积算子设计;反褶积算子应用;GAP=20MS;地表一致性反褶积;基本原理和计算方法; 将各个频谱写成指数形式： W(f)= Aweiφw S(f)= Aseiφs R(f)=Areiφr G(f)= Ageiφg M(f)= Ameiφm 其中，字母A表示振幅谱，字母φ表示相位谱。 这样一来，子波复频谱的振幅谱和相位谱可写为： Aw=AsArAgAm (振幅谱) Φw=φs+φr+φg+φm (相位谱) 对振幅谱两端取对数，得： LnAw=LnAs+LnAr+LnAg+LnAm （对数振幅谱） 当炮点和检波点变化时，对数振幅谱方程将变成一系列方程。方程个数一般大于未知数个数，但其中独立的方程个数一般小于未知数个数，因此该方程组无确定解。为了求解该方程组，用最小平方法，即要求输入谱与分解谱的误差能量达到最小。于是可得出四个分量的迭代计算公式，再用Gauss-Seidel 法计算出子波对数振幅谱的炮点、检波点、共中心点和共偏移距分量。 ;高斯-赛德尔迭代法 为书写方便，将对数振幅谱表达式中的符号稍加改变。 记：Aw=A, A