Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “injection rate”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2022

Estimation of China Fusion Engineering Test Reactor performance and burning fraction in different pellet fueling scenarios by a multi-species radial transport model

Tritium self-sufficiency in future deuterium–tritium fusion reactors is a crucial challenge. As an engineering test reactor, the China Fusion Engineering Test Reactor requires a burning fraction of 3% for the goal to test the accessibility to the future fusion plant. To self-consistently simulate burning plasmas with profile changes in pellet injection scenarios and to estimate the corresponding burning fraction, a one-dimensional multi-species radial transport model is developed in the BOUT++ framework. Here several pellet-fueling scenarios are then tested in the model. The results show that the increased fueling depth improves the burning fraction by particle confinement improvement and fusion power increase. Nevertheless, by increasing the depth, the pellet cooling-down may significantly lower the temperature in the core region. Taking the density perturbation into consideration, the reasonable parameters of the fueling scenario in these simulations are estimated as pellet radius ${r}_{{\rm{p}}}=3\,{\rm{mm}},$ injection rate $=\,4\,\mathrm{Hz},$ and pellet injection velocity $=\,1000\mbox{-}2000\,{\rm{m}}\,{{\rm{s}}}^{-1}$ without drift or $450\,{\rm{m}}\,{{\rm{s}}}^{-1}$ with high-field-side drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A micromechanical hydro-mechanical-damage coupled model for layered rocks considering multi-scale structures

An anisotropic Hydro-Mechanical-Damage coupled model was established in the framework of thermodynamics for saturated layered rocks, which contain arbitrarily-distributed cracks of smaller scales and orientated bedding planes of much larger scales. Here, the proposed model considers the anisotropic damage growth and frictional sliding of cracks and shear sliding and dilatancy of bedding planes, as well as the interaction between cracks and bedding planes. Numerical implementation of the nonlinear coupled model was developed based on TOUGHREACT, which is a well-established simulator for multiphase fluid flow and reactive transport analysis. The proposed model was first validated with laboratory water injection test results on rock blocks, and then used to study the effects of injection rates, in situ stresses and bedding planes on fluid injection processes at the field scale, demonstrating the importance to consider the variation of multi-scale structures for better understanding the coupled Hydro-Mechanical responses in layered rocks.

58 GEOSCIENCES↗
DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Synergies between H, He and radiation damage in dual and triple ion irradiation of candidate fusion blanket materials

Three ferritic/martensitic alloys were studied to understand the synergistic effect between single ion beam (Fe 2+ ), dual ion beam (Fe 2+ +He 2+ and Fe 2+ +H + ), and triple ion beam (Fe 2+ +He 2+ +H + ) irradiations on cavity evolution. A commercial alloy, F82H, a castable nanostructured alloy, CNA3, and a high purity model alloy, Fe8Cr2W, were irradiated at 400°C to 600°C to a damage level of 50 dpa at a damage rate of 1 × 10 –3 dpa/s with He and H injection rates of 10 and 40 appm/dpa, respectively. Post-irradiation characterization via bright field transmission electron microscopy and high-angle annular dark-field scanning transmission electron microscopy was performed on all irradiated conditions to characterize the cavity size distribution and determine the effects of H/He injection on cavity microstructure. In all three alloys, hydrogen co-injection with helium resulted in an increased cavity number density and maximum cavity size, producing an increase in swelling over that from helium injection alone. Swelling in F82H appears to peak between 450°C and 500°C. At 600°C, swelling was minimal and cavities of high density and small size were confined to grain boundaries and dislocations while at 400°C, swelling is also low with a nearly homogeneous, high density, distribution of very small cavities throughout. Swelling was least in the commercial alloy F82H due to the high sink strength. The CNA3 alloy underwent dissolution of precipitates that lowered the sink strength and resulted in higher swelling than F82H, but less than the model alloy. Electron energy loss spectroscopy (EELS) elemental mapping revealed hydrogen forming a halo-like structure about the periphery of the cavities and helium residing within the cavities themselves. Finally, this observation suggests that hydrogen reduces the surface energy of helium-filled cavities which results in both increased cavity number density and cavity size in triple beam irradiation over dual beam irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Deuterium retention characteristics during lithium powder and granule injection in EAST

