Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constrained control”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

THE CORRELATION BETWEEN H{alpha} EMISSION AND VISUAL MAGNITUDE DURING LONG-TERM VARIATIONS IN CLASSICAL Be STARS

H{alpha} equivalent widths and UBV magnitudes are calculated for Be star disk models that grow in size and/or density with time. We show that these simple models are consistent with the known Be star classes of positive and inverse correlations between long-term variations in H{alpha} and V magnitude as identified by Harmanec. We support the conclusion of Harmanec that the distinction is controlled by the inclination of the disk to the line of sight. We demonstrate that the strength of these correlations, particularly those of an inverse correlation where the system becomes fainter as the H{alpha} emission strength grows, is strongly influenced by the scale height of the inner Be star disk and the extent of the gravitational darkening of the central B star due to rapid rotation. This dependence may allow coordinated spectroscopic and photometric observations to better constrain these poorly known Be star parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Screening for hepatitis D and PEG-Interferon over Tenofovir enhance general hepatitis control efforts in Brazil

Hepatitis D virus (HDV), which requires the presence of hepatitis B virus (HBV), is a deadly yet neglected disease that rapidly leads to liver cancer and disease-induced mortality. This co-dependence creates complex transmission dynamics that make it difficult to predict the efficacy of interventions aimed at HBV and/or HDV control in endemic regions, such as certain municipalities of Brazil, where up to 65% of HBV-infected persons are co-infected. We created a mathematical model that captures the joint transmission dynamics of HBV and HDV, incorporating mother-to-child, sexual and household transmission. With an aim to minimize the number of total infections and disease-induced mortality in 2027, we then determined optimal strategies for Brazil and its sub-regions under a constrained budget, which was dynamically allocated among HBV and HDV screening, HBV and HDV treatment, HBV newborn and adult vaccination, and awareness programs. Three treatment options were considered, namely: Tenofovir, PEGylated-Interferon, and nucleic acid polymers (NAP). The additional cost of HDV screening and the use of a more expensive PEGylated-Interferon are offset by not wasting resources on treating co-infected persons with Tenofovir. The introductory price of NAP treatment must be less than $16,000 per course to become competitive with Tenofovir and PEGylated-Interferon in Brazil. Additional screening for HDV is beneficial, even in a low HBV and HDV endemic regions of Brazil. We recommend PEGylated-Interferon, wherever possible, for both HBV and HDV. If PEGylated-Interferon is not available in abundance, PEGylated-Interferon for co-infections and 4-year Tenofovir treatment for mono-infections is recommended.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1988

The Cretaceous-Tertiary boundary marine extinction and global primary productivity collapse

The extinction of marine phyto-and zoo-plankton across the K-T boundary has been well documented. Such an event may have resulted in decreased photosynthetic fixation of carbon in surface waters and a collapse of the food chain in the marine biosphere. Because the vertical and horizontal distribution of the carbon isotopic composition of total dissolved carton (TDC) in the modern ocean is controlled by the transfer of organic carbon from the surface to deep reservoirs, it follows that a major disruption of the marine biosphere would have had a major effect on the distribution of carbon isotopes in the ocean. Negative carbon isotope excursions have been identified at many marine K-T boundary sequences worldwide and are interpreted as a signal of decreased oceanic primary productivity. However, the magnitude, duration and consequences of this productivity crisis have been poorly constrained. On the basis of planktonic and benthic calcareous microfossil carbon isotope and other geochemical data from DSDP Site 577 located on the Shatsky Rise in the north-central Pacific, as well as other sites, researchers have been able to provide a reasonable estimate of the duration and magnitude of this event.

Zachos, J. C.↗
NASA NTRS1990

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗
NASA NTRS1993

Heliocentric zoning of the asteroid belt by aluminum-26 heating

