反褶积模型课件.pptVIP

* 反褶积模型 当地下为水平层状介质，地震波垂直入射 条件下，无噪声的地震记录表示为 x i = r i * b i = b i * r i 式中，r i 为反射系数序列，b i 为地震子波。 地震记录的形成可用滤波过程表示 r i x i b i 在滤波方程中, 已知x i , 未知r i , b i 从理论上讲, 由一个已知函数x i , 求解两个 未知函数r i 、b i , 这样的方程无解。但通过一定 的假设条件， 可以从滤波方程中近似估计地震 子波b i , 解出反射系数r i 。 反褶积概念 求解一个滤波器a i , 它作用于地震记录 x i 后的输出为反射系数r i , 这个过程恰好是 地震记录形成过程(褶积)的反过程, 所以称 之为反褶积。 a i 为反褶积因子。 r i b i x i x i r i a i 地震记录形成过程 反褶积过程 反褶积过程的输出： r i = x i * a i = r i *b i * a i 满足上述方程的条件是 b i * a i = ? i ?函数 ?函数的定义： 1, i = 0 ? i = 0, i ? 0 ?函数的频谱与Z变换： ?(f) = 1， ?(Z) = 1， ?函数的性质： y i * ? i = y i { 反滤波 当地震子波b i已知时，方程 b i * a i = ? i 的Z变换为 B(Z)A(Z) = ?(Z) = 1 则有 A(Z) = 1? B(Z) 反滤波 反滤波因子a i为子波b i的逆，b i 是有 限长的，准确的a i是无限长的，而我们只 能取有限长度。 设计有限长度反滤波因子a i时，我们 应用最小二乘法(最小平方)的原理。 关于反褶积的一个基本认识 从通过简化(无噪、子波已知)得到的反 滤波方程中，不能求解准确的反褶积因子。 任何一种反褶积方法都不能精确求解反 射系数，只能近似求解。 子波的自相关函数 设反射系数是白噪声，则地震记录的自 相关等于子波的自相关。 地震记录自相关的频谱 R x x(f ) = ?X(f )? = ?R(f )B(f ) ? = ?B(f )? = R b b(f ) 2 2 2 白噪声的振幅谱和相位谱 反射系数 振幅谱 相位谱 反射系数 振幅谱 相位谱 反射系数 振幅谱 相位谱 合成地震记录及其自相关(左) 子波及其自相关(右) 合成地震记录及其自相关(左) 子波及其自相关(右) * * * * *