20 December 2025, Volume 21 Issue 6
    

  • Select all
    |
  • Wei Lingxiang, Guo Dongjun, Wu Yanhua, Liang Sufang, Chen Zhilong
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1851-1860. https://doi.org/10.20174/j.JUSE.2025.06.01
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The utilization scale of underground space in China will continue to grow for a long period of time in the future, which can provide sufficient and stable space and facility support for a variety of emission reduction and sink enhancement technologies, and has great significance to the achievement of the goal of carbon peaking and carbon neutrality. In order to reveal the research hotspots and trend of China's underground space in response to the goal of carbon peaking and carbon neutrality, 944 academic papers published in the CNKI database before Jan., 2024 were used as data samples, and a visual analysis was conducted based on the VOSviewer software to analyze the number of annual publications, the distribution of journals, keywords, and institutions of publications in this field. The results of the study show that the research literature on the response of China's underground space to the goal of carbon peaking and carbon neutrality has shown fluctuating growth, with a significant increase in the number of literature in the past two years; the research content mainly focuses on the low-carbon effect of the subway, carbon utilization and storage in underground, and underground energy storage, etc. The research trend of China's underground space in response to the goal of carbon peaking and carbon neutrality is analyzed from four aspects, namely, spatial compactness, energy diversity, functional integration, and construction phasing.
  • Yu Hongbin, Chen Zhilong, Zhang Jian, Zhang Mengyu
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1861-1873. https://doi.org/10.20174/j.JUSE.2025.06.02
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    At present, many integration construction problems such as the disconnection of integrated development of urban rail stations and surrounding areas, insufficient integration of stations and cities, insufficient accessibility, and low comprehensive benefits are becoming more and more prominent, and at the same time, the mechanism for supervising and implementing the assessment of urban rail stations is not yet perfect. Based on the framework of urban rail station assessment system, this paper established an integrated assessment index system for urban rail stations from eight dimensions, including land function, traffic organisation, spatial environment, commercial vitality, economic benefits, safety resilience, underground space and implementation effectiveness, and took 32 rail stations in Beijing as the objects of study, and carried out the assessment and classification of the rail stations in different radial scales by using entropy Weight-TOPSIS method. Evaluation and classification of different radiation scales of rail stations were carried out by using entropy Weight-TOPSIS method. The results show that there are spatial differences in the influencing factors of rail station integration construction in different radiation scales of rail stations, and the assessment results are divided into four categories by using K-means clustering, and classification optimisation suggestions are put forward.
  • Li Xiaodan, Li Jing, Liu Zhiping, Song Yaobing, Hao Yifan
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1874-1883. https://doi.org/10.20174/j.JUSE.2025.06.03
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Urban land scarcity and increasing student numbers have limited the development of old university campuses. This study proposes a new concept of integrating above-ground and underground space for the study of the renovation design of old university campuses and applies this new concept to the Xueyuan Road Campus of China University of Mining and Technology (Beijing). The study was carried out with a focus on four aspects: functional zoning integration, traffic optimization, emergency evacuation of underground space, and space quality improvement. Adjustment of the functionality and zoning of the campus by analyzing the demolition of some of the existing old buildings, taking into account the surrounding environment and the above and below-ground conditions of the campus after the demolition, to separate pedestrian and motor vehicle traffic on campus. In this design, the large-scale gymnasium is underground beneath the playground where the upper space cannot be opened. An innovative "functional area-secondary safe area-safe area (sunken space)" evacuation system is proposed, and its feasibility was verified using the simulation software. Sunken corridors, sunken courtyards, and skylights were used to solve the problems of lighting in underground spaces. The methodology described in this paper effectively increased the area of various buildings on the campus and significantly improved the quality of the campus space. Our results provide a plan for the renovation of the Xueyuan Road Campus of China University of Mining and Technology (Beijing), and provide an effective method for the renovation of old university campuses.
  • Zheng Yingren, Abi Erdi, Yang Yingming
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1884-1887. https://doi.org/10.20174/j.JUSE.2025.06.04
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The traditional algorithm is unable to accurately calculate the ultimate bearing capacity of heavy soil foundation, whereas numerical limit analysis methods can do so precisely. It can be seen that numerical limit analysis can solve geotechnical engineering problems more complex than traditional analytical methods. The incremental load method is used, which applies the principle of failure of engineering materials from quantitative to qualitative changes in stress. As the foundation load gradually increases, the stress on the soil foundation transitions from elastic to plastic, and finally reaches the ultimate failure state, obtaining the corresponding ultimate load value. The ultimate bearing capacity of the foundation under strip foundation load was obtained by combining engineering examples, and compared with the traditional empirical results of Terzaghi, Vesic, and Chen. The results show that the numerical limit analysis method can save 40%, 20%, and similar foundation engineering quantities compared to Terzaghi, Vesic, and Chen empirical methods, respectively. At present, both Terzaghi formula and Vesic formula are applied in practical projects, while Chen formula has not been applied yet.
  • Shan Yao, Dong Yacheng, Wu Yaojie, Chen Zhining, Yao Xiping
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1888-1898. https://doi.org/10.20174/j.JUSE.2025.06.05
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To elucidate the variations in non-limit state earth pressure for limited soils mass under the rotation at the top of a rigid retaining wall (RT) mode, and to address the incomplete consideration of wall displacement, horizontal shear stress, and other influencing factors in existing earth pressure calculation methods, a numerical calculation method for passive earth pressure of a finite soil mass was developed within the framework of the differential element method. This method integrates considerations of wall displacement, soil strength parameters (cm-φm), soil arching effects, and the impact of horizontal shear stress. By combining triaxial unloading tests and stress Mohr circles, the relationship between soil strength parameters (cm-φm) under non-limit states and wall displacement was derived. Based on the principles of soil arching effect, stress analysis of the finite soil body was conducted, and static equilibrium equations for differential soil elements were established, further creating a numerical iterative format for passive earth pressure within a depth range. The proposed method was validated for its rationality and accuracy through comparisons with existing research. The results show that the lateral constraints under the RT mode help enhance the soil's shear resistance; Non-linear distribution of passive earth pressures along depth in finite soils; the magnitude of earth pressure decreases with a reduction in soil strength parameters (cm-φm), a decrease in the initial slip surface angle, a reduction in the bottom of wall displacement ratio, and an increase in the width-to-depth ratio. The method of this paper can provide theoretical guidance for the economic design of retaining structures in practical engineering.
  • Li Tiantian, Ha Da, Wang Shilei
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1899-1906. https://doi.org/10.20174/j.JUSE.2025.06.06
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In practice engineering, it is commonly using pumping tests to determine hydraulic parameters, but the previous of estimated parameters cannot be proved due to the lack of unified standard and hard to choose a suitable algorithm. Meanwhile, due the limitation of independent algorithm, in some difficult hydraulic models, the engineers have to adjust parameters continually to meet the convergence criterion which cost lots of time and the accuracy of results cannot be proved. In this study, a series of pumping tests has been conducted in Dongli district of Tianjin city, China. With considering storage effect in the pumping wells and observation wells which has important impact in the early experiment, the effective ration has been proposed to eliminate the effect of well storage. Based on the local hydrogeological environment and measured drawdown, an unsteady leaky seepage model, called Hantush-Jacob model, is used to estimate hydraulic parameters. By combination of measured drawdown and analytical solution,three methods, extend Kalman filter algorithm, LM algorithm and Genetic algorithm, are all adopted to estimate hydraulic parameters in this study. According to the result, three algorithm all shows its limitation in the calculation. Thus, a combination of Genetic algorithm and LM algorithm has been proposed in this study. Based on the result, the combination algorithm shows its superiority in searching the optimal solution accurately and avoiding the defect of the three separate algorithm. According to the comparison, the measured drawdown fit well with the calculated drawdown based on the estimated hydraulic parameter. Furthermore, compared with the single algorithms used in previously published literature, the combined GA-LM algorithm has demonstrated its unique advantages in the inversion of hydrological parameters.
  • Liu Zhongxian, Wei Xia, Jin Liguo, Zhou Tao, Wang Ying
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1907-1915. https://doi.org/10.20174/j.JUSE.2025.06.07
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Discussing the development of the indirect boundary element method to the simulation of seismic dynamic interaction between river valleys and tunnels under Rayleigh wave incidence. On the basis of verifying the accuracy of the proposed method, the influences of tunnel buried depth and valley width-depth ratio on the surface displacement amplitude and circumferential stress of the lining tunnel under the action of the Rayleigh wave with different frequencies are quantitatively analyzed. The results show that: With the increase of tunnel depth, the displacement above the tunnel gradually decreases, and the amplification effect of surface displacement gradually weakens. The influence of the ratio of width to depth of the valley on the displacement amplitude in the valley and on the right side of the valley is more significant. Under shallow burial conditions, the circumferential stress amplitude coefficient of the tunnel near the river valley can reach more than 10.0, showing a significant dynamic stress amplification effect. Under η=2.0 incidence, the amplitude-amplitude-amplitude-amplitude-is significantly affected by the valley width-to-depth ratio, but the amplitude-amplitude-effect is not obvious when the buried depth d=3.0a. The seismic design of the tunnel along the river should focus on the amplitude effect of tunnel stress under shallow buried low frequency.
  • Zhao Guojun, Guo Chengchao, Wang Fuming, Wang Haibo
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1916-1926. https://doi.org/10.20174/j.JUSE.2025.06.08
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to study the crack-resistance performance of fiber-reinforced composite (FRP) reinforced tunnel cracked lining, based on the peridynamics theory, a peridynamics simulation method for FRP reinforced concrete structures is established by introducing the interface weakening criteria. Four-point bending tests of FRP reinforced concrete beams are conducted to verify the simulation method. The problem of FRP reinforced tunnel cracked lining is modeled and analyzed to explore the influence of FRP reinforcement width and layers on the reinforcement performance of various cracked linings. The results indicate that: Interfacial debonding is the main failure mode of FRP reinforced cracked lining. The width and number of layers of FRP reinforcement are positively correlated with the reinforcement performance. Nevertheless, when the reinforcement width exceeds 2/3 of the arch bottom and the number of layers exceeds 2, the structural reinforcement performance does not significantly improve. In addition, the distribution and depth of cracks have a significant impact on the reinforcement performance of FRP, and the reinforcement effect of mid-span cracks is significantly better than that of eccentric cracks. Moreover, when the height ratio of mid span cracks is 1/5, the global damage index of FRP reinforcement lining decreases by 40.3%, significantly improving the crack resistance performance of the structure.
  • Zhang Yueming, Wu Shuguang, Huang Kai, Lan Yuanxin, Yang Kaicheng
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1927-1934. https://doi.org/10.20174/j.JUSE.2025.06.09
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Gravel piles are widely used in engineering, and their bearing capacity is an important indicator for evaluating the performance of gravel pile projects. Based on the cavity expansion theory, the calculation formula of bearing capacity of single gravel pile is derived by introducing the extended Lade-Duncan failure criterion, and considering the superimposed effect of stress expansion, the calculation formula of bearing capacity of single gravel pile under square layout scheme is derived. The laboratory scale model test shows that the maximum swelling depth of stone column is about twice the depth of pile diameter, which provides a basis for theoretical calculation. Finally, the proposed calculation method is compared and verified by field tests and existing research results. The results show that: The calculation results of the existing bearing capacity calculation formula of the crushed stone column are conservative, and the bearing capacity obtained by introducing the calculation formula of the extended Lade-Duncan failure criterion is closer to the measured value. This study can provide reference for related engineering and research.
  • Guo Nan, Jia Shining, Yang Xiaohui, Chen Zhenghan
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1935-1945. https://doi.org/10.20174/j.JUSE.2025.06.10
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to further investigate the deformation and strength characteristics of transversely isotropic unsaturated remolded loess under complex stress paths, three sets of 39 consolidation shear tests were conducted on transversely isotropic unsaturated soil using different experimental instruments (unsaturated true triaxial apparatus and unsaturated triaxial apparatus). The results show that: In true triaxial shear tests, the stress-strain relationship curve of the specimen exhibits strain hardening characteristics and has a hyperbolic shape, and the stress also increases with the increase of suction force.As the shearing process continues, the moisture content of the true triaxial specimen decreases continuously. Under the same b value and net confining pressure conditions, the greater the suction force, the greater the slope of the moisture content axial strain relationship curve. The true triaxial specimen, K0CD test specimen, and K0′CD test specimen consistently exhibit a state of shear shrinkage during the shear test. The magnitude of the failure stress of the test specimens under three different initial conditions is highest for the true triaxial test specimen, followed by the K0CD test specimen, and lowest for the K0′CD test specimen. The initial tangent modulus of the K0CD test specimen is the highest, followed by the K0′CD test specimen, and the true triaxial test specimen is the smallest. The research results are of great significance for promoting the study of mechanical properties of transversely isotropic unsaturated loess, improving relevant theories and criteria, and guiding engineering practice.
  • Fang Tao, Xu Yinan, Zhou Peng, Fang Wenyuan, Xu Jialei
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1946-1955. https://doi.org/10.20174/j.JUSE.2025.06.11
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The space around the proposed excavation in complex environments becomes increasingly limited, the semi-infinite space assumed by the classical earth pressure theory is inconsistent with actual engineering, resulting in a lack of reasonable calculation methods for active earth pressure in excavations with limited soil. Taking the flexible retaining piles and sandy soil of limited soil excavations as the research objects, using the single-row pile cantilever retaining mode and the bottom ends of the retaining piles fixed, four excavations model tests under different soil widths were carried out. To explore the distribution characteristics of active earth pressure and horizontal displacement of retaining piles along the retaining piles. The results show that the horizontal displacement of the retaining structure on each soil width side approximates a triangular displacement pattern, and the earth pressure gradually approximates a D-shaped distribution as the excavation depth increases, with the maximum value appearing at approximately 0.8H depth. The force analysis of the soil differential unit was performed through the limit equilibrium method, and the active earth pressure distribution of the finite soil under the limit state was deduced. Combined with the friction angle development model that considers the influence of displacement, the non-limit active earth pressure distribution and resultant force of the finite soil were deduced. and the height of the point of action. The effectiveness of the calculation method in this paper is verified by comparing the earth pressure test values with the theoretical values in this paper. Through parameter analysis, it was found that the earth pressure is positively related to the aspect ratio and negatively related to the internal friction angle. The critical aspect ratio ncr = 0.5 to distinguish between finite soil and semi-infinite soil.
  • Li Yanxia, Mou Chunmei, Yang Jin, Pan Qiang, Zhang Bo
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1956-1967. https://doi.org/10.20174/j.JUSE.2025.06.12
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The triaxial test based on photogrammetry can realize the real-time measurement of the local deformation of the specimen, but the refraction of light leads to the error in the measurement of specimen deformation. The introduction of the refraction correction factor (κ) can not only solve this problem, but also improve the test efficiency. In order to investigate the trend of the refraction correction coefficient in the axial and radial deformation measurements of different specimens, and to analyze the difference between the stress-strain curves measured by applying the refraction correction coefficient and those of the conventional triaxial tests, the photogrammetric-based consolidation and drainage triaxial tests were carried out on red clay and sandy clay, respectively. The results show that: (1) In the axial deformation measurement, the refraction correction coefficients of the different local regions of the red clay and the sandy clay converged to 1, and the refraction amplification effect on the axial deformation measurement of the specimen is negligible; (2) In the radial deformation measurement, the refraction correction coefficients of the red clay and sandy soil are 0.812 and 0.757, respectively. The refraction amplification effect has a greater impact on the radial deformation measurement of the specimen, and the refraction correction coefficients of the different specimens correspond to different refraction correction coefficients; (3) The change trends of the stress-strain curves measured by conventional triaxial and the application of refraction correction coefficient. The stress-strain curves measured by conventional triaxial and by applying the refraction correction coefficient have the same trend, and the bias stress measured by applying the refraction correction coefficient photogrammetry is relatively small in the middle and late stages of the test.
  • Zhong Xiaochun, Zhang Wenbin, Zhou Zhi, Huang Siyuan, Yun Qiang
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1968-1975. https://doi.org/10.20174/j.JUSE.2025.06.13
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Earth pressure balance shield tunneling is prone to spewing during excavation under conditions of large burial depth and high water pressure. Polymers are generally used as emergency modifiers to deal with spewing. Due to the high cost of polymers and the need for specialized mixing equipment on the ground to mix them into aqueous solutions, there are few cases of their application in soil conditioning. In order to explore the adaptability of water-soluble polymer in improving the formation of medium coarse sand, polymer soil conditioning was carried out for medium coarse sand with different fine-grained contents. Pressure permeability tests, critical spewing pressure tests, and flowability tests were carried out. The results show that: The permeability coefficient of improved soil is greatly affected by changes in fine particle content, while the effect of polymer addition on its permeability coefficient is relatively gentle; The flowability and critical surge pressure of the improved soil increase linearly with the increase of fine particle content; Exploring the adaptability of water-soluble polymers to the improvement of medium coarse sand residue from three aspects: permeability coefficient, critical spewing pressure, and flowability, and determining that polymers are suitable for medium coarse sand formations with a fine particle content of up to 4%; Comparing the economic benefits of the polymer with commonly used bentonite slurry, the results show that the polymer not only has a small impact on construction efficiency but also has certain cost advantages.
  • Bao Weixing, Wu Qian, Wu Qian, Qin Chuan, Hou Tianqi
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1976-1987. https://doi.org/10.20174/j.JUSE.2025.06.14
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to study the influence of different water content and salt content on the water-salt migration mechanism of Ili salinized loess. The water-salt migration test of saline loess under the action of the freeze-thaw cycle was carried out, and the temperature, moisture, and salt fields and their cohesion change rules were analyzed. The study shows that under the action of freezing and thawing, the temperature field, moisture field and salt field of saline loess show cyclic changes. With the increase of the number of freezing and thawing cycles, the water content of the soil column as a whole is increased, showing a pattern of high water content at the two ends and low water content in the middle. The change of salts in the soil layer under freezing and thawing is consistent with the overall change of moisture, and the salts in each soil layer show a homogeneous distribution before the unfreezing and thawing. Under the effect of temperature potential, the water carries the salts to migrate upward, resulting in higher salinity in the middle and upper part of the soil column. The increase of water content and salinity has a promoting effect on the migration of water and salts. In addition, the increase in water content will have a lubricating effect between the soil particles, resulting in the cohesion of the sample with low water content being greater than that of the sample with high water content. With the increase of salt content, the agglomerates formed by the salt crystals wrapped around the soil particles play the role of skeleton, which makes the cohesion of the specimen with high salt content larger than that of the specimen with low salt content, the value of internal friction angle of salinized loess after freeze-thawing shows a wave distribution.
  • Chen Lichao, Zhang Chaopeng, Zhang Diankun, Lü Shuaifeng, Wang Shengwei
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1988-1996. https://doi.org/10.20174/j.JUSE.2025.06.15
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to reveal the influence mechanism of acid-rock reaction in different bedding directions on fracture properties of Marine shale, SNDB samples were used to obtain the load-displacement relationship and fracture load of samples under three-point bending loading, and the fracture toughness(KIC) and fracture energy of shale were calculated. Combined with network fracturing requirements, acid injection modes of marine shale reservoirs under different depths and stress states were initially discussed. The results show that: (1) The fracture process of the sample under three-point loading includes the stages of compaction, elastic compression, fracture, and fracture expansion. The bearing capacity of the sample is rapidly lost after fracture, and the shale is brittle. (2) In parallel, vertical, and shear stratification directions, the fracture toughness of shale decreases with the increase of acid concentration, KIC decreases from 1.12 to 0.44 MPa· m0.5, and fracture energy decreases from 211.50 to 56.00 J·m-2. The fracture performance of Marine shale is significantly controlled by bedding orientation and acid concentration. (3) According to the acid filtration mechanism, the acid intrusion mode of parallel bedding shale is fracture intrusion, and the intrusion depth is large; The vertical bedding shale is a porous intrusion with shallow intrusion depth. The shear stratification sample is of mixed type. (4) Shallow shale fracturing needs to inhibit fracture propagation along the horizontal bedding direction, and acid injection is not appropriate; In deep shale fracturing, acid injection is suggested to enhance the horizontal expansion trend of fractures and improve the complexity of fractures. Due to the limitation of vertical fractures in ultra-deep shale fracturing, it is suggested that acid injection should be cautious to promote vertical fracture expansion and realize the expansion of ultra-deep reservoirs.
  • Ma Lina, Li Shiping, Qiao Danyang
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 1997-2004. https://doi.org/10.20174/j.JUSE.2025.06.16
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Through the indoor unconfined compressive strength test, focusing on the influence of different stone content and fine-grained material proportion on its compressive strength. Combined with the XTDIC analysis system to obtain the strain and displacement cloud diagrams of the sample, and the qualitative and quantitative analyses is carried out. The results show that under static loading, the increase in strain leads to an accelerated rate of initial stress growth and a significant increase in peak stress, and the maximum increase of adjacent peak stress is 22.2%, after reaching the peak, the strain continues to increase while the rate of stress reduction slows down, and the residual stress rises; when the content of stone powder in the specimen is fixed and the stone content varies, the growth rate of initial stress is accelerated and leads to the elevation of peak stress. From the stress-strain curve analysis, it can be seen that the specimen shows the best mechanical properties and the smallest deformation when the stone content is 70% and the stone powder accounts for 80%; the specimen shows the worst mechanical properties when the stone content is 30% and the stone powder accounts for 20%. The effect of fractal dimension on the unconfined compressive strength after compaction was significant, and with the decrease of fractal dimension after compaction (the increase of stone content), the unconfined compressive strength increases first and then decreases.
  • Chen Yongqi, Xiao Yao, Deng Huafeng, Zhu Wenxi, Huang Xiaoyun
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2005-2016. https://doi.org/10.20174/j.JUSE.2025.06.17
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to enhance the microbial improvement effect of red mudstone filler,calcium lignosulfonate and basalt fiber were used as admixtures to cooperate with microorganisms to improve the red mudstone filler. The MICP mixing method was used to test the physical and mechanical properties of the red mudstone filler. The effects of different amounts of calcium lignosulfonate and basalt fiber on the unconfined compressive strength, failure characteristics, no-load expansion rate and disintegration resistance of the red mudstone filler treated by MICP mixing and their changing laws were studied through relevant experiments. The mechanism of calcium lignosulfonate-basalt fiber and microorganisms in the improvement of red mudstone filler was revealed. The results show that: (1) Compared with the MICP group, when the content of calcium lignosulfonate and basalt fiber is 3% and 0.2% respectively, the synergistic improvement effect of the sample is the best, the longitudinal wave velocity is increased by 39.49%, the unconfined compressive strength is increased by 134.55%, and the peak strain at failure is increased by 32.34%. The calcium lignosulfonate-fiber and MICP synergistic effect can significantly improve the compactness, strength, and deformation resistance of red mudstone filler. (2) After the addition of calcium lignosulfonate and basalt fiber, the no-load expansion rate of the red mudstone filler is reduced by 53.84%~73.07%, and the disintegration of the optimal content group sample is improved from moderate disintegration to no disintegration. (3) The mechanism of calcium lignosulfonate and basalt fiber synergistic with microorganisms to improve the red mudstone filler mainly includes three aspects: the cementing and filling effect of lignin polymer, the reinforcement effect of basalt fiber, and the synergistic promotion effect of calcium lignosulfonate and basalt fiber on MICP. The research results can provide new ideas for the improvement of red mudstone filler materials in high-speed railway subgrades.
  • Ding Weiwei, Huang Guangli
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2017-2025. https://doi.org/10.20174/j.JUSE.2025.06.18
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    During the long-term operation of salt cavern gas storage, the duration of pressure changes within the reservoir is relatively long, and the frequency of pressure changes is relatively low. The deformation process of the reservoir's surrounding rock is a creep deformation process under low-frequency cyclic loading. To deeply study the mechanical properties of salt rock under long-term creep conditions, three sets of uniaxial creep tests on salt rock were conducted. The test results show that: (1) For damaged salt rock samples, the deformation process is divided into three stages: decelerated deformation, stable deformation, and accelerated deformation. The corresponding stress-strain curve exhibits a trend of "sparse"-"dense"-"sparse". (2) Among the two influencing factors of loading rate and stress upper limit, the stress upper limit has a more significant impact on the mechanical properties of salt rock. Under the same stress amplitude, as the number of cycles increases, the impact of different loading rates on salt rock deformation gradually decreases. (3) During the loading stage, the deformation rate of the salt rock samples gradually decreases over time until the stress increases to 60% of the stress amplitude, at which point the strain rate starts to increase. (4) In the constant stress upper limit stage, the strain rate gradually decreases over time, but significantly increases when the sample is close to failure. (5) Throughout the test process, the strain rates during different cycles are nearly equal in both the loading stage and the constant stress upper limit stage. In the final cycle before sample failure, the strain rate in the constant stress upper limit stage increases noticeably. (6) In the unloading stage, the deformation rate increases as the load decreases, with the elastic strain of the salt rock recovering in this stage, showing a negative deformation rate. (7) During the constant stress lower limit stage, the strain rate is close to 0 and fluctuates in a negative value state.
  • Fan Xiang, Hong Ming
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2026-2039. https://doi.org/10.20174/j.JUSE.2025.06.19
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The stress of soft and hard interbedded mass is complex. Therefore, pre-fabricated specimens of interbedded surrounding rock tunnel in different thickness ratio of soft and hard interbedded mass were prefabricated to study the deformation and failure characteristics around the interbedded rock tunnel through experiments. The stress and strain distribution around the interbedded rock tunnel were analyzed and compared by numerical simulation. The results show that: When the relatively soft layer is thick, the damage of the tunnel is more affected by the compression of the relatively soft layer, and the axial compression deformation of the tunnel is larger. When the relatively hard layer is thick, the damage of the tunnel is more affected by the fracture of the relatively hard layer. When the interlayer thickness is thin, the thinner rock layer makes the deformation and failure of rock mass around the tunnel more continuous. (1) The overall failure mode of the test specimen under the condition of 10 mm relatively hard layer and 10 mm relatively soft layer is strong hole failure mode, in which the tensile failure is more serious at the vault and arch bottom, and pressure-shear composite crack is easily occurs at the arch shoulder and arch foot. (2) The overall failure mode of the test specimen under the condition of 10 mm relatively hard layer and 20 mm relatively soft layer is strong compression failure mode, in which the compression-shear failure is obvious at the arch waist, and the compression-shear composite crack is easily occurs at the arch shoulder and arch foot. (3) The overall failure mode of the test specimen under the condition of 20 mm relatively hard layer and 10mm relatively soft layer is fracture failure mode, in which the vault obviously sinks, and the stress level near the arch waist is high, but the failure degree of the arch waist is small.
  • Luo Hongchuan, Ma Baosong, Zheng Lining, Bi Jingjie, Zhao Yahong
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2040-2050. https://doi.org/10.20174/j.JUSE.2025.06.20
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The elderly water supply pipelines rehabilitated with cured-in-place-pipe (CIPP) technology can withstand internal water pressure loads. However, the existing pipeline's hole defects provide seepage channels for groundwater around the pipe, leading to a significant increase in local static water external pressure and causing instability and failure of the lining. Through static water external pressure tests and finite element analysis of two full-scale pipeline groups, DN 600 and DN 1000, this study investigates the deformation patterns and instability characteristics of the cured lining under the influence of static water external pressure. It explores the influence of lining DR values and ellipticity on the buckling behavior of the inner lining under interface bonding conditions. The results indicate that: The inner lining undergoes five stages under static water external pressure: uplift, bending, ultimate external pressure buckling, local fiber fracture, and lining fracture. The bonding strength of the interface significantly affects the buckling of the lining. When repairing water supply pipes using CIPP technology, it is essential to ensure the durability of the bonding performance and bonding effect between the CIPP lining and the existing pipeline.
  • Peng Bo, Wang Zhonghao, Mao Long, Wu Xiaoquan, Yang Jiayi
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2051-2059. https://doi.org/10.20174/j.JUSE.2025.06.21
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Rock anchorage is a new type anchorage structure of suspension bridge with complex stress mechanisms. In order to study the deformation characteristics and bearing capacity of the suspension bridge rock anchorage, the 1∶10 scaled model tests on site are carried out based on the Xihoumen Rail-cum-Road Bridge rock anchorage project. By sequentially conducting tests on 1.2P prestressing condition loading, 1P design load condition loading, and overload condition loading, the deformation distribution law of the rock anchorage under various working conditions are analyzed and the bearing characteristics of the rock anchorage are studied. According to the research, the deformation of rock anchorage mainly occurs in the two rock anchorage bodies and the rock mass area between them under various working conditions. And the deformation distribution curve is characterized by a bimodal shape centered on the rock mass between the anchorages.Under the design load condition, the maximum deformation of the actual bridge rock anchorage is 4.17 mm. Moreover, and the safety stability coefficient of the actual bridge rock anchorage is 11. Both deformation and safety stability coefficient of the suspension bridge anchorage are meet the requirements of the specifications. The research results may provide the experimental basis for the safety evaluation of the rock anchorage of the Xihoumen Rail-cum-Road Bridge, and also provide the technical reference for the research and application of rock anchorage related issues.
  • Wang Xuesong, Guo Lianjun, Zhang Jiuyang, Yan Dayang, Xu Zhenyang
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2060-2070. https://doi.org/10.20174/j.JUSE.2025.06.22
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In the process of blasting rock breaking, the existence of joints and fissures in rock mass affects the blasting crushing effect. Some finite element numerical simulation software can realize the blasting rock breaking simulation of jointed rock mass, while it is difficult to realize the quantitative description of cracks. A numerical model of jointed rock mass was established using LS-DYNA software to analyze the obstruction and reflection effects of joint surfaces on explosion stress waves. The images of blast-induced cracks were binarized, and crack extraction from the binary images was achieved using the Hough transform method. Statistical results of the crack images were subsequently analyzed through the maximum likelihood estimation method. The results show that: The existence of joints seriously hinders the propagation of blasting stress wave, and the peak stress is reduced by about 60%. The rock mass between blast hole and joint is subjected to stronger stress superposition due to reflection, and the crushing effect is more remarkable; In the initial stage of explosion, the occurrence probability of small cracks is large, and at the end, it is dominated by long cracks. In this process, many short cracks expand into long cracks, and the distribution of crack angles is relatively average. The distribution of long cracks with joint angles of 45° and 60° is more uniform, and the extreme value is not obvious. When the joints are 15° and 30°, the crack development in the vertical direction is obvious. The Hough transform method realizes the quantitative result statistics of crack images in the numerical simulation results.
  • Liang Shihua, Yu Yang, Feng Deluan
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2071-2093. https://doi.org/10.20174/j.JUSE.2025.06.23
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    A large area of deep coastal soft soil layer is widely distributed in the coastal area. It has the remarkable characteristics of 'three-high, two-low, and one-long characteristic' with high water content, high compressibility, high organic matter content, low bearing capacity, low permeability, and long deformation stability time. Solidification treatment is one of the most used and most effective methods. An in-depth understanding of the solidification mechanism, strength characteristics, and durability of solidified soft soil is crucial for ensuring the long-term safety and stability of solidified coastal soft soil foundation. The strength and durability of the cement solidified coastal soft soil are reviewed. The multi-scale coupling solidification mechanism of cement on coastal soft soil is clarified. The shortcomings of the current research on solidified coastal soft soil are discussed, and the potential research suggestions with cutting-edge are discussed, which provides new thoughts and ideas for the follow-up scientific research and engineering application.
  • Pu Rumin, Zhu Pengcheng, Li Xiaokun, Liang Ninghui, Yang Fang
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2094-2103. https://doi.org/10.20174/j.JUSE.2025.06.24
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Mudstone exhibits strong hydrophilicity, leading to volumetric expansion upon exposure to water, which adversely affects the construction and operation of tunnel engineering projects. To investigate the impact of rock expansion forces on the mechanical characteristics of secondary lining in tunnels, the Heijing Tunnel in Chuxiong Prefecture, Yunnan Province, is utilized as a case study. Based on the theory of stratum structure method, a comprehensive analysis was conducted on the mechanical characteristics and safety of secondary lining under both the influence and absence of rock expansion forces. Additionally, an economically viable design scheme for anti-expansion secondary lining is proposed. The research results indicate that: Variations in the rise-to-span ratio of the invert primarily affect the secondary lining in the region below the tunnel's springline (including the side walls, arch foot, and invert). An increase in the rise-to-span ratio can optimize the stress state of the structure, thereby enhancing its safety. Under the influence of expansion forces, the safety of the secondary lining will diminish, and the most critical structural positions will shift in response to variations in the rise-to-span ratio of the arch. When the expansion force is within 100 kPa, the structural safety and crack requirements can be met by increasing the reinforcement ratio. However, when the expansion force exceeds 100 kPa, adopting a rise-to-span ratio of 0.12 significantly mitigates the adverse effects of the expansion force on the structure, allowing for a reduction in the reinforcement ratio and a corresponding decrease in costs.
  • Huang Dawei, Chen Kai, Chen Houhong, Chen Yongqing, Yu Jue
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2104-2111. https://doi.org/10.20174/j.JUSE.2025.06.25
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The attitude control of the shield machine is mainly realized by the telescopic pushing jack and the control of the elongation of the over-excavation cutter. Due to the complexity of the curve line of the transition curve section, the attitude of the shield machine is often difficult to control. In view of the above problems, based on the over-excavation mode and plane attitude analysis of the shield machine, the reasonable over-excavation mode of the shield machine is determined to be the inner half-circle over-excavation, and the difference between the inner and outer pushing strokes and the maximum elongation of the over-excavation cutter is calculated. In the calculation process, a two-dimensional plane is used to simplify the jacking jack of the shield machine, and only the influence of the jacks at the left and right ends on the attitude of the shield machine is considered. Finally, the calculation formulas of the jacking stroke difference and the maximum elongation of the over-excavation cutter of the shield machine in the straight-line transition curve section and the transition curve section are obtained. Taking the construction data of a project into the above formula, it is concluded that the push stroke difference between the inner and outer sides of the straight line into the transition curve section and the maximum elongation of the over-excavation cutter is in a quadratic parabola relationship with the tunneling distance of the shield machine, and the transition curve section. The research can provide a scientific basis for the plane attitude and over-excavation control of the shield machine when the shield machine is over-gentle curve.
  • Liao Yi, Zhang Ao, Tao Tiejun, Liu Enlong, Xu Taolong
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2112-2121. https://doi.org/10.20174/j.JUSE.2025.06.26
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In response to the blockage situation in the discharge pipeline of the slurry circulation system in the slurry shield machine, a two-way coupling method based on Computational Fluid Dynamics (CFD) and Discrete Element Method (DEM) is used to simulate the transportation process of slurry and stone particles interacting in corrugated pipelines. The evolution of the particulate phase and the fluid phase over time was studied. The results show that: The corrugated pipeline achieves an accelerated effect on the converging pipe flow field with inward corrugations. The flow velocity of the slurry significantly increases when passing through the corrugated structure, and the high-speed flow area expands. More stone particles at the bottom of the pipeline change from sliding friction with the pipe wall to rolling friction, reducing the frictional wear between the stones and the pipeline, and the slurry at the bottom tends to flow upwards, giving the stone particles an upward velocity component, distributing them in a position with stronger flow field effects. Near the corrugated pipe wall, the turbulent kinetic energy increases with the flow of the slurry, while the friction loss between the slurry and the pipe wall and the thermal energy dissipation caused by the turbulent flow significantly reduce the total pressure. Increasing the amplitude of the corrugated function or decreasing the period is beneficial in alleviating the blockage of stone particles, but an excessively small period can cause individual stone particles to be retained for extended periods. Using a corrugated acceleration structure in pipe sections where stone particles tend to accumulate can effectively prevent clogging.
  • Cao Zhensheng, Xu Haibin, Hu Xiaomin, Zhang Guangkui, Liu Keyuan
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2122-2132. https://doi.org/10.20174/j.JUSE.2025.06.27
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Tunnels crossing fault zones are common in railway engineering, water conservancy engineering and highway engineering. Engineering geological disasters induced by crossing fault fracture zones account for more than 50% of the total number of tunnel disasters in China. Ensuring the safety and stability of tunnels crossing fault zones is the top priority in the development of tunnel engineering in China. Based on the engineering background of Tabaiyi Tunnel, this paper explores the effective support measures of surrounding rock under the influence of fault fracture zone. In order to solve the problem of large deformation of soft rock in Tabaiyi Tunnel, this paper first explores the main controlling factors of large deformation by carrying out on-site point load test, in-situ stress test and indoor mineral composition analysis test. Then, according to the deformation mechanism of Tabaiyi Tunnel, a high pre-tightening force and long-short NPR anchor net support scheme is proposed. According to the excavation compensation theory, the high pre-tightening force of NPR anchor cable is used to compensate the stress of tunnel surrounding rock. Through numerical simulation and field monitoring, it is shown that the tunnel can effectively reduce the large deformation of surrounding rock under the influence of fault fracture zone under the high pre-tightening force and long-short NPR anchor net support scheme. The research results can provide reference for tunnel support crossing fault zone.
  • Jiang Chuanbin, Yuan Peng, Shang Xueyi, Chen Jie, Yuan Qiang
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2133-2140. https://doi.org/10.20174/j.JUSE.2025.06.28
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Distributed acoustic sensing (DAS) technology plays an increasingly important role in underground engineering safety monitoring where the identification of intrusion events and source localization are key issues. In response to the issues of poor location stability of single method-based P-wave arrival time data and inconsideration of P-wave arrival time data uncertainty in source localization, this study proposed an intrusion event Bayesian source localization method for distributed acoustic sensing that adapted the selection of P-wave arrival time data and considered its uncertainty. Specifically, we established a source localization objective function for intrusion events that attenuates with distance, developed a P-wave arrival time data adaptive selection and source preliminary localization method, and then proposed a Bayesian source localization method that considered the uncertainty of P-wave arrival time data, providing the source localization results and its uncertainty. Tests on intrusion event source localization in the drainage tunnel show that: (1) Sensor data closer to the intrusion event location have significantly greater localization weight, while the weights of sensor data from further distances tend to zero; (2) For sensors closer to the propagation distance, the adaptively selected P-wave arrival time is consistent with the theoretical propagation time; (3) The Bayesian source localization considering P-wave arrival time uncertainty converges quickly, with the source localization results showing a circular distribution of uncertainty in the horizontal direction and an elliptical distribution in the vertical direction, with an intrusion event localization error within 1 meter; (4) The average location errors of rock fall, water release, human walking, and manual knocking tests are 0.98 m, 0.81 m, 1.82 m and 0.35 m, respectively. In summary, this study provides a new approach to the adaptive selection of P-wave arrival time data and the analysis of uncertainty in source localization results.
  • Hu Wei, Tang Yin, Tan Xinrong, Yü Yü
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2141-2147. https://doi.org/10.20174/j.JUSE.2025.06.29
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    China preliminarily completes the transition of limit state design method for railway tunnel currently and establishes the limit state design method for railway tunnel structure under self weight and surrounding rock load. However, detailed research are not conducted on the limit state method for railway tunnel under seismic load. In order to improve the limit state design system of railway tunnel, clarify the target reliability index of railway tunnel structure for seismic resistance and establish the limit state method for railway tunnel for seismic resistance, the paper calibrated the reliability of railway tunnel lining and open cut tunnel in current general reference maps for seismic resistance by calibration method, and put forward the target reliability index of railway tunnel for seismic resistance based on calibration results. The results show that: The reliability index of railway tunnel structure for seismic resistance is mainly related to the surrounding rock level and structure type, specifically, the reliability index of the lining is significantly lower than that of the open cut tunnel, and the reliability index of the structure under Grade V surrounding rock is significantly lower than that under Grade IV; After considering the weight coefficients of various structures under different surrounding rock levels, the average reliability index of railway tunnel for seismic resistance is 6.81, which is greater than target reliability index of 4.2 for railway tunnel structure under self weight and surrounding rock pressure, the reliability level of the current structure for seismic resistance is relatively high; Based on the calibration results of the current general reference maps of railway tunnel for seismic resistance and the target reliability index in current code, the target reliability index of railway tunnel for seismic resistance is determined to be 4.7 for Level 1, 4.2 for Level 2, and 3.7 for Level 3.
  • Xie Cheng, Zhao Zizhen, Liu Jun, Xiao Shiguo, Dai Tianyi
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2148-2158. https://doi.org/10.20174/j.JUSE.2025.06.30
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Petroleum pipelines across slopes can produce obvious mechanical responses to seismic landslide thrust. In order to rationally calculate internal forces and displacements of the pipeline due to the seismic landslide of the soil mass surrounding it, a mechanical analysis method of the pipeline is proposed based on the Timoshenko beam considering lateral shear deformation on the Winkler foundation with the normal trapezial earth pressure on it, in which the seismic earth pressure is determined using the transfer coefficient method combined with the pseudo-static method. The formulas of internal forces and displacement of the pipeline are accordingly derived, and the material determination approach is put forward according to the control of its safe stress state. Analysis results of an example show that the pipeline deflection increases with the distance to the slide mass boundary, and the maximum bending moment and shear force on the pipeline are nearly located in the boundary, which is consistent with that by the numerical simulation via FLAC3D. The distribution mode of the internal forces is not almost varied with seismic coefficients, pipeline length in the slide mass, and cross-section of the pipeline. However, the maximum bending moment and shear force are nonlinearly increasing with the horizontal seismic coefficient, slightly increasing with the cross-section, and marginally decreasing with the increase of the pipeline length in the slide mass.
  • Ma Xingdong, Zhang Shishu, Li Changyou, Yuan Guoqing, Jia Xiang
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2159-2166. https://doi.org/10.20174/j.JUSE.2025.06.31
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The deep buried large underground powerhouse cavern group of Shuangjiangkou Hydropower Station is arranged in the left bank mountain, and the geological lithology is brittle hard rock granite. The measured ground stress reaches 37.82 MPa, which is a high to extremely high stress zone. Research has found that the relaxation depth of underground power plants ranges from 2.2 m to 4.4 m, and the overall damage to surrounding rock shows an increasing trend with layered blasting excavation. However, the relaxation depth variation is small, and the average fragmentation variation of the relaxation layer is less than 10%. The lag effect of surrounding rock fracture is shortened. The GEDZ distribution in high stress brittle hard rock caverns can be evaluated using acoustic detection as the main method, supplemented by innovative semi quantitative indicators such as "rupture morphology" (new cracks, residual blasting hole half hole rate, rock burst failure depth), and quantitative indicators such as permeability to reasonably evaluate the distribution of excavation damage areas in surrounding rock. This can better guide on-site design and construction, and provide reference for similar hydropower underground power plants and Sichuan Tibet Railway projects in the future.
  • Yang Yunhao, Zhang Shunli, Yuan Guoqing
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2167-2179. https://doi.org/10.20174/j.JUSE.2025.06.32
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The temporal and spatial characteristics and mechanisms of perimeter rock deformation in the three large cavern chambers of Shuangjiangkou Hydropower Station are analyzed, and with the plant cavern as the focus of the study, the monitoring feedback design for monitoring during the construction period was carried out using perimeter rock deformation, support load monitoring data as well as drilling acoustic wave detection data, and by adopting the inversion method of the mechanical parameter that took into account the dynamic evolution of the relaxation zone around the cavern. The results show that: The spatial characteristics of rockmass deformation result from superposition of unloading effect, the cavern group effect and high wall effect with unloading effect as a main factor which result in unsymmetrical deformation of powerhouse; Thanks to the good quality of granite rockmass, the high ratio of rock strength to in-situ stress and scarcity of weak structural plane, the largest peripheral displacement of powerhouse is 28mm; For surrounding rockmass that has been just exposed or merely experienced 2~3 underlying slicing excavation, the only component of total displacement is excavation induced rebounding deformation; For rockmass that has experienced more than three underlying slicing excavation, time related deformation due to progression degradation of rockmass can be found in total displacement but occupy a small proportion and converge very quickly after the excavation of powerhouse has been finished; In case of time effect existed properly modeling of damaged zone around cavern is necessary and numerical calculation based on such model can give a good prediction of deformation increment. There exists no abnormal deformation in surrounding rockmass when the intensity of rock bolt support is decreased at some important construction stages according to prediction of future deformation, which means that feedback design of support is reasonable.
  • Xiong Zhongming, Zhang Jiabin, Chen Xuan, Liang Jiawei
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2180-2191. https://doi.org/10.20174/j.JUSE.2025.06.33
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In cities with dense ground fissure distributions, linear infrastructures like underground utility tunnel (UUT) inevitably intersect with ground fissures. The vertical displacement between the hanging wall and footwall poses a significant threat to the safety of UUT in ground fissure site. For this reason, taking the UUT in Xi'an Happiness Forest Belt as the prototype structure, 3D finite element models that considered the soil-structure interaction were developed. The strain distribution, damage characteristics of the UUT, and the stress variations in the site soil under different dislocation magnitudes were analyzed. In addition, a damage index was introduced to establish the classification criteria (DS0-DS3) for the safety evaluation of the UUT in order to assess the structural damage. The results indicate that: The tensile damage areas in the UUT are concentrated in the upper level compartment of footwall, the wall at the fissure location, and the lower level compartment of hanging wall; the compressive damage areas are concentrated in the wall at the fissure location; the main damage in the UUT is primarily controlled by tensile damage; the weak areas of the UUT are the bottom slab of the lower level large compartment in hanging wall, the side wall of the lower level large compartment at the fissure location, and the top slab of the upper level large compartment in footwall. When the vertical displacement reaches 300 mm, the plastic zone of the upper soil of the hanging wall UUT develops upwards and runs through the soil. Most of the damaged areas of the UUT reach DS3, causing large-scale concrete peeling and structure failure.
  • Liu Yuyang, Zhang Yixiao, Ma Desheng, Zhang Haonan, Chen Mingjie
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2192-2202. https://doi.org/10.20174/j.JUSE.2025.06.34
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Crystal disease is very common in mountain tunnel drainage system, and the establishment of quantitative evaluation method of crystal risk has great practical engineering significance. Based on the analysis of the chemical crystallization process and physical plugging process in the tunnel drainage tube, seven solution parameters were identified as the key judgment indicators for identifying crystallization risk. The influence tests of single factor crystallization and orthogonal crystallization were designed and carried out, and the influence rules and degree ranking of each judgment index on crystallization were determined respectively. The index grading criteria were determined according to the change rate of calcium ion concentration. The index weights were calculated by the improved analytic hierarchy process (AHP), the crystallization risk evaluation system was established, and the comparison with the field situation was verified. The results show that the crystallization plugging process of mountain tunnel includes three stages: calcium osmosis in the first branch, drainage crystallization reaction and physical precipitation in the tube. In the drainage solution parameters, Ca2+, HCO-3, pH, SO2-4, flow rate, Mg2+ and temperature can be used as indicators to judge the crystallization risk. The first five have a positive promoting effect on crystallization, and the last two have a negative inhibiting effect. The relative influence degree of crystallization reaction is pH > Ca2+ > flow rate > HCO-3 > SO2-4 > temperature > Mg2+ from large to small, and the corresponding weight varies between 0.319 and 0.034. Based on the quantitative curve of calcium ion concentration change rate, the scoring criteria can be divided into 4~6 intervals, and the scoring criteria can be determined by curve slope interpolation. The engineering application shows that the crystallization risk grade evaluation system of tunnel drainage system can objectively reflect the crystallization risk grade of target tunnel.
  • Su Peidong, Shi Shutong, Li Yougui
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2203-2216. https://doi.org/10.20174/j.JUSE.2025.06.35
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In response to the unsatisfactory treatment effect and high construction risk of shallow buried tunnels crossing high gas formations, based on the theory of fluid porous media, the coal seam gas displacement and extraction technology is applied to the construction treatment of shallow buried tunnels in high gas formations. The performance of gas displacement in the formation was analyzed using theoretical analysis and numerical simulation methods. The influence range and time of gas displacement were determined. Four factors that may affect displacement performance, including negative pressure of drainage, gas injection pressure, formation permeability, and borehole diameter, were studied. Finally, the numerical simulation results of gas displacement in the formation were demonstrated through on-site experiments. The results show that under the condition of negative pressure PC=50 kPa and injection pressure PZ=2.0 MPa, when the distance between injection and extraction holes is 12 m, gas injection for 40 min is used, and the displacement effect achieved by pumping for 90 min is better. The smaller the suction negative pressure, the better the displacement performance, but the impact of suction negative pressure on the displacement process is relatively small; The increase of gas injection pressure can significantly increase gas drainage capacity and reduce displacement time; There is a nearly linear positive correlation between the decrease in formation gas, gas flow rate from drainage holes, and permeability; For the commonly used borehole diameter, increasing the borehole diameter has little effect on the effectiveness of gas displacement control in the formation. Due to factors such as geological conditions, equipment limitations, and construction disturbances, there are certain differences between on-site experiments and numerical simulations of gas content monitoring data. However, the two have a high degree of consistency in terms of displacement effect and time progress.
  • Zhai Yue, Han Zhulin, Qu Lu, Han Ce
    Chinese Journal of Underground Space and Engineering. 2025, 21(6): 2217-2226. https://doi.org/10.20174/j.JUSE.2025.06.36
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to enhance the safety and efficiency of pedestrian-vehicle evacuation in underground parking lot fires, a background field model based on cellular automata is constructed using Python. The model incorporates fire field data of the PyroSim simulation to create a dynamic coupling model of fire and pedestrian-vehicle evacuation. A simulation study is conducted in the underground parking lot of a commercial complex in Xi'an. The results show that, compared to pedestrian evacuation, the Required Safety Egress Time (RSET) for pedestrian-vehicle evacuation is extended by 20% to 130%. Safe evacuation is only achievable under specific conditions of pedestrian density and parking space occupancy. The arch effect is observed at the exit when pedestrian density is high, which reduces the evacuation efficiency of both pedestrians and vehicles. The probability of vehicle evacuation failure is approximately 22 times higher than that of pedestrians. Even if pedestrians successfully reach the exit, 15% to 25% still suffer from varying degrees of injury, with the proportion of serious injuries and near-death cases positively correlated with the pedestrian evacuation density. These findings can serve as a reference for optimizing fire protection design in underground parking lots and formulating evacuation plans for both people and vehicles.