Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluids”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Prominin-2 expression increases protrusions, decreases caveolae and inhibits Cdc42 dependent fluid phase endocytosis

Highlights: •Prominin-2 expression induced protrusions that co-localized with lipid raft markers. •Prominin-2 expression decreased caveolae, caveolar endocytosis and increased pCav1. •Prominin-2 expression inhibited fluid phase endocytosis by inactivation of Cdc42. •These endocytic effects can be reversed by adding exogenous cholesterol. •Caveolin1 knockdown restored fluid phase endocytosis in Prominin2 expressing cells. -- Abstract: Background: Membrane protrusions play important roles in biological processes such as cell adhesion, wound healing, migration, and sensing of the external environment. Cell protrusions are a subtype of membrane microdomains composed of cholesterol and sphingolipids, and can be disrupted by cholesterol depletion. Prominins are pentaspan membrane proteins that bind cholesterol and localize to plasma membrane (PM) protrusions. Prominin-1 is of great interest as a marker for stem and cancer cells, while Prominin-2 (Prom2) is reportedly restricted to epithelial cells. Aim: To characterize the effects of Prom-2 expression on PM microdomain organization. Methods: Prom2-fluorescent protein was transfected in human skin fibroblasts (HSF) and Chinese hamster ovary (CHO) cells for PM raft and endocytic studies. Caveolae at PM were visualized using transmission electron microscopy. Cdc42 activation was measured and caveolin-1 knockdown was performed using siRNAs. Results: Prom2 expression in HSF and CHO cells caused extensive Prom2-positive protrusions that co-localized with lipid raft markers. Prom2 expression significantly decreased caveolae at the PM, reduced caveolar endocytosis and increased caveolin-1 phosphorylation. Prom2 expression also inhibited Cdc42-dependent fluid phase endocytosis via decreased Cdc42 activation. Effects on endocytosis were reversed by addition of cholesterol. Knockdown of caveolin-1 by siRNA restored Cdc42 dependent fluid phase endocytosis in Prom2-expressing cells. Conclusions: Prom2 protrusions primarily localize to lipid rafts and recruit cholesterol into protrusions and away from caveolae, leading to increased phosphorylation of caveolin-1, which inhibits Cdc42-dependent endocytosis. This study provides a new insight for the role for prominins in the regulation of PM lipid organization.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Adsorption of simple square-well fluids in slit nanopores: Modeling based on Generalized van der Waals partition function and Monte Carlo simulation

A model based on the Generalized van der Waals partition function is derived to predict the adsorption of square-well fluid in slit pores of any size, the walls of which also have square-well potential. The space inside the pore is divided into several regions based on the extent of the attractive regions generated by the walls. Closed-form expressions of the chemical potentials of the confined fluid in different regions in the pore are obtained. The densities of fluid in different regions are calculated by equalizing the chemical potentials of fluid in those regions to that of the bulk phase. To examine the accuracy of the model, the Grand Canonical Monte Carlo (GCMC) simulation is also conducted. Here we find that the model well captures the effects of the bulk conditions and the properties of adsorbate and pore on the density of adsorbate in the pore. The model is also shown to be able to predict the adsorption of real gases in various activated carbons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Selective transport of light vs. heavy rare earth elements by sulfate/bisulfate complexes in hydrothermal fluids

Here, this study explores the transport of rare earth elements (REE) in acidic sulfate-bearing hydrothermal fluids and the implications for the fractionation of light/heavy REE in critical mineral deposits. The speciation of Nd (light REE) and Yb (heavy REE) sulfate complexes were determined via in situ Raman spectroscopy using fused SiO 2 capillary cells up to 300 ºC at saturated water vapor pressure and in a hydrothermal diamond anvil cell up to 500 ºC and 540 MPa. The REE monosulfate (REESO 4 + ) and REE disulfate (REE(SO 4 ) 2 - ) species are stable below 150 to 250 ºC but become less stable at higher temperatures, particularly the light REE, due to the decreased solubility of REE sulfate solids. At higher pressure, these REE sulfate complexes display an increased stability field up to 400 ℃. The REE bisulfate complex (REEHSO 4 2+ ) was identified with a wide stability field below 400 ºC for the HREE in acidic Yb 2 (SO 4 ) 3 -bearing solutions, whereas in Nd 2 (SO 4 ) 3 -bearing solutions, the LREE bisulfate complex is restricted to below 100 ℃. These results suggest that bisulfate is a previously unrecognized selective ligand for heavy REE transport at low temperature in acidic oxidized hydrothermal fluids. Such fluids are responsible for hydrothermal alteration in many REE deposits. Prediction of phase stabilities across pressure, temperature, and composition space (P-T- x ) is crucial for predicting the role of aqueous REE sulfate complexes in the mobilization of REE in crustal fluids.

