基于预条件共轭梯度法的直流电阻率三维有限元正演研究 - 岩土工程学报.pdf

32卷 12期 岩 土 工 程 学 报 Vol.32 No.12 2010年 .12月 Chinese Journal of GeotechnicalEngineering Dec. 2010 于预条件共轭梯度法的直流电阻率三维有限元 正演研究 * 斌，李术才，李树忱，聂利超 (山东大学岩土与结构工程研究中心，山东 济南 250061) 摘 要：对直流电阻率法勘探而言，对其反演结果的精度和速度的要求越来越高，这就需要提出一套优化有限元数值 正演速度和精度的计算方案。设计了系数矩阵的一维非零元素压缩存储模式，设置了索引数组以便按照行号和列号对 元素 行索引，与变带宽存储模式相比，其内存占用量明显减小。为提高正演计算速度，利用预条件共轭梯度法（PCG） 求解有限元中的大型稀疏线性方程组。在PCG 法中，将雅可比迭代中的对角阵作为预处理矩阵，与其它预处理矩阵相 比，其具有求逆方便、无需存储空间的特点，使得大型线性方程组的求解速度大大提高。另外，在直流电阻率三维正 演中，采用了异常电位法，提高了电源点附近的解的精度。以二层地层的电阻率勘探为例，初步验证了计算方案的实 用性。利用上述方案，重点对隧道含水断层的电阻率法超前探测 行了有限元数值正演，并 行了相应的物理模型试 验。对比显示，数值正演结果与试验数据基本一致，且数值正演的速度和精度均显著提高。 关键词：直流电阻率勘探；3D FEM ；PCG；预条件矩阵；一维非零元素压缩存储；含水构造超前探测；模型试验 中图分类号：TU45 文献标识码：A 文章编号：1000 –4548(20 10) 12 –1846 –08 作者简介：刘 (1983 – )，男，山东高唐人，讲师，主要从事不良地质体超前预报与治理、电阻率法与激电法的正 反演等方面的研究。E-mail: liubin0635@163.com 。 3D FEM numerical forward modeling of direct current electrical resistivity based on PCG algorithm * LIU Bin, LI Shu-cai, LI Shu-chen , NIE Li-chao (Geotechnical Structural Engineering Research Center of Shandong University, Jinan 250061, China) Abstract: For the direct current electrical resistivity prospecting, the requirements for the precision and the speed of inversion are higher and higher. Thus, an optimization scheme is proposed to improve the precision and the speed of FEM numerical modeling. 1D compression storage mode