三维频域电磁响应的数值模拟中迭代算法的计算效率的研究.pdfVIP

三维频域电磁响应的数值模拟中迭代算法的计算效率研究 沈金松 (石油大学盆地与油藏研究中心北京102200) I．引言 尽管过去的十多年中发表了许多三维电磁正演模拟算法，但对离散化大型稀疏线性方程组解法的计算 效率的研究极少。事实上．在三维电磁正演模拟和反演计算中，离散化线性方程组的有效求解已成为整个 处理系统的瓶颈。本文通过交错网格有限差分的离散化问题，讨论了定常松弛迭代算法(SOR／和非定常迭 代的Krylov子空间迭代算法的计算效率和收敛稳定性等问题。利用数值模拟考察了迭代计算中消除伪解 的散度校正方法和加速收敛的不完全乔累斯基分解预处理手续。文中所用的交错网格有限差分将三维介质 离散成正六面体单元，每个单元为一个均匀电性体．将离散电场分量定义在正六面体单元边的中点．离散 磁场分量定义在正六面体单元每个侧面面元的中心。文中还将总场分离成背景场和二次场，背景场用解析 方法求取，二次场用有限差分实现。 2．麦克斯韦(Maxwell)方程的离散化及其选代解法 假设场的时间变化为一，其中F／=i，。=2”，是圆频率，频率域中麦克斯韦方程为： V×E…icoaHMp和vxH=iⅢ(￡一i-o)E+』， (1) 埘’ rrXl0 式中．，·是磁导率，,uo=4 H／m；￡一6DEr是介电常数；e。是相对介电常数；口是电导率；E和H是 电磁场；』。和M。分别是电流密度和磁流密度。 以磁偶极子(．『p≠o和M，=o)为例导出其离散化形式。由式(1)得到二阶电场矢量的方程： V×叮xE (2) k2,E=，。 其中，层一0∥(e～i旦)为介质的复渡数。 tO 为丁消除源点附近的奇异性，将总场￡分离成背景场E。和二次场p，二次场满足如下方程： 可×V×f产t；E3=(睇一薛)‘ (3J 式中，Fp是均匀全空间或半空间的背景场；背景介质的参数为唧和e。。 上述方程用交错网格有限差分法离散得到线性方程，其系数阵既非正定，也不是Hermilian矩阵，传 统的共轭梯度算法不能直接应用。针对该系数阵的特殊结构，选用了松弛迭代法(SOR)和适用于非Her 顶处理和散度校正加速收敛手续。 次迭代得到的解与高斯一赛德尔迭代得到的解加权平均。用矩阵形式表示有： 』：”一(D ∞L)_[甜U十(1一埘)D]J“”+Ⅲ(D一ⅢL)一B(4) 式中，f)是离散化方程系数阵的对角元索组成的对角阵，L和u分别是系数阵的下三角元素组和上三角 元豢组成的矩阵，B是右端向量，。是松驰因子。 松弛迭代解法中，松驰因子“在(o，2)范围内时，迭代算法保证收敛，而且m值的大小决定了迭代算 法的收敛特性。 (21 稳定陆。所用的方法为：正交双共轭梯度方法CoBiCG．广义乘积型双共轭梯度方法GPBiCG，加阻尼的 ———408—— 理建立的。Krylov子空问迭代解法是在共轭梯度迭代解法基础上发展起来的。最初是针对正定行列为系 数阵的线性程组提出的。其基本思路是由线性方程组构造一个二次泛数： F(x)=—昙(x，AX)一(x，B) (5) L 通过构造适当的解空间和选择适当的搜索方向，使第k步搜索方向Pt与先前的搜索方向P㈠tm! 样上述n个残差构成了一个n维的非零向量空间，要使第n+1步的残差^+-与前面所有的残差向量垂直， 必有h+。一0，即对n元线性方程组须n步迭代即可得到真解。 3．迭代解法的加速收敛手续 对于Krylov子空间迭代解法，基于不完全分解的预处理是加速收敛的有效手段。通过预处理使系数 矩阵接近于单位阵，改善其谱分布特性，使共轭梯度迭代方法快速收敛。 仪利用系数阵的对角块进行不完全乔累斯基分