常见的反褶积_学习基础.ppt

预测反褶积的原理和计算方法 主要内容 什么叫预测 预测的条件 预测滤波 预测反褶积 预测反褶积的计算 Omega系统的预测反褶积模块及其主要参数 什么叫预测? 预测就是根据过去和现在已发生的事实判定将来会出现的情况。 在数学上，对一个时间函数的预测是指该函数某一点的值用其前面若干个值的线性组合表示出来。这种预测称为线性预测。 预测的条件 并非所有事物都可线性预测。函数x(t)可线性预测的条件是: x(t)为平稳随机过程，即它的统计特征：数学期望 和方差 是与时间无关的量，且自相关函数rxx(τ)只与时差τ有关 。 我们认为地震记录满足以上条件，因而可做预测。 预测滤波 在地震勘探中，我们认为地震记录是平稳随机过程，因而可以预测。 根据地震记录褶积模型的假设，地震记录x(t)由地震子波b(t)和地层反射系数g(t)的褶积构成： 我们先假定b(t)为一物理可实现的最小相位信号， g(t)为白噪序列。在时刻(t+τ),地震记录的振幅值可表示为： 在右端第二项中，令j=s-τ,上式变为： 记 设b(t)的反信号为a(t)，有a(t)*x(t))=a(t)*b(t)*g(t)=δ(t)*g(t)=g(t) 因为b(t)为一物理可实现的最小相位信号，因此有：当t0时，a(t)=0 将 g(t) =a(t)*x(t)带入x’(t+τ),得： 令s=j+k,上式变为： 再令 得到： 上式为一褶积表达式，它说明：x’(t+τ)是c(s)对x(t)的过去和现在值的滤波结果，称它为x(t+τ)的预测值，c(s)称为预测滤波因子。实际值与预测值的差 e(t+τ)=x(t+τ)-x’ (t+τ) 称为预测误差。 τ 叫做预测间隙、预测步长或预测距离。 预测反褶积 因为： 所以预测误差为： 特别地，当τ=1时，有： e(t+1)=b(0)g(t+1) 该式表明，当预测距离等于1时，预测误差与反射系数只差一个常数因子，因而可视为反射系数。于是，只要在预测滤波中输出预测误差就达到预测反褶积的目的，这就是预测反褶积。 但通常不用τ=1这种理想情形，而是令τ为大于1的某个数。 当τ=1时，预测反褶积就是脉冲反褶积。 预测反褶积的计算 预测反褶积计算的关键是求得预测滤波因子c(s)。由于子波未知，不能用公式 直接计算。 可用最小平方法。 最小平方法的数学模如下： 输入信号： x(t) 设预测滤波因子： c(s)=[c(0),c(1),…,c(m)] 期望输出： x(t+τ) (τ0) 预测输出： 预测误差： ? 误差总能量： 预测反褶积处理模块 在Omega系统中，预测反褶积处理由以下三个模块完成，即： 1、预测反褶积谱分析（PRD_DCN_SPCTRL_ANL) ； 2、预测反褶积算子设计(PER_DCN_OPR_DESIGN)： 3、反褶积算子应用(DCN_OPR_APPLY)。 预测反褶积谱分析 功 能：计算给定时窗的自相关函数。该自相关函数用于构造计算预测滤波算 子的方程： 主要参数：1、确定时窗 的参数（起始时间、时窗长度）： 根据资料情况和处理目的确定。 2、自相关长度： 可根据算子长度确定：自相关半长度不能小于算子长度，因 为方程右端用到了自相关函数的值，所以自相关半长度至少 要等于算子长度加预测距离。 预测反褶积算子设计 功 能：计算预测反褶积算子。每个道的各个时窗都有自己的预测反褶积 算子。如果需要，你也可以先将谱分析输出的自相关函数按某种 方式（如炮集）进行叠加，然后设计统一的算子。 输 入：自相关函数文件 输 出：1、预测反褶积算子；2、估算子波=预测反褶积算子的反信号。 主要参数： 1、算子长度：这里指的是预测滤波算子的长度。 预测反褶积算子长度=预测滤