Lattice modeling of excavation damage in argillaceous clay formations: Influence of deformation and strength anisotropy

This paper presents modeling of mechanical anisotropy in argillaceous rocks using an irregular lattice modeling approach, namely the rigid-body-spring network. To represent the mechanical anisotropy, new schemes are implemented in the modeling framework. The directionality of elastic deformation is resolved by modifying the element formulation with anisotropic elastic properties. The anisotropy of strength and failure characteristics is facilitated by adopting orientation-dependent failure criteria into the failure model. The verification of the improved modeling procedures is performed against theoretical model predictions for unconfined compression tests with various bedding orientations. Furthermore, excavation damage and fracturing processes in rock formations are simulated for different geomechanical configurations, such as rock anisotropy and tectonic heterogeneity. The simulated excavation damage characteristics are realistic and comparable with the actual field observation at a tunnel located in an argillaceous clay formation. The simulation results provide insights into the excavation damage zone phenomena with an explicit representation of fracturing processes

Disposal R&D in the Used Fuel Disposition Campaign: A Discussion of Opportunities for Active International Collaboration

For DOE's Used Fuel Disposition Campaign (UFDC), international collaboration is a beneficial and cost-effective strategy for advancing disposal science with regards to multiple disposal options and different geologic environments. While the United States disposal program focused solely on Yucca Mountain tuff as host rock over the past decades, several international programs have made significant progress in the characterization and performance evaluation of other geologic repository options, most of which are very different from the Yucca Mountain site in design and host rock characteristics. Because Yucca Mountain was so unique (e.g., no backfill, unsaturated densely fractured tuff), areas of direct collaboration with international disposal programs were quite limited during that time. The decision by the U.S. Department of Energy to no longer pursue the disposal of high-level radioactive waste and spent fuel at Yucca Mountain has shifted UFDC's interest to disposal options and geologic environments similar to those being investigated by disposal programs in other nations. Much can be gained by close collaboration with these programs, including access to valuable experience and data collected over recent decades. Such collaboration can help to efficiently achieve UFDC's long-term goals of conducting 'experiments to fill data needs and confirm advanced modeling approaches' (by 2015) and of having a 'robust modeling and experimental basis for evaluation of multiple disposal system options' (by 2020). This report discusses selected opportunities of active international collaboration, with focus on both Natural Barrier System (NBS) and Engineered Barrier System (EBS) aspects and those opportunities that provide access to field data (and respective interpretation/modeling) or allow participation in ongoing field experiments. This discussion serves as a basis for the DOE/NE-53 and UFDC planning process for FY12 and beyond

Using pressure and volumetric approaches to estimate CO2 storage capacity in deep saline aquifers

Various approaches are used to evaluate the capacity of saline aquifers to store CO2, resulting in a wide range of capacity estimates for a given aquifer. The two approaches most used are the volumetric “open aquifer” and “closed aquifer” approaches. We present four full-scale aquifer cases, where CO2 storage capacity is evaluated both volumetrically (with “open” and/or “closed” approaches) and through flow modeling. These examples show that the “open aquifer” CO2 storage capacity estimation can strongly exceed the cumulative CO2 injection from the flow model, whereas the “closed aquifer” estimates are a closer approximation to the flow-model derived capacity. An analogy to oil recovery mechanisms is presented, where the primary oil recovery mechanism is compared to CO2 aquifer storage without producing formation water; and the secondary oil recovery mechanism (water flooding) is compared to CO2 aquifer storage performed simultaneously with extraction of water for pressure maintenance. This analogy supports the finding that the “closed aquifer” approach produces a better estimate of CO2 storage without water extraction, and highlights the need for any CO2 storage estimate to specify whether it is intended to represent CO2 storage capacity with or without water extraction

Reduced Order Models for Prediction of Groundwater Quality Impacts from CO2 and Brine Leakage

