When a porous rock is subjected to overall compressive loading,either increasing pore pressure or decreasing confining pressure could result in rock failure.The stress path and the applied pressure change rate may aff...When a porous rock is subjected to overall compressive loading,either increasing pore pressure or decreasing confining pressure could result in rock failure.The stress path and the applied pressure change rate may affect the initiation and propagation of fractures within brittle materials.Understanding the physical mechanisms leading to failure is crucial for underground engineering applications and geo-energy exploration and storage.We conducted triaxial compression experiments on porous Bentheim sandstone samples at different stress paths and pressure change rates.First,at a constant confining pressure of 35 MPa and pore pressure of 5 MPa,intact cylindrical samples were axially loaded up to about 85%of the peak strength.Subsequently,the axial piston position was fixed,and then either the pore pressure was increased or the confining pressure was decreased at two different rates(0.5 MPa/min or 2 MPa/min),leading to final catastrophic failure.The mechanical results revealed that samples subjected to higher rates of decreasing effective confining pressure exhibited larger stress drop rates,higher slip rates,higher total breakdown work,higher rates of acoustic emissions(AEs)before failure,and higher post-failure AE decay rates.In contrast,the applied stress path did not significantly affect rock failure characteristics.Comparison of located AE events with post-mortem microstructures of deformed samples shows a good agreement.The AE source type determined from the P-wave first-motion polarity shows that shear failure dominated the fracture process when approaching failure.Gutenberg-Richter b-values revealed a significant decrease before failure in all tests.Our results indicate that,in contrast to the stress path,the rate of effective stress change strongly affects fracturing behavior and AE rate changes.展开更多
To investigate the damage evolution caused by stress-driven and sub-critical crack propagation within the Beishan granite under multi-creep triaxial compressive conditions,the distributed optical fiber sensing and X-r...To investigate the damage evolution caused by stress-driven and sub-critical crack propagation within the Beishan granite under multi-creep triaxial compressive conditions,the distributed optical fiber sensing and X-ray computed tomography were combined to obtain the strain distribution over the sample surface and internal fractures of the samples.The Gini and skewness(G-S)coefficients were used to quantify strain localization during tests,where the Gini coefficient reflects the degree of clustering of elements with high strain values,i.e.,strain localization/delocalization.The strain localization-induced asymmetry of data distribution is quantified by the skewness coefficient.A precursor to granite failure is defined by the rapid and simultaneous increase of the G-S coefficients,which are calculated from strain increment,giving an earlier warning of failure by about 8%peak stress than those from absolute strain values.Moreover,the process of damage accumulation due to stress-driven crack propagation in Beishan granite is different at various confining pressures as the stress exceeds the crack initiation stress.Concretely,strain localization is continuous until brittle failure at higher confining pressure,while both strain localization and delocalization occur at lower confining pressure.Despite the different stress conditions,a similar statistical characteristic of strain localization during the creep stage is observed.The Gini coefficient increases,and the skewness coefficient decreases slightly as the creep stress is below 95%peak stress.When the accelerated strain localization begins,the Gini and skewness coefficients increase rapidly and simultaneously.展开更多
A multi-stage stress relaxation test was performed on a granodiorite sample to understand the deformation process prior to the macroscopic failure of brittle rocks,as well as the transient response during stress relax...A multi-stage stress relaxation test was performed on a granodiorite sample to understand the deformation process prior to the macroscopic failure of brittle rocks,as well as the transient response during stress relaxation.Distributed optical fiber sensing was used to measure strains across the sample surface by helically wrapping the single-mode fiber around the cylindrical sample.Close agreement was observed between the circumferential strains obtained from the optical fibers and the extensometer.The reconstructed full-field strain contours show strain heterogeneity from the crack closure phase,and the strains in the later deformation phase are dominantly localized within the former high-strain zone.The Gini coefficient was used to quantify the degree of strain localization and shows an initial increase during the crack closure phase,a decrease during the linear elastic phase,and a subsequent increase during the post-yielding phase.This behavior corresponds to a process of initial localization from an imperfect boundary condition,homogenization,and eventual relocalization prior to the macroscopic failure of the sample.The transient strain rate decay during the stress relaxation phase was quantified using the p-value in the“Omori-like"power law function.A higher initial stress at the onset of relaxation results in a lower p-value,indicating a slower strain rate decay.As the sample approaches macroscopic failure,the lowest p-value shifts from the most damaged zone to adjacent areas,suggesting stress redistribution or crack propagation in deformed crystalline rocks under stress relaxation conditions.展开更多
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).展开更多
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.展开更多
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.展开更多
文摘When a porous rock is subjected to overall compressive loading,either increasing pore pressure or decreasing confining pressure could result in rock failure.The stress path and the applied pressure change rate may affect the initiation and propagation of fractures within brittle materials.Understanding the physical mechanisms leading to failure is crucial for underground engineering applications and geo-energy exploration and storage.We conducted triaxial compression experiments on porous Bentheim sandstone samples at different stress paths and pressure change rates.First,at a constant confining pressure of 35 MPa and pore pressure of 5 MPa,intact cylindrical samples were axially loaded up to about 85%of the peak strength.Subsequently,the axial piston position was fixed,and then either the pore pressure was increased or the confining pressure was decreased at two different rates(0.5 MPa/min or 2 MPa/min),leading to final catastrophic failure.The mechanical results revealed that samples subjected to higher rates of decreasing effective confining pressure exhibited larger stress drop rates,higher slip rates,higher total breakdown work,higher rates of acoustic emissions(AEs)before failure,and higher post-failure AE decay rates.In contrast,the applied stress path did not significantly affect rock failure characteristics.Comparison of located AE events with post-mortem microstructures of deformed samples shows a good agreement.The AE source type determined from the P-wave first-motion polarity shows that shear failure dominated the fracture process when approaching failure.Gutenberg-Richter b-values revealed a significant decrease before failure in all tests.Our results indicate that,in contrast to the stress path,the rate of effective stress change strongly affects fracturing behavior and AE rate changes.
基金supported by the National Natural Science Foundation of China(Grant No.52339001).
文摘To investigate the damage evolution caused by stress-driven and sub-critical crack propagation within the Beishan granite under multi-creep triaxial compressive conditions,the distributed optical fiber sensing and X-ray computed tomography were combined to obtain the strain distribution over the sample surface and internal fractures of the samples.The Gini and skewness(G-S)coefficients were used to quantify strain localization during tests,where the Gini coefficient reflects the degree of clustering of elements with high strain values,i.e.,strain localization/delocalization.The strain localization-induced asymmetry of data distribution is quantified by the skewness coefficient.A precursor to granite failure is defined by the rapid and simultaneous increase of the G-S coefficients,which are calculated from strain increment,giving an earlier warning of failure by about 8%peak stress than those from absolute strain values.Moreover,the process of damage accumulation due to stress-driven crack propagation in Beishan granite is different at various confining pressures as the stress exceeds the crack initiation stress.Concretely,strain localization is continuous until brittle failure at higher confining pressure,while both strain localization and delocalization occur at lower confining pressure.Despite the different stress conditions,a similar statistical characteristic of strain localization during the creep stage is observed.The Gini coefficient increases,and the skewness coefficient decreases slightly as the creep stress is below 95%peak stress.When the accelerated strain localization begins,the Gini and skewness coefficients increase rapidly and simultaneously.
基金support of her postdoctoral research at the GFZ Helmholtz Centre for Geosciences.P.Pan acknowledges the financial support of the National Natural Science Foundation of China(Grant No.52339001)H.Hofmann and Y.Ji acknowledge the financial support of the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(contract number VH-NG-1516).
文摘A multi-stage stress relaxation test was performed on a granodiorite sample to understand the deformation process prior to the macroscopic failure of brittle rocks,as well as the transient response during stress relaxation.Distributed optical fiber sensing was used to measure strains across the sample surface by helically wrapping the single-mode fiber around the cylindrical sample.Close agreement was observed between the circumferential strains obtained from the optical fibers and the extensometer.The reconstructed full-field strain contours show strain heterogeneity from the crack closure phase,and the strains in the later deformation phase are dominantly localized within the former high-strain zone.The Gini coefficient was used to quantify the degree of strain localization and shows an initial increase during the crack closure phase,a decrease during the linear elastic phase,and a subsequent increase during the post-yielding phase.This behavior corresponds to a process of initial localization from an imperfect boundary condition,homogenization,and eventual relocalization prior to the macroscopic failure of the sample.The transient strain rate decay during the stress relaxation phase was quantified using the p-value in the“Omori-like"power law function.A higher initial stress at the onset of relaxation results in a lower p-value,indicating a slower strain rate decay.As the sample approaches macroscopic failure,the lowest p-value shifts from the most damaged zone to adjacent areas,suggesting stress redistribution or crack propagation in deformed crystalline rocks under stress relaxation conditions.
基金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).
基金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.
基金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.