Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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.

At least 1,045 records · Page 58

The Fluids Integrated Rack and Light Microscopy Module Integrated Capabilities

The Fluids Integrated Rack (FIR), a facility class payload, and the Light Microscopy Module (LMM), a subrack payload, are scheduled to be launched in 2005. The LMM integrated into the FIR will provide a unique platform for conducting fluids and biological experiments on ISS. The FIR is a modular, multi-user scientific research facility that will fly in the U.S. laboratory module, Destiny, of the International Space Station (ISS). The first payload in the FIR will be the Light Microscopy Module (LMM). The LMM is planned as a remotely controllable, automated, on-orbit microscope subrack facility, allowing flexible scheduling and control of fluids and biology experiments within the FIR. Key diagnostic capabilities for meeting science requirements include video microscopy to observe microscopic phenomena and dynamic interactions, interferometry to make thin film measurements with nanometer resolution, laser tweezers for particle manipulation, confocal microscopy to provide enhanced three-dimensional visualization of structures, and spectrophotometry to measure photonic properties of materials. The LMM also provides experiment sample containment for frangibles and fluids. This paper will provide a description of the current FIR and LMM designs, planned capabilities and key features. In addition a brief description of the initial five experiments planned for LMM/FIR will be provided.

Motil, Susan M.↗

Isolated Fluid Oxygen Drop Behavior in Fluid Hydrogen at Rocket Chamber Pressures

A model has been developed for the behavior of an isolated fluid drop of a single compound immersed into another compound in finite, quiescent surroundings at supercritical conditions. The model is based upon fluctuation theory which accounts for both Soret and Dufour effects in the calculation of the transport matrix relating molar and heat fluxes to the transport properties and the thermodynamic variables.

Fluid↗

Potential of a New Lunar Surface Radiator Concept for Hot Lunar Thermal Environments

The optimum radiator configuration in hot lunar thermal environments is one in which the radiator is parallel to the ground and has no view to the hot lunar surface. However, typical spacecraft configurations have limited real estate available for top-mounted radiators, resulting in a desire to use the spacecraft s vertically oriented sides. Vertically oriented, flat panel radiators will have a large view factor to the lunar surface, and thus will be subjected to significant incident lunar infrared heat. Consequently, radiator fluid temperatures will need to exceed approx.325 K (assuming standard spacecraft radiator optical properties) in order to provide positive heat rejection at lunar noon. Such temperatures are too high for crewed spacecraft applications in which a heat pump is to be avoided. A recent study of vertically oriented radiator configurations subjected to lunar noon thermal environments led to the discovery of a novel radiator concept that yielded positive heat rejection at lower fluid temperatures. This radiator configuration, called the Upright Lunar Terrain Radiator Assembly (ULTRA), has exhibited superior performance to all previously analyzed concepts in terms of heat rejection in the lunar noon thermal environment. A key benefit of the ULTRA is the absence of louvers or other moving parts and its simple geometry. Analysis of the ULTRA for a lunar extravehicular activity (EVA) portable life support system (PLSS) is shown to provide moderate heat rejection, on average, at all solar incident angles assuming an average radiator temperature of 294 K, whereas prior concepts exhibited insignificant heat rejection or heat absorption at higher incident angles. The performance of the ULTRA for a lunar lander is also discussed and compared to the performance of a vertically oriented, flat panel radiator at various lunar latitudes.

Ochoa, Dustin A.↗

Anomalous and Ungrouped Carbonaceous Chondrites in the US Antarctic Meteorite Collection and their Potential Relevance to Ryugu and Bennu

