胡祥云,,
彭荣华
中国地质大学(武汉)地球物理与空间信息学院, 武汉 430074
基金项目: 国家自然科学基金（41630317，41704133）和国家重点专项（2016YFC0600201-7）联合资助
详细信息
作者简介: 韩骑, 女, 1990年生, 博士研究生, 研究方向为大地电磁模拟和地热解释.E-mail:hanqi@cug.edu.cn
通讯作者: 胡祥云, 男, 1966年生, 教授, 博士生导师, 主要从事电磁法理论与应用研究.E-mail:xyhu@cug.edu.cn
中图分类号: P631收稿日期:2019-05-24
修回日期:2020-01-01
上线日期:2020-08-05
Spherical magnetotelluric modeling based on non-uniform source
HAN Qi,HU XiangYun,,
PENG RongHua
Institute of Geophysics and Geomatics, China University of Geosciences, Wuhan 430074, China
More Information
Corresponding author: HU XiangYun,E-mail:xyhu@cug.edu.cn
MSC: P631--> Received Date: 24 May 2019
Revised Date: 01 January 2020
Available Online: 05 August 2020
摘要
摘要:大地电磁测深理论与数据处理解释均假定平面电磁波垂直入射地下空间，但随着研究尺度的逐渐扩大，使得因地球弧度产生的影响难以忽略.此时，传统笛卡尔坐标体系及平面波场源不再适用于大尺度的大地电磁数据正反演解释.本文提出并实现了一种基于球坐标系的大地电磁交错网格有限差分三维正演算法，并对电场进行极向-环向分解，结合球谐函数和贝塞尔函数构建了可取代平面波的场源模型.首先利用经纬度信息构建三维地电模型，将场源设置于模型空间正上方，然后通过直接求解球坐标系下麦克斯韦方程来获得大地电磁响应.在此基础上，本文设计了球坐标下具有不同分辨率的多个三维地电模型，阐述了由球体模型到笛卡尔模型的转换方法，详细对比了两种坐标体系在计算效率、所求得的电场和视电阻率方面的差异.结果表明二者差异度主要与电性横向分布和地图投影方法有关，与周期并不存在明显的单调递增关系.
关键词: 大地电磁/
球坐标/
非均一场源/
三维模拟
Abstract:Traditionally, magnetotellurics (MT) rests on an assumption that plane wave is incident vertically into flat Earth. With more and more magnetotellurics data coming in and the expansion of research area, Cartesian approximation and plane wave assumption may no longer be suitable for the large-scale modeling due to the inevitable Earth curve. We have developed a kind of magnetotellurics forwarding code in spherical coordinates based on staggered grid finite difference method. In addition, the electric fields are represented by poloidal-toroidal decomposition, then solved based on spherical harmonic and Bessel function to construct the source model that can be used to replace plane wave. The spherical model is built directly through geographic information like latitude and longitude, and the forward response will be obtained by solving the Maxwell's equations in terms of spherical format. Electric field components defined on curving cell edges are the primary EM field that is computed iteratively utilizing a quasi-minimum residual scheme, with a divergence correction applied periodically.
Taking into account the versatility and convenience, this spherical forward code is developed on the open source ModEM which provides a ready framework for rapid implementation of new application or algorithms. Large scale three-dimensional conductive models with different resolution are built and transformed to their corresponding Cartesian models by Lambert Conformal Conic map projection which can maintain the shape of original spherical model. We compare their forward response such as electric fields and apparent resistivity computed from spherical and Cartesian coordinates respectively. Under the combined effect of conductivity distribution and the unavoidable deformation, electric fields and apparent resistivities show stronger difference at the place where conductivity contrasts. All the forward calculations are implemented on the blade server in the form of frequency parallel scheme, spherical model usually takes twice time than Cartesian. According to real bathometry data, we make the China continent and its adjacent sea model to research the forward difference between spherical and Cartesian coordinates.
Key words:Magnetotellurics/
Spherical coordinates/
Non-uniform source/
Three-dimensional modeling
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15560
