[发明专利]一种用于由卫星磁异常反演居里等温面埋深的方法

摘要

本发明公开了一种用于由卫星磁异常反演居里等温面埋深的方法，首先，选择全球初始CPD分布，代入到CPD埋深；其次，由CPD分布，经ESD正演计算得到全球全波段的岩石圈磁场(LMF)分布，再其次，采用高通滤波滤去计算岩石圈磁场中的长波长，得到研究区域内中短波长的岩石圈磁场；接下来，将滤波后的岩石圈磁场与观测数据或者模型进行对比，如果二者小于设定的反演误差，则输出CPD埋深；否则，进行非线性共轭梯度(NCGM)反演，得到新的CPD；将这一迭代过程不断循环，直到得到理想的输出结果。本发明解决了如何由卫星岩石圈磁异常观测或者模型获取区域性居里面埋深，使得反演结果可靠，具有更高分辨率和更高精度，反演算法快速收敛。

主权项

1.一种用于由卫星磁异常反演居里等温面埋深的方法，其特征在于，包括如下步骤：首先，选择全球初始CPD分布，代入到CPD埋深；其次，由CPD分布，经ESD正演计算得到全球全波段的岩石圈磁场(LMF)分布，在ESD正演过程中，任一场点位置处的LMF近似为其所对应的正下方地壳内一定区域内的ESD叠加效果，考虑磁场三分量对磁异常的贡献，以及球面几何的影响，计算过程中采用DA公式；再其次，由于LMF的长波长部分被主磁场所掩盖，一般的LMF模型给出的都是中短波长部分，采用高通滤波滤去计算岩石圈磁场中的长波长，得到研究区域内中短波长的岩石圈磁场；接下来，在研究区域内，将滤波后的岩石圈磁场与观测数据或者模型进行对比，如果二者小于设定的反演误差，则输出CPD埋深；否则，进行非线性共轭梯度(NCGM)反演，得到新的CPD；最后，将这一迭代过程不断循环，直到得到理想的输出结果。