Variations in petrology among meteorites attest to a strong heating event early in solar system history, but the heat source has remained unresolved. Aluminum-26 has been considered the most likely high-energy, short-lived radionuclide (half-life 0.72 million years) since the discovery of its decay product - excess Mg-26 - in Allende CAI's. Furthermore, observation of relict Mg-26 in an achondritic clast and in feldspars within ordinary chondrites (3,4) provided strong evidence for live Al-26 in meteorite parent bodies and not just in refractory nebular condensates. The inferred amount of Al-26 is consistent with constraints on the thermal evolution of both ordinary and carbonaceous chondrite parent objects up to a few hundred kilometers in diameter. Meteorites can constrain the early thermal evolution of their parent body locations, provided that a link can be established between asteroid spectrophotometric signature and meteorite class. Asteroid compositions are heliocentrically distributed: objects thought to have experienced high metamorphic or even melting temperatures are located closer to the sun, whereas apparently unaltered or mildly heated asteroids are located farther away. Heliocentric zoning could be the result of Al-26 heating if the initial amount of the radionuclide incorporated into planetesimals was controlled by accretion time, which in turn varies with semimajor axis. Analytic expressions for planetary accretion may be integrated to given the time, tau, required for a planetesimal to grow to a specified radius: tau varies as a(sup n), where n = 1.5 to 3 depending on the assumptions about variations in the surface density of the planetesimal swarm. Numerical simulations of planetesimal accretion at fixed semimajor axis demonstrate that variations in accretion time among small planetesimals can be strongly nonlinear depending on the initial conditions and model assumptions. The general relationship with semimajor axis remains valid because it depends only on the initial orbit properties and distribution of the planesimal swarm. In order to demonstrate the basic dependence of thermal evolution on semimajor axis, we parameterized accretion time across the asteroid belt according to tau varies as a(sup n) and calculated the subsequent thermal history. Objects at a specified semimajor axis were assumed to have the same accretion time, regardless of size. We set the initial Al-26/Al-27 ratio = 6 x 10(exp -5) and treated n and tau(sub 0) at a(sub 0) = 3 AU as adjustable parameters. The thermal model included temperature-dependent properties of ice and rock (CM chondrite analog) and the thermodynamic effects of phase transitions.

Grimm, R. E.↗
NASA NTRS2016

Metabolic Potential and Activity in Fluids of the Coast Range Ophiolite Microbial Observatory, California, USA

Metabolic potential and activity associated with hydrogen and carbon monoxide were characterized in fluids sampled from the the Coast Range Ophiolite Microbial Observatory (CROMO). CROMO consists of two clusters of science-dedicated wells drilled to varying depths up to 35m in the actively serpentinizing, Jurassic-age Coast Range Ophiolite of Northern California, along with a suite of pre-existing monitoring wells at the same site. Consistent with the fluid chemistry observed in other serpentinizing systems, CROMO fluids are highly alkaline, with pH up to 12.5, high in methane, with concentrations up 1600 micromolar, and low in dissolved inorganic carbon (DIC), with concentrations of 10's to 100's of micromolar. CROMO is conspicuous for fluid H2 concentrations that are consistently sub-micromolar, orders of magnitude lower than is typical of other systems. However, higher H2 concentrations (10's -100's of micromolar) at an earlier stage of fluid chemical evolution are predicted by, or consistent with: thermodynamic models for fluid chemistry based on parent rock composition equivalent to local peridotite and with water:rock ratio constrained by observed pH; the presence of magnetite at several wt% in CROMO drill cores; and concentrations of formate and carbon monoxide that would require elevated H2 if formed in equilibrium with H2 and DIC. Calculated Gibbs energy changes for reaction of H2 and CO in each of several metabolisms, across the range of fluid composition encompassed by the CROMO wells, range from bioenergetically feasible (capable of driving ATP synthesis) to thermodynamically unfavorable. Active consumption relative to killed controls was observed for both CO and H2 during incubation of fluids from the pre-existing monitoring wells; in incubations of freshly cored solids, consumption was only observed in one sample set (corresponding to the lowest pH) out of three. The specific metabolisms by which H2 and CO are consumed remain to be determined.

Hydrogen↗
NASA NTRSDate not supplied

Towards Erecting Straighter Lightweight Towers on the Moon Using Deployable Guy Wires

