Bayesian inversion of laboratory experiments of transport through limestone fractures

In this study, a novel experimental setup is proposed for which a column filled with glass beads and parallelepiped-shaped limestone beams is used to reconstruct a multiple fracture limestone media. The proposed setup produces asymmetric breakthrough curves (BTCs) that are consistent with the shape expected from the past field and lab-scale studies. Three transport experiments have been conducted under fast, medium, and slow flow velocity conditions. The research focuses on parameter and state estimation using Bayesian inference via Markov Chain Monte Carlo (MCMC) sampler, investigating the degree to which three models of transport through fractured media can reproduce the experimental results under the three flow conditions. The first transport model, named ADE, is based on the equivalent porous medium (EPM) approach and corresponds to the linear advection dispersion equation (ADE). The second model, named FOMIM (first-order mobile immobile), is based on the mobile/immobile approach and uses the dual porosity model with a linear first-order transfer between mobile and immobile regions. The third model, named NLMIM (non-linear mobile-immobile), uses a nonlinear transfer function between these two regions. The results of the three models show that almost all the unknown model input parameters can be well-estimated with narrow confidence intervals using the MCMC method. With respect to state estimation, the ADE model fails to reproduce correctly the tail of the BTCs observed under slow and medium flow conditions. The FOMIM model improves the tailing of the BTCs, but significant discrepancies remain between simulated and measured concentrations. The NLMIM model with velocity-dependent parameters is the only model that captures BTCs under all three conditions of slow, medium, and fast flow velocities.