With two different carbonaceous asteroid sample return missions in full swing, attention has focused on what connections can be made between the asteroid samples and the wide range of carbonaceous chondrite meteorites in worldwide collections. The US Antarctic meteorite collection contains nearly 1000 carbonaceous chondrites of various types including many in well-established groups as well as ungrouped, unusual or anomalous groups [1]. Some of the latter have been included in, or are possibly related to, recently proposed new carbonaceous chondrite classifications – CA and CY chondrites [2,3]. In addition to these, there are numerous ungrouped samples that have properties intermediate between established groups (like CM and CO; [4]), distinct from any other groups [5], or have anomalous properties that might be attributable to parent body processes such as heating, fluid interaction, or impacts [6]. Some CM anomalous or ungrouped chondrites share spectral features with Ryugu, which has a small hydration peak arguably due to hydrated minerals left after either impact heating or shock in carbonaceous chondrites [7]. PCA 91008, PCA 02012, GRO 95566, and LEW 85311 are all CMs that have experienced heating or metamorphism that may be due to impacts, solar radiation, or radiogenic decay [6]. These samples all have low H contents, C/H (bulk), and low 17O [6,8] and may hold clues to understanding the mineralogy of Ryugu, or interpretations of its spectral properties. WIS 91600, on the other hand, appears to be related to several other highly altered CM [6], and shares properties with the newly proposed CY chondrites [6]. An understanding of this grouplet will also aid in the interpretation of Bennu samples which have strong hydration features. In addition, ungrouped carbonaceous chondrites may provide valuable insights into the aqueous alteration potentially recorded in Bennu and Ryugu samples; such as the relatively moderate aqueous alteration recognized and dated at 4-5 myr after CAIs in MAC 88107 (C2-ungrouped) by [9] to the extensive aqueous alteration apparent in MIL 090292 (C1-ungrouped) [10]. Finally, LON 94101/94102 is a brecciated CM chondrite with numerous lithologies. Its appearance is similar to some of the brecciated lithologies visible at the surface of Bennu and Ryugu [11]. Although CM chondrites with multiple lithologies are common (e.g., [12]), the lithologies are often difficult to resolve at the hand specimen scale and only after some detailed e-beam characterization are the subtle lithologic differences evident (e.g., [13]). LON 94102 contains visibly distinct clasts at the hand specimen scale, relatively rare for CM breccias. The largest CM chondrites by mass from the U.S. Antarctic meteorite collection may provide insights into the extent of brecciation and heterogeneity within more typical CM chondrite-like source asteroids. For example, recent curation CT scans [14] show possible clasts within in a subsplit of ALH 83100 (CM1/2) which is the largest CM chondrite or CM chondrite pairing group in the U.S. Antarctic collection with an original mass of 3.019 kg. Similar work on additional subsplits or meteorites may aid our understanding of heterogeneity within CM chondrites from Antarctica. Initial classification of CM chondrites in U.S. Antarctic meteorite collection includes preliminary pairing when petrographically similar CM chondrites that have been previously found in the same field area. The two largest CM chondrite preliminary pairing groups by mass are the ALH 83102 (CM2) pairing group with an original mass of 2.554 kg and the EET 96005 (CM2) pairing group with an original mass of 1.125 kg. However, preliminary pairing groups assigned at classification—intrinsically— do not include stones that are petrographically distinct. Rigorous pairing group studies of CM chondrites from specific field sites are needed to investigate if there is unrecognized heterogeneity in CM chondrites from Antarctica that may represent common asteroid impactors (pre-atmospheric entry asteroids/meteoroids). Detailed pairing studies have the potential to recognize initial stones with multiple lithologies and investigate if there are stones that only sample one of those respective lithologies in the collection. These groups of heated CMs, extensively hydrated CMs, intermediate between CM and CO chondrites, and brecciated samples are all potentially relevant to Bennu and Ryugu samples, where heating, hydration, and brecciation have all come into play. These bodies might also be comprised of material intermediate to CM and CO chondrites, or at least distinct from the well-established CC groups. These small and unusual groups of carbonaceous chondrites may help to unlock new information about early solar system processes and aid in the understanding of the evolution of these carbonaceous asteroids.

meteorites↗

Models of the earth's core

Combined inferences from seismology, high-pressure experiment and theory, geomagnetism, fluid dynamics, and current views of terrestrial planetary evolution lead to models of the earth's core with five basic properties. These are that core formation was contemporaneous with earth accretion; the core is not in chemical equilibrium with the mantle; the outer core is a fluid iron alloy containing significant quantities of lighter elements and is probably almost adiabatic and compositionally uniform; the more iron-rich inner solid core is a consequence of partial freezing of the outer core, and the energy release from this process sustains the earth's magnetic field; and the thermodynamic properties of the core are well constrained by the application of liquid-state theory to seismic and labroatory data.

Stevenson, D. J.↗

The preparation of new perfluoroether fluids exhibiting excellent thermal-oxidative stabilities