This paper reports our static stability test findings for a simple guy wire system to correct the natural lateral deflections of an 8.5m tall, compact deployable composite tower intended to support exploration of lunar permanently shadowed regions by nearby robotic assets. Deployable composite booms with microgravity flight heritage are currently being investigated at NASA Langley Research Center (LaRC) and Massachusetts Institute of Technology(MIT)’s Space Resources Workshop for their potential to be vertically deployed in the lunar gravity field, in support of NASA’s Artemis campaign. These applications include vertical solar arrays and the provision of elevated lines-of-sight to science or engineering payloads on landers and rovers, in support of nearby or distant crewed or robotic assets exploring scientifically interesting and hard to reach areas. Useful elevated payloads include radio repeaters, remote sensing and imaging, navigation and power beaming systems. However, while these lightweight booms have an excellent height to mass ratio, they typically exhibit slight axial curvature upon deployment resulting in appreciable lateral dead-load deflection of the tip mass relative to the tower base. This static deflection increases with tower height and tip mass, not only constraining the value delivered by the tower but also endangering its integrity. To develop a competitive, lightweight deployable composite boom tower, a capability to correct static deflections during and after deployment may be required. This paper presents a pathfinder deployable guy wire stability system for the MIT / LaRC self-erecting composite boom lunar tower that provides real time measurements, maintains tension passively, and can serve as a reconfigurable platform to test new guy wire components, configurations and control algorithms. Using a validated, calibrated photogrammetry system, the natural lateral deflection of the boom tip relative to the boom base at different deployed heights in Earth’s gravity field was recorded. With real-time tension measurements it was found that guy wires can significantly reduce the tip deflection of a deployable composite boom under dead load. Specifically, we found that (1) control capability is greatest where it is needed most, i.e. for the lever arm closest to being opposite the direction of deflection, and (2) for a tower height of at least 8.5 m and arm length of at least 60 cm, a solution of differential tension in all three arms exists and, in principle, provides sufficient control capability to correct or significantly reduce boom deflections. We also found that natural deflections occur almost entirely out-of-plane of the seams of the boom cross-section, which was expected, and that the natural boom tip lateral deflection under dead load upon deployment was ~5% of boom deployed length, unexpectedly exceeding the manufacturing acceptance specification of 1%.Ongoing and future collaborative work between LaRC and MIT includes the further investigation of the unexpected lateral deflection, testing of alternative guy wire system designs at higher tensions and higher deployed heights, as well as trade studies of costs and benefits of an optimized integrated guywire system compared to other types of static stability solutions.

lunar towers↗
DOE OSTI2015

Taxonomy for Common-Cause Failure Vulnerability and Mitigation

Applying current guidance and practices for common-cause failure (CCF) mitigation to digital instrumentation and control (I&C) systems has proven problematic, and the regulatory environment has been unpredictable. The potential for CCF vulnerability inhibits I&C modernization, thereby challenging the long-term sustainability of existing plants. For new plants and advanced reactor concepts, concern about CCF vulnerability in highly integrated digital I&C systems imposes a design burden that results in higher costs and increased complexity. The regulatory uncertainty in determining which mitigation strategies will be acceptable (e.g., what diversity is needed and how much is sufficient) drives designers to adopt complicated, costly solutions devised for existing plants. To address the conditions that constrain the transition to digital I&C technology by the US nuclear industry, crosscutting research is needed to resolve uncertainty, demonstrate necessary characteristics, and establish an objective basis for qualification of digital technology for nuclear power plant (NPP) I&C applications. To fulfill this research need, Oak Ridge National Laboratory is investigating mitigation of CCF vulnerability for nuclear-qualified applications. The outcome of this research is expected to contribute to a fundamentally sound, comprehensive basis to qualify digital technology for nuclear power applications. This report documents the development of a CCF taxonomy. The basis for the CCF taxonomy was generated by determining consistent terminology and establishing a classification approach. The terminology is based on definitions from standards, guides, and relevant nuclear power industry technical reports. The classification approach is derived from identified classification schemes focused on I&C systems and key characteristics, including failure modes. The CCF taxonomy provides the basis for a systematic organization of key systems aspects relevant to analyzing the potential for CCF vulnerability and the suitability of mitigation techniques. Development of an effective CCF taxonomy will help to provide a framework for establishing the objective analysis and assessment capabilities desired to facilitate rigorous identification of fault types and triggers that are the fundamental elements of CCF.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Shape transition of endotaxial islands growth from kinetically constrained to equilibrium regimes

