Time-lapse inversion of one-dimensional magnetotelluric data

Abstract

We present a new tool for modelling time-lapse magnetotelluric (MT) data, an emerging technique for monitoring changes in subsurface electrical resistivity. Time-lapse MT data have been acquired in various settings, including sites of hydraulic fracturing, dewatering and sequestration. It has been shown in other geophysical techniques that the most effective way to model time-lapse data is with simultaneous inversion, which uses information from all timesteps to produce models with higher accuracy and fewer artefacts. We introduce this method to model time-lapse 1D MT data. As with a standard MT inversion, our routine penalises spatial roughness at each time-step, however we also introduce temporal regularisation. The inversion is simple to apply, requiring only the ratio between regularisation parameters and the desired level of misfit from the user. The algorithm is tested on both synthetic data, and a case study. We find that in the synthetic example our inversion successfully retrieves the main characteristics of the test model and introduces minimal artefacts, even in the presence of significant noise. We also test the effect of changing the ratio of regularisation parameters. In the case study, we produce an easily interpretable model that compares favourably with previous inversions of the synthetic data. We conclude that time-lapse modelling of 1D MT data can be a valuable tool for imaging subsurface change.

One‐Dimensionl Inversion of Magnetotelluric Field Data by SIS and Occam's Inversion Technique: A Case Study from the Indian Himalaya

“…The modeling and inversion for 1D structure are relatively simple and routinely used by numerous MT practitioners. [5][6][7][8][9][10][11][12] 1D inversion of MT data is also used to construct an initial guess model for 2D and 3D inversion. [13][14][15] Linearized least square inversion and minimum norm inversion are the basic approaches for 1D inversion.…”