Effect of data error correlations on trans-dimensional MT Bayesian inversions

Real magnetotelluric (MT) data errors are commonly correlated, but MT inversions routinely neglect such correlations without an investigation on the impact of this simplification. This paper applies a hierarchical trans-dimensional (trans-D) Bayesian inversion to examine the effect of correlated MT data errors on the inversion for subsurface geoelectrical structures, and the model parameterization (the number of conductivity interfaces) is treated as an unknown. In the inversion considering error correlations, the data errors are parameterized by the first-order autoregressive (AR(1)) process, which is included as an unknown in the inversion. The data information itself determines the AR(1) parameter. The trans-D inversion applies the reversible-jump Markov chain Monte Carlo algorithm to sample the trans-D posterior probability density (PPD) for the model parameters, model parameterization and AR(1) parameters, accounting for the uncertainties of the model dimension and data error correlation in the uncertainty estimates of the conductivity profile. In the inversion ignoring the correlation, we neglect the correlation effect by turning off the AR(1) parameter. Then the correlation effect on the MT inversion can be examined upon comparing the posterior marginal conductivity profiles from the two inversions. Further investigation is then carried out for a synthetic case and a real MT data example. The results indicate that for strong correlation cases, neglecting error correlations can significantly affect the inversion results.

Large-scale 3D inversion of marine magnetotelluric data: Case study from the Gemini prospect, Gulf of Mexico

Three-dimensional magnetotelluric (MT) inversion is an emerging technique for offshore hydrocarbon exploration. We have developed a new approach to the 3D inversion of MT data, based on the integral equation method. The Tikhonov regularization and physical constraint have been used to obtain a stable and reasonable solution of the inverse problem. The method is implemented in a fully parallel computer code. We have applied the developed method and software for the inversion of marine MT data collected by the Scripps Institution of Oceanography (SIO) in the Gemini prospect, Gulf of Mexico. The inversion domain was discretized into 1.6 million cells. It took nine hours to complete 51 iterations on the 832-processor cluster with a final misfit between the observed and predicted data of 6.2%. The inversion results reveal a resistive salt structure, which is confirmed by a comparison with the seismic data. These inversion results demonstrate that resistive geoelectrical structures like salt domes can be mapped with reasonable accuracy using the 3D inversion of marine MT data.