Computations of permeability of large rock images by dual grid domain decomposition
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Grain size in martian sedimentary rocks can be constrained using point-to-point chemical variabilities in Laser Induced Breakdown Spectroscopy (LIBS) data from the ChemCam instrument on the Mars Science Laboratory (MSL) Curiosity rover. The diameter of each point ablated by the ChemCam laser is in the range of medium to coarse sand in size. Thus, rocks with grains significantly smaller than the laser spot size produce bulk rock compositions at each LIBS point and low point-to-point chemical variability among LIBS points. In contrast, analyses of rocks with grains about the size of the spot or larger contain contributions from individual grains at each point and often have high point-to-point chemical variability. Here the Gini index, a statistical parameter, was used to calculate the point-to-point chemical variability in major-element oxide compositions derived from the ChemCam LIBS data. First, the total range of each LIBS major-element oxide composition was normalized from 0 to 1 across all LIBS observations. Then the Gini index was calculated for each oxide in each LIBS observation. Finally, the Gini indices of each oxide were averaged to derive a Gini index mean score, G MEAN , for each LIBS observation. A correlation between G MEAN and grain size was validated using sedimentary rocks of various grain sizes from the Yellowknife Bay formation and the Pahrump Hills member of the Murray formation in Gale crater. Overall, finer-grained rocks had smaller G MEAN than coarser-grained rocks. To calibrate G MEAN to grain size, grain size estimates based on visual assessment of high-resolution images were compared to G MEAN values for the same targets to create a calibrated scale. This calibrated scale was used to infer the grain size of rocks with unknown grain size. Overall, the grain sizes predicted for rocks with unknown grain size overlapped with those of known grain size from the same units and/or bedrock targets. The grain sizes inferred using the G MEAN based on ChemCam LIBS data are complimentary to those determined from images and both techniques can be used to improve interpretations of the depositional environments of rocks analyzed by Curiosity and future Mars missions with LIBS, such as the Mars 2020 rover.
The surface configuration of pristine layered oxide cathode particles for Li-ion batteries significantly affects the electrochemical behavior, which is generally considered to be a thin rock-salt layer in the surface. Unfortunately, aside from its thin nature and spatial location on the surface, the true structural nature of this surface rock-salt layer remains largely unknown, creating the need to understand its configuration and the underlying mechanisms of formation. Using scanning transmission electron microscopy, we have found a correlation between the surface rock-salt formation and the crystal facets on pristine LiNi0.80Co0.15Al0.05O2 primary particles. It is found that the originally (01 $ \overline{4}\ $) and (003) surfaces of the layered phase result in two kinds of rock-salt reconstructions: the (002) and (111) rock-salt surfaces, respectively. Stepped surface configurations are generated for both reconstructions. The (002) configuration is relatively flat with monoatomic steps while the (111) configuration shows significant surface roughening. Both reconstructions reduce the ionic and electronic conductivity of the cathode, leading to a reduced electrochemical performance.
Despite the numerous studies on changes within the reservoir following CO2 injection and the effects of CO2 release into overlying aquifers, little or no literature is available on the effect of CO2 release on rock between the storage reservoirs and subsurface. To address this knowledge gap, relevant rock materials, temperatures and pressures were used to study mineralogical and elemental changes in this intermediate zone. After rocks reacted with CO2, liquid analysis showed an increase of major elements (e.g., Ca, and Mg) and variable concentrations of potential contaminants (e.g., Sr and Ba); lower concentrations were observed in N2 controls. In experiments with As/Cd and/or organic spikes, representing potential contaminants in the CO2 plume originating in the storage reservoir, most or all of these contaminants were removed from the aqueous phase. SEM and Mössbauer spectroscopy results showed the formation of new minerals and Fe oxides in some CO2-reacted samples, indicating potential for contaminant removal through mineral incorporation or adsorption onto Fe oxides. These experiments show the interactions between the CO2-laden plume and the rock between storage reservoirs and overlying aquifers have the potential to affect the level of risk to overlying groundwater, and should be considered during site selection and risk evaluation.
Explore the source record for details and available documents.
Accurate prediction of strength in rocks with distinct bedding planes requires knowledge of the bedding plane orientation relative to the load direction. Thermal softening adds complexity to the problem since it is known to have significant influence on the strength and strain localization properties of rocks. In this paper, we utilize a recently proposed thermoplastic constitutive model appropriate for rocks exhibiting transverse isotropy in both the elastic and plastic responses to predict their strength and strain localization properties. Recognizing that laboratory-derived strengths can be influenced by material and geometric inhomogeneities of the rock samples, we consider both stress-point and boundary-value problem simulations of rock strength behavior. Both plane strain and 3D loading conditions are considered. Results of the simulations of the strength of a natural Tournemire shale and a synthetic transversely isotropic rock suggest that the mechanical model can reproduce the general U-shaped variation of rock strength with bedding plane orientation quite well. We show that this variation could depend on many factors, including the stress loading condition (plane strain versus 3D), degree of anisotropy, temperature, shear-induced dilation versus shear-induced compaction, specimen imperfections, and boundary restraints, among others.
Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.
Explore the source record for details and available documents.
In order to make predictions of the potential 37 Ar signal from an underground nuclear explosion, a value for the 37 Ar emanation fraction is needed. Here, a system was developed to quantify the emanation fraction from samples ranging in size from powder to small rocks. Assorted rock samples were irradiated with Watt fission spectrum neutrons to produce radioactive argon gas via activation of calcium. The argon produced was extracted on a custom apparatus and cryogenically trapped. The extracted gas was then processed to separate the argon and load it into a proportional counter for counting. The emanation fraction of argon gas from the rock samples was then calculated by comparing the measured and predicted activities.
Rock slope commonly fails due to the shear failure of rock fractures. Shear strength of rock fractures are reduced substantially once the fracture surfaces are mismatched or opened. We propose a new criterion to predict the shear strength of rock fractures in different opening states. The degree of interlocking representing the true asperity contact area is incorporated into the modified model of Saeb and Amadei. The effect of fracture opening on asperity dilation and degradation is separately considered. The transitional stress that is a critical parameter involved in the model is analytically determined based on energy consideration. The new model is validated with experimental results from direct shear tests on synthetic fractures with regular‐shaped asperities. Good agreement between the analytical solution and the experimental data confirms the capacity of the proposed model. Therefore, the model has great potential for assessing the stability of rock slopes where fractures are often opened due to stress relief and engineering disturbances.
We present a reservoir geophysics study, including rock physics modeling and seismic inversion, of a carbon dioxide sequestration site in Southwestern Wyoming, namely the Rock Springs Uplift, and build a petrophysical model for the potential injection reservoirs for carbon dioxide sequestration. Our objectives include the facies classification and the estimation of the spatial model of porosity and permeability for two sequestration targets of interest, the Madison Limestone and the Weber Sandstone. The available dataset includes a complete set of well logs at the location of the borehole available in the area, a set of 110 core samples, and a seismic survey acquired in the area around the well. The proposed study includes a formation evaluation analysis and facies classification at the well location, the calibration of a rock physics model to link petrophysical properties and elastic attributes using well log data and core samples, the elastic inversion of the pre-stack seismic data, and the estimation of the reservoir model of facies, porosity and permeability conditioned by seismic inverted elastic attributes and well log data. In particular, the rock physics relations are facies-dependent and include granular media equations for clean and shaley sandstone, and inclusion models for the dolomitized limestone. The permeability model has been computed by applying a facies-dependent porosity-permeability relation calibrated using core sample measurements. Finally, the study shows that both formations show good storage capabilities. The Madison Limestone includes a homogeneous layer of high-porosity high-permeability dolomite; the Weber Sandstone is characterized by a lower average porosity but the layer is thicker than the Madison Limestone.
Not provided.
At >7 km depths in the Tarim Basin, hydrocarbon reservoirs in Ordovician rocks of the Yijianfang Formation contain large cavities ( c . 10 m or more), vugs, fractures and porous fault rocks. Although some Yijianfang Formation outcrops contain shallow (formed near surface) palaeokarst features, cores from the Halahatang oilfield lack penetrative palaeokarst evidence. Outcrop palaeokarst cavities and opening-mode fractures are mostly mineral filled but some show evidence of secondary dissolution and fault rocks are locally highly ( c . 30%) porous. Cores contain textural evidence of repeated formation of dissolution cavities and subsequent filling by cement. Calcite isotopic analyses indicate depths between c . 220 and 2000 m. Correlation of core and image logs shows abundant cement-filled vugs associated with decametre-scale fractured zones with open cavities that host hydrocarbons. A Sm–Nd isochron age of 400 ± 37 Ma for fracture-filling fluorite indicates that cavities in core formed and were partially cemented prior to the Carboniferous, predating Permian oil emplacement. Repeated creation and filling of vugs, timing constraints and the association of vugs with large cavities suggest dissolution related to fractures and faults. In the current high-strain-rate regime, corroborated by velocity gradient tensor analysis of global positioning system (GPS) data, rapid horizontal extension could promote connection of porous and/or solution-enlarged fault rock, fractures and cavities. Supplementary material: Stable isotopic analyses and the velocity gradient tensor and principal direction and magnitude calculation are available at https://doi.org/10.6084/m9.figshare.c.4946046 Thematic collection: This article is part of the The Geology of Fractured Reservoirs collection available at: https://www.lyellcollection.org/cc/the-geology-of-fractured-reservoirs
To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy's Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.
To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy’s Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.
Time-dependent rock deformation caused by the initiation and growth of fractures leads to the weakening of the rock mass. Understanding the fracturing mechanisms involved in the time-dependent behavior in brittle rocks is very important and to achieve this goal, a systematic series of three types of experiments was performed on double-flawed prismatic Barre granite specimens under unconfined compression. The first series aimed to identify the failure mechanism in the short-term failure mode under monotonic loading, the, second series involved multistage relaxation (constant strain) experiments to analyze the damage at different strain levels, and the third series explored the fracture propagation under multistage creep (constant load) experiments. The spatial and temporal evolution of cracking mechanisms were evaluated using the acoustic emission (AE) and two-dimensional digital image correlation (2D-DIC) techniques to observe the whole crack growth process as well as the accumulated inelastic strain at the specified region of interest. Results suggest that in the case of multistage creep experiments, the time to failure was less compared to the multistage relaxation, when loaded above the crack damage threshold (CD) estimated from the monotonic testing. Furthermore, the frequency magnitude distribution of the AE events generated in the three loading conditions followed the Gutenberg Richter model. A relatively lower b-value was obtained for the creep experiments, indicative of high energy AE events and faster crack growth. In addition, the AE and DIC results also revealed high evolution of tensile cracks at-different stages of creep and relaxation compared to shear and mixed-mode cracks.
Explore the source record for details and available documents.
Abstract not provided