为实现井—地激发极化数据的快速反演,分析三维有限元正演四面体网格剖分形成的系数矩阵的元素规律,应用MSR非零元素压缩存储和SSOR-PCG方程求解技术,提高了正演执行效率。在正演基础上,利用Jacobian矩阵与正演方程的简单关系,采用共轭梯度法求解最小二乘目标函数的模型修改量,避开Jacobian矩阵的直接计算和存储,实现了三维井-地激电数据的分钟级快速反演。模型正演合成数据的反演结果表明,能基本再现已知的模型。
Based on an analysis of coefficient matrix formed from three-dimensional FEM(finite element method) forward with tetrahedral cells division,the authors used MSR(Modified Sparse Row) memory method and SSOR-PCG for solving equations to realize the fast forward of borehole-surface IP surveying.The inversion technique is based on the least-squares and CG(conjugate gradient),and it makes use of the simple relationship between Jacobian matrix A and the forward equation,which needs only to solve and multiplying a vector and avoids direct calculation and access of matrix.The 3-D fast inversion with minutes of PC time is achieved.The inversion results of three-dimensional borehole-surface forward synthetic data reproduce the known models approximately.