Graphical abstract: - Highlights: • All Fe{sub 13}Ge{sub 8} islands will grow into Ge(0 0 1) substrate at temperatures from 350 to 675 °C. • Shape transition occurred from kinetically constrained to equilibrium regime. • All endotaxial islands can be clarified into two types. • The mechanisms of endotaxial growth and shape transition have been rationalized. - Abstract: A comprehensive study of Fe grown on Ge(0 0 1) substrates has been conducted at elevated temperatures, ranging from 350 to 675 °C. All iron germinide islands, with the same Fe{sub 13}Ge{sub 8} phase, grow into the Ge substrate with the same epitaxial relationship. Shape transition occurs from small square islands (low temperatures), to elongated orthogonal islands or orthogonal nanowires (intermediate temperatures), and then finally to large square orthogonal islands (high temperatures). According to both transmission electron microscopy (TEM) and atomic force microscopy (AFM) investigations, all islands can be defined as either type-I or type-II. Type-I islands usually form at kinetically constrained growth regimes, like truncated pyramids. Type-II islands usually appear at equilibrium growth regimes forming a dome-like shape. Based on a simple semi-quantitative model, type-II islands have a lower total energy per volume than type-I, which is considered as the dominant mechanism for this type of shape transition. Moreover, this study not only elucidates details of endotaxial growth in the Fe–Ge system, but also suggests the possibility of controlled fabrication of temperature-dependent nanostructures, especially in materials with dissimilar crystal structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Forming a database to study reversed magnetic shear from the National Spherical Torus eXperiment using machine learning

Achieving a long-lived reversed magnetic shear (RMS) target plasma in the National Spherical Torus eXperiment Upgrade will require developing various sustainment scenarios. To help with the ongoing plasma control efforts, the development of a new analysis for the motional Stark effect (MSE) diagnostic using a machine learning algorithm, namely, MSE-ML, is described. MSE-ML will be used to identify patterns during RMS discharges, some of which suffer magnetohydrodynamic (MHD) events resulting in current redistribution and monotonic q-profiles. A database consisting of q and magnetic shear profiles is being constructed primarily based on the existing National Spherical Torus eXperiment data with equilibrium reconstructions constrained by the magnetic field pitch angle profile measured using the multi-channel MSE diagnostic. An unsupervised k-means clustering of the data is developed to study the RMS formation as a function of time. The initial clustering from the q-profiles shows significant differences in both amplitude and the duration of the RMS period. As a goal, the clustering results that detect and distinguish shots with substantial and sustained RMS are to be used as a preprocessing step in a supervised algorithm to identify the underlying conditions that lead to long-lasting improved confinement with RMS. Another aim of the MSE-ML study is to identify precursors of RMS-destroying MHD events in either derived data such as the q-profile or directly measured data such as the magnetic field pitch angle profile.

Uzun-Kaymak, I. U. (ORCID:0000000276251493)↗
DOE OSTI2014

A grid of MHD models for stellar mass loss and spin-down rates of solar analogs

Stellar winds are believed to be the dominant factor in the spin-down of stars over time. However, stellar winds of solar analogs are poorly constrained due to observational challenges. In this paper, we present a grid of magnetohydrodynamic models to study and quantify the values of stellar mass loss and angular momentum loss rates as a function of the stellar rotation period, magnetic dipole component, and coronal base density. We derive simple scaling laws for the loss rates as a function of these parameters, and constrain the possible mass loss rate of stars with thermally driven winds. Despite the success of our scaling law in matching the results of the model, we find a deviation between the 'solar dipole' case and a real case based on solar observations that overestimates the actual solar mass loss rate by a factor of three. This implies that the model for stellar fields might require a further investigation with additional complexity. Mass loss rates in general are largely controlled by the magnetic field strength, with the wind density varying in proportion to the confining magnetic pressure B {sup 2}. We also find that the mass loss rates obtained using our grid models drop much faster with the increase in rotation period than scaling laws derived using observed stellar activity. For main-sequence solar-like stars, our scaling law for angular momentum loss versus poloidal magnetic field strength retrieves the well-known Skumanich decline of angular velocity with time, Ω{sub *}∝t {sup –1/2}, if the large-scale poloidal magnetic field scales with rotation rate as B{sub p}∝Ω{sub ⋆}{sup 2}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Trailhead for quantum simulation of SU(3) Yang-Mills lattice gauge theory in the local multiplet basis

