Sensitivity Study of the Reactive Transport Model for CO2 Injection into the Utsira Saline Formation Using 3D Fluid Flow Model History Matched with 4D Seismic

AbstractIn the UK, disposal of packaged intermediate-level radioactive waste may involve waste packages being placed in a geological disposal facility (GDF) and surrounded by a cementitious backfill. Cracking of the backfill could occur due to a number of mechanisms, and this could affect the post-closure performance of the GDF.This work has assessed potential cracking in the backfill during the backfilling and early post-closure period of GDF vaults with an open crown space in a higher strength rock. From the comprehensive range of processes considered, three were identified as potentially causing cracking: (1) during backfilling, plastic settlement under solid horizontal surfaces could result in horizontal gaps beneath waste packages; (2) within days of backfilling, early-age thermal contraction of the backfill could result in primarily vertical cracks; (3) over a number of years, expansion of waste packages could result in large horizontal cracks.A groundwater flow model incorporating a representation of the cracks was used to calculate flows through a backfilled GDF vault, and through the cracks themselves. Including cracks increased the flow rate significantly. A reactive transport model was used to estimate the evolution of the pore water chemistry as groundwater flows through the cracked backfill. Calcite and brucite were predicted to precipitate, with brucite subsequently dissolving. Calcite build-up could seal some cracks.

Download Full-text

Inhomogeneous rock compositions and varying dissolution rates affect evolution and shape of leaching zones in potash seams

10.5194/egusphere-egu21-1196 ◽

2021 ◽

Author(s):

Svenja Steding ◽

Thomas Kempka ◽

Axel Zirkler ◽

Michael Kühn

Keyword(s):

Fluid Flow ◽

Reactive Transport ◽

Transport Model ◽

Chemical Species ◽

Simulation Environment ◽

Dissolution Rates ◽

Reactive Transport Model ◽

Hazard Potential ◽

Transport Simulation ◽

Salt Deposits

Salt deposits host an important industrial raw material and provide storage capacities for energy and nuclear waste. However, leaching zones can seriously endanger the development and utilisation of salt deposits for these purposes, especially if these occur in potash seams. Their increased solubility enables even NaCl-saturated solutions, if present, to deeply penetrate these seams. The resulting salt dissolution processes generate fluid flow paths and affect the mechanical rock integrity. To model the timely evolution of leaching zones and to assess their hazard potential, a reactive transport model has been developed, taking into account not only the complex dissolution and precipitation behaviour of potash salts, but also the resulting porosity and permeability changes as well as density-driven chemical species transport. Additionally, the model makes use of an approach to describe transport and chemical reactions at the interface between impermeable (dry) salt rocks and permeated leaching zones (Steding et al., 2021). In the present study, we focus on the effect of heterogeneity of the mineral distribution within potash seams and on the influence of mineral- and saturation-dependent dissolution rates.The applied reactive transport model is based on a coupling of the geochemical module PHREEQC (Parkhurst & Appelo, 2013) with the TRANSport Simulation Environment (Kempka, 2020) as well as the newly developed extension of an interchange approach (Steding et al., 2021). A numerical model has been developed and applied to simulate the leaching process of a carnallite-bearing potash seam due to natural density-driven convection. The results show that both, the mineral composition and dissolution rate of the original salt rock, strongly influence the shape and evolution of the leaching zone (Steding et al., 2021).In nature, strong variations of the mineralogy occur within potash seams with random or stratified distributions. Furthermore, dissolution rates depend on the mineral itself as well as on its saturation state. Both may considerably influence the growth rate of a leaching zone. Therefore, the reactive transport model has been extended by mineral- and saturation-dependent dissolution rates. A scenario analysis has been undertaken to compare the impact of homogeneous and heterogeneous rock compositions. For that purpose, the carnallite content in the potash seam was varied from 5 to 25&#160;wt.&#160;% including different stratifications and random distributions. The simulations were classified by means of the P&#233;clet and Damk&#246;hler numbers, and the long-term behaviour as well as hazard potential are discussed.&#160;References:Parkhurst, D.L.; Appelo, C.A.J. (2013). Description of Input and Examples for PHREEQC Version 3 - a Computer Program for Speciation, Batch-reaction, One-dimensional Transport, and Inverse Geochemical Calculations. In Techniques and Methods; Publisher: U.S. Geological Survey; Book 6, 497 ppKempka, T. (2020). Verification of a Python-based TRANsport Simulation Environment for density-driven fluid flow and coupled transport of heat and chemical species. Adv. Geosci. 54, 67&#8211;77. Steding, S.; Kempka, T.; Zirkler, A.; K&#252;hn, M. (2021). Spatial and temporal evolution of leaching zones within potash seams reproduced by reactive transport simulations. Water 13, 168.

Download Full-text