Development and application of novel high‐efficiency composite ultrafine cement grouts for roadway in fractured surrounding rocks

The fractured surrounding rocks of roadways pose major challenges to safe mining.Grouting has often been used to reinforce the surrounding rocks to mitigate the safety risks associated with fractured rocks.The aim of this study is to develop highly efficient composite ultrafine cement(CUC)grouts to reinforce the roadway in fractured surrounding rocks.The materials used are ultrafine cement(UC),ultrafine fly ash(UF),ultrafine slag(US),and additives(superplasticizer[SUP],aluminate ultrafine expansion agent[AUA],gypsum,and retarder).The fluidity,bleeding,shrinkage,setting time,chemical composition,microstructure,degree of hydration,and mechanical property of grouting materials were evaluated in this study.Also,a suitable and effective CUC grout mixture was used to reinforce the roadway in the fractured surrounding rock.The results have shown that the addition of UF and US reduces the plastic viscosity of CUC,and the best fluidity can be obtained by adding 40%UF and 10%US.Since UC and UF particles are small,the pozzolanic effect of UF promotes the hydration reaction,which is conductive to the stability of CUC grouts.In addition,fine particles of UC,UF,and US can effectively fill the pores,while the volumetric expansion of AUA and gypsum decreases the pores and thus affects the microstructure of the solidified grout.The compressive test results have shown that the addition of specific amounts of UF and US can ameliorate the mechanical properties of CUC grouts.Finally,the CUC22‐8 grout was used to reinforce the No.20322 belt roadway.The results of numerical simulation and field monitoring have indicated that grouting can efficaciously reinforce the surrounding rock of the roadway.In this research,high‐performance CUC grouts were developed for surrounding rock reinforcement of underground engineering by utilizing UC and some additives.

Surrounding rock pressure in the tunnel portal section through moraine under freeze-thaw action

Moraines,characterized by the accumulation of rock and soil debris transported by glacial activity,present unique challenges for tunnel construction,particularly in portal sections,due to prevailing geographical and climatic conditions that facilitate freeze-thaw action.Despite these challenges,there is a dearth of studies investigating the influence of freeze-thaw action and water content on the mechanical properties of moraines,and no research on calculating surrounding rock pressure in moraine tunnels subjected to freeze-thaw conditions.In this study,direct shear tests under freeze-thaw cycles were conducted to examine the effects of freeze-thaw cycles and water content on the mechanical properties of frozen moraine.A comprehensive parameter K,integrating the number of freeze-thaws and water content,was introduced to model cohesion c.Drawing on Terzaghi Theory,we propose an improved algorithm for calculating surrounding rock pressure at the portal section of moraine tunnels.Using a tunnel as a case study,surrounding rock pressure was calculated under various conditions to validate the Improved Algorithm's efficacy.The results show that:(1)Strength loss exhibits a linear trend with the number of freeze-thaw cycles at water content levels of 4%and 8%,while at 12%water content,previous freeze-thaw cycles induce more significant damage to the soil.(2)Moraine saturation peaks between 8%and 12%water content.Following repeated freeze-thaw cycles,moraine shear strength initially increases before decreasing with varying water content.(3)The internal friction angle of moraine experiences slight reductions with prolonged freeze-thaw cycles,but both freeze-thaw cycles and water content significantly influence cohesion.(4)Vertical surrounding rock pressure increases after the initial freeze-thaw cycle,particularly with higher water content,although freeze-thaw cycles have minimal effect on it.(5)Freeze-thaw cycles lead to a substantial increase in lateral surrounding rock pressure,necessitating reinforced support structures at the arch wall,arch waist,and arch foot in engineering projects to mitigate freeze-thaw effects.This study provides a foundation for designing and selecting tunnel support structures in similar geological conditions.

Sub-Homogeneous Peridynamic Model for Fracture and Failure Analysis of Roadway Surrounding Rock