58 GEOSCIENCES↗
DOE OSTI2018

A novel algorithm of immersed moving boundary scheme for fluid-particle interactions in DEMLBM

This paper presents an efficient and accurate Immersed Moving Boundary (IMB) algorithm for solving fluid-particle interactions in the framework of the Lattice Boltzmann Method (LBM). Although the IMB scheme has been widely employed in many fluid-particle coupling problems in a wide range of applications, the algorithm of its implementation, especially in identifying both fluid and solid boundary nodes for particles, is seldom reported. Besides, the computational cost of handling fluid-particle coupling is very expensive. To provide a bridge between theory and application and improve the computing efficiency of IMB, a novel boundary tracing procedure and an efficient method for computing the solid nodal ratio using Gaussian quadrature are proposed in this paper. Both accuracy and efficiency of the proposed algorithm are examined by two benchmark tests. It is also found that the IMB scheme are more efficient and stable compared to another widely used the Immersed Boundary Method (IBM) in LBM.

42 ENGINEERING↗
DOE OSTI2017

Special issue of Computers and Fluids in honor of Cecil E. (Chuck) Leith

Here, this special issue of Computers and Fluids is dedicated to Cecil E. (Chuck) Leith in honor of his research contributions, leadership in the areas of statistical fluid mechanics, computational fluid dynamics, and climate theory. Leith's contribution to these fields emerged from his interest in solving complex fluid flow problems--even those at high Mach numbers--in an era well before large scale supercomputing became the dominant mode of inquiry into these fields. Yet the issues raised and solved by his research effort are still of vital interest today.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Sr and O isotopes in western Aleutian seafloor lavas: Implications for the source of fluids and trace element character of arc volcanic rocks

High Mg# andesites and dacites (Mg# = molar Mg/Mg + Fe) from western Aleutian seafloor volcanoes carry high concentrations of Sr (>1000 ppm) that is unradiogenic ( 87 Sr/ 86 Sr < 0.7029) compared to lavas from emergent volcanoes throughout the arc (200–800 ppm Sr, 87 Sr/ 86 Sr >0.7030). Data patterns in plots of 87 Sr/ 86 Sr vs Y/Sr and Nd/Sr imply the existence of an eclogite-melt source component – formed by partial melting of MORB eclogite in the subducting Pacific Plate – which is most clearly expressed in the compositions of western Aleutian andesites and dacites (Nd/Sr and Y/Sr < 0.02) and which dominates the source budget for Sr in volcanic rocks throughout the arc. When viewed in combination with inversely correlated ε Nd and 87 Sr/ 86 Sr, these patterns rule out aqueous fluids as an important source of Sr because mixtures of fluids from altered oceanic crust with depleted mantle and sediment produce compositions with 87 Sr/ 86 Sr higher than in common Aleutian rocks. The unradiogenic nature of Sr in the western Aleutian andesite–dacite end-member may be understood if H 2 O required to drive melting of the subducting oceanic crust is transported in fluids containing little Sr. Mass balance demonstrates that such fluids may be produced by dewatering of serpentinite in the mantle section of the subducting plate. If the eclogite-melt source component is present throughout the Aleutian arc, melting of the subducting plate must extend into minimally altered parts of the sheeted dike section or upper gabbros, at depths >2 km below the paleo-seafloor. Oxygen isotopes in western Aleutian seafloor lavas, which fall within a narrow range of MORB-like values (δ 18 O=5.1–5.7), are also consistent with this model. These results indicate that the subducting Pacific lithosphere beneath the Aleutian arc is significantly hotter than indicated my most thermal models.

58 GEOSCIENCES↗
DOE OSTI2023

Uncertain characterization of reservoir fluids due to brittleness of equation of state regression

