8第七章面波.docVIP

第七章 面波 至此，我们的论述局限于体波，即在全空间地震波动方程的解。然而，当介质中存在自由表面时，可能有其他的解，称之为面波。沿地球表面传播的面波 有两类：Rayleigh波和Love波。对横向均匀的模型，Rayleigh波径向偏振（），存在于任何自由表面，而Love波横向偏振，且要求速度随深度（或球面几何结构）增大。面波通常是远距离记录中最强的波至，对地球的浅层结构和震源的低频特性提供最好的约束。面波与体波有许多方面的差别——面波传播比较慢，振幅随距离的衰减通常小得多，速度与频率有很强的依赖关系。 7.1 Love波 7.1.1 Love波的形成条件 设有均匀弹性半空间，上面覆盖一弹性层，层厚为h，用这样的模型来简单描述地壳覆盖在上地幔的情况。取x,y在自由表面上（z=0），z轴垂直向下，令层中横波速度为，密度为，令半空间中横波速度为，密度为，且有，Love波在层中为v1。半空间中为v2，为简化分析，仍考虑平面波的情况，并令波沿x方向传播，由于考虑的是SH型面波。因此振动应垂直于x轴且平行于分界面，即振动应沿y方向。则V1V2应满足波动方程： 我们解上面的第一个方程：根据分离变量法，并且与x,z,t有关，我们可设，其中，代入上面的式子，可得： 该方程的特征方程为，其解为，根据微分方程理论，，其中A,B为常数。同理可以得到。因要满足的收敛条件，将中第二项去掉，上边方程的解可以写为： 其中 k为波数，c为面波速度。上边的式子应满足自由表面边界条件及层与半空间的层面的边界条件 因此有： 消去B得： 这是一个齐次方程组，要使得该方程有非零解，需使得行列式的系数为零，即： 其中，为双曲正切，其定义为，下面为x为实数的图像为： %tanh(x)函数 plot([-10:0.1:10],tanh([-10:0.1:10])) hold on plot([0 0],ylim,k); xlabel(x);ylabel(y); title(y=tanh(x)) 当x为实数时，由于和为正值，找不到函数值对应的情况。当和均为虚数时，方程右边为实数。而为虚数，也不满足方程。因此只有为虚数，而为实数时才能满足方程，即方程有非零解的条件为： （7.10） 上式规定了Love波在层里传播的频散曲线。注意，在与之间变化（没有过临界反射）。由于正切函数的周期性，对每个值，有多个值。最小的值定义了基阶模式，第二个最小的是第一个高阶模式，等等。方程（7.10）没有的解析解，必须用数值法确定的值（见练习7.1）。 简正振型 解上述频散方程得： 这就给出了与的对应关系。要使得方程有实根，应有：。对于满足该条件的值，有无穷多个与之对应，而每一个值对应于一种Love波。当n=0时成为基阶Love波（Fundamental mode），n=n称为n阶Love波(overtone)。 由于，由上式可知： 截止频率 对n=0的基阶Love波，当时，,因而，因此有波长或。当时，，此时没有截止频率。 对n阶Love波，当时，，频率的最小值为；波长的最大值为；这就是n阶Love波的最低频率和最长波长，我们称之为截止频率；当时，，。 位移分布 下面讨论简正振型Love波的位移分布，根据前面的式子，消去B得： 并且根据前面的方程，A=B,得： 所以有 令 根据频散方程得： 可见在覆盖层中振幅为，令其为零，则可得到： ，即，得到：，在此深度上，Love波的振幅为零。由于在0~之间，因此基阶Love波在覆盖层无节点，一阶Love波有一个节面，二阶Love波有二个节面。 下面我们模拟地壳厚度为30km的频散曲线和Love波振幅随深度的分布： H=30; %地壳的厚度为30km vs1=3.6;vs2=4.6; %地壳和地幔的速度，单位km/s %den1=2.9;den2=3.38; %地壳和地幔的密度，单位kg/m3 miu21=1.8; n=0; c=vs1:0.01:vs2; c2=c.*c; sqc1=sqrt(c2/vs1/vs1-1); sqc2=sqrt(1-c2/vs2/vs2); atann=atan(miu21.*sqc2./sqc1); omiga=c./(H*sqc1).*(atann+n*pi); figure(1) plot(2*pi./omiga,c) hold on omiga=c./(H*sqc1).*(atann+1*pi); plot(2*pi./omiga,c) omiga=c./(H*sqc1).*(atann+2*pi); plot(2*pi./omiga,c) xlabel(Period/s); ylabel(Velocity/km.s^-^1) figure(2) c=4.0; c2=