A series of low molecular weight perfluoroalkylethers (PFAE) were synthesized by direct fluorination. Viscosity-temperature properties and oxidation stabilities were determined. Viscosity-temperature correlations indicated that increases in branching and increases in the size of the branching substituent caused a deterioration in viscometric properties (i.e., an increase in ASTM slope). In addition, increasing the ratio of carbon to oxygen in these compounds also increased the ASTM slope. Preliminary oxidation stability tests indicated that highly branched PFAE fluids. (i.e., those containing quaternary carbons) may be less stable than either those containing a single trifluoromethyl pendant group or those containing no branching at all.

Jones, W. R., Jr.↗

Adaptive Mesh Refinement Large Eddy Simulation of the Supercritical Carbon Dioxide Round Turbulent Jet

Supercritical carbon dioxide (sCO2) is of interest to a range of engineering problems, including carbon capture, utilization, and storage (CCUS) as well as advanced cycles for power generation. Non-ideal variations in physical properties of sCO2 impact the physics of these systems. In this study, we simulate turbulent sCO2 jets to gain a better understanding of these physics.We use a second order finite volume method with adaptive mesh refinement as implemented in the first-principles simulation code PeleC to perform a Large Eddy Simulation (LES) of three turbulent jets of sCO2. Additionally, we use the Soave-Redlich-Kwong equation of state to close the system and examine the impact of a cubic equation of state on the turbulent flow physics. We look at velocity and Reynolds stress profiles at different downstream locations for three cases in which the temperature of the jet andthat of the ambient fluid differ in order to capture the effects of widely varying thermal properties in the pseudocritical region. These results are then contrasted with established theory for ideal gas jets.

adaptive mesh refinement↗

Imaging of a fluid injection process using geophysical data — A didactic example

In many subsurface industrial applications, fluids are injected into or withdrawn from a geologic formation. It is of practical interest to quantify precisely where, when, and by how much the injected fluid alters the state of the subsurface. Routine geophysical monitoring of such processes attempts to image the way that geophysical properties, such as seismic velocities or electrical conductivity, change through time and space and to then make qualitative inferences as to where the injected fluid has migrated. The more rigorous formulation of the time-lapse geophysical inverse problem forecasts how the subsurface evolves during the course of a fluid-injection application. Using time-lapse geophysical signals as the data to be matched, the model unknowns to be estimated are the multiphysics forward-modeling parameters controlling the fluid-injection process. Properly reproducing the geophysical signature of the flow process, subsequent simulations can predict the fluid migration and alteration in the subsurface. The dynamic nature of fluid-injection processes renders imaging problems more complex than conventional geophysical imaging for static targets. This work intents to clarify the related hydrogeophysical parameter estimation concepts.

58 GEOSCIENCES↗

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗

A hybrid machine-learning approach for analysis of methane hydrate formation dynamics in porous media with synchrotron CT imaging

Fast multi-phase processes in methane hydrate bearing samples pose a challenge for quantitative micro-computed tomography study and experiment steering due to complex tomographic data analysis involving time-consuming segmentation procedures. This is because of the sample's multi-scale structure, which changes over time, low contrast between solid and fluid materials, and the large amount of data acquired during dynamic processes. Here, a hybrid approach is proposed for the automatic segmentation of tomographic data from time-resolved imaging of methane gas-hydrate formation in sandy granular media, which includes a deep-learning 3D U-Net model. To prepare a training dataset for the 3D U-Net, a technique to automate data labeling based on sample-specific information about the mineral matrix immobility and occasional fluid movement in pores is proposed. Automatic segmentation allowed for studying properties of the hydrate growth in pores, as well as dynamic processes such as incremental flow and redistribution of pore brine. Results of the quantitative analysis showed that for typical gas-hydrate stability parameters (100 bar methane pressure, 7°C temperature) the rate of formation is slow (less than 1% per hour), after which the surface area of contact between brine and gas increases, resulting in faster formation (2.5% per hour). Hydrate growth reaches the saturation point after 11 h of the experiment. Finally, the efficacy of the proposed segmentation scheme in on-the-fly automatic data analysis and experiment steering with zooming to regions of interest is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytical Study and Experimental Verification of Shear-Thinning Ink Flow in Direct Ink Writing Process

