Strong seismic excitation and fault dislocation are likely to occur simultaneously in high-intensity seismic zones,causing severe damage to tunnels crossing active fault zones.This paper aims to develop a novel analyt...Strong seismic excitation and fault dislocation are likely to occur simultaneously in high-intensity seismic zones,causing severe damage to tunnels crossing active fault zones.This paper aims to develop a novel analytical solution to determine the longitudinal mechanical responses of tunnels subjected to the combined effects of seismic waves and strike-slip faulting.Adopting the elastic springbeam model,the seismic waves are modelled as shear horizontal(SH)waves and the fault dislocation follows an S-shaped pattern;the superposition principle for free-fielddisplacements caused by both effects is assumed.In addition,the transmission and reflectionof seismic waves at the fault-rock geological interface and the tangential contact conditions at the tunnel-rock interface are considered.The analytical model is validated against numerical simulations,confirmingits accuracy in calculating tunnel responses.Moreover,a parametric study is conducted to evaluate the impact of key factors,including fault displacement,fault zone width,fault dip angle,earthquake frequency,rock conditions,tunnel lining stiffness,and tangential contact conditions,on tunnel responses.Compared with each effect alone,the combined effects of seismic waves and strike-slip faulting significantlychange the tunnel deformation and internal forces,leading to increased tunnel responses,especially within the fault zone and near the fault-rock interfaces.Depending on specificparameters,tunnel responses can be classifiedinto seismic-dominated,faulting-dominated,and seismic-faulting coupled responses on the basis of the relative contributions of each effect.The proposed analytical solution can be applied to quickly predict the longitudinal mechanical behaviour of tunnels under such combined effects in engineering applications.展开更多
The widely distributed loess deposits in the Yellow River Basin exhibit unique engineering geological characteristics.The variations in their thickness and stratigraphic structure significantly amplify ground motion p...The widely distributed loess deposits in the Yellow River Basin exhibit unique engineering geological characteristics.The variations in their thickness and stratigraphic structure significantly amplify ground motion parameters,directly influencing the regional seismic hazard risk level.This study methodically conducted on-site studies and observations of building collapses and damages resulting from seismic amplification effects,using the Wenchuan M_(S)8.0 earthquake as a case study.Comprehensive experimental and numerical simulation studies were carried out.A large-scale shaking table test was performed,and numerical models for 14 different loess sites types were established.Various types of seismic waves were incorporated into these models for systematic numerical simulation calculations.The research reveals the mechanisms by which loess deposit thickness and stratigraphic structure in the Yellow River Basin affect seismic ground motion amplification.The results indicate that as the epicentral distance increases,the peak ground motion shows a marked attenuation trend,with the horizontal component attenuating substantially faster than the vertical component.As the overlying loess layer thickness increases from 50 to 100 m,the seismic intensity may escalate by 3−4 degrees,and the peak acceleration may amplify by 1.5−2.2 times.With the augmentation of loess deposit thickness and the proliferation of soil layers,both the peak acceleration response spectrum and the characteristic period demonstrate an upward tendency,exhibiting slight fluctuations contingent upon the seismic wave type.展开更多
The northern segment of the North-South Seismic Belt is characterized by intense crustal deformation,well-developed active tectonics,and frequent occurrences of strong earthquakes.Therefore,conducting a Probabilistic ...The northern segment of the North-South Seismic Belt is characterized by intense crustal deformation,well-developed active tectonics,and frequent occurrences of strong earthquakes.Therefore,conducting a Probabilistic Seismic Hazard Analysis(PSHA)for this region is of significant importance for supporting seismic fortification in major engineering projects and formulating disaster prevention and mitigation policies.In this study,a composite seismic source model was constructed by integrating data on historical earthquakes,active faults,and paleoseismicity.Furthermore,a logic tree framework was employed to quantify epistemic uncertainties,enabling a systematic seismic hazard assessment of the region.To more accurately characterize the spatial heterogeneity of seismic activity,improvements were made to both the Circular Spatial Smoothing Model(CSSM)with a fixed radius and the Adaptive Spatial Smoothing Model(ASSM),with full consideration given to the spatiotemporal completeness of historical earthquake magnitudes.Regarding the CSSM,for scenarios involving small sample sizes in earthquake catalogs,the cross-validation method proposed in this study demonstrated higher robustness than the maximum likelihood method in determining the optimal correlation distance.Performance evaluation results indicate that while both models effectively characterize seismic activity,the ASSM exhibits superior overall predictive performance compared to the CSSM,owing to its ability to adaptively adjust the smoothing radius according to seismic density.Significant discrepancies were observed in the Peak Ground Acceleration(PGA)results calculated with a 10%probability of exceedance in 50 years across different combinations of seismic source models.The single spatially smoothed point-source model yielded a maximum PGA of approximately 0.52 g,with high-value areas concentrated near historical epicenters,thereby significantly underestimating the hazard associated with major fault zones.When combined with the simple fault-source model,the maximum PGA increased to 0.8 g,with high-value zones exhibiting a zonal distribution along faults;however,the risk remained underestimated for faults with low slip rates that are nevertheless approaching their recurrence cycles.Following the introduction of the time-dependent characteristic fault-source model,local PGA values for faults in the middle-to-late stages of their recurrence cycles increased by a factor of 2 to 7 compared to the single model.These results demonstrate that the characteristic fault-source model reasonably delineates the time-dependence of large earthquake recurrence,thereby providing a more accurate assessment of imminent seismic risks.By comprehensively applying the improved spatially smoothed pointsource model,the simple fault-source model,and the characteristic fault-source model,the following faults within the region were identified as having high seismic hazard:the Huangxianggou,Zhangxian,and Tianshui segments of the Xiqinling northern edge fault;the Maqin-Maqu segment of the Dongkunlun fault;the Longriqu fault;the Maoergai fault;the Elashan fault;the Riyueshan fault;the eastern segment of the Lenglongling fault;the Maxianshan segment of the Maxianshan northern Margin fault;and the Maomaoshan-Jinqianghe segment of the Laohushan-Maomaoshan fault.As these faults are located within seismic gaps or are approaching the recurrence periods of large earthquakes,they should be prioritized for current and future seismic monitoring as well as disaster prevention and mitigation efforts.展开更多
Topological insulators with localized edge or interface states have been extensively studied,particularly in phononic crystals and related fields;however,their application in seismic metamaterials remains largely unex...Topological insulators with localized edge or interface states have been extensively studied,particularly in phononic crystals and related fields;however,their application in seismic metamaterials remains largely unexplored.To address this gap,we designed a topological seismic metamaterial,where the topological interface is formed by joining the ends of two distinct one-dimensional periodic lattices.The first full-scale field experiment confirms the existence of topological interface states,which exhibit pronounced localization characteristics and induce a resonant amplification effect of 7.2 dB on the total energy of seismic surface waves.This study provides the first experimental validation for the implementation of topological principles in the design of seismic metamaterials,enabling novel approaches to high-sensitivity seismic detection and efficient energy localization for wave control.展开更多
The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velo...The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velocity models are available for the LMS fault zone,high-resolution velocity models are lacking.Therefore,a dense array of 240 short-period seismometers was deployed around the central segment of the LMS fault zone for approximately 30 days to monitor earthquakes and characterize fine structures of the fault zone.Considering the large quantity of observed seismic data,the data processing workflow consisted of deep learning-based automatic earthquake detection,phase arrival picking,and association.Compared with the earthquake catalog released by the China Earthquake Administration,many more earthquakes were detected by the dense array.Double-difference seismic tomography was adopted to determine V_(p),V_(s),and V_(p)/V_(s)models as well as earthquake locations.The checkerboard test showed that the velocity models have spatial resolutions of approximately 5 km in the horizontal directions and 2 km at depth.To the west of the Yingxiu–Beichuan Fault(YBF),the Precambrian Pengguan complex,where most of earthquakes occurred,is characterized by high velocity and low V_(p)/V_(s)values.In comparison,to the east of the YBF,the Upper Paleozoic to Jurassic sediments,where few earthquakes occurred,show low velocity and high V_(p)/V_(s)values.Our results suggest that the earthquake activity in the LMS fault zone is controlled by the strength of the rock compositions.When the high-resolution velocity models were combined with the relocated earthquakes,we were also able to delineate the fault geometry for different faults in the LMS fault zone.展开更多
The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismi...The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismic wave simulation with a data-constrained finite-fault rupture model.The constraint is implemented by identifying the optimal ground motion models(GMMs)through a scoring system that selects the best-fit GMMs to mid-and far-field China Earthquake Networks Center(CENC)seismic network data;and applying the optimal GMMs to refine the rupture model parameters for near-fault intensity field simulation.The simulated near-fault seismic intensity field reproduces seismic intensities collected from Myanmar’s sparse seismic network and concentrated in≥Ⅷintensity zones within 50 km of the projected fault plane;and identifies abnormal intensity regions exhibiting≥Ⅹintensity along the Meiktila-Naypyidaw corridor and near Shwebo that are attributed to soft soil amplification effects and near-fault directivity.This framework can also be applied to post-earthquake assessments in other similar regions.展开更多
Studying the inner structure of intraplate earthquakes originating in aseismic areas,which are poorly covered by seismic networks or as historical earthquakes is usually the only way to get knowledge about their sourc...Studying the inner structure of intraplate earthquakes originating in aseismic areas,which are poorly covered by seismic networks or as historical earthquakes is usually the only way to get knowledge about their source mechanisms,which is partially essential for a deeper understanding of intraplate geodynamics.The epicentral zones of earthquakes are situated in hard-toreach areas,so,using active seismic methods for such purposes is unreasonable or even impossible because of high cost and logistical difficulties.We propose a novel approach that combines diverse passive seismic methods,which allows us to get sufficient information about geological environment structure for such task solutions.As an example,we investigated the inner structure of platform earthquake epicentral zone originated up north of Russia.We used four passive seismic methods:microseismic sounding method,passive seismic interferometry,HVSR method,and microseismic activity method.We show that passive seismic data,recoded in the same installation and processed by these different methods,can provide sufficient information about structure of studied environment,needed to explain source mechanism.In sum,the hypocenter zone is presented by intersection of vertical faults and a lateral fractured zone in the middle crust.Results were confirmed by comparison with results by active seismic methods.展开更多
This article aims to enhance seismic hazard assessment methods for Kazakhstan’s seismotectonic conditions.It combines probabilistic seismic hazard analysis(PSHA),ground motion simulation,sitespecific geological and g...This article aims to enhance seismic hazard assessment methods for Kazakhstan’s seismotectonic conditions.It combines probabilistic seismic hazard analysis(PSHA),ground motion simulation,sitespecific geological and geotechnical data analysis,and seismic scenario analysis to develop Probabilistic General Seismic Zoning(GSZ)maps for Kazakhstan and Probabilistic Seismic Microzoning maps for Almaty.These maps align with Eurocode 8 principles,incorporating seismic intensity and engineering parameters like peak ground acceleration(PGA).The new procedure,applied in national projects,has resulted in GSZ maps for the country,seismic microzoning maps for Almaty,and detailed seismic zoning maps for East Kazakhstan.These maps,part of a regulatory document,guide earthquake-resistant design and construction.They offer a comprehensive assessment of seismic hazards,integrating traditional Medvedev-Sponheuer-Karnik(MSK-64)intensity scale points with quantitative parameters like peak ground acceleration.This innovative approach promises to advance methods for quantifying seismic hazards in specific regions.展开更多
A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the ...A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the studied area,and(ii)to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate.Despite the limited number of stations and the short-period char-acter of the seismic sensors,the empirical Green's functions were retrieved for all station pairs using two years of passive data.Both group and phase velocities were derived,the former using the widespread frequency-time analysis,the latter through the analysis of the real part of the cross-spectra.The main advantage of combining these two methods is a more accurate identification of higher modes,resulting in a reduction of ambiguity during picking and data interpretation.Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths.The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method,for this specific case.The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast.Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs.While the average model can be used for a more accurate location of the local natural seismicity,the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth.Finally,the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field.Errors as large as 4%have been found,related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources,and to the very short time widows.It was shown that using a one-year time window these errors arereduced to 0.3%.展开更多
The Chilean subduction zone is one of the most seismically active regions globally,characterized by extensive intermediatedepth seismicity in the slab.In this study,we construct a new earthquake catalog for northern C...The Chilean subduction zone is one of the most seismically active regions globally,characterized by extensive intermediatedepth seismicity in the slab.In this study,we construct a new earthquake catalog for northern Chile using seismic waveforms assembled for the period of 2014-2019,from which 320,070 P-wave and 232,907 S-wave first arrivals are obtained for 25,763 earthquakes.Grid search location method NonLinLoc is applied to determine initial earthquake locations and double-difference location method is used to improve relative event locations.The distribution of earthquakes exhibits distinct patterns to the north and south of 21°S.There are many more earthquakes deeper than~150 km to the south of 21°S,while relatively fewer to the north.The intraslab earthquakes shallower than~80 km generally reveal a distinct double seismic zone,and the gap between the two seismic planes disappears at a depth of approximately~80 km,followed by a concentration of seismicity in the depth range of~80-150 km.In the deeper slab,there exist several seismicity clusters with distinct earthquake activities down to~300 km.These characteristics shown in slab seismicity are likely caused by different mechanisms and can be helpful for understanding the subduction process.展开更多
The quantitative assessment of seismic hazard of Uzbekistan has been examined,and new maps of seismic zoning has been developed.Quasihomogeneous seismological areas and seism-generating zones are considered as seismic...The quantitative assessment of seismic hazard of Uzbekistan has been examined,and new maps of seismic zoning has been developed.Quasihomogeneous seismological areas and seism-generating zones are considered as seismic sources,based on the analysis of seismotectonic data.The seismological parametrization of seismic sources has been carried out,including determination of parameters of earthquake’s reoccurrence for different power levels,seismic potential and a predominant type motion in the source of occurring earthquakes.The reoccurrence parameters of seismic sources were both determined by directly from Gutenberg-Richter dependence law and using summation and distribution methods.Setting were conducted separately for a sample of strong(M≥5)and weak(M<5)earthquakes.The seismic potential Mmax of seismic active zones was determined by seismological and seismotectonic methods.The predominant type of movement in the faulting for each seismic source is reverse fault.Regional regularities of seismic intensity attenuation with distance of different power levels are determined.The seismic hazard is expressed in points of a macro seismic scale,in velocities and accelerations of ground motions,and it is characterized by calculated seismic intensity with the set probability(P=0.9,P=0.95,P=0.98 and P=0.99)not exceed within 50 years in the constructed maps.In seismic zoning map,the seismic intensity was in average soil conditions.The developed complex of maps considers a number of uncertainties of the input parameters in relation to both incompleteness of initial seismological and seismotectonic data,as well the probabilistic nature of seismic process and seismic intensity.The factors of uncertainty of the input parameters are taken into account by constructing a logic tree.The constructed maps of the general seismic zoning are intended for imple menting antiseismic actions in Uzbekistan.展开更多
Earthquakes pose significant perils to the built environment in urban areas.To avert the calamitous aftermath of earthquakes,it is imperative to construct seismic resilient cities.Due to the intricacy of the concept o...Earthquakes pose significant perils to the built environment in urban areas.To avert the calamitous aftermath of earthquakes,it is imperative to construct seismic resilient cities.Due to the intricacy of the concept of urban seismic resilience(USR),its assessment is a large-scale system engineering issue.The assessment of USR should be based on the notion of urban seismic capacity(USC)assessment,which includes casualties,economic loss,and recovery time as criteria.Functionality loss is also included in the assessment of USR in addition to these criteria.The assessment indicator system comprising five dimensions(building and lifeline infrastructure,environment,society,economy,and institution)and 20 indicators has been devised to quantify USR.The analytical hierarchy process(AHP)is utilized to compute the weights of the criteria,dimensions,and indicators in the urban seismic resilience assessment(USRA)indicator system.When the necessary data for a city are obtainable,the seismic resilience of that city can be assessed using this framework.To illustrate the proposed methodology,a moderate-sized city in China was selected as a case study.The assessment results indicate a high level of USR,suggesting that the city possesses strong capabilities to withstand and recover from potential future earthquakes.展开更多
Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mi...Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mitigation. Key controls on seismicity are defined in terms of fault and fracture strength, second-order frictional response and stability, and competing fluid-driven mechanisms for arrest. We desire to constrain maximum event magnitudes in triggered earthquakes by relating pre-existing critical stresses to fluid injection volume to explain why some recorded events are significantly larger than anticipated seismic moment thresholds. This formalism is consistent with several uncharacteristically large fluid injection-triggered earthquakes. Such methods of reactivating fractures and faults by hydraulic stimulation in shear or tensile fracturing are routinely used to create permeability in the subsurface. Microearthquakes (MEQs) generated by such stimulations can be used to diagnose permeability evolution. Although high-fidelity data sets are scarce, the EGS-Collab and Utah FORGE hydraulic stimulation field demonstration projects provide high-fidelity data sets that concurrently track permeability evolution and triggered seismicity. Machine learning deciphers the principal features of MEQs and the resulting permeability evolution that best track permeability changes – with transfer learning methods allowing robust predictions across multiple eological settings. Changes in permeability at reactivated fractures in both shear and extensional modes suggest that permeability change (Δk) scales with the seismic moment (M) of individual MEQs as Δk∝M. This scaling relation is exact at early times but degrades with successive MEQs, but provides a method for characterizing crustal permeability evolution using MEQs, alone. Importantly, we quantify for the first time the role of prestress in defining the elevated magnitude and seismic moment of fluid injection-triggered events, and demonstrate that such MEQs can also be used as diagnostic in quantifying permeability evolution in the crust.展开更多
The occurrence time and magnitude of injection-induced seismicity are influenced by engineering factors,such as wellhead pressure,injection location,injection volume,and injection rate.Understanding the relationship b...The occurrence time and magnitude of injection-induced seismicity are influenced by engineering factors,such as wellhead pressure,injection location,injection volume,and injection rate.Understanding the relationship between injection operations and seismic magnitude is of great significance for optimizing industrial production and reducing earthquake disasters.Numerical simulation of hydromechanical coupling is a crucial method for studying injection-induced seismicity.However,few studies have explored the risk management measures for injection-induced seismicity from the perspective of engineering.How seismic magnitudes can be reduced through reasonable adjustments to injection operations in engineering remains unclear.Therefore,in this study,a 3D hydro-mechanical coupling model involving multiple faults and injection wells was established based on the geological background and well location of Fox Creek,Canada.Different injection schemes under multi-well and multi-fault conditions were studied,and a traffic light system was used to simulate and control the magnitudes under a multi-well injection scheme.Specifically,we simulated injection scenarios involving up to three wells and analyzed the response of five faults.We compared the maximum moment magnitude of different scenarios by controlling the same injection volume.The results revealed the effect and advantage of the multi-well scheme in reducing seismic magnitude.To reduce the risk of induced seismicity,utilizing far-fault operational wells to compensate for the effects of near-fault operational wells proves to be an efficient and cost-effective method,with potential for wide practical applications.展开更多
A seismic swarm occurred southeast of King George Island,South Shetland Islands,Antarctica,between August 2020 and February 2021.This work intends to parameterize seismic events recorded by seismic station AM.R4DE2 fr...A seismic swarm occurred southeast of King George Island,South Shetland Islands,Antarctica,between August 2020 and February 2021.This work intends to parameterize seismic events recorded by seismic station AM.R4DE2 from 15 September to 31 October 2020.Using the localization methodology with a single station,the record of the entire period was analyzed manually to determine the local magnitude,hypocentral distance,epicentral distance,backazimuth,and location of the epicenter for each event.We could parameterize 6362 events,although we estimate the occurrence to be around 20000 for the period.The results suggest a magmatic origin for the swarm,supporting previous studies.Seismicity exhibited a southeastward migration away from King George Island,as indicated by a progressive increase in epicentral distance over time.Most events were classified as volcanic and volcano-tectonic,supporting a magmatogenesis hypothesis linked to the opening of Bransfield Ridge.展开更多
Mining-related seismicity poses significant challenges in underground coal mining due to its complex rupture mechanisms and associated hazards.To bridge gaps in understanding these intricate processes,this study emplo...Mining-related seismicity poses significant challenges in underground coal mining due to its complex rupture mechanisms and associated hazards.To bridge gaps in understanding these intricate processes,this study employed a multi-local seismic monitoring network,integrating both in-mine and local instruments at overlapping length scales.We specifically focused on a damaging local magnitude(ML)2.6 event and its aftershocks that occurred on 10 September 2022 in the vicinity of the 3308 working face of the Yangcheng coal mine in Shandong Province,China.Moment tensor(MT)inversion revealed a complex cascading rupture mechanism:an initial moment magnitude(M_(w))2.2 normal fault slip along the DF60 fault in an ESEeWNW direction,transitioning to a M_(w)3.0 event as the FD24 and DF60 faults unclamped.The scale-independent self-similarity and stress heterogeneity of mining-related seismicity were investigated through source parameter calculations,providing valuable insights into the driving mechanism of these seismic sequences.The in-mine network,constrained by its low dynamic changes,captured only the nucleation phase of the DF60 fault.Furthermore,standard decomposition of the MT solution from the seismic network proved inadequate for accurately identifying the complex nature of the rupture.To enhance safety and risk management in mining environments,we examined the implications of source reactivation within the cluster area post-stress-adjustment.This comprehensive multiscale analysis offers crucial insights into the complex rupture mechanisms and hazards associated with mining-related seismicity.The results underscore the importance of continuous multi-local network monitoring and advanced analytical techniques for improved disaster assessment and risk mitigation in mining operations.展开更多
The widespread variation of focal depths and fault plane solutions observed in the Hindukush region depicts continuous deformation along the Indian-Eurasian collision zone.For period of twelve years i.e.from 2010 to 2...The widespread variation of focal depths and fault plane solutions observed in the Hindukush region depicts continuous deformation along the Indian-Eurasian collision zone.For period of twelve years i.e.from 2010 to 2022,a total of 89 intermediate-depth earthquakes of magnitude(Mw)≥5.5 of the Hindukush Region were considered,relocated using both regional and tele seismic data with 90 per cent confidence limits of less than 20 km.Two distinct seismic activity clusters:First one at a deeper depth and second at a shallower depth having different P-axes were observed that verifies the internal structure and geometry of Hindukush zone as suggested in previous studies.Beneath the Hindukush collision zone,there exists a complex pattern of deformation,arising from a combination of compression,tension,shearing and necking states due to an unusual and a rare case of subduction that is not from oceanic plate.The Hindukush seismic zone extends from 70 to 300 km depth and mostly strikes east-west and then turns northeast.The relocated seismicity by merging data of seismic network close to Hindukush along with international data shows that the Hindukush zone may be divided vertically into upper and lower slabs separated by a gap at about 150 km depth at which strike and dip directions change sharply with significant structural changes.Seismicity rate is higher in the lower part of Hindukush,having large magnitude events in a small volume below 180 km forming complex pattern of source mechanisms.Contrary in upper part seismicity rate is lower and scattered.The Global CMT(Global Centroid-Moment-Tensor Project)source mechanisms of intermediate depth earthquakes have a systematic pattern of reverse faulting with the vertical T-axes,while shallow events do not have such pattern.The vertical T-axes of the intermediate-depth events may be attributed to negative buoyancy caused by subduction of the cold and denser slab.展开更多
With the advancement of fracturing technologies in deeper and more geologically complex formations,fault reactivation and induced seismicity have attracted increasing attention.The increasing frequency and magnitude o...With the advancement of fracturing technologies in deeper and more geologically complex formations,fault reactivation and induced seismicity have attracted increasing attention.The increasing frequency and magnitude of these events underscore the need for a robust understanding of the governing physical mechanisms.Elevated pore pressure,modified fault-loading conditions,and aseismic slip are widely acknowledged as the primary drivers.Recent studies have explored these mechanisms under varying factors,including fluid properties,rock ductility,poroelastic responses,and evolving fault stress states,thereby offering critical insights into model refinement.Probabilistic forecasting approaches,which combine statistical analyses of historical data with real-time monitoring,are being increasingly adopted in seismic risk assessments.In parallel,machine learning techniques are employed to process large seismic datasets and identify key patterns.However,their predictive capabilities remain limited by geological heterogeneity,subsurface complexity,and scarce observational data.Moreover,fluid–rock interactions further complicate the development of universally applicable models,thereby constraining the generalizability of mitigation strategies.This review synthesizes the current understanding of induced seismicity mechanisms,evaluates the prevailing prediction and mitigation methods,and identifies major challenges and future research directions.Advancements in these areas are essential to enhancing seismic risk management and ensuring the safe,sustainable development of deep-subsurface energy resources.展开更多
Accurate characte rization of the fault system is crucial for the exploration and development of fractu red reservoirs.The fault characterization technique based on multi-azimuth and multi-attribute fusion is a hotspo...Accurate characte rization of the fault system is crucial for the exploration and development of fractu red reservoirs.The fault characterization technique based on multi-azimuth and multi-attribute fusion is a hotspot.In this way,the fault structures of different scales can be identified and the characterization details of complex fault systems can be enriched by analyzing and fusing the fault-induced responses in multi-azimuth and multi-type seismic attributes.However,the current fusion methods are still in the stage of violent information stacking in utilizing fault information of multi-azimuth and multi-type seismic attributes,and the fault or fracture semantics in multi-type attributes are not fully considered and utilized.In this work,we propose a physic-guided multi-azimuth multi-type seismic attributes intelligent fusion method,which can mine fracture semantics from multi-azimuth seismic data and realize the effective fusion of fault-induced abnormal responses in multi-azimuth seismic coherence and curvature with the cooperation of the deep learning model and physical knowledge.The fused result can be used for multi-azimuth comprehensive characterization for multi-scale faults.The proposed method is successfully applied to an ultra-deep carbonate field survey.The results indicate the proposed method is superior to self-supervised-based,principal-component-analysis-based,and weighted-average-based fusion methods in fault characterization accuracy,and some medium-scale and microscale fault illusions in multi-azimuth seismic coherence and curvature can be removed in the fused result.展开更多
Ground response analysis and determination of site-specific ground motion parameters are necessary for evaluating seismic loads to enable sustainable design of aboveground and underground structures,particularly in de...Ground response analysis and determination of site-specific ground motion parameters are necessary for evaluating seismic loads to enable sustainable design of aboveground and underground structures,particularly in deep overburden sites.This study investigates the influence of bedrock interface conditions and depth of soil deposits on obtained site-specific ground motion parameters.Employing the one-dimensional seismic response analysis program SOILQUAKE,the ground responses of five representative soil profiles and 1050 case studies are calculated considering three different site models of seismic input interfaces.The analysis employs the actual bedrock interface with a shear wave velocity of 760 m/s as the reference input bedrock interface.The results illustrate that the selection of the bedrock interface condition significantly affects the seismic response on the ground surface of deep overburden sites.Specifically,the ground surface acceleration response spectra at longer periods are notably smaller compared to those at the actual bedrock site.This may present a challenge for designing long-period high-rise buildings situated in deep overburden sites.It is recommended to select a seismic input bedrock interface closely approximating the actual bedrock depth when conducting seismic response analyses for deep overburden sites.展开更多
基金supported by the National Natural Science Foundation of China(No.41941018)Shanghai Gaofeng Discipline Construction Funding.
文摘Strong seismic excitation and fault dislocation are likely to occur simultaneously in high-intensity seismic zones,causing severe damage to tunnels crossing active fault zones.This paper aims to develop a novel analytical solution to determine the longitudinal mechanical responses of tunnels subjected to the combined effects of seismic waves and strike-slip faulting.Adopting the elastic springbeam model,the seismic waves are modelled as shear horizontal(SH)waves and the fault dislocation follows an S-shaped pattern;the superposition principle for free-fielddisplacements caused by both effects is assumed.In addition,the transmission and reflectionof seismic waves at the fault-rock geological interface and the tangential contact conditions at the tunnel-rock interface are considered.The analytical model is validated against numerical simulations,confirmingits accuracy in calculating tunnel responses.Moreover,a parametric study is conducted to evaluate the impact of key factors,including fault displacement,fault zone width,fault dip angle,earthquake frequency,rock conditions,tunnel lining stiffness,and tangential contact conditions,on tunnel responses.Compared with each effect alone,the combined effects of seismic waves and strike-slip faulting significantlychange the tunnel deformation and internal forces,leading to increased tunnel responses,especially within the fault zone and near the fault-rock interfaces.Depending on specificparameters,tunnel responses can be classifiedinto seismic-dominated,faulting-dominated,and seismic-faulting coupled responses on the basis of the relative contributions of each effect.The proposed analytical solution can be applied to quickly predict the longitudinal mechanical behaviour of tunnels under such combined effects in engineering applications.
基金supported by the Earthquake Science and Technology Spark Plan Project(No.XH23041C)The Natural Science Foundation of Gansu Province(No.22JR11RA090)Gansu Lanzhou Geophysics National Observation and Research Station(No.2021Y14).
文摘The widely distributed loess deposits in the Yellow River Basin exhibit unique engineering geological characteristics.The variations in their thickness and stratigraphic structure significantly amplify ground motion parameters,directly influencing the regional seismic hazard risk level.This study methodically conducted on-site studies and observations of building collapses and damages resulting from seismic amplification effects,using the Wenchuan M_(S)8.0 earthquake as a case study.Comprehensive experimental and numerical simulation studies were carried out.A large-scale shaking table test was performed,and numerical models for 14 different loess sites types were established.Various types of seismic waves were incorporated into these models for systematic numerical simulation calculations.The research reveals the mechanisms by which loess deposit thickness and stratigraphic structure in the Yellow River Basin affect seismic ground motion amplification.The results indicate that as the epicentral distance increases,the peak ground motion shows a marked attenuation trend,with the horizontal component attenuating substantially faster than the vertical component.As the overlying loess layer thickness increases from 50 to 100 m,the seismic intensity may escalate by 3−4 degrees,and the peak acceleration may amplify by 1.5−2.2 times.With the augmentation of loess deposit thickness and the proliferation of soil layers,both the peak acceleration response spectrum and the characteristic period demonstrate an upward tendency,exhibiting slight fluctuations contingent upon the seismic wave type.
基金supported by the National Key R&D Program of China(No.2022YFC3003502).
文摘The northern segment of the North-South Seismic Belt is characterized by intense crustal deformation,well-developed active tectonics,and frequent occurrences of strong earthquakes.Therefore,conducting a Probabilistic Seismic Hazard Analysis(PSHA)for this region is of significant importance for supporting seismic fortification in major engineering projects and formulating disaster prevention and mitigation policies.In this study,a composite seismic source model was constructed by integrating data on historical earthquakes,active faults,and paleoseismicity.Furthermore,a logic tree framework was employed to quantify epistemic uncertainties,enabling a systematic seismic hazard assessment of the region.To more accurately characterize the spatial heterogeneity of seismic activity,improvements were made to both the Circular Spatial Smoothing Model(CSSM)with a fixed radius and the Adaptive Spatial Smoothing Model(ASSM),with full consideration given to the spatiotemporal completeness of historical earthquake magnitudes.Regarding the CSSM,for scenarios involving small sample sizes in earthquake catalogs,the cross-validation method proposed in this study demonstrated higher robustness than the maximum likelihood method in determining the optimal correlation distance.Performance evaluation results indicate that while both models effectively characterize seismic activity,the ASSM exhibits superior overall predictive performance compared to the CSSM,owing to its ability to adaptively adjust the smoothing radius according to seismic density.Significant discrepancies were observed in the Peak Ground Acceleration(PGA)results calculated with a 10%probability of exceedance in 50 years across different combinations of seismic source models.The single spatially smoothed point-source model yielded a maximum PGA of approximately 0.52 g,with high-value areas concentrated near historical epicenters,thereby significantly underestimating the hazard associated with major fault zones.When combined with the simple fault-source model,the maximum PGA increased to 0.8 g,with high-value zones exhibiting a zonal distribution along faults;however,the risk remained underestimated for faults with low slip rates that are nevertheless approaching their recurrence cycles.Following the introduction of the time-dependent characteristic fault-source model,local PGA values for faults in the middle-to-late stages of their recurrence cycles increased by a factor of 2 to 7 compared to the single model.These results demonstrate that the characteristic fault-source model reasonably delineates the time-dependence of large earthquake recurrence,thereby providing a more accurate assessment of imminent seismic risks.By comprehensively applying the improved spatially smoothed pointsource model,the simple fault-source model,and the characteristic fault-source model,the following faults within the region were identified as having high seismic hazard:the Huangxianggou,Zhangxian,and Tianshui segments of the Xiqinling northern edge fault;the Maqin-Maqu segment of the Dongkunlun fault;the Longriqu fault;the Maoergai fault;the Elashan fault;the Riyueshan fault;the eastern segment of the Lenglongling fault;the Maxianshan segment of the Maxianshan northern Margin fault;and the Maomaoshan-Jinqianghe segment of the Laohushan-Maomaoshan fault.As these faults are located within seismic gaps or are approaching the recurrence periods of large earthquakes,they should be prioritized for current and future seismic monitoring as well as disaster prevention and mitigation efforts.
基金supported by the National Natural Science Foundation of China(Grant No.11974044)。
文摘Topological insulators with localized edge or interface states have been extensively studied,particularly in phononic crystals and related fields;however,their application in seismic metamaterials remains largely unexplored.To address this gap,we designed a topological seismic metamaterial,where the topological interface is formed by joining the ends of two distinct one-dimensional periodic lattices.The first full-scale field experiment confirms the existence of topological interface states,which exhibit pronounced localization characteristics and induce a resonant amplification effect of 7.2 dB on the total energy of seismic surface waves.This study provides the first experimental validation for the implementation of topological principles in the design of seismic metamaterials,enabling novel approaches to high-sensitivity seismic detection and efficient energy localization for wave control.
基金supported by the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology under Grant 2024yjrc64the National Key R&D Program of China under Grant 2018YFC1504102。
文摘The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velocity models are available for the LMS fault zone,high-resolution velocity models are lacking.Therefore,a dense array of 240 short-period seismometers was deployed around the central segment of the LMS fault zone for approximately 30 days to monitor earthquakes and characterize fine structures of the fault zone.Considering the large quantity of observed seismic data,the data processing workflow consisted of deep learning-based automatic earthquake detection,phase arrival picking,and association.Compared with the earthquake catalog released by the China Earthquake Administration,many more earthquakes were detected by the dense array.Double-difference seismic tomography was adopted to determine V_(p),V_(s),and V_(p)/V_(s)models as well as earthquake locations.The checkerboard test showed that the velocity models have spatial resolutions of approximately 5 km in the horizontal directions and 2 km at depth.To the west of the Yingxiu–Beichuan Fault(YBF),the Precambrian Pengguan complex,where most of earthquakes occurred,is characterized by high velocity and low V_(p)/V_(s)values.In comparison,to the east of the YBF,the Upper Paleozoic to Jurassic sediments,where few earthquakes occurred,show low velocity and high V_(p)/V_(s)values.Our results suggest that the earthquake activity in the LMS fault zone is controlled by the strength of the rock compositions.When the high-resolution velocity models were combined with the relocated earthquakes,we were also able to delineate the fault geometry for different faults in the LMS fault zone.
基金Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No.2023C01National Natural Science Foundation of China under Grant No.52478570Distinguished Young Scholars Program of the Natural Science Foundation of Heilongjiang Province,China under Grant No.JQ2024E002。
文摘The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismic wave simulation with a data-constrained finite-fault rupture model.The constraint is implemented by identifying the optimal ground motion models(GMMs)through a scoring system that selects the best-fit GMMs to mid-and far-field China Earthquake Networks Center(CENC)seismic network data;and applying the optimal GMMs to refine the rupture model parameters for near-fault intensity field simulation.The simulated near-fault seismic intensity field reproduces seismic intensities collected from Myanmar’s sparse seismic network and concentrated in≥Ⅷintensity zones within 50 km of the projected fault plane;and identifies abnormal intensity regions exhibiting≥Ⅹintensity along the Meiktila-Naypyidaw corridor and near Shwebo that are attributed to soft soil amplification effects and near-fault directivity.This framework can also be applied to post-earthquake assessments in other similar regions.
基金funded by the Russian Federation Ministry of Science and Higher Education according to the research project FUUW-2025-0033 of the FECIAR UrB RASthe Arkhangelsk branch of the Territorial Fund of Geological Information for the North-Western Federal District for providing geological and geophysical reports for the objects under consideration。
文摘Studying the inner structure of intraplate earthquakes originating in aseismic areas,which are poorly covered by seismic networks or as historical earthquakes is usually the only way to get knowledge about their source mechanisms,which is partially essential for a deeper understanding of intraplate geodynamics.The epicentral zones of earthquakes are situated in hard-toreach areas,so,using active seismic methods for such purposes is unreasonable or even impossible because of high cost and logistical difficulties.We propose a novel approach that combines diverse passive seismic methods,which allows us to get sufficient information about geological environment structure for such task solutions.As an example,we investigated the inner structure of platform earthquake epicentral zone originated up north of Russia.We used four passive seismic methods:microseismic sounding method,passive seismic interferometry,HVSR method,and microseismic activity method.We show that passive seismic data,recoded in the same installation and processed by these different methods,can provide sufficient information about structure of studied environment,needed to explain source mechanism.In sum,the hypocenter zone is presented by intersection of vertical faults and a lateral fractured zone in the middle crust.Results were confirmed by comparison with results by active seismic methods.
基金The work was carried out in the framework of earmarked funding“Assessment of seismic hazard of territories of Kazakhstan on modern scientific and methodological basis”,programme code number F.0980.Source of funding-Ministry of Science and Higher Education of the Republic of Kazakhstan.
文摘This article aims to enhance seismic hazard assessment methods for Kazakhstan’s seismotectonic conditions.It combines probabilistic seismic hazard analysis(PSHA),ground motion simulation,sitespecific geological and geotechnical data analysis,and seismic scenario analysis to develop Probabilistic General Seismic Zoning(GSZ)maps for Kazakhstan and Probabilistic Seismic Microzoning maps for Almaty.These maps align with Eurocode 8 principles,incorporating seismic intensity and engineering parameters like peak ground acceleration(PGA).The new procedure,applied in national projects,has resulted in GSZ maps for the country,seismic microzoning maps for Almaty,and detailed seismic zoning maps for East Kazakhstan.These maps,part of a regulatory document,guide earthquake-resistant design and construction.They offer a comprehensive assessment of seismic hazards,integrating traditional Medvedev-Sponheuer-Karnik(MSK-64)intensity scale points with quantitative parameters like peak ground acceleration.This innovative approach promises to advance methods for quantifying seismic hazards in specific regions.
基金This study is part of the project 2021RUAPON-REACT EU-Finanziamento PON“Ricerca e Innovazione”20142020,grant n.19-G-12543-2,funded by the Italian Ministry of University and Research(MUR)This study was developed in the frame of“The Geosciences for Sustainable Development”project(Budget Minis-tero dell'Universita e della Ricerca-Dipartimenti di Eccellenza 2023-2027,code n.C93C23002690001).
文摘A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the studied area,and(ii)to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate.Despite the limited number of stations and the short-period char-acter of the seismic sensors,the empirical Green's functions were retrieved for all station pairs using two years of passive data.Both group and phase velocities were derived,the former using the widespread frequency-time analysis,the latter through the analysis of the real part of the cross-spectra.The main advantage of combining these two methods is a more accurate identification of higher modes,resulting in a reduction of ambiguity during picking and data interpretation.Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths.The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method,for this specific case.The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast.Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs.While the average model can be used for a more accurate location of the local natural seismicity,the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth.Finally,the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field.Errors as large as 4%have been found,related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources,and to the very short time widows.It was shown that using a one-year time window these errors arereduced to 0.3%.
基金supported by National Natural Science Foundation of China(No.41861134009)the PIA ANID(AFB180004)the ANID(PCI PII-180003).
文摘The Chilean subduction zone is one of the most seismically active regions globally,characterized by extensive intermediatedepth seismicity in the slab.In this study,we construct a new earthquake catalog for northern Chile using seismic waveforms assembled for the period of 2014-2019,from which 320,070 P-wave and 232,907 S-wave first arrivals are obtained for 25,763 earthquakes.Grid search location method NonLinLoc is applied to determine initial earthquake locations and double-difference location method is used to improve relative event locations.The distribution of earthquakes exhibits distinct patterns to the north and south of 21°S.There are many more earthquakes deeper than~150 km to the south of 21°S,while relatively fewer to the north.The intraslab earthquakes shallower than~80 km generally reveal a distinct double seismic zone,and the gap between the two seismic planes disappears at a depth of approximately~80 km,followed by a concentration of seismicity in the depth range of~80-150 km.In the deeper slab,there exist several seismicity clusters with distinct earthquake activities down to~300 km.These characteristics shown in slab seismicity are likely caused by different mechanisms and can be helpful for understanding the subduction process.
基金Institute of Seismology,Academy of Science of the Republic of Uzbekistan Grant number is FA-F8-007
文摘The quantitative assessment of seismic hazard of Uzbekistan has been examined,and new maps of seismic zoning has been developed.Quasihomogeneous seismological areas and seism-generating zones are considered as seismic sources,based on the analysis of seismotectonic data.The seismological parametrization of seismic sources has been carried out,including determination of parameters of earthquake’s reoccurrence for different power levels,seismic potential and a predominant type motion in the source of occurring earthquakes.The reoccurrence parameters of seismic sources were both determined by directly from Gutenberg-Richter dependence law and using summation and distribution methods.Setting were conducted separately for a sample of strong(M≥5)and weak(M<5)earthquakes.The seismic potential Mmax of seismic active zones was determined by seismological and seismotectonic methods.The predominant type of movement in the faulting for each seismic source is reverse fault.Regional regularities of seismic intensity attenuation with distance of different power levels are determined.The seismic hazard is expressed in points of a macro seismic scale,in velocities and accelerations of ground motions,and it is characterized by calculated seismic intensity with the set probability(P=0.9,P=0.95,P=0.98 and P=0.99)not exceed within 50 years in the constructed maps.In seismic zoning map,the seismic intensity was in average soil conditions.The developed complex of maps considers a number of uncertainties of the input parameters in relation to both incompleteness of initial seismological and seismotectonic data,as well the probabilistic nature of seismic process and seismic intensity.The factors of uncertainty of the input parameters are taken into account by constructing a logic tree.The constructed maps of the general seismic zoning are intended for imple menting antiseismic actions in Uzbekistan.
基金supported by the National Key R&D Program of China(No.2023YFC3805100)the National Natural Science Foundation of China(Nos.52222811 and 52494963)。
文摘Earthquakes pose significant perils to the built environment in urban areas.To avert the calamitous aftermath of earthquakes,it is imperative to construct seismic resilient cities.Due to the intricacy of the concept of urban seismic resilience(USR),its assessment is a large-scale system engineering issue.The assessment of USR should be based on the notion of urban seismic capacity(USC)assessment,which includes casualties,economic loss,and recovery time as criteria.Functionality loss is also included in the assessment of USR in addition to these criteria.The assessment indicator system comprising five dimensions(building and lifeline infrastructure,environment,society,economy,and institution)and 20 indicators has been devised to quantify USR.The analytical hierarchy process(AHP)is utilized to compute the weights of the criteria,dimensions,and indicators in the urban seismic resilience assessment(USRA)indicator system.When the necessary data for a city are obtainable,the seismic resilience of that city can be assessed using this framework.To illustrate the proposed methodology,a moderate-sized city in China was selected as a case study.The assessment results indicate a high level of USR,suggesting that the city possesses strong capabilities to withstand and recover from potential future earthquakes.
基金Derek Elsworth acknowledges the support from a Gledden Visiting Fellowship from the Institute of Advanced Studies at the University of Western Australia,Australia,and the G.Albert Shoemaker Endowment at Pennsylvania State University,USA.
文摘Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mitigation. Key controls on seismicity are defined in terms of fault and fracture strength, second-order frictional response and stability, and competing fluid-driven mechanisms for arrest. We desire to constrain maximum event magnitudes in triggered earthquakes by relating pre-existing critical stresses to fluid injection volume to explain why some recorded events are significantly larger than anticipated seismic moment thresholds. This formalism is consistent with several uncharacteristically large fluid injection-triggered earthquakes. Such methods of reactivating fractures and faults by hydraulic stimulation in shear or tensile fracturing are routinely used to create permeability in the subsurface. Microearthquakes (MEQs) generated by such stimulations can be used to diagnose permeability evolution. Although high-fidelity data sets are scarce, the EGS-Collab and Utah FORGE hydraulic stimulation field demonstration projects provide high-fidelity data sets that concurrently track permeability evolution and triggered seismicity. Machine learning deciphers the principal features of MEQs and the resulting permeability evolution that best track permeability changes – with transfer learning methods allowing robust predictions across multiple eological settings. Changes in permeability at reactivated fractures in both shear and extensional modes suggest that permeability change (Δk) scales with the seismic moment (M) of individual MEQs as Δk∝M. This scaling relation is exact at early times but degrades with successive MEQs, but provides a method for characterizing crustal permeability evolution using MEQs, alone. Importantly, we quantify for the first time the role of prestress in defining the elevated magnitude and seismic moment of fluid injection-triggered events, and demonstrate that such MEQs can also be used as diagnostic in quantifying permeability evolution in the crust.
基金funded by the Joint Funds of the National Natural Science Foundation of China(U23A20671)the Major Project of Inner Mongolia Science and Technology(2021ZD0034)the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering(No.Z021003)。
文摘The occurrence time and magnitude of injection-induced seismicity are influenced by engineering factors,such as wellhead pressure,injection location,injection volume,and injection rate.Understanding the relationship between injection operations and seismic magnitude is of great significance for optimizing industrial production and reducing earthquake disasters.Numerical simulation of hydromechanical coupling is a crucial method for studying injection-induced seismicity.However,few studies have explored the risk management measures for injection-induced seismicity from the perspective of engineering.How seismic magnitudes can be reduced through reasonable adjustments to injection operations in engineering remains unclear.Therefore,in this study,a 3D hydro-mechanical coupling model involving multiple faults and injection wells was established based on the geological background and well location of Fox Creek,Canada.Different injection schemes under multi-well and multi-fault conditions were studied,and a traffic light system was used to simulate and control the magnitudes under a multi-well injection scheme.Specifically,we simulated injection scenarios involving up to three wells and analyzed the response of five faults.We compared the maximum moment magnitude of different scenarios by controlling the same injection volume.The results revealed the effect and advantage of the multi-well scheme in reducing seismic magnitude.To reduce the risk of induced seismicity,utilizing far-fault operational wells to compensate for the effects of near-fault operational wells proves to be an efficient and cost-effective method,with potential for wide practical applications.
文摘A seismic swarm occurred southeast of King George Island,South Shetland Islands,Antarctica,between August 2020 and February 2021.This work intends to parameterize seismic events recorded by seismic station AM.R4DE2 from 15 September to 31 October 2020.Using the localization methodology with a single station,the record of the entire period was analyzed manually to determine the local magnitude,hypocentral distance,epicentral distance,backazimuth,and location of the epicenter for each event.We could parameterize 6362 events,although we estimate the occurrence to be around 20000 for the period.The results suggest a magmatic origin for the swarm,supporting previous studies.Seismicity exhibited a southeastward migration away from King George Island,as indicated by a progressive increase in epicentral distance over time.Most events were classified as volcanic and volcano-tectonic,supporting a magmatogenesis hypothesis linked to the opening of Bransfield Ridge.
基金funded by the National Natural Science Foundation of China(Grant No.51574225)Shandong Energy Group(Grant No.SNKJ2022BJ03-R28)for Caiping Lu+1 种基金the Research Team on MonitoringActivity Mechanisms of Unnatural Earthquakes of Shandong Earthquake Agency(Grant No.TD202301)for Chengyu Liu.
文摘Mining-related seismicity poses significant challenges in underground coal mining due to its complex rupture mechanisms and associated hazards.To bridge gaps in understanding these intricate processes,this study employed a multi-local seismic monitoring network,integrating both in-mine and local instruments at overlapping length scales.We specifically focused on a damaging local magnitude(ML)2.6 event and its aftershocks that occurred on 10 September 2022 in the vicinity of the 3308 working face of the Yangcheng coal mine in Shandong Province,China.Moment tensor(MT)inversion revealed a complex cascading rupture mechanism:an initial moment magnitude(M_(w))2.2 normal fault slip along the DF60 fault in an ESEeWNW direction,transitioning to a M_(w)3.0 event as the FD24 and DF60 faults unclamped.The scale-independent self-similarity and stress heterogeneity of mining-related seismicity were investigated through source parameter calculations,providing valuable insights into the driving mechanism of these seismic sequences.The in-mine network,constrained by its low dynamic changes,captured only the nucleation phase of the DF60 fault.Furthermore,standard decomposition of the MT solution from the seismic network proved inadequate for accurately identifying the complex nature of the rupture.To enhance safety and risk management in mining environments,we examined the implications of source reactivation within the cluster area post-stress-adjustment.This comprehensive multiscale analysis offers crucial insights into the complex rupture mechanisms and hazards associated with mining-related seismicity.The results underscore the importance of continuous multi-local network monitoring and advanced analytical techniques for improved disaster assessment and risk mitigation in mining operations.
文摘The widespread variation of focal depths and fault plane solutions observed in the Hindukush region depicts continuous deformation along the Indian-Eurasian collision zone.For period of twelve years i.e.from 2010 to 2022,a total of 89 intermediate-depth earthquakes of magnitude(Mw)≥5.5 of the Hindukush Region were considered,relocated using both regional and tele seismic data with 90 per cent confidence limits of less than 20 km.Two distinct seismic activity clusters:First one at a deeper depth and second at a shallower depth having different P-axes were observed that verifies the internal structure and geometry of Hindukush zone as suggested in previous studies.Beneath the Hindukush collision zone,there exists a complex pattern of deformation,arising from a combination of compression,tension,shearing and necking states due to an unusual and a rare case of subduction that is not from oceanic plate.The Hindukush seismic zone extends from 70 to 300 km depth and mostly strikes east-west and then turns northeast.The relocated seismicity by merging data of seismic network close to Hindukush along with international data shows that the Hindukush zone may be divided vertically into upper and lower slabs separated by a gap at about 150 km depth at which strike and dip directions change sharply with significant structural changes.Seismicity rate is higher in the lower part of Hindukush,having large magnitude events in a small volume below 180 km forming complex pattern of source mechanisms.Contrary in upper part seismicity rate is lower and scattered.The Global CMT(Global Centroid-Moment-Tensor Project)source mechanisms of intermediate depth earthquakes have a systematic pattern of reverse faulting with the vertical T-axes,while shallow events do not have such pattern.The vertical T-axes of the intermediate-depth events may be attributed to negative buoyancy caused by subduction of the cold and denser slab.
基金supported in part by the National Key Research and Development Project of China(No.2022YFC3004602)in part by the National Natural Science Foundation of China(Nos.52121003 and 52442406)。
文摘With the advancement of fracturing technologies in deeper and more geologically complex formations,fault reactivation and induced seismicity have attracted increasing attention.The increasing frequency and magnitude of these events underscore the need for a robust understanding of the governing physical mechanisms.Elevated pore pressure,modified fault-loading conditions,and aseismic slip are widely acknowledged as the primary drivers.Recent studies have explored these mechanisms under varying factors,including fluid properties,rock ductility,poroelastic responses,and evolving fault stress states,thereby offering critical insights into model refinement.Probabilistic forecasting approaches,which combine statistical analyses of historical data with real-time monitoring,are being increasingly adopted in seismic risk assessments.In parallel,machine learning techniques are employed to process large seismic datasets and identify key patterns.However,their predictive capabilities remain limited by geological heterogeneity,subsurface complexity,and scarce observational data.Moreover,fluid–rock interactions further complicate the development of universally applicable models,thereby constraining the generalizability of mitigation strategies.This review synthesizes the current understanding of induced seismicity mechanisms,evaluates the prevailing prediction and mitigation methods,and identifies major challenges and future research directions.Advancements in these areas are essential to enhancing seismic risk management and ensuring the safe,sustainable development of deep-subsurface energy resources.
基金sponsorship of the National Natural Science Foundation of China (42430809, 42030103)the Basic Research Funds for Northeast Petroleum University in Heilongjiang Province (2025GPL-01)。
文摘Accurate characte rization of the fault system is crucial for the exploration and development of fractu red reservoirs.The fault characterization technique based on multi-azimuth and multi-attribute fusion is a hotspot.In this way,the fault structures of different scales can be identified and the characterization details of complex fault systems can be enriched by analyzing and fusing the fault-induced responses in multi-azimuth and multi-type seismic attributes.However,the current fusion methods are still in the stage of violent information stacking in utilizing fault information of multi-azimuth and multi-type seismic attributes,and the fault or fracture semantics in multi-type attributes are not fully considered and utilized.In this work,we propose a physic-guided multi-azimuth multi-type seismic attributes intelligent fusion method,which can mine fracture semantics from multi-azimuth seismic data and realize the effective fusion of fault-induced abnormal responses in multi-azimuth seismic coherence and curvature with the cooperation of the deep learning model and physical knowledge.The fused result can be used for multi-azimuth comprehensive characterization for multi-scale faults.The proposed method is successfully applied to an ultra-deep carbonate field survey.The results indicate the proposed method is superior to self-supervised-based,principal-component-analysis-based,and weighted-average-based fusion methods in fault characterization accuracy,and some medium-scale and microscale fault illusions in multi-azimuth seismic coherence and curvature can be removed in the fused result.
基金supported by the National Natural Science Foundation of China(Nos.52408435,52278384)。
文摘Ground response analysis and determination of site-specific ground motion parameters are necessary for evaluating seismic loads to enable sustainable design of aboveground and underground structures,particularly in deep overburden sites.This study investigates the influence of bedrock interface conditions and depth of soil deposits on obtained site-specific ground motion parameters.Employing the one-dimensional seismic response analysis program SOILQUAKE,the ground responses of five representative soil profiles and 1050 case studies are calculated considering three different site models of seismic input interfaces.The analysis employs the actual bedrock interface with a shear wave velocity of 760 m/s as the reference input bedrock interface.The results illustrate that the selection of the bedrock interface condition significantly affects the seismic response on the ground surface of deep overburden sites.Specifically,the ground surface acceleration response spectra at longer periods are notably smaller compared to those at the actual bedrock site.This may present a challenge for designing long-period high-rise buildings situated in deep overburden sites.It is recommended to select a seismic input bedrock interface closely approximating the actual bedrock depth when conducting seismic response analyses for deep overburden sites.
