用遗传算法确定三维横向不均匀介质中的近震震源位置.PDF

第 19 卷　第2 期 西　北　地　震　学　报 Vol.19　No.2 1997 年6 月 NORTHWESTERN SEISMOLOGICAL JOURNAL June 1997 用遗传算法确定三维横向不均匀介质中的 近震震源位置 万永革 (国家地震局防灾技术高等专科学校, 北京东燕郊　101601) 李清河 (国家地震局兰州地震研究所, 兰州　730000) 李鸿吉　丁志峰 (国家地震局地球物理研究所, 北京　100081) 摘要　介绍了 遗传算法在三维横向不均匀介质中确定震源位置的研究思路并 简述了三维横向不均匀介质中走时计算的基本原理.在遗传算法的迭代过程中, 考虑 到初始迭代阶段有较大的参数搜索空间而在最后迭代阶段有较小的参数空间, 提出 了遗传算法中变异概率随迭代次数下降的函数形式并对其进行了分析.结果表明, 随 迭代次数指数下降的变异概率使得算法收敛加快, 而且反演得到的震源位置的各台 站走时残差比原定震源位置的明显减小. 主题词:地震定位　横向不均 一性　遗传算法　变异概率 1　引 准确的震源位置是研究地震活动图象、地震与活动构造的关系以及工程地震等问题的基 础.可靠的地震定位需要区域内准确的速度模型和高效稳定的反演方法. 大多数定位算法基于盖革法, 此方法通过对起始假想震中的迭代微扰逐渐趋于最优解.每 [1] 一步迭代过程需要计算假想震源到台站的走时的偏导数 .这种方法在初始模型适当给定的 情况下可以得到合理解.但当约束不够严格时, 常需要某种形式的数值阻尼来避免不稳定性. 这种情况下得到的震源位置强烈依赖于阻尼的程度.在三维横向不均匀介质中计算地震波走 时的偏导数需要花费较大的代价.另外, 定位问题的非线性常使解依赖于起算模型, 从而使得 算法限于局部极值. 不依赖于求导信息的定位算法可以避免上述限制.Sa bridge等(1986)提出了这种算法. 他们通过分离时间和空间坐标运用一系列网格搜索可达到合理的效率.然而, 对于许多问题进 行可靠的解常需要大量的拟合差函数估计.所以, 当计算三维横向不均匀介质中的走时时, 这 种方法并不实用. 遗传算法是一种新型的完全非线性搜索方法, 近来在人工智能领域内得到广泛关注.在地 本论文为硕士论文的一部分 收稿日期:1996-06-02 第一作者简介:万永革, 男, 1967 年2 月生, 助教, 主要从事地壳结构和震源位置反演研究. 　　　　　　　　　　　　　　　　　　西　北　地　震　学　报　　　　　　　　　　　　　　　　　　第 19 卷 8 球物理学中已解决了许多复杂的最优化问题[2 ～4] .其主要优点是全局搜索, 不用求导数, 不依 赖于初始模型的选取, 普遍适用于各种优化问题.本文中, 我们描述了如何用遗传算法在三维 横向不均匀介质中进行高效稳健的地震定位, 并对遗传算法的基本运行机制进行了修改以使 得避免陷入局部极值的能力有所改善, 收敛的速度有所增加. 2　地震定位问题 假如在时刻t , 在(x , y , z )处发生了一次地震, 那么在台网中观测到的是一组到时.用 0 0 0 0 这些数据求地震的发震时刻和震源位置(t , x , y , z ). 0 0 0 0 为了用台站位置为(xk, y k, z k)(k =1, 2, …, )的一组到时 τk 作地震定位, 我们假定一 ＊ ＊ ＊ 地壳模型, 用它能计算出由位于(X ,