Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid injection”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Reactive Transport Modeling of Shale–Fluid Interactions after Imbibition of Fracturing Fluids

Injection of hydraulic fracturing fluid (HFF) into shale formations for unconventional oil/gas production results in chemical reactions in the shale matrix. Our recent experimental study determined the depths to which different types of reactions between the shale matrix and the HFF extended. In the present study, we built continuum-scale reactive transport models to understand the coupling of chemical reactions and the transport of aqueous species in these shale–HFF systems. Calibration of the model with our previous experimental results reveals that it takes hours to months for the shale matrix to completely neutralize the imbibed acids, depending primarily upon the carbonate content of the shale. Both the HFF pH and pore pH affect the location of barite precipitation, resulting in unique barite precipitation profiles extending millimeters into calcite-rich Eagle Ford shale but only tens of micrometers into low-carbonate Marcellus shale. In addition, dissolved oxygen and extracted bitumen are key to reproducing the experimental observation of Fe(III) (oxyhydr)oxide formation in the shale matrix as a result of pyrite dissolution in the shales. A comparison between the modeling results of porosity in the present study to experimentally measured permeabilities in our previous study suggests that chemical reactions occurring at a greater depth than the observable reaction zone might have impacted permeability. Our model serves as a benchmark for efficiently modeling water–rock interactions in similar systems where bulk rock samples react with a solution in batch reactors. Important reactive transport processes were ascertained via modeling, which allows for quantitative prediction of shale–HFF interactions in shale matrices given the shale and HFF compositions.

58 GEOSCIENCES↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗
DOE OSTI2024

Modeling injection-induced fault slip using long short-term memory networks

Stress changes due to changes in fluid pressure and temperature in a faulted formation may lead to the opening/shearing of the fault. This can be due to subsurface (geo)engineering activities such as fluid injections and geologic disposal of nuclear waste. Such activities are expected to rise in the future making it necessary to assess their short- and long-term safety. Here, a new machine learning (ML) approach to model pore pressure and fault displacements in response to high-pressure fluid injection cycles is developed. The focus is on fault behavior near the injection borehole. To capture the temporal dependencies in the data, long short-term memory (LSTM) networks are utilized. To prevent error accumulation within the forecast window, four critical measures to train a robust LSTM model for predicting fault response are highlighted: (i) setting an appropriate value of LSTM lag, (ii) calibrating the LSTM cell dimension, (iii) learning rate reduction during weight optimization, and (iv) not adopting an independent injection cycle as a validation set. Several numerical experiments were conducted, which demonstrated that the ML model can capture peaks in pressure and associated fault displacement that accompany an increase in fluid injection. The model also captured the decay in pressure and displacement during the injection shut-in period. Further, the ability of an ML model to highlight key changes in fault hydromechanical activation processes was investigated, which shows that ML can be used to monitor risk of fault activation and leakage during high pressure fluid injections.

58 GEOSCIENCES↗
DOE OSTI2021

Multi-fluid renewable geo-energy systems and methods

A geo-energy production system and method extracts thermal energy from a reservoir formation, and stores either thermal waste heat or excess heat in a storage zone of the reservoir formation. A compressed fluid injection injects an unheated, compressed working fluid into the storage zone. A fluid injection well injects a working fluid laden with thermal waste heat or excess heat into the storage zone. The storage zone is located below a caprock layer and above a native brine zone of the reservoir formation and is partially circumscribed by a hot brine storage zone. The compressed working fluid assists with a withdrawal of pressurized brine residing below and/or to the sides of the storage zone. A compressed CO2, N2, or air production well helps to remove compressed working fluid from the storage zone for use in power production.

Buscheck, Thomas A.↗
DOE OSTI2022

Operational Forecasting of Induced Seismicity

