[发明专利]一种接收函数和重力联合反演地壳厚度和波速比的方法

摘要

本发明属于地球动力学研究技术领域，具体涉及一种接收函数和重力反演地壳厚度和波速比的方法。该方法包括如下步骤读入接收函数H‑κ叠加反演的H与κ及区域重力异常数据；选择待处理台站，设置窗口大小及网格间距；将窗口范围内的所有H、κ及重力数据网格化；求取密度参数及噪声参数；采用重力反演的似然估计法计算H‑κ似然谱；将接收函数H‑κ叠加图与重力H‑κ似然图归一化并点乘，得到新的H‑κ谱，从中拾取最优的H和κ作为联合反演结果；选择下一个待处理台站继续前述步骤。本发明通过重力似然图与接收函数叠加图相互加权实现最优H和κ的拾取，为解决复杂情况下接收函数反演的多解性问题提供有效途径。

主权项

一种获取地壳厚度和波速比的方法，所述方法用于降低接收函数H‑κ叠加法的不确定性，提高获取地壳厚度H和波速比κ的精度和效率，其特征在于，所述方法基于接收函数和重力进行联合反演，所述方法首先利用布格重力异常反演获取地壳厚度H初值，进而获得波速比κ初值；然后，对选择的任意一个待反演地震台站为中心进行单一台站联合反演，得到单一台站的联合反演的H‑κ图，从而获取该地震台站的最佳估计的地壳厚度H和波速比κ值；最后，对研究区域内所有待反演站台进行上述单一台站联合反演过程，以获得研究区域内所有待反演地震台站的地壳厚度H和波速比κ；所述联合反演的H‑κ图中地壳厚度H和波速比κ拾取简单；所述方法能够为地壳增厚和地壳变形机制模式的研究提供相对可靠的地壳厚度和波速比数值。