During the construction of underground caverns, the stability of deep underground cavern excavation, which affects the safety and sustainable development of such projects, is a hot issue. First, based on the mechanical properties of surrounding rock in deep tunnels, the strain-softening behavior, damage, and heterogeneity of rock masses are analyzed. Then, a strain-softening model of heterogeneous jointed rock mass that considers statistical damage (SSD) is developed and implemented through FLAC3D simulation software. Finally, the SSD is applied to a deep roadway in the Jinchuan mining area, and a comparative analysis of the computation results of the Mohr–Coulomb (MC) model and the strain-softening (SS) model are carried out. The numerical results are compared with the field-monitoring results, which show that the SSD model simulated the behavior of the surrounding rocks well. The results show that the deformations of the roof and floor are larger, which may serve as a reference for the support pattern of deep roadways.

The phenomenon of toppling deformation and failure is common in slopes with anti-dip structures, especially in soft metamorphic rock slopes. This paper aims to explore the instability mechanism of anti-dip layered soft metamorphic rock landslides. Taking the slope of a mining area in the southern Qinling Mountains of China as a geological prototype, a large-scale centrifuge model test and a numerical simulation based on the combined finite and discrete element method (FDEM) were performed. The deformation and failure process, failure mode, and failure path of the slope under rainfall and excavation conditions were simulated. The results show that both the physical centrifuge model test and the new numerical model test can simulate the instability process of anti-dip layered soft metamorphic rock slopes, and the phenomena simulated by the two methods are also very close. Rainfall mainly weakens the mechanical properties of rock, while the excavation at the slope toe mainly changes the stress field distribution and provides space for slope deformation, both of which accelerate the instability of the anti-dip soft metamorphic rock slope. The failure process of an anti-dip layered soft rock slope can be described as follows: bending of the rock layer–tensile fracture along the layer–flexural toppling and cracking perpendicular to the rock layer–extension and penetration of the tensile fracture surface–sliding and instability of the slope.

Gas flow mechanisms and apparent permeability are important factors for predicating gas production in shale reservoirs. In this study, an apparent permeability model for describing gas multiple flow mechanisms in nanopores is developed and incorporated into the COMSOL solver. In addition, a dynamic permeability equation is proposed to analyze the effects of matrix shrinkage and stress sensitivity. The results indicate that pore size enlargement increases gas seepage capacity of a shale reservoir. Compared to conventional reservoirs, the ratio of apparent permeability to Darcy permeability is higher by about 1–2 orders of magnitude in small pores (1–10 nm) and at low pressures (0–5 MPa) due to multiple flow mechanisms. Flow mechanisms mainly include surface diffusion, Knudsen diffusion, and skip flow. Its weight is affected by pore size, reservoir pressure, and temperature, especially pore size ranging from 1 nm to 5 nm and reservoir pressures below 5 MPa. The combined effects of matrix shrinkage and stress sensitivity induce nanopores closure. Therefore, permeability declines about 1 order of magnitude compare to initial apparent permeability. The results also show that permeability should be adjusted during gas production to ensure a better accuracy.

Mixing calculations have been widely applied to identify sources of groundwater recharge, but these calculations have assumed that the concentrations of end-members are well known. However, the end-members of water remain unclear and are not easily available in practical applications. To better determine end-members and mixing ratios, an end-member mixing analysis combining multivariate statistical methods was used on a large, complex water chemistry dataset collected from the Shashandao gold mine in China. Multivariate statistical methods, including principal component analysis (PCA) and hierarchical cluster analysis (HCA), were applied to determine the specific end-members (these two methods verified each other). On the basis of the identified end-members, a maximum likelihood method was then used to estimate the mixing ratios of the water sources. The combined method proposed in this study can help to identify more accurate end-members and deal with uncertainty in end-member concentrations, and it can also adjust the concentrations until the optimal mixing ratios for the calculation are obtained. This method can be a powerful tool for groundwater management and in predicting water inrush in mining operations.

Submarine mine water inrush has become a problem that must be urgently solved in coastal gold mining operations in Shandong, China. Research on water in subway systems introduced classifications for the types of mine groundwater and then established the functions used to identify each type of water sample. We analyzed 31 water samples from −375 m underground using multivariate statistical analysis methods. Cluster analysis combined with principle component analysis and factor analysis divided water samples into two types, with one type being near the F3 fault. Principal component analysis identified four principle components accounting for 91.79% of the total variation. These four principle components represented almost all the information about the water samples, which were then used as clustering variables. A Bayes model created by discriminant analysis demonstrated that water samples could also be divided into two types, which was consistent with the cluster analysis result. The type of water samples could be determined by placing Na⁺ and CHO₃⁻ concentrations of water samples into Bayes functions. The results demonstrated that F3, which is a regional fault and runs across the whole Xishan gold mine, may be the potential channel for water inrush, providing valuable information for predicting the possibility of water inrush and thus reducing the costs of the mining operation.

Cracks are critical for the deformation and failure of rock masses, but the effects of real cracks are rarely considered when evaluating the stability and safety of practical engineering. This paper presents a study on the application of fractures in the Sanshandao Gold Mine. Field investigation and statistical analysis methods were adopted to obtain the distribution laws of the cracks. Laboratory tests, MATLAB programming, and simulation using the software, GDEM (Gdem Technology, Beijing, China, Co., Ltd.), were employed to study the mechanical behaviors of rock masses with real fractures after excavation. The main results are as follows: (1) Three sets of highly discrete cracks were developed in the study area. Their inclination and dip can be approximately considered to follow a Gaussian distribution or uniform distribution. They had close ties to the three faults developed in the mining area. (2) Compared with the model that did not consider cracks and the model processed by the equivalence idea, the surrounding rock deformation caused by excavation of the model that considered real cracks was larger than that of the former and smaller than that of the latter. However, its influence range was larger than that of the other two models. The results show that it is reasonable to use three sets of discrete cracks to characterize the fracture distribution of the surrounding rock. In the evaluation of roadway stability, it is not advisable to use the equivalence method to deal with all the cracks. Considering a part of the cracks that are compatible with the size of the calculation model, a relatively accurate evaluation can be obtained in terms of the deformation, failure, and permeability changes of the surrounding rock.