Site-Scale Saturated Zone Flow Model

The purpose of this model report is to document the components of the site-scale saturated-zone flow model at Yucca Mountain, Nevada, in accordance with administrative procedure (AP)-SIII.lOQ, ''Models''. This report provides validation and confidence in the flow model that was developed for site recommendation (SR) and will be used to provide flow fields in support of the Total Systems Performance Assessment (TSPA) for the License Application. The output from this report provides the flow model used in the ''Site-Scale Saturated Zone Transport'', MDL-NBS-HS-000010 Rev 01 (BSC 2003 [162419]). The Site-Scale Saturated Zone Transport model then provides output to the SZ Transport Abstraction Model (BSC 2003 [164870]). In particular, the output from the SZ site-scale flow model is used to simulate the groundwater flow pathways and radionuclide transport to the accessible environment for use in the TSPA calculations. Since the development and calibration of the saturated-zone flow model, more data have been gathered for use in model validation and confidence building, including new water-level data from Nye County wells, single- and multiple-well hydraulic testing data, and new hydrochemistry data. In addition, a new hydrogeologic framework model (HFM), which incorporates Nye County wells lithology, also provides geologic data for corroboration and confidence in the flow model. The intended use of this work is to provide a flow model that generates flow fields to simulate radionuclide transport in saturated porous rock and alluvium under natural or forced gradient flow conditions. The flow model simulations are completed using the three-dimensional (3-D), finite-element, flow, heat, and transport computer code, FEHM Version (V) 2.20 (software tracking number (STN): 10086-2.20-00; LANL 2003 [161725]). Concurrently, process-level transport model and methodology for calculating radionuclide transport in the saturated zone at Yucca Mountain using FEHM V 2.20 are being carried out in the model report, ''Site-Scale Saturated Zone Transport'', MDL-NBS-HS-000010 Rev 01 (BSC 2003 [162419]). The velocity fields are calculated by the flow model, described herein, independent of the transport processes, and are then used as inputs to the transport model. Justification for this abstraction is presented in the model report, ''Saturated Zone Flow and Transport Model Abstraction'', MDL-NBS-HS-000021 (BSC 2003 [164870])

Models and methods summary for the FEHMN application

This models and methods summary provides a detailed description of the mathematical models and numerical methods employed by the finite element heat and mass tranfer code (FEHMN) application

Mathematical Formulation Requirements and Specifications for the Process Models

The Advanced Simulation Capability for Environmental Management (ASCEM) is intended to be a state-of-the-art scientific tool and approach for understanding and predicting contaminant fate and transport in natural and engineered systems. The ASCEM program is aimed at addressing critical EM program needs to better understand and quantify flow and contaminant transport behavior in complex geological systems. It will also address the long-term performance of engineered components including cementitious materials in nuclear waste disposal facilities, in order to reduce uncertainties and risks associated with DOE EM's environmental cleanup and closure activities. Building upon national capabilities developed from decades of Research and Development in subsurface geosciences, computational and computer science, modeling and applied mathematics, and environmental remediation, the ASCEM initiative will develop an integrated, open-source, high-performance computer modeling system for multiphase, multicomponent, multiscale subsurface flow and contaminant transport. This integrated modeling system will incorporate capabilities for predicting releases from various waste forms, identifying exposure pathways and performing dose calculations, and conducting systematic uncertainty quantification. The ASCEM approach will be demonstrated on selected sites, and then applied to support the next generation of performance assessments of nuclear waste disposal and facility decommissioning across the EM complex. The Multi-Process High Performance Computing (HPC) Simulator is one of three thrust areas in ASCEM. The other two are the Platform and Integrated Toolsets (dubbed the Platform) and Site Applications. The primary objective of the HPC Simulator is to provide a flexible and extensible computational engine to simulate the coupled processes and flow scenarios described by the conceptual models developed using the ASCEM Platform. The graded and iterative approach to assessments naturally generates a suite of conceptual models that span a range of process complexity, potentially coupling hydrological, biogeochemical, geomechanical, and thermal processes. The Platform will use ensembles of these simulations to quantify the associated uncertainty, sensitivity, and risk. The Process Models task within the HPC Simulator focuses on the mathematical descriptions of the relevant physical processes

Sizing of a hot dry rock reservoir from a hydraulic fracturing experiment

Hot dry rock (HDR) reservoirs do not lend themselves to the standard methods of reservoir sizing developed in the petroleum industry such as the buildup/drawdown test. In a HDR reservoir the reservoir is created by the injection of fluid. This process of hydraulic fracturing of the reservoir rock usually involves injection of a large volume (5 million gallons) at high rates (40BPM). A methodology is presented for sizing the HDR reservoir created during the hydraulic fracturing process. The reservoir created during a recent fracturing experiment is sized using the techniques presented. This reservoir is then investigated for commercial potential by simulation of long term power production. 5 refs., 7 figs

Reservoir modeling of the Phase II Hot Dry Rock System