The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel function to represent accurately the spatial decline of long-range force.Additionally,modifications to the traditional bondbased PD model are made.By considering the micro-structure of coal-rock materials within a uniform discrete model,heterogeneity characterized by bond random pre-breaking is introduced.This approach facilitates the proposal of a novel model capable of handling the random distribution characteristics of material heterogeneity,rendering the PD model suitable for analyzing the deformation and failure of heterogeneous layered coal-rock mass structures.The established numerical model and simulation method,termed the sub-homogeneous PD model,not only incorporates the support effect but also captures accurately the random heterogeneous micro-structure of roadway surrounding rock.The simulation results obtained using this model show good agreement with field measurements from the Fucun coal mine,effectively validating the model’s capability in accurately reproducing the deformation and failure mode of surrounding rock under bolt-supported(anchor cable).The proposed subhomogeneous PD model presents a valuable and effective simulation tool for studying the deformation and failure of roadway surrounding rock in coal mines,offering new insights and potential advancements.

Numerical Simulation of Surrounding Rock Deformation and Grouting Reinforcement of Cross-Fault Tunnel under Different Excavation Methods

Tunnel construction is susceptible to accidents such as loosening, deformation, collapse, and water inrush, especiallyunder complex geological conditions like dense fault areas. These accidents can cause instability and damageto the tunnel. As a result, it is essential to conduct research on tunnel construction and grouting reinforcementtechnology in fault fracture zones to address these issues and ensure the safety of tunnel excavation projects. Thisstudy utilized the Xianglushan cross-fault tunnel to conduct a comprehensive analysis on the construction, support,and reinforcement of a tunnel crossing a fault fracture zone using the three-dimensional finite element numericalmethod. The study yielded the following research conclusions: The excavation conditions of the cross-fault tunnelarray were analyzed to determine the optimal construction method for excavation while controlling deformationand stress in the surrounding rock. The middle partition method (CD method) was found to be the most suitable.Additionally, the effects of advanced reinforcement grouting on the cross-fault fracture zone tunnel were studied,and the optimal combination of grouting reinforcement range (140°) and grouting thickness (1m) was determined.The stress and deformation data obtained fromon-site monitoring of the surrounding rock was slightly lower thanthe numerical simulation results. However, the change trend of both sets of data was found to be consistent. Theseresearch findings provide technical analysis and data support for the construction and design of cross-fault tunnels.

Heterogeneous information phase space reconstruction and stability prediction of filling body–surrounding rock combination

Traditional research believes that the filling body can effectively control stress concentration while ignoring the problems of unknown stability and the complex and changeable stress distribution of the filling body–surrounding rock combination under high-stress conditions.Current monitoring data processing methods cannot fully consider the complexity of monitoring objects,the diversity of monitoring methods,and the dynamics of monitoring data.To solve this problem,this paper proposes a phase space reconstruction and stability prediction method to process heterogeneous information of backfill–surrounding rock combinations.The three-dimensional monitoring system of a large-area filling body–surrounding rock combination in Longshou Mine was constructed by using drilling stress,multipoint displacement meter,and inclinometer.Varied information,such as the stress and displacement of the filling body–surrounding rock combination,was continuously obtained.Combined with the average mutual information method and the false nearest neighbor point method,the phase space of the heterogeneous information of the filling body–surrounding rock combination was then constructed.In this paper,the distance between the phase point and its nearest point was used as the index evaluation distance to evaluate the stability of the filling body–surrounding rock combination.The evaluated distances(ED)revealed a high sensitivity to the stability of the filling body–surrounding rock combination.The new method was then applied to calculate the time series of historically ED for 12 measuring points located at Longshou Mine.The moments of mutation in these time series were at least 3 months ahead of the roadway return dates.In the ED prediction experiments,the autoregressive integrated moving average model showed a higher prediction accuracy than the deep learning models(long short-term memory and Transformer).Furthermore,the root-mean-square error distribution of the prediction results peaked at 0.26,thus outperforming the no-prediction method in 70%of the cases.

Assessing the range of blasting-induced cracks in the surrounding rock of deeply buried tunnels based on the unified strength theory

Blasting-induced cracks in the rock surrounding deeply buried tunnels can result in water gushing and rock mass collapse,posing significant safety risks.However,previous theoretical studies on the range of blasting-induced cracks often ignore the impact of the in-situ stress,especially that of the intermediate principal stress.The particle displacement−crack radius relationship was established in this paper by utilizing the blasthole cavity expansion equation,and theoretical analytical formulas of the stress−displacement relationship and the crack radius were derived with unified strength theory to accurately assess the range of cracks in deep surrounding rock under a blasting load.Parameter analysis showed that the crushing zone size was positively correlated with in-situ stress,intermediate principal stress,and detonation pressure,whereas negatively correlated with Poisson ratio and decoupling coefficient.The dilatancy angle-crushing zone size relationship exhibited nonmonotonic behavior.The relationships in the crushing zone and the fracture zone exhibited opposite trends under the influence of only in-situ stress or intermediate principal stress.As the in-situ stress increased from 0 to 70 MPa,the rate of change in the crack range and the attenuation rate of the peak vibration velocity gradually slowed.

