The Bayesian method is a powerful tool to estimate the resistivity distribution and associate uncertainty from time-domain electromagnetic (TDEM) data. As the forward simulation of the TDEM method is computationally expensive and a large number of samples are needed to globally explore the model space, the full Bayesian inversion of TDEM data is limited to layered models. To make high-dimensional Bayesian inversion tractable, we propose a divide-and-conquer strategy to speed up the Bayesian inversion of TDEM data. First, the full datasets and model spaces are divided into disjoint batches based on the coverage of the sources so that independent and highly efficient Bayesian subsampling can be conducted. Then, the samples from each subsampling procedure are combined to get the full posterior. To obtain an asymptotically unbiased approximation to the full posterior, a kernel density product method is used to reintegrate samples from each subposterior. The model parameters and their uncertainty are estimated from the full posterior. The proposed method is tested on synthetic examples and applied to a field dataset acquired with a large fixed-loop configuration. The 2D section from the Bayesian inversion revealed several mineralized zones, one of which matches well with the information from a nearby drill hole. The field example shows the ability of Bayesian inversion to infer reliable resistivity and uncertainty.
