利用ICCG迭代技术加快电阻率三维正演计算.docVIP

利用ICCG迭代技术加快电阻率三维正演计算　摘要　一般而言，有限差分法求解点源三维地电场正问题所形成的大型稀疏线性方程组Ax＝b，直接解法的计算效率极低。本文从系数矩阵A的不完全Cholesky分解及矩阵特征值的特点等角度，说明了不完全Cholesky共轭梯度(ICCG)迭代技术可大大提高电阻率三维正演速度的内在原因。结合矩阵A的稀疏存储模式，使得内存需求也大大减少。　　关键词　电阻率法勘探　三维正演　预处理　不完全的Cholesky共轭梯度法　　中国图书资料分类法分类号　P631.34 3-D RESISTIVITY FORWARD CALCULATION ACCELERATEDBY ICCG ITERATION TECHNIQUE Wu Xiaoping　Xu Guoming　(University of Science and Technology of China) 　Abstract　For the large sparse linear equations:Ax＝b，which are formed from the finite difference method used to solve the 3-D forward problem of geoelectrical field,in general,the computational efficiency with direct method is quite slow.In this paper,from the view point of the characteristics of the incomplete Cholesky decomposition of matrix A and its eigenvalue,the internal cause of greatly increased 3-D resistivity forward speed using the incomplete Cholesky conjugate gradient (ICCG) iteration technique is explained .Introducing the row-indexed sparse storage mode to store matrix A,the internal storage demand is greatly decreased.　Keywords　resistivity prospecting;3-D forward;incomplete Cholesky conjugate gradient;sparse storage mode 1　引言　　随着近年来高密度电测技术的应用，传统的直流电阻率法又展现了新的活力，成为浅层水文、工程、环境、考古等与人类社会生活密切相关的探测领域的重要手段［1，2，3］。而浅层地质目标多表现为复杂的三维电性结构，对解释带来很大困难。典型的如大地电磁测深(MT)中的静位移效应［4］，其根本的解决方法就是三维反演解释，这也是高精度电法发展的趋势［5］。　　众所周知，正演是反演的基础和前提。而大多数情况下，地电三维正演只能通过有限差分［6］、有限单元［7，8，9］等求得数值解，因而探寻高效快速的三维正演计算就显得尤为重要。基于有限差分的数值计算方法，笔者曾引入不完全Cholesky共轭梯度(ICCG)方法［10，11］，同时结合系数矩阵的稀疏存储模式，求解二、三维地电场的正问题，在提高计算速度和减少内存需求上均获一定成功，取得初步效果。本文从讨论系数矩阵A的不完全Cholesky分解及矩阵特征值的特点出发，探讨了利用ICCG迭代技术加快电阻率三维正演计算的内在原因，希望能为ICCG方法在地球物理三维正演计算中的更好应用和发展起到抛砖引玉的作用。　　2　点源三维地电场的有限差分计算　　设点电源的电流强度为I，位于坐标点(x0，y0，z0)处，则其产生的点源三维地电场电位φ(x，y，z)满足微分方程： 其边界条件为：  式中　σ(x，y，z)——地下介质的电导率；　　　　　n——边界外法线方向的坐标变量；　　　　　θ——n和r的夹角；　　　　　n——边界外法线方向的单位矢量；　　　　　r——源点到边界上点的向径。　　上述边值问题用有限差分的数值计算方法求解［6，8］，可获得复杂三维结构上的地电场分布。只是需对整个研究区域进行Nx×Ny×Nz的三维网格剖分，未知节点个数太多，最后形成一大型稀疏线性方程组 Aφ＝S。(2) 其中A为容量矩阵，是一大型稀疏对称正定带状矩阵，有如下形式： 这里Ct、Cb0、Cl、Cr、Cf、Cba、Cp分别是节点(i，j，