AbstractA careful assessment of the risk associated with geologic CO2 storage is critical to the deployment of large-scale storage projects. A potential risk is the deterioration of groundwater quality caused by the leakage of CO2 and brine leakage from deep subsurface reservoirs. In probabilistic risk assessment studies, numerical modeling is the primary tool employed to assess risk. However, the application of traditional numerical models to fully evaluate the impact of CO2 leakage on groundwater can be computationally complex, demanding large processing times and resources, and involving large uncertainties. As an alternative, reduced order models (ROMs) can be used as highly efficient surrogates for the complex process-based numerical models.In this study, we represent the complex hydrogeological and geochemical conditions in a heterogeneous aquifer and subsequent risk by developing and using two separate ROMs. The first ROM is derived from a model that accounts for the heterogeneous flow and transport conditions in the presence of complex leakage functions for CO2 and brine. The second ROM is obtained from models that feature similar, but simplified flow and transport conditions, and allow for a more complex representation of all relevant geochemical reactions. To quantify possible impacts to groundwater aquifers, the basic risk metric is taken as the aquifer volume in which the water quality of the aquifer may be affected by an underlying CO2 storage project. The integration of the two ROMs provides an estimate of the impacted aquifer volume taking into account uncertainties in flow, transport and chemical conditions. These two ROMs can be linked in a comprehensive system level model for quantitative risk assessment of the deep storage reservoir, wellbore leakage, and shallow aquifer impacts to assess the collective risk of CO2 storage projects

A new Eulerian-Lagrangian finite element simulator for solute transport in discrete fracture-matrix systems

Fracture network simulators have extensively been used in the past for obtaining a better understanding of flow and transport processes in fractured rock. However, most of these models do not account for fluid or solute exchange between the fractures and the porous matrix, although diffusion into the matrix pores can have a major impact on the spreading of contaminants. In the present paper a new finite element code TRIPOLY is introduced which combines a powerful fracture network simulator with an efficient method to account for the diffusive interaction between the fractures and the adjacent matrix blocks. The fracture network simulator used in TRIPOLY features a mixed Lagrangian-Eulerian solution scheme for the transport in fractures, combined with an adaptive gridding technique to account for sharp concentration fronts. The fracture-matrix interaction is calculated with an efficient method which has been successfully used in the past for dual-porosity models. Discrete fractures and matrix blocks are treated as two different systems, and the interaction is modeled by introducing sink/source terms in both systems. It is assumed that diffusive transport in the matrix can be approximated as a one-dimensional process, perpendicular to the adjacent fracture surfaces. A direct solution scheme is employed to solve the coupled fracture and matrix equations. The newly developed combination of the fracture network simulator and the fracture-matrix interaction module allows for detailed studies of spreading processes in fractured porous rock. The authors present a sample application which demonstrate the codes ability of handling large-scale fracture-matrix systems comprising individual fractures and matrix blocks of arbitrary size and shape

A new Lagrangian-Eulerian finite element method for modeling contaminant transport in fractured porous formations

Fracture network simulators have been extensively used in the past for obtaining a better understanding of flow and transport processes in fractured rock. However, most of these models do not account for fluid or solute exchange between the fractures and the porous matrix, although diffusion into the matrix pores can have a major impact on the spreading of contaminants. In the present paper a new finite element code TRIPOLY is introduced which combines a powerful Lagrangian-Eulerian approach for solving flow and transport in networks of discrete fractures with an efficient method to account for the diffusive interaction between the fractures and the adjacent matrix blocks. The code is capable of handling large-scale fracture-matrix systems comprising individual fractures and matrix blocks of arbitrary size, shape, and dimension

A Turbulent Transport Network Model in MULTIFLUX Coupled with TOUGH2

A new numerical method is described for the fully iterated, conjugate solution of two discrete submodels, involving (a) a transport network model for heat, moisture, and airflows in a high-permeability, air-filled cavity; and (b) a variably saturated fractured porous medium. The transport network submodel is an integrated-parameter, computational fluid dynamics solver, describing the thermal-hydrologic transport processes in the flow channel system of the cavity with laminar or turbulent flow and convective heat and mass transport, using MULTIFLUX. The porous medium submodel, using TOUGH2, is a solver for the heat and mass transport in the fractured rock mass. The new model solution extends the application fields of TOUGH2 by integrating it with turbulent flow and transport in a discrete flow network system. We present demonstrational results for a nuclear waste repository application at Yucca Mountain with the most realistic model assumptions and input parameters including the geometrical layout of the nuclear spent fuel and waste with variable heat load for the individual containers. The MULTIFLUX and TOUGH2 model elements are fully iterated, applying a programmed reprocessing of the Numerical Transport Code Functionalization model-element in an automated Outside Balance Iteration loop. The natural, convective airflow field and the heat and mass transport in a representative emplacement drift during postclosure are explicitly solved in the new model. The results demonstrate that the direction and magnitude of the air circulation patterns and all transport modes are strongly affected by the heat and moisture transport processes in the surrounding rock, justifying the need for a coupled, fully iterated model solution such as the one presented in the paper

