3D inverse modelling of groundwater flow at a fractured site using a stochastic continuum model with multiple statistical populations

 3D groundwater flow at the fractured site of Aspö (Sweden) is simulated. The aim was to characterise the site as adequately as possible and to provide measures on the uncertainty of the estimates. A stochastic continuum model is used to simulate both groundwater flow in the major fracture planes and in the background. However, the positions of the major fracture planes are deterministically incorporated in the model and the statistical distribution of the hydraulic conductivity is modelled by the concept of multiple statistical populations; each fracture plane is an independent statistical population. Multiple equally likely realisations are built that are conditioned to geological information on the positions of the major fracture planes, hydraulic conductivity data, steady state head data and head responses to six different interference tests. The experimental information could be reproduced closely. The results of the conditioning are analysed in terms of ensemble averaged average fracture plane conductivities, the ensemble variance of average fracture plane conductivities and the statistical distribution of the hydraulic conductivity in the fracture planes. These results are evaluated after each conditioning stage. It is found that conditioning to hydraulic head data results in an increase of the hydraulic conductivity variance while the statistical distribution of log hydraulic conductivity, initially Gaussian, becomes more skewed for many of the fracture planes in most of the realisations.