科研之友
微信
新浪微博
Facebook
分享链接
摘要
地球物理勘探中，地下复杂地质体电磁场三维数值模拟计算量大、存储要求高.针对这一问题，本文提出一种空间-波数域三维电磁场数值模拟方法，该方法利用电磁场积分方程为卷积的特点，将积分沿水平方向进行二维傅里叶变换，将电磁场三维空间域卷积问题转换为不同波数相对独立的多个一维积分问题，由此计算量和存储需求大大减少，易于实现并行计算.采用有限单元法中的形函数进行一维数值积分，一维积分离散为多个单元积分之和，每个单元采用二次形函数表征电流变化，可得出单元积分的解析表达式；保留垂向为空间域，优势之一在于可根据实际情况合理调整单元疏密程度，准确模拟任意复杂地形和导电率异常，兼顾计算精度与计算效率，优势之二是用形函数拟合求得积分的解析解，计算精度和效率高；最后引入压缩算子，采用迭代求得电磁场的数值解.本文方法充分利用不同波数之间一维积分高度并行性、一维形函数积分的高精度及快速傅里叶变换的高效性，实现电磁场三维高效高精度数值模拟.设计模型将本文方法数值模拟结果和软件INTEM3D的数值模拟结果对比，验证了方法的正确性；设计高阻和低阻异常体，研究了异常电导率与背景电导率差异对迭代收敛速度的影响；改变计算规模，随着计算网格增多，算法耗时与存储呈近似线性增长；设计复杂模型与目前主流数值模拟方法对比，本文算法速度快1个数量级以上，且计算规模越大，算法效率的优势越明显.研究结果表明，本文提出的空间-波数域三维电磁场数值模拟方法理论和方法正确，计算效率高，对计算机存储要求低，算法高度并行，适合任意复杂条件大规模三维电磁场高效、高精度数值模拟.
- 单位
