检索规则说明:AND代表“并且”;OR代表“或者”;NOT代表“不包含”;(注意必须大写,运算符两边需空一格)
检 索 范 例 :范例一: (K=图书馆学 OR K=情报学) AND A=范并思 范例二:J=计算机应用与软件 AND (U=C++ OR U=Basic) NOT M=Visual
机构地区:[1]海底科学与探测技术教育部重点实验室中国海洋大学海洋地球科学学院,青岛266100 [2]青岛海洋科学与技术国家实验室海洋矿产资源评价与探测技术功能实验室
出 处:《地球物理学报》2017年第12期4874-4886,共13页Chinese Journal of Geophysics
基 金:国家自然科学基金重点项目(41130420);国家自然科学青年基金项目(41604063)联合资助
摘 要:本文实现了可控源电磁(CSEM)场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性:一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.We present a 3D finite-element (FE) algorithm for the CSEM forward problem based on the secondary field formulation. In this method, the electric fields are approximated by edge elements rather than nodal elements, so that spurious solutions can be eliminated. The computational domain is subdivided into unstructured tetrahedral elements, which allows using complex model geometries including bathymetry and tilted electrical interfaces. To further improve the accuracy of the FE solution, adaptive grid refinement guided by an a posteriori error estimator is performed iteratively.To demonstrate the performance of the presented algorithm, we first benchmark FE solutions against analytical solutions for the 1D canonical model, and show that the errors of FE solutions decrease rapidly at the first few grid refinement iterations, and after 10 refinement iterations a highly accurate solution is achieved, with the relative amplitude error of electromagnetic fields close to 1 percent and the error in phase less than 1°. Then we compare FE solutions and finite volume (FV) solutions for the 3D canonical model, yielding an excellent agreement with each other. Finally, forward responses of a slightly complicated model in which a 3D sloping conductor is presented is adequately modeled by the FE algorithm proposed in this paper.The synthetic numerical examples indicate that the presented FE algorithm is very effective, well suitable to handle complex models with unstructured meshes employed, and highly precise modeling results can be obtained by performing the adaptive grid refinement.
关 键 词:可控源电磁 三维正演 自适应矢量有限元 非结构网格
分 类 号:P631[天文地球—地质矿产勘探]
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在载入数据...
正在链接到云南高校图书馆文献保障联盟下载...
云南高校图书馆联盟文献共享服务平台 版权所有©
您的IP:216.73.216.200