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.

Tensile Fractures and in situ Stress Measurement Data Constraints on Cretaceous-Present Tectonic Stress Field Evolution of the Tanlu Fault Zone in Shandong Province,North China Craton

Tectonic stress fields are the key drivers of tectonic events and the evolution of regional structures.The tectonic stress field evolution of the Tanlu fault zone in Shandong Province,located in the east of the North China Craton(NCC),may have preserved records of the NCC’s tectonic history.Borehole television survey and hydraulic fracturing were conducted to analyze the paleo and present tectonic stress fields.Three groups of tensile fractures were identified via borehole television,their azimuths being NNW-SSE,NW-SE and NE-SW,representing multiple stages of tectonic events.Hydraulic fracturing data indicates that the study region is experiencing NEE-SWW-oriented compression and nearly-N-Soriented extension,in accordance with strike-slip and compression.Since the Cretaceous,the orientation of the extensional stress has evolved counterclockwise and sequentially from nearly-NW-SE-oriented to NE-SW-oriented and even nearly N-S-oriented,the stress state having transitioned from strike-slip-extension to strike-slip-compression,in association with the rotating and oblique subduction of the Pacific Plate beneath the NCC,with the participation of the Indian Plate.

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.

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.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

Comprehensive assessment on dynamic roof instability under fractured rock mass conditions in the excavation disturbed zone

The damage process of fractured rock mass showed that the fracture in rocks induced roof collapse in Yangchangwan Coal Mine, China. The rock mass was particularly weak and fractured. There occurred 6 large-scale dynamical roof falls in the excavation disturbed zone （EDZ） with the collapsing volume of 216 m^3. First, the field detailed geological environment, regional seismic dynamics, and dynamic instability of roadways were generally investigated. Second, the field multiple-index monitoring measurements for detecting the deep delamination of the roof, convergence deformation, bolt-cable load, acoustic emission （AE） characteristic parameters, total AE events, AE energy-releasing rate, rock mass fracture, and damage were arranged. Finally, according to the time-space-strength relations, a quantitative assessment of the influence of rock-mass damage on the dynamic roof instability was accomplished.

Mechanical mechanism of overlying strata breaking and development of fractured zone during close-distance coal seam group mining

This study mainly investigates the mechanical mechanism of overlying strata breaking and the development of fractured zones during close-distance coal seam group mining in the Gaojialiang coal mine.First,a mechanical model for the second"activation"of broken overlying strata is established,and the related mechanical"activation"conditions are obtained.A recursive formula for calculating the separation distance of overlying strata is deduced.Second,a height determining method for predicting the height of fractured zones during close-distance coal seam group mining is proposed based on two values,namely,the separation distance and ultimate subsidence value of overlying strata.This method is applied to calculate the fractured zone heights in nos.20107 and 20307 mining faces.The calculated results are almost equal to the field observation results.Third,a modified formula for calculating the height of a waterflowing fractured zone is proposed.A comparison of the calculated and observed results shows that the errors are small.The height determining method and modified formula not only build a theoretical foundation for water conservation mining at the Gaojialiang coal mine,but also provide a reference for estimating the height of water-flowing fractured zones in other coal mines with similar conditions.

Overburden fracture evolution laws and water-controlling technologies in mining very thick coal seam under water-rich roof

Considering the danger of water inrush in mining very thick coal seam under water-rich roof in Majialiang Coal Mine,the universal discrete element(UDEC)software was used to simulate the overburden fracture evolution laws when mining 4#coal seam.Besides,this study researched on the influence of face advancing length,speed and mining height on the height of the water flowing fractured zones(HWFFZ),and analyzed the correlation of face advancing length and change rules of aquifer water levels and goaf water inflow.Based on those mentioned above,this research proposed the following water-controlling technologies:draining the roof water before mining,draining goaf water,reasonable advancing speed and mining thickness.These water-controlling technologies were successfully used in the feld,thus ensured safely mining the very thick coal seam under water-rich roof.

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.

Fractal Feature of Western Fracture Zone in Xikuangshan Antimony Mine and its Geological Significance

In Xikuangshan antimony ore field, the western fracture zone is a composite of major fault, F 75 , and its secondary faults, such as F 71 , F 72 and F 3 etc.. On plane, the fracture zone scatters from southwest to northeast, and concentrates from upper to deeper level on profile. All ore bodies exist in the carbonate of footwall of the major fault or that of the footwall of its secondary faults. From 480 m and 320 m to 120 m level, the fractal dimensional number of the fault system decreases from 1 482 2 and 1 448 6 to 1 339 2, which indicates the form of fracture zone becoming more simple at deeper level. And in five sub ranges, the Ⅲ and Ⅳ sub ranges are the known area, and the Ⅰ, Ⅱ and Ⅴ sub ranges are unknown. The fractal studies of the western fracture zone in these sub ranges show that the fractal dimensional numbers of the Ⅰ and Ⅱ, being 1 201 5 and 1 278 0, respectively, are smaller than that of the Ⅲ and Ⅳ, being 1 475 9 and 1 576 9, respectively; and that of the Ⅴ, being 1 571 2, keeps with that of the Ⅲ, Ⅳ sub ranges. So mineralization is not well in Ⅰ and Ⅱ sub ranges, and Ⅴ sub range is the best to benefit mineralization.

Height of fractured zone inside overlying strata under high-intensity mining in China

