论文部分内容阅读
基于波动方程的层析成像方法被认为是研究地球内部结构和动力学过程的新一代成像方法。首先,本文利用伴随方法、Born近似以及格林函数等发展了基于波动方程的地震层析成像和震源定位联合反演方法。针对该方法需要求解正演波场和伴随波场来计算层析成像的敏感核的问题,本研究提出了具有四阶精度的高效率、低存储需求、低数值频散、易并行计算的三维近似解析中心差分方法,分析了其稳定性和数值频散关系,并针对不同的介质模型给出了地震波传播的波场模拟结果。GPU/CPU异构并行加速试验表明该方法具有优良的并行加速效果。随后,利用推导的差分方法在一个地壳-Moho面-地幔模型中计算了不同地震波震相,包括直达波、反射波、多次反射波以及首波的波形和走时敏感核。结果表明,无论是波形敏感核还是走时敏感核,都主要分布在以地震波传播几何路径为中心的第一菲涅尔带区域内;后续震相的敏感核不仅存在于第一菲涅尔带,也包含反射点或散射点的等时线。基于波动方程层析成像敏感核包含等时线是其能天然解决多路径问题的原因之一。在应用方面,利用基于波动方程的走时层析成像对北京及其周边区域地壳结构进行层析成像研究。成像结果显示该区域地壳的强非均匀性与地震的发生存在关联性。唐山地震震中位置下方的低速异常可能表明流体的存在;唐山地震及其余震的发生与流体对孕震层的弱化作用以及扩散、运移有关。在对美国南加利福尼亚Landers地震区域的层析成像研究中,尝试了同时利用基于波动方程的层析成像方法和Moho反射震相对该区域整个地壳进行成像。检测板测试结果表明Moho反射震相能够显著提高下地壳的层析成像结果分辨率。高分辨率成像结果揭示了近地表的非均匀性与地表地质构造有关;地震多发生于高速、低泊松比异常区域或者高低速异常交界区域。成像结果揭示的中下地壳中沿圣安德烈斯断层东南端靠近索尔顿海槽存在明显的低速异常,可能是岩浆扩张对地壳入侵和驱动作用导致下地壳存在流体所致;其北端的高速异常为含流体区域的边界。在Big Bear地震震中下方存在一个高速P波,低速S波,高泊松比异常,该异常的形成可能与受上地幔中残留的古法拉隆板块下沉过程中对部分熔融的软流圈物质驱动所形成的低速异常对下地壳侵袭作用有关。