3D Modeling of Saline Groundwater Flow and Transport in a Flooded Salt Mine in Stassfurt, Germany

Abstract

The abandoned, flooded deep salt mine in Stassfurt, Germany has been extensively studied. The mine causes serious uncontrolled salt leaching and mass displacement effects underground that have led to significant regional-scale subsidence. Collapsing mined areas and further large-scale subsidence pose a potential danger to the city of Stassfurt. Comprehensive modeling of the three dimensional (3-D) variable density flow dynamics is required at both large and local scale to assess the stability of the cap rocks and the potential of subsidence due to salt leaching processes, and to better control (minimize) salt leaching. The impact of different engineering solutions and appropriate remedial strategies (e.g. pumping schemes) on the subsurface spatio-temporal distribution of salinity is of specific concern. Coupled variable density flow and transport simulations were performed using FEFLOW, a finite-element simulator, which is capable of both porous media (Darcy-type) flow computations and lower-dimensional discrete elements for laminar or turbulent flow. The geometry of the mine workings is incorporated in suitable detail and resolution in the developed 3-D model. Computational results were successfully compared with field observations and geophysical measurements. Different model scenarios provide better insight into flow and mass transport processes in the study area, which could provide significant information to quantify subsidence potential of the mine site.

Zusammenfassung

Das stillgelegte und geflutete Untertagesalzbergwerk Staßfurt wurde ausführlich untersucht. Beträchtliche unkontrollierte Salzauslaugung und untertägige Gesteinsdeformationen haben in der Umgebung des Bergwerks zu signifikanten Bodenabsenkungen geführt. Für die Stadt Staßfurt bedeuten diese Verbrüche und die großflächigen Absenkungen eine potentielle Gefahr. Um die Salzauslaugung besser kontrollieren und minimieren zu können sowie um die Stabilität des Deckgebirges und die Gefahr von Absenkungen durch Auslaugung abzuschätzen, sind umfangreiche klein- und großräumige 3D-Modellierungen der sich ändernden Dichteströmungsdynamik erforderlich. Die Auswirkungen unterschiedlicher technischer Maßnahmen und geeigneter Sanierungsvarianten (z.B. Pumpenprojekte) auf die räumliche und zeitliche Verteilung der Salinität untertage sind daher von besonderer Bedeutung. Gekoppelte Transport- und veränderliche Dichteströmungsmodelle wurden mit dem Finite Elemente Simulator FEFLOW durchgeführt. Dieses Programm kann sowohl Strömungen in porösen Medien (Darcy-Strömung) als auch laminare oder turbulente Strömungen in diskreten Elementen berechnen. Das 3D-Model beinhaltet die Geometrie der Grubenhohlräume in geeigneter Detailtreue und Auflösung. Die Berechnungsergebnisse zeigen eine gute Übereinstimmung mit den Geländebeobachtungen und geophysikalischen Messwerten. Unterschiedliche Modellierungsszenarien ermöglichen im Untersuchungsgebiet einen besseren Einblick in die Prozesse des Strömungs- und Massentransports. Sie liefern damit wichtige Informationen, um mögliche Absenkungen des Grubengeländes zu quantifizieren.

Resumen

La mina de sal inundada y abandonada en Stassfurt, Alemania, ha sido extensamente estudiada. La mina causa una lixiviación seria e incontrolada de sal y efectos de desplazamiento subterráneo de masa que provoca un significativo hundimiento. Las áreas colapsadas de la mina y un posterior hundimiento a gran escala, plantean un potencial riesgo para la ciudad de Stassfurt. Un modelado comprensivo en tres dimensiones (3-D) de la dinámica de flujo, es requerido a gran escala y a escala local para relevar la estabilidad de las rocas y el potencial de hundimiento debido a los procesos de lixiviación de sal, y para un mejor control (minimización) de la lixiviación de sal. El impacto de diferentes soluciones ingenieriles y apropiadas estrategias de remediación (por ejemplo, esquemas de bombeo) sobre la distribución espacio-temporal de la salinidad debajo de la superficie, es la preocupación específica. Simulaciones que acoplen el flujo de densidad y el transporte, fueron realizadas usando un simular finito, FEFLOW, que es capaz de computar tanto flujos en medios porosos (Darcy-tipo) y elementos discretos para flujos laminares o turbulentos. La geometría de los trabajos en la mina es incorporada en adecuado detalle y resolución en el modelo 3-D. Los resultados computacionales fueron exitosamente comparados con observaciones de campo y mediciones geofísicas. Diferentes escenarios para el modelo proporcionan una mejor compresión de los procesos de transporte de masa y de flujo en el área de estudio, que podría generar significativa información para cuantificar el potencial de hundimiento del sitio de la mina.

抽象

德国斯塔斯弗德 (Stassfurt) 地区的废弃充水盐矿已被广泛研究. 该盐矿难以控制的盐溶滤作用和置换作用已经引发显著的区域地面沉降. 矿区塌陷及大范围地面沉降对斯塔斯弗德 (Stassfurt) 城市安全构成潜在危胁. 为了评价区域和局部基岩稳定性以及盐溶滤引起的地面沉降、有效控制盐溶滤作用 (使其最小化), 需要建立矿区地下水的三维可变密度水动力学模型. 不同的工程措施和修复方案 (如抽水计划等) 对地下水盐分时空分布的影响是明显的. 有限元模拟软件FEFLOW不仅能够完成多孔介质达西流计算, 而且能够进行层流和紊流低维离散元计算. 本文利用FEFLOW耦合模拟了变密度流和物质运移过程. 在3-D模型中, 以合理的详细程度和分辨率考虑了工作面几何形状. 模拟计算结果与现场观测和地球物理探测结果吻合. 不同情景下模型的建立及讨论有助于更好地理解研究区水流和溶质运移过程, 能够为定量评价矿区地面降深提供了重要依据.