Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

Formation mechanism and height calculation of the caved zone and water-conducting fracture zone in solid backfill mining

To study the heights of the caved zone and water-conducting fracture zone in backfill mining,the failure mechanism of strata during backfill mining was analyzed,and a method for determining the heights of the two zones was proposed based on key strata theory.The movement and failure regularity of the strata above the backfilling panel were revealed through numerical simulation.Considering the geologic conditions of the CT101 backfilling panel,the height of the fracture zone was determined using the proposed method along with empirical calculation,numerical simulation,and borehole detection.The results of the new calculation method were similar to in situ measurements.The traditional empirical formula,which is based on the equivalent mining height model,resulted in large errors during calculation.The findings indicate the reliability of the new method and demonstrate its significance for creating reference data for related studies.

Height prediction of water-flowing fracture zone with a geneticalgorithm support-vector-machine method

Prediction of the height of a water-flowing fracture zone(WFFZ)is the foundation for evaluating water bursting conditions on roof coal.By taking the Binchang mining area as the study area and conducting an in-depth study of the influence of coal seam thickness,burial depth,working face length,and roof category on the height of a WFFZ,we proposed that the proportion of hard rock in different roof ranges should be used to characterise the influence of roof category on WFFZ height.Based on data of WFFZ height and its influence index obtained from field observations,a prediction model is established for WFFZ height using a combination of a genetic algorithm and a support-vector machine.The reliability and superiority of the prediction model were verified by a comparative study and an engineering application.The results show that the main factors affecting WFFZ height in the study area are coal seam thickness,burial depth,working face length,and roof category.Compared with multiple-linear-regression and back-propagation neural-network approaches,the height-prediction model of the WFFZ based on a genetic-algorithm support-vector-machine method has higher training and prediction accuracy and is more suitable for WFFZ prediction in the mining area.

Early Pliocene Paleoceanography of the Vityaz Fracture Zone,Central Indian Ridge

Planktic foraminifera from the calcareous substrate of a ferromanganese crust in the Vityaz Fracture Zone （VFZ）, Central Indian Ridge were studied to reconstruct the early Pliocene paleoceanography of this region. Eleven species of planktic foraminifera were encountered, among them Globorotalia menardii, Neogloboquadrina dutertrei, Globigerina bulloides and Globigerinoides tuber are prominent. Predominance of N. dutertrei in the top 3 cm of the carbonate substrate is attributed to an influx of fresh water which eventually triggered their productivity by increasing the nutrient level. The presence of G. bulloides and G. menardii in significant proportions in deeper layers suggests the prevalence of open ocean upwelling. The bulk chemical compositions of the substrate at different depth intervals indicates higher enrichment of trace metals in the upper sections which could have been supplied through oceanic water by the chemical weathering of terrestrial matter during the peak of Pliocene Asian monsoon. Thus, it is concluded that during the early Pliocene the biogenic components of the substrate were distinctly contributed by both upwelling and productivity triggered by an influx of fresh water originating from the intensification of the Asian monsoon during the early Pliocene Period.

Response of underground pipeline through fault fracture zone to random ground motion

It is assumed that a pipeline is laid through a vertical fault fracture zone, and is excited by seismic ground motion modelled as stationary stochastic process. For horizontal incidence of waves, the cross-PSD （Power Spectral Density） function is developed using wave propagation theory, while for vertical incidence of waves the cross-PSD function is composed by auto-PSD model, coherence model and site response model. As the seismic input, the eross-PSD function is used to calculate the the axial and lateral seismic responses of underground pipeline through the fracture zone. The results show that the incident directions of seismic waves, width and soil property of the fracture zone have great influence on underground pipeline. It is suggested that the flexible joints with appropriate stiffness should be added into the pipeline near the interfaces between the fracture zone and the surrounded media.

Stresses and Shear Fracture Zone of Jinshazhou Tunnel Surrounding Rock in Rich Water Region

Field evidence has shown that large-scale and unstable discontinuous planes in the rock mass surrounding tunnels in rich water region are probably generated after excavation. The tunnel surrounding rock was divided into three zones, including elastic zone, plastic damage zone and shear fracture zone for assessing the stability of the tunnel surrounding rock. By local hydrogeology, the stresses of surrounding rock of Jinshazhou circular tunnel was analyzed and the stress solutions on the elastic and plastic damage zones were obtained by applying the theories of fluid-solid coupling and elasto-plastic damage mechanics. The shear fracture zone generated by joints was studied and its range was determined by using Molar-Coulomb strength criterion. Finally, the correctness of the theoretical results was validated by comparing the scopes of shear fracture zones calculated in this paper with those from literature.

Imaging Rock Density Distribution beneath Liwa Fracture Zone in the Southern Part of the Great Sumatran Fault System, Indonesia

