Anisotropic shearing mechanism of Kangding slate:Experimental investigation and numerical analysis

The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.

Effects of seepage pressure on the mechanical behaviors and microstructure of sandstone

Surrounding rocks of underground engineering are subjected to long-term seepage pressure,which can deteriorate the mechanical properties and cause serious disasters.In order to understand the impact of seepage pressure on the mechanical property of sandstone,uniaxial compression tests,P-wave velocity measurements,and nuclear magnetic resonance(NMR)tests were conducted on saturated sandstone samples with varied seepage pressures(i.e.0 MPa,3 MPa,4 MPa,5 MPa,6 MPa,7 MPa).The results demonstrate that the mechanical parameters(uniaxial compressive strength,peak strain,elastic modulus,and brittleness index),total energy,elastic strain energy,as well as elastic strain energy ratio,decrease with increasing seepage pressure,while the dissipation energy and dissipation energy ratio increase.Moreover,as seepage pressure increases,the micro-pores gradually transform into meso-pores and macro-pores.This increases the cumulative porosity of sandstone and decreases P-wave velocity.The numerical results indicate that as seepage pressure rises,the number of tensile cracks increases progressively,the angle range of microcracks is basically from 50-120to 80-100,and as a result,the failure mode transforms to the tensile-shear mixed failure mode.Finally,the effects of seepage pressure on mechanical properties were discussed.The results show that decrease in the effective stress and cohesion under the action of seepage pressure could lead to deterioration of strength behaviors of sandstone.

A coupled cryogenic thermo-hydro-mechanical model for frozen medium:Theory and implementation in FDEM

This paper presents the development of a coupled modeling approach to simulate cryogenic thermo-hydro-mechanical(THM)processes associated with a freezing medium,which is then implemented in the combined finite-discrete element method code(FDEM)for multi-physics simulation.The governing equations are deduced based on energy and mass conservation,and static equilibrium equations,considering water/ice phase change,where the strong couplings between multi-fields are supplemented by critical coupling parameters(e.g.unfrozen water content,permeability,and thermal conductivity).The proposed model is validated against laboratory and field experiments.Results show that the cryogenic THM model can well predict the evolution of strongly coupled processes observed in frozen media(e.g.heat transfer,water migration,and frost heave deformation),while also capturing,as emergent properties of the model,important phenomena(e.g.latent heat,cryogenic suction,ice expansion and distinct three-zone distribution)caused by water/ice phase change at laboratory and field scales,which are difficult to be all revealed by existing THM models.The novel modeling framework presents a gateway to further understanding and predicting the multi-physical coupling behavior of frozen media in cold regions.

Experimental study on the influences of cutter geometry and material on scraper wear during shield TBM tunnelling in abrasive sandy ground

When shield TBM tunnelling in abrasive sandy ground,the rational design of cutter parameters is critical to reduce tool wear and improve tunnelling efficiency.However,the influence mechanism of cutter parameters on scraper wear remains unclear due to the lack of a reliable test method.Geometry and material optimisation are often based on subjective experience,which is unfavourable for improving scraper geological adaptability.In the present study,the newly developed WHU-SAT soil abrasion test was used to evaluate the variation in scraper wear with cutter geometry,material and hardness.The influence mechanism of cutter parameters on scraper wear has been revealed according to the scratch characteristics of the scraper surface.Cutter geometry and material parameters have been optimised to reduce scraper wear.The results indicate that the variation in scraper wear with cutter geometry is related to the cutting resistance,frictional resistance and stress distribution.An appropriate increase in the front angle(or back angle)reduces the cutting resistance(or frictional resistance),while an excessive increase in the front angle(or back angle)reduces the edge angle and causes stress concentration.The optimal front angle,back angle and edge angle for quartz sand samples areα=25°,β=10°andγ=55°,respectively.The wear resistance of the modelled scrapers made of different metal materials is related to the chemical elements and microstructure.The wear resistances of the modelled scrapers made of 45#,06Cr19Ni10,42CrMo4 and 40CrNiMoA are 0.569,0.661,0.691 and 0.728 times those made of WC-Co,respectively.When the alloy hardness is less than 47 HRC(or greater than 58 HRC),scraper wear decreases slowly with increasing alloy hardness as the scratch depth of the particle asperity on the metal surface stabilizes at a high(or low)level.However,when the alloy hardness is between 47 HRC and 58 HRC,scraper wear decreases rapidly with increasing alloy hardness as the scratch depth transitions from high to low levels.The sensitive hardness interval and recommended hardness interval for quartz sand are[47,58]and[58,62],respectively.The present study provides a reference for optimising scraper parameters and improving cutterhead adaptability in abrasive sandy ground tunnelling.

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.

