基于PARDISO直接求解器的三维自然电位正反演

近年来自然电位法在海底硫化物资源的勘探和评价中发挥了重要作用。本文开展的是基于PARDISO直接求解器的3D自然电位正反演算法研究。首先，利用有限体积法离散自然电位控制方程，采用PARDISO直接求解器提高正演计算的效率，通过数值解与解析解对比，验证了正演算法的可靠性。其次，在3D反演算法中考虑了地形因素，同时将最小支撑约束与深度加权加入目标泛函中，理论模型数据的反演结果很好地恢复了矿体的结构。最后，利用该算法对室内沙箱实验获得的自然电位数据进行反演，结果显示得出的电流密度异常与金属棒的位置基本一致。因此，本文提出的反演算法在未来大规模自然电位数据反演中具有重要作用。