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.
Clays have many positive attributes leading to their use in waste isolation scenarios. One negative of their use is the formation of clay colloids that can transmit otherwise strongly sorbing contaminants long distances. In this work, we introduce a material combining a layered double hydroxide and a montmorillonite clay. This material is functionally characterized in an advective column under conditions where clay colloid formation is extensive. The material is shown to retard clay colloids once formed, and to help limit the initial formation of clay colloids. Adding this material to a waste barrier will help to limit the erosion of the clays in the barrier and limit the colloidal transport of contaminants.
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). CO2 and H2O, either preexisting or introduced, are two major components that coexist with shale gas (predominately CH4) during hydrofracturing and gas extraction. Note that liquid or supercritical CO2 has been suggested as an alternative fluid for subsurface fracturing such that CO2 enhanced gas recovery can also serve as a CO2 sequestration process. Limited data indicate that CO2 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 CH4 in shale matrix. Therefore, fundamental understanding of CH4-CO2-H2O behavior and their interactions in shale nanopores is of great importance for gas production and the related CO2 sequestration. This project focuses on the systematic study of CH4-CO2-H2O 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.
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.
One of the major challenges of simulating flow and transport in the far field of a geologic repository in crystalline host rock is related to reproducing the properties of the fracture network over the large volume of rock with sparse fracture characterization data. Various approaches have been developed to simulate flow and transport through the fractured rock. The approaches can be broadly divided into Discrete Fracture Network (DFN) and Equivalent Continuum Model (ECM). The DFN explicitly represents individual fractures, while the ECM uses fracture properties to determine equivalent continuum parameters. We compare DFN and ECM in terms of upscaled observed transport properties through generic fracture networks. The major effort was directed on making the DFN and ECM approaches similar in their conceptual representations. This allows for separating differences related to the interpretation of the test conditions and parameters from the differences between the DFN and ECM approaches. The two models are compared using a benchmark test problem that is constructed to represent the far field (1 × 1 × 1 km3) of a hypothetical repository in fractured crystalline rock. The test problem setting uses generic fracture properties that can be expected in crystalline rocks. The models are compared in terms of the: 1) effective permeability of the domain, and 2) nonreactive solute breakthrough curves through the domain. The principal differences between the models are mesh size, network connectivity, matrix diffusion and anisotropy. We demonstrate how these differences affect the flow and transport. We identify the factors that should be taken in consideration when selecting an approach most suitable for the site-specific conditions.
Jove Colon, Carlos F.; Wang, Yifeng; Hadgu, Teklu H.; Zheng, Liange; Rutqvist, Jonny; Xu, Hao; Kim, Kunhwi; Voltolini, Marco; Cao, Xiaoyuan; Fox, Patricia; Nico, Peter S.; Caporuscio, Florie A.; Norskog, Katherine E.; Zavarin, Mavrik; Wolery, Thomas J.; Atkins-Duffin, Cindy; Jerden, James; Gattu, Vineeth K.; Ebert, William; Buck, Edgar C.; Wittman, Richard S.
The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key advances in experimental and modeling aspects of chemical and physical phenomena towards the long-term safety assessment of nuclear waste disposition in deep clay/shale/argillaceous rock. Experimental activities on clay barrier interactions with fluids and radionuclides provide the much needed knowledge to evaluate engineered barrier system (EBS) performance. Thermal-Hydrological-Mechanical-Chemical (THMC) model development of clay provides a rigorous simulation platform to assess the complex dynamic behavior of engineered and natural barrier materials in response to coupled process phenomena induced by heat-generating nuclear waste. This report describes the ongoing disposal R&D efforts on the advancement and refinement of coupled THMC process models, hydrothermal experiments and geochemical modeling of on barrier material (clay/metal) interactions, spent fuel and canister material degradation, radiolytic phenomena and UO2 degradation, and thermodynamic database development. These play an important role to the evaluation of sacrificial zones as part of the EBS exposure to thermally-driven chemical and transport processes. Clay-zeolite phase equilibria play a key role in the mineralogical transformations of clay barrier conducive to loss in swelling properties but also in controlling H20 uptake/release through hydration/dehydration reactions. The result is volume changes can affect the interface / bulk phase porosities, transport, and the mechanical (stress) state of the bentonite barrier. Characterization studies on barrier samples (bentonite/cement) from controlled tests at underground research laboratories (URLs) provide key insights into barrier materials interactions at EBS interfaces. Spent fuel degradation modeling coupled with canister and cladding corrosion effects demonstrate the strong influence of H2 generation on the source term.
The objective of the Crystalline Disposal R&D control account is to advance our understanding of long-term disposal of used fuel in crystalline rocks and to develop necessary experimental and computational capabilities to evaluate various disposal concepts in such media.
Active participation in international R&D is crucial for achieving the Spent Fuel Waste Science & Technology (SFWST) long-term goals of conducting “experiments to fill data needs and confirm advanced modeling approaches” and of having a “robust modeling and experimental basis for evaluation of multiple disposal system options” (by 2020). DOE’s Office of Nuclear Energy (NE) has developed a strategic plan to advance cooperation with international partners. The international collaboration on the evaluation of crystalline disposal media at Sandia National Laboratories (SNL) in FY17 focused on the collaboration through the Development of Coupled Models and their Validation against Experiments (DECOVALEX-2019) project. The DECOVALEX project is an international research and model comparison collaboration, initiated in 1992, for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. SNL has been participating in three tasks of the DECOVALEX project: Task A. Modeling gas injection experiments (ENGINEER), Task C. Modeling groundwater recovery experiment in tunnel (GREET), and Task F. Fluid inclusion and movement in the tight rock (FINITO).
The experimental breeder reactor (EBR-II) used fuel with a layer of sodium surrounding the uranium-zirconium fuel to improve heat transfer. Disposing of this EBR-II used 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 decided to treat the EBR-II sodium-bonded used fuel in an electrorefiner (ER), which produces a metallic waste, mostly from the cladding. The salt remaining in the ER contains most of the actinides and fission products. Two baseline 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 are being investigated that may reduce the complexity. For example, performance assessments show that both mined repositories in salt and deep boreholes in basement crystalline rock can easily accommodate the ER salt waste without treating it to form a ceramic waste form. Hence the focus of a direct disposal option, as described herein, is now on the feasibility of packaging the ER salt waste in the near term such that it can be transported to a repository in the future without repackaging. A vessel for direct disposal of ER salt waste has been previously proposed, designed, and a prototype manufactured based on desirable features for use in the hot cell. The reported analysis focused on the feasibility of transporting this proposed vessel and whether any issues would suggest that a smaller or larger size is more appropriate. Specifically, three issues are addressed (1) shielding necessary to reduce doses to acceptable levels; (2) the criticality potential and the ease which it can be shown to be inconsequential, and (3) temperatures of the containers in relation to acceptable cask limits. The generally positive results demonstrate that direct disposal of ER in the proposed packaging is feasible without the need to secure funding to modify the facility.
The migration and trapping of supercritical CO2 (scCO2) in geologic carbon storage is strongly dependent on the geometry and wettability of the pore network in the reservoir rock. During displacement, resident fluids may become trapped in the pits of a rough pore surface forming an immiscible two-phase fluid interface with the invading fluid, allowing apparent slip flow at this interface. We present a two-phase fluid dynamics model, including interfacial tension, to characterize the impact of mineral surface roughness on this slip flow. We show that the slip flow can be cast in more familiar terms as a contact-angle (wettability)-dependent effective permeability to the invading fluid, a nondimensional measurement which relates the interfacial slip to the pore geometry. The analysis shows the surface roughness-induced slip flow can effectively increase or decrease this effective permeability, depending on the wettability and roughness of the mineral surfaces. Configurations of the pore geometry where interfacial slip has a tangible influence on permeability have been identified. The results suggest that for large roughness features, permeability to CO2 may be enhanced by approximately 30% during drainage, while the permeability to brine during reimbibition may be enhanced or diminished by 60%, depending on the contact angle with the mineral surfaces and degrees of roughness. For smaller roughness features, the changes in permeability through interfacial slip are small. As a result, a much larger range of effective permeabilities are suggested for general fluid pairs and contact angles, including occlusion of the pore by the trapped phase.
The migration and trapping of supercritical CO2 (scCO2) in geologic carbon storage is strongly dependent on the geometry and wettability of the pore network in the reservoir rock. During displacement, resident fluids may become trapped in the pits of a rough pore surface forming an immiscible two-phase fluid interface with the invading fluid, allowing apparent slip flow at this interface. We present a two-phase fluid dynamics model, including interfacial tension, to characterize the impact of mineral surface roughness on this slip flow. We show that the slip flow can be cast in more familiar terms as a contact-angle (wettability)-dependent effective permeability to the invading fluid, a nondimensional measurement which relates the interfacial slip to the pore geometry. The analysis shows the surface roughness-induced slip flow can effectively increase or decrease this effective permeability, depending on the wettability and roughness of the mineral surfaces. Configurations of the pore geometry where interfacial slip has a tangible influence on permeability have been identified. The results suggest that for large roughness features, permeability to CO2 may be enhanced by approximately 30% during drainage, while the permeability to brine during reimbibition may be enhanced or diminished by 60%, depending on the contact angle with the mineral surfaces and degrees of roughness. For smaller roughness features, the changes in permeability through interfacial slip are small. A much larger range of effective permeabilities are suggested for general fluid pairs and contact angles, including occlusion of the pore by the trapped phase.