Maintaining local interactions in the quantum simulation of gauge field theories relegates most states in the Hilbert space to be unphysical—theoretically benign, but experimentally difficult to avoid. Reformulations of the gauge fields can modify the ratio of physical to gauge-variant states often through classically preprocessing the Hilbert space and modifying the representation of the field on qubit degrees of freedom. This paper considers the implications of representing SU(3) Yang-Mills gauge theory on a lattice of irreducible representations in both a global basis of projected global quantum numbers and a local basis in which controlled-plaquette operators support efficient time evolution. Classically integrating over the internal gauge space at each vertex (e.g., color isospin and color hypercharge) significantly reduces both the qubit requirements and the dimensionality of the unphysical Hilbert space. Initiating tuning procedures that may inform future calculations at scale, the time evolution of one and two plaquettes are implemented on one of IBM’s superconducting quantum devices, and early benchmark quantities are identified. The potential advantages of qudit environments, with either constrained two-dimensional hexagonal or one dimensional nearest-neighbor internal state connectivity, are discussed for future large-scale calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Chance-Constrained AC Optimal Power Flow for Distribution Systems With Renewables

This paper focuses on distribution systems featuring renewable energy sources (RESs) and energy storage systems, and presents an AC optimal power flow (OPF) approach to optimize system-level performance objectives while coping with uncertainty in both RES generation and loads. The proposed method hinges on a chance-constrained AC OPF formulation where probabilistic constraints are utilized to enforce voltage regulation with prescribed probability. A computationally more affordable convex reformulation is developed by resorting to suitable linear approximations of the AC power-flow equations as well as convex approximations of the chance constraints. The approximate chance constraints provide conservative bounds that hold for arbitrary distributions of the forecasting errors. An adaptive strategy is then obtained by embedding the proposed AC OPF task into a model predictive control framework. Finally, a distributed solver is developed to strategically distribute the solution of the optimization problems across utility and customers.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Physical, resource supply, and biological controls on nutrient processing along the river continuum

Nutrient impairment has led to damages to US surface and groundwater systems in excess of 100 billion dollars per year. Therefore, there is a strong need to develop methods to predict the transport, uptake, and export of nutrients along fluvial networks. We present results that are based on a data-driven mechanistic understanding of three factors that largely control nutrient uptake and export: 1) interactions between transport-related processes (mass transfer to metabolically active zones), 2) resource supply dynamics (nutrient concentration, stoichiometric constraints, etc.), and 3) biological controls (microbial community structure and function). Our results were generated from column experiments conducted along the Jemez River-Rio Grande continuum, which spans four orders of magnitude in mean annual discharge, more than 2000 m in altitude, and more than 500 km of stream longitude. Two resource supply injections were performed on each of the columns, i.e., a nitrate only addition, followed by a stoichiometrically ‘balanced’ 106Carbon:16Nitrogen:1Phosphorus addition. We quantified NO3-N uptake kinetics while constraining three variables: stream order, sediment type and type of injection (N vs stoichiometrically ‘balanced’ C:N:P). Following the laboratory nutrient uptake experiments, the columns were destructively sampled and the contents were homogenized to collect subsamples for DNA sequencing. Amplicon analysis was carried out as described by the Earth Microbiome Protocol for 16s and ITS sequencing.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

Inventions Utilizing Microfluidics and Colloidal Particles

