地震层析算法网格尺度适应性改进.doc

地震层析算法网格尺度适应性改进 张小磊，张东 (武汉大学 物理科学与技术学院，武汉430072 摘 要：常规地震走时层析成像算法由于很少考虑网格尺度对反演结果的影响，容易导致计算结果受到网格尺度很大的制约，即只有在某一个特定的网格尺度之下，算法的计算结果才是比较可靠的。针对这一问题，本文提出了通过改进的插值算法、结合层内与层间处理的平滑算法、以及计算区域边缘递推算法，以提高算法的网格尺度适应性，使不同尺度之下的反演结果，都能够与实际情况相符，并使其变化规律保持一致。实验结果表明，不仅提高了层析成像对网格尺度的的适应性，也提高了精度，更加符合实际地质结构。 关键词尺度插值平滑递推 分类： 631.4 文献标志码： The improvement of the adaptability to grid scale for seismic travel-time tomography Zhang Xiaolei, Zhang Dong (Wuhan University, School of Physics and Technology, Hubei Wuhan 430072, China) Abstract:The effect of grid scale on the inversion results was rarely considered in conventional seismic traveltime inversion algorithm, resulting in the calculation results was seriously restricted by grid scale, and the inversion results of conventional algorithm were more reliable only in a certain scale. In order to solve this problem, this paper proposed an improved method to improve the grid-scale adaptability of the algorithm and make the inversion results of different scales match with the actual situation, also keep the change rule consistent by the improving interpolation and smoothing algorithm and using recursive algorithm to correct the speed value of edge area. The numerical results showed that this method would not only improve the grid-scale adaptability of the tomography results, but also improve the accuracy of the inversion, in addition make the results more in line with the actual geological structure. Keywords:grid-scale adaptability; interpolation; smoothing; recursive 0 引言 地震走时反演在速度成像、风化层速度与层厚探测以及层析静校正中都起着重要的作用，已广泛应用于地震勘探资料处理与地质结构调查[。4-6]与实际计算结果分析都表明，常规地震走时层析成像算法，其离散网格单元的尺度选择为一个值一般情况下，网格尺度根据、的以及要求的层析深度来确定4,6]，选择合适的网格尺度可以得出较好的层析结果。但由于条件的限制，在实际中，选择网格尺度进行层析因此解决传统算法对网格尺度依赖的问题。 尺度arametrization）问题，6,7]。2003年Hu-Wei Zhou等人在走时层析，网格尺度[8-10]，然后按照设定的权值取网格尺度计算结果的均值作为最终结果方法，。E. Kissling[4]等人提出三网格技术，认为以较小的网格尺度进行射线追踪，以最适合的尺度描述地下介质的速度结构，以较大的尺度进行反演计算是最佳的模型参数化方法。2002年赵连峰等人提出了交错网格技术[1]，正演采用小网格尺度，反演采用，；2005年马德堂等人[12]， 针对层析结果对网格尺度依赖性[13]，作者通过算法研究与模型实验研究发现，对网格尺度适应性影响较大的因