On Water Flow in Hot Fractured Rock -- A Sensitivity Study on theImpact of Fracture-Matrix Heat Transfer

Dual-continuum models have been widely used in modeling flowand transport in fractured porous rocks. Among many other applications,dual-continuum approaches were utilized in predictive models of thethermal-hydrological conditions near emplacement tunnels (drifts) atYucca Mountain, Nevada, the proposed site for a radioactive wasterepository in the U.S. In unsaturated formations such as those at YuccaMountain, the magnitude of mass and heat exchange between the twocontinua fracture network and matrix is largely dependent on the flowcharacteristics in the fractures, because channelized finger-type flowstrongly reduces the interface area between the matrix surfaces and theflowing liquid. This effect may have important implications, for example,during the time period that the fractured rock near the repository driftswould be heated above the boiling point of water. Depending on themagnitude of heat transfer from the matrix, water percolating down thefractures will either boil off in the hot rock region above drifts or maypenetrate all the way to the drift walls and possibly seep into the opencavities. In this paper, we describe a sensitivity analysis using avariety of approaches to treat fracture-matrix interaction in athree-dimensional dual-continuum setting. Our simulation example is alaboratory heater experiment described in the literature that providesevidence of rapid water flow in fractures, leading to drift seepagedespite above-boiling conditions in the adjacent fractured rock. Theexperimental finding can only be reproduced when the interface area forheat transfer between the matrix and fracture continua is reduced toaccount for flow channeling

Pressure control for managing and optimizing adjacent subsurface operations in large scale CCS

Injecting CO2 in to the subsurface for safe storage of CO2 the pressure propagates far away from the injection point and this can be a potential problem if the overpressure extents to neighbouring subsurface activities or potential leakage pathways. For structural closure trap configurations the CO2 plume is captured within the local structural closure but the pressure footprint is on a more regional scale. This rise the question on, how large the storage complex needs to be for any individual storage operations and how large an area monitoring activities have to cover. The EC CCS guidance document addresses the issues with statements on competitions between subsurface operations but returns no absolute values. Pressure modelling of CO2 injection process with state of the art reservoir simulation tools is challenges by use of realistic model boundary conditions in order to model a realistic pressure level. Combined use of models on a site scale and on a regional scale can instruct how boundary conditions are set-up for a site scale model. Pressure management through pressure release wells could be an option to mitigate undesirable over-pressure developments. For local structural closures the pressure release wells can be placed outside the closure hereby mitigate the overpressure without introducing a potential leakage by drilling inside the trap. The paper addresses the issue of selecting model boundary conditions and modelling mitigation of pressure development by use of a large regional model with local structural traps in the Bunter Sandstone Formation in the UK Southern North Sea

Natural convection in tunnels at Yucca Mountain and impact on drift seepage

The decay heat from radioactive waste that is to be disposed in the once proposed geologic repository at Yucca Mountain (YM) will significantly influence the moisture conditions in the fractured rock near emplacement tunnels (drifts). Additionally, large-scale convective cells will form in the open-air drifts and will serve as an important mechanism for the transport of vaporized pore water from the fractured rock in the drift center to the drift end. Such convective processes would also impact drift seepage, as evaporation could reduce the build up of liquid water at the tunnel wall. Characterizing and understanding these liquid water and vapor transport processes is critical for evaluating the performance of the repository, in terms of water-induced canister corrosion and subsequent radionuclide containment. To study such processes, we previously developed and applied an enhanced version of TOUGH2 that solves for natural convection in the drift. We then used the results from this previous study as a time-dependent boundary condition in a high-resolution seepage model, allowing for a computationally efficient means for simulating these processes. The results from the seepage model show that cases with strong natural convection effects are expected to improve the performance of the repository, since smaller relative humidity values, with reduced local seepage, form a more desirable waste package environment