PSV波的展平变换和其一级修正数值解.docVIP

PSV波的展平变换及其一级修正数值解 地球物理系2000级 黄慧 摘要 我们考虑一个球状分层忽略中理想的地球模型，通过展平变换，将球状介质转换成平层介质，以便快速有效的利用现有的程序算出球状介质下PSV波的本征位移。同时，对展平变换的结果进行一级修正，并和没有修正过的结果作比较。 介绍 展平变换简化了球状介质中波的传播问题，它最早应用于地球周围空间中无线电波的传播(Price 1953; Koo and Katzin 1960)。当地球介质弯曲的影响变得重要时，分析弹性波的传播也遇到了类似需要展平的问题。对于SH波，已经有精确的展平变换公式(Biswas and Knopoff 1970)；而对于PSV波，Biswas(1972)和Bhttacharya(1997)分别给出了不同的一阶近似公式。本文直接从球状介质和平层介质的方程出发，推导了一个一般意义上的展平变换公式。 在计算PSV波的本征位移时，Biswas(1972)和Bhttacharya(1997)直接忽略了一阶小量。本文对计算结果进行一级修正，并和以前未进行修正过的结果作比较。 基本方程 2.1平层弹性介质运动方程（位移表示形式） 弹性介质中地震波的运动方程如下： (1) ,,分别为拉梅参数和介质密度；为位移，即位移在i方向上的分量。对于Rayleigh波，在直角坐标系下可以设定： 则方程(1)化为 (2) 2.2球状介质的弹性动力学方程 在球坐标系下弹性介质的运动方程为： …... (3) 其中为体积膨胀因子，为应变分量，为三个位移分量)，为拉梅系数，为密度。以上方程忽略了重力项。用球谐函数来展开位移： (4) 其中为n阶球谐函数，是振动的频率。将（4）代入（3），有： (5a) 其中，N=n(n+1)，应用琼斯（Jeans）公式：n+1/2=ka，有N=。注意该公式是在“自由振荡高频近似下”成立。不考虑重力因素下的边界条件为： (5b) 2.3用展平变换公式对球状介质的方程进行展平变换 对球状介质下的弹性波运动方程进行改造，采用如下公式（其中m为常数，a是地球半径）： (6) 将公式（6）代入方程（5），同时令 (7) 得到： (8) 注意到在忽略了于1/a成正比的项以后，方程（2）于方程（8）在形式上一样。有对应关系： 称公式（6）与公式（7）为展平变换公式。两个公式意味着两次变换：一是用公式（7）把球状介质的弹性参数转化为平层介质参数；二是在计算出平层介质的位移后，用（6）式变换得到球状介质的位移。Biswas(1976)在其文章中把m取为1；Bhattacharya(1997)则取m=3。 求解方程 3.1对方程进行降阶 对于半空间分层均匀的介质，考虑它的第j层，运动方程如下： (9) 其中 取中间函数 (10) 可以证明，，正好对应第j层的应力场。联立(9),(10)两个方程，可以得到： (11) 其中 现在问题就转化为求解方程（11）这个一阶常微分非齐次矢量方程。 3.2方程的0级近似解 由于a=6371，所以在方程（11）中，相对于其他各项，是小量，忽略，得到0级近似方程，如下： (12) 这个方程现在已经很成熟和快速的程序计算其数值解(Chen,1993)。 3.3方程的一级修正解 将已求得的0级近似解代入，求解方程。 (13) 在形式上有 (14) 其中 注意到在数值计算中，在z方向上是离散化的，由一阶常微分方程的线性性质，可以分别求对应一个离散点的，再对所有的求和。为了保证积分效果一致，单个离散点对应的方程为： 其中为该离散点的深度，为两个离散点之间的距离。 计算结果与讨论 取平层介质模型如下： Thickness of each layer(km) () 1 20 2.8 3.50 6.0 2 10 2.9 3.65