A new numerical model is presented to simulate fracture initiation and propagation in geological structures. This model is based on the recent amalgamation of established failure and fracture mechanics theory, which h...A new numerical model is presented to simulate fracture initiation and propagation in geological structures. This model is based on the recent amalgamation of established failure and fracture mechanics theory, which has been implemented to the finite difference FLAC code as a constitutive FISH userdefined-model. Validation of the model has been studied on the basis of comparing the transitional failure modes in rock. It is shown that the model is capable of accurately simulating fracture distributions over entire brittle to ductile rock phases. The application of the model during longwall retreat simulation highlighted several caving characteristics relevant to varying geological condition. The distribution and behaviour of modelled fractures were both realistic and shown to provide an enhanced post failure analysis to geological structures in FLAC. Moreover, the model introduces new potential insight towards the failure analysis of more complicated problems. This is best suited towards improving safety and efficiency in mines through the prediction of various key fractures and caving characteristics of geological structures.展开更多
The relatively high stress probably leads to generation of a fractured or even instable area around a working coalface. Also, the generated weak area often evolves into an easy-infiltrating field of water/gas to great...The relatively high stress probably leads to generation of a fractured or even instable area around a working coalface. Also, the generated weak area often evolves into an easy-infiltrating field of water/gas to greatly increase probability of accident occurrence. To reveal the distribution of high stress around working faces, we put forward the mode-I-crack compression model. In this model, the goaf following a working face is regarded as a mode-I crack in an infinite plate, and the self-gravity of overlaying strata is transformed into an uniform pressure applied normal to the upper edge of the model crack. Solving this problem is based on the Westergaard complex stress function. For comparison, the software RFPA-2D is also employed to simulate the same mining problem, and furthermore extendedly to calculate the stress interference induced by the simultaneous advances of two different working faces. The results show that, the area close to a working face or the goaf tail has the maximum stress, and the stress is distributed directly proportional to the square root of the advance and inversely proportional to the square root of the distance to the working face. The simultaneous advances of two neighboring working faces in different horizontals can lead to extremely high resultant stress in an interference area.展开更多
Shear wave splitting is studied based on the digital waveforms of three seismic stations DJS, SZD and WUJ, which were set up after the Jiujiang-Ruichang MS5.7 earthquake of November 26, 2005 around the epicenter area....Shear wave splitting is studied based on the digital waveforms of three seismic stations DJS, SZD and WUJ, which were set up after the Jiujiang-Ruichang MS5.7 earthquake of November 26, 2005 around the epicenter area. The result shows that the time delays of slow shear waves of the DJS station, which is not far from the epicenter and where the distribution of faults is complex near the station, are relatively larger and the polarization directions of fast shear waves are not concentrated; the predominant polarization direction of fast shear waves of WUJ station, with single fault distributed nearby, has a difference of 35° to the strike of the fault and is inconsistent with the direction of regional principal compressive stress. The predominant polarization direction of fast shear waves of SZD station with no faults nearby is in accordance with regional principal compressive stress. There is no obvious regular relationshipship between the delay time and the focal depth.展开更多
Shallow fissures, being the main infiltration paths of fluid on the surface of a slope, played an important role in the whole process of a landslide. However, the spatial distribution characteristics of fissures in th...Shallow fissures, being the main infiltration paths of fluid on the surface of a slope, played an important role in the whole process of a landslide. However, the spatial distribution characteristics of fissures in the slope are difficult to be determined. In this study, we attempted to characterize the variation pattern of slope fissures along depth in the Wenchuan earthquake area in Sichuan Province by combining engineering geological investigation, geomorphologic analysis and geophysical investigation. The geophysical methods that were used in this study include Multichannel Analysis of Surface Wave(MASW), Ground Penetrating Radar(GPR) and Electrical Resistivity Tomography(ERT). The results suggested that geophysical parameters(shear wave velocity, electromagnetic signals attenuation and resistivity) could provide valuable information for the spatial network of shallow fissures. Through the verification by engineering geological survey and geophysical sensitivity analysis, this work highlighted that MASW was the most appropriate technique to delineate the propagation of shallow fissures in a gravel soil slope.展开更多
Based on seismic data from the regional network for the last 34 years, we analyzed the present fault behavior of major fault zones around the Mabian area, southern Sichuan, and identified the risky fault-segments for ...Based on seismic data from the regional network for the last 34 years, we analyzed the present fault behavior of major fault zones around the Mabian area, southern Sichuan, and identified the risky fault-segments for potential future. The method of analysis is a combination of activity background of historical strong earthquakes mainly show ~ ( 1 ) The spatial distribution of b-values strong and large earthquakes in the spatial distribution of b-values with and current seismicity. Our results indicates significant heterogeneity in the studied area, which reflects the spatial difference of cumulative stress levels along various fault zones and segments. (2) Three anomalously low b-value areas with different dimensions were identified along the Mabian-Yanjin fault zone. These anomalies can be asperities under relatively high cumulated stress levels. Two asperities are located in the north of Mabian county, in Lidian town in western Muchuan county, and near Yanjin at the south end of the fault zone. These two areas represent potential large earthquake seismogenic sites around the Mabian area in the near future. Besides them, the third relatively smaller asperity is identified at southern Suijiang, as another potential strong- earthquake source. (3) An asperity along the southwestern segment of the Longquanshan fault zone indicates the site of potential moderate-to-strong earthquakes. (4) The asperity along the segment between Huangmu town in Hanyuan county and Longchi town in Emeishan city on Jinkouhe-Meigu fault has potential for a moderate-strong earthquake.展开更多
The influences of subsurface cracks,distributing along the axial direction,on the rolling contact fatigue(RCF)faliure in a bearing ring are investigated.A realistic three-dimensional model of the bearing ring containi...The influences of subsurface cracks,distributing along the axial direction,on the rolling contact fatigue(RCF)faliure in a bearing ring are investigated.A realistic three-dimensional model of the bearing ring containing three subsurface cracks is used to evaluate the fatigue crack propagation based on stress intensity factor(SIF)calculations.Moreover,the distributions of the subsurface cracks along the axial direction are varied to study their effects on RCF.The results provide valuable guidelines for enhanced understanding of RCF in bearings.展开更多
On 25 April 2015,an M_w 7.8 earthquake occurred on the Main Himalaya Thrust fault with a dip angle of^7° about77 km northwest of Kathmandu,Nepal.This Nepal Gorkha event is the largest one on the Himalayan thrust ...On 25 April 2015,an M_w 7.8 earthquake occurred on the Main Himalaya Thrust fault with a dip angle of^7° about77 km northwest of Kathmandu,Nepal.This Nepal Gorkha event is the largest one on the Himalayan thrust belt since 1950.Here we use the compressive sensing method in the frequency domain to track the seismic radiation and rupture process of this event using teleseismic P waves recorded by array stations in North America.We also compute the distribution of static shear stress changes on the fault plane from a coseismic slip model.Our results indicate a dominant east-southeastward unilateral rupture process from the epicenter with an average rupture speed of ~3 km s^(-1).Coseismic radiation of this earthquake shows clear frequency-dependent features.The lower frequency(0.05-0.3 Hz) radiation mainly originates from large coseismic slip regions with negative coseismic shear stress changes.In comparison,higher frequency(0.3-0.6 Hz) radiation appears to be from the down-dip part around the margin of large slip areas,which has been loaded and presents positive coseismic shear stress changes.We propose an asperity model to interpret this Nepal earthquake sequence and compare the frequency-dependent coseismic radiation with that in subduction zones.Such frequency-dependent radiation indicates the depth-varying frictional properties on the plate interface of the Nepal section in the main Himalaya thrust system,similar to previous findings in oceanic subduction zones.Our findings provide further evidence of the spatial correlation between changes of static stress status on the fault plane and the observed frequency-dependent coseismic radiation during large earthquakes.Our results show that the frequency-dependent coseismic radiation is not only found for megathrust earthquakes in the oceanic subduction environment,but also holds true for thrust events in the continental collision zone.展开更多
文摘A new numerical model is presented to simulate fracture initiation and propagation in geological structures. This model is based on the recent amalgamation of established failure and fracture mechanics theory, which has been implemented to the finite difference FLAC code as a constitutive FISH userdefined-model. Validation of the model has been studied on the basis of comparing the transitional failure modes in rock. It is shown that the model is capable of accurately simulating fracture distributions over entire brittle to ductile rock phases. The application of the model during longwall retreat simulation highlighted several caving characteristics relevant to varying geological condition. The distribution and behaviour of modelled fractures were both realistic and shown to provide an enhanced post failure analysis to geological structures in FLAC. Moreover, the model introduces new potential insight towards the failure analysis of more complicated problems. This is best suited towards improving safety and efficiency in mines through the prediction of various key fractures and caving characteristics of geological structures.
基金Projects 50774083 and 40811120546 supported by the National Natural Science Foundation of ChinaNCET-07-0803 by the Program for New Century Ex-cellent Talents in University 2005CB221502 by the National Basic Research Program of China
文摘The relatively high stress probably leads to generation of a fractured or even instable area around a working coalface. Also, the generated weak area often evolves into an easy-infiltrating field of water/gas to greatly increase probability of accident occurrence. To reveal the distribution of high stress around working faces, we put forward the mode-I-crack compression model. In this model, the goaf following a working face is regarded as a mode-I crack in an infinite plate, and the self-gravity of overlaying strata is transformed into an uniform pressure applied normal to the upper edge of the model crack. Solving this problem is based on the Westergaard complex stress function. For comparison, the software RFPA-2D is also employed to simulate the same mining problem, and furthermore extendedly to calculate the stress interference induced by the simultaneous advances of two different working faces. The results show that, the area close to a working face or the goaf tail has the maximum stress, and the stress is distributed directly proportional to the square root of the advance and inversely proportional to the square root of the distance to the working face. The simultaneous advances of two neighboring working faces in different horizontals can lead to extremely high resultant stress in an interference area.
基金sponsored by the China Spark Program of Earthquake Science and Technology(XH12027)the Three-Combination Topics of China Earthquake Administration of"Research on the Crustal Medium Anisotropy in the Jiujiang-Ruichang Earthquake Area"the Special Fund of Seismic Industry Research(201008007)
文摘Shear wave splitting is studied based on the digital waveforms of three seismic stations DJS, SZD and WUJ, which were set up after the Jiujiang-Ruichang MS5.7 earthquake of November 26, 2005 around the epicenter area. The result shows that the time delays of slow shear waves of the DJS station, which is not far from the epicenter and where the distribution of faults is complex near the station, are relatively larger and the polarization directions of fast shear waves are not concentrated; the predominant polarization direction of fast shear waves of WUJ station, with single fault distributed nearby, has a difference of 35° to the strike of the fault and is inconsistent with the direction of regional principal compressive stress. The predominant polarization direction of fast shear waves of SZD station with no faults nearby is in accordance with regional principal compressive stress. There is no obvious regular relationshipship between the delay time and the focal depth.
基金financially supported by the National Basic Research program(973 program)of China(Grant No.2013CB733201)the Key Program of the Chinese Academy of Sciences(KZZD-EW-05-01)the“Hundred Talents”program of Chinese Academy of Sciences for supporting the research
文摘Shallow fissures, being the main infiltration paths of fluid on the surface of a slope, played an important role in the whole process of a landslide. However, the spatial distribution characteristics of fissures in the slope are difficult to be determined. In this study, we attempted to characterize the variation pattern of slope fissures along depth in the Wenchuan earthquake area in Sichuan Province by combining engineering geological investigation, geomorphologic analysis and geophysical investigation. The geophysical methods that were used in this study include Multichannel Analysis of Surface Wave(MASW), Ground Penetrating Radar(GPR) and Electrical Resistivity Tomography(ERT). The results suggested that geophysical parameters(shear wave velocity, electromagnetic signals attenuation and resistivity) could provide valuable information for the spatial network of shallow fissures. Through the verification by engineering geological survey and geophysical sensitivity analysis, this work highlighted that MASW was the most appropriate technique to delineate the propagation of shallow fissures in a gravel soil slope.
基金This research is supported by the National Key Basic Research 973 Project(Grant No.:2008CB425701)the Special Project M7.0~8.0 of China Earthquake Administration
文摘Based on seismic data from the regional network for the last 34 years, we analyzed the present fault behavior of major fault zones around the Mabian area, southern Sichuan, and identified the risky fault-segments for potential future. The method of analysis is a combination of activity background of historical strong earthquakes mainly show ~ ( 1 ) The spatial distribution of b-values strong and large earthquakes in the spatial distribution of b-values with and current seismicity. Our results indicates significant heterogeneity in the studied area, which reflects the spatial difference of cumulative stress levels along various fault zones and segments. (2) Three anomalously low b-value areas with different dimensions were identified along the Mabian-Yanjin fault zone. These anomalies can be asperities under relatively high cumulated stress levels. Two asperities are located in the north of Mabian county, in Lidian town in western Muchuan county, and near Yanjin at the south end of the fault zone. These two areas represent potential large earthquake seismogenic sites around the Mabian area in the near future. Besides them, the third relatively smaller asperity is identified at southern Suijiang, as another potential strong- earthquake source. (3) An asperity along the southwestern segment of the Longquanshan fault zone indicates the site of potential moderate-to-strong earthquakes. (4) The asperity along the segment between Huangmu town in Hanyuan county and Longchi town in Emeishan city on Jinkouhe-Meigu fault has potential for a moderate-strong earthquake.
基金supported by the National Basic Research Program of China(Grant No.2011CB706605)State Key Program of National Natural Science Foundation of China(Grant No.51135007)+1 种基金Innovative Research Groups of the National Natural Science Foundation of Hubei Province(Grant No.2011CDA12)the Fundamental Research Funds for the Central Universities(Grant Nos.2012-Ia-017,2013-IV-014)for the support given to this research
文摘The influences of subsurface cracks,distributing along the axial direction,on the rolling contact fatigue(RCF)faliure in a bearing ring are investigated.A realistic three-dimensional model of the bearing ring containing three subsurface cracks is used to evaluate the fatigue crack propagation based on stress intensity factor(SIF)calculations.Moreover,the distributions of the subsurface cracks along the axial direction are varied to study their effects on RCF.The results provide valuable guidelines for enhanced understanding of RCF in bearings.
基金supported by National Natural Science Foundation of China(Grant Nos.41374055,41225010)Chinese University of Hong Kong Direct Grant for Research(Grant No. 3132771)+1 种基金HKSAR Research Grant Council ECS(Grant No.2191093) and GRF(Grant No.2130509)the Fundamental Research Funds for the Central Universities(Grant No.WK2080000053)
文摘On 25 April 2015,an M_w 7.8 earthquake occurred on the Main Himalaya Thrust fault with a dip angle of^7° about77 km northwest of Kathmandu,Nepal.This Nepal Gorkha event is the largest one on the Himalayan thrust belt since 1950.Here we use the compressive sensing method in the frequency domain to track the seismic radiation and rupture process of this event using teleseismic P waves recorded by array stations in North America.We also compute the distribution of static shear stress changes on the fault plane from a coseismic slip model.Our results indicate a dominant east-southeastward unilateral rupture process from the epicenter with an average rupture speed of ~3 km s^(-1).Coseismic radiation of this earthquake shows clear frequency-dependent features.The lower frequency(0.05-0.3 Hz) radiation mainly originates from large coseismic slip regions with negative coseismic shear stress changes.In comparison,higher frequency(0.3-0.6 Hz) radiation appears to be from the down-dip part around the margin of large slip areas,which has been loaded and presents positive coseismic shear stress changes.We propose an asperity model to interpret this Nepal earthquake sequence and compare the frequency-dependent coseismic radiation with that in subduction zones.Such frequency-dependent radiation indicates the depth-varying frictional properties on the plate interface of the Nepal section in the main Himalaya thrust system,similar to previous findings in oceanic subduction zones.Our findings provide further evidence of the spatial correlation between changes of static stress status on the fault plane and the observed frequency-dependent coseismic radiation during large earthquakes.Our results show that the frequency-dependent coseismic radiation is not only found for megathrust earthquakes in the oceanic subduction environment,but also holds true for thrust events in the continental collision zone.