Prediction model of surrounding rock deformation in doublecontinuous-arch tunnel based on the ABC-WNN

Purpose–The wavelet neural network(WNN)has the drawbacks of slow convergence speed and easy falling into local optima in data prediction.Although the artificial bee colony(ABC)algorithm has strong global optimization ability and fast convergence speed,it also has the drawbacks of slow speed while finding the optimal solution and weak optimization ability in the later stage.Design/methodology/approach–This article uses an ABC algorithm to optimize the WNN and establishes an ABC-WNN analysis model.Based on the example of the Jinan Yuhan underground tunnel project,the deformation of the surrounding rock of the double-arch tunnel crossing the fault fracture zone is predicted and analyzed,and the analysis results are compared with the actual detection amount.Findings–The comparison results show that the predicted values of the ABC-WNN model have a high degree of fitting with the actual engineering data,with a maximum relative error of only 4.73%.On this basis,the results show that the statistical features of ABC-WNN are the lowest,with the errors at 0.566 and 0.573,compared with the single back propagation(BP)neural network model and WNN model.Therefore,it can be derived that the ABC-WNN model has higher prediction accuracy,better computational stability and faster convergence speed for deformation.Originality/value–This article uses firstly the ABC-WNN for the deformation analysis of double-arch tunnels.This attempt laid the foundation for artificial intelligence prediction in deformation analysis of multiarch tunnels and small clearance tunnels.It can provide a new and effective way for deformation prediction in similar projects.

Study on creep deformation and energy development of underground surrounding rock under four‐dimensional support

There is an urgent need to develop optimal solutions for deformation control of deep high‐stress roadways,one of the critical problems in underground engineering.The previously proposed four‐dimensional support(hereinafter 4D support),as a new support technology,can set the roadway surrounding rock under three‐dimensional pressure in the new balanced structure,and prevent instability of surrounding rock in underground engineering.However,the influence of roadway depth and creep deformation on the surrounding rock supported by 4D support is still unknown.This study investigated the influence of roadway depth and creep deformation time on the instability of surrounding rock by analyzing the energy development.The elastic strain energy was analyzed using the program redeveloped in FLAC3D.The numerical simulation results indicate that the combined support mode of 4D roof supports and conventional side supports is highly applicable to the stability control of surrounding rock with a roadway depth exceeding 520 m.With the increase of roadway depth,4D support can effectively restrain the area and depth of plastic deformation in the surrounding rock.Further,4D support limits the accumulation range and rate of elastic strain energy as the creep deformation time increases.4D support can effectively reduce the plastic deformation of roadway surrounding rock and maintain the stability for a long deformation period of 6 months.As confirmed by in situ monitoring results,4D support is more effective for the long‐term stability control of surrounding rock than conventional support.

Surrounding Rock Control Technology of Strong Dynamic Pressure Roadway in Hudi Coal Industry

Aiming at the problems of large deformation and difficult maintenance of deep soft rock roadway under the influence of high ground stress and strong dynamic pressure, taking the surrounding rock control of 1105 lane in Hudi Coal Industry as an example, the deformation characteristics and surrounding rock control measures of deep soft rock roadway are analyzed and discussed by means of geological data analysis, roadway deformation monitoring, rock crack drilling and field test. The results show that the main causes of roadway deformation are high ground stress, synclinal tectonic stress, advance mining stress, roadway penetration and surrounding rock fissure development. Based on the deformation characteristics and mechanism of lane 1105, the supporting countermeasures of “roof synergic support, layered grouting, anchor cable beam support, closed hardening of roadway surface” are proposed, which can provide reference for the control of deep roadway surrounding rock under similar conditions.

Surrounding Rock Failure Mechanism and Control Technology of Gob-Side Entry with Triangle Coal Pillar at Island Longwall Panel in 15 m Extra-Thick Coal Seams

