复杂地表地震勘探中的高精度近地表速度建模方法.docVIP

复杂地表地震勘探中的高精度近地表速度建模方法 　　[摘 要]在勘察内蒙阿拉善某高放废物地质处置预选场址目标花岗岩体完整性与稳定性的地震勘探过程中，利用基于MSFM（Multistencils Fast Marching Methods）射线追踪算法的地震初至波走时反演方法进行了近地表高精度速度建模研究。获得了准确、精细的近地表速度模型成果。该成果不但可被用于解决对于陆上复杂地表地震勘探来说至关重要的层析静校正问题，而且可作为地震解释工作中重要的参考资料 [关键词]高放废物地质处置; 复杂地表; 近地表速度建模; 射线追踪; 速度反演 中图分类号：P5 文献标识码：A 文章编号：1009-914X（2016）23-0013-02 根据以往的探测经验可知，近地表速度不准确会严重影响地下地质结构的地震成像效果，对于地表条件复杂的工区更是如此。在勘察内蒙阿拉善某高放废物地质处置预选场址目标花岗岩体的过程中，为了解决对于陆上工区地震勘探至关重要的层析静校正问题，利用基于MSFM（Multistencils Fast Marching Methods， 多模板快速步进）射线追踪算法的地震初至波走时反演技术，建立了沿地震测线准确、精细的近地表速度模型。该成果不但被用于有效解决了后续的层析静校正问题，而且是最终地震地质解释工作中精确、可靠的参考资料 1 方法原理 复杂地表问题通常指包含剧烈起伏的崎岖地表和强速度变化的近地表速度结构。在选择适用于复杂地表初至波走时反演的射线追踪算法时，通常需要考虑计算精度、计算效率、以及算法对于模型复杂度的适应能力三个方面[1] 1.1 MSFM走时计算公式 MSFM（Multistencils Fast Marching Methods）多模板快速步进算法是由Hassouna等提出的对于经典FMM方法[2]的改进算法[3]。其原理是在标准FMM算法的基础上通过坐标旋转产生多个FMM计算模板，使对角方向邻点也能够参与计算。该算法相对标准FMM方法在计算精度方面有显著提高，并且更容易推广到处理高维高阶问题 其计算公式如下（为了简化讨论，仅列举横向和纵向空间步长相等时的情况）： ① 相邻模板（S1）： 对于方向导数的一阶近似，可通过下式求解： （1） 对于方向导数的二阶近似，可通过下式求解： （2） ② 对角模板（S2）： 对于方向导数的一阶近似，可通过下式求解： （3） 对于方向导数的二阶近似，可通过下式求解： （4） 在上述公式（1）～公式（4）中，为待求结果（也就是理论计算走时），为地下介质空间慢度的函数，为空间步长，为邻域范围内的最小到达时 1.2 算法实现过程 算法的实现需要利用Dijkstra算法思想[4]，即首先将计算网格点划分成三个子集，分别为到达时已确定的{已知点子集（known）}，已经被计算但最终到达时待定的{窄带子集（narrow_band）}和尚未被计算的{远点子集（far）}： （1） 首先将炮点初始化为{已知点子集}，到达时为零;并根据MSFM模板计算炮点所有在{远点子集}中的邻点的到达时，并将这些点移入{窄带子集}. （2） 选取{窄带子集}中走时最小的网格点移入{已知点子集}，并利用模板计算该点所有在{窄带子集}和{远点子集}内的邻点的到达时;若邻点在{远点子集}中，则将该点移入{窄带子集};若在{窄带子集}中，则根据到达时大小选择更新或者保持原到达时（取较小者） （3） 重复实施步骤（2）直至所有网格点移入{已知点子集} 1.3 Model94模型测试 为了检验MSFM射线追踪算法的计算精度及其对于复杂模型的适应能力，本次研究首先设计了针对Model94复杂地表模型的测试试验。具体试验流程为：① 给定观测系统，利用时间2阶、空间20阶交错网格有限差分法[5]求解声波方程正演出一套地震记录。② 利用2阶MSFM算法计算在相同观测系统下各检波点位置的波前面到达时（初至波到时）。③ 将MSFM计算结果投影到有限差分法地震模拟记录上，对比通过MSFM方法计算得到的初至波到达时（走时）与模拟地震录中初至波的吻合程度 按照上述流程共进行了等间隔分布的277炮数据的数值测试试验。其中位于模型正中的第139号单炮数据数值测试结果如图1所示 可以看到，就Model94模型复杂度的地质问题来说MSFM算法的计算精度是可以满足要求的。图中显示的波前面到时计算结果与通过高阶交错网格有限差分法得到的模拟地震记录中的初至波基本吻合 1.4 初至波走时反演速度建模实验 相对于全波形反演（FWI）来说，射线类走时反演算法具有计算效率高、在应用中稳健可靠、