Hydrogeophysical data integration through Bayesian Sequential Simulation with log-linear pooling. Issue 3 (6th February 2020)

• Record Type: Journal Article
• Authors: Nussbaumer, Raphaël; Mariethoz, Grégoire; Gloaguen, Erwan; Holliger, Klaus
• Is Part Of: Geophysical journal international. Volume 221:Issue 3(2020)
• Page Start: 2184
• Page End: 2200
• Publication Date: 2020-02-06
• Subjects: Hydrogeophysics -- Electrical resistivity tomography (ERT) -- Statistical methods; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggaa072
• Languages: English
• ISSNs: 0956-540X

Abstract

SUMMARY: Bayesian sequential simulation (BSS) is a geostastistical technique, which uses a secondary variable to guide the stochastic simulation of a primary variable. As such, BSS has proven significant promise for the integration of disparate hydrogeophysical data sets characterized by vastly differing spatial coverage and resolution of the primary and secondary variables. An inherent limitation of BSS is its tendency to underestimate the variance of the simulated fields due to the smooth nature of the secondary variable. Indeed, in its classical form, the method is unable to account for this smoothness because it assumes independence of the secondary variable with regard to neighbouring values of the primary variable. To overcome this limitation, we have modified the Bayesian updating with a log-linear pooling approach, which allows us to account for the inherent interdependence between the primary and the secondary variables by adding exponential weights to the corresponding probabilities. The proposed method is tested on a pertinent synthetic hydrogeophysical data set consisting of surface-based electrical resistivity tomography (ERT) data and local borehole measurements of the hydraulic conductivity. Our results show that, compared to classical BSS, the proposed log-linear pooling method using equal constant weights for the primary and secondary variables enhances the reproduction of the spatial statistics of the stochastic realizations, while maintaining a faithful correspondence with the geophysical data. Significant additional improvements can be achieved by optimizing the choice of these constant weights. We also explore a dynamic adaptation of the weights during the course of the simulation process, which provides valuable insights into the optimal parametrization of the proposed log-linear pooling approach. The results corroborate the strategy of selectively emphasizing the probabilities of the secondary and primary variables at the very beginning and for the remainder of the simulation process, respectively.