Simulation of non-Gaussian transmissivity fields honoring piezometric data and integrating soft and secondary information

The conditional probabilities (CP) method implements a new procedure for the generation of transmissivity fields conditional to piezometric head data capable to sample nonmulti-Gaussian random functions and to integrate soft and secondary information. The CP method combines the advantages of the self-calibrated (SC) method with probability fields to circumvent some of the drawbacks of the SC method-namely, its difficulty to integrate soft and secondary information or to generate non-Gaussian fields. The SC method is based on the perturbation of a seed transmissivity field already conditional to transmissivity and secondary data, with the perturbation being function of the transmissivity variogram. The CP method is also based on the perturbation of a seed field; however, the perturbation is made function of the full transmissivity bivariate distribution and of the correlation to the secondary data. The two methods are applied to a sample of an exhaustive non-Gaussian data set of natural origin to demonstrate the interest of using a simulation method that is capable to model the spatial patterns of transmissivity variability beyond the variogram. A comparison of the probabilistic predictions of convective transport derived from a Monte Carlo exercise using both methods demonstrates the superiority of the CP method when the underlying spatial variability is non-Gaussian.