Experimental and numerical investigation on failure characteristics and mechanism of coal with different water contents

Water injection, as a widely used technique to prevent coal burst, can restrain the fractured coal seam and released the energy storage. In this study, laboratory tests were frstly carried out on standard coal specimens with fve diferent water contents (i.e., 0%, 0.6%, 1.08%, 1.5%, 2.0%, and 2.3%). The failure mode, fragment size, and energy distribution characteristics of coal specimens were investigated. Experimental results show that strength, elastic strain energy, dissipated energy, brittleness index, as well as impact energy index decrease with increasing water content. Besides, the failure mode transitions gradually from splitting ejection to tensile-shear mixed failure mode as water content increases, and average fragment size shows positively related to water content. Moreover, scanning electron microscope tests results indicate that water in coal sample mainly causes the mineral softening and defects increase. Furthermore, a numerical model containing roadway excavation was established considering the water on coal burst prevention. Modelling results revealed that water injection can reduce degree of coal burst and ejection velocity of coal blocks, while it will raise up the depth of crack zone and surface displacement of roadway. Combined with laboratory tests and numerical results, the micro mechanism, energy mechanism, and engineering signifcance of water injection on coal burst prevention were fnally analyzed.

Role of intrinsic defects on carbon adsorbent for enhanced removal of Hg^(2+)in aqueous solution

Carbon is a normally used adsorbent for removal of heavy metal ion in aqueous solutions,but the efficient adsorbent needs intensive modification by heteroatom doped or supported noble metals that cause severe pollution and easy leaching of active components during use.In this paper,the role of intrinsic defects on Hg^(2+)adsorption for carbon adsorbent was investigated.The maximum adsorbing capacity of defectrich carbon has been improved up to 433 mg·g^(-1)which is comparable to most of the modified carbon adsorbents via supported metal chloride or noble metal components.The basicity is increased with the content of defective sites and the strong chemical bonding can be formed via electron transformation between the defect sites with adsorbed Hg^(2+).The present study gives a direction to explore cheap and easily scale-up high-performance mercury adsorbents by simply tuning the intrinsic defective structure of carbon without the necessity to support metal or other organic compounds.

An efficient 3D cell-based discrete fracture-matrix flow model for digitally captured fracture networks

Complex hydraulic fracture networks are critical for enhancing permeability in unconventional reservoirs and mining indus-tries.However,accurately simulating the fluid flow in realistic fracture networks(compared to the statistical fracture net-works)is still challenging due to the fracture complexity and computational burden.This work proposes a simple yet efficient numerical framework for the flow simulation in fractured porous media obtained by 3D high-resolution images,aiming at both computational accuracy and efficiency.The fractured rock with complex fracture geometries is numerically constructed with a cell-based discrete fracture-matrix model(DFM)having implicit fracture apertures.The flow in the complex fractured porous media(including matrix flow,fracture flow,as well as exchange flow)is simulated with a pipe-based cell-centered finite volume method.The performance of this model is validated against analytical/numerical solutions.Then a lab-scale true triaxial hydraulically fractured shale sample is reconstructed,and the fluid flow in this realistic fracture network is simu-lated.Results suggest that the proposed method achieves a good balance between computational efficiency and accuracy.The complex fracture networks control the fluid flow process,and the opened natural fractures behave as primary fuid pathways.Heterogeneous and anisotropic features of fluid flow are well captured with the present model.

Real-time rock mass condition prediction with TBM tunneling big data using a novel rock-machine mutual feedback perception method

