一种新的线性地震定位方法.docVIP

一种新的线性单事件地震定位方法 地球物理系98级 田玥 摘要 本文介绍了一种完全线性的单事件地震定位方法。由各台站P波和S波的到时差与台站和震源距离的关系得到一组线性方程，采用奇异值分解来求解方程组，能较好地确定震中；利用台站与震源的距离可事先求出的特点，分别对震源深度和发震时刻单独求解，二者完全解耦，结果较理想。在选取方程组时，每个台站与其它所有台站配对相减，做法最简单而且结果精度最高；如果对到时读数进行扰动作为虚拟台站的到时，可进一步提高定位精度。本文分别讨论了到时读数误差、台站平面分布和台站深度对定位精度的影响，并且定义了定量描述台站平面分布的无量纲量q值。与经典的Geiger方法相比，本文方法不需要对走时求偏导数，不需要迭代，并且解决了无法得到震源深度有效信息的问题，但适用范围较狭小。 Abstract In this article we introduce a new linear method for single event seismic location. With the relationship between the differential travel time of direct P and S waves of a station and the hypocenter distance, we derive a set of linear equations. By using the SVD technique, those linear equations can be solved effectively. Moreover, we investigate the effects of data noise and configuration of stations-hypocenter on the solution, and find that direct solution of the derived linear equations give rise to a satisfactory precision for epicenter location, but a poor one for focal depth. Making use of the fact that the hypocenter distance can be calculated beforehand, we get the focal depth and the original time respectively and independently in a satisfactory precision. . If we add random noise perturbations to the observed arrivals and then take them as arrivals of more dummy stations, a still higher precision can be get with more information. To quantitatively assess the influence of the configuration of stations-hypocenter on epicenter location, we define a“q”. Compared with the classic Geiger’s method, iterations and differential calculations are not required in our method. 一、引 言 地震定位是地震学中最经典、最基本的问题之一，对于研究诸如地震活动构造、地球内部结构、震源的几何构造等此类地震学中的基本问题有重要意义。此外， 基于快速准确的地震定位的地震速报，对于震后的减灾、救灾工作也是至关重要的。因此，地震学家一直在不断改进或提出新的地震定位方法。 现行的线性地震定位方法大多是源于Geiger的经典方法[1]，将问题归结为求一目标函数的极小值，并通过对走时求偏导数做一阶近似，将非线性问题转化为线性问题。在此基础上，产生了一系列地震定位方法：HYPO71[2]，HYPOINVERSE[3]，HYPOCENTER[4]，联合定位（JHD）[5]，相对定位[6]等。 这里，我们提出一种不同于Geiger的方法，它通过各台站P波和S波的到时差与台站和震源距离的关系得到完全线性的方程组，由此确定震中，并且分别对震源深度和发震时刻单独求解，不需要对走时求偏导数，求解过程不需要迭代。 二、基本原理 设震源坐标为