地震定位研究综述-2017年11月.docVIP

地震定位研究综述 作 者 田 玥 指导教师 陈晓非 摘要 综述了各种地震定位方法的基本原理，重点介绍了Geiger的经典方法以及在此基础上建立的各种线性方法：联合定位法，相对定位法，和最新的双重残差法；对每一种方法的应用情况，尤其是国内的工作做了总结；同时也指出了各种方法的特点，并进行了相应的比较。此外，还简要介绍了空间域的定位方法和各种非线性定位方法。 关键词：震定位;线性定位;非线性定位 目 录 引言 1 1经典定位方法 1 1.1经典方法 2 1.2各种改进方法 3 2经典定位方法 4 2.1震源位置与台站校正的联合反演(JED， JHD) 4 2.2 震源位置与速度结构的联合反演(SSH) 5 2.3 相对定位法（主事件定位法，ATD） 5 3空间域内的定位方法——台偶时差法 6 4非线性定位方法 7 4.1 牛顿法 7 4.2 全局搜索方法 7 4.3 Bayesian方法 8 5最新定位方法 8 5.1 EHB方法 8 5.2 双重残差法（DDA） 8 结论 9 致谢 9 参考文献 10 引言 地震定位是地震学中最经典、最基本的问题之一，对于研究诸如地震活动构造、地球内部结构、震源的几何构造等此类地震学中的基本问题有重要意义。此外， 基于快速准确的地震定位的地震速报，对于震后的减灾、救灾工作也是至关重要的。因此，地震学家一直在不断改进或提出新的定位方法。 地震定位问题的提法如下：根据台站对地震到时的观测资料，来确定震源的空间坐标和发震时刻，有时还给出对解的评价。早期的地震定位方法以几何作图法为主[1]。近三十年来由于计算机技术的飞速发展和广泛应用，基于科学计算和计算机技术的智能化数值自动定位方法也得到了迅速发展，并业已成为当前地震定位的主流方法。我国最初的地震定位工作由李善邦先生于1930年在北京鹫峰地震台开创，1953年开始采用多台站大规模观测数据确定震中，现在大多使用国际流行的定位方法。 本文只介绍目前广泛使用的计算机定位方法，重点介绍Geiger的经典方法以及在此基础上建立的各种线性方法：联合定位法，相对定位法，和双重残差法，并且重点总结了国内的有关工作。 1经典定位方法 1.1经典方法 现行的线性定位方法大都源于1912年Geiger提出的经典方法[2]： 设n个台站的观测到时为 求震源及发震时刻，使得目标函数 (1) 最小。其中为到时残差 ， (2) 为震源到第i个台站的计算走时。 使目标函数取极小值也即 ， (3) 其中. 为方便，记 ， (4) 则由(3)式，在真解附近任意试探解及其校正矢量满足 . (5) 也即 . (6) 由的定义可得公式(6)的具体表达式 . (7) 若偏离真解不大，则和较小，可忽略二阶导数项，(7)式被简化为线性最小二乘解： . (8) 以矩阵形式表示，上式为 ， (9) 其中 . 若二阶导数项不可忽略，则(7)式给出非线性最小二乘解 . (10) 通常各台站的到时数据具有不同的精度，如果不加以区别，则具有较低精度的数据将严重干扰结果的精度，这一问题可以通过引入加权目标函数来解决。设各台站到时残差的方差为，引入加权目标函数 ， (11) 按照上述同样的步骤，通过求(11)式的极小值，得到如下加权线性最小二乘解 ， (12) 其中为加权方差矩阵：. 由方程(9), (10), 或(12)求得后，以作为新的尝试点，再求解相应方程。如此反复迭代，直至足够小（或满足一定的循环结束条件），此时即得估计解。 1.2各种改进方法 直到20世纪70年代，随着计算机的迅速兴起，Geiger的思想才被广泛用于地震定位工作。Lee等人连续给出了HYPO71, HYPO78~81系列程序[3]，至今仍被普遍使用，我国的赵仲和参与了80、81版本程序的研制。Backus和Gilbert提出新的反演理论后，Klein提出HYPOINVERSE