Equations of state (EoS) play a central role in modeling the phase equilibrium of fluid mixtures. Their parameterization involves fitting a model to experimental data, i.e., solving a nonlinear, non-convex, multivariate optimization problem. The latter requires one to select design variables, domains of definition for each variable, and weights assigned to individual measurements. We demonstrate that subjective choices of an optimization algorithm and an initial guess also impact the regression process. Consequently, EoS predictions are fundamentally uncertain even after the EoS tuning to a limited set of experimental data points. We demonstrate this observation for two hydrocarbon reservoir fluids, in which five properties of the heaviest carbon fraction are treated as design variables. While all the optimization algorithms and initial guesses match experimental data for the gas and liquid properties, the resulting EoS parameterizations lead to dramatically different predictions of the fluid’s thermophysical behavior in the unsampled pressure and temperature regions. In conclusion, we propose the probabilistic treatment of design variables to quantify the predictive uncertainty of the resulting fluid models.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

"Hidden" hydrothermal technical potential & technoeconomics: Revealing permeability & fluids with more data

Historical hydrothermal estimates have largely relied on temperature or heat flow estimates ignoring the need for natural flowing fluids. More accurate hydrothermal estimates require some indication of permeability and fluids that naturally exist in the subsurface. This paper describes a novel approach that includes proxies of permeability and fluids in hydrothermal estimates by leveraging the relatively data-rich Great Basin. Specifically, nameplate capacities (megawatts) of operating geothermal plants, negative (0 megawatt) locations and 48 geophysical and geologic features are used to used in eXtreme Gradient Boosting (XGBoost) regression to make hydrothermal capacity predictions. Additionally, this work inputs the XGBoost-based hydrothermal predictions into the Renewable Energy Potential (reV) model to quantify technical capacity, its uncertainty and techno-economics. Compared to historical hydrothermal estimates, these predictions adhere to the 37 operating geothermal plants and negative locations. We present a method for subsampling the negative sites to bring the labels into balance that uses the geologic domain knowledge to proportionally represent negatives. Overall, the distributions of the hydrothermal technical capacity and the site levelized cost of energy are respectively much tighter, lower and more accurate than the previous estimates for the Great Basin, as they include geological and geophysical surrogates for permeability and fluids. Percentile (50th and 90th, median and high estimate, respectively) models provide bookends for these metrics.

13 HYDRO ENERGY↗
DOE OSTI2015

Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters

We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Metabolomic changes in follicular fluid induced by soy isoflavones administered to rats from weaning until sexual maturity

Female Wistar rats at 21 days of age were treated with one of three concentrations of soy isoflavones (SIF) (50, 100 or 200 mg/kg body weight, orally, once per day) from weaning until sexual maturity (3 months) in order to evaluate the influence of SIF on ovarian follicle development. After treatment, the serum sex hormone levels and enumeration of ovarian follicles of the ovary were measured. The metabolic profile of follicular fluid was determined using HPLC-MS. Principal component analysis (PCA) and partial least-squares-discriminant analysis (PLS-DA) was used to identify differences in metabolites and reveal useful toxic biomarkers. The results indicated that modest doses of SIF affect ovarian follicle development, as demonstrated by decreased serum estradiol levels and increases in both ovarian follicle atresia and corpora lutea number in the ovary. SIF treatment-related metabolic alterations in follicular fluid were also found in the PCA and PLS-DA models. The 24 most significantly altered metabolites were identified, including primary sex hormones, amino acids, fatty acids and metabolites involved in energy metabolism. These findings may indicate that soy isoflavones affect ovarian follicle development by inducing metabolomic variations in the follicular fluid. - Highlights: ► Modest doses of soy isoflavones (SIF) do affect ovarian follicle development. ► SIF treatment-related metabolic alterations in follicular fluid were found. ► The 24 most significantly altered metabolites were identified.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Reduction of hexavalent chromium by fasted and fed human gastric fluid. II. Ex vivo gastric reduction modeling

