Multiple-point statistical simulation for hydrogeological models: 3-D training image development and conditioning strategies

TL;DR

This paper introduces a novel workflow for 3D multiple-point statistics (MPS) modeling in hydrogeology, emphasizing realistic training image development and data integration for large-scale, detailed geological simulations.

Contribution

It presents a new strategy for creating 3D training images and integrating diverse data sets, improving large-scale hydrogeological modeling accuracy.

Findings

Effective 3D TIs can significantly influence simulation outcomes.

Iterative TI development enhances model realism and reliability.

Combining geological and geophysical data improves model fidelity.

Abstract

Most studies on the application of geostatistical simulations based on multiple-point statistics (MPS) to hydrogeological modelling focus on relatively fine-scale models and on the estimation of facies-level structural uncertainty. Less attention is paid to the input data and the construction of Training Images (TIs). E.g. even though the TI should capture a set of spatial geological characteristics, the majority of the research still relies on 2D or quasi-3D training images. Here, we demonstrate a novel strategy for 3D MPS modelling characterized by (i) realistic 3D TIs and (ii) an effective workflow for incorporating a diverse group of geological and geophysical data sets. The study covers 2810 km^2 in southern Denmark. MPS simulations are performed on a subset of the geological succession (the lower to middle Miocene sediments) which is characterized by relatively uniform structures and dominated by sand and clay. The simulated domain is large and each of the geostatistical realizations contains approximately 45 x 10^6 voxels with size 100 m x 100 m x 5 m. Data used for the modelling include water well logs, seismic data, and a previously published 3D geological model. We apply a series of different strategies for the simulations based on data quality and develop a novel method to effectively create observed spatial trends. The TI is constructed as a relatively small 3D voxel model covering an area of 90 km^2. We use an iterative training image development strategy and find that even slight modifications in the TI create significant changes in simulations. Thus, this study shows how to include both the geological environment and the type and quality of input information in order to achieve optimal results from MPS modelling. We present a practical workflow to build the TI and effectively handle different types of input information to perform large-scale geostatistical modelling