psv波的展平变换及其1级修正数值解.docVIP

PAGE PAGE 130PSV波的展平变换及其一级修正数值解 地球物理系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)()1202.83.506.02102.93.656.33153.52.504.94254.03.806.55203.24.507.662.45.609.34.1对不同