Operational Forecasting of Induced Seismicity (ORION) is an open-source toolkit developed as part of the US Department of Energy SMART and NRAP programs. ORION is designed to help both expert and non-expert users understand the seismic hazard at a site. This tool analyzes how the seismic hazard evolves during injection of fluid, and can be used to identify potential mitigation strategies for observed or forecasted seismic activity. ORION is built using Python, and has an optional GUI where users can easily configure model parameters, evaluate temporal and spatial forecasts, and inspect key data (seismicity, fluid injection, pressure models, etc.).

Kroll, Kayla↗
DOE OSTI2026

Role of Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes for Enhanced Geothermal Systems

The ability to sustain high energy extraction efficiency from Enhanced Geothermal Systems (EGS) is affected by the ability to measure, characterize, and predict the effect of different perturbations that arise from fluid injection rates, fluid temperature, and shut-in conditions on existing and stimulated fracture systems throughout a subsurface reservoir’s lifecycle. The difficulty arises from limited knowledge of the interaction of fluid and temperature driven fractures with frictional interfaces that govern local deformation and frictional behavior in the surrounding rock. Thermo-poro-mechanical coupled processes tend to dominate this interaction and strongly affect the permeability and local stress distributions that impact efficiency and prevent optimal stimulation conditions over multi-decadal time frames.

Pyrak-Nolte, Laura [Purdue Univ., West Lafayette, ↗
DOE OSTI2020

More Than 40 yr of Potentially Induced Seismicity Close to the San Andreas Fault in San Ardo, Central California

Evidence for fluid–injection–induced seismicity is rare in California hydrocarbon basins, despite widespread injection close to seismically active faults. We investigate a potential case of injection–induced earthquakes associated with San Ardo oilfield operations that began in the early 1950s. The largest potentially induced events occurred in 1955 (M L 5.2) and 1985 (M w 4.5) within ~6 km from the oilfield. We analyze Synthetic Aperture Radar interferometric images acquired by Sentinel–1A/B satellites between 2016 and 2020 and find surface deformation of up to 1.5 cm/yr, indicating pressure–imbalance in parts of the oilfield. Fluid injection in San Ardo is concentrated within highly permeable rocks directly above the granitic basement at a depth of ~800 m. Seismicity predominantly occurs along basement faults at 6–13 km depths. Seismicity and wastewater disposal wells are spatially correlated to the north of the oilfield. Temporal correlations are observed over more than 40 yr with correlation coefficients of up to 0.71 for seismicity within a 24 km distance from the oilfield. Such large distances have not previously been observed in California but are similar to the large spatial footprint of injection in Oklahoma. Furthermore, the San Ardo seismicity shows anomalous clustering with earthquakes consistently occurring at close spatial proximity but long interevent times. Similar clustering has previously been reported in California geothermal fields and may be indicative of seismicity driven by long–term, spatially persistent external forcing.

58 GEOSCIENCES↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗
DOE OSTI2023

Geomechanical in situ testing of fault reactivation in argillite repositories

Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.

Birkholzer, Jens T.↗
DOE OSTI2023

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗
DOE OSTI2021

Geochemical reactions and alteration of pore architecture in saturated shale after injection of stimulation fluid

Pore architecture regulates fluid flow between unconventional shale reservoir and hydraulically-induced fractures. Imbibition of stimulation fluid may change this architecture and alter hydrocarbon flow. Hydrothermal experiments were conducted at reservoir conditions (125 °C, 45 MPa) to test two hypotheses: 1) Shale, not stimulation fluid, dominates the geochemistry of an unconventional reservoir containing formation water; and 2) Mineral dissolution and precipitation induced by stimulation fluid is transient and manifests across micro-, meso-, and macro-scales of pore architecture. Experiments replicated a shut-in well in the Cretaceous Baxter Shale, Green River Basin, Wyoming USA. Stimulation fluid was injected into one experiment after formation water (I = 0.35 mol/kg, pH = 6.4) and core reacted for 48 days. This novel approach equilibrated rock and formation water and saturated pores with formation water before introducing stimulation fluid. The second experiment served as a control. Trends of aqueous calcium, silica and aluminum in the injection experiment suggest transient dissolution of calcite and feldspar and/or quartz as well as clay and barite precipitation; mineralogic evidence was limited to calcite dissolution and barite precipitation. The rock maintained reducing conditions (Eh = +0.08 to -0.16 V) despite injection of oxidizing stimulation fluid (Eh = 1.1 V). pH of the stimulation fluid-formation water mixture evolved from 2.3 to pre-injection values (~6) within 24 h. The results indicate that mineral dissolution and precipitation manifest in macropores with barely detectable alteration to micro- and mesopores. Formation water or organic matter in pores may have inhibited access of stimulation fluid to micro- and mesopores

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

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗
DOE OSTI2023

