Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Methane (CH4) and carbon dioxide (CO2), the two major components generated from kerogen maturation, are stored dominantly in nanometer-sized pores in shale matrix as (1) a compressed gas, (2) an adsorbed surface species and/or (3) a species dissolved in pore water (H2O). In addition, supercritical CO2 has been proposed as a fracturing fluid for simultaneous enhanced oil/gas recovery (EOR) and carbon sequestration. A mechanistic understanding of CH4-CO2-H2O interactions in shale nanopores is critical for designing effective operational processes. Using molecular simulations, we show that kerogen preferentially retains CO2 over CH4 and that the majority of CO2 either generated during kerogen maturation or injected in EOR will remain trapped in the kerogen matrix. The trapped CO2 may be released only if the reservoir pressure drops below the supercritical CO2 pressure. When water is present in the kerogen matrix, it may block CH4 release. However, the addition of CO2 may enhance CH4 release because CO2 can diffuse through water and exchange for adsorbed methane in the kerogen nanopores.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Shale is characterized by the predominant presence of nanometer-scale (1-100 nm) pores. The behavior of fluids in those pores directly controls shale gas storage and release in shale matrix and ultimately the wellbore production in unconventional reservoirs. Recently, it has been recognized that a fluid confined in nanopores can behave dramatically differently from the corresponding bulk phase due to nanopore confinement (Wang, 2014). CO₂ and H₂O, either preexisting or introduced, are two major components that coexist with shale gas (predominately CH₄) during hydrofracturing and gas extraction. Note that liquid or supercritical CO₂ has been suggested as an alternative fluid for subsurface fracturing such that CO₂ enhanced gas recovery can also serve as a CO₂ sequestration process. Limited data indicate that CO₂ may preferentially adsorb in nanopores (particularly those in kerogen) and therefore displace CH4 in shale. Similarly, the presence of water moisture seems able to displace or trap CH₄ in shale matrix. Therefore, fundamental understanding of CH₄-CO₂-H₂O behavior and their interactions in shale nanopores is of great importance for gas production and the related CO₂ sequestration. This project focuses on the systematic study of CH₄-CO₂-H₂O interactions in shale nanopores under high-pressure and high temperature reservoir conditions. The proposed work will help to develop new stimulation strategies to enable efficient resource recovery from fewer and less environmentally impactful wells.

Kerogen plays a central role in hydrocarbon generation in an oil/gas reservoir. In a subsurface environment, kerogen is constantly subjected to stress confinement or relaxation. The interplay between mechanical deformation and gas adsorption of the materials could be an important process for shale gas production but unfortunately is poorly understood. Using a hybrid Monte Carlo/molecular dynamics simulation, we show here that a strong chemo-mechanical coupling may exist between gas adsorption and mechanical strain of a kerogen matrix. The results indicate that the kerogen volume can expand by up to 5.4% and 11% upon CH4 and CO2 adsorption at 192 atm, respectively. The kerogen volume increases with gas pressure and eventually approaches a plateau as the kerogen becomes saturated. The volume expansion appears to quadratically increase with the amount of gas adsorbed, indicating a critical role of the surface layer of gas adsorbed in the bulk strain of the material. Furthermore, gas uptake is greatly enhanced by kerogen swelling. Swelling also increases the surface area, porosity, and pore size of kerogen. Our results illustrate the dynamic nature of kerogen, thus questioning the validity of the current assumption of a rigid kerogen molecular structure in the estimation of gas-in-place for a shale gas reservoir or gas storage capacity for subsurface carbon sequestration. The coupling between gas adsorption and kerogen matrix deformation should be taken into consideration.

Experimental hydrology data from the Mizunami Underground Research Laboratory in Central Japan have been used to develop a site-scale fracture model and a flow model for the study area. The discrete fracture network model was upscaled to a continuum model to be used in flow simulations. A flow model was developed centered on the research tunnel, and using a highly refined regular mesh. In this study development and utilization of the model is presented. The modeling analysis used permeability and porosity fields from the discrete fracture network model as well as a homogenous model using fixed values of permeability and porosity. The simulations were designed to reproduce hydrology of the modeling area and to predict inflow of water into the research tunnel during excavation. Modeling results were compared with the project hydrology data. Successful matching of the experimental data was obtained for simulations based on the discrete fracture network model.