Real-time perception of rock mass information is of great importance to efficient tunneling and hazard prevention in tunnel boring machines(TBMs).In this study,a TBM-rock mutual feedback perception method based on data mining(DM) is proposed,which takes 10 tunneling parameters related to surrounding rock conditions as input features.For implementation,first,the database of TBM tunneling parameters was established,in which 10,807 tunneling cycles from the Songhua River water conveyance tunnel were accommodated.Then,the spectral clustering(SC) algorithm based on graph theory was introduced to cluster the TBM tunneling data.According to the clustering results and rock mass boreability index,the rock mass conditions were classified into four classes,and the reasonable distribution intervals of the main tunneling parameters corresponding to each class were presented.Meanwhile,based on the deep neural network(DNN),the real-time prediction model regarding different rock conditions was established.Finally,the rationality and adaptability of the proposed method were validated via analyzing the tunneling specific energy,feature importance,and training dataset size.The proposed TBM-rock mutual feedback perception method enables the automatic identification of rock mass conditions and the dynamic adjustment of tunneling parameters during TBM driving.Furthermore,in terms of the prediction performance,the method can predict the rock mass conditions ahead of the tunnel face in real time more accurately than the traditional machine learning prediction methods.

A real-time prediction method for tunnel boring machine cutter-head torque using bidirectional long short-term memory networks optimized by multi-algorithm

Based on data from the Jilin Water Diversion Tunnels from the Songhua River(China),an improved and real-time prediction method optimized by multi-algorithm for tunnel boring machine(TBM)cutter-head torque is presented.Firstly,a function excluding invalid and abnormal data is established to distinguish TBM operating state,and a feature selection method based on the SelectKBest algorithm is proposed.Accordingly,ten features that are most closely related to the cutter-head torque are selected as input variables,which,in descending order of influence,include the sum of motor torque,cutter-head power,sum of motor power,sum of motor current,advance rate,cutter-head pressure,total thrust force,penetration rate,cutter-head rotational velocity,and field penetration index.Secondly,a real-time cutterhead torque prediction model’s structure is developed,based on the bidirectional long short-term memory(BLSTM)network integrating the dropout algorithm to prevent overfitting.Then,an algorithm to optimize hyperparameters of model based on Bayesian and cross-validation is proposed.Early stopping and checkpoint algorithms are integrated to optimize the training process.Finally,a BLSTMbased real-time cutter-head torque prediction model is developed,which fully utilizes the previous time-series tunneling information.The mean absolute percentage error(MAPE)of the model in the verification section is 7.3%,implying that the presented model is suitable for real-time cutter-head torque prediction.Furthermore,an incremental learning method based on the above base model is introduced to improve the adaptability of the model during the TBM tunneling.Comparison of the prediction performance between the base and incremental learning models in the same tunneling section shows that:(1)the MAPE of the predicted results of the BLSTM-based real-time cutter-head torque prediction model remains below 10%,and both the coefficient of determination(R^(2))and correlation coefficient(r)between measured and predicted values exceed 0.95;and(2)the incremental learning method is suitable for realtime cutter-head torque prediction and can effectively improve the prediction accuracy and generalization capacity of the model during the excavation process.

Impact of stress on solute transport in a fracture network: A comparison study

This paper compares numerical modeling of the effect of stress on solute transport （advection and matrix diffusion） in fractured rocks in which fracture apertures are correlated with fracture lengths. It is mainly motivated by the performance and safety assessments of underground radioactive waste repositories. Five research teams used different approaches to model stress/deformation, flow and transport pro- cesses, based on either discrete fracture network or equivalent continuum models. The simulation results derived by various teams generally demonstrated that rock stresses could significantly influence solute transport processes through stress-induced changes in fracture apertures and associated changes in per- meability. Reasonably good agreement was achieved regarding advection and matrix diffusion given the same fracture network, while some observed discrepancies could be explained by different mechanical or transport modeling approaches.

ASYMPTOTIC PROPERTIES OF A BRANCHING RANDOM WALK WITH A RANDOM ENVIRONMENT IN TIME

We consider a branching random walk in an independent and identically distributed random environment ξ=(ξn) indexed by the time. Let W be the limit of the martingale Wn=∫e^-txZn(dx)/Eξ∫e^-txZn(dx), with Zn denoting the counting measure of particles of generation n, and Eξ the conditional expectation given the environment ξ. We find necessary and sufficient conditions for the existence of quenched moments and weighted moments of W, when W is non-degenerate.

Failure criterion for soft rocks considering intermediate principal stress