Abstract Direct ink writing (DIW) process is a facile additive manufacturing technology to fabricate three-dimensional (3D) objects with various materials. Its versatility has attracted considerable interest in academia and industry in recent years. As such, upsurging endeavors are invested in advancing the ink flow behaviors in order to optimize the process resolution and the printing quality. However, so far, the physical phenomena during the DIW process are not revealed in detail, leaving a research gap between the physical experiments and its underlying theories. Here, we present a comprehensive analytical study of non-Newtonian ink flow behavior during the DIW process. Different syringe-nozzle geometries are modeled for the comparative case studies. By using the computational fluid dynamics (CFD) simulation method, we reveal the shear-thinning property during the ink extrusion process. Besides, we study the viscosity, shear stress, and velocity fields, and analyze the advantages and drawbacks of each syringe-nozzle model. On the basis of these investigations and analyses, we propose an improved syringe-nozzle geometry for stable extrusion and high printing quality. A set of DIW printing experiments and rheological characterizations are carried out to verify the simulation studies. The results developed in this work offer an in-depth understanding of the ink flow behavior in the DIW process, providing valuable guidelines for optimizing the physical DIW configuration toward high-resolution printing and, consequently, improving the performance of DIW-printed objects.

Engineering↗

Materials Assessment for FLiBe Fusion Blankets

This evaluation was prepared at the request of DOE-FES and expands on an earlier evaluation prepared by Wilson, Wiffen and Keiser in 2015 (attached as Appendices A and B). Additional useful information for materials in fission molten salt reactor systems that use FLiBe can be found in a prior NRC report [Busby 2019]. The purpose of this report is to review the considerations that will determine the fusion blanket structural materials that will be suitable for use in a FLiBe breeder-coolant fusion power plant (FPP). Identifying a fusion-relevant structural material that is sufficiently compatible with FLiBe is a key issue to be addressed. The constraints and trade-offs made in structural materials selection include the usual power generation requirements such as mechanical properties, fabricability and durability, which assumes reasonable compatibility with the working fluid. For an FPP, additional constraints include neutron economy (breeding ratio), resistance to irradiation damage and low activation (i.e. waste disposal) characteristics.

36 MATERIALS SCIENCE↗

Produced Fluid Induced Mineralogy and Elemental Alterations of Caney Shale, Southern Oklahoma

ABSTRACT This study involves batch reactor experiments and subsequent analyses of samples from Caney Shale in the Ardmore Basin of South-Central Oklahoma. Samples include mainly rock cores and cuttings recovered from two wells respectively drilled vertically through and horizontally across the Caney Shale. Mineralogical compositions are obtained by X-Ray Diffraction (XRD) measurements whilst microstructure and elemental distribution are acquired by Scanning Electron Microscopy/Energy Dispersive Spectroscopy (SEM/EDS) respectively. Batch experiments are then conducted using selected rock samples and produced fluid from the Caney Formation. Deionized water is also reacted with some samples to serve as standard. Experiments are conducted at 95°C and ambient pressure for 7 and 30 days to assess the geochemical rock-fluid interactions. Results show rock mineralogical compositions are predominantly quartz, feldspar, carbonates, and clay with minor pyrite. Post-experimental mineralogical changes observed in samples include increased amorphous entities especially within the clay portions of XRD plots and dissolution of feldspar and carbonate minerals and formation of new mineral phases, mostly clays and salts. These are corroborated by EDS elemental analyses which show decreased elemental compositions. The implications of reactions mentioned above include but not limited to, scale formation, clay fines migration and shale softening all of which pose significant permeability impairment on formation over time. INTRODUCTION Shale reservoirs account for a large share of unconventional reservoirs in the world (Lyu et al., 2015). However, ultra-low permeability and high clay compositions pose significant challenges when producing from these reservoirs (Dawuda and Srinivasan, 2022, 2023). Producing from these reservoirs therefore requires horizontal drilling and hydraulic fracturing technologies which have proven their efficacy in generating substantial permeability in reservoirs to ensure production (Fujian et al., 2019; Liu et al., 2018). Even after expensive horizontal drilling and hydraulic fracturing, geochemical reactions between engineered fluids and formation leads to fracture constriction and adversely impact petrophysical properties (permeability and porosity) of the reservoir. These technologies are therefore under constant development and improvement in various aspects to ensure fine tuning for specific reservoirs. Under present conditions, much of the hydrocarbon reserves in unconventional shale reservoirs are left unproduced due to rapid decline in permeability following resumption to production after hydraulic fracturing. It is therefore essential to understand the range of geochemical reactions that cause rapid depletion of permeability after hydraulic fracturing and apply these to each shale reservoir to ensure substantial recovery rates.

Awejori, G. A.↗