各向异性地球介质中地震波传播特震研究-北京大学教务部.doc

北京大学政学者论文集（2003年） 各向异性地球介质中地震波传播特震研究 PAGE  PAGE 331 各向异性地球介质中地震波传播特震研究 地球物理系00级 黄俊伟 摘要 通过方程的演绎，论文较完整地论述了一套研究地震波在分层各向异性半空间介质模型中的传播特性的理论方法。在理论上涵盖了地震波传播介质的一般情况。文章既保留了原方法的理论框架，又在此基础上做了些改进。已有的理论方法主要包括分层半空间介质的特性、各种各向异性介质模型下的地震波方程和均匀各向同性介质模型中地震波的传播特性。我们将这些理论成果结合起来，按均匀各向同性介质模型中地震波传播特性的理论框架重新论述地震波在分层、各向异性、半空间介质中的波动方程。以时间域中地震波的位移为结论。 Abstract In this article, through the deduction of principle functions, we derive a set of theory concerning seismic propagation in Multi-layered Anisotropic Half-space Media, which, theoretically, contains the general conditions about wave propagation. Preserving the original theoretical frame, we achieve some improvements. The original theories contain the characteristics of Multi-layered Half-space media, seismic functions in various anisotropic models, and the characteristics of wave propagation in homogeneous media. We combine all those results and re-deduct the wave propagation functions in Multi-layered Anisotropic Half-space media with the seismic wave displacement in time domain as our ultimate goal for seismogram synthesis. 一、引言 随着分析地震数据方法的发展，各向同性地球介质模型已逐渐不能满足研究的需要。而各向异性的特性在地球介质中几乎普遍存在，它对地震数据的影响必须考虑在内。其中一个最紧迫的要求是能对各向异性介质合成地震图。 Fryer 和 Frazer[6][7] 在1984年描述了波在纵向分层任意各向异性介质中传播的数学特性， 并在原则上给出完整的合成如何计算。Kennett [10]对分层介质的特性已有深入地研究。陈晓非教授[3]在1999年提出了对成层半空间均匀各向同性介质中传播的地震波的合成。各种理论各有所长，但随着对地震数据的增长及各向异性介质模型趋于一般化，计算的效率与各向异性特性必须考虑在内。 本文只在已有的理论工作的基础上推出成层半空间各向异性介质中的地震波波动方程，重点针对三维介质，为数值实现提供全部方程。 二、基本原理 本文采用图一所示的地球介质模型：分布于半空间的N层各向异性均匀介质，各层交 界面水平， 和 为第j层的上下界面。最顶层为自由表面；第j层的弹性参数由表示。一任意点源置于第s层，认为s(N+1)，即最底层无源。 自由表面 源 图一 1 根据弹性力学基本原理，分层介质中的波动方程为： 其中，(j,s)为介质层数编号 分别取 考虑到介质弹性常数的对称性，我们可以采用对应关系[14]：“11~1、22~2、33~3、44~4、12,21 ~6、13,31~5、23,32~4” 即可与6x6矩阵元素一一对应。于是得单斜介质的刚度矩阵： （3） 共13个参数，值由所在层介质物性决定。 对时间作Fourier 变换，在域中方程为： 相应的频率域中的边条件为： (i) (ii) 对 (iii) 三、二维情况 力沿着y轴方向延伸到正负无穷，因此相应的波场沿y轴为常量，即与y坐标无关，但依赖于x和z坐标。为此，如下选择基函数：