Doctoral Thesis
Refine
Year of publication
- 2023 (1)
Document Type
- Doctoral Thesis (1) (remove)
Language
- English (1)
Has Fulltext
- yes (1)
Is part of the Bibliography
- no (1) (remove)
Keywords
- Hydrogeologie (1) (remove)
Institute
Underground hard coal mining operations irreversibly disrupt the pre-existing mechanical equilibrium of the geological media. The employment of high-recovery methods modifies the stress field of the sedimentary sequence, generating movement and faulting of the rock layers above and below mined seams. These new fracture zones do affect the original conditions of the hydrogeological system by modifying flow pathways and increasing the permeability of the rock sequence. Moreover, the surface area of rock exposed to air and water is increased, conditioning the water-rock interaction. Despite this rather clear conceptualization, flow and reactive transport processes in fractured overburdens are rarely modeled simultaneously. Discrete setups that consider fractures and porous matrix require extensive characterization of both media, which is impractical for regional case studies. As a result, most post-mining models explicitly ignore fracture structures by employing the equivalent porous approach or even both media with lumped parameter models. However, omitting either medium represents a delicate simplification, considering that mining-related fractures control the rate and direction of water flow within moderately permeable but relatively highly porous rock sequences.
In this dissertation, the specific contribution of fractured and matrix continua to the transient discharge and water quality of a post-mining coal zone is quantified and evaluated. For this purpose, dual and multiple interacting continua models are employed to simulate fluid flow and reactive mass transport in fractured and variable water-saturated rock sequences. The effectiveness of the models is evaluated by simulating the origin, generation and transport of acid mine drainage (i.e., water with elevated concentrations of hydrogen, iron, sulfate and chloride) within the shallow overburden of the Ibbenbüren Westfield. Compared to other coal districts in Germany, this area is strongly delimited by the local geology and topography, resulting in a well-defined hydrogeological system to test the models. Petrographic and chemical analyses performed on core samples from the area show the strong influence of mining-derived fractures on the water-rock interaction within the Carboniferous sequence. The presence of oxidized pyrite along with amorphous iron hydroxide phases in weathering fronts on both sides of the fractures demonstrates the exchange of solutes and gases between the fractured and the porous matrix media.
Based on the previous evidence, the TOUGHREACT software is employed to characterize flow and reactive transport processes in the Westfield. However, each of the two processes is simulated at separate stages to have more control in the adjustment of sensitive parameters for which little information is available. For the flow component, a dual continuum model, with Richard’s equations is used to characterize the unsaturated water flow in both fractured and matrix media. Under this approach, the model adequately reproduces the bimodal flow behavior of the discharges measured in the mine drainage for the years 2008 and 2017. Simulation results show how the fractured continuum generates intense discharge events during the winter months while the rock matrix controls smooth discharge limbs in summer, when water is slowly released back to the fractures. With the flow component calibrated, the second part of the study incorporates the geochemical processes into the model based on actual data from the rock samples. Their simulation requires extending the two-continuum setup to a multiple continua model with five nested block strings: one for the fractures and four for the rock matrix. This further subdivision prevents under-representations of kinetic reactions with short equilibrium length scales and numerical instabilities due to lack of chemical and flow gradients. As a result, the new multiple continua model provides good agreement with respect to long- and short-term concentrations and discharge trends measured in the mine drainage. The flow of oxygen and meteoric water through the fractured continuum leads to a high and steady release of hydrogen, iron and sulfate ions derived from pyrite oxidation in the matrix continua closest to the fractures. Moreover, high chloride concentrations result from the mixing and gradual release of relatively immobile solutes in the matrix as they interact with percolating water in the fracture. Both findings are equally congruent with the reactive pyrite oxidation and iron hydroxide precipitation fronts identified in the fractured core samples.
In the end, the multiple continua models, the simulation procedure and the results of the benchmark and sensitivity analysis scenarios developed for the Westfield pave the way for the application of the approach in other mining zones. The first candidate emerges in the Ibbenbüren Eastfield, where a coupled elemental-isotopic approach included in this thesis has confirmed that water-conducting fracture zones are primary elements for solute generation and transport in the first 300 meters of the overburden. In the latter case, calibration and verification of the models can be complemented with measurements of δ34S in sulfates and δ18O, δ2H, and Tritium in water.