The height of fractured zone(HFZ) at the high-intensity longwall mining panel plays a vital role in the safety analysis of coal mining under bodies of water. This paper described definitions of the highintensity mining. The processes of overburden failure transfer(OFT) were analyzed, which were divided into the development stage and the termination stage. Through theoretical analysis, the limited suspension-distance and the limited overhanging distance were proposed to judge the damage of each stratum. Mechanical models of strata suspended integrity and overhanging stability were established.A theoretical method to predict the HFZ at the high-intensity longwall mining panel was put forward based on the processes of OFT. Taking a high-intensity longwall mining panel(No. 11915 panel) as an example, the theoretical method proposed, the engineering analogy and the empirical formulas in the Regulation were used to predict the HFZ. The results show that the theoretical result is consistent with the engineering analogies' result and empirical formulas' result. The rationality and reliability of the theoretical method proposed is verified.

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.

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.

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.

INFLUENCE OF LOCAL BRITTLE ZONE ON THE FRACTURETOUGHNESS OF HIGH-STRENGTH LOW-ALLOYEDMULTIPASS WELD METALS

In this paper, toughness properties and microstructurc of low-alloyed multipass welds with yield strength above 700MPa have 6een studied using the weld thermal simulation and throughout thickness CTOD fracture mechanics tests. Impact testing of thermal simulated specimens showed that the primary weld metal and the fine gmmed weld metal had good toughness, while the coarse grained weld metal had the lowest toughness value as the local brittle zone (LBZ) in multipass weld metals. Cleavage fracture in CTOD testing of thick multipass weld metals was initiated from martensite-austenite (MA) phases in the LBZ. MA phases were distributed at the prior austenite grain boundaries and around ferrite grains. As the size of the local brittle zone along the fatigue crack front increases, CTOD frncture toughness of multipass weld metals decreases. The weakest link theory was used to evaluate effect of the local brittle zone on fracture toughness of thick multipass weld metals. The estimated curves agree well with the eaperimental data.

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.

Fractured zone height of longwall mining and its effects on the overburden aquifers

As mining depth becomes deeper and deeper,the possibility of undermining overburden aquifers is increasing.It is very important for coal miners to undertake studies on the height of fractured zone during longwall mining and the effects of longwall mining on the underground water while mining under surface water bodies and underground aquifers.In order to study this problem,piezometers for monitoring underground water levels were installed above the longwall panels in an American coalmine.Large amounts of pre-mining,during mining and post-mining monitoring data were collected.Based on the data,the heights of fractured zones were obtained and the effects of longwall mining on the underground water were studied.The results demonstrate that when the piezometer monitoring wells had an interburden thickness of less than 72.7 m,the groundwater level decreased immediately to immeasurable levels and the wells went dry after undermining the face of longwall.The height of the fractured zone is 72.7-85.3 m in the geological and mining conditions.The results also show that the calculated values of fractured zones by the empirical formulae used in China are smaller than the actual results.Therefore,it is not always safe to use them for analyses while mining under water bodies.

Investigation of FRP and SFRC Technologies for Efficient Tunnel Reinforcement Using the Cohesive Zone Model

Amid urbanization and the continuous expansion of transportation networks,the necessity for tunnel construction and maintenance has become paramount.Addressing this need requires the investigation of efficient,economical,and robust tunnel reinforcement techniques.This paper explores fiber reinforced polymer(FRP)and steel fiber reinforced concrete(SFRC)technologies,which have emerged as viable solutions for enhancing tunnel structures.FRP is celebrated for its lightweight and high-strength attributes,effectively augmenting load-bearing capacity and seismic resistance,while SFRC’s notable crack resistance and longevity potentially enhance the performance of tunnel segments.Nonetheless,current research predominantly focuses on experimental analysis,lacking comprehensive theoretical models.To bridge this gap,the cohesive zone model(CZM),which utilizes cohesive elements to characterize the potential fracture surfaces of concrete/SFRC,the rebar-concrete interface,and the FRP-concrete interface,was employed.A modeling approach was subsequently proposed to construct a tunnel segment model reinforced with either SFRC or FRP.Moreover,the corresponding mixed-mode constitutive models,considering interfacial friction,were integrated into the proposed model.Experimental validation and numerical simulations corroborated the accuracy of the proposed model.Additionally,this study examined the reinforcement design of tunnel segments.Through a numerical evaluation,the effectiveness of innovative reinforcement schemes,such as substituting concrete with SFRC and externally bonding FRP sheets,was assessed utilizing a case study from the Fuzhou Metro Shield Tunnel Construction Project.

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.

Analysis of fracture process zone in brittle rock subjected to shear-compressive loading

An analytical expression for the prediction of shear-compressive fracture process zone(SCFPZ) is derived by using a proposed local strain energy density criterion, in which the strain energy density is separated into the dilatational and distortional strain energy density, only the former is considered to contribute to the brittle fracture of rock in different loading cases. The theoretical prediction by this criterion shows that the SCFPZ is of asymmetric mulberry leaf in shape, which forms a shear-compression fracture kern. Dilatational strain energy density along the boundary of SCFPZ reaches its maximum value. The dimension of SCFPZ is governed by the ratio of K_Ⅱ to (K_Ⅰ.) The analytical results are then compared with those from literatures and the tests conducted on double edge cracked Brazilian disk subjected to diametrical compression. The obtained results are useful to the prediction of crack extension and to nonlinear analysis of shear-compressive fracture of brittle rock.