The Mizunami Underground Research Laboratory is located in the Tono area (Central Japan). Its main purpose is providing a scientific basis for the research and development of technologies needed for deep geological disposal of radioactive waste in fractured crystalline rocks. The current work is focused on the experiments in the research tunnel (500 m depth). The collected tunnel and borehole data were shared with the participants of DEvelopment of COupled models and their VALidation against EXperiments (DECOVALEX) project. This study describes how these data were used to (1) develop the fracture model of the granite rocks around the research tunnel and (2) validate the model.

An ongoing Sandia National Laboratories’ (SNL) research study is evaluating a potential design of an injectable sacrificial material (SM) system that could contain and cool corium ejected from a reactor vessel lower head failure during a potential severe accident involving melting fuel at a commercial light water nuclear reactor (LWR). An injectable system could be installed at any existing LWR, without significant modification to the cavity or to the drywell pedestal region of the plant. The conceptual design under consideration is a passive system. The SM is being optimized to quickly cool the corium mixture while creating gas to form porosity in the solid, such that subsequent water flooding can penetrate the structure and provide additional cooling. The SM would form a barrier and limit corium-concrete interactions. This three-year project takes a joint experimental and computational approach. In this paper, we will first discuss the success of our small-scale experiments conducted on the interactions between the surrogate corium material (SCM) and SM, used to evaluate the injectable concept. A larger experimental study, currently underway, will further validate the injectable concept, with a focus on accurately measuring interactions. This paper details the modeling study and its progress, including modeling the experiments on a surrogate system and extending the model to bench-scale corium flow from validation experiments. The project’s modeling studies will use the SNL engineering code suite SIERRA Mechanics to understand the interaction of injectable SM and molten corium and predict corium spreading. Spreading is modeled using a level set method to track the front in conjunction with a pressure-stabilized finite element method on the fully three-dimensional mass, momentum, and energy conservation equations. Using this diffuse-interface method, the corium spreading front can be tracked and an appropriate pseudo-solidification viscosity models can be implemented to accurately model the corium spreading physics. Finally, an injectable SM delivery system is discussed along with its deployment to the six-common commercial LWR designs currently operating in the United States. At the end of this project, a simplified model based on SIERRA simulations will be developed for implementation into MELCOR, a severe reactor analysis code, developed at SNL for the U.S. Nuclear Regulatory Commission. This will allow us to demonstrate the ability of the injectable SM system to mitigate the ex-vessel corium spreading, provide containment and negate the release of radionuclides.

The porosity of clay aggregates is an important property governing chemical reactions and fluid flow in low-permeability geologic formations and clay-based engineered barrier systems. Pore spaces in clays include interlayer and interparticle pores. Under compaction and dewatering, the size and geometry of such pore spaces may vary significantly (sub-nanometer to microns) depending on ambient physical and chemical conditions. Here we report a molecular dynamics simulation method to construct a complex and realistic clay-like nanoparticle aggregate with interparticle pores and grain boundaries. The model structure is then used to investigate the effect of dewatering and water content on micro-porosity of the aggregates. The results suggest that slow dewatering would create more compact aggregates compared to fast dewatering. Furthermore, the amount of water present in the aggregates strongly affects the particle-particle interactions and hence the aggregate structure. Detailed analyses of particle-particle and water-particle interactions provide a molecular-scale view of porosity and texture development of the aggregates. The simulation method developed here may also aid in modeling the synthesis of nanostructured materials through self-assembly of nanoparticles.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The experimental breeder reactor (EBR-II) used fuel with a layer of sodium surrounding the uranium-zirconium fuel to improve heat transfer. Disposing of EBR-II fuel in a geologic repository without treatment is not prudent because of the potentially energetic reaction of the sodium with water. In 2000, the US Department of Energy (DOE) decided to treat the sodium-bonded fuel with an electrorefiner (ER), which produces metallic uranium product, a metallic waste, mostly from the cladding, and the salt waste in the ER, which contains most of the actinides and fission products. Two waste forms were proposed for disposal in a mined repository; the metallic waste, which was to be cast into ingots, and the ER salt waste, which was to be further treated to produce a ceramic waste form. However, alternative disposal pathways for metallic and salt waste streams may reduce the complexity. For example, performance assessments show that geologic repositories can easily accommodate the ER salt waste without treating it to form a ceramic waste form. Because EBR-II was used for atomic energy defense activities, the treated waste likely meets the definition of transuranic waste. Hence, disposal at the Waste Isolation Pilot Plant (WIPP) in southern New Mexico, may be feasible. This report reviews the direct disposal pathway for ER salt waste and describes eleven tasks necessary for implementing disposal at WIPP, provided space is available, DOE decides to use this alternative disposal pathway in an updated environmental impact statement, and the State of New Mexico grants permission.

