有限频体波走时层析成像及其在华北地区的应用

张风雪

（中国地震局地球物理研究所，北京100081）

有限频体波走时层析成像及其在华北地区的应用

张风雪

（中国地震局地球物理研究所，北京100081）

自层析成像被引入到研究地球内部结构以来，地震层析成像发展迅速，现已成为研究地球深部结构的有效方法之一。按层析成像的理论依据，可以把地震层析成像分为射线层析成像和有限频层析成像。本文旨在掌握一种有限频灵敏度算核的计算方法，进而构建有限频层析成像的反演计算，并用有限频层析成像的方法研究华北克拉通上地幔的速度结构及其动力学意义。

本文的主要研究内容及结论如下：

在射线层析成像方面，本文改进了FMM（Fast Marching Method）射线追踪方法，这种方法可以得到整个空间的完整走时场，并且求得的射线路径是最速路径。本文在算法上对FMM方法求取射线路径的过程进行了改进，使其求解速度更快。

有限频层析成像是随着宽频带数字地震仪的发展而兴起的一种新方法，该方法克服了射线理论中“无限高频”假设所带来的弊端，考虑更接近于真实情况的地震波频率衰减、波散射及波前复原等情况。本文使用改进的FMM射线追踪方法构建了走时表，并在此走时表的基础上，采用查表算法计算一维速度模型下的灵敏度算核，结果表明：查表算法可以获得较高的计算效率，为下一步开展有限频层析成像提供了有力的支持。

本文讲述了波形的多道互相关法和自适应叠加法拾取相对走时残差的基本原理，实现了两种方法的图形界面操作，简化了操作过程，方便结果的对比，剔除差别较大的数据，从而提高了相对走时残差的可信度。

在实际地震资料方面，本文从华北地震科学台阵的190套宽频带地震仪，在2006年10月—2009年3月期间所记录的震中距在30°～90°的地震事件中选取多频段的相对走时残差，进行P波和S波的速度结构的成像反演，得到如下结论：

（1）研究区域西部的低速异常延伸至约400 km，此低速异常反映了深部存在的热物质，本文的成像结果发现，此低速异常没有贯穿地幔转换带到达下地幔，所以，本文认为，此低速异常并不是来源下地幔，至于这个低速异常的形成原因不能单从速度结构一方面去推断，应结合其他参数（热值、密度、Q值、物质成份等）进行综合的分析。

（2）区域地震层析成像研究发现，太平洋俯冲板块前缘已经抵达约120°E处。本文研究发现，在地幔转换带中也有一个高速异常体，推断其为太平洋俯冲板块的前缘。

（3）华北盆地为裂谷型盆地，其速度异常的分布较为复杂，在100～250 km范围内存在岩石圈拆沉的高速异常残留体。

（4）与古老的克拉通岩石圈厚度相比，华北地区的岩石圈正在减薄，华北平原地区表现为高速异常与低速异常交互分布的复杂现象，而燕山地区则为高速异常，表明华北平原地区发生了明显的侵蚀与破坏，燕山地区则没有发生侵蚀与破坏过程或此过程发生的现象不明显。

射线理论；有限频理论；层析成像；华北；克拉通；走时残差

（作者电子信箱，张风雪：zhangfengxue336＠163.com）

P315.3；

A；

10.3969/j.issn.0235-4975.2012.01.014