Three-dimensional magnetotelluric inversion under topographic relief based on the limited-memory quasi-Newton algorithm (L-BFGS)

Yu Hui,Deng JuZhi,Chen Hui,Chen Xiao,Wang XianXiang,Zhang ZhiYong,Ye YiXin,Chen ShuNi
DOI: https://doi.org/10.6038/cjg2019M0258
2019-01-01
Abstract:We present an algorithm of three-dimensional (3D) magnetotelluric (MT) inversion under topographic relief using a limited-memory quasi-Newton optimization based on Broyden-Fletcher-Goldfarb-Shanno (BFGS) formula in this study. Firstly, we establish the approximate expression and calculation of the inverse of Hessian matrix for 3D MT inversion based on iterative minimization of a classical Tikhonov regularized penalty function. A covariance matrix is designed to keep the resistivity of air unchanged and to impose the smoothness of the model while implementing the inversion with topographic relief. In addition, the cooling method for the regularization factor and a strategy to search for the step length based on the Wolf condition are adopted in this algorithm to improve the stability of the inversion. Secondly, we make a 3D inversion of synthetic data which are generated utilizing a series of conductivity structure models with topographic relief and compare the results with that of the well-known 3D MT inversion program (ModEM) to verify the correctness and reliability of the algorithm in this paper. Finally, we use this 3D algorithm to invert observed MT data with topographic relief from a mountainous survey area in South China, yielding an electrical structure model. It can be roughly divided into three layers, i. e. a high-resistivity basement, a conductive layer and the unconformity surface in the middle and the sedimentary cover with high-resistivity. The applicability of this algorithm is further verified.
What problem does this paper attempt to address?