Fuel retention remains a critical challenge for magnetic confinement fusion devices. As a low-Z material, lithium plays a key role in tokamak wall conditioning and impurity control. In this study, the fuel retention behaviors associated with two lithium forms—powder and granule—are systematically examined using the gas balance method. Both forms of lithium injection significantly enhance fuel retention, shifting the wall behavior from net outgassing to net absorption. For lithium powder injection, compared with the reference discharge, a suppression efficiency exceeding 40% is achieved relative to the number of injected lithium atoms. Moreover, discharges with different injection rates show a monotonic increase in retained fuel with the injected lithium quantity. In contrast, repeated granule injections exhibit a pronounced cumulative effect, with the retention ratio varying from 0.16 to 0.65 over successive discharges, indicating progressive enhancement of wall absorption and deuterium retention. These results advance the understanding of wall behavior modification induced by solid material injection and provide insights for developing effective wall-conditioning strategies using low-Z materials in ITER and future fusion reactors, while the expected high retention may ultimately exclude the use of Li wall conditioning in future fusion devices.

Fuel retention↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
DOE OSTI2024

Carbon storage cost modeling for the offshore Gulf of Mexico

Groundbreaking for geologic carbon storage (GCS) projects in the offshore Gulf of Mexico is imminent, and there is great interest in utilizing this region for GCS projects. Offshore saline reservoirs provide a significant and accessible resource for GCS. However, conducting GCS in the offshore environment will pose distinct challenges pertaining to site selection, operations, infrastructure use, and monitoring compared to operating onshore that ultimately affect technoeconomic assessment of offshore GCS projects. Carbon storage and transport costs are critical to project developers looking to deploy carbon storage in the offshore environment. We present CO2_S_COM_Offshore, a model developed by the National Energy Technology Laboratory (NETL) as a screening-level offshore saline GCS cost modeling tool. Based on NETL’s widely used CO2_S_COM cost model for onshore saline CS, CO2_S_COM_Offshore enables technoeconomic analysis of GCS in offshore areas. This model comprehensively incorporates multiple facets of offshore GCS projects, from regional evaluation and site selection to permitting, transport, operations, monitoring, site closure, and decommissioning. In general, the model can explore the cost implications for potential offshore GCS project(s) by enabling the user to change several project operational and financial attribute configurations. Key inputs include offshore storage formation options, CO2 injection rate and duration, infrastructure types, monitoring intensity, project financing, and post-injection site care duration. Supporting cost algorithms within CO2_S_COM_Offshore were compiled utilizing S&P Global’ s QUE$TORTM cost estimation software alongside a variety of open-source scientific literature. In addition to reviewing key model components, we discuss several sensitivity analyses, input variabilities, and results on analysis of break-even CO2 price required by a project based on different regulation/policy and operational scenarios for the offshore Gulf of Mexico. These results indicate the value of modeling offshore GCS specifically, and the potential of offshore GCS within a decarbonization value chain. Presented at the 41st USAEE/IAEE North American Conference, 3-6 November 2024, Baton Rouge, LA, United States.

Mark-Moser, Mackenzie K.↗
DOE OSTI2025

Time-Lapse Electromagnetic Methods for Monitoring Plume Development in a Carbon Storage Reservoir

Conference presentation at International Meeting for Applied Geoscience & Energy (IMAGE), Houston, Texas, August 25–28, 2025. The Energy & Environmental Research Center (EERC) is leading applied research on electromagnetic (EM) monitoring methods at an active carbon storage site in North Dakota. Injection operations at the site began in February 2024, with a permitted injection rate of up to 2.7 million tonnes of CO 2 annually using six injection wells. CO 2 is captured on-site and injected into the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer located at a depth of approximately 1800 meters. The EERC led acquisition of multiple active- and passive-source EM techniques between August and October of 2024 to provide a thorough understanding of the resistivity profile at the site.

02 PETROLEUM↗
DOE OSTI2023

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗
DOE OSTI2021

Preferential Direction of Electron Transfers at a Dye–Metal Oxide Interface with an Insulating Fluorinated Self-Assembled Monolayer and MgO

Slowing non-desirable electron transfer reactions at metal oxide-dye interfaces is important for many technologies. Specifically, after an interfacial photoinduced charge separation event at a metal oxide-dye interface, it is critically important to limit the rate of electron transfer reactions back to sensitizers and to limit electron transfer reactions between the electrolyte and the metal oxide. Ruthenium-based dyes at metal oxide interfaces are widely used in many fields; however, these dyes often have poor surface insulation resulting in fast recombination kinetics with transition metal-based redox shuttles in an electrolyte. This work explores two semiconductor surface modification strategies designed to minimize recombination events of electrons in TiO 2 with oxidized redox shuttles using a fluorinated siloxane insulator (PFTS) and a metal oxide insulator (MgO) with a well-known Ru dye, B11. Additionally, the influence of these treatments on the rate and duration of photoinduced interfacial charge separation at the TiO 2 -dye interface was examined. The TiO 2 -dye-redox shuttle systems were studied via dye-sensitized solar cell current-voltage curve, incident photon-to-current conversion efficiency, small-modulated photovoltage transient, time-correlated single photon counting, and transient absorption spectroscopy measurements. MgO was found to decrease the rate of the electron transfer reaction from the metal oxide to the electrolyte, decrease the rate of the electron transfer reaction to the oxidized dye from TiO 2 , increase the electron transfer reaction rate from a reduced redox shuttle to an oxidized dye, and decrease the electron injection rate from the photoexcited dye to TiO2. Interestingly, PFTS was found to desirably improve on these rates relative to MgO or untreated TiO 2 . Here, a model based on electrostatic interactions is presented to explain the exceptional behavior of PFTS with DFT computational analysis of PFTS supporting this model.

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

