Geostatistical Simulation of Hydrofacies Heterogeneity of the West Thessaly Aquifer Systems in Greece

Integrating geological properties, such as relative positions and proportions of different hydrofacies, is of highest importance in order to render realistic geological patterns. Sequential indicator simulation (SIS) and Plurigaussian simulation (PS) are alternative methods for conceptual and deterministic modeling for the characterization of hydrofacies distribution. In this work, we studied the spatial differentiation of hydrofacies in the alluvial aquifer system of West Thessaly basin in Greece. For this, we applied both SIS and PS techniques to an extensive set of borehole data from that basin. Histograms of model versus experimental hydrofacies proportions and indicative cross sections were plotted in order to validate the results. The PS technique was shown to be more effective in reproducing the spatial characteristics of the different hydrofacies and their distribution across the study area. In addition, the permeability differentiations reflected in the PS model are in accordance to known heterogeneities of the aquifer capacity.