第四章射线理论.docVIP

第四章 射线理论：走时 地震的射线理论与光学的射线理论相似，在过去100年中被广泛地用于解释地震资料。由于它简明和适合于广阔的领域，因此直到今天仍被广泛的应用。这些应用包括多数的地震定位算法、体波震源机制的确定及壳幔速度结构的反演。射线理论直观、易懂，程序的编制简单、有效。与更完整的解法比较，射线理论较直截了当地给出三维的速度模型。然而，射线理论也有一些重要的极限。它是高频的近似，对长周期或陡的速度梯度的介质，则有欠缺。它不容易对“非几何”效应，例如首波和绕射波做出预测。由于在层里多次的反射，因此必须对射线的几何形状作详细的说明，否则就难以研究混响和共振效应。 在这一章，我们只涉及地震波到时，而把对振幅和其他细节的考虑推迟到后面。然而尽管如此，这些内容对解决许多问题仍然是很用的。现代地震学的部分重要研究只用了走时资料。许多射线理论的理论基础来自于短时距方程（参见附录3）。然而，因为多数研究不需要这些结果，因此这里没有进行叙述。 波动方程向射线理论的过度 首先阐明在什么条件下波动方程能够过渡到几何地震学，即几何地震学在何种条件下能反映真实波动在空间传播的情况。 在第三章我们已经得到，如果在一个波长范围内的相对变化相对很小，即可以看成常数时，等价的P波和S波方程可以表示为： 我们统一写为： 取为纵波方程，为横波方程。一般情况下，该方程的解为： 这里讨论一般的谐和体波，式中r为坐标原点（或震源）至波面的法向距离，为振幅，而则表示波动在空间传播的相位。 讨论波前的形状及其运动情况即可描述地震波的传播过程。因为波前是一个运动着的曲面，在该面上每一点于相同时刻处在同一个相位。波前方程可以由来确定。令 =0 即是的初始值。上式表明随着时间t的增加，r/c也相应增加。波动随时间向外传播。当波在介质中传播时，有一系列时间满足上式，即具有零相位的等相位面在空间中随着波的传播而连续分布，数目是无限的，t也是空间的函数： 称为特征函数，它确定波沿射线的走时，所谓射线就是处处与等位面垂直的线。于是波动方程可表述为： 则 可得： 根据实部和实部相等，虚部和虚部相等可得： 我们讨论第一个方程： 得到： 当较小，而很小时， 即 由于 而为两个向量相乘的模，故有， 这就是特征方程或哈密顿方程，又称时间场方程、程函方程。它是一个具有纯粹几何形像的波面方程，通过它，波动地震学就过渡为几何地震学了。 上式具有重要的物理意义，如果介质的波速参数c已知，利用边界条件或初始条件，就可以求得时间场，从而可知任何时刻波前的空间位置，也就求得了地震波传播的全部情况，而用不着求波动方程的解。因此上式是几何地震学中最基本的公式。 但我们要记住，从波动地震学过渡到几何地震学的两个根本条件：（1）波长趋近于零，即只对高频适用；（2）不能趋于无限，对于球面波而言，，当r趋于零时，趋于无穷，即球心或聚焦点处几何地震学不适用。 上式还可以表示成向量形式： 其中为沿波传播方向的单位向量，从上式计算沿任意方向的线积分，得到从点S1到S2的波的传播时间为： 这是因为沿梯度的方向其值总是最小的。因此，上式告诉我们，两点之间波沿射线传播的时间最短，即 为变分符号，这就是著名的费马原理，它说明沿射线传播的时间与沿其他路径的时间相比为一个极值。 4.1 斯奈尔定律 考虑一个以均匀的速度在介质里传播，并与水平地表面相交的平面波（图4.1）。沿射线路径，波阵面在之间的距离为。射线相对于垂直方向的角度叫做入射角。把与波阵面在地表面分开的距离联系起来： 图4.1 在水平的表面平面波的入射。射线相对于垂直方向的角度叫做入射角。 （4.1） 因为，故有： （4.2） 或 ： （4.3） 这里u是慢度（，是速度），p叫做射线参数。注意，如果界面是地表自由表面，我们通过计算波阵面到达两个不同台站的时间，就可以直接测量p。射线参数p表示波阵面在水平方向的视慢度，这正是为什么有时把p叫做射线的水平慢度。 现在考虑一个下行波到达了一个上、下层分别为均匀层的界面的情况及其在下一层产生的透射波（图4.2）。如果沿射线以均匀的时间间隔画出波阵面，那么在不同的层里，射线就被分成不同的间距，我们可以看到，为了保持在分界面两边波阵面同步，在分界面，射线的角度必然发生变化。 图4.2 两个均匀半空间水平界面两边的平面波。 下一层较高的速度使波阵面间距更大。 在上一层有较低的速度（）和相应的较大的慢度（）的情况下。可以根据每一层里的慢度和射线相对于垂直方向的角度来表示射线参数： （4.4） 注意，这只不过是几何光学的斯奈尔定律在地震学上的说法。