c语言地震子波6,地震子波反褶积程序

地震子波反褶积

盲解卷积算法

1. 原理：

几个重要概念：

1.1 褶积模型假设：

假设1：地层是由具有常速的水平层组成；

假设2：震源产生一个平面压缩波(P波)，法向入射到层边界上，在这种情况下，不产生剪切波(S波)；

假设3：震源波形在地下传播过程中不变，即它是稳定的； 假设4：噪音成分是零；

假设5：震源波形是已知的；

假设6：反射系数序列是一个随机过程。这意味着地震记录具有地； 震子波的特征，即它们的自相关和振幅谱是相似的；

假设7：地震子波是最小相位的，因此，它有一个最小相位的逆。

1.2 反滤波

如果定义滤波算子为f (t),则f (t) 与已知地震记录x(t)的褶积得到一个对地层脉冲响应e(t)的估计

e(t) = f (t) x(t)； (1) x(t) = w(t)* f (t)* x(t)； (2) δ(t) = w(t)* f (t)； (3)

1f(t) (t)* (4) w(t)

用流程图表示为：