利用symmetric successive over-relaxation(SSOR)预处理共轭梯度法的三维位场数据反演

摘要

随着位场数据测量在现代地球物理勘探中的发展,及广泛的应用前景,位场数据的定量解释越来越受到人们的重视.近年来,三维反演技术在位场数据处理等相关领域也得到各个地球物理学家的重视.文本中，首次引用了一种新的预处理共轭梯度法（SSOR），这种方法可以很好的解决大数据量的位场数据反演方法：本文提出的预处理方法的预处理矩阵可以直接从系数矩阵中获得，同时不会有任何附近的计算时间的消耗。不需要过多的计算机内存来存储系数矩阵的逆。由于SSOR预处理矩阵本身就是系数矩阵的一部分，在反演算法中，它可以更加有效的应用在矩阵向量的乘法中，进而可以有效的减少反演的迭代次数。并通过简单的地球物理正反演模型试验，与传统的对角线预处理方法，以及不完全Choleshky分解方法预处理等预处理方法进行对比，从理论模型试验中可以发现这种新的SSOR预处理方法的反演迭代次数要明显的少于其他的预处理方法，说明SSOR预处理反演方法可以很好的提高位场数据反演计算效率。