水平层状海洋-地球模型中地震面波振型解.pdfVIP

中目村荸(D辑) 第3l卷第9期 SCIENCElNCHINA(SeriesD)2001年9月 水平层状海洋．地球模型中地震面波振型解4 陈蔚天 陈晓非” 北京大学地球物理学系．北京10087 摘要 提出了一种求解水平层状海洋一地球模型中面波振型问题的新算法这种算法 具有理论上的简洁和计算上高效特性以及在高频情况下数值计算上的稳定性，从而彻 底消除了在高频情况下数值计算上的精度失真问题，为研究海洋一地球模型中面渡传播 及相关问题提供7一个有力的工具 关键词 面波振型解 水平层状海洋-地球模型 地震面波 含有液体层的水平层状介质中的地震面波问题是地震学中的一个经典问题，对于研究地 球的结构、确定震源参数以及合成理论地震图都有重要应用【】。】．这方面最具代表性的工作是 计算任意多层液、固体层状介质中地震面波问题的Haskell矩阵方法f4．51以及传播矩阵方法[61 由于Haskell矩阵与传播矩阵方法在高频情况下存在数值计算的不稳定性问题，人们又提出了 一些改进方法07～191，但均未从根本上解决这一问题Kennett等人从完全不同的角度出发，采用 反射一透射系数，提出了一种数值计算上稳定的算法来求解层状介质中的面波振型问题…’ 陈晓非在Luco和Apsel的广义反射．透射系数方法f121基础上提出了一种系统、有效的算 法来求解固体层状介质中的面波振型问题…】．这种算法不仅具有理论E的简洁和计算上的高 效特性，而且彻底消除r高频情况下数值计算的不稳定性问题．然而由于只考虑了固体情况， 陈的算法不能用于计算海洋一地球模型中的面波振型问题本文将推广陈所提出的算法，使其 也适用于包含多层液体层的层状海洋。地球模型． 1基本方程及其通解 由于海洋中不同深度的温度、压力和盐度的不同，液体部分可分为不同层．固体部分根据 不同深度密度p和弹性常数A，“也划分为不同层本文所考虑的水平层状海洋一地球模型可 参照文献【13】中的图1，所不同的是第1至第工层为液体在本文中，我们将根据面波振型是齐 次弹性动力学方程(不含源项)在给定边界条件下的非零解这一基本原则”…，系统地导出水平 层状海洋一地球模型中面波振型解 在频率域上，均匀各向同性的多层海洋一地球模型中的齐次弹性动力学方程可表述如下． 在液体层中： 一P(』’甜2u【J’扛，∞)=Ⅳ川VVH‘7’(x，∞)，J=1，2，…．厶 (1) 在固体层中： (2) 一p(J’∞2Ⅱ(J’(工，∞)=(∥，’+2“(川)vv“‘，1(j，co)一“‘，’VxVxu‘，’(工，∞)， 2001．01．20收辅，2001一02—27收修改稿 t国家杰出青年科学基金(DY和国家重点基础规划发展研究(G1998040702)资助项目 ”通讯联系人 第9期 陈蔚天等：水平层状海洋一地球模型中地震面波振型解 713 J=L+1，L+2，…，N．其中P¨’，Ⅳ【川，孙。’分 别为第J层的密度和Lame常数，在液体层中 “‘J’为零，由于液体层中的sH波场恒为零， 而固体层中SH渡的解与纯固体层中的SH 波的解一致，因此本文将只考虑海洋一地球 模型中P．SV渡场的解． 为了有效地求解方程(1)和(2)，我们利 用横向波函数矢量基R：‘(r，0)，s?(r．0)和 r，(r，8)将波场展开”l，如此一来，方程(” 和(2)被转化为如下关于由展开系数函数 相速度#ran·s—o 图l周期lS时久期函数与波慢度倒数的关系以及 “R(z)，“s(z)，。R(z)和盯5(z)所定义的位 其根的分布 移．应力矢量的一阶常微分方程组如下．