ELM Suppression by Boron Powder Injection and Comparison with Lithium Powder Injection on EAST

Type I edge-localized modes (ELMs) in the Experimental Advanced Superconducting Tokamak (EAST) were completely suppressed via boron powder injection into the X-point region of an upper-single null configuration over a wide range of operating conditions (2.8 < P aux < 7.5 MW, 3.8 × 10 19 < n e < 6 × 10 19 m -3 , RF-only and RF + NBI heating scenarios, both grad-B drift directions, and even He ion majority plasmas) (Sun et al. in Nucl. Fusion, 2020). A window of edge B concentration for stable long pulse operation was identified: too low and ELMs return, too high and the discharge suffers radiative collapse. The injection of boron powder above the minimum for ELM suppression coincided with the occurrence of an edge harmonic oscillation detected in magnetics (both on the high-field side and low-field side), in AXUV diodes near the upper X-point, divertor Dα emission, and in a range of other diagnostics (Diallo et al. in: Proceedings of 2020 IAEA fusion energy conference, 2021). Furthermore, no harmonic oscillation was observed when ELMs were present, and stored energy was slightly increased at constant density during ELM suppression. Core tungsten emission during ELM suppression either increased or decreased relative to ELMy H-mode, but the W emission was maintained at acceptable levels. The threshold B injection rate was measured for several conditions, and found to increase with heating power. Li powder injection into comparable discharges also resulted in a short phase of ELM suppression, but density and stored energy both decreased due to the strong pumping effect of lithium; no edge harmonic oscillation was observed with Li injection, indicating that the ELM suppression mechanisms differ. The new set of B-seeded, ELM-suppressed discharges exhibited certain characteristics of quiescent H-mode (Burrell et al. in Phys Plasmas 8:2153, 2001), but did not require high shear, counter beams, etc. The wide operating window and compatibility with RF-only discharges paves the way for future experiments targeting long pulse H-mode discharges with complete ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Scaling of Turbulence and Microphysics in a Convection–Cloud Chamber of Varying Height

Abstract The convection–cloud chamber enables measurement of aerosol and cloud microphysics, as well as their interactions, within a turbulent environment under steady‐state conditions. Increasing the size of a convection–cloud chamber, while holding the imposed temperature difference constant, leads to increased Rayleigh, Reynolds and Nusselt numbers. Large–eddy simulation coupled with a bin microphysics model allows the influence of increased velocity, time, and spatial scales on cloud microphysical properties to be explored. Simulations of a convection–cloud chamber, with fixed aspect ratio and increasing heights of H = 1, 2, 4, and (for dry conditions only) 8 m are performed. The key findings are: Velocity fluctuations scale as H 1/3 , consistent with the Deardorff expression for convective velocity, and implying that the turbulence correlation time scales as H 2/3 . Temperature and other scalar fluctuations scale as H −3/7 . Droplet size distributions from chambers of different sizes can be matched by adjusting the total aerosol injection rate as the horizontal cross‐sectional area (i.e., as H 2 for constant aspect ratio). Injection of aerosols at a point versus distributed throughout the volume makes no difference for polluted conditions, but can lead to cloud droplet size distribution broadening in clean conditions. Cloud droplet growth by collision and coalescence leads to a broader right tail of the distribution compared to condensation growth alone, and this tail increases in magnitude and extent monotonically as the increase of chamber height. These results also have implications for scaling within turbulent, cloudy mixed‐layers in the atmosphere, such as fog layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Facile Synthesis of Palladium-Based Nanocrystals with Different Crystal Phases and a Comparison of Their Catalytic Properties

A relatively unexplored aspect of noble-metal nanomaterials is polymorphism, or their ability to crystallize in different crystal phases. Here, in this work, a method is reported for the facile synthesis of Ru@Pd core–shell nanocrystals featuring polymorphism, with the core made of hexagonally close-packed (hcp)-Ru while the Pd shell takes either an hcp or face-centered cubic (fcc) phase. The polymorphism shows a dependence on the shell thickness, with shells thinner than ≈1.4 nm taking the hcp phase whereas the thicker ones revert to fcc. The injection rate provides an experimental knob for controlling the phase, with one-shot and drop-wise injection of the Pd precursor corresponding to fcc-Pd and hcp-Pd shells, respectively. When these nanocrystals are tested as catalysts toward formic acid oxidation, the Ru@Pd hcp nanocrystals outperform Ru@Pd fcc in terms of both specific activity and peak potential. Density functional theory calculations are also performed to elucidate the origin of this performance enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Humate Amendment Injection Viability Testing – Supplementary Batch Mixture and Soil Column Studies