Microseismic Monitoring of Laboratory Hydraulic Fracturing Experiments in Granitic Rocks for Different Fracture Propagation Regimes

While hydraulic fracturing (HF) is a widely employed process, the underlying fracturing processes are still heavily contested. The attributes of the HF generated fracture network can exhibit substantial variation when dealing with specific HF propagation regimes encountered in the field. In this study, HF experiments were performed on true-triaxially loaded Barre granite cubes, with microseismic monitoring, to identify and characterize the fracturing mechanisms associated with different viscosity injection fluids. Utilizing fluids with high (oil/1450 cP) and low (water/1 cP) viscosity represented two key HF propagation regimes: viscosity- and toughness-dominated. The experiments conducted with oil involved higher breakdown pressures, larger fluid volumes, and slower fracture propagation speeds. Furthermore, the frequency–magnitude distribution (b value) for all experiments (1.9–2.3) was similar to those encountered for large-scale operations. Slightly larger b values were encountered during the initiation phase (2.4–2.7) relative to the fracture propagation and post-fracturing phases (1.9–2.2). Techniques such as polarity and moment-tensor inversion were utilized to characterize the source mechanisms. For the HF experiments with oil, tensile fractures were most dominant (92%) in the initiation phase compared to fracture propagation and post-fracturing phases (70–75%). Similar tensile fracturing dominancy was not observed with water, attributable to fluid permeation and leak-off. Regardless of the injection fluid or classification criteria employed, tensile fractures were the dominant type consistently, with fewer occurring in water experiments but the specific ratio of crack types varied with different source mechanism criteria employed.

58 GEOSCIENCES↗
DOE OSTI2021

Injection Parameters That Promote Branching of Hydraulic Cracks

Fluid injection into rock formations can either produce complex branched hydraulic fractures, create simple planar fractures, or be dominated by porous diffusion. Currently, the optimum injection parameters to create branched fractures are unknown. We conducted repeatable hydraulic fracturing experiments using analog-rock samples with controlled heterogeneity to quantify the fluid parameters that promote fracture branching. A large range of injection rates and fluid viscosities were used to investigate their effects on induced fracture patterns. Paired with a simple analytical model, our results identify the threshold at which fracture transitions from an isolated planar crack to branched cracks when closed natural fractures exist. These results demonstrate that this transition can be controlled by injection rate and fluid viscosity. In relation to the field practices, the present model predicts slickwater and lower viscosity fluid injections promote fracture branching, with the Marcellus shale used as an example.

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 OSTI2022

The 2017 Pohang, South Korea, Mw 5.4 main shock was either natural or triggered, but not induced

Understanding the cause of the November 2017 Pohang main shock is of considerable importance to the future of the geothermal industry because of its large magnitude compared to prior expectations based on case histories of other projects involving underground fluid injection. Of the three possibilities – induced, triggered or natural, “induced” can be ruled out based on the disproportionately large seismic moment of the main shock. Furthermore, whether natural or triggered, the source of seismic hazard at Pohang was tectonic strain accumulation, not fluid injection. Arguably, the most timely indicator of seismic hazard and risk in the environs of Pohang was the September 2016 Mw 5.4 Gyeongju earthquake, which was natural and located about 40 km south of Pohang along the same active fault system.

58 GEOSCIENCES↗