Several related inventions pertain to families of devices that utilize microfluidics and/or colloidal particles to obtain useful physical effects. The families of devices can be summarized as follows: (1) Microfluidic pumps and/or valves wherein colloidal-size particles driven by electrical, magnetic, or optical fields serve as the principal moving parts that propel and/or direct the affected flows. (2) Devices that are similar to the aforementioned pumps and/or valves except that they are used to manipulate light instead of fluids. The colloidal particles in these devices are substantially constrained to move in a plane and are driven to spatially order them into arrays that function, variously, as waveguides, filters, or switches for optical signals. (3) Devices wherein the ultra-laminar nature of microfluidic flows is exploited to effect separation, sorting, or filtering of colloidal particles or biological cells in suspension. (4) Devices wherein a combination of confinement and applied electrical and/or optical fields forces the colloidal particles to become arranged into three-dimensional crystal lattices. Control of the colloidal crystalline structures could be exploited to control diffraction of light. (5) Microfluidic devices, incorporating fluid waveguides, wherein switching of flows among different paths would be accompanied by switching of optical signals.

Marr, David W.↗
NASA NTRS2007

Estimating the Inertia Matrix of a Spacecraft

A paper presents a method of utilizing some flight data, aboard a spacecraft that includes reaction wheels for attitude control, to estimate the inertia matrix of the spacecraft. The required data are digitized samples of (1) the spacecraft attitude in an inertial reference frame as measured, for example, by use of a star tracker and (2) speeds of rotation of the reaction wheels, the moments of inertia of which are deemed to be known. Starting from the classical equations for conservation of angular momentum of a rigid body, the inertia-matrix-estimation problem is formulated as a constrained least-squares minimization problem with explicit bounds on the inertia matrix incorporated as linear matrix inequalities. The explicit bounds reflect physical bounds on the inertia matrix and reduce the volume of data that must be processed to obtain a solution. The resulting minimization problem is a semidefinite optimization problem that can be solved efficiently, with guaranteed convergence to the global optimum, by use of readily available algorithms. In a test case involving a model attitude platform rotating on an air bearing, it is shown that, relative to a prior method, the present method produces better estimates from few data.

Acikmese, Behcet↗
DOE OSTI2012

Microelectromechanical (MEM) thermal actuator

Microelectromechanical (MEM) buckling beam thermal actuators are disclosed wherein the buckling direction of a beam is constrained to a desired direction of actuation, which can be in-plane or out-of-plane with respect to a support substrate. The actuators comprise as-fabricated, linear beams of uniform cross section supported above the substrate by supports which rigidly attach a beam to the substrate. The beams can be heated by methods including the passage of an electrical current through them. The buckling direction of an initially straight beam upon heating and expansion is controlled by incorporating one or more directional constraints attached to the substrate and proximal to the mid-point of the beam. In the event that the beam initially buckles in an undesired direction, deformation of the beam induced by contact with a directional constraint generates an opposing force to re-direct the buckling beam into the desired direction. The displacement and force generated by the movement of the buckling beam can be harnessed to perform useful work, such as closing contacts in an electrical switch.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Permeability-controlled migration of induced seismicity to deeper depths near Venus in North Texas

Abstract Migration of seismic events to deeper depths along basement faults over time has been observed in the wastewater injection sites, which can be correlated spatially and temporally to the propagation or retardation of pressure fronts and corresponding poroelastic response to given operation history. The seismicity rate model has been suggested as a physical indicator for the potential of earthquake nucleation along faults by quantifying poroelastic response to multiple well operations. Our field-scale model indicates that migrating patterns of 2015–2018 seismicity observed near Venus, TX are likely attributed to spatio-temporal evolution of Coulomb stressing rate constrained by the fault permeability. Even after reducing injection volumes since 2015, pore pressure continues to diffuse and steady transfer of elastic energy to the deep fault zone increases stressing rate consistently that can induce more frequent earthquakes at large distance scales. Sensitivity tests with variation in fault permeability show that (1) slow diffusion along a low-permeability fault limits earthquake nucleation near the injection interval or (2) rapid relaxation of pressure buildup within a high-permeability fault, caused by reducing injection volumes, may mitigate the seismic potential promptly.

58 GEOSCIENCES↗