NMO.rar

1、动校正简介

       多次叠加技术是地震勘探发展中的一次革命，简单说就是利用动校正速度（通过拾取速度谱得到）将CMP道集（经过静校正、去噪等处理）拉平，然后叠加起来提高资料的信噪比。

      动校正是地震勘探中一个必不可少的处理步骤，动校正这个名词翻译自英文Normal Moveout，准确来讲应该称为正常时差校正。根据精度的不同、介质的不同等，动校正有不同的实现方式，这里通过最简单的双曲线动校正来说明动校正是如何实现的。

2、动校正理论

      对于一个CMP道集，动校正做的工作是：对指定偏移距与反射旅行时(x,t)的振幅值，根据动校正公式校正到垂直入射的反射旅行时t0。动校正公式（公式1）为：

很明显，对NMO道集中指定(x,to)，利用下式（公式2）可以计算得到CMP道集中对应的反射旅行时t：

     实际中动校正公式是如何实施的呢？这里给出一个简化的实现步骤：（1）首先将一个二维矩阵作为NMO后的道集，矩阵的所有元素赋值为0，需要注意的是NMO后的道集的反射旅行时为t0；（2）在NMO道集中(x,t0)处，利用公式(2)计算CMP道集中对应的t；（3）将反射旅行时t处的振幅值作为NMO道集(x,t0)处的振幅值；（4）循环进行（2-3）步。上述步骤不精确的地方在哪呢？就在第（3）步，按参考文献[1]中所说，应该通过内插得到(x,t0)处的振幅值。

3、编程实现

      实现动校正的输入文件两个：CMP道集和动校正速度文件。动校正速度可以通过拾取速度谱或其他速度分析方法得到。输出文件为：动校正后的CMP道集。

      这里通过一个模拟数据进行NMO，NMO的Fortran程序及对应的输入输出文件可以通过附件下载。这里简单说明一下程序中的关键：（1）需要计算动校正速度，由于我们这里的模拟数据时三层层状模型，可以通过均方根计算公式

得到两个反射界面对应的均方根速度，公式中 $\Delta \tau _{i}$ 为第i层的双程旅行时；（2）要进行NMO，动校正速度文件必须是规则的，也就是说必须内插得到其他采样点的均分根速度（用商业软件也是这样实现的，速度谱通过只拾取数个或数十个动校正速度，其他采样点的需要内插）；（3）下一步就可以利用动校正速度对输入的CMP道集进行NMO了。结果如下图所示：

       这里就可以思考一下动校正拉伸的来源了，我认为就是由于动校正速度是通过内插得到的（除了拾取的速度谱点）。

4、参考文献

[1] Lenardo Uieda. Step-by-step NMO correction[J]. The Leading Edge, 2017.

[2] Yilmaz O. Seismic data analysis: Processing, inversion, and interpretation of seismic data.SEG, 2001.

5、注意事项

     附件中的输入输出文件是带道头的segy格式。