论文部分内容阅读
直流电阻率法和激发极化法在浅层水文、找矿、工程、环境等领域应用广泛,而浅层的地质目标体一般表现为三维电性结构,一维和二维的数据反演不能完整的恢复三维电性结构,为了避免在反演中三维体不均匀性的影响,需要采用三维数据观测和反演。在找矿方面,基于矿山中钻井的相对位置,我们可以将井中装置与地表装置相结合,即形成地表、地井、井地和井间四种观测方式相结合的复合采集方式。由于三维地球物理问题不仅是欠定的、非线性的,在进行三维反演时还存在耗时长的特点,而且线性反演算法往往使得反演陷入局部极小,所以本文采用非线性共轭梯度反演算法来规避使用线性反演时所带来陷入局部极小的问题,并且在该反演算法中加入MPI和GPU并行算法来提高计算效率。针对直流电阻率法和激发极化法的正演问题,本文采用有限差分法求解点电源三维地电场。针对正演算法中的解方程,采用不完全Cholesky分解的共轭梯度法和不完全LU分解的稳定双共轭梯度法;针对大型系数矩阵,本文采用行压缩存储的格式,不仅减小了内存消耗,还提高了计算的效率。反演算法采用非线性共轭梯度法,非线性共轭梯度法的优势就在于不用求雅克比矩阵,而是直接求取雅克比矩阵或者其转置和一个向量的乘积。该方法不仅避免了存储雅克比矩阵,还大大减少了计算量,每次反演迭代计算一次正演,两次“拟正演”。在并行计算方面,由于MPI含有将多个独立任务分配给多进程同时计算的特性,即可利用MPI进行多个点电源地电场的计算。而在计算正演和“拟正演”过程中,都存在解方程的过程,反演过程的大部分时间都消耗在解方程中,因此,利用GPU快速运算的特性,在解方程时,调用基于GPU开发的CULA Sparse函数库参与解方程运算。CULA Sparse函数库包含多种解方程方法和预处理方法,在本文中,我们采用不完全LU分解的共轭梯度解方程。最后,分别设计单一采集装置和多个将地表、地井、井地和井间相结合的装置,通过模型算例验证实现的非线性共轭梯度并行算法的有效性,并说明了多装置相结合采集的数据参与反演比使用单一装置采集的数据参与反演的结果更好。