第八章 地震定位及地球速度反演3.docVIP

第八章 走时数据的反演 89 8.1 球层介质的速度反演 90 8.1.1 拐点法（古登堡法） 90 8.1.2 Hergloz – Wiechert方法 91 8.2 单台地震定位及Inglada法定位 100 8.2.1单台定位法 100 8.2.2三站求取震中位置 105 8.2.3 采用三台站求解地震位置的计算机算法 106 8.3 地震定位结果与台站分布的关系 109 8.4 迭代定位方法 113 8.4.1 迭代定位算法的基本思路 113 8.4.2迭代方法的具体实现 114 8.5 相对定位方法 131 第八章 走时数据的反演 在前两章我们根据已知的速度结构对射线追踪和走时曲线的计算问题进行探讨，推导了波速只随深度变化的一维速度模型的射线追踪的表达式。三维结构的射线追踪虽然较复杂，但它遵循的原则是类似的。现在我们来研究根据观测到时数据得到速度结构和地震位置的问题。前面两章是已知速度结构和地震位置，如何计算地震波走时的问题，而本章是根据地震波走时得到地震结构和震源位置的问题，是前面两章的反问题。在这个问题的研究中，地震学家通常把问题进行简化：第一种简化是对于震源位置已经清楚的地震，根据地震波走时数据求解速度结构，第二种简化是已经得出了该地区的速度结构，根据地震波到时求解地震位置。限于篇幅，对于第一种简化，本书仅讲授球层介质一维速度结构的反演。对于第二种简化，我们讲授地震定位的基本方法。 8.1 球层介质的速度反演 8.1.1 拐点法（古登堡法） 这个方法是利用走时曲线的拐点处与从震源处水平方向射出的射线相对应，据此可求出震源处的波速。 根据第七章球层介质中的Snell定律有： 其中h代表震源处的参数，都是常数，从震源向不同方向射出的射线，值不同，相应的射线参数p也不同，当=90(时，达到极大值，即射线参数也达到极大值，因此有 根据球层介质中的Bendorff定律（7-5-3）式，，有 而对应于走时曲线的拐点。因此走时曲线的拐点与从震源处水平射出的射线相对应。 这样，我们只要求得某地震的震源深度及有观察分析归纳得到其走时曲线，找出走时曲线的拐点，并确定该点的走时曲线的斜率，则 而对应的M点有，因此有 其中，h为震源深度，为走时曲线拐点M点的视速度（参看第七章球层介质的Benndorf定律），R为地球半径。 此方法原理清楚、方法简单、计算方便。但由于地震发生的深度大约在0-700km的范围内，采用这种方法只能求出0-700km处的波速。由于拐点不易找准，得到的精度较差，它又要求对应每个深度的地震就要总结出一条走时曲线，因而资料分析工作量很大。 8.1.2 Hergloz – Wiechert方法 考虑球面介质，速度随半径r 变化[v( r )]，根据7.2节中得出震中距的表达式有 其中=，为射线转折点距地心的距离，p=，R为地球半径。改变积分变数，则 ， 其中。 在此介绍一积分式， ， 其中，为转折半径为的射线参数，所以式表示从转折半径为0至到转折半径为的所有射线的积分。将式作用于震中距表达式式的两边，得到 方程的左边变为，考虑到,可以用分部积分的形式给出 式的第一项代入积分上下限，由于cosh-1(1)=0，且在p=，半径为R，Δ=0，因此第一项消失。仅剩下第二项，可写成： 的右侧为一双重积分式，改变积分变数，得到： 这里采用了不定积分公式： 整理变为 与式结合，得出 令，可以得到： 这就是射线最深点的半径的表达式。此表达式可以根据观测走时曲线t()得出。如果t()为一平滑曲线，我们可由其切线斜率得到p()