Molecular structures of kerogen control hydrocarbon production in unconventional reservoirs. Significant progress has been made in developing model representations of various kerogen structures. These models have been widely used for the prediction of gas adsorption and migration in shale matrix. However, using density functional perturbation theory (DFPT) calculations and vibrational spectroscopic measurements, we here show that a large gap may still remain between the existing model representations and actual kerogen structures, therefore calling for new model development. Using DFPT, we calculated Fourier transform infrared (FTIR) spectra for six most widely used kerogen structure models. The computed spectra were then systematically compared to the FTIR absorption spectra collected for kerogen samples isolated from Mancos, Woodford and Marcellus formations representing a wide range of kerogen origin and maturation conditions. Limited agreement between the model predictions and the measurements highlights that the existing kerogen models may still miss some key features in structural representation. A combination of DFPT calculations with spectroscopic measurements may provide a useful diagnostic tool for assessing the adequacy of a proposed structural model as well as for future model development. This approach may eventually help develop comprehensive infrared (IR)-fingerprints for tracing kerogen evolution.

Shale is characterized by the predominant presence of nanometer-scale (1-100 nm) pores. The behavior of fluids in those pores directly controls shale gas storage and release in shale matrix and ultimately the wellbore production in unconventional reservoirs. Recently, it has been recognized that a fluid confined in nanopores can behave dramatically differently from the corresponding bulk phase due to nanopore confinement (Wang, 2014). CO₂ and H₂O, either preexisting or introduced, are two major components that coexist with shale gas (predominately CH₄) during hydrofracturing and gas extraction. Note that liquid or supercritical CO₂ has been suggested as an alternative fluid for subsurface fracturing such that CO₂ enhanced gas recovery can also serve as a CO₂ sequestration process. Limited data indicate that CO₂ may preferentially adsorb in nanopores (particularly those in kerogen) and therefore displace CH₄ in shale. Similarly, the presence of water moisture seems able to displace or trap CH₄ in shale matrix. Therefore, fundamental understanding of CH₄-CO₂-H₂O behavior and their interactions in shale nanopores is of great importance for gas production and the related CO₂ sequestration. This project focuses on the systematic study of CH₄-CO₂-H₂O interactions in shale nanopores under high-pressure and high temperature reservoir conditions. The proposed work will help to develop new stimulation strategies to enable efficient resource recovery from fewer and less environmentally impactful wells.

Abstract not provided.

Quantitative understanding of reactivity and stability for a chemical species is fundamental to chemistry. The concept has undergone many changes and additions throughout the history of chemistry, stemming from the ideas such as Lewis acids and bases. For a given complexing ligand (Lewis base) and a group of isovalent metal cations (Lewis acids), the stability constants of metal-ligand (ML) complexes can simply correlate to the known properties of metal ions [ionic radii (rMn+), Gibbs free energy of formation (ΔG°f, Mn+), and solvation energy (ΔG°s, Mn+)] by 2.303RT log KML = (α∗MLΔG°f, Mn+ - β∗MLrMn+ + γ∗MLΔG°s, Mn+ - ∗ML), where the coefficients (α∗ML, β∗ML, γ∗ML, and intercept δ∗ML) are determined by fitting the equation to the existing experimental data. Coefficients β∗ML and γ∗ML have the same sign and are in a linear relationship through the origin. Gibbs free energies of formation of cations (ΔG°f, Mn+) are found to be natural indices for the softness or hardness of metal cations, with positive values corresponding to soft acids and negative values to hard acids. The coefficient α∗ML is an index for the softness or hardness of a complexing ligand. Proton (H+) with the softness index of zero is a unique acid that has strong interactions with both soft and hard bases. The stability energy resulting from the acid-base interactions is determined by the term α∗MLΔG°f, Mn+; a positive product of α∗ML and ΔG°f, Mn+ indicates that the acid-base interaction between the metal cation and the complexing ligand stabilizes the complex. The terms β∗MLrMn+ and γ∗MLΔG°s, Mn+, which are related to ionic radii of metal cations, represent the steric and solvation effects of the cations. The new softness indices proposed here will help to understand the interactions of ligands (Lewis bases) with metal cations (Lewis acids) and provide guidelines for engineering materials with desired chemical reactivity and selectivity. The new correlation can also enhance our ability for predicting the speciation, mobility, and toxicity of heavy metals in the earth environments and biological systems.

Abstract not provided.

Abstract not provided.