The Phase II system has been created with a series of hydraulic fracturing experiments at the Fenton Hill Hot Dry Rock site. Experiment 2032, the largest of the fracturing operations, involved injecting 5.6 million gallons (21,200m/sup 3/) of water into wellbore EE-2 over the period December 6-9, 1983. The experiment has been modeled using geothermal simulator FEHM developed at Los Alamos National Laboratory. The modeling effort has produced strong evidence of a large highly fractured reservoir. Two long term heat extraction schemes for the reservoir are studied with the model

Simulation of heat transfer in the unsaturated zone

Heat transfer can play an important role in fluid flow near the emplacement site of high-level nuclear waste. The effects on far- field flow can be important in understanding net moisture fluxes above the repository zone. The convection in the unsaturated zone at the Yucca Mountain site was responsible for this movement. If this is so, then the convection could provide a mechanism for drying the rock above the repository zone and thus provide a buffer for heavy rainfall events. In addition, the convection would increase the movement of gaseous radionuclides such as {sup 14}CO{sub 2}, tritiated water vapor, and {sup 129}I (Weeks, 1987). Because of the complexity of the problem, numerical models were required to calculate gas flow and vapor transport at the site. Kipp previously modeled this problem using the code HST3D. This code represents the flow of a single-phase fluid with both heat- and mass-transfer effects included. Water density and partial pressure effects are accounted for by the virtual temperature method. In this paper, the problem was simulated using the code FEHMN, a finite-element heat- and mass-transfer code being developed for the Yucca Mountain Project. The work described in this paper was done in preparation of the upcoming problem to be formulated for the Performance Assessment Calculation Exercise. 5 refs., 9 figs., 1 tab

Non-Darcy flow in geothermal reservoirs

The effects of non-Darcy flow laws are investigated for two geothermal reservoir types: multiphase and Hot Dry Rock (HDR). Long-term thermal behavior is emphasized as short-term pressure transient behavior is addressed in the oil field literature. Comparisons of Darcy and non-Darcy flow laws are made

Calibration of the Site-Scale Saturated Zone Flow Model

The purpose of the flow calibration analysis work is to provide Performance Assessment (PA) with the calibrated site-scale saturated zone (SZ) flow model that will be used to make radionuclide transport calculations. As such, it is one of the most important models developed in the Yucca Mountain project. This model will be a culmination of much of our knowledge of the SZ flow system. The objective of this study is to provide a defensible site-scale SZ flow and transport model that can be used for assessing total system performance. A defensible model would include geologic and hydrologic data that are used to form the hydrogeologic framework model; also, it would include hydrochemical information to infer transport pathways, in-situ permeability measurements, and water level and head measurements. In addition, the model should include information on major model sensitivities. Especially important are those that affect calibration, the direction of transport pathways, and travel times. Finally, if warranted, alternative calibrations representing different conceptual models should be included. To obtain a defensible model, all available data should be used (or at least considered) to obtain a calibrated model. The site-scale SZ model was calibrated using measured and model-generated water levels and hydraulic head data, specific discharge calculations, and flux comparisons along several of the boundaries. Model validity was established by comparing model-generated permeabilities with the permeability data from field and laboratory tests; by comparing fluid pathlines obtained from the SZ flow model with those inferred from hydrochemical data; and by comparing the upward gradient generated with the model with that observed in the field. This analysis is governed by the Office of Civilian Radioactive Waste Management (OCRWM) Analysis and Modeling Report (AMR) Development Plan ''Calibration of the Site-Scale Saturated Zone Flow Model'' (CRWMS M&O 1999a)

Hydro-thermo-mechanical response of a fractured rock block

Hydro-thermo-mechanical effects in fractured rocks are important in many engineering applications and geophysical processes. Modeling these effects is made difficult by the fact that the governing equations are nonlinear and coupled, and the problems to be solved are three dimensional. In this paper we describe a numerical code developed for this purpose. The code is finite element based to allow for complicated geometries, and the time differencing is implicit, allowing for large time steps. The use of state-of-the-art equation solvers has resulted in a practical code. The code is capable of fully three dimensional simulations, however, in this paper we consider only the case of two dimensional heat and mass flow coupled to one dimensional deformation. Partial verification of the code is obtained by comparison with published semianalytical results. Several examples are presented to demonstrate the effects of matrix expansion, due to pore pressure and heating, on fracture opening due to fluid injection. 16 refs., 11 figs

An efficient, three dimensional, fully-coupled hydro-thermo- mechanical simulator

Hydro-thermo-mechanical effects in fractured rocks are important in many oilfield processes. Modeling these effects is made difficult by the fact that the governing equations are nonlinear and coupled, and the problems to be solved are three dimensional. In this paper the authors describe a numerical code developed for this purpose. The code is finite element based to allow for complicated geometries, and the time differencing is implicit, allowing for large time steps. The use of state-of-the-art equation solvers has resulted in a practical code. An example is presented to demonstrate the effects of matrix expansion, due to pore pressure and heating, on fracture opening due to fluid injection, and induced stress changes at a distant well bore. 19 refs., 9 figs., 3 tabs