We have imaged rock density distribution beneath Liwa fracture zone in the southern part of the the Sumatran Fault Zone by modelling and inverting Bouguer gravity data in two-and three-dimensional environments, respectively. The purpose of this study is aimed to figure out the subsurface distribution of rock densities associated with subsurface basement structure representing the evidence of trans-tensional tectonic product in the SF. In the gravity modeling, to eliminate distortions to the measured gravity values before modelling and inverting the data, Bouguer anomalies obtained in field measurements are reduced to the horizontal plane of z = +800 m as a representation of the average elevation in Liwa. For the inversion, we used algorithm implementing depth-and minimum volume weighting parameters in order to obtain a smooth model with better vertical resolution. The two-dimensional models show clearly surface topography of the basement rocks and the presence of normal faults. The reduced Bouguer anomaly of +800 m elevation shows the presence of structural lineaments extending primarily in a northwest-southeast direction, parallel to Sumatran Fault Zone and older graben faults showing a negative flower structure. From the three-dimensional inversion, the model illustrates an increase of density contrast, lower values being found near the surface and higher values in the deeper parts. The lower density contrast of 0.15 to 0.3 g/cm3 found in the rock groups at depths of 2 km and less is characteristic of relatively homogeneous and poorly compacted rocks. Rocks with moderate to high density contrast (>1.0 g/cm3) are recognized at depths of over 2 km. This model suggests a change of basement morphology as a function of depth, and delineates structural lineaments extending in northwest-southeast direction. This study supports the previous thought that Liwa area is underlain by graben structures, formed by trans-tensional tectonic events. Higher-density Tertiary volcanic breccia and lower-density Quaternary volcanic products of the Ranau Formation form the basement rocks and the overlying younger sediments, respectively.

Unconfined compressive strength and failure behaviour of completely weathered granite from a fault zone

Understanding the strength characteristics and deformation behaviour of the tunnel surrounding rock in a fault zone is significant for tunnel stability evaluation.In this study,a series of unconfined compression tests were conducted to investigate the mechanical characteristics and failure behaviour of completely weathered granite(CWG)from a fault zone,considering with height-diameter(h/d)ratio,dry densities(ρd)and moisture contents(ω).Based on the experimental results,a regression mathematical model of unconfined compressive strength(UCS)for CWG was developed using the Multiple Nonlinear Regression method(MNLR).The research results indicated that the UCS of the specimen with a h/d ratio of 0.6 decreased with the increase ofω.When the h/d ratio increased to 1.0,the UCS increasedωwith up to 10.5%and then decreased.Increasingρd is conducive to the improvement of the UCS at anyω.The deformation and rupture process as well as final failure modes of the specimen are controlled by h/d ratio,ρd andω,and the h/d ratio is the dominant factor affecting the final failure mode,followed byωandρd.The specimens with different h/d ratio exhibited completely different fracture mode,i.e.,typical splitting failure(h/d=0.6)and shear failure(h/d=1.0).By comparing the experimental results,this regression model for predicting UCS is accurate and reliable,and the h/d ratio is the dominant factor affecting the UCS of CWG,followed byρd and thenω.These findings provide important references for maintenance of the tunnel crossing other fault fractured zones,especially at low confining pressure or unconfined condition.

A method for predicting the water-flowing fractured zone height based on an improved key stratum theory

In the process of using the original key stratum theory to predict the height of a water-flowing fractured zone(WFZ),the influence of rock strata outside the calculation range on the rock strata within the calculation range as well as the fact that the shape of the overburden deformation area will change with the excavation length are ignored.In this paper,an improved key stratum theory(IKS theory)was proposed by fixing these two shortcomings.Then,a WFZ height prediction method based on IKS theory was established and applied.First,the range of overburden involved in the analysis was determined according to the tensile stress distribution range above the goaf.Second,the key stratum in the overburden involved in the analysis was identified through IKS theory.Finally,the tendency of the WFZ to develop upward was determined by judging whether or not the identified key stratum will break.The proposed method was applied and verified in a mining case study,and the reasons for the differences in the development patterns between the WFZs in coalfields in Northwest and East China were also fully explained by this method.

Failure responses of rock tunnel faces during excavation through the fault-fracture zone

It is essential to cast light on the construction risks in tunnel excavations through the fault-fracture zone(FFZ).This study adopts the material point method(MPM)to simulate the failure responses of a rock tunnel face during excavation through the FFZ.A numerical study was conducted to compare a physical model test and validate the feasibility of using the MPM in simulating tunnel face failure.One hundred ninety numerical simulation cases were constructed to represent a rock tunnel excavation project with different site con-figurations.The simulation results suggest that the cohesion and the friction angle significantly influence failure responses.The tunnel cover depth can magnify the failure responses,and the FFZ thickness significantly affects the mobilized rock mass volume when the FFZ consists of a weak rock mass.The numerical simulation results suggest three deformation patterns:face bulge,partial failure,and slide collapse.The failure responses can be characterized by stress arch,slip surface,angle of reposing,and influence range.The insights suggested by the face failure responses during excavation through the FFZ can aid field engineers in determining the scope of possible damage,and in establishing emergency measures to minimize losses if such failure occurs.