Taking the return air roadway of Tashan 8204 isolated island working face as the background, the evolution law of the stress field in the surrounding rock of the widened coal pillar area roadway during the mining period of the isolated island working face is obtained through numerical simulation. The hazardous area of strong mine pressure under different coal pillar widths is determined. Through simulation, it is known that when the width of the coal pillar is less than 20 m, there is large bearing capacity on the coal side of the roadway entity. The force on the side of the coal pillar is relatively small. When the width of the coal pillar ranges from 25 m to 45 m, the vertical stress on the roadway and surrounding areas is relatively high. Pressure relief measures need to be taken during mining to reduce surrounding rock stress. When the width of the coal pillar is greater than 45 m, the peak stress of the coal pillar is located in the deep part of the surrounding rock, but it still has a certain impact on the roadway. It is necessary to take pressure relief measures to transfer the stress to a deeper depth to ensure the stability of the triangular coal pillar during the safe mining period of the working face. This provides guidance for ensuring the stability of the triangular coal pillar during the safe mining period of the working face.

Forward prediction for tunnel geology and classification of surrounding rock based on seismic wave velocity layered tomography

Excavation under complex geological conditions requires effective and accurate geological forward-prospecting to detect the unfavorable geological structure and estimate the classification of surround-ing rock in front of the tunnel face.In this work,a forward-prediction method for tunnel geology and classification of surrounding rock is developed based on seismic wave velocity layered tomography.In particular,for the problem of strong multi-solution of wave velocity inversion caused by few ray paths in the narrow space of the tunnel,a layered inversion based on regularization is proposed.By reducing the inversion area of each iteration step and applying straight-line interface assumption,the convergence and accuracy of wave velocity inversion are effectively improved.Furthermore,a surrounding rock classification network based on autoencoder is constructed.The mapping relationship between wave velocity and classification of surrounding rock is established with density,Poisson’s ratio and elastic modulus as links.Two numerical examples with geological conditions similar to that in the field tunnel and a field case study in an urban subway tunnel verify the potential of the proposed method for practical application.

Failure mechanism and control countermeasures of surrounding rock at deep large section chamber intersection in the Wanfu Coal Mine

In order to solve the problem of large deformation at the intersection of deep large section soft rock roadway,this paper takes the intersection of kilometer-deep roadway in the Wanfu Coal Mine as an engineering example and applies Negative Poisson’s Ratio(NPR)steel anchor cable in roadway support for the first time.By combining numerical simulation indoor test,theoretical analysis and field test,the deformation mechanism of surrounding rock at the intersection of deep-buried roadway was analyzed,and the control strategy with micro NPR steel anchor cable as the core was put forward.Through numerical simulation,the numerical analysis model of roadway intersection with different intersection angles and excavation sequence was constructed,and the impact of two key variables of rake angle and excavation sequence on the stability of surrounding rock at roadway intersection was studied.The optimal dip angle is 90°and the optimal excavation sequence was determined as pump house-pump house passage-substation.The mechanical properties of the micro-NPR steel anchor cable were studied through the static tensile test in the laboratory.The results showed that the micro-NPR steel anchor cable showed high constant resistance,uniform tensile,no yield platform,and no obvious necking phenomenon during breaking.Through theoretical derivation,it was calculated that the vertical stress of roadway intersection is 45 MPa,and the bearing capacity of superposed arch composed of micro NPR steel anchor cable is 1257 kN,which is enough to guarantee the overall stability of intersection.Support application test and monitoring were carried out on site,and it was verified that the combined support strategy of short and long micro NPR steel anchor cable has a good control effect on large deformation of surrounding rock at intersection,which provides a new support material and support means for the safety and stabilization control of surrounding rock at intersection.

Frost heaving force considering synchronous damage to tunnel lining and surrounding rocks under freeze—thaw cycles

