地球AK135模型中P震相的有限频走时灵敏度算核的计算.pdfVIP

2011 年 6 月 APPLIED GEOPHYSICS （应用地球物理） 第 8 卷，第 2 期 地球AK135模型中P震相的有限频走时灵敏度算核的计算//The computation of a finite-frequency travel time sensitive kernel for P-waves in the AK135 earth model, Zhang Feng-Xue, Wu Qing-Ju, Pan Jia-Tie, Zhang Guang-Cheng, and Feng Qiang-Qiang, APPLIED GEOPHYSICS, 2011, 8(2), P.158 - 163. DOI:10.1007/s11770-011-0280-6 地球 AK135 模型中 P 震相的有限频走时灵敏度算核的计算 1 1 1 1 1 张风雪 ，吴庆举 ，潘佳铁 ，张广成 ，冯强强 （1. 中国地震局地球物理研究所，北京 100081 ） 摘要：有限频走时层析成像是近年发展起来的一种新方法，这种新方法的一个主要过程是走时灵敏度算核的计 算。求解灵敏度算核要多次用到同一散射点的走时，多次地求解同一走时是相当耗时的任务，如果介质为均匀 或速度线性变化等简单模型，散射点的走时可以用解析公式快速地求出，从而灵敏度算核的计算耗时相对较少。 然而各种地球模型中，介质速度大多为分层模型，从解析公式中得到走时信息就比较繁锁。为了提高计算效率， 本文采用查表算法研究地球分层速度模型中的 P 震相有限频走时灵敏度算核的计算，选用的速度模型是地球 AK135 模型，用查表算法求解走时，节省了约 50%的计算时间。在相同的速度模型下，与已有结果的对比，本 文所用的查表算法，能在兼顾精度的前提下，以较小的存储要求换取较高的计算效率，这对提高有限频走时层 析成像算法的速度具有一定的参考价值。 关键词：AK135 ，有限频，灵敏度算核，P 震相 射效应，其结果在震相走时上也会反映出来。另 引言 外，随着传播距离的增加，波前由于非均匀介质 发生延迟或提前的现象，也将随着波前复原逐渐 自1977 年 Aki 和 Lee （1976）将医学层析成 消失（Nolet 和 Dahlen, 2000; Hung 等, 2001 ）。 像应用于地震学研究以来，利用地震波的走时层 针对以上问题，Dahlen 和 Hung 等提出了“香 析成像研究取得了显著地进展，它在全球地震学 蕉——甜甜圈” （Banana-Doughnut Kernel ）理论 和地震勘探中的应用已经取得了一批重要成果 （Dahlen 等, 2000; Hung 等, 2000 ），基于此理论 （Zhao 等, 1994; Huang 和 Zhao, 2004, 2006 ）。传 的层析成像又称有限频（finite-frequency ）走时层 统走时层析成像反演基于射线理论，射线理论中 析成像，它与射线理论走时层析成像的主要区别 一般假定地震波的频带是无限高频的，按照费马 是:射线理论将地震波看作是无限高频的，用一条 原理，即地震波沿着耗时最短的射线路径传播， 狭窄的几何射线路径来代表地震波的传播路径， 其到时在不同观测台站的差异完全取决于射线路 地震波的走时主要受这条射线路径上速度结构变 径上的速度结构，即观测时间仅受连接震源和接 化的影响，而射线路径以外的速度扰动对走时没 收点的大圆下方射线路径上介质速度结构的影 有影响；然而有限频理论认为，地