Relationship between Remotely Sensed Vegetation Change and Fracture Zones Induced by the 2008 Wenchuan Earthquake, China

The Wenchuan earthquake triggered cascading disasters of landslides and debris flows that caused severe vegetation damage. Fracture zones can affect geodynamics and spatial pattern of vegetation damage. A segment tracing algorithm method was applied for identifying the regional fracture system through lineament extractions from a shaded digital elevation model with 25 m mesh for southern Wenchuan. Remote sensing and geographic information system techniques were used to analyze the spatiotemporal vegetation pattern. The relationship between vegetation type identified from satellite images and lineament density was used to characterize the distribution patterns of each vegetation type according to fracture zones. Broad-leaved forest, mixed forest, and farmland persist in areas with moderate lineament density. Deciduous broad-leaved and coniferous forest persists in less frac- tured areas. Shrub and meadow seem to be relatively evenly distributed across all lineament densities.Meadow, farmland, and shrub persist in the fractured areas. Changes of spatial structure and correlation between vegetation patterns before and after the earthquake were examined using semivariogram analysis of normalized difference vegetation indices derived from Landsat enhanced thematic mapper images. The sill values of the semivariograms show that the spatial heterogeneity of vegetation covers increased after the earthquake. Moreover, the anisotropic behaviors of the semivariograms coincide with the vegetation changes due to the strikes of fracture zones.

Seismic response of tunnel near fault fracture zone under incident SV waves

This study investigated the impact of a non-causative fault on the dynamic response of a nearby lined tunnel under the incidence of plane SV waves using the indirect boundary element method.The effects of several critical parameters,such as the incident frequency,the inclination degree of the fault,the distance between the fault and the tunnel on the hoop stress of the lined inner and outer walls,were explored intensively.The numerical results indicated that the non-causative fault could significantly change the hoop stress distribution of inner and outer surfaces of the tunnels.In general,for the vertically incident seismic waves,when the tunnel was located in the foot wall(under the fault),the hoop stress within the tunnel was significantly greater than that of the tunnels in the non-fault half space,with an amplification factor of up to 117%.The amplification effect became more pronounced as the fault dip angle increased.However,when the tunnel was located in the hanging wall(above the fault),the non-causative fault could produce a significant shielding effect on the dynamic response of the tunnel under high frequency wave incidence,with the reduction of hoop stress being up to 81%.For lowfrequency waves,though,the fault could lead to an increase of the hoop stress of the tunnel of up to 152%.The research results will provide a reference for the seismic design and safety protection of underground structures in non-causative fault sites.

Predicting the height of water-flow fractured zone during coal mining under the Xiaolangdi Reservoir

It is very important to determine the extent of the fractured zone through which water can flow before coal mining under the water bodies.This paper deals with methods to obtain information about overburden rock failure and the development of the fractured zone while coal mining in Xin'an Coal Mine.The risk of water inrush in this mine is great because 40%of the mining area is under the Xiaolangdi reservoir.Numerical simulations combined with geophysical methods were used in this paper to obtain the development law of the fractured zone under different mining conditions.The comprehensive geophysical method described in this paper has been demonstrated to accurately predict the height of the water-flow fractured zone.Results from the new model, which created from the results of numerical simulations and field measurements,were successfully used for making decisions in the Xin'an Coal Mine when mining under the Xiaolangdi Reservoir.Industrial scale experiments at the number 11201,14141 and 14191 working faces were safely carried out.These achievements provide a successful background for the evaluation and application of coal mining under large reservoirs.

Structural failure mechanism and strengthening method of fracture plugging zone for lost circulation control in deep naturally fractured reservoirs

Focused on the lost circulation control in deep naturally fractured reservoirs, the multiscale structure of fracture plugging zone is proposed based on the theory of granular matter mechanics, and the structural failure pattern of plugging zone is developed to reveal the plugging zone failure mechanisms in deep, high temperature, high pressure, and high in-situ stress environment. Based on the fracture plugging zone strength model, key performance parameters are determined for the optimal selection of loss control material(LCM). Laboratory fracture plugging experiments with new LCM are carried out to evaluate the effect of the key performance parameters of LCM on fracture plugging quality. LCM selection strategy for fractured reservoirs is developed. The results show that the force chain formed by LCMs determines the pressure stabilization of macro-scale fracture plugging zone. Friction failure and shear failure are the two major failure patterns of fracture plugging zone. The strength of force chain depends on the performance of micro-scale LCM, and the LCM key performance parameters include particle size distribution, fiber aspect ratio, friction coefficient, compressive strength, soluble ability and high temperature resistance. Results of lab experiments and field test show that lost circulation control quality can be effectively improved with the optimal material selection based on the extracted key performance parameters of LCMs.