Groundwater in the Lost Lake Aquifer Zone (LLAZ) in the Southern Sector of the M-Area Hazardous Waste management Facility (HWMF) is contaminated with chlorinated ethenes, including trichloroethylene (TCE) and tetrachloroethylene (PCE). Treatment of contaminated groundwater with humic acid is being evaluated as a potential corrective action for these volatile organic compounds (VOCs) in the LLAZ in Southern Sector (SRNS, 2019). Pilot scale injections of Huma-K brand humic acid for groundwater treatment were previously performed in M-Area between 2017 and 2020 (Amidon, 2023). Groundwater was extracted from the Lost Lake Aquifer Zone (LLAZ), then a solution of humate was mixed with recovered groundwater and subsequently re-injected into the same aquifer unit. Operations of the historical humate pilot testing were discontinued in 2020 due to COVID-19 restrictions, low injection rates, and recirculation fluid spillage. No further field scale or laboratory research to investigate humate amended groundwater injection was performed at the SRS until 2024. At the request of Area Completion Projects (ACP), the Savannah River National Laboratory (SRNL) performed supplementary testing to help determine the viability of further injection of humate as a remedial option for VOCs within the LLAZ in the Southern Sector of the M-Area HWMF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Simple Statistical Models for Predicting Overpressure Due to CO2 and Low-Salinity Waste-Fluid Injection into Deep Saline Formations

Deep saline aquifers have been used for waste-fluid disposal for decades and are the proposed targets for large-scale CO2 storage to mitigate CO2 concentration in the atmosphere. Due to relatively limited experience with CO2 injection in deep saline formations and given that the injection targets for CO2 sometimes are the same as waste-fluid disposal formations, it could be beneficial to model and compare both practices and learn from the waste-fluid disposal industry. In this paper, we model CO2 injection in the Patterson Field, which has been proposed as a site for storage of 50 Mt of industrial CO2 over 25 years. We propose general models that quickly screen the reservoir properties and calculate pressure changes near and far from the injection wellbore, accounting for variable reservoir properties. The reservoir properties we investigated were rock compressibility, injection rate, vertical-to-horizontal permeability ratio, average reservoir permeability and porosity, reservoir temperature and pressure, and the injectant total dissolved solids (TDS) in cases of waste-fluid injection. We used experimental design to select and perform simulation runs, performed a sensitivity analysis to identify the important variables on pressure build-up, and then fit a regression model to the simulation runs to obtain simple proxy models for changes in average reservoir pressure and bottomhole pressure. The CO2 injection created more pressure compared to saline waste-fluids, when similar mass was injected. However, we found a more significant pressure buildup at the caprock-reservoir interface and lower pressure buildup at the bottom of the reservoir when injecting CO2 compared with waste-fluid injection.

Ansari, Esmail↗
DOE OSTI2020

Reply to Comment by Peterie Et Al. on “Accelerated Fill‐Up of the Arbuckle Group Aquifer and Links to U.S. Midcontinent Seismicity”

Abstract Peterie et al. question one observation in our paper: associating pressure increases to injection volumes at distances of up to 25 km from an injection well. In this reply, we show that the comment misunderstands our analysis and the evidence that led to this conclusion. We also show that gauge‐depth‐corrected pressures, used by the authors to produce statewide pressure maps, are discrepant with the static fluid level data, provided in our original compilation and analysis. The discrepancies are a result of the pressure correction method employed, which naïvely substitutes formation pressure for bottomhole pressure to calculate wellbore fluid density. Their linearly interpolated pressure maps, based on sparse data, contain interpolation and extrapolation artifacts that contradict injection trends in the state, the Theis solution, and the superposition principle. We reiterate that pressure and static fluid level increases in Class I wells existed prior to 2013, most notably in central Kansas, where recent earthquakes are cited in the comment as evidence of a pressure plume emanating from the Kansas‐Oklahoma border, 90 km away. We show that the space‐time pattern of seismicity in this area is inconsistent with a northward propagating pressure plume and, instead, seismicity appears to be centered on and near a cluster of high‐rate injection wells, two of which are among the highest rate wells in the state. These observations, along with recent M4 + earthquakes during continued decreases in wastewater injection in southern Kansas and northern Oklahoma, question the usefulness of the comment for understanding and managing societally significant earthquakes.

Ansari, Esmail↗