二维平层半空间波场中震相的分离和数值实现-北京大学教务部.doc

二维平层半空间波场中震相的分离和数值实现 地球与空间科学学院00级 尚学峰 摘 要 本文讨论的是二维的均匀的分层平层半空间介质模型。由广义射线的方法可以得到全波场的解，在此基础上，按照到时和路径的不同，可把全波场分解为各个独立震相。每个独立震相在时间域的贡献可以用Cagniard De-Hoop方法通过简单的数学变换操作后直接在时间域上作Laplace褶积即可得到，从而避免了复杂的Fourier反变换。理论上，Cagniard De-Hoop方法在处理二维问题时是精确和完全的，但是在数值上，远场和理论解误差较小而近场时个别点上误差较大。 关键词：平层半空间介质，广义射线，独立震相,Cagniard De-Hoop方法 导 言 由于地球的横向不均匀性一般远远小于纵向不均匀性，计算半空间层状介质的合成理论地震图，对于我们分析地震数据和认识地球内部结构以及震源动力学都有着重要的意义。上个世纪七十年代以来，很多研究工作完善了理论地震图的计算。其中离散波数法（Bouchon and Aki，1977；Bouchon 1979,1981））域上的全波场解 域上的全波场解的基本方程 考虑如下的平层半空间介质模型： 如图1-1所示，半空间中共有N层横向各向同性的均匀介质，各层间界面水平，和之间称为第j层，每一层的性质由、和描述。最上层为自由表面。任一点源（单力点源、双力偶源、爆炸源）位于第s层中。 第j层中任意一点在域上满足弹性动力学方程： ， (1-1-1) 当 相应的边界条件为： 1)自由表面条件 , 当 ; (1-1-1a) (2)边界上位移、应力连续条件 ; (1-1-1b) (3)无穷远边条件 . (1-1-1c) 其中， 为第j层任一水平面上的应力。 在二维问题中，我们引入一组完备的正交基： (1-1-2) 利用上式的完备性,将位移场和源项展开为： , (1-1-3) 和 . (1-1-4) 结合、和的正交归一性，得到： , (1-1-5) (1-1-6) (1-1-7) 引入应力函数： , for SH波 (1-1-8) for P-SV波 (1-1-9) 利用中间函数(1-1-8)和(1-1-9)，我们可以将方程(1-1-5)、 (1-1-6) 和(1-1-7)统一的写为： (1-1-10) 其中，SH波有： , (1-1-11a) , (1-1-11b) , (1-1-11c) P-SV波有： , (1-1-12a) , (1-1-12b) . (1-1-12c)相应的边界条件变为： (1)自由边条件 for SH 波, (1-1-13a) for P-SV 波; (1-1-13b) (2)连续性边条件 for SH 波, (1-1-14a) for P-SV 波; (1-1-14b) (3)无穷远边条件 for SH 波, (1-1-15a) for P-SV 波 . (1-1-15b) 所以SH波对应的方程和P-SV波对应的方程形式上是类似的。只要求出 (1-1-10)的通解，即可求出、和，进而求出在域上的总波场。 SH波和P-SV波的各类反透射系数 考虑方程(1-1-11),根据常微分方程的理论，形如(1-1-11)的非齐次的线性常微分方程组的通解可写为：