多重网格法在CSAMT三维正演中的计算效率分析.pdfVIP

复数线性方程组阍。而多重网格算法是已被证明了的近乎最优的求解大型线性方程组的方法哺1，该方法已在直流 电阻率三维数值模拟中取得了较好的效果H1。因此本文采用多重网格算法实现基于交错采样有限差分法的CSA硼I 三维数值模拟。通过与两种传统迭代方法比较证明多重网格算法能够大大提高CSAgT三维数值模拟的计算速度。 l基于交错网格采用有限差分的CsAMT三维数值模拟 根据可控源音频大地电磁法的电磁场响应特征，将电磁场分解成一次场和二次场的叠加。其中一次场一般为 在均匀半空间(或水平层状介质)下有限长电偶源的电磁响应，可以利用快速汉克尔变换实现快速、准确的求取 全空间电磁响应。而二次场方程根据总场的MAxVWELL方程组减去一次场的MAXWELL方程组得到，其表达式为： f秘4棚=护4·豳 l和4棚=肌棚4·豳 次场。 首先对求解区域在直角坐标系中采用三组平行x，Y，Z坐标轴平面剖分，对于每个六面体单元采用YEE单元(交 错网格采样)定义方式。对式(1)中进行离散可以得到关于磁场和电磁二次场离散方程，消去电场二次场电场分 量可以得到关于磁场二次场三分量的离散线性方程组。同时在分析磁场三个分量衰减特性的基础上，提出了研究 区域足够大时，取各边界的二次场值为零的简洁边界条件，降低了CSAMT三维数值模拟的难度。 2 CSAMT三维数值模拟的多重网格法 多重网格法是求解偏微分方程离散化所得的大型线性方程组的最有效的方法之一，基本原理为通过把迭代分 散在许多不同间距的网格中进行，不同频率的误差分量在不同间距的网格中被削减掉从而实现快速求解。因此在 CSAMT三维数值模拟中需要研究多重网格法的粗化策略，限制算子、延拓算子、光滑迭代等关键技术。 本论文中采用标准粗化策略(在各个剖分方向上剖分数减半)对求解单元进行粗化：在各层网格上采用 Guass-Seidel迭代方法对离散后得到的关于磁场的大型线性方程组进行光滑迭代，并且通过V循环的方式将各层 网格上的求解残差通过限制和延拓方式实现粗网格校正。 由于插值和限制算子是互逆的。可以通过矩阵的互逆运算可以实现相互转换。限制算子通过加权平均的方式 推导了三个方向上的磁场限制过程的表达式，其x方向的表达式为： ● 2， 2， 2上 (2) H，(I’J，K)=÷∑艺∑W：w二w：H沁m，n) 1 1=2f—I_一2 J一2H一2t一2 其中： I垒纽 ，；2i一1 ，纷=2j一2 些±=2 刀：2k一2 △22¨+△z2I一2 ， 吖：{她≮+蝇， ’矿： 辨=2j-I嵋= l 刀=2k—l ’I 坚f f-2i ”_一 【△也卜I+蝇， ，玎=2j 竺±=! 刀：2七 ％一一。‰一一 Az2t—I+Az2t一2 3 CSAgT维正演实例及分析 为验证CSAgT三维数值模拟中多重网格法(啪)的优势，从计算时间、误差衰减速度、迭代次数等角度比较 了№法与两种常用的迭代方法(不完全LU分解共轭梯度法(ILU—cG)和不完全LU分解双共轭梯度迭代算法 (ILu—BIcG))的优劣。所用模型、网格剖分及计算频率均与散度校正研究时的相同。在采用RRCM散度校正的基 础上分别采用ILU-CG、ILU—BIC6、gG法进行迭代求解，淞法在三个方向上的最小网格单元为2，多重网格层数 为5，采用V循环方式进行迭代。 曲线。由图可以看出ILU-CG迭代的前lO次迭代误差曲线出现震荡。但第lO次