There was an evident increase in the number of earthquakes in the Xinfengjiang Reservoir from June to July 2014 after the landing of Typhoon Hagibis.To understand the spatial and temporal evolution of this microseismi...There was an evident increase in the number of earthquakes in the Xinfengjiang Reservoir from June to July 2014 after the landing of Typhoon Hagibis.To understand the spatial and temporal evolution of this microseismicity,we built a high-precision earthquake catalog for 2014 and relocated 2275 events using recently developed methods for event picking and catalog construction.Seismicity occurred in the southeastern part of the reservoir,with the preferred fault plane orientation aligned along the Heyuan Fault.The total seismic energy peaked when the typhoon passed through the reservoir,and seismicity correlated with typhoon energy.In contrast,a limited seismic response was observed during the later Typhoon Rammasun.Combining data regarding the water level in the Xinfengjiang Reservoir and seismicity frequency changes in the Taiwan region during these two typhoon events,we suggest that typhoon activity may increase microseism energy by impacting fault stability around the Xinfengjiang Reservoir.Whether a fault can be activated also depends on how close the stress accumulation is to its failure point.展开更多
Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operation...Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.展开更多
Monitoring seismicity in real time provides significant benefits for timely earthquake warning and analyses.In this study,we propose an automatic workflow based on machine learning(ML)to monitor seismicity in the sout...Monitoring seismicity in real time provides significant benefits for timely earthquake warning and analyses.In this study,we propose an automatic workflow based on machine learning(ML)to monitor seismicity in the southern Sichuan Basin of China.This workflow includes coherent event detection,phase picking,and earthquake location using three-component data from a seismic network.By combining Phase Net,we develop an ML-based earthquake location model called Phase Loc,to conduct real-time monitoring of the local seismicity.The approach allows us to use synthetic samples covering the entire study area to train Phase Loc,addressing the problems of insufficient data samples,imbalanced data distribution,and unreliable labels when training with observed data.We apply the trained model to observed data recorded in the southern Sichuan Basin,China,between September 2018 and March 2019.The results show that the average differences in latitude,longitude,and depth are 5.7 km,6.1 km,and 2 km,respectively,compared to the reference catalog.Phase Loc combines all available phase information to make fast and reliable predictions,even if only a few phases are detected and picked.The proposed workflow may help real-time seismic monitoring in other regions as well.展开更多
Earthquakes triggered by dynamic disturbances have been confirmed by numerous observations and experiments.In the past several decades,earthquake triggering has attracted increasing attention of scholars in relation t...Earthquakes triggered by dynamic disturbances have been confirmed by numerous observations and experiments.In the past several decades,earthquake triggering has attracted increasing attention of scholars in relation to exploring the mechanism of earthquake triggering,earthquake prediction,and the desire to use the mechanism of earthquake triggering to reduce,prevent,or trigger earthquakes.Natural earthquakes and large‐scale explosions are the most common sources of dynamic disturbances that trigger earthquakes.In the past several decades,some models have been developed,including static,dynamic,quasi‐static,and other models.Some reviews have been published,but explosiontriggered seismicity was not included.In recent years,some new results on earthquake triggering have emerged.Therefore,this paper presents a new review to reflect the new results and include the content of explosion‐triggered earthquakes for the reference of scholars in this area.Instead of a complete review of the relevant literature,this paper primarily focuses on the main aspects of dynamic earthquake triggering on a tectonic scale and makes some suggestions on issues that need to be resolved in this area in the future.展开更多
Based on the seismic data recorded by the China Earthquake Networks Center(CENC) in the Luxian area from January 2009 to October 2021,the 3D V_P,V_S, V_P/V_S structures and seismic locations of the area are obtained b...Based on the seismic data recorded by the China Earthquake Networks Center(CENC) in the Luxian area from January 2009 to October 2021,the 3D V_P,V_S, V_P/V_S structures and seismic locations of the area are obtained by joint inversion using the V_P/V_S model consistency-constrained double-difference tomography method(tomoDDMC).The earthquakes in the study area are mainly concentrated at a depth of 2-6 km,and the focal depth is generally shallow.The Ms 6.0 Luxian earthquake occurred at the transition zone of high-and low-velocity anomalies and the aftershock sequence was distributed along the edge of the low-V_P zone.A small number of foreshocks occurred on the west side of the M_S 6.0 Luxian earthquake,while most of the aftershocks were distributed on the east side of the M_S 6.0 Luxian earthquake.The aftershock sequence consisted of three seismic bands with different trends,and the overall distribution was in a NWW direction,which was inconsistent with the spatial distribution of the main active faults nearby.In addition,the spatiotemporal distribution of earthquakes and the variation of b-values are closely related to the industrial water injection activities in the study area,reflecting the activation of pre-existing hidden faults under certain tectonic and stress environments leading to seismic activities in the area.展开更多
Buckling-restrained braces(BRBs)have shown their capability to provide building structures with stiffness,strength,and ductility.Estimating the seismic drifts of buckling-restrained braced frames(BRBFs)is an important...Buckling-restrained braces(BRBs)have shown their capability to provide building structures with stiffness,strength,and ductility.Estimating the seismic drifts of buckling-restrained braced frames(BRBFs)is an important design step to control structural and non-structural damage.In current practice of seismic design,the estimation of seismic drifts of BRBFs is performed by using empirical calculations that are independent upon either the type of the structural system or the design level of seismicity.In these empirical calculations,the seismic drifts are estimated by amplifying the reduced elastic drifts obtained under design lateral loading with a displacement amplification factor(DAF).The value of DAF is considered equal to the product of the response modification factor R and the inelastic displacement ratioρ.The goal of the current research is to assess the value ofρfor low-to mid-rise BRBFs designed under low and high levels of seismicity.This goal has been achieved by conducting a series of elastic and inelastic time-history analyses pertaining to an ensemble of earthquake records on 3-,6-and 9-story BRBFs.The results indicate that theρ-ratio increases with an increase in design seismic intensity and an increase in experienced inelasticity.The range ofρfor low seismicity designs ranges from 0.63 to 0.9,while for high seismicity designs this range stretches from 0.83 to 1.29.It has been found that the consideration of a generalρ-ratio of 1.0 is a reasonable estimation for the design of the BRBFs considered in this study.展开更多
A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distri...A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distribution and magnitude of hydraulic fracturing-induced seismicity. An integrated machine learning-based investigation is conducted to systematically evaluate multiple factors that contribute to induced seismicity. Feature importance indicates that a distance to fault, a distance to basement, minimum principal stress, cumulative fluid injection, initial formation pressure, and the number of fracturing stages are among significant model predictors. Our seismicity prediction map matches the observed spatial seismicity, and the prediction model successfully guides the fracturing job size of a new well to reduce seismicity risks. This study can apply to mitigating potential seismicity risks in other seismicity-frequent regions.展开更多
We present a detailed catalog of 13671 earthquakes in the Eastern Tennessee Seismic Zone(ETSZ)that spans January 1,2005 to July 31,2020.We apply a matched filter detection technique on over 15 years of continuous data...We present a detailed catalog of 13671 earthquakes in the Eastern Tennessee Seismic Zone(ETSZ)that spans January 1,2005 to July 31,2020.We apply a matched filter detection technique on over 15 years of continuous data,resulting in arguably the most complete catalog of seismicity in the ETSZ yet.The magnitudes of newly detected events are determined by computing the amplitude ratio between the detections and templates using a principal component fit.We also compute the b-value for the new catalog and comparatively relocate a subset of newly detected events using XCORLOC and hypoDD,which shows a more defined structure at depth.We find the greatest concentration along and to the east of the New York-Alabama Lineament,as defined by the magnetic anomaly,supporting the argument that this feature likely is related to the generation of seismicity in the ETSZ.We examine seismicity in the vicinity of the Watts Bar Reservoir,which is located about 5 km from the epicenter of the M_(W) 4.4 December 12,2018 Decatur,Tennessee earthquake,and find possible evidence for reservoir modulated seismicity in this region.We also examine seismicity in the entire ETSZ to search for a correlation between shallow earthquakes and seasonal hydrologic changes.Our results show limited evidence for hydrologicallydriven shallow seismicity due to seasonal groundwater levels in the ETSZ,which contradicts previous studies hypothesizing that most intraplate earthquakes are associated with the dynamics of hydrologic cycles.展开更多
Induced seismicity is strongly related to various engineering projects that cause anthropogenic in-situ stress change at a great depth.Hence,there is a need to estimate and mitigate the associated risks.In the past,va...Induced seismicity is strongly related to various engineering projects that cause anthropogenic in-situ stress change at a great depth.Hence,there is a need to estimate and mitigate the associated risks.In the past,various simulation methods have been developed and applied to induced seismicity analysis,but there is still a fundamental diference between simulation results and feld observations in terms of the spatial distribution of seismic events and its frequency.The present study aims to develop a method to simulate spatially distributed on-fault seismicity whilst reproducing a complex stress state in the fault zone.Hence,an equivalent continuum model is constructed,based on a discrete fracture network within a fault damage zone,by employing the crack tensor theory.A fault core is simulated at the center of the model as a discontinuous plane.Using the model,a heterogeneous stress state with stress anomalies in the fault zone is frst simulated by applying tractions on the model outer boundaries.Subsequently,the efective normal stress on the fault plane is decreased in a stepwise manner to induce slip.The simulation result is validated in terms of the b-value and other seismic source parameters,hence demonstrating that the model can reproduce spatially and temporally distributed on-fault seismicity.Further analysis on the parameters shows the variation of frequency-magnitude distribution before the occurrence of large seismic events.This variation is found to be consistent with feld observations,thus suggesting the potential use of this simulation method in evaluating the risk for seismic hazards in various engineering projects.展开更多
Post shut‐in seismic events in enhanced geothermal systems(EGSs)occur predominantly at the outer rim of the co‐injection seismic cloud.The concept of postinjection fracture and fault closure near the injection well ...Post shut‐in seismic events in enhanced geothermal systems(EGSs)occur predominantly at the outer rim of the co‐injection seismic cloud.The concept of postinjection fracture and fault closure near the injection well has been proposed and validated as a mechanism for enhancing post shut‐in pressure diffusion that promotes seismic hazard.This phenomenon is primarily attributed to the poro‐elastic closure of fractures resulting from the reduction of wellbore pressure after injection termination.However,the thermal effects in EGSs,mainly including heat transfer and thermal stress,may not be trivial and their role in postinjection fault closure and pressure evolution needs to be explored.In this study,we performed numerical simulations to analyze the relative importance of poro‐elasticity,heat transfer,and thermo‐elasticity in promoting postinjection fault closure and pressure diffusion.The numerical model wasfirst validated against analytical solutions in terms offluid pressure diffusion and against heatedflow‐through experiments in terms of thermal processes.We then quantified and distinguished the contribution of each individual mechanism by comparing four different shut‐in scenarios simulated under different coupled conditions.Our results highlight the importance of poro‐elastic fault closure in promoting postinjection pressure buildup and seismicity,and suggest that heat transfer can further augment the fault closure‐induced pressure increase and thus potentially intensify the postinjection seismic hazard,with minimal contribution from thermo‐elasticity.展开更多
The paper is devoted to analysis of hydrogeological, geomagnetic and seismic response to the two great remote geophysical events, 2022 Tonga volcano eruption and 2020-2023 Türkiye earthquakes in Georgia (Caucasus...The paper is devoted to analysis of hydrogeological, geomagnetic and seismic response to the two great remote geophysical events, 2022 Tonga volcano eruption and 2020-2023 Türkiye earthquakes in Georgia (Caucasus). The geophysical observation system in Georgia, namely, water level stations in the network of deep wells, atmospheric pressure and the geomagnetic sensors of the Dusheti Geophysical Observatory (DGO) as well as seismic data in Garni Observatory (Armenia) respond to the Tonga event by anomalies in the time series. These data show that there are two types of respond: infrasound disturbances in atmospheric pressure and seismic waves in the Earth generated by the eruption. After Tonga eruption January 15 at 04:21 UTC three groups of N-shaped waveforms were registered in the water level corresponding to the global propagation characteristics of the N-shaped waveform of infrasound signals on the barograms generated by eruption at the distance ~15,700 km: they were identified as the Lamb wave, a surface wave package running in the atmosphere with a velocity around ~314 m/s. The paper also presents the WL reactions to three strong EQs that occur in Türkiye 2020-2023, namely Elazığ, Van and Türkiye-Syria EQs. WL in Georgian well network reacts to these events by anomalies of different intensity, which points to the high sensitivity of hydrosphere to remote (several hundred km) strong EQs. The intensity and character of WL reactions depend strongly on the local hydrogeological properties of rocks, surrounding the well.展开更多
基金Strategic Priority Research Program(B)of the Chinese Academy of Sciences(No.XDB42020304)National Natural Science Foundation of China(No.42074059).
文摘There was an evident increase in the number of earthquakes in the Xinfengjiang Reservoir from June to July 2014 after the landing of Typhoon Hagibis.To understand the spatial and temporal evolution of this microseismicity,we built a high-precision earthquake catalog for 2014 and relocated 2275 events using recently developed methods for event picking and catalog construction.Seismicity occurred in the southeastern part of the reservoir,with the preferred fault plane orientation aligned along the Heyuan Fault.The total seismic energy peaked when the typhoon passed through the reservoir,and seismicity correlated with typhoon energy.In contrast,a limited seismic response was observed during the later Typhoon Rammasun.Combining data regarding the water level in the Xinfengjiang Reservoir and seismicity frequency changes in the Taiwan region during these two typhoon events,we suggest that typhoon activity may increase microseism energy by impacting fault stability around the Xinfengjiang Reservoir.Whether a fault can be activated also depends on how close the stress accumulation is to its failure point.
基金funded by the joint fund of the National Key Research and Development Program of China(No.2021YFC2902101)National Natural Science Foundation of China(Grant No.52374084)+1 种基金Open Foundation of National Energy shale gas R&D(experiment) center(2022-KFKT-12)the 111 Project(B17009)。
文摘Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.
基金the financial support of the National Key R&D Program of China(2021YFC3000701)the China Seismic Experimental Site in Sichuan-Yunnan(CSES-SY)。
文摘Monitoring seismicity in real time provides significant benefits for timely earthquake warning and analyses.In this study,we propose an automatic workflow based on machine learning(ML)to monitor seismicity in the southern Sichuan Basin of China.This workflow includes coherent event detection,phase picking,and earthquake location using three-component data from a seismic network.By combining Phase Net,we develop an ML-based earthquake location model called Phase Loc,to conduct real-time monitoring of the local seismicity.The approach allows us to use synthetic samples covering the entire study area to train Phase Loc,addressing the problems of insufficient data samples,imbalanced data distribution,and unreliable labels when training with observed data.We apply the trained model to observed data recorded in the southern Sichuan Basin,China,between September 2018 and March 2019.The results show that the average differences in latitude,longitude,and depth are 5.7 km,6.1 km,and 2 km,respectively,compared to the reference catalog.Phase Loc combines all available phase information to make fast and reliable predictions,even if only a few phases are detected and picked.The proposed workflow may help real-time seismic monitoring in other regions as well.
基金supported by the National Natural Science Foundation of China(NSFC grants No.12172036,51774018)the Program for Changjiang Scholars and Innovative Research Team in University(PCSIRT,IRT_17R06)+2 种基金the Russian Foundation for Basic Research,Grant Number 20‐55‐53032Russian State Task number 1021052706247‐7‐1.5.4the Government of Perm Krai,research project No.С‐26/628.
文摘Earthquakes triggered by dynamic disturbances have been confirmed by numerous observations and experiments.In the past several decades,earthquake triggering has attracted increasing attention of scholars in relation to exploring the mechanism of earthquake triggering,earthquake prediction,and the desire to use the mechanism of earthquake triggering to reduce,prevent,or trigger earthquakes.Natural earthquakes and large‐scale explosions are the most common sources of dynamic disturbances that trigger earthquakes.In the past several decades,some models have been developed,including static,dynamic,quasi‐static,and other models.Some reviews have been published,but explosiontriggered seismicity was not included.In recent years,some new results on earthquake triggering have emerged.Therefore,this paper presents a new review to reflect the new results and include the content of explosion‐triggered earthquakes for the reference of scholars in this area.Instead of a complete review of the relevant literature,this paper primarily focuses on the main aspects of dynamic earthquake triggering on a tectonic scale and makes some suggestions on issues that need to be resolved in this area in the future.
基金supported by the Fundamental Research Funds for the Institute of Earthquake Forecasting, China Earthquake Administration(CEAIEF20220211 and CEAIEF20220401)。
文摘Based on the seismic data recorded by the China Earthquake Networks Center(CENC) in the Luxian area from January 2009 to October 2021,the 3D V_P,V_S, V_P/V_S structures and seismic locations of the area are obtained by joint inversion using the V_P/V_S model consistency-constrained double-difference tomography method(tomoDDMC).The earthquakes in the study area are mainly concentrated at a depth of 2-6 km,and the focal depth is generally shallow.The Ms 6.0 Luxian earthquake occurred at the transition zone of high-and low-velocity anomalies and the aftershock sequence was distributed along the edge of the low-V_P zone.A small number of foreshocks occurred on the west side of the M_S 6.0 Luxian earthquake,while most of the aftershocks were distributed on the east side of the M_S 6.0 Luxian earthquake.The aftershock sequence consisted of three seismic bands with different trends,and the overall distribution was in a NWW direction,which was inconsistent with the spatial distribution of the main active faults nearby.In addition,the spatiotemporal distribution of earthquakes and the variation of b-values are closely related to the industrial water injection activities in the study area,reflecting the activation of pre-existing hidden faults under certain tectonic and stress environments leading to seismic activities in the area.
文摘Buckling-restrained braces(BRBs)have shown their capability to provide building structures with stiffness,strength,and ductility.Estimating the seismic drifts of buckling-restrained braced frames(BRBFs)is an important design step to control structural and non-structural damage.In current practice of seismic design,the estimation of seismic drifts of BRBFs is performed by using empirical calculations that are independent upon either the type of the structural system or the design level of seismicity.In these empirical calculations,the seismic drifts are estimated by amplifying the reduced elastic drifts obtained under design lateral loading with a displacement amplification factor(DAF).The value of DAF is considered equal to the product of the response modification factor R and the inelastic displacement ratioρ.The goal of the current research is to assess the value ofρfor low-to mid-rise BRBFs designed under low and high levels of seismicity.This goal has been achieved by conducting a series of elastic and inelastic time-history analyses pertaining to an ensemble of earthquake records on 3-,6-and 9-story BRBFs.The results indicate that theρ-ratio increases with an increase in design seismic intensity and an increase in experienced inelasticity.The range ofρfor low seismicity designs ranges from 0.63 to 0.9,while for high seismicity designs this range stretches from 0.83 to 1.29.It has been found that the consideration of a generalρ-ratio of 1.0 is a reasonable estimation for the design of the BRBFs considered in this study.
基金This research has been made possible by contributions from the Natural Sciences and Engineering Research Council(NSERC)/Energi Simulation Industrial Research Chair in Reservoir Simulation and the Alberta Innovates(iCore)Chair in Reservoir ModelingThis research was supported by the Science Foundation of China University of Petroleum,Beijing(No.2462023BJRC001)the National Natural Science Foundation of China Joint Fund Key Support Project(No.U19B6003).
文摘A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distribution and magnitude of hydraulic fracturing-induced seismicity. An integrated machine learning-based investigation is conducted to systematically evaluate multiple factors that contribute to induced seismicity. Feature importance indicates that a distance to fault, a distance to basement, minimum principal stress, cumulative fluid injection, initial formation pressure, and the number of fracturing stages are among significant model predictors. Our seismicity prediction map matches the observed spatial seismicity, and the prediction model successfully guides the fracturing job size of a new well to reduce seismicity risks. This study can apply to mitigating potential seismicity risks in other seismicity-frequent regions.
基金supported by USGS NHERP grant G20AP00039Matched Filter detection was run on the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation (NSF) grant number ACI-1548562it used the Bridges system, which is supported by NSF award number ACI-1445606, at the Pittsburgh Supercomputing Center (PSC).
文摘We present a detailed catalog of 13671 earthquakes in the Eastern Tennessee Seismic Zone(ETSZ)that spans January 1,2005 to July 31,2020.We apply a matched filter detection technique on over 15 years of continuous data,resulting in arguably the most complete catalog of seismicity in the ETSZ yet.The magnitudes of newly detected events are determined by computing the amplitude ratio between the detections and templates using a principal component fit.We also compute the b-value for the new catalog and comparatively relocate a subset of newly detected events using XCORLOC and hypoDD,which shows a more defined structure at depth.We find the greatest concentration along and to the east of the New York-Alabama Lineament,as defined by the magnetic anomaly,supporting the argument that this feature likely is related to the generation of seismicity in the ETSZ.We examine seismicity in the vicinity of the Watts Bar Reservoir,which is located about 5 km from the epicenter of the M_(W) 4.4 December 12,2018 Decatur,Tennessee earthquake,and find possible evidence for reservoir modulated seismicity in this region.We also examine seismicity in the entire ETSZ to search for a correlation between shallow earthquakes and seasonal hydrologic changes.Our results show limited evidence for hydrologicallydriven shallow seismicity due to seasonal groundwater levels in the ETSZ,which contradicts previous studies hypothesizing that most intraplate earthquakes are associated with the dynamics of hydrologic cycles.
文摘Induced seismicity is strongly related to various engineering projects that cause anthropogenic in-situ stress change at a great depth.Hence,there is a need to estimate and mitigate the associated risks.In the past,various simulation methods have been developed and applied to induced seismicity analysis,but there is still a fundamental diference between simulation results and feld observations in terms of the spatial distribution of seismic events and its frequency.The present study aims to develop a method to simulate spatially distributed on-fault seismicity whilst reproducing a complex stress state in the fault zone.Hence,an equivalent continuum model is constructed,based on a discrete fracture network within a fault damage zone,by employing the crack tensor theory.A fault core is simulated at the center of the model as a discontinuous plane.Using the model,a heterogeneous stress state with stress anomalies in the fault zone is frst simulated by applying tractions on the model outer boundaries.Subsequently,the efective normal stress on the fault plane is decreased in a stepwise manner to induce slip.The simulation result is validated in terms of the b-value and other seismic source parameters,hence demonstrating that the model can reproduce spatially and temporally distributed on-fault seismicity.Further analysis on the parameters shows the variation of frequency-magnitude distribution before the occurrence of large seismic events.This variation is found to be consistent with feld observations,thus suggesting the potential use of this simulation method in evaluating the risk for seismic hazards in various engineering projects.
文摘Post shut‐in seismic events in enhanced geothermal systems(EGSs)occur predominantly at the outer rim of the co‐injection seismic cloud.The concept of postinjection fracture and fault closure near the injection well has been proposed and validated as a mechanism for enhancing post shut‐in pressure diffusion that promotes seismic hazard.This phenomenon is primarily attributed to the poro‐elastic closure of fractures resulting from the reduction of wellbore pressure after injection termination.However,the thermal effects in EGSs,mainly including heat transfer and thermal stress,may not be trivial and their role in postinjection fault closure and pressure evolution needs to be explored.In this study,we performed numerical simulations to analyze the relative importance of poro‐elasticity,heat transfer,and thermo‐elasticity in promoting postinjection fault closure and pressure diffusion.The numerical model wasfirst validated against analytical solutions in terms offluid pressure diffusion and against heatedflow‐through experiments in terms of thermal processes.We then quantified and distinguished the contribution of each individual mechanism by comparing four different shut‐in scenarios simulated under different coupled conditions.Our results highlight the importance of poro‐elastic fault closure in promoting postinjection pressure buildup and seismicity,and suggest that heat transfer can further augment the fault closure‐induced pressure increase and thus potentially intensify the postinjection seismic hazard,with minimal contribution from thermo‐elasticity.
文摘The paper is devoted to analysis of hydrogeological, geomagnetic and seismic response to the two great remote geophysical events, 2022 Tonga volcano eruption and 2020-2023 Türkiye earthquakes in Georgia (Caucasus). The geophysical observation system in Georgia, namely, water level stations in the network of deep wells, atmospheric pressure and the geomagnetic sensors of the Dusheti Geophysical Observatory (DGO) as well as seismic data in Garni Observatory (Armenia) respond to the Tonga event by anomalies in the time series. These data show that there are two types of respond: infrasound disturbances in atmospheric pressure and seismic waves in the Earth generated by the eruption. After Tonga eruption January 15 at 04:21 UTC three groups of N-shaped waveforms were registered in the water level corresponding to the global propagation characteristics of the N-shaped waveform of infrasound signals on the barograms generated by eruption at the distance ~15,700 km: they were identified as the Lamb wave, a surface wave package running in the atmosphere with a velocity around ~314 m/s. The paper also presents the WL reactions to three strong EQs that occur in Türkiye 2020-2023, namely Elazığ, Van and Türkiye-Syria EQs. WL in Georgian well network reacts to these events by anomalies of different intensity, which points to the high sensitivity of hydrosphere to remote (several hundred km) strong EQs. The intensity and character of WL reactions depend strongly on the local hydrogeological properties of rocks, surrounding the well.