To extend previous models of hexavalent chromium [Cr(VI)] reduction by gastric fluid (GF), ex vivo experiments were conducted to address data gaps and limitations identified with respect to (1) GF dilution in the model; (2) reduction of Cr(VI) in fed human GF samples; (3) the number of Cr(VI) reduction pools present in human GF under fed, fasted, and proton pump inhibitor (PPI)-use conditions; and (4) an appropriate form for the pH-dependence of Cr(VI) reduction rate constants. Rates and capacities of Cr(VI) reduction were characterized in gastric contents from fed and fasted volunteers, and from fasted pre-operative patients treated with PPIs. Reduction capacities were first estimated over a 4-h reduction period. Once reduction capacity was established, a dual-spike approach was used in speciated isotope dilution mass spectrometry analyses to characterize the concentration-dependence of the 2nd order reduction rate constants. These data, when combined with previously collected data, were well described by a three-pool model (pool 1 = fast reaction with low capacity; pool 2 = slow reaction with higher capacity; pool 3 = very slow reaction with higher capacity) using pH-dependent rate constants characterized by a piecewise, log-linear relationship. These data indicate that human gastric samples, like those collected from rats and mice, contain multiple pools of reducing agents, and low concentrations of Cr(VI) (< 0.7 mg/L) are reduced more rapidly than high concentrations. The data and revised modeling results herein provide improved characterization of Cr(VI) gastric reduction kinetics, critical for Cr(VI) pharmacokinetic modeling and human health risk assessment. - Highlights: • SIDMS allows for measurement of Cr(VI) reduction rate in gastric fluid ex vivo • Human gastric fluid has three reducing pools • Cr(VI) in drinking water at < 0.7 mg/L is rapidly reduced in human gastric fluid • Reduction rate is concentration- and pH-dependent • A refined PK model is used to characterize inter-individual variability in Cr(VI) gastric reduction capacity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Tracking interface and common curve dynamics for two-fluid flow in porous media

Pore-scale studies of multiphase flow in porous medium systems can be used to understand transport mechanisms and quantitatively determine closure relations that better incorporate microscale physics into macroscale models. Multiphase flow simulators constructed using the lattice Boltzmann method provide a means to conduct such studies, including both the equilibrium and dynamic aspects. Moving, storing, and analyzing the large state space presents a computational challenge when highly-resolved models are applied. We present an approach to simulate multiphase flow processes in which in-situ analysis is applied to track multiphase flow dynamics at high temporal resolution. We compute a comprehensive set of measures of the phase distributions and the system dynamics, which can be used to aid fundamental understanding and inform closure relations for macroscale models. The measures computed include microscale point representations and macroscale averages of fluid saturations, the pressure and velocity of the fluid phases, interfacial areas, interfacial curvatures, interface and common curve velocities, interfacial orientation tensors, phase velocities and the contact angle between the fluid-fluid interface and the solid surface. Test cases are studied to validate the approach and illustrate how measures of system state can be obtained and used to inform macroscopic theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Probing Carbon Mineralization Mechanisms in Pore and Bulk Fluids by Harnessing Architected Calcium Silicates

The ability to synthesize materials with well-controlled pore structures gives us unprecedented control over probing fluid interactions with reactive interfaces and advancing calibrated insights into coupled chemo-morphological interactions. One of the primary challenges in developing crystalline silicate materials lies in achieving ordered pore structures. Existing approaches of producing amorphous mesoporous metal silicates via sol–gel methods and heat-treatment of these materials to produce crystalline phases cause the pore structures in the amorphous phases to collapse. To overcome this challenge, carbon coating of amorphous mesoporous calcium silicate particles is carried out to retain the pore structure, while the material is heated to produce crystalline calcium silicate with calcium sulfate inclusions. The pore diameter in these materials is about 3.9 nm, with a surface area and a pore volume of 28.75 m 2 /g and 0.092 cm 2 /g, respectively. The mechanisms of carbon mineralization are investigated by reacting architected calcium silicate with 1 M Na 2 CO 3 and monitoring the evolution in the structural phases using operando wide-angle X-ray scattering (WAXS) measurements. Formation of stable calcium carbonate polymorph or calcite and metastable calcium carbonate polymorph or vaterite in pore and bulk fluids, respectively, resulting from the reaction between Na 2 CO 3 and CaSiO 3 , are noted. The mechanisms associated with carbon mineralization are delineated using ReaxFF molecular dynamics (MD) simulations. The surface dissolution reaction is initiated by 2H + ions that replace a Ca 2+ ion in the Ca–silicate matrix. Ca 2+ ions in the solution initially react with water to form calcium hydroxide and eventually form calcium (bi)carbonate. A slow and gradual increase in the formation of sodium silicate in the solution resulting from the reactions of silicic acid or the silicon dioxide reaction with sodium hydroxide is noted. When carbon mineralization occurs in environments bearing interfacial fluids, calcite is the dominant calcium carbonate polymorph, as determined using experiments with pore fluids and molecular-scale simulations. In conclusion, these studies provide fundamental insights into the mechanisms underlying the carbon mineralization of calcium silicate informed by experiments and molecular-scale simulations.

Calcium↗
DOE OSTI2024

Molecular Modes Elucidate the Nuclear Magnetic Resonance Relaxation of Viscous Fluids

