KH方法合成理论地震图(最新修改)0409.docVIP

横向不均匀介质中Kirchhoff-Helmholtz积分合成地震图 李红光1， 吴庆举1 1 中国地震局地球物理研究所，北京 100081 摘 要 本文将Kirchhoff-Helmholtz积分方法应用于横向非均匀介质的合成地震图计算。与反射率方法及二维有限差分方法进行了比较。表明KH积分方法能很好的模拟反射波震相，且精度较高。KH积分方法能够计算复杂界面的反射波，且其计算效率明显大于有限差分。 关键词 合成地震图，KH方法，反射率方法，有限差分方法 1 引言 合成地震图是地震波形解释的一个重要工具。自Helmberger提出了广义射线合成地震图计算方法，并将其应用于实际地震波形资料的对比和解释工作后，水平层状介质合成地震图的计算方法得到了迅猛发展，产生了广义射线、反射率法、WKBJ近似及全波理论等合成地震图计算方法。对于区域不均匀性介质，发展了计算横向不均匀介质中地震波响应的方法，大致可以分为三类[1]：直接求解弹性动力学波动方程及边界条件的数值方法，例如有限差分方法[2]，有限元方法[3]；将非均匀介质分解成参考介质和微扰介质，把参考介质中的无微扰场跟散射场相加得到非均匀介质中的总波场的扰动法[4]；高频渐近似方法。各种方法都有其一定的局限性：有限差分方法虽然能很好的模拟复杂介质地震波传播过程，但由于有限差分方法非常耗时，限制了其在地壳、上地幔尺度的地震波波形模拟和解释工作中的广泛应用。扰动法只适用于横向非均匀性相对较弱条件下的波形模拟工作。因此，目前复杂介质中波形模拟应用较为广泛的是高频渐近似方法。随着地震台网的密集，特别是现代大型地震科学台阵观测技术的问世，使地壳上地幔横向不均匀性的调查和研究成为可能。所以我们迫切需要一种能够快速计算横向不均匀介质的合成地震图方法。 当介质的非均匀性比地震波波长大得多的时候，高频渐近方法效果最佳。最典型的高频近似方法是传统的射线理论方法，但当存在奇点的时候，传统射线理论方法不能正常工作。为了克服奇点带来的困难，人们陆续发展了Kirchhoff-Helmholtz积分方法[5]、高斯光束方法[6]、Maslov方法[7]、WKBJ[8]等近似方法。KH积分方法是一种非常好的计算横向非均匀介质合成地震图的方法，KH积分是把反射界面上每个点都看作一个点源，认为每个点源对地震波振幅都有一定的贡献，把反射界面上每个点的贡献相加就得到了地震波响应。KH方法不仅计算效率高，而且还能计算介质中存在焦点和散射点的地震波响应。 KH积分方法最早被应用在声波传播中[9,10]。Hilterman、Trorey和Berryhill用KH技术模拟了水平分层介质近垂直反射。Hilterman＆Larsen，Berryhill，Carter＆Frazer，Deregowski＆Brown分别对此方法进行了拓展,并把它应用在速度横向变化的情况。 2 KH积分方法的原理 如图1所示，V是一个2D或3D弹性介质体，令为V的边界，是V边界上任一点的法线向量。是作用于空间x1处的体力，且在V和上为零。分别是体力产生的位移和应力张量，它们满足频率域弹性动力学方程。同样，力是作用在空间x2处的另一体力，分别是体力产生位移和应力张量，它们也满足方程。 根据Gauss散度定理得到： （1） 图1 一个带有边界的体积单元V，是边界的外法线方向， 我们希望计算出力在点处引起的位移。 Figure.1 A volume V in(E2 or E3) with boundary .The vector is the outward-pointing unit normal to .We wish to calculate the displacement at x2 due to a force f1 at x1. 因为对于任意向量u和二阶张量有：。公式（1）中等号右边部分就变为： （2） 由于都是与四阶弹性张量c相关的，而弹性张量c又是对称的，所以：。那么积分公式（2）的最后两项就互相抵消了。公式（2）剩余部分，我们分别用、代替、。因为在V和上为零，积分公式（2）就化简为：。根据以上的化简我们最终得出：如果在V和上为零，那么有： （3） 我们可以根据实际需要来选择力的表达式。如果接收点处的检波器测量的是沿方向的位移，那么我们选择，从而： 公式（1）简化为：