基于非结构网格的三维大地电磁法有限内存拟牛顿反演研究

作者：曹晓月,殷长春,张博,黄鑫,刘云鹤,蔡晶

摘要：Traditional 3D Magnetotelluric(MT)forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures.By contrast,unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods,such as the high computational cost for improving model accuracy and the difficulty of inverting with topography.In this paper,we used the limited-memory quasi-Newton(L-BFGS)method with an unstructured finite-element grid to perform 3D MT inversions.This method avoids explicitly calculating Hessian matrices,which greatly reduces the memory requirements.After the fi rst iteration,the approximate inverse Hessian matrix well approximates the true one,and the Newton step(set to 1)can meet the sufficient descent condition.Only one calculation of the objective function and its gradient are needed for each iteration,which greatly improves its computational effi ciency.This approach is well-suited for large-scale 3D MT inversions.We have tested our algorithm on data with and without topography,and the results matched the real models well.We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.

发文机构：College of Geo-exploration Science and Technology Department of Earth Sciences

关键词：Magnetotelluric(MT)3DinversionUNSTRUCTUREDfinite-elementMETHODQUASI-NEWTONMETHODL-BFGS

分类号： O[理学]