The Bloembergen, Purcell, and Pound (BPP) theory of nuclear magnetic resonance (NMR) relaxation in fluids dating back to 1948 continues to be the linchpin in interpreting NMR relaxation data in applications ranging from characterizing fluids in porous media to medical imaging (MRI). The BPP theory is founded on assuming molecules are hard spheres with 1 H– 1 H dipole pairs reorienting randomly; assumptions that are severe in light of modern understanding of liquids. Nevertheless, it is intriguing to this day that the BPP theory was consistent with the original experimental data for glycerol, a hydrogen-bonding molecular fluid for which the hard-sphere-rigid-dipole assumption is inapplicable. To better understand this incongruity, atomistic molecular simulations are used to compute 1 H NMR T 1 relaxation dispersion (i.e., frequency dependence) in two contrasting cases: glycerol, and a (non hydrogen-bonding) viscosity standard. At high viscosities, simulations predict distinct functional forms of T 1 for glycerol compared to the viscosity standard, in agreement with modern measurements, yet both in contrast to BPP theory. The cause of these departures from BPP theory is elucidated, without assuming any relaxation models and without any free parameters, by decomposing the simulated T 1 response into dynamic molecular modes for both intramolecular and intermolecular interactions. The decomposition into dynamic molecular modes provides an alternative framework to understand the physics of NMR relaxation for viscous fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A Broadband Laboratory Study of the Seismic Properties of Cracked and Fluid-Saturated Synthetic Glass Media

In this paper, for better understanding of frequency dependence (dispersion) of seismic wave velocities caused by stress-induced fluid flow, broadband laboratory measurements were performed on a suite of synthetic glass media containing both equant pores and thermal cracks. Complementary forced oscillation, resonant bar, and ultrasonic techniques provided access to millihertz-hertz frequencies, ~1 kHz frequency, and ~1 MHz frequency, respectively. The wave speeds or effective elastic moduli and associated dissipation were measured on samples under dry, argon- or nitrogen-saturated, and water-saturated conditions in sequence. The elastic moduli, in situ permeability, and crack porosity inferred from in situ X-ray computed tomography all attest to strong pressure-induced crack closure for differential (confining-minus-pore) pressures <30 MPa, consistent with zero-pressure crack aspect ratios <4 × 10 -4 . The low permeabilities of these materials allow access to undrained conditions, even at subhertz frequencies. The ultrasonically measured elastic moduli reveal consistently higher shear and bulk moduli upon fluid saturation—diagnostic of the saturated-isolated regime. For a glass rod specimen, containing cracks but no pores, saturated-isolated conditions apparently persist to subhertz frequencies—requiring in situ aspect ratios (minimum/maximum dimension) <10 -5 . In marked contrast, the shear modulus measured at subhertz frequencies on a cracked glass bead specimen of 5% porosity, is insensitive to fluid saturation, consistent with the Biot-Gassmann model for the saturated-isobaric regime. Finally, the measured dispersion of the shear modulus approaches 10% over the millihertz-megahertz frequency range for the cracked and fluid-saturated media—implying that laboratory ultrasonic data should be used with care in the interpretation of field data.

58 GEOSCIENCES↗
DOE OSTI2020

Complexity of Fault Rupture and Fluid Leakage in Shale: Insights From a Controlled Fault Activation Experiment

Abstract We observed rupture growth caused by controlled fluid injections at 340‐m depth within a fault zone in the low‐permeability Opalinus Clay in the Mont Terri Underground Research Laboratory (Switzerland). The rupture mechanisms were evaluated using measurements of the three‐component borehole wall displacements and fluid pressure in two sections of the fault zone and located horizontally 3 m apart from each other. One section was set across a secondary segment of the fault and used for stepwise fluid injection intended to trigger rupture growth. The other section was set across the principal shear zone of the fault for monitoring. After stepwise pressure increase up to 5.95 MPa at injection, rupture initiated as slip activation, followed by an overall opening of the fault planes connected to the injection. After 19 s of continued injection, displacements arrived at the monitoring point on the principal shear zone. These displacements are about 2.4 times larger than in the secondary fault segment. Overall, the displacements corresponded to a normal fault activation. About 9 s after the displacement front arrived, a strong pressure increase of 4.17 MPa was measured at the monitoring point, indicating a hydraulic connection had formed along the initially very low permeability fault planes between the injection and the monitoring points. Our analyses highlight that the fault activation is consistent with the state of stress but that injection pressure must be close to the normal stress acting on the fault for permeability to be generated and for fluid leakage to occur.

Guglielmi, Yves↗
DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