| 大地电磁法三维交错采样有限差分数值模拟 |
| |
| 引用本文： | 谭捍东,余钦范,John Booker,魏文博.大地电磁法三维交错采样有限差分数值模拟[J].地球物理学报,2003,46(5):705-711. |
| |
| 作者姓名： | 谭捍东 余钦范 John Booker 魏文博 |
| |
| 作者单位： | 1.中国地质大学地球物理与信息技术学院，北京 100083 2 中国科学院地质与地球物理研究所，北京 100029 3 University of Washington，Seattle，USA |
| |
| 基金项目： | 国家自然科学基金项目 (4 0 0 740 19)，教育部重大项目 (重大 0 2 11)资助 |
| |
| 摘 要： | 系统地论述了大地电磁三维交错采样有限差分数值模拟算法实现过程中交错网格剖分、积分公式离散化、边界条件、方程组求解、三维张量阻抗的计算等内容. 由于提出了简洁的边界条件，采用了解大型系数矩阵方程组的双共轭梯度稳定解法，所实现的三维交错采样有限差分数值模拟算法具有迭代收敛稳定、计算精度高、速度快等特点. 通过两个理论模型的计算结果检验了算法的正确性和计算精度. 所实现的三维交错采样有限差分数值模拟算法为研究三维反演问题奠定了基础.
|
| 关 键 词： | 大地电磁法 交错采样有限差分法 三维正演 边界条件 双共轭梯度稳定解法 |
| 文章编号： | 0001-5733(2003)05-0705-07 |
| 收稿时间： | 2002-12-19 |
| 修稿时间： | 2003-5-27 |
MAGNETOTELLURIC THREE-DIMENSIONAL MODELING USING THE STAGGERED-GRID FINITE DIFFERENCE METHOD |
| |
| TAN HANDONG.MAGNETOTELLURIC THREE-DIMENSIONAL MODELING USING THE STAGGERED-GRID FINITE DIFFERENCE METHOD[J].Chinese Journal of Geophysics,2003,46(5):705-711. |
| |
| Authors: | TAN HANDONG |
| |
| Institution: | 1.School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China 2 Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China 3 University of Washington, Seattle, USA |
| |
| Abstract: | The crucial problems of 3 T5BZ]D forward modeling using the staggered-grid finite difference method are described in detail in this paper. They are staggered-grid, discretization of integrated form of Maxwell equation, boundary condition, solving linear algebra equations and calculating D forward modeling using the staggered-grid finite difference method are described in detail in this paper. They are staggered-grid, discretization of integrated form of Maxwell equation, boundary condition, solving linear algebra equations and calculating 3 T5BZ]D tensor impedance. Giving the more explicit boundary conditions and using Bi-conjugate gradients stabilized method to solve the linear algebra equation with large coefficient matrix, we get a fast algorithm of high-precision to calculate effectively electrical and magnetic fields in the whole space. This has been proved by comparing the 3D forward modeling solutions with analytic solution to abutting quarter-spaces, and 2D forward modeling to 2D prism using the 2D finite element method. This efficient algorithm has setup basis for research of D tensor impedance. Giving the more explicit boundary conditions and using Bi-conjugate gradients stabilized method to solve the linear algebra equation with large coefficient matrix, we get a fast algorithm of high-precision to calculate effectively electrical and magnetic fields in the whole space. This has been proved by comparing the 3D forward modeling solutions with analytic solution to abutting quarter-spaces, and 2D forward modeling to 2D prism using the 2D finite element method. This efficient algorithm has setup basis for research of 3D inversion . |
| |
| Keywords: | Magnetotellurics Staggered-grid finite difference method Three-dimensional forward modeling Boundary condition Bi-conjugate gradients stabilized method |
| 本文献已被 CNKI 维普 万方数据 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
