Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

Abstract

Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.

Résumé

Les processus de sédimentation génèrent fréquemment des hétérogénéités tridimensionnelles dans les paramètres hydrogéologiques en subsurface. Les approches stochastiques basées sur des variogrammes ne sont souvent pas en mesure de décrire cette hétérogénéité de manière satisfaisante dans des environnements géologiques si complexes. La présente étude expose comment des analyses géostatistiques à points multiples peuvent être mises en œuvre pour une application hydrogéologique réaliste, afin de déterminer l’impact des hétérogénéités géologiques complexes sur l’écoulement et le transport souterrains. Cette approche est mise en application sur un aquifère de Belgique, qui démontre une hétérogénéité sédimentaire et une anisotropie complexes. Une image expérimentale est construite sur la base des données géologiques et hydrogéologiques de terrain. Les statistiques à points multiples sont empruntées à cette image afin de simuler l’occurrence des faciès perméables, tandis que la variabilité de la perméabilité entre les faciès est simulée en utilisant les méthodes géostatistiques standard, basées sur les variogrammes. Les transmissivités simulées sont utilisées comme données entrantes dans un modèle d’écoulement et de transport souterrain, afin de simuler l’effet des hétérogénéités sédimentaires à petite échelle sur la migration d’un panache de polluant. Les résultats montrent que les hétérogénéités à petite échelle ont un effet significatif sur le transport de polluant dans l’aquifère étudié. L’incertitude sur la distribution spatiale des faciès et des transmissivités intrafaciès génère une incertitude significative sur la distribution calculée des concentrations. La comparaison avec les techniques standard basées sur les variogrammes démontre que les analyses géostatistiques à points multiples permettent une reproduction plus fidèle des horizons argileux peu perméables et discontinus qui influencent le transport de solutés.

Resumen

Los procesos sedimentológicos a menudo dan como resultado heterogeneidad tridimensional compleja en el subsuelo de los valores de los parámetros hidrogeológicos. Aproximaciones estocásticas basadas en variogramas no son a menudo capaces de describir la heterogeneidad en tales ambientes geológicos complejos. Este trabajo muestra de que manera la geoestadística de múltiples puntos puede ser aplicada en una aplicación realista para determinar el impacto de una heterogeneidad geológica compleja en el flujo y transporte de aguas subterráneas. La aproximación es aplicada a un acuífero real en Bélgica que exhibe una complejo heterogeneidad y anisotropía sedimentaria. Una imagen de ensayo se construyó basada en datos de hidrogeológicos y geológicos de campo. La estadística de múltiples puntos es tomada de esta imagen de ensayo para simular la ocurrencia de hidrofacies, mientras que la variabilidad de la permeabilidad de intrafacies está simulada usando métodos geoestadísticos convencionales basados en variogramas. Los ensayos simulados de la conductividad hidráulica son usados como entrada al modelo de flujo y transporte para investigar el efecto de la heterogeneidad sedimentaria de pequeña escala en la migración de la pluma contaminante. Los resultados muestran que la heterogeneidad sedimentaria de pequeña escala tiene un efecto significativo sobre el transporte de contaminantes en el acuífero estudiado. La incertidumbre en la distribución espacial de las facies y la distribución de la conductividad hidráulica de intrafacies da como resultado en una significativa incertidumbre sobre el cálculo de la distribución de la concentración. La comparación con las técnicas standard basadas en variogramas muestra que la geoestadística de múltiples puntos permite una mejor reproducción de cortinas de arcilla de forma irregular y baja permeabilidad que influyen en el transporte de soluto.

摘要

沉积过程往往会导致水文地质参数在地下空间的复杂的非均质性。基于方差图的随机方法往往不能描述这种复杂地质环境的非均质性。本次工作展示了如何将多点统计地质学应用于实际水文地质问题, 确定复杂的地质非均质性对地下水流和溶质运移的影响。该方法被用于比利时的一个实际的具有复杂的沉积非均质性和各向异性含水层中。根据地质和水文地质实测资料绘制了训练图像。从这个训练图像中得出多点统计数值用于模拟水相分布, 亚相渗透率变化率利用传统的基于方差图的地质统计方法进行模拟。模拟出的水力传导度作为地下水流和溶质运移模型的输入值, 用来研究小范围内沉积造成的非均质性对污染物晕运移的影响。结果显示, 小范围的沉积造成的非均质性对研究目的含水层中的污染物运移有重要影响。空间相分配和亚相水力传导系数分配的不确定性导致计算出的浓度分布的显著不确定性。同标准的基于方差图方法对比, 多点统计地质学能够更好地再现影响溶质运移的不规则形态的低渗透率粘土隔层。

Resumo

Os processos sedimentológicos resultam frequentemente em heterogeneidades subsuperficiais tridimensionais complexas dos valores dos parâmetros hidrogeológicos. Abordagens estocásticas baseadas na análise variográfica são frequentemente incapazes de descrever a heterogeneidade nesses ambientes geológicos complexos. O presente trabalho evidencia como a geoestatística multiponto pode ser utilizada numa aplicação hidrogeológica realista para determinar o impacte da heterogeneidade geológica complexa no fluxo e transporte hídrico subterrâneo. A abordagem é aplicada a um aquífero real na Bélgica que apresenta uma heterogeneidade e anisotropia sedimentar complexa. É construída uma imagem de treino baseada em dados de campo de natureza geológica e hidrogeológica. As estatísticas multiponto são retiradas desta imagem de treino para simular a ocorrência de hidrofácies, enquanto a variabilidade da permeabilidade intrafácies é simulada utilizando métodos geoestatísticos baseados na análise variográfica convencional. Os resultados da simulação da condutividade hidráulica são utilizados como “input” no modelo de simulação de fluxo e de transporte para investigar o efeito de heterogeneidades sedimentares de pequena escala na migração de plumas de contaminação. Os resultados mostram que as heterogeneidades sedimentares de pequena escala têm um efeito significativo sobre o transporte de contaminantes no aquífero estudado. A incerteza na distribuição espacial das fácies e na distribuição da condutividade hidráulica intrafácies resulta numa incerteza significativa na distribuição das concentrações calculadas. A comparação com técnicas clássicas de análise variográfica mostra que a geoestatística multiponto permite uma melhor reprodução das barreiras argilosas irregulares de baixa permeabilidade que influenciam o transporte em solução.