In areas with seasonal freezing,when the tunnel lining concrete is saturated with water infiltrating the interior,the lining and the surrounding rocks will simultaneously freeze.However,the current calculation of the frost heaving force fails to consider the synchronous damage to the lining and surrounding rocks under freeze-thaw cycles.Therefore,as per the elastic calculation model of the frost heaving force and model of steady-state heat transfer of circular tunnels,this study introduces the frost heaving rate of lining and surrounding rocks.First,the analytical solution of frost heaving force is obtained for simultaneous frost heaving of lining and surrounding rocks under any steady-state temperature field.Then,based on the fracture theory and meso-damage mechanics,the damage variables of lining and surrounding rocks under freeze-thaw cycles are extracted,representing their elastic modulus and porosity.Finally,the formula of frost heaving force for synchronous damage to the lining and surrounding rocks at any steady-state temperature field is obtained.The calculation results demonstrate that the lower the temperature inside the lining,the greater the frost heaving force.With the increasing number of freeze-thaw cycles,frost heaving force tends to gradually increase initially,reaching a peak value at 85 freeze-thaw cycles,decreasing to 80%of the peak value at 140 cycles before reaching a constant value.The lining participates in frost heaving,increasing the frost heaving force.The initial increase rate of frost heaving force is 15.7%.Changing the fitting coefficients s1 and s2 of the lining and surrounding rocks can effectively control the magnitude of the frost heaving force in the tunnels.

Heterogeneity identification method for surrounding rock of large-section rock tunnel faces based on support vector machine

Purpose–The purpose of the study is to quickly identify significant heterogeneity of surrounding rock of tunnel face that generally occurs during the construction of large-section rock tunnels of high-speed railways.Design/methodology/approach–Relying on the support vector machine(SVM)-based classification model,the nominal classification of blastholes and nominal zoning and classification terms were used to demonstrate the heterogeneity identification method for the surrounding rock of tunnel face,and the identification calculation was carried out for the five test tunnels.Then,the suggestions for local optimization of the support structures of large-section rock tunnels were put forward.Findings–The results show that compared with the two classification models based on neural networks,the SVM-based classification model has a higher classification accuracy when the sample size is small,and the average accuracy can reach 87.9%.After the samples are replaced,the SVM-based classification model can still reach the same accuracy,whose generalization ability is stronger.Originality/value–By applying the identification method described in this paper,the significant heterogeneity characteristics of the surrounding rock in the process of two times of blasting were identified,and the identification results are basically consistent with the actual situation of the tunnel face at the end of blasting,and can provide a basis for local optimization of support parameters.

Rock burst risk evaluation based on equivalent surrounding rock strength

On-site investigations consistently show that the rock burst inherent to coal seams varies greatly with coal seam thickness.In this study,impact factors related to coal seam thickness and surrounding rock strength were analyzed and a corresponding rock burst risk assessment method was constructed.The model reflects the influence of coal seam thickness on the stress distribution of surrounding rock at the roadway.Based on the roadway excavation range,a stress distribution model of surrounding roadway rock is established and the influence of coal seam thickness on rock burst risk is analyzed accordingly.The proposed rock burst risk assessment method is based on the equivalent surrounding rock strength and coal seam bursting liability.The proposed method was tested in a 3500 mining area to find that it yields rock burst risk assessment results as per coal seam thickness that are in accordance with real-world conditions.The results presented here suggest that coal seam thickness is a crucial factor in effective rock burst risk assessment.

Stress distribution and surrounding rock control of mining near to the overlying coal pillar in the working face

The occurrence of overlying coal pillar(OCP)exerts a strong effect on the stress and strain distribution of the surrounding rock in the stope.In this paper,the stress distribution characteristics are analyzed via the numerical calculation with the account of OCP presence or absence.In addition,this study revealed the joint effect of side pressure relief area of the goaf and stress concentration in OCP on the final stress distribution.Furthermore,the rules of abutment stress distribution affected by three influencing factors,namely horizontal-vertical distances between OCP and working face and buried depth of OCP,are analyzed.The functional model linking the peak stress of surrounding rock with the above influencing factors is developed.The field application of the above results proved that the rib spalling and deformation of a 2.95 m-high and 5.66 m-wide roadway could be efficiently controlled by rationally adjusting working states of the support,and adopting the hydraulic prop coordinated with the p type metal beam and anchor cable to strengthen the surrounding rock of working face and roadway,respectively.The proposed measures are considered appropriate to satisfy the safe operation requirements.

Numerical analysis of deformation and failure characteristics of deep roadway surrounding rock under static-dynamic coupling stress

