The Macao Science Satellite-1(known as MSS-1)is the first scientific exploration satellite that was designed to measure the Earth's low latitude magnetic field at high resolution and with high precision by collect...The Macao Science Satellite-1(known as MSS-1)is the first scientific exploration satellite that was designed to measure the Earth's low latitude magnetic field at high resolution and with high precision by collecting data in a near-equatorial orbit.Magnetic field data from MSS-1's onboard Vector Fluxgate Magnetometer(VFM),collected at a sample rate of 50 Hz,allows us to detect and investigate sources of magnetic data contamination,from DC to relevant Nyquist frequency.Here we report two types of artificial disturbances in the VFM data.One is V-shaped events concentrated at night,with frequencies sweeping from the Nyquist frequency down to zero and back up.The other is 5-Hz events(ones that exhibit a distinct 5 Hz spectrum peak);these events are always accompanied by intervals of spiky signals,and are clearly related to the attitude control of the satellite.Our analyses show that VFM noise levels in daytime are systematically lower than in nighttime.The daily average noise levels exhibit a period of about 52 days.The V-shaped events are strongly correlated with higher VFM noise levels.展开更多
INTRODUCTION.On May 1st,2024,around 2:10 a.m.,a catastrophic collapse occurred along the Meilong Expressway near Meizhou City,Guangdong Province,China,at coordinates 24°29′24″N and 116°40′25″E.This colla...INTRODUCTION.On May 1st,2024,around 2:10 a.m.,a catastrophic collapse occurred along the Meilong Expressway near Meizhou City,Guangdong Province,China,at coordinates 24°29′24″N and 116°40′25″E.This collapse resulted in a pavement failure of approximately 17.9 m in length and covering an area of about 184.3 m^(2)(Chinanews,2024).展开更多
Source properties and stress fields are critical to understand fundamental mechanisms for fluid-induced earthquakes.In this study,we identify the focal mechanism solutions(FMSs)of 360 earthquakes with local magnitude ...Source properties and stress fields are critical to understand fundamental mechanisms for fluid-induced earthquakes.In this study,we identify the focal mechanism solutions(FMSs)of 360 earthquakes with local magnitude M_(L)≥1.5 in the Changning shale gas field from January 2016 to May 2017 by fitting three-component waveforms.We then constrain the directions of the maximum horizontal stress(σ_(H_(max)))for four dense earthquake clusters using the stress tensor inversion method.The stress drops of 121 earthquakes with M_(L)≥1.5 are calculated using the spectral ratio method.We examine the spatiotemporal heterogeneity of stress field,and discuss the cause of non-double-couple(non-DC)components in seismicity clusters.Following the Mohr-Coulomb criterion,we estimate the fluid overpressure thresholds from FMS for different seismic clusters,providing insights into potential physical mechanisms for induced seismicity.The FMS results indicate that shallow reverse earthquakes,with steep dip angles,characterize most events.The source mechanisms of earthquakes with M_(L)≥1.5 are dominated by DC components(>70%),but several earthquakes with M_(L)>3.0 and the microseismic events nearby during injection period display significant non-DC components(>30%).Stress inversion results reveal that the σ_(H_(max)) direction ranges from 120°to 128°.Stress drops of earthquakes range between 0.10 and 64.49 MPa,with high values occurring on reverse faults situated at a greater distance from the shale layer,accompanied by a moderate rotation(≤25°)in the trend of σ_(H_(max)).The seismic clusters close to the shale layer exhibit low fluid overpressure thresholds,prone to being triggered by high pore-pressure fluid.The integrated results suggest that the diffusion of high pore pressures is likely to be the primary factor for observed earthquakes.The present results are expected to offer valuable insights into the origin of anomalous seismicity near the shale gas sites.展开更多
This study aims to understand the effect of injection rate on injection-induced fracture activation in granite.We performed water injection-induced slip tests on samples containing either a smooth or a rough fracture ...This study aims to understand the effect of injection rate on injection-induced fracture activation in granite.We performed water injection-induced slip tests on samples containing either a smooth or a rough fracture at four different injection rates under undrained conditions and monitored the acoustic emission(AE)signals during the tests.Experimental results reveal that the critical activation fluid pressure is related to the injection rate,pressure diffusion rate,stress state,and fracture roughness.For the smooth fracture,as the injection rate increases,the critical activation fluid pressure increases significantly,while the injection rate has little effect on the critical activation fluid pressure of the rough fracture.The quasi-static slip distance of fractures decreases as the injection rate increases,with rough fractures exhibiting a greater overall slip distance compared to smooth fractures.The number of AE events per unit sliding distance increases with the injection rate,while the global b value decreases.These results indicate that higher injection rates produce more large-magnitude AE events and more severe slip instability and asperity damage.We established a linkage between fluid injection volume,injection rate,and AE events using the seismogenic index(Σ).The smooth fracture exhibits a steadily increasingΣwith the elapse of injection time,and the rate of increase is higher at higher injection rates;while the rough fracture is featured by a fluctuatingΣ,signifying the intermittent occurrence of large-magnitude AE events associated with the damage of larger fracture asperities.Our results highlight the importance of fracture surface heterogeneity on injection-induced fracture activation and slip.展开更多
In low-permeability geothermal reservoirs,hydro-shearing of pre-existing natural fractures plays a crucial role in improving connectivity between injection and production wells,thereby enhancing heat extraction effici...In low-permeability geothermal reservoirs,hydro-shearing of pre-existing natural fractures plays a crucial role in improving connectivity between injection and production wells,thereby enhancing heat extraction efficiency.This process increases fracture conductivity through dilation caused by injectioninduced slip;however,it also carries the risk of inducing seismic events,posing significant challenges for geothermal operations.This study employs a coupled hydro-mechanical numerical model based on the boundary element method to simulate hydro-shearing under two distinct fluid injection scenarios:(1)monotonic injection and(2)cyclic injection regulated by a traffic light system(TLS).The model assesses the effectiveness of these injection regimes in enhancing fracture conductivity while mitigating seismic hazards.Results indicate that monotonic injection frequently triggers a cascade of seismic events,disrupting pressure and stress distributions on nearby faults and resulting in complex seismic and aseismic interactions.In contrast,TLS-regulated cyclic injection,when carefully managed,promotes stable slip behavior and improves fracture conductivity.This approach proves particularly effective over extended durations during the simultaneous stimulation of two parallel faults.However,in multi-stage stimulation scenariosdwhere natural fractures are stimulated sequentiallydTLS-based cyclic injection,while more efficient at enhancing conductivity,may increase seismicity risk with prolonged application,thereby limiting its safe operational window.展开更多
In subsurface projects where the host rock is of low permeability,fractures play an important role in fluid circulation.Both the geometrical and mechanical properties of the fracture are relevant to the permeability o...In subsurface projects where the host rock is of low permeability,fractures play an important role in fluid circulation.Both the geometrical and mechanical properties of the fracture are relevant to the permeability of the fracture.To evaluate this relationship,we numerically generated self-affine fractures reproducing the scaling relationship of the power spectral density(PSD)of the measured fracture surfaces.The fractures were then subjected to a uniform and stepwise increase in normal stress.A fast Fourier transform(FFT)-based elastic contact model was used to simulate the fracture closure.The evolution of fracture contact area,fracture closure,and fracture normal stiffness were determined throughout the whole process.In addition,the fracture permeability at each step was calculated by the local cubic law(LCL).The influences of roughness exponent and correlation length on the fracture hydraulic and mechanical behaviors were investigated.Based on the power law of normal stiffness versus normal stress,the corrected cubic law and the linear relationship between fracture closure and mechanical aperture were obtained from numerical modeling of a set of fractures.Then,we derived a fracture normal stiffness-permeability equation which incorporates fracture geometric parameters such as the root-mean-square(RMS),roughness exponent,and correlation length,which can describe the fracture flow under an effective medium regime and a percolation regime.Finally,we interpreted the flow transition behavior from the effective medium regime to the percolation regime during fracture closure with the established stiffness-permeability function.展开更多
Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensiona...Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensional(3D)benchmark model of a laboratory experiment on graywacke to examine the dynamic hydraulic fracturing process under a polyaxial stress state.In the numerical model,injection pressures after breakdown(postbreakdown)are varied to study the impact on fracture growth.The fluid pressure front and crack front are identified in the numerical model to analyze the dynamic relationship between fluid diffusion and fracture propagation.Following the hydraulic fracturing test,the polyaxial stresses are rotated to investigate the influence of the stress field rotation on the fracture slip behavior and permeability.The results show that fracture propagation guides fluid diffusion under a high postbreakdown injection pressure.The crack front runs ahead of the fluid pressure front.Under a low postbreakdown injection pressure,the fluid pressure front gradually reaches the crack front,and fluid diffusion is the main driving factor of fracture propagation.Under polyaxial stress conditions,fluid injection not only opens tensile fractures but also induces hydroshearing.When the polyaxial stress is rotated,the fracture slip direction of a fully extended fracture is consistent with the shear stress direction.The fracture slip direction of a partly extended fracture is influenced by the increase in shear stress.Normal stress affects the permeability evolution by changing the average mechanical aperture.Shear stress can induce shearing and sliding on the fracture plane,thereby increasing permeability.展开更多
Purpose–This study solves the key problem that the static level monitoring is susceptible to temperature interference and affects the accuracy in slope instability/deformation monitoring.The purpose is to develop a r...Purpose–This study solves the key problem that the static level monitoring is susceptible to temperature interference and affects the accuracy in slope instability/deformation monitoring.The purpose is to develop a reliable temperature compensation method for the system,improve the accuracy of slope stability monitoring and provide support for improving the safety and safety monitoring of engineering spoil slope and other projects.Design/methodology/approach–Combined with theoretical analysis and experimental verification,the temperature compensation method is explored.The working principle of the hydrostatic leveling monitoring system is analyzed and the data processing formula,the temperature error calculation formula and the calculation formula for eliminating the error settlement value are derived.The temperature compensation method is established and verified by the field test of the engineering spoil slope which is disturbed by a debris flow.Findings–The experimental results show that this method can reduce the error of the static level monitoring system by about 40%.The field test shows that the fluctuation of slope settlement monitoring value is reduced after temperature compensation and the monitoring value is consistent with the actual situation,which has certain practicability.Originality/value–The originality of this study is to derive a theoretical formula for quantifying/eliminating temperature errors in static leveling and to establish a practical temperature compensation method.The accuracy of the system is improved,which provides a reference for slope stability monitoring under complex environment(especially railway geotechnical engineering)and promotes the development of precision monitoring technology.展开更多
0 INTRODUCTION Synthetic Aperture Radar(SAR)remote sensing,particularly with the C-band Sentinel-1 mission,has been widely used for landslide displacement analysis due to its high spatial resolution and revisit freque...0 INTRODUCTION Synthetic Aperture Radar(SAR)remote sensing,particularly with the C-band Sentinel-1 mission,has been widely used for landslide displacement analysis due to its high spatial resolution and revisit frequency(Zhou et al.,2024;Dai et al.,2021).However,in densely vegetated or humid mountainous regions such as the Three Gorges Reservoir(TGR),C-band signals suffer from temporal decorrelation,limiting their effectiveness for landslide monitoring.展开更多
In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events cau...In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events caused by fault slip and their potential effects on rockbursts.Through Bayesian inversion,it is determined that the sources near fault FQ14 have a significant shear component.Additionally,we analyzed the stress and displacement fields of high-energy events,along with the hypocenter distribution of aftershocks,which aided in identifying the slip direction of the critically stressed fault FQ14.We also performed forward modeling to capture the complex dynamics of fault slip under varying friction laws and shear fracture modes.The selection of specific friction laws for fault slip models was based on their ability to accurately replicate observed slip behavior under various external loading conditions,thereby enhancing the applicability of our findings.Our results suggest that the slip behavior of fault FQ14 can be effectively understood by comparing different scenarios.展开更多
基金supported by the National Key R&D Program of China(Grant2022YFF0503700)the National Natural Science Foundation of China(42474200 and 42174186)。
文摘The Macao Science Satellite-1(known as MSS-1)is the first scientific exploration satellite that was designed to measure the Earth's low latitude magnetic field at high resolution and with high precision by collecting data in a near-equatorial orbit.Magnetic field data from MSS-1's onboard Vector Fluxgate Magnetometer(VFM),collected at a sample rate of 50 Hz,allows us to detect and investigate sources of magnetic data contamination,from DC to relevant Nyquist frequency.Here we report two types of artificial disturbances in the VFM data.One is V-shaped events concentrated at night,with frequencies sweeping from the Nyquist frequency down to zero and back up.The other is 5-Hz events(ones that exhibit a distinct 5 Hz spectrum peak);these events are always accompanied by intervals of spiky signals,and are clearly related to the attitude control of the satellite.Our analyses show that VFM noise levels in daytime are systematically lower than in nighttime.The daily average noise levels exhibit a period of about 52 days.The V-shaped events are strongly correlated with higher VFM noise levels.
基金supported by the National Natural Science Foundation of China(Nos.42371094,41907253)partially supported by the Interdisciplinary Cultivation Program of Xidian University(No.21103240005)the Postdoctoral Fellowship Program of CPSF(No.GZB20240589)。
文摘INTRODUCTION.On May 1st,2024,around 2:10 a.m.,a catastrophic collapse occurred along the Meilong Expressway near Meizhou City,Guangdong Province,China,at coordinates 24°29′24″N and 116°40′25″E.This collapse resulted in a pavement failure of approximately 17.9 m in length and covering an area of about 184.3 m^(2)(Chinanews,2024).
基金financially supported by the National Natural Science Foundation of China(Grant Nos.U20A20266 and 12302503)Scientific and technological research projects in Sichuan province(Grant No.2024NSFSC0973).
文摘Source properties and stress fields are critical to understand fundamental mechanisms for fluid-induced earthquakes.In this study,we identify the focal mechanism solutions(FMSs)of 360 earthquakes with local magnitude M_(L)≥1.5 in the Changning shale gas field from January 2016 to May 2017 by fitting three-component waveforms.We then constrain the directions of the maximum horizontal stress(σ_(H_(max)))for four dense earthquake clusters using the stress tensor inversion method.The stress drops of 121 earthquakes with M_(L)≥1.5 are calculated using the spectral ratio method.We examine the spatiotemporal heterogeneity of stress field,and discuss the cause of non-double-couple(non-DC)components in seismicity clusters.Following the Mohr-Coulomb criterion,we estimate the fluid overpressure thresholds from FMS for different seismic clusters,providing insights into potential physical mechanisms for induced seismicity.The FMS results indicate that shallow reverse earthquakes,with steep dip angles,characterize most events.The source mechanisms of earthquakes with M_(L)≥1.5 are dominated by DC components(>70%),but several earthquakes with M_(L)>3.0 and the microseismic events nearby during injection period display significant non-DC components(>30%).Stress inversion results reveal that the σ_(H_(max)) direction ranges from 120°to 128°.Stress drops of earthquakes range between 0.10 and 64.49 MPa,with high values occurring on reverse faults situated at a greater distance from the shale layer,accompanied by a moderate rotation(≤25°)in the trend of σ_(H_(max)).The seismic clusters close to the shale layer exhibit low fluid overpressure thresholds,prone to being triggered by high pore-pressure fluid.The integrated results suggest that the diffusion of high pore pressures is likely to be the primary factor for observed earthquakes.The present results are expected to offer valuable insights into the origin of anomalous seismicity near the shale gas sites.
基金supported by the National Key Research and Development Program of China(Grant No.2023YFB2390300)the National Natural Science Foundation of China(Grant No.42172292)Yinlin Ji is grateful for the support by the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(Contract No.VH-NG-1516).
文摘This study aims to understand the effect of injection rate on injection-induced fracture activation in granite.We performed water injection-induced slip tests on samples containing either a smooth or a rough fracture at four different injection rates under undrained conditions and monitored the acoustic emission(AE)signals during the tests.Experimental results reveal that the critical activation fluid pressure is related to the injection rate,pressure diffusion rate,stress state,and fracture roughness.For the smooth fracture,as the injection rate increases,the critical activation fluid pressure increases significantly,while the injection rate has little effect on the critical activation fluid pressure of the rough fracture.The quasi-static slip distance of fractures decreases as the injection rate increases,with rough fractures exhibiting a greater overall slip distance compared to smooth fractures.The number of AE events per unit sliding distance increases with the injection rate,while the global b value decreases.These results indicate that higher injection rates produce more large-magnitude AE events and more severe slip instability and asperity damage.We established a linkage between fluid injection volume,injection rate,and AE events using the seismogenic index(Σ).The smooth fracture exhibits a steadily increasingΣwith the elapse of injection time,and the rate of increase is higher at higher injection rates;while the rough fracture is featured by a fluctuatingΣ,signifying the intermittent occurrence of large-magnitude AE events associated with the damage of larger fracture asperities.Our results highlight the importance of fracture surface heterogeneity on injection-induced fracture activation and slip.
基金the financial support of the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(Contract number VHNG-1516).
文摘In low-permeability geothermal reservoirs,hydro-shearing of pre-existing natural fractures plays a crucial role in improving connectivity between injection and production wells,thereby enhancing heat extraction efficiency.This process increases fracture conductivity through dilation caused by injectioninduced slip;however,it also carries the risk of inducing seismic events,posing significant challenges for geothermal operations.This study employs a coupled hydro-mechanical numerical model based on the boundary element method to simulate hydro-shearing under two distinct fluid injection scenarios:(1)monotonic injection and(2)cyclic injection regulated by a traffic light system(TLS).The model assesses the effectiveness of these injection regimes in enhancing fracture conductivity while mitigating seismic hazards.Results indicate that monotonic injection frequently triggers a cascade of seismic events,disrupting pressure and stress distributions on nearby faults and resulting in complex seismic and aseismic interactions.In contrast,TLS-regulated cyclic injection,when carefully managed,promotes stable slip behavior and improves fracture conductivity.This approach proves particularly effective over extended durations during the simultaneous stimulation of two parallel faults.However,in multi-stage stimulation scenariosdwhere natural fractures are stimulated sequentiallydTLS-based cyclic injection,while more efficient at enhancing conductivity,may increase seismicity risk with prolonged application,thereby limiting its safe operational window.
基金supported by the China Postdoctoral Science Foundation Funded Project(Grant No.2023M740385)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20233326)the support by the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(Contract No.VH-NG-1516).
文摘In subsurface projects where the host rock is of low permeability,fractures play an important role in fluid circulation.Both the geometrical and mechanical properties of the fracture are relevant to the permeability of the fracture.To evaluate this relationship,we numerically generated self-affine fractures reproducing the scaling relationship of the power spectral density(PSD)of the measured fracture surfaces.The fractures were then subjected to a uniform and stepwise increase in normal stress.A fast Fourier transform(FFT)-based elastic contact model was used to simulate the fracture closure.The evolution of fracture contact area,fracture closure,and fracture normal stiffness were determined throughout the whole process.In addition,the fracture permeability at each step was calculated by the local cubic law(LCL).The influences of roughness exponent and correlation length on the fracture hydraulic and mechanical behaviors were investigated.Based on the power law of normal stiffness versus normal stress,the corrected cubic law and the linear relationship between fracture closure and mechanical aperture were obtained from numerical modeling of a set of fractures.Then,we derived a fracture normal stiffness-permeability equation which incorporates fracture geometric parameters such as the root-mean-square(RMS),roughness exponent,and correlation length,which can describe the fracture flow under an effective medium regime and a percolation regime.Finally,we interpreted the flow transition behavior from the effective medium regime to the percolation regime during fracture closure with the established stiffness-permeability function.
基金supported by the Knowledge Innovation Program of Wuhan-Basic Research (Grant No.2022010801010159)support from the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES (Contract number VH-NG-1516)supported by the Swedish Radiation Safety Authority (Project SSM2020-2758).
文摘Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensional(3D)benchmark model of a laboratory experiment on graywacke to examine the dynamic hydraulic fracturing process under a polyaxial stress state.In the numerical model,injection pressures after breakdown(postbreakdown)are varied to study the impact on fracture growth.The fluid pressure front and crack front are identified in the numerical model to analyze the dynamic relationship between fluid diffusion and fracture propagation.Following the hydraulic fracturing test,the polyaxial stresses are rotated to investigate the influence of the stress field rotation on the fracture slip behavior and permeability.The results show that fracture propagation guides fluid diffusion under a high postbreakdown injection pressure.The crack front runs ahead of the fluid pressure front.Under a low postbreakdown injection pressure,the fluid pressure front gradually reaches the crack front,and fluid diffusion is the main driving factor of fracture propagation.Under polyaxial stress conditions,fluid injection not only opens tensile fractures but also induces hydroshearing.When the polyaxial stress is rotated,the fracture slip direction of a fully extended fracture is consistent with the shear stress direction.The fracture slip direction of a partly extended fracture is influenced by the increase in shear stress.Normal stress affects the permeability evolution by changing the average mechanical aperture.Shear stress can induce shearing and sliding on the fracture plane,thereby increasing permeability.
基金funded by the Scientific Research Project of China Academy of Railway Sciences Group Co.,Ltd(No.2024YJ332 and No.2024QT005)Scientific Research Special Project of China State Railway Group Co.,Ltd(No.TICSTR-2024-Ⅳ-007).
文摘Purpose–This study solves the key problem that the static level monitoring is susceptible to temperature interference and affects the accuracy in slope instability/deformation monitoring.The purpose is to develop a reliable temperature compensation method for the system,improve the accuracy of slope stability monitoring and provide support for improving the safety and safety monitoring of engineering spoil slope and other projects.Design/methodology/approach–Combined with theoretical analysis and experimental verification,the temperature compensation method is explored.The working principle of the hydrostatic leveling monitoring system is analyzed and the data processing formula,the temperature error calculation formula and the calculation formula for eliminating the error settlement value are derived.The temperature compensation method is established and verified by the field test of the engineering spoil slope which is disturbed by a debris flow.Findings–The experimental results show that this method can reduce the error of the static level monitoring system by about 40%.The field test shows that the fluctuation of slope settlement monitoring value is reduced after temperature compensation and the monitoring value is consistent with the actual situation,which has certain practicability.Originality/value–The originality of this study is to derive a theoretical formula for quantifying/eliminating temperature errors in static leveling and to establish a practical temperature compensation method.The accuracy of the system is improved,which provides a reference for slope stability monitoring under complex environment(especially railway geotechnical engineering)and promotes the development of precision monitoring technology.
基金supported by the National Natural Science Foundation of China(Nos.42371094,41907253)the Fundamental Research Funds for the Central Universities(No.B250201054).
文摘0 INTRODUCTION Synthetic Aperture Radar(SAR)remote sensing,particularly with the C-band Sentinel-1 mission,has been widely used for landslide displacement analysis due to its high spatial resolution and revisit frequency(Zhou et al.,2024;Dai et al.,2021).However,in densely vegetated or humid mountainous regions such as the Three Gorges Reservoir(TGR),C-band signals suffer from temporal decorrelation,limiting their effectiveness for landslide monitoring.
基金the Graduate Innovation Program of China University of Mining and Technology,the Fundamental Research Funds for the Central Universities(Grant No.2023WLKXJ017)the Postgraduate Research&Practice Innovation Program of Jiangsu Province(Grant No.KYCX23_2776)the Shandong Energy Group(Grant No.SNKJ2022BJ03-R28)。
文摘In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events caused by fault slip and their potential effects on rockbursts.Through Bayesian inversion,it is determined that the sources near fault FQ14 have a significant shear component.Additionally,we analyzed the stress and displacement fields of high-energy events,along with the hypocenter distribution of aftershocks,which aided in identifying the slip direction of the critically stressed fault FQ14.We also performed forward modeling to capture the complex dynamics of fault slip under varying friction laws and shear fracture modes.The selection of specific friction laws for fault slip models was based on their ability to accurately replicate observed slip behavior under various external loading conditions,thereby enhancing the applicability of our findings.Our results suggest that the slip behavior of fault FQ14 can be effectively understood by comparing different scenarios.
