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

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

A Comparative Review of Hydrologic Issues Involved in Geologic Storage of CO2 and Injection Disposal of Liquid Waste

The paper presents a comparison of hydrologic issues and technical approaches used in deep-well injection and disposal of liquid wastes, and those issues and approaches associated with injection and storage of CO{sub 2} in deep brine formations. These comparisons have been discussed in nine areas: (1) Injection well integrity; (2) Abandoned well problems; (3) Buoyancy effects; (4) Multiphase flow effects; (5) Heterogeneity and flow channeling; (6) Multilayer isolation effects; (7) Caprock effectiveness and hydrogeomechanics; (8) Site characterization and monitoring; and (9) Effects of CO{sub 2} storage on groundwater resources There are considerable similarities, as well as significant differences. Scientifically and technically, these two fields can learn much from each other. The discussions presented in this paper should help to focus on the key scientific issues facing deep injection of fluids. A substantial but by no means exhaustive reference list has been provided for further studies into the subject

Coupled In-Rock and In-Drift Hydrothermal Model Stuudy For Yucca Mountain

A thermal-hydrologic-natural-ventilation model is configured for simulating temperature, humidity, and condensate distributions in the coupled domains of the in-drift airspace and the near-field rockmass in the proposed Yucca Mountain repository. The multi-physics problem is solved with MULTIFLUX in which a lumped-parameter computational fluid dynamics model is iterated with TOUGH2. The solution includes natural convection, conduction, and radiation for heat as well as moisture convection and diffusion for moisture transport with half waste package scale details in the drift, and mountain-scale heat and moisture transport in the porous and fractured rock-mass. The method provides fast convergence on a personal computer computational platform. Numerical examples and comparison with a TOUGH2 based, integrated model are presented

Modeling Single Well Injection-Withdrawal (SWIW) Tests for Characterization of Complex Fracture-Matrix Systems

The ability to reliably predict flow and transport in fractured porous rock is an essential condition for performance evaluation of geologic (underground) nuclear waste repositories. In this report, a suite of programs (TRIPOLY code) for calculating and analyzing flow and transport in two-dimensional fracture-matrix systems is used to model single-well injection-withdrawal (SWIW) tracer tests. The SWIW test, a tracer test using one well, is proposed as a useful means of collecting data for site characterization, as well as estimating parameters relevant to tracer diffusion and sorption. After some specific code adaptations, we numerically generated a complex fracture-matrix system for computation of steady-state flow and tracer advection and dispersion in the fracture network, along with solute exchange processes between the fractures and the porous matrix. We then conducted simulations for a hypothetical but workable SWIW test design and completed parameter sensitivity studies on three physical parameters of the rock matrix - namely porosity, diffusion coefficient, and retardation coefficient - in order to investigate their impact on the fracture-matrix solute exchange process. Hydraulic fracturing, or hydrofracking, is also modeled in this study, in two different ways: (1) by increasing the hydraulic aperture for flow in existing fractures and (2) by adding a new set of fractures to the field. The results of all these different tests are analyzed by studying the population of matrix blocks, the tracer spatial distribution, and the breakthrough curves (BTCs) obtained, while performing mass-balance checks and being careful to avoid some numerical mistakes that could occur. This study clearly demonstrates the importance of matrix effects in the solute transport process, with the sensitivity studies illustrating the increased importance of the matrix in providing a retardation mechanism for radionuclides as matrix porosity, diffusion coefficient, or retardation coefficient increase. Interestingly, model results before and after hydrofracking are insensitive to adding more fractures, while slightly more sensitive to aperture increase, making SWIW tests a possible means of discriminating between these two potential hydrofracking effects. Finally, we investigate the possibility of inferring relevant information regarding the fracture-matrix system physical parameters from the BTCs obtained during SWIW testing

Forecast of thermal-hydrological conditions and air injection test results of the single heater test at Yucca Mountain

The heater in the Single Heater Test (SHT) in alcove 5 of the Exploratory Studies Facility (ESF) was turned on August 26, 1996. A large number of sensors are installed in the various instrumented boreholes to monitor the coupled thermal-hydrological-mechanical-chemical responses of the rock mass to the heat generated in the single heater. In this report the authors present the results of the modeling of both the heating and cooling phases of the Single Heater Test (SHT), with focus on the thermal-hydrological aspect of the coupled processes. Also in this report, the authors present simulations of air injection tests will be performed at different stages of the heating and cooling phase of the SHT

Sensitivity study of CO2 storage capacity in brine aquifers withclosed boundaries: Dependence on hydrogeologic properties

In large-scale geologic storage projects, the injected volumes of CO{sub 2} will displace huge volumes of native brine. If the designated storage formation is a closed system, e.g., a geologic unit that is compartmentalized by (almost) impermeable sealing units and/or sealing faults, the native brine cannot (easily) escape from the target reservoir. Thus the amount of supercritical CO{sub 2} that can be stored in such a system depends ultimately on how much pore space can be made available for the added fluid owing to the compressibility of the pore structure and the fluids. To evaluate storage capacity in such closed systems, we have conducted a modeling study simulating CO{sub 2} injection into idealized deep saline aquifers that have no (or limited) interaction with overlying, underlying, and/or adjacent units. Our focus is to evaluate the storage capacity of closed systems as a function of various reservoir parameters, hydraulic properties, compressibilities, depth, boundaries, etc. Accounting for multi-phase flow effects including dissolution of CO{sub 2} in numerical simulations, the goal is to develop simple analytical expressions that provide estimates for storage capacity and pressure buildup in such closed systems

Estimating maximum sustainable injection pressure duringgeological sequestration of CO2 using coupled fluid flow andgeomechanical fault-slip analysis

This paper demonstrates the use of coupled fluid flow andgeomechanical fault slip (fault reactivation) analysis to estimate themaximum sustainable injection pressure during geological sequestration ofCO2. Two numerical modeling approaches for analyzing faultslip areapplied, one using continuum stress-strain analysis and the other usingdiscrete fault analysis. The results of these two approaches to numericalfault-slip analyses are compared to the results of a more conventionalanalytical fault-slip analysis that assumes simplified reservoirgeometry. It is shown that the simplified analytical fault-slip analysismay lead to either overestimation or underestimation of the maximumsustainable injection pressure because it cannot resolve importantgeometrical factors associated with the injection induced spatialevolution of fluid pressure and stress. We conclude that a fully couplednumerical analysis can more accurately account for the spatial evolutionof both insitu stresses and fluid pressure, and therefore results in amore accurate estimation of the maximum sustainable CO2 injectionpressure