In actual production,deep coal mine roadways are often under typical static-dynamic coupling stress(SDCS)conditions with high ground stress and strong dynamic disturbances.With the increasing number of disasters and accidents induced by SDCS conditions,the safe and efficient production of coal mines is seriously threatened.Therefore,it is of great practical significance to study the deformation and failure characteristics of the roadway surrounding rock under SDCS.In this paper,the effects of different in-situ stress fields and dynamic load conditions on the surrounding rock are studied by numerical simulations,and the deformation and failure characteristics are obtained.According to the simulation results,the horizontal stress,vertical stress and dynamic disturbance have a positive correlation with the plastic failure of the surrounding rock.Among these factors,the influence of the dynamic disturbance is the most substantial.Under the same stress conditions,the extents of deformation and plastic failure of the roof and ribs are always greater than those of the floor.The effect of horizontal stresses on the roadway deformation is more notable than that of vertical stresses.The results indicate that for the roadway under high-stress conditions,the in-situ stress test must be strengthened first.After determining the magnitude of the in-situ stress,the location of the roadway should be reasonably arranged in the design to optimize the mining sequence.For roadways that are strongly disturbed by dynamic loads,rock supports(rebar/cable bolts,steel set etc.)that are capable of maintaining their effectiveness without failure after certain dynamic loads are required.The results of this study contribute to understanding the characteristics of the roadway deformation and failure under SDCS,and can be used to provide a basis for the support design and optimization under similar geological and geotechnical circumstances.

Macro/meso failure behavior of surrounding rock in deep roadway and its control technology

Engineering disasters occur frequently and violently with the increase in mining depth, which is mostly due to insufficient study on the failure mechanism of the deep rock mass. In this paper, theoretical and experimental researches on the failure behaviors and deformation control of deep surrounding rock in recent years were reviewed. Macro/meso failure mechanism of deep rock or coal-rock combined body under different loading conditions have been systematically investigated. Stress gradient failure theory of surrounding rock, uniform strength support in the deep roadway, and the analogous hyperbola movement model of overlying strata were preliminary established and a combined grouting control technology for surface and underground was proposed. Abovementioned achievements are expected to offer theoretical bases and technical supports for the exploitation of China's deep mineral resources in the future.

Domino instability effect of surrounding rock-coal pillars in a room-and-pillar gob

To discuss the domino instability effect and large area roof falling and roof accidents of surrounding rockcoal pillars in a room-and-pillar gob,the equilibrium equation for a roof-coal pillar-floor system with the influence of mining floor was developed based on the engineering conditions of the surrounding rock in a room-and-pillar gob in the 3^(-2)coal seam of Tanggonggou mine.The conditions of system instability and the relationship between system stability and system stiffness were analyzed from an energetic point of view.Numerical simulation using the discrete element software UDEC was also carried out to simulate conditions causing the domino effect on surrounding rock-coal pillars in a 3^(-2)room-and-pillar gob.The results show that:if we want the system to destabilize,the collective energy in roof-and-floor must be larger than that in the coal pillar.When the stiffness of the coal pillars and the roof-and-floor are both greater than zero,the system is stable.When the stiffness of the coal pillars is negative but the summed stiffness of the coal pillars and roof-and-floor is larger than or equal to zero,the system is statically destroyed.When the sum of the coal pillars and the roof-floor stiffness is negative,the system suffers from severe damages.For equal advance distances of the coal mining face,the wider coal pillars can reduce the probability of domino type instability.Conversely,the smaller width pillars can increase the instability probability.Domino type instability of surrounding rock-coal pillars is predicted to be unlikely when the width of coal pillars is not less than 8 m.

Surrounding rock control of gob-side entry driving with narrow coal pillar and roadway side sealing technology in Yangliu Coal Mine

Gob-side entry driving can increase coal recovery ratio, and it is implied in many coal mines. Based on geological condition of 10416 working face tailentry in Yangliu Coal Mine, the surrounding rock deformation characteristics of gob-side entry driving with narrow coal pillar is analysed, reasonable size of coal pillar and reasonable roadway excavation time after mining are achieved. Surrounding rock control technology and effective roadway side sealing technology are proposed and are taken into field practice. The results showed that a safer and more efficient mining of working face can be achieved. In addition, results of this paper also have important theoretical significance and valuable reference for surrounding rock control technology of gob-side entry driving with narrow coal pillar under special geological condition.