大地震破裂成像的时间域和频率域反投影方法研究及应用

摘要

从过去到现在，地震一直是困扰着人类的重大自然灾难，这其中，以大地震造成的灾害最为严重，因此研究大地震的孕育和发生过程及其产生的灾害特征一直是地震学研究的核心问题。为了更好地理解和研究大地震，科研人员提出了多种地球物理学方法，以更好地从不同角度研究大地震。其中，基于地震学的反投影方法可以利用区域高密度台阵或全球台阵来还原地震破裂过程中能量辐射的时空分布特征，并以此为基础研究破裂过程。
　　本文从信号传输理论出发，从信号传播的视角展示了反投影的基本原理和处理过程，自然地推导出了反投影的优缺点以及所依赖的数据特征。在理论推导基础上，介绍了传统时间域和频率域反投影方法的数据处理过程，并介绍了常见的各种处理流程上的变化及后续发展。
　　为了验证多台阵时间域反投影方法的有效性，本文用多台阵反投影方法处理了北美、澳洲和欧洲三个台阵的波形数据，研究了2015年4月25日尼泊尔Mw7.8级地震的破裂过程和高频能量辐射特征。结果显示，高频能量辐射位于高滑移区的北部，破裂从震中开始沿断层走向向东东南方向传播，速度约为2.8km/s，各个台站所观测的结果基本一致。将反投影显示的高频辐射与余震分布和滑移分布结合起来，本文最终得出结论，受深度控制的温压条件控制着同震高频和低频辐射的分布特征，而破裂区东部与西部的应力和结构不均匀控制着破裂的开始与结束。
　　在传统反投影方法中，通常只考虑基于一维模型的理论走时和从初始震中到台站下方的三维结构差异，而将源区结构视为均匀和各向同性。为了考察源区结构对反投影的影响，本文提出了一种较为系统的基于余震的源区走时校正方法，并用这种方法处理和研究了2011年日本Tohoku Mw9.0级地震的破裂和高频能量释放过程。本文通过对比进行源区结构校正和没有进行源区结构校正的反投影结果发现，校正后的高频辐射位置向靠近海沟远离海岸的方向有所移动，离滑移区更近。这说明Tohoku地震破裂区确实存在着结构的横向非均匀性，也说明了本文提出的校正方法的有效性。这种校正方法可以结合其他先验信息和假设，具有很强的适应性。
　　除了时间域的多台阵反投影方法，本文还研究了频率域的多台阵稀疏反演（压缩感知）反投影方法。这种方法结合基于l2,1范数的稀疏优化，将频率域压缩感知反投影方法从单台阵推广到了多台阵，发展了一种新的基于多观测量的压缩感知反投影算法。这种方法只假设不同台阵观测的能量辐射来自相同的稀疏位置，而对不同台阵观测到的辐射源的振幅和相位没有限制，从而合理有效地利用了多台阵观测反演得到更为准确的震源信息。为了验证这种方法的有效性，本文进行了数学和模拟数据测试，测试结果都表明方法的有效性和可靠性。
　　今后，将多台阵的频率域压缩反投影和源区走时校正结合起来，可以有效地提高反投影结果的真实性、可靠性和准确性，为认识大地震破裂的物理过程提供更好的约束。

目录

声明

摘要

第1章 引言

1．1 大地震的危害及研究的困难

1．2 目前对大地震的主要研究手段及认识。

1．3 本文结构

2．1 概论

2．2 理论基础

2．3 传统反投影

2．4 频率域反投影

2．5 基于稀疏优化的频率域反投影

2．6 多台阵反投影

2．7 反投影方法总结

第3章 用多台阵反投影方法研究2015尼泊尔Mw 7．8 地震的破裂过程

3．1 概述

3．2 数据与方法

3．3 结果

3．4 讨论

3．4．1 多阶段破裂过程

3．4．2 同震高频辐射、滑移、余震和结构尺度的关系

3．5 总结

第4章 利用余震进行走时校正的反投影方法：以2011年日本Tohoku Mw 9．0 地震为例

4．1 概述

4．2 数据与方法

4．3 结果

4．4 讨论

4．4．1 绝对走时与相对走时

4．4．2 与其他余震走时校正方法的比较

4．4．3 可行性、可靠性及提升

4．5 总结

第5章 基于稀疏优化的多台阵频率域反投影

5．1 概述

5．2 方法

5．3 模拟测试

5．3．1 数学测试

5．3．2 用Ricker子波生成的模拟波形的测试

5．3．3 用有限断层生成的模拟波形的测试

5．4 讨论

5．4．1 网格间距

5．4．2 对不同震源机制的适应性

5．4．3 如何理解观测到的相位变化

5．4．4 模拟测试与处理真实数据之间的距离

5．5 总结

6．1 总结

6．2 存在的问题

6．2．1 深度震相的影响

6．2．2 二维结果与三维真实破裂

6．2．3 如何理解多台阵结果的相似与差异

6．3 展望

参考文献

附录

致谢