| 基于矢量有限元的带地形大地电磁三维反演研究 |
| |
| 引用本文： | 顾观文,李桐林.基于矢量有限元的带地形大地电磁三维反演研究[J].地球物理学报,2020,63(6):2449-2465. |
| |
| 作者姓名： | 顾观文 李桐林 |
| |
| 作者单位： | 1. 吉林大学, 地球探测科学与技术学院, 长春 130026;2. 防灾科技学院, 地球科学学院, 三河 065201;3. 河北省地震动力学重点实验室, 三河 065201 |
| |
| 基金项目： | 中央高校创新团队项目（ZY20180102）和自然资源部"十五"重点科技项目（20010211）联合资助. |
| |
| 摘 要： | 研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.
|
| 关 键 词： | 大地电磁 矢量有限元法 PARDISO 三维反演 带地形 |
| 收稿时间： | 2019-08-02 |
Three-dimensional magnetotelluric inversion with surface topography based on the vector finite element method |
| |
| GU GuanWen,LI TongLin.Three-dimensional magnetotelluric inversion with surface topography based on the vector finite element method[J].Chinese Journal of Geophysics,2020,63(6):2449-2465. |
| |
| Authors: | GU GuanWen LI TongLin |
| |
| Institution: | 1. College of Geo-Exploration Sciences and Technology, Jinlin University, Changchun 130026, China;2. Schoolof Earth Sciences, Instituteof Disaster Prevention, Sanhe 065201, China;3. Hebei Key Laboratory of Earthquake Dynamics, Sanhe 065201, China |
| |
| Abstract: | In this paper, the three-dimension (3D) Magnetotelluric (MT) inversion algorithm with topography based on the vector finite element method is studied and the 3D inversion program code is developed. In the 3D forward numerical modeling of the MT field, a fast forward modeling scheme based on Parallel Direct Sparse Solver (PARDISO) without divergence correction is adopted. For typical models with topography, under the conditions of medium-scale calculations, comparison shows that the PARDISO method is more than 10 times faster in calculation speed than Bi-Conjugate Gradient (BICG). The correctness of the 3D forward modeling program with topography is also confirmed by tests on a theoretical model and the comparison with the results in previous studies using the finite element method. In 3D inversion, an inversion code of MT with topography based on the Conjugate Gradient (CG) method is compiled in this study to avoid directly calculation of the Jacobian matrix, thus the problem existing in calculation of the Jacobian matrix is turned into solving two "quasi-forward" problems, and the fast forward solution of PARDISO is further applied to the solution of the "quasi-forward" modeling to improve the efficiency of inversion calculation. The developed algorithm is used to carry out 3D inversion of synthetic data of a series of geo-electric models with topography. The inversion results can well reproduce the electrical structure of the theoretical model, which verifies the correctness and reliability of the 3D inversion algorithm developed in this paper. Finally, the algorithm is used to invert the MT measurement data in a mine area. The resultant 3D electrical structure from the inversion clearly reveals geoelectric characteristics of the study area. Data interpretation is conducted successfully by combining the inversion results with geological data available, further verifying the effectiveness of the algorithm proposed in this work. |
| |
| Keywords: | Magnetotelluric Vector finite element method Parallel Direct Sparse Solver (PARDISO) 3D Inversion Topography |
| 本文献已被 CNKI 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
