罗威^1,2,,
王绪本^1,,,
王堃鹏¹,
张刚³,
李德伟¹
1. 成都理工大学地球物理学院, 成都 610059
2. 四川省蜀通岩土工程公司, 成都 610084
3. 西南科技大学环境与资源学院, 四川绵阳 621010

基金项目: 国家高技术研究发展计划项目（2014AA06A612）和国家自然科学基金项目（41674078，41704105）联合赞助


详细信息

作者简介: 罗威, 男, 1988年生, 博士研究生, 主要从事电磁法的正反演研究.E-mail:luovitas@qq.com

通讯作者: 王绪本, 男, 1956年生, 教授, 博士生导师, 主要从事地球物理电磁勘探理论的研究.E-mail:wxb@cdut.edu.cn

中图分类号: P631

收稿日期:2018-07-15
修回日期:2019-08-01
上线日期:2019-10-05

Three-dimensional forward modeling of the magnetotelluric method in spherical coordinates

LUO Wei^1,2,,
WANG XuBen^1,,,
WANG KunPeng¹,
ZHANG Gang³,
LI DeWei¹
1. College of Geophysics, Chengdu University of Technology, Chengdu 610059, China
2. Sichuan Shutong Geotechnical Engineering Company, Chengdu 610084, China
3. School of Environment and Resource, Southwest University of Science and Technology, Mianyang Sichuan 621010, China


More Information

Corresponding author: WANG XuBen,E-mail:wxb@cdut.edu.cn

MSC: P631

--> Received Date: 15 July 2018
Revised Date: 01 August 2019
Available Online: 05 October 2019


摘要
摘要:大地电磁正演理论研究热点一直以来主要集中在如何提高计算效率和精度，但在剖面足够长、探测深度足够大的情况下，传统的笛卡尔坐标系数值模拟方式难以准确拟合地球曲率形态.本文研究了基于球坐标系的三维大地电磁正演，推导了交错网格有限差分三维正演公式，与一维解析解和三维标准模型测试对比，验证了正演算法的正确性.通过理论模型计算，对比分析球坐标和笛卡尔坐标系正演结果表明：球坐标系模拟更合理，避免了传统笛卡尔坐标拉伸投影所引入的误差，可代替目前的笛卡尔坐标模拟方法.基于球坐标和笛卡尔坐标系的三维大地电磁正演响应值随着频率变低差异越明显.球坐标和笛卡尔坐标计算结果差异度与频率、模型结构和电阻率有关.本文模型计算结果在数万秒周期处已出现接近10%的差异，对于较大尺度的长周期大地电磁，地球曲率的影响不能忽略.
关键词: 三维大地电磁/
球坐标/
地球曲率/
有限差分

Abstract:Theoretical research on forward modeling of the magnetotelluric (MT) method has focused on how to improve the calculation efficiency and precision. When the profile is long enough and detection depth is sufficiently deep, it is difficult to fit the curvature of the earth accurately by numerical simulation in a Cartesian coordinate system. This paper studies three-dimensional MT forward modeling with a staggered grid based in a spherical coordinate system. We deduce the formula for staggered-grid finite difference three-dimensional forward modeling and compare it with of one-dimensional analytical solution and the DTM1 three-dimensional modeling to verify correctness of the algorithm. Calculation on a theoretical model shows that the modeling in the spherical coordinates is more reasonable, which can avoid errors caused by the conventional Cartesian coordinate system, so can replace the Cartesian coordinates modeling method. The differences of 3D MT forward modeling responses in theses two coordinate systems are related with frequency, model structure and resistivity. For example, such differences become more notable with lowering frequency. They are close to 10% at the period of tens of thousands of seconds. It implies that the effect of earth curvature on large-scale long-period MT cannot be ignored, so MT forward modeling on large scales should be conducted in spherical coordinates.
Key words:3-D magnetotelluric/
Spherical coordinates/
Earth curvature/
Finite difference



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15196