Estimating porosity distribution of a heterogeneous alluvial aquifer by GPU-accelerated 3D conditional stochastic inversion of common-offset GPR reflection data

Three-dimensional (3D) high-resolution and accurate characterization of the heterogeneity in an alluvial aquifer is critical to predicting the groundwater flow and contaminant transport. In this study, we extend the previous 2D conditional stochastic inversion method with respect to common-offset ground-penetrating radar (GPR) data to 3D for the estimation of 3D porosity distribution in a heterogeneous alluvial aquifer. First, 3D stochastic re-alizations of subsurface properties can be obtained using the fast Fourier transform moving average (FFT-MA) method, and the realizations are conditioned to borehole porosity measurements available in the 3D survey area, as well as to the geostatistical parameters derived from the borehole logs and the processed 3D GPR data. Next, the realizations are constantly modified and regenerated via a localized simulated annealing optimization strategy, which ultimately makes their corresponding synthetic data offer an acceptable fit to the GPR data. To accelerate the computational speed of 3D inversion, we design and implement the multi-GPUs multi-node par-allelization of the 3D conditional stochastic inversion algorithm. The proposed 3D inversion algorithm is posi-tively verified through the application to field GPR data collected at the Boise Hydrogeophysical Research Site near Boise, Idaho, USA. Our results indicate that the proposed 3D inversion method has the potential to effec-tively recover the fine-scale porosity distribution of a heterogeneous alluvial aquifer.