论文部分内容阅读
地下各种具有不同磁性差异的岩(矿)石会在其周围产生不同的磁性异常场,磁力勘探是通过分析解释观测所得到的磁异常数据,用以勘察有用矿产及寻找地下地质构造的一种较为常用地球物理方法,具有广泛的应用。磁测数据的采集简单易行、数据采集效率高,且不受地域的影响。磁异常反映了地下磁性特征的变化,包括磁性异常体的分布和磁化强度。通过有效的反演方法,结合已知的地质资料,就可以对磁性异常体的大小、形状、磁化强度等信息进行估算。 本文在不考虑磁性异常体的剩磁及退磁影响的情况下,对二度磁性异常体进行了磁化强度成像研究。在实际工作中,由于二度体和三度体在剖面上可能有相同的异常形态,而本文所介绍的方法中,磁异常数据为一维(剖面)数据,所有地质体都当作二度体的来解释,部分解释结果与实际情况不一致是在所难免的,若要进行更为细致和准确的解释,那么建议采用平面异常数据,而对于位场反演问题来说,三维反演与二维反演相比,并没有本质上的差异,实现三维反演仅需要做细节上的调整便可,反演方法并没有太大的改变。 在研究了二度磁异常体正演公式的基础上,本文给出了二维磁化强度成像的基本公式。在构造目标函数时,基于Occam’s反演方法,在目标函数中加入了光滑度项及最小模型项来对模型进行约束,同时为了克服反演结果富集在地表附近的现象,采用深度加权函数对空间单元赋予不同的权重,得到比较符合实际的结果。本文还讨论了如何选择最优的正则化参数,通过选取合理的正则化参数,使反演的结果较为合理。在求解方程组时,采用了计算速度快、节省内存且收敛速度快的共轭梯度算法。 通过对多组理论模型的正演计算及其磁化强度成像结果的对比分析,讨论了反演算法的分辨率和可靠性,并对在反演中所选择的参数进行了总结和归纳。本文的反演算法不依赖于初始模型设置,成像效果较好。将该反演方法应用于航磁剖面,反演的结果为构造的解释提供了一定依据。