Structural formation and evolution mechanisms of fracture plugging zone

A coupled CFD-DEM method is used to simulate the formation process of fracture plugging zone.A photo-elastic system characterizing mesoscale force chain network developed by our own is used to model the pressure evolution in fracture plugging zone to reveal the evolution mechanism of the structure of fracture plugging zone.A theoretical basis is provided for improving the lost circulation control effect in fractured reservoirs and novel methods are proposed for selecting loss control materials and designing loss control formula.CFD-DEM simulation results show that bridging probability is the key factor determining the formation of fracture plugging zone and fracture plugging efficiency.Critical and absolute bridging concentrations are proposed as the key indexes for loss control formula design.With the increase of absolute bridging concentration,the governing factor of bridging is changed from material grain size to the combination of material grain size and friction force.Results of photo-elastic experiments show that mesoscale force chain network is the intrinsic factor affecting the evolution of pressure exerting on the fracture plugging zone and determines the macroscopic strength of fracture plugging zone.Performance parameters of loss control material affect the force chain network structure and the ratio of stronger force chain,and further impact the stability and strength of fracture plugging zone.Based on the study results,the loss control formula is optimized and new-type loss control material is designed.Laboratory experiments results show that the fracture plugging efficiency and strength is effectively improved.

In-situ observations of damage-fracture evolution in surrounding rock upon unloading in 2400-m-deep tunnels

The damage-fracture evolution of deep rock mass has obvious particularity,which is revealed in 2400-mdeep tunnels by field tests.The evolution of the excavation damaged zone depth is consistent with that of the fractured zone depth.The ratio of the excavation damaged zone depth to the excavation fractured zone depth is greater than 2.0 in a rock mass with both high strength and good integrity,but less than1.5 in a rock mass with lower strength or poor integrity.Zonal disintegration in a rock mass with high strength and fair integrity is more likely to occur when it contains more than two groups of primary fractures in damaged zones.Fractures develop outward in zonal disintegration but are totally different from the single-zone fracture,in which the fractures develop inward,and it is the starting position of the fractured zone when the excavation surface of the middle pilot is 7–9 m close to the pre-set borehole and it stops after the excavation surface of the baseplate is 11–14 m away.The most intense evolution occurs around 2–4 m from the pre-set borehole in the sidewall expansion stage.The research results provide a reference for the monitoring scheme and support design of CJPL-Ⅲin its future construction.

ROCK FRACTURE PROCESS ZONE AND ITS INFLUENCE ON MODEⅡFRACTURE BEHAVIOR

On the basis of that rock material usually has a larger fracture process zone,a new fracture criterion which is different from that of linear elastic fracture theory was presented.On this basis,the fracture behavior and influencing factors under modeⅡor compressive shear loading were investigated.

Height Detection and Analysis of Water Flowing Fractured Zone of Coal Face

Taking 91105 working face as the research object, the observation method of water flowing fractured zone and the layout of mining holes were determined by analyzing the field geological structure. It was shown that the fractured zone height and the ratio given by the measured method were 52.33 and 12.46, respectively. By the numerical simulation method with the software of UDEC, the fractured zone height and the ratio were 42.5 and 10.12. By comparison of measured height data and UDEC numerical simulation, there were some differences between the measured height and the calculated results of UDEC numerical simulation method. The method of simulation can be used as the technical basis for the design of waterproof coal pillar in the future.

A fracture model for assessing tensile mode crack growth resistance of rocks

Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict tensile mode(mode I)failure loads of cracked rocks.The basic principle of the model is to estimate the reference crack corresponding to the fracture process zone(FPZ)based on the maximum normal strain(MNSN)ahead of the crack tip,and then use the effective crack to calculate the fracture toughness.We emphasize that the non-singular stress/strain terms should be considered in the description of the MNSN.In this way,the FPZ,non-singular terms and the biaxial stress state at the crack tip are simul-taneously considered.The principle of the model is explicit and easy to apply.To verify the proposed model,laboratory experiments were performed on a rock material using six groups of specimens.The model predicted the specimen geometry dependence of the measured fracture toughness well.More-over,the potential of the model in analyzing the size effect of apparent fracture toughness was discussed and validated through experimental data reported in the literature.The model was demonstrated su-perior to some commonly used fracture models and is an excellent tool for the safety assessment of cracked rock structures.

Fracture development around deep underground excavations: Insights from FDEM modelling

Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method （FDEM） code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone （EDZ） around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.