Hybrid Monte Carlo 1-D joint inversion of LOTEM and MT

Abstract

Abstract In this paper two different geophysical electro-magnetic techniques, long-offset time domain EM and the magnetotelluric method, are jointly inverted using a variation of Bayesian inversion called Hybrid Monte-Carlo. Long-offset time domain EM is a method that records Earth response to a controlled source square wave signal, while the magnetotelluric method records Earth response to naturally occurring electro-magnetic signals. We obtain better understanding of subsurface resistivity structures because the information they provide is complementary and joint inversion reduces the set of permissible inversion models. In most situations where long-offset time domain EM is recorded magnetotelluric data is also collected, since receiver instrumentation is shared. Hybrid Monte Carlo uses the gradient of the forward model to enable larger jumps through the solution space than are permissible using the standard Metropolis-Hastings algorithm, avoiding the pitfall of random walk behavior and giving faster convergence. The method is successfully applied to a survey area near the edge of Wuhan China. We obtain inversion results in a highly urbanized area, where the magnetotelluric method is typically compromised by background electro-magnetic noise.