The significant differences between hard rocks(more brittle)and soft rocks(more ductile)may suggest the use of different failure criteria.A strength criterion for soft rocks that includes intermediate principal stress was proposed.The new criterion includes two independent parameters:the uniaxial compressive strength(σ_(ci)),which can be obtained from common laboratory tests or indirectly estimated by alternative index tests in the laboratory or field;and f(joint),which is used to characterize the rock mass quality and can be easily estimated.The authors compared the predictive capabilities of the new criterion with other criteria using the database of soft rocks under two conditions:with and without triaxial data.For the estimation of triaxial and true-triaxial strengths,the new criterion generally produced a better fit.The proposed criterion is practical for an approximate first estimation of rock mass strength,even without triaxial data,as it balances accuracy(lower prediction error)and simplicity(fewer independent parameters).

Theoretical and numerical studies of crack initiation and propagation in rock masses under freezing pressure and far-field stress

Water-bearing rocks exposed to freezing temperature can be subjected to freezeethaw cycles leading tocrack initiation and propagation, which are the main causes of frost damage to rocks. Based on theGriffith theory of brittle fracture mechanics, the crack initiation criterion, propagation direction, andcrack length under freezing pressure and far-field stress are analyzed. Furthermore, a calculation methodis proposed for the stress intensity factor （SIF） of the crack tip under non-uniformly distributed freezingpressure. The formulae for the crack/fracture propagation direction and length of the wing crack underfreezing pressure are obtained, and the mechanism for coalescence of adjacent cracks is investigated.In addition, the necessary conditions for different coalescence modes of cracks are studied. Using thetopology theory, a new algorithm for frost crack propagation is proposed, which has the capability todefine the crack growth path and identify and update the cracked elements. A model that incorporatesmultiple cracks is built by ANSYS and then imported into FLAC3D. The SIFs are then calculated using aFISH procedure, and the growth path of the freezing cracks after several calculation steps is demonstratedusing the new algorithm. The proposed method can be applied to rocks containing fillings such asdetritus and slurry. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Using benzene carboxylic acids to prepare zirconium-based catalysts for the conversion of biomass-derived furfural

Benzene carboxylic acid (BCAs) are common and useful chemical blocks, which can be derived from the abundant low rank coals (LRCs) via oxidative degradation. In this work, we proposed a novel strategy to utilize BCAs as raw materials to prepare catalysts with transition metal zirconium, and the prepared catalysts were applied into the conversion of the renewable biomass resources. Typical model BCAs in the oxidative products of LRCs, including pyromellitic acid, trimesic acid (TMSA), trimellitic acid, and benzoic acid, were used as the block to construct the Zr- BCAs catalysts. The chemoselective conversion of furfural into furfuryl alcohol (FAL). an important reaction in the biomass conversion chain, is chosen to evaluate the activity of the catalysts. The preparation conditions of the catalysts and experiment factors during the reaction were systematically investigated. The prepared catalysts were characterized by SEM, TEM, XRD and TG-DTG. The results showed that the prepared catalysts were efficient for the conversion of furfural into FAL, among which Zr-TMSA gave the highest activity. Zr-TMSA could be recycled for ten times without obvious deactivation, indicating an excellent stability. The strategy proposed in this work may be beneficial for the value-added utilization of both LRCs and biomass resources.

Attenuation characteristics of impact-induced seismic wave in deep tunnels:An in situ investigation based on pendulum impact test

The radiated seismic energy is an important index for the intensity assessment of microseismic(MS)events and the early warning of dynamic disasters.However,the energy of MS signals is significantly attenuated due to the heterogeneity and viscous damping of rock media.Therefore,the study on attenuation characteristics of MS signals in underground engineering has practical significance for the accurately estimation of radiated seismic energy.Based on a pendulum impact test facility and MS monitoring system,an in situ investigation was carried out to explore attenuation characteristics at a deep tunnel.The results show that the seismic energy and peak particle velocity(PPV)attenuation are exponentially related to the propagation distance.The attenuation coefficient of energy is larger than that of PPV.With the increase in the input impact-energy,the seismic energy attenuation coefficient decreases as a power function.An empirical relationship between energy attenuation coefficient and wave impedance of rock mass was established in this scenario.Moreover,the time-frequency characteristics and energy distribution laws of impact-induced signals were investigated by the continuous wavelet transform(CWT)and wavelet packet analyses,respectively.The dominant frequency of signals decreases gradually as the propagation distance increases.Based on the energy attenuation characteristics,a new method was proposed to calculate the released source energy of MS events in the field.This study can provide an insight into energy attenuation characteristics of seismic waves and references for attenuation correction in seismic energy calculation.

A Numerical Study of Wave Propagation and Cracking Processes in Rock-Like Material under Seismic Loading Based on the Bonded-Particle Model Approach

An earthquake is usually followed by a considerable number of aftershocks that play a significant role in earthquake-induced landslides,During the aftershock,the cracking process in rocks becomes more complex because of the formation of faults.In order to investigate the effects of seismic loading on the cracking processes in a specimen containing a single flaw,a numerical approach based on the bonded-particle model(BPM)was adopted to study the seismic loading applied in two orthogonal directions.The results reveal that no transmission and reflection phenomena were observable in the small specimens(76 mm×152 mm)because they were considerably smaller than the wavelength of the P-wave.Furthermore,under seismic loading,the induced crack was solely tensile in nature.Repeated axial seismic loading did not induce crack propagation after the first axial seismic loading.Cracks began to propagate only when the seismic loading direction was changed from axial to lateral,and then back to axial,ultimately resulting in the failure of the specimen.

Transient Electro-Osmotic and Pressure Driven Flows through a Microannulus

Flow behavior of transient mixed electro-osmotic and pressure driven flows (EOF/PDF) through a microannulus is investigated based on a linearized Poisson-Boltzmann equation and Navier-Stokes equation. A semi-analytical solution of EOF velocity distribution as functions of relevant parameters is derived by Laplace transform method. By numerical computations of inverse Laplace transform, the effects of inner to outer wall zeta potential β, the normalized pressure gradient Ω and the inner to outer radius ratio α on transient EOF velocity are presented.

Evolvement behavior of microstructure and H_2O adsorption of lignite pyrolysis

The effect of pyrolysis on the microstructure and moisture adsorption of lignite was investigated with low field nuclear magnetic resonance spectroscopy. Changes in oxygen-containing groups were analyzed by Fourier transform infrared spectroscopy （FTIR）, and H20 adsorption mechanism on the surface of lignite pyrolysis was inferred. Two major changes in the pore structure of lignite char were observed as temperature increased in 105-200 ℃ and 500-700 ℃. Pyrolysis temperature is a significant factor in removing carboxyl and phenolic hydroxyl from lignite. Variation of ether bond content can be divided into three stages; the content initially increased, then decreased, and finally increased. The equilibrium adsorption ratio, content of oxygen-containing groups, and variation of pore volume below 700° were closely correlated with each other. The amount of adsorbed water on char pyrolyzed at 700 ℃ increased. Moreover, the adsorption capacity of the lignite decreased, and the adsorption state changed.

A thermo-hydro-mechano-chemical formulation for modeling water transport around a ventilated tunnel in an argillaceous rock

Argillaceous rocks are being considered as potential host rocks for deep geological disposal. For the research work in DECOVALEX-2011, 5 participant research teams performed simulations of a labora- tory drying test and a ventilation experiment for Mont Terri underground laboratory built in argillaceous rock formation. Our study starts with establishing a coupled thermo-hydro-mechano-chemical （THMC） processes model to simulate water transport in rock around the ventilated tunnel. Especially in this THMC formulation, a three-phase and two-constituent hydraulic model is introduced to simulate the processes which occur during tunnel ventilation, including desaturation/resaturation in the rock, phase change and air/rock interface, and to explore the Opalinus clay parameter set, It can be found that water content evolution is very sensitive to intrinsic permeability, relative permeability and capillary pressure in clay rock. Water loss from surrounding rock is sensitive to the change of permeability in clay which is induced by excavation damaged zone. Chemical solute transport in the rock near ventilation experiment tunnel is simulated based on the coupled THMC formulation. It can be estimated that chemical osmotic flow has little significance on water flow modeling. Comparisons between simulation results from 5 teams and experimental observations show good agreement. It increases the confidence in modeling and indi- cates that it is a good start for fully THMC understanding of the moisture transportation and mechanical behavior in argillaceous rock.