论文部分内容阅读
合理截取半圆形计算区域，采取局部加密的λ-等级网格，结合对称行索引存贮格式（CSR）及并行稀疏直接求解器PARDISO，提出一种高效、高精度的2.5D直流电阻率法有限元正演方案，并编制相应的Fortran程序，对具有解析解的3个典型地电模型进行计算与分析。研究结果表明：圆形截断边界不仅便于在径向方向上采取λ-等级网格剖分，而且能大大简化有限元模拟中单元刚度矩阵的计算；结构化的等级网格避开了通常非结构化网格有限元计算时繁琐的网格剖分及总体刚度阵的集成过程，且能在不增大问题规模的前提下，显著提高2.5D直流电