直流电阻率三维正演的预条件共轭梯度法研究.pdfVIP

·710· 中国地球物理2011 (19)地质调查和矿产勘查地球物理 直流电阻率三维正演的预条件共轭梯度法研究 3一Ddirectcurrent forward the resistivitymodelingby method conjugategradient precondition 杨金凤+邓居智陈辉 JuzhiChenHui Jinfeng+Deng Vang 东华理工大学放射性地质与勘查技术国防重点科学实验室 江西抚州 344000 1．引言 直流电阻率法因其具有灵活性高、勘探费用低等特点被广泛的应用于矿产资源勘查、地质调查等领域。 如何快速、准确的求解离散后的线性方程组是直流电阻率三维正演的核心问题，对三维反演等后续工作具 有重要意义。本文基于有限差分的数值计算方法，引入对称超松弛预条件共轭梯度(SSOR．CG)法，研究直 流电阻率三维有限差分正演问题，并通过异常场法来消除总场中由点电源导致的奇异性。 2．方法原理 为了提高计算精度和计算速度以及避免由源奇异性引起的差分计算复杂性，采用解析计算和数值计算 相结合的方法，即异常场法，推导出不含电流源项的异常场方程。对求解区域进行六面体剖分，然后在网 格点上利用七点差分格式对变形后的积分方程进行离散得到节点电位满足的线性方程组。边界条件采用第 一类边界条件和混合边界条件。 对该方程组的系数矩阵采用对称对角行存储格式进行存储，即用两个一维数组存储。其中一个实型数 组存储矩阵上三角中的非零元，另一个整型数组存储上一个数组中的元素在矩阵中的位置索引。然后采用 对称超松弛预条件共轭梯度(SSOR．CG)法求解。预处理技术的关键是通过合适的预调件矩阵的选择，来降 低系数矩阵的条件数。本算法首先将系数矩阵分解成一个上三角矩阵和下三角矩阵及单位矩阵，然后将预 条件矩阵定义为含松弛因子的一个上三角阵与其转置的乘积。对于松弛因子的选择，根据一定的迭代精度 试算了不同模型最终确定了一个最佳的松弛因子。对于该算法的另一关键步是通过高斯消去法先顺带后回 带最终可求出预处理矩阵与任一列向量的乘积，该计算过程仅涉及少量非零元素的操作，均可以非常简单、 快速地完成。 3．模型验证及结论 首先设计了水平层状介质模型和垂直接触模型进行试算，结果表明采用对称超松弛预条件共轭梯度法 求解有限差分正演中的大型线性方程组，与高斯直接法相比速度显著得到提高，而且内存需求大大减少， 验证了该算法的可行性与有效性。然后设计了典型的高阻(低阻)正方体异常模型进行计算，结果表明将 解析计算和数值计算相结合并利用对称超松弛预条件共轭梯度法求解该正演计算中的异常电位，可以消除 点电源奇异性影响，在保持快速计算的基础上，显著提高了计算精度。 本研究由“东华理工大学研究生创新基金(DYCAl0015)”资助。 参考文献 M P removal：arefinementof 【1】LowryT’AllenB，ShiveN．Singularity resistivitymodeling 1989，54(6)：766．--77．