接收函数近邻反演方法的改进和对海拉尔台下地壳速度结构的研究.pdfVIP

维普资讯 第22卷第 3期 中 国科 学 院研 究 生 院 学报 Vo1．22 No．3 2005年5月 JournaloftheGraduateSchooloftheChineseAcademyofSciences Mav 2005 文章编号：1002．1175(2005)03．0322—07 接收函数近邻反演方法的改进和 对海拉尔台下地壳速度结构的研究* 沈旭章 周蕙兰 (中国科学院研究生院计算地球动力学重点实验室 ，北京 100049) (2004年7月2813收稿；2004年 11月2531收修改稿) ShenXZ，ZhouHL．Improvementoninversionmethodofneighborhoodalgorithm withreceiverfunctionsandresearchofcrust structureunderstationHLR．JournaloftheGraduateSchoolofhteChineseAcademyofSciences，2005，22(3)：322—328 摘 要 数值试验结果指出，用接收函数反演地壳速度结构时，在正演和反演过程中使用不匹 配的地壳模型参数化方式，将对结果产生不良影响；同时发现，在处理实际资料时将接收函数 归一化，可能会丢弃介质信息．在此基础上，改进 了接收函数和近邻算法反演地壳速度结构的 程序，并用来反演了海拉尔台下的地壳速度结构． 关键词 接收函数，近邻算法，地壳结构，模型参数化 中图分类号 P313．2 1 引言 近年来 ，随着观测数据的积累和高性能大规模并行计算技术的发展，耗时多但不依赖初始模型的非 线性反演方法的研究和应用，有了更多的发展 ． 在用接收函数研究台站下方地壳上地幔结构的非线性反演技术中，全局最优解的寻找方法、算法的 搜索效率以及多解 的存在都是人们关注的重要问题 ．针对这些问题，人们提出了不同的方法 ．刘启元 发展了接收函数复谱 比的非线性反演方法 ，并用之研究了中国 10个数字地震台站下面的速度结 构 ；Shibutani和Sambridge则使用了遗传算法(GeneticAlgorithm，GA) ，周蕙兰和杨毅提出了剥壳遗传 算法 ，将接收函数的反演深度扩展到 了700km左右；Sambridge提 出了近邻算法 (Neighborhood Algorithm，NA)[73，将其用于接收函数的非线性反演中，并将它与模拟退火法(SimulationAlgorithm，SA)和 遗传算法进行了比较．相对于常用的线性反演方法，这些非线性反演方法的最大优点是不需要选用一个 初始模型，而只需给出模型参数的搜索范围，从而消除了线性反演中对初始模型的依赖．国内外许多学 者使用这些非线性反演方法和接收函数已经得到了不少区域性的研究结果 ． 本文在使用NA方法和接收函数反演地壳速度结构时，发现对地壳模型不合理的参数化方式和接 收函数的归一处理，都会对反演结果有不 良影响．为此，本文在数值试验的基础上对反演程序进行了改 进，并用海拉尔台站的实际资料反演了它下面的速度结构． 2 地壳结构的NA反演方法 NA这种近年来兴起的非线性反演方法 ，可以在最优解的搜索过程中，充分利用前一次迭代的信 * 国家 自然科学基金 资助 十通讯联系人 ，E．mail：zhoutd@gscas．ac．ca 维普资讯 第 3期 沈旭章 ，周蕙兰：接收函数近邻反演方法的改进和对海拉尔台下地壳速度结构的研究 323 一∞，LI互一 1f 息；而且，只需要改变2个参数，就可以在模型空间中产生许多复杂的搜索．同时，该方法又具有5 天然5 的 4