基于ILU-BICGSTAB算法的起伏地形MT二维正演模拟  

2D MT forward modeling with topography using the ILU-BICGSTAB algorithm

作　　者：郭胜男[1] 熊彬[1] 罗天涯[1] 吴延强[1] 陈欣[1] 

机构地区：[1]桂林理工大学地球科学学院,桂林541006

出　　处：《物探化探计算技术》2017年第3期306-312,共7页Computing Techniques For Geophysical and Geochemical Exploration

基　　金：国家自然科学基金(41164004;41674075);广西自然科学基金(2016GXNSFA038004;2013GXNSFAA019277);广西壮族自治区桂林市"漓江学者"(2013)专项;桂林理工大学研究生创新项目(BS201601)

摘　　要：针对起伏地形对异常体的探测效果的影响,利用有限单元法,基于MTLAB平台编写了大地电磁二维正演程序,为更好地模拟实际地形,采用矩形网格中进一步三角细化的剖分方式,并利用不完全LU分解的预条件BICGSTAB算法提高计算速度与精度,最后通过建立多种地电模型分析讨论了起伏地形对大地电磁场的影响,验证了该方法的实用性与准确性。In allusion to the issue that complex topography effects on magnetotelluric sounding, using finite element method to compute 2D magnetotelluric responses on MATLAB platform. We used rectangular grid further refines the triangulation mesh mode in order to better simulate the actual topography. The BICGSTAB algorithm with incomplete LU decomposition for preconditioning was used to improve calculating speed and precision. This approach was verified through the calculations of va-rious geoelectric model. Finally, we analysis topographic responses in magnetotelluric.

关 键 词：大地电磁 有限单元 起伏地形 正演 

分 类 号：P631.4[天文地球—地质矿产勘探]