Transport Modeling in a Finite Fractured Rock Domain

Abstract

ABSTRACTThe Waste Isolation Pilot Plant (WIPP) is a U.S. Department of Energy facility intended to demonstrate the safe disposal of transuranic nuclear waste. The WIPP is located in the thick halite beds of the Salado Formation in the Delaware Basin of southeastern New Mexico. Overlying the repository is the Culebra Dolomite Member of the Rustler Formation, a dolomitic unit containing accessory minerals such as calcite, gypsum, and clays. The Culebra is the predominant continuous water-bearing unit in the area. Part of the scientific demonstration of the safety of WIPP involves understanding and estimating what might happen should the WIPP inadvertently be intruded by exploration boreholes long after the WIPP has been decommissioned. Should such a breach occur, it is possible that dissolved transuranic waste could be transported up a borehole and into the Culebra, through which it could eventually reach the accessible environment. The rates and mechanisms for actinide transport in Culebra Dolomite are being investigated as a component of repository performance assessment.The two primary mechanisms for delaying solute release from a double porosity system are called physical and chemical retardation, both of which retard solute migration relative to bulk water flow. Physical retardation occurs when solutes diffuse out of the advective transport regions (fractures) and into diffusive transport regions (matrix). Results of hydrologic field tests suggest that the Culebra behaves as a double porosity medium. Chemical retardation occurs when chemical interactions between dissolved species and mineral surfaces remove solutes from the aqueous phase. Some fracture surfaces in the Culebra are clay-lined, and, because clays can have a large potential for chemical interactions with solutes (e.g., ion exchange or adsorption), it is possible that fracture linings alone may cause significant retardation of any actinides that could reach the Culebra. This report examines these two features in Culebra transport: delimiting the parameter space where a fully two-dimensional model is needed, and examining the importance of clay linings during transport.Numerical simulations with the FMT (Fracture Matrix Transport) finite difference model in a finite, idealized double porosity system were performed to: 1) delineate parameter values for which a set of two coupled one-dimensional equations is insufficient to describe behavior and a fully two-dimensional representation is needed, 2) compare chemical retardation by ion exchange for fractures with and without a clay lining, and 3) compare physical and chemical retardation by ion exchange as a function of fracture velocity. The simulations differ from others in that the matrix is finite, and thus can become saturated with tracers introduced into the domain. The modeling suggests that the chemical retardation caused by clay linings may not be important in a double porosity system. Over a wide range, a single parameter group essentially determines the complexity of the model.