复杂电性结构大地电磁二维响应特征  

Magnetotelluric Responses of Two-Dimensional Complex Conductivity Structures

作　　者：罗天涯[1] 熊彬[1] 蔡红柱 陈欣[1] 刘云龙[1] 兰怀慷 李祖强[1] 梁卓[1] Luo Tianya Xiong Bin Cai Hongzhu Chen Xin Liu Yunlong Lan Huaikang Li Zuqiang Liang Zhuo(College of Earth Sciences, Guilin University of Technology, Guilin 541006, Guangxi ,China College of Mines ＆ Earth Sciences, University of Utah, Salt Lake City 84112, UT , USA)

机构地区：[1]桂林理工大学地球科学学院,广西桂林541006 [2]犹他大学矿产和地球科学学院,美国犹他州盐湖城84112

出　　处：《吉林大学学报（地球科学版）》2017年第1期215-223,共9页Journal of Jilin University：Earth Science Edition

基　　金：国家自然科学基金项目(40974077;41164004;41674075);广西自然科学基金项目(2016GXNSFA380004;2013GXNSFAA019277);桂林市"漓江学者"专项经费资助项目(2013005);桂林理工大学研究生创新项目(BS201601)~~

摘　　要：为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。To improve the discretization technique,adaptive unstructured quadrilateral grid finite element（FE）procedure is presented for solving weighted residual formulation of magnetotelluric（MT）fields in 2-D conductivity structures.The finer the meshes are,the more accurate the numericalsolutions are.However,this is at the cost of time-consuming and large memory requirement.Adaptive mesh refinement generally provides an optimization of both run time and accuracy since in the method meshes are only refined where the electrical parameters vary intensely,while coarsen in the slow electrical resistivity change region.Due to the fact that regular grids have no potential to incorporate arbitrary topography and complex geologic structures,we introduced COMSOL Multiphysics software package to build arbitrary geometries and subdivide them with free quadrilaterals,then apply meshes information to our algorithm.Comparison with analytical solutions and previous numerical solutions demonstrate a high-accuracy and flexibility of numerical simulations based on unstructured grid,a new technology to discretize the computating domain.

关 键 词：非结构化网格 大地电磁 有限单元法 正演 

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