CSAMT三维交错网格有限差分数值模拟中的散度校正方法研究.pdfVIP

第十届中国国际地球电磁学术讨论会201l f叮Ⅳ。·d／=伊。·dS (1) bndr=J弛棚。·岱 次场。 首先对求解区域在直角坐标系中采用三组平行X，Y，Z坐标轴平面剖分，对于每个六面体单元采用YEE单元(交 错网格采样)定义方式。对式(1)中进行离散可以得到关于磁场和电磁二次场离散方程，消去电场二次场电场分 量可以得到关于磁场二次场三分量的离散线性方程组。同时在分析磁场三个分量衰减特性的基础上，提出了研究 区域足够大时，取各边界的二次场值为零的简洁边界条件，降低了CSA-iT三维数值模拟的难度。 2磁场散度校正 由于麦克斯韦方程组磁场散度为零，因此总场和背景场的磁场强度散度也应为零，即： V·H=O (2) V．H‘=0 (3) 两式相减得： V·(Ⅳ·H6)=V·H4=O (4) 即二次场的磁场强度散度也为零。 然而，在CS／WT三维数值模拟过程中由于计算机的截断误差、迭代求解的误差积累、模型电阻率的突变、 不等距剖分等因素的影响，很难保证每个求解单元内的磁场散度为零，从而影响计算结果的准确性和精度。为避 免上述现象发生，在CS／u{T三维数值模拟求解过程中常采用磁场散度校正的方法对迭代求得的磁场值进行修正。 本章从磁场散度为零的等式出发，提出了直接磁场散度校正(DDcM)、磁场散度残差复校正(RPcM)和磁场散度 残差实校正(RRCM)三种磁场散度校正方法。在详细讨论磁场散度校正技术的基础上，通过对比确定了一种适用 于CSAMT三维数值模拟的最佳磁场散度校正方法。 3散度校正计算效率分析 (NDCM)的计算效率、适用性及优缺点，设计了一个典型低阻异常体模型，从数值模拟的计算速度、迭代次数、 误差衰减曲线等进行比较。模型中低阻异常体大小为500米×500米×500米，电阻率为1Q·m的，顶部埋深为 200米，围岩电阻率为100Q·m。将求解区域剖分成32×32×24个网格单元，采用指数延拓增加8层空气层，计 算频率为lOHz。采用ILU-BICG方法求解离散后形成的线性方程组。 图中可以看出采用NDCM方案时，开始20次迭代能够较快收敛，但第20次以后的迭代，收敛速度明显变慢，第 x10—6)，这说 50次迭代后，误差曲线出现震荡，并且最大精度只能达到1×10-4，未能达到设定的计算精度(1 明NDCM求解方案不收敛。造成这种现象的主要原因是散度校正公式破坏了迭代求解的对称性，从而使方程求解 发散。RPcM和RRCM方案误差衰减曲线前半段几乎一样。这是因为两种校正方法的原理相同。但是，由于RPCM 求解复数型磁场散度方程，相对于RRCM迭代增加了虚部的计算机截断误差积累，因此RRCM比RPCM的迭代次数 少，收敛速度更快。表1为采用不同散度校正方法进行数值模拟所用的计算时间(求二次场时解线性方程组的时 次迭代时间最短，为0．365秒，而RPCM单次迭代时间最长，这是因为NDCM单次迭代不需要求解磁场散度残差方 程，而RPCM需要求解复数型的磁场散度残差，相对于RRCi4求解实数型的残差方案计算量更大。 从上面的分析可以看出，在不破坏方程求解对称性的情况下，散度校正可以明显提高计算速度和精度，加快 收敛速度。在散度校正中，RRCM方案相对于RPCM方案计算量更小、收敛速度更快、迭代次数更少。因此在CSAMT 三维数值模拟中应采用RRcM的散度校正方法。 为了验证RRCM散度校正方案的适用性和稳定性，对不同网格剖分和不同频点数时RRCM方案误差衰减曲线和 迭代时间进行了比较。共设计了8×8×6、16×16X 数值模拟的计算实例，仍用前面的典型低阻异常模型。图2为不同剖分数下采用RRCh{散度校正时数值模拟的误 19 ———————卫!堂堕!!坐堕竺!!!竺堂堂!!罂!：!燮业兰!! 差衰减曲线，从圈中rⅡ以看出-四种剖分方案计算误*部能快速