20 June 2026, Volume 22 Issue 3
    

  • Select all
    |
  • Yin Ziyuan, Wang Xinran, Cui Pengjie
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 753-760. https://doi.org/10.20174/j.JUSE.2026.03.01
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The renewal of the area around high-density urban rail transit stations is an important topic for the revitalization of existing urban public spaces. The pedestrian network and the infrastructure around it have a profound impact on the aggregation of passenger flow at rail transit stations. As urban rail transit lines continue to expand across China, a diminishing marginal effect is observed in the incremental passenger flow at various stations. To ensure a more effective alignment between the development intensity, land use nature, public facility density, and passenger flow intensity in urban renewal projects surrounding rail transit stations, and to fully leverage the agglomeration effects of spatial elements within these areas, this paper takes the urban renewal of the surrounding region of the Hi-Tech Central Station in Shenzhen as a case study. Firstly, based on the empirical data of well-functioning passenger flow and station area development intensity, the Gaoxinzhong Station is identified as an arrival-type station among the categories of arrival, departure, and balanced stations. Secondly, the quantitative relationship between the planned passenger flow and anticipated development intensity at the Hi-Tech Central Station is analyzed, revealing that the urban renewal capacity of the station area fails to meet expectations. Finally, a quantitative assessment is conducted on the ground and underground building space capacity, road network structure, and other relevant factors of the surrounding plots around the Hi-Tech Central Station, and a comprehensive evaluation is made of the potential induced passenger flow aggregation capacity of each plot, providing a design basis for achieving the planned capacity.
  • Xue Minghui, Wang Qingyan
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 761-770. https://doi.org/10.20174/j.JUSE.2026.03.02
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    With the transition towards the trend of multifunctional urban development, various types of public spaces have become increasingly integrated and polysemous. As a unique form of urban public space and an integral part of underground commercial streets, entrance spaces not only serve as vital nodes connecting above ground and underground areas but also as focal points accommodating diverse urban activities. However, many underground commercial street entrances in Northeast China face challenges such as aging facilities, outdated functions, and low spatial utilization rates, resulting in insufficient vitality and an inability to meet the diverse needs of modern urban life. Implementing scientific and reasonable strategies for renewal and revitalization to fully unlock the potential of these entrance spaces has become essential for improving the quality of underground commercial streets and enhancing the overall vibrancy of urban public spaces. A typological approach was adopted to categorize existing underground commercial street entrances in Northeast China. From the perspective of polysemy, it establishes an evaluation system based on three dimensions: environmental relationality, perceptual process, and functional ambiguity. By distributing questionnaires to experts and residents, it identifies and compares the influencing factors of entrance space polysemy and assigns values to polysemy indicators. Based on the evaluation results, the entrances are classified into categories of polysemy, and typical underground commercial street entrance spaces catering to diverse activity demands are summarized. Finally, optimization strategies for the polysemy of these typical entrance spaces are proposed to promote the renewal of existing underground commercial street entrances, meet the diverse needs of citizens, and enhance the vitality of urban underground spaces.
  • Sun Tao, Yang Yukui, Liu Mingyu, Wu Jiao, Tan Jingwen
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 771-779. https://doi.org/10.20174/j.JUSE.2026.03.03
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In recent years, the scale and speed of underground space construction have grown rapidly, and the development of deep underground spaces has become increasingly frequent. Under the background of the national spatial planning system, the urban underground space planning of Guangzhou city adheres to the people-centered approach, bottom line thinking, coordinated development, and safety, and explores the system and methods for preparing the underground space planning of super large cities. Firstly, the safety bottom line, resource protection bottom line, and control bottom line range are proposed as areas where underground space is not suitable for construction; Priority should be given to safety and comprehensive resource assessment should be carried out to determine suitable construction areas for underground spaces and general construction areas; To coordinate the comprehensive utilization of underground space, functional composite three-dimensional layering, graded zoning and classification guidance; To achieve deep integration of above ground and underground facilities, infrastructure+ is proposed based on the characteristics of underground facilities; To build a safe and resilient underground space, identify disaster risks, classify and stage safety prevention measures, and coordinate above ground and underground disaster prevention and reduction facilities; To improve the efficiency of planning and control, and assist in the sustainable development of underground spaces, it is proposed to strengthen the guarantee mechanism.
  • Wang Zhichao, Lin Yanghao, Tian Yinghui, Zhang Chunhui, Luo Guangcai
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 780-787. https://doi.org/10.20174/j.JUSE.2026.03.04
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To overcome the limitations of the conventional super-subloading surface constitutive model in capturing soil stress-induced anisotropy and the nonlinear degradation of shear modulus at small strains, this study presents an enhanced model. By incorporating the g(θ) method and integrating classical small-strain stiffness theory, the proposed model offers improved representation of both anisotropic behavior and the nonlinear shear modulus reduction under small-strain conditions. The new model was subsequently applied to predict triaxial shear test results for Shanghai soft clay, Fukakusa clay, and Hefei slightly expansive clay, as well as to simulate the deep excavation of the Huifu Road Station in Hefei Metro. The results demonstrate that the proposed model effectively captures the high initial shear modulus and its nonlinear attenuation under small-strain conditions, unifies the application of four yield criteria (von Mises, Mohr-Coulomb, Matsuoka-Nakai, and Lade-Duncan) to characterize stress-induced anisotropy, and successfully describes both the structural shear-shrinkage softening of soft clay and the shear-dilation softening of overconsolidated soil. The improved model not only effectively characterizes complex mechanical behaviors of natural soil but also accurately predicts deformation patterns of retaining piles during excavation processes.
  • Mei Yuan, Liu Ziyang, Zhou Dongbo, Wang Yifei, Zhang Yuhang
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 788-799. https://doi.org/10.20174/j.JUSE.2026.03.05
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    With the acceleration of urbanization, underground space development has become an important way to relieve the pressure of urban traffic. In composite formation, the construction of a double-line large-diameter mud-water shield tunnel faces complicated geological conditions and construction problems. Based on the measured surface settlement data of a double-line tunnel project in Hangzhou, this study adopts a mathematical method to perform regression analysis and introduces the correction coefficient α of the maximum surface settlement and the correction coefficient β of the width of the settlement trough to correct the double-line Peck formula. Meanwhile, combining the superposition principle, an applicable correction method for the superposed section of the double-line tunnel in composite formation is proposed. In order to accurately predict the ground settlement law of double-line large-diameter mud-water shield tunnel construction in composite formation. The results show that there is a great difference between the measured value and the predicted value of the Peck formula of the double-line large-diameter mud shield before correction. Examples show that the superposition section correction method can effectively improve the fitting degree of the superposition section prediction curve when the center distance L of the two-track tunnel is greater than (i left +i right). By comparing the corrected coefficient range ofthe Peck formula with different diameters, it can be seen that the large-diameter shield is better than the small-diameter shield in vertical settlement control, and the distribution range of the corrected coefficient of sedimentation tank width is larger than that of the small-diameter shield. The research can provide a reference for the design and construction of large-diameter shield tunnels in China.
  • Pan Jingtao, Guo Zhuang'ao, Zhao Dan, Fan Yuzhong, Qin Gengmu
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 800-812. https://doi.org/10.20174/j.JUSE.2026.03.06
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Understanding the formation and evolution of plastic zones is fundamental to controlling the deformation and instability of surrounding rock in deep roadways. However, existing plastic zone theories primarily focus on analytical solutions for circular roadway radii under uniform stress fields. To address this limitation, based on the Drucker-Prager strength criterion and butterfly-shaped failure theory, a single-factor analysis method is adopted to investigate the variation trends and sensitivity of the plastic zone in the surrounding rock of the 1201 return airway in a deep mine. Furthermore, numerical simulations based on FLAC3D are conducted to determine the plastic zone expansion morphology, range, and failure characteristics of different roadway cross-sections. The results show that: Under a non-uniform stress field with a lateral pressure coefficient λ=2.5, the maximum plastic zone depths at the roof, shoulder corners, and floor of a rectangular roadway are 5.3 m, 6.4 m, and 5.0 m, respectively, which are on average 27.5% larger than those of a circular roadway (4.2 m, 5.1 m, 4.0 m). The plastic zone depth in a straight-wall three-centered arch-shaped roadway (4.8 m, 5.7 m, 4.5 m) is 15.2% greater than that of the circular section, demonstrating the significant influence of cross-sectional shape on plastic zone distribution. Field measurements are also employed to validate the numerical results. The study reveals that when the intermediate principal stress coefficient m=0.7, the plastic zone radius reaches its minimum value. However, for m>0.7, stress redistribution leads to a radius rebound, exhibiting an interval effect. Additionally, increasing the surrounding rock strength parameters can reduce the plastic zone radius by approximately 35% to 48%. Field observations confirm that the measured plastic zone depths at the shoulder corners (5.5 m~5.6 m) have an error margin of less than 10% compared to numerical simulations, and the fracture distribution is highly consistent with the butterfly-shaped failure theory. These findings provide valuable insights for analyzing plastic zones in deep roadways and optimizing excavation cross-section designs.
  • Ran Jiaju, Xie Zhitian, Niu Xiaokai, Guo Jinpeng, Song Lihong
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 813-821. https://doi.org/10.20174/j.JUSE.2026.03.07
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To investigate the ultimate bearing capacity of shield tunnel segments during service, a full-scale experimental investigation was conducted. Using three rings of service-aged segments as test subjects, the experimental testing content and loading scheme were formulated. The segments underwent three-stage loading to simulate static loading tests under overburden pressure, with monitoring of overall convergence deformation, joint misalignment, longitudinal joint opening, and segment cracks during the loading process. The stress distribution, failure modes, and ultimate bearing capacity of the middle ring segments were analyzed. The results indicate that under a simulated burial depth of 15 m, the overall structural deformation of the segments was minimal. The deformation of the segments was divided into three stages, with rapid increases in overall convergence deformation when the relative displacement between the crown and invert exceeded 2.83‰D. The failure cracks in the segments were primarily concentrated around the bolt holes and longitudinal joint edges.
  • Chang Lei, Rong Chuanxin, Cheng Hua, Cai Haibing, Cao Yi
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 822-836. https://doi.org/10.20174/j.JUSE.2026.03.08
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In the process of freezing method construction of underground engineering in Huainan and Huaibei mining areas, the deep frozen soil is subjected to loading and unloading in different directions at the same time, showing a dynamic change from the initial original stress state to the increasing stress on one side and the decreasing stress on the other side. The uneven distribution of stress leads to different stress paths of the soil in front of the working face, so that the soil is in a three-way unequal stress state, which in turn affects its mechanical properties. Based on this, in order to further study the mechanical behavior of frozen soil under true triaxial stress state, this study analyzed the effects of different temperatures, initial stress states, and loading and unloading rates on the strength and deformation characteristics of frozen calcareous clay using a self-developed true triaxial apparatus for frozen soil. The test results show that under the same test conditions, when the loading rate increases from 0.4 mm/min to 0.8 mm/min, the slope of the generalized shear stress q and the generalized shear strain curve γ of the frozen calcareous clay gradually increases, that is, the slope of the stress-strain curve of the frozen calcareous clay increases with the increase of the loading rate, and the peak stress corresponding to the strain softening group also gradually increases. In addition, with the increase of initial stress and loading rate, the stress-strain curve of frozen calcareous clay changes from strain softening to strain hardening as a whole. Under different stress paths, the strengthening effect of loading rate on frozen calcareous clay gradually decreases with the increase of loading rate, and satisfies the power function relationship. Under the same freezing temperature T and loading rate v, the triaxial strength of frozen calcareous clay conforms to the Drucker-Prager criterion in the principal stress space, and the relationship expressions of α and K with freezing temperature T and loading rate v are given respectively. The research results can provide a theoretical basis for deep freezing shaft sinking.
  • Zeng Guojian, Li Kegang, Li Mingliang, Qiao Zhiqiang, Gan Yanken
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 837-849. https://doi.org/10.20174/j.JUSE.2026.03.09
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to explore the influence mechanism of joint position on rockburst of straight wall arched roadway under dynamic and static loads, the limestone of Huize Lead-Zinc Mine in Yunnan Province was selected and processed into cube specimens with straight wall arched holes. The joints were prefabricated, and the spacing between the joints and the straight wall was 5 mm, 10 mm, 15 mm, and 20 mm, and the jointless specimens were set as the control. Based on the true triaxial rockburst test system, the rockburst simulation tests of roadways with different joint positions were carried out. Combined with video, acoustic emission monitoring, and fractal theory, the influence of joint position on rockbursts was analyzed from the aspects of rockburst failure process, rockburst pit morphology, acoustic emission signal, and debris characteristics. The results show that: (1) The rockburst failure stage of all specimens were same, and the depth of the rockburst pit decreases with the increase of the distance between the joint and the straight wall. (2) The absolute energy of acoustic emission and the fractal dimension of debris reflect the rockburst intensity. Both of them and the rockburst intensity decrease with the increase of the distance between the joint and the tunnel wall. (3) The farther the joint is from the tunnel wall, the lower the energy consumption of rockburst, and the more obvious the flake characteristics of debris.
  • Guo Yanwei, Yang Tao, Yang Minghui, Song Muyuan
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 850-859. https://doi.org/10.20174/j.JUSE.2026.03.10
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The influence of tunnel excavation on the overlying pipeline cannot be ignored when the distance between the pipeline and tunnel is minor. For the special condition of tunnel excavation through the existing pipeline, the transparent soil and PIV technologies were adopted to conduct a comparative study on the deformation regularities of continuous and non-continuous pipelines. The dynamic change of soil displacement with tunnel excavation was monitored, and the influence of bending stiffness, depth and section length of pipelines on their deformation was analyzed by the control variable method. The results show that: The settlement patterns of both continuous and non-continuous pipelines are basically consistent with the Gaussian settlement curve. The maximum settlement value of continuous pipeline decreases with the increase of flexural stiffness and increases with the increase of pipeline burial depth. The maximum settlement value of discontinuous pipeline increases first and then decreases with the increase of pipeline section length. In addition, the relative position between the non-continuous pipeline and tunnel has a significant impact on the pipeline settlement. When the non-continuous pipeline interface is directly above the tunnel axis, its maximum settlement value is significantly higher than those value when the tunnel axis is directly below the pipeline joint center, indicating that the tunnel crossing position close to the non-continuous pipeline interface is the most dangerous situation which should be avoided as soon as possible in design.
  • Cheng Gong, Zhou Changbing, Wang Yunmin, Li Xiaoshuang
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 860-869. https://doi.org/10.20174/j.JUSE.2026.03.11
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Regarding the influence of cyclic loading and dry-wet cycling coupling effects on the evolution of mechanical properties of dolomite, triaxial compression tests are conducted under different dry-wet cycle counts and cyclic loading upper limits. Combined with macro- and microstructural analysis, the strength degradation, deformation characteristics, and damage evolution mechanisms of dolomite are systematically investigated. The results indicate that: As the number of dry-wet cycles increases, the mass of dolomite decreases, while the water content shows a logarithmic increasing trend. With the increase in dry-wet cycles, the peak compressive strength and elastic modulus of dolomite gradually decrease, and the damage factor exhibits a logarithmic increase. Different upper-limit cyclic loadings have a certain promoting effect on the dry-wet degradation. The elastic modulus of the specimens shows varying trends with the increase in the upper limit of cyclic loading. Dry-wet cycling has no significant impact on the failure mode of dolomite, which remains shear failure or mixed shear-tensile failure. Dry-wet cycling leads to the dissolution of internal soluble salts and colloids in dolomite, as well as the initiation and propagation of cracks, while cyclic loading and unloading at different upper limits accelerate crack propagation and the destruction of cementing materials between particles. When these two factors act together, the deterioration rate of the rock's physical properties is faster than under either factor alone.
  • Duan Xiaoxiao, Yang Dengke
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 870-882. https://doi.org/10.20174/j.JUSE.2026.03.12
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To explore the effect of water content on the mechanical properties and fracture characteristics of fractured rock in a freeze-thaw environment, five prefabricated oval double-hole red sandstones with different saturations were prepared, and CT scanning and uniaxial compression tests were carried out after freeze-thaw cycles. By analyzing the mesoscopic evolution law of rock during freeze-thaw cycle and the macroscopic failure characteristics after freeze-thaw cycles and loading, the relationship between macroscopic failure characteristics and meso-evolution was explored. The results show that: (1) After freeze-thaw cycles, microcracks begin to develop at the tips of existing fracture in rocks with saturation ranging from 60% to 100%, and the total volume of microcracks increases with the increase of saturation. In particular, connected microcracks appear in the rock bridge area of completely saturated rocks after 120 freeze-thaw cycles. (2) With the increase in the number of freeze-thaw cycles, the peak strength of rocks with saturation of 60% to 100% decreases linearly. The elastic modulus of rocks with saturation of 60% to 90% decreases linearly, while the fitting curve of elastic modulus of completely saturated rock is “upper convex”. (3) The failure mode of rocks with different saturations is tensile-shear failure in the early stage of freeze-thaw cycle. With the increase in the number of freeze-thaw cycles, the failure mode of completely saturated rock changes from tensile-shear to tensile. (4) After freezing-thawing, the degree of particle cementation is relatively weak in the position where the microcracks develop, and the macroscopic fractures are easy to occur in this position. The evolution process of microcracks affects the macroscopic failure mode. There is a positive linear correlation between the percentage reduction in peak strength and elastic modulus and the percentage increase in fractal dimension of microcracks. This study can provide a theoretical basis for the prediction of freeze-thaw resistance of rock mass engineering in cold regions.
  • Wu Mengjun, Zhong Zuliang, Xu Miao, Hu Xuebing, Cao Peng
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 883-890. https://doi.org/10.20174/j.JUSE.2026.03.13
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The internal energy evolution process and mechanism of high-performance fiber-reinforced shotcrete single-layer lining during layered construction remain unclear. To investigate the influence of layered construction on the energy evolution characteristics of the single-layer lining structure, this study conducted uniaxial compression strength tests on specimens with six different thickness ratios (the thickness ratio of the second layer to the first layer of fiber-reinforced shotcrete). The results indicate that: The uniaxial compression strength of fiber-reinforced shotcrete is significantly affected by the thickness ratio. The peak strength reaches its maximum value of 43.32 MPa when the thickness ratio is 1, with an optimal stress distribution. When the thickness ratio exceeds 4, debonding and crushing phenomena occur due to the excessive thinness of the lower layer. During the uniaxial compression process, the total input energy of the specimens increases with the thickness ratio. The elastic energy storage capacity is optimal at a thickness ratio of 1, while the dissipated energy exhibits significant fluctuations at a thickness ratio of 5, reflecting concentrated energy release due to local instability. Based on the energy dissipation coefficient analysis, the energy dissipation rate decreases during the compaction stage, indicating a significant reduction in structural stability.
  • Wei Jianxin, Zheng Hong, Liang Ninghui, Huang Yizhan, Wu Mengjun
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 891-899. https://doi.org/10.20174/j.JUSE.2026.03.14
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to study the mechanical properties and reinforcement mechanism of special-shaped cross-section polypropylene fiber reinforced shotcrete, compressive, splitting tensile, and flexural tests were conducted on benchmark shotcrete and special-shaped cross-section polypropylene fiber reinforced shotcrete, and SEM micro scanning tests were carried out. The results showed that the addition of special-shaped cross-section polypropylene fibers can improve the compressive strength, splitting tensile strength, and flexural strength of sprayed concrete. The most significant improvement was observed when the fiber content was 6 kg/m3, the compressive strength at 1 day, 3 days, and 28 days, splitting tensile strength at 28 days, and flexural strength at 28 days were increased by 13.9%, 5.8%, 5.1%, 26.8%, and 8.1%, respectively, compared to the benchmark concrete group. The toughness of sprayed concrete is enhanced by adding special-shaped cross-section polypropylene fibers, and the failure mode changes from brittle to ductile. Special-shaped cross-section polypropylene fibers can reduce the proportion of large pores in the concrete matrix, improve the compactness of the matrix, and reduce the width and number of microcracks at the interface between coarse aggregates and matrix, provide a bridging effect to absorb destructive energy, and enhance the strength and toughness of sprayed concrete. The research results can provide a theoretical basis for the application of special-shaped cross-section polypropylene fiber reinforced shotcrete.
  • Zhou Ting, Xie Jianbin, Xu Yao, Yang Le, Jia Ronggu
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 900-909. https://doi.org/10.20174/j.JUSE.2026.03.15
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Red clay, characterized by high moisture content and a high plasticity index, is widely distributed across the Yunnan Plateau and commonly serves as a predominant foundation soil for infrastructure development. The Punching-compaction dynamic compaction (PCDC) technique, notable for its deep penetration capacity and for forming granular column-composite foundations through gravel-clay mixtures in compaction holes, is widely adopted to reinforce underconsolidated red clay sites. However, its dynamic behavior remains insufficiently studied. Laboratory model tests were conducted to evaluate the effects of filler ratios, compaction energy (8 000~10 000 kN·m), hammer angles, and blow counts on cumulative penetration depth, soil acceleration, and peak pressure. Results revealed that cumulative penetration depth correlates logarithmically with blow counts, with 10 000 kN·m achieving 10% greater penetration than lower energy levels. Higher energy induced stronger vertical/horizontal accelerations, while flat-bottom hammers combined with high energy optimized efficiency. Gravel-soil mixtures enhance column density during densification, with increased gravel content improving soil compactness.
  • Jia Rui, Hao Daiheng, Qiao Xiaoli
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 910-919. https://doi.org/10.20174/j.JUSE.2026.03.16
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Sand production in the reservoir is a key issue restricting the safe exploitation of natural gas hydrate, and reservoir permeability is a key parameter affecting the production capacity of hydrate exploitation. The effects of particle size distribution of reservoir, hydrate saturation, depressurization amplitude and pore size of sand prevention layer on reservoir sand production and permeability during hydrate depressurization exploitation were investigated by laboratory model tests. The correlation between the reservoir sand production and change in reservoir permeability was analyzed. The results show that: (1) Factors influencing reservoir sand production: the lower the clay content, the larger the depressurization amplitude, and the greater the pore size of sand prevention layer were, the larger the amount of reservoir sand production. With the same pore size of sand prevention layer, the amount of reservoir sand production has approximately a linear relationship with the average flow rate of water. (2) Factors influencing reservoir permeability: the reservoir permeability generally increased with time, and finally gradually tended to a stable value. The greater the hydrate saturation, and the greater the clay content were, the smaller the initial value of permeability. The greater the depressurization amplitude, and the larger the pore size of sand prevention layer were, the greater the ratio of final value to initial value of permeability. (3) The correlation between the sand production and variation in permeability: with the same hydrate saturation, the more the amount of reservoir sand production was, the relatively greater the increase amount of reservoir permeability. With the same hydrate saturation and particle size distribution of reservoir, the more the amount of reservoir sand production was, the relatively greater the increase rate of reservoir permeability.
  • Huang Haoyong, Zhong Liqing, Zeng Bo, Xu Ersi, Zhang Fengshou
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 920-927. https://doi.org/10.20174/j.JUSE.2026.03.17
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    With the exploration and development of shale gas resources in southern Sichuan, casing deformation related to fracture and fault shear slip occurs frequently, severely restricting the construction and production of deep shale gas. Therefore, it is significant to determine the friction coefficient of shale gouge under the high-temperature and high-pressure conditions of deep geologic environments from the perspective of fracture/fault friction. In-situ shale samples from three wells (A, B, and C) at different depths of the Longmaxi Formation were used to prepare fault gouges for friction shear tests. The clay mineral content of the samples ranged from 10% to 45%. The results of high-temperature, high-pressure triaxial shear friction tests show that the friction coefficient of the deep shale is between 0.52 and 0.70, and it decreases linearly as the clay mineral content increases. Based on the mineral composition of the Longmaxi shale, the distribution of the friction coefficient was predicted, revealing that the friction coefficient increases with relative depth within a single well. Finally, the slip amount of fractures/faults in the Longmaxi shale was calculated from the stress drop observed in the shear tests. This calculated slip corresponds well with actual casing deformation measurements. These results further illustrate the relationship between fracture/fault slip and casing deformation and provide a corresponding calculation method.
  • Cai Zhi, Xia Peikai, Zhang Tingzhong, Zhang Rui, Shen Zhifu
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 928-936. https://doi.org/10.20174/j.JUSE.2026.03.18
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In geotechnical site investigation, the distance between boreholes often makes the inter-borehole soil layer inferences rely on human experience under conditions of limited in situ tests and sampling. How to conduct intelligent and reliable soil stratigraphic division is a crucial research direction in the current information-oriented development of geotechnical engineering. This paper presents an intelligent soil stratigraphic layer division method by combining Bayesian Compressed Sensing (BCS), Support Vector Machine (SVM) classification, Gaussian Mixture Model (GMM), and Hidden Markov Random Field (HMRF) model. The application flowchart and soil layer division results are presented by taking the Nanjing Yangtze River floodplain ground as an example. The study shows that: BCS can reliably extend the blow count data from standard penetration test (SPT) for subsequent soil layer division; SVM classification can intelligently learn soil boundaries in the two-dimensional space of SPT blow count versus test depth, achieving an initial soil stratigraphic division; based on this, the preliminary optimization of soil layer division can be realized by using the GMM by considering the probability distribution of soil characteristic parameters; finally, the secondary optimization of soil layer division can be realized by using the HMRF model by incorporating spatial correlation constraints (i.e., adjacent points tending to be the same soil type). Combining the four methods can intelligently and automatically divide soil layers, and can significantly improve the accuracy of overall soil layer division and soil layer boundary identification.
  • Liu Shi, Yu Huagui, Wang Jianning, Guo Hong, Yan Shilong
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 937-947. https://doi.org/10.20174/j.JUSE.2026.03.19
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Injecting CO2 into underground formations for permanent storage is one of the key technologies for achieving carbon neutrality goals. Wellblock 38 in Wuqi Oilfield is the core area of the first domestic full-process CCUS demonstration project-the Yanchang Petroleum CCUS Demonstration Project, which has been injecting gas since 2014, now for a decade. The study of the migration state of CO2 after injection into the formation and the dynamic changes in the temperature and pressure of the storage layer are crucial for the safety assessment of carbon storage projects. This study, based on core, geological and test data from the CCUS demonstration area in Wuqi Oilfield, utilized the TOUGH2/ECO2N software platform to construct a three-dimensional numerical model to simulate the plume state of CO2 after injection into the formation and its impact on the temperature and pressure system of the formation. The results show that: CO2 tends to migrate to the upper part of the storage layer. In the initial stage of CO2 injection, the temperature around the well decreased sharply from 59.94 ℃ to 48.2 ℃ (a drop of 11.74 ℃), and the pressure increased from 12.6 MPa to a peak of 15.6 MPa (an increase of 23.81%). However, after 10 years of injection cessation, the temperature gradually returned to near the original value (60.2℃), and the pressure dropped to 12.8MPa (1.59% higher than the initial value). The temperature and pressure changes in the formation were more significant in the initial stage and gradually slowed down, eventually reaching a stable state. The impact of the injected CO2 fluid on the overall temperature and pressure of the formation was relatively limited. This study quantified the temperature and pressure response thresholds for CO2 storage in low-permeability oil reservoirs, revealing the coupling mechanism between the migration of supercritical CO2 and the energy balance of the formation, providing theoretical support for the design of storage schemes and the optimization of monitoring systems in similar reservoirs.
  • Sun Yanyang, Chen Hai, Cheng Kun, Lojain Suliman, Zhou Xiaohan
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 948-954. https://doi.org/10.20174/j.JUSE.2026.03.20
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Clarifying the settlement characteristics of small clearance tunnel excavation in urban areas is the premise of ensuring construction safety. In this paper, based on a rail transit subway station in Chongqing, the settlement characteristics of small clear distance tunnels (clear distance of 0.57 m) under different excavation sequences were explored by combining on-site monitoring and numerical simulation, and the optimal tunnel excavation spacing was determined. The results show that: (1) The removal of the temporary support will aggravate the deformation of the surrounding rock, and the deformation increment of the surrounding rock is about 1.43 mm~4.00 mm; (2) The main tunnel is excavated first, and then the final settlement curve of the No. 6 entrance and exit tunnel is "W" shaped, and the No. 6 entrance and exit tunnel are excavated first, and then the main tunnel is excavated, and the final settlement curve is "U" shaped, and the influence of the latter tunnels on each other is strong. (3) With the increase of the excavation spacing, the settlement curve changes from "U" shape to "W" shape, and the influence of tunnels caused by excavation is weakening. The surface load can lead to greater settlement, and when the load is large enough to a certain extent, it will affect the form of the settlement curve. (4) Taking the maximum settlement of the main tunnel vault as the control index, the effectiveness coefficient k of the excavation spacing is defined, and the optimal tunnel excavation spacing to meet the on-site construction conditions is 1.5 m under the premise of ensuring construction safety.
  • Yan Wei, Yue Bin, Xiong Tianfang, Geng Jia, Cheng Xuesong
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 955-964. https://doi.org/10.20174/j.JUSE.2026.03.21
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Shield method is widely used in tunnel construction due to its high construction speed, high economic efficiency and good safety characteristics, and the application of communication channels will also become increasingly widespread. However, the research for small-diameter shield tunnel stiffness efficiency is relatively limited. Based on a shield method contact channel project in Tianjin, with ABAQUS to establish a variety of tunnel and segment numerical calculation model, study the size of small diameter tunnel stiffness efficiency and influencing factors. The study found that: (1) The stiffness efficiency of the contact channel is not different from that of ordinary tunnel, the longitudinal stiffness efficiency is around 0.07, the transverse stiffness efficiency is between 0.5 and 0.6, the transverse stiffness improvement has a positive relationship on the longitudinal stiffness; (2) The segment assembly method has little influence on the efficiency of longitudinal stiffness, and extremely asymmetric deformation under the joint assembly; (3) The tension of the bolt, which has some influence on the efficiency of the longitudinal stiffness, but has little influence on the efficiency of transverse stiffness.
  • Chen Wuyuan, Wei Yongpeng, Guan Zhewei, Xue Chenxi, Ma Longxiang
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 965-972. https://doi.org/10.20174/j.JUSE.2026.03.22
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    At present, there is a lack of theoretical and quantitative analysis on the influence of shield tunnel posture on the upward floating phenomenon of shield tunnel segments under synchronous grouting, which makes it difficult to provide a reasonable shield tunnel posture adjustment plan for construction sites. By considering the influence of shield machine attitude on the shield-tail segments, based on the Euler beam theory, the Pasternak two-parameter foundation model, and other related mechanical theories, an analysis model for the influence of shield tunnel segment uplift under synchronous grouting is established, and the uplift phenomenon, deformation, and internal force responses of shield tunnel segments under different shield attitudes are analyzed. The results show that: Reasonable adjustment of the shield attitude can control the uplift phenomenon of segments during construction to a certain extent. The vertical position of the shield machine can be appropriately adjusted downward, and the pitch angle can be adjusted so that the tail of the shield machine rotates downward, which can alleviate the uplift phenomenon of the segments and reduce their internal forces. Unreasonable adjustment of the shield attitude will not only aggravate the uplift phenomenon of the segments but also significantly increase their internal forces. In actual projects, in order to control the uplift phenomenon, the vertical position of the shield machine can be appropriately adjusted downward, and the pitch angle of the shield machine can be appropriately adjusted to keep the tail of the shield machine rotating downward. After adjusting the shield attitude according to the research conclusions of this paper, the uplift phenomenon of the segments in the referenced project was alleviated. The research results can be used to guide the quantitative adjustment of shield tunneling posture when controlling the upward movement of construction pipe segments.
  • Li Ruixue, Miao Yichen, Lu Jiale, Su Hexian, Pan Wen
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 973-986. https://doi.org/10.20174/j.JUSE.2026.03.23
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Energy-absorbing anchor support is a key technology for controlling large deformations in high-stress surrounding rock. However, the relationship between the geometric parameters of conventional anchor energy-consuming components and their energy-consuming performance is unclear, which limits the optimization of these parameters. This paper designs a cone-head-sleeve damping structure as the core energy-consuming component of energy-absorbing anchors. Optimizing its cone angle (α) and sleeve wall thickness (t) parameters improves anchor support performance. Through indoor experiments and numerical simulations, we reveal the resistance characteristics and energy dissipation law of the damping structure during the slip process and establish a design method for optimizing the parameters. The results show that the SCATW (small cone angle and thick wall) design principle should be followed when using conventional materials to fabricate the damping structure. When the cone angle is between 0.5° and 5° and the sleeve wall thickness is between 4 and 7 mm, the damping structure has greater stabilizing resistance and better energy dissipation performance. A method of optimizing the geometric parameters to enhance energy dissipation is provided for engineering applications. with a 2° cone angle as the neutral point. In other words, when the cone angle α is between 0.5° and 2° or between 2° and 5°, respectively, increasing the cone angle and wall thickness improves the energy dissipation of the damping structure.
  • Wang Tiehang, Lü Siqin, Luo Yang, Yin Hongbo, Ruan Jiabin
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 987-994. https://doi.org/10.20174/j.JUSE.2026.03.24
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to enhance the stability of the coal mine shaft lifting vessel, the structural form of four-corner guide rails is often used to replace the interrupted part of the end tank passage. However, the current design method of four-corner guide rails is based on the design method of the end tank passage. However, due to the great differences in the structural forms and constraints between four-corner guide rails and the end tank passage, the reliability of the force of four-corner guide rails calculated by referring to the design method of the end tank passage is doubtful. In order to solve this problem, the acceleration time history curve of the four-corner guide rails under loading and unloading conditions is obtained through the field test of the force of the four-corner guide rails. By analyzing the time-history curve, the force of the four-corner guide rails is determined, and its characteristics and influencing factors are proved. After further regression analysis, the calculation method of horizontal force and vertical force of four-corner guide rails of coal mine shaft is established considering the lifting terminal load level. The calculation method based on the field test results is reliable and has important significance for improving the stability, safety and economy of the four-corner guide rails.
  • Li Xiaoyong, Yuan Feng, Li Zhi, Zhao Dongyang, Wen Yuanping
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 995-1002. https://doi.org/10.20174/j.JUSE.2026.03.25
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The special spatial structure of karst tunnels is due to the existence of beaded caves, which leads to the complexity of the stress-strain field of the surrounding rock, and the influence on the surrounding rock in the process of traversing the construction is unclear. Taking the Zhongliangshan Tunnel of Chongqing Rail Transit Line 27 is taken as the background of the project, and a three-dimensional numerical model is constructed based on the existing energy field visualization and analysis system to investigate the influence of the different locations of the caves on the stability of the tunnel surrounding rock. The results show that: The right side of the tunnel cavern has the most significant influence on the stability of the surrounding rock, and the vertical displacement of its arch foot is up to 50%, and the plastic zone is connected with the cavern area; when the cavern is located in the upper and lower part of the tunnel and the spacing is 5 m, the plastic zone is not connected; The energy evolution of the surrounding rock shows that the elastic strain energy around the support structure is released in a negative value, and the right cavern results in the formation of the dissipative energy density through the band; Based on the plastic zone through the criterion of the bead-like cavern safety protection thickness, it is proposed to be a safe and safe protection against the breakout. Based on the plastic zone penetration criterion, it is proposed that the thickness of the safety anti-protrusion of bead-like cavern should be more than 5 m.
  • Li Jichao, Guo Dong, Zhao Shiwei, Qiu Yunjun, Deng Feng
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1003-1012. https://doi.org/10.20174/j.JUSE.2026.03.26
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Large-diameter shield tunneling technology is widely applied in urban rail transit construction. However, its advancement faces significant challenges, including large excavation cross-sections, great burial depths, high thrust and torque requirements, and pronounced disturbances to rock and soil masses. This paper proposes a real-time monitoring method for the surrounding rock stress during shield tunneling, enabling continuous surveillance of stress variations. Based on the test results, the study analyzes the variation patterns of additional stress in the surrounding rock when operating under balanced conditions in an Earth Pressure Balance (EPB) shield. Furthermore, the paper discusses the magnitude and extent of influence of this additional stress. The results indicate that: The stress in the surrounding rock at the tunnel face gradually increases as the shield cutterhead approaches, reaching its peak value when the cutterhead arrives at the monitoring point. This peak value is close to the theoretical value of the additional contact pressure at the tunnel face. Both the circumferential and radial stresses around the tunnel exhibit a trend of initial increase, followed by a decrease, and finally stabilization. Specifically, the circumferential stress reaches its maximum value before the cutterhead arrives at the monitoring point, whereas the radial stress peaks when the cutterhead is directly at the monitoring point. Analysis of the variation pattern of additional stress reveals that the influence range of shield tunneling on the surrounding rock ahead is approximately 0.8 times the cutterhead diameter. Furthermore, both the maximum total additional stress and the residual total additional stress in the surrounding rock show a linear relationship with the distance from the monitoring point to the tunnel wall. Based on the theoretical value of the additional contact pressure at the tunnel face, calculation formulas for the maximum and residual total additional stresses are proposed. These stress variation patterns were further verified using data from shield tunneling in other sections. The research findings can provide a basis for controlling the stability of surrounding rock during large-diameter shield tunneling in urban rail transit projects.
  • Zhang Pei, Yang Chengru, Du Xiuli, Li Qianqian
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1013-1022. https://doi.org/10.20174/j.JUSE.2026.03.27
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The artificial freezing method is an effective technique for tunnel construction in water-rich sandy cobble strata. It is of great engineering significance to clarify the construction mechanical response when tunneling in frozen sandy cobble strata. Considering the complexity sandy cobble strata, the uniaxial compression indoor tests and numerical tests of frozen sandy cobble soil under different rock contents were conducted. The concrete plastic damage constitutive model was employed to describe the mechanical property of frozen sandy cobble soil, a numerical model and an analysis method for simulating tunnel excavation in artificially frozen sandy cobble soil were established. Then, the tunnel excavation simulations in frozen strata with different rock contents were carried out, and the construction mechanical responses of frozen strata and lining structure were analyzed. The results show that rock content has a significant impact on the stress-strain curve of frozen sandy cobble soil. The center of the lining bottom plate, the junction between the lining bottom plate and the straight wall are the dangerous locations. The plastic zone width in the center of the bottom plate increases gradually with increasing rock content. However, the stress distribution pattern of the lining is not affected by the rock content. The surface settlement curves under different rock contents all can be described by Peck's empirical formula. With the increase of rock content, the maximum value of ground surface settlement and the uplift of bottom plate all increase linearly, while the vertical displacement of arch shoulder and straight wall monitoring point increase firstly and then decrease.
  • Duan Xianbiao, Huang Yuanlou, Tang Jixiang, Zhang Qingzhao, Li Zongzheng
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1023-1033. https://doi.org/10.20174/j.JUSE.2026.03.28
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To investigate the mechanical characteristics of mud-water balance steel pipe jacking under deep-buried high-hydrostatic pressure strata, field monitoring was conducted on the pipe-soil contact pressure, pore water pressure, and pipe strain in the advanced drainage pipe jacking project at Wangjiazai Tunnel in Yunnan Province. The monitoring contact pressures were systematically compared with the design values in domestic and international standards. Additionally, numerical simulations were used to study the influence of secondary grouting pressure on the pipe. The results show that: (1) By comparing the monitoring contact pressure values under non-operational conditions with the calculated values in the standards, the earth pressure calculated using the calculation model and calculation parameters specified in the Japanese super large diameter PC method and Chinese T/CECS 246-2020 demonstrated superior agreement with the monitoring values. (2) During the jacking process, secondary grouting pressure significantly affected the contact pressure of the pipe, especially after entering the strong water-rich zone, where a rapid increase in grouting pressure led to sharp changes in the contact pressure on the grouting side. (3) Circumferential strain monitoring showed an alternating tension-compression pattern between opposite sides of the pipe, and axial strain deviated from the zero point due to the Poisson's effect and negative frictional forces, indicating an eccentric compression state during the jacking process. (4) After validation, numerical simulations were used to study the impact of secondary grouting pressure on the pipe. As the grouting pressure increased, the horizontal position yielded first in unilateral grouting, while in the bilateral grouting, the top and bottom of the pipe yielded first. The pressure control boundaries for unilateral and bilateral grouting were found to be 2.5 MPa and 0.6 MPa, respectively.
  • Lu Chenghua
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1034-1042. https://doi.org/10.20174/j.JUSE.2026.03.29
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Locking steel pipe pile cofferdams are widely used in deepwater bridge foundation construction in China. However, limited studies have examined their deformation characteristics under sloped foundation conditions. This study investigates the construction of a steel pipe pile cofferdam for the Xin'anjiang Bridge foundation in Jiande City, Zhejiang Province. A three-dimensional finite element model of the entire construction process is developed using Midas Civil and validated against on-site measurements. The analysis explores the deformation behavior of the cofferdam during construction and evaluates the effects of foundation inclination, construction techniques, and water flow velocity on lateral deformation. Finally, the deformation of the structure was effectively controlled by reasonably optimizing the size of the retaining structure. The results indicate that: Lateral displacement varies significantly due to differences in constraint conditions; Compared to the layered lowering method, the overall lowering approach mitigates cumulative deformation through the synergistic effect of each layer of purlins; Furthermore, the influence of the foundation inclination on cofferdam stability is more pronounced than that of water flow velocity; By optimizing the size of the maintenance system and balancing deformation control with material efficiency, cofferdam performance can be effectively improved. The findings of this study provide theoretical guidance and practical data for deepwater bridge foundation construction under similar geological conditions.
  • Tu Jinguang, Zeng Hongbo, Shi Zhou, Xiao Zhonglin, Xie Xiongyao
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1043-1055. https://doi.org/10.20174/j.JUSE.2026.03.30
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The stratum in coastal area is weak and rich in water with low bearing capacity. During the engineering construction by frozen method, the deformation and temperature fields of the structure and soil respond violently, which can easily lead to large-scale construction disturbances. Temperature effect parameters of soil and groundwater are obtained through field and indoor tests, a 3D thermo-mechanical coupled numerical analysis model for the whole process of freezing excavation thawing is established, the development and diffusion laws of the frozen curtain ring during the construction process is explored, the deformation and temperature evolution characteristics of soil and support structures are revealed, and the research results are finally applied and verified through engineering practice. The results show that: During the channel excavation, the soil within a range of 2.75~7.16 m from the channel center is compressed and compacted, and there are sharp temperature fluctuations. Within 15~25 days of heating and thawing, grouting reinforcement can be carried out on the soil within a depth range of 15.8 m below the ground surface. Especially within a range of 3.99 m from the channel center, the soil temperature rise faster than other regions. During the soil freezing process, the soil on the inner side of the freezing tube freezes faster than the outer side, and the soil on both sides of the channel freezes faster than the upper and lower sides. During the heating and thawing process, the upper and lower soil in the channel thaw faster than the two sides, and the thawing of the upper and lower soil in the channel mainly occurs during the continuous heating stage. During the thawing process, the deformation of the upper and lower parts of the primary support and secondary lining is significantly greater than that of the two sides, and there will be significant stress concentration at the corners of the steel rib bottom and the top and lower parts of the reinforcement mesh. Pre-embedded grouting pipes can be used to inject grout into the upper and lower soil 20 to 40 days after thawing begins, in order to slow down the development of structural displacement and stress concentration.
  • Liu Zhiqiang, Liu Guofeng, Chen Xueqi, Duan Shuqian, Pei Shufeng
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1056-1067. https://doi.org/10.20174/j.JUSE.2026.03.31
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Rock spalling is a common phenomenon of local rock mass failure in large underground cavern projects under high geostress, which seriously threatens the stability of engineering and construction safety. Relying on the underground powerhouse projects on both the left and right banks of the Baihetan Hydropower Station in Southwest China, a numerical evaluation method for the depth of rock spalling based on case back analysis is proposed through the use of numerical simulation, parameter inversion, field testing, and case investigation. Firstly, by collecting and organizing field data, the typical distribution pattern of spalling during the excavation of the roof arch and sidewalls of the Baihetan left and right bank powerhouses is statistically analyzed. Secondly, by taking rock displacement and loosening zone depth as target variables, key rock mechanics parameters in different rock spalling segmentations are obtained through the combination of a large number of field test results and the genetic-neural network algorithm (GA-ANN). Thirdly, the numerical simulation and the evaluation index of rock fracture damage (RFD) are used to analyze the range and depth of the brittle failure zone of the surrounding rock after excavation based on the hard rock degradation model (RDM) applicable to deep rock engineering, and the results are compared with the observed rock spalling damage on-site. The results show that more than 91% of the rock spalling depths correspond to RFD thresholds ranging from 1.35 to 1.50. Finally, the threshold was used to analyze the excavation of the layer Ⅲ of the underground powerhouses on the left and right banks of the Baihetan Hydropower Station. The results showed that the predicted accuracy of the rock spalling failure depth reached over 88%, proving the good applicability of the model in practical engineering, indicating that the RFD can be used to effectively predict the depth of rock spalling failure. This study can provide important support for predicting the depth of rock mass failure in deep underground engineering.
  • Long Haiping
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1068-1078. https://doi.org/10.20174/j.JUSE.2026.03.32
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to investigate the deformation law of pile foundation caused by the construction of a tunnel shield that obliquely penetrating into the pile foundation of adjacent existing bridges, the influence of the relative positional relationship between the bridge pile foundation and the shield tunnel is considered by taking the circular cross-section pile as an example. Using the Kerr beam foundation model, differential evolutionary algorithm and energy variational method, a prediction model for the deformation of adjacent piles caused by the diagonal tunneling of the shield tunnel is established. The engineering applicability of the prediction model is verified by the monitoring results at the construction site, and the discrepancy between the prediction results of different methods and the measured results is analyzed in terms of error. The results show that the model can predict the horizontal and vertical displacements of pile foundations caused by the construction of pile foundations at different locations, and can be used in two working conditions, positive underpassing and diagonal underpassing, with a wide range of engineering applicability, and the maximum error between the prediction results and the monitoring results does not exceed the 20% engineering error limit, which is of good applicability to the project. Under the condition of Kerr beam foundation model, the model in this paper also considers the relative position relationship between pile and tunnel, and effectively sets the shear layer parameter c. Compared with other methods, the prediction accuracy is significantly improved, which can provide theoretical guidance for the deformation control of adjacent piles caused by similar tunnel shield tunnel underpassing construction.
  • Wang Yanhua, Dai Fanfei, He Jianqun, Li Lichen, Cai Jihua
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1079-1089. https://doi.org/10.20174/j.JUSE.2026.03.33
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Based on the on-site data of vertical pipe jacking construction, a numerical model of vertical pipe jacking was established by using the particle flow code PFC2D. The failure process of the overlying marine soft soil layer and the jacking force were analyzed. The numerical model was verified to be reasonable by comparing the jacking force obtained from the field measurement and theoretical calculation. Then, the influence laws of the overburden depth, pipe diameter and pipe spacing on the failure process of the overlying marine soft soil layer were discussed from a mesoscopic perspective. The results show that: (1) During the pipe jacking process, the failure of the overlying soil layer mainly goes through the following processes: Compression at the bottom of the soil layer leads to shear failure → oblique cracks appear around the pipe → cracks expand to form voids → redistribution of the soil layer above the cracks → filling of the voids → new cracks are generated; (2) In the numerical simulation of vertical jacking, the failure surface of the soil layer has obvious nonlinear characteristics, and the error between the simulated maximum jacking force and the measured value is within 9%, which verifies the rationality of the numerical model; (3) With the increase of the overburden depth, the top failure range of the overlying marine soft soil layer increases linearly; (4) When jacking through a layer of silty clay interbedded with silt, the failure range of the surface soil is in a power function relationship with the pipe diameter, and its influence on the surface failure range is much smaller than that of the overburden depth; (5) The failure range of the overlying soil caused by jacking a single pipe was determined, and a reference formula for the minimum pipe spacing in pipe layout was given. The above numerical simulation results have reference value for the optimization of vertical jacking construction technology in the sea area and the layout of pipe spacing.
  • Ding Ziwei, Zhang Chenchen, Sun Wenliang, Dong Yunjun, Yu Maoqing
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1090-1102. https://doi.org/10.20174/j.JUSE.2026.03.34
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Layered sandstone is widely distributed in the surrounding rock environment of underground space engineering. The method based on cluster analysis is used to study the failure modes of layered sandstone under cyclic loading and identify failure signals, which can analyze the failure behavior of rock mass under complex stress conditions. Cyclic loading experiments were conducted on layered sandstone samples to monitor the acoustic emission signals of rock failure. The signals were classified using K-means++clustering algorithm, revealing the failure mode of layered sandstone under cyclic loading and identifying precursor signals related to rock failure. The results show that: The macroscopic failure of uniaxial loaded specimens is mainly shear failure, with a mixed shear tensile failure observed in specimens loaded three times per level, and tensile failure observed in specimens loaded five times per level; When the sample is loaded three times per level, micro cracks and large-scale cracks develop together in the early and middle stages. In the later stage, the stable expansion of large-scale cracks is the main trend, and ultimately the main cracks are formed by the penetration of micro cracks. When the sample is loaded five times per level, the large-scale cracks in the sample develop stably and quickly penetrate when approaching failure; The contour coefficients of cluster analysis under different loading conditions were 0.81 and 0.93, respectively, and three types of signals were classified: Type I signal was mainly tensile signal, accounting for 39.2% and 37.9% of the two groups of samples; Type II signal was mainly mixed failure, accounting for 43.3% and 57.3% of the two groups of samples; Type III signal was mainly shear failure, accounting for 17.4% and 4.8% of the two groups of samples, and the frequency of occurrence increased with the failure of adjacent rock samples. The research provides a scientific basis for the failure mode and failure prediction of layered sandstone.
  • Ning Xue, Ma Yongchi, Chen Kun
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1103-1114. https://doi.org/10.20174/j.JUSE.2026.03.35
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In order to explore the causes of subway construction accidents and prevent their occurrence, an analysis was conducted on 18 serious and above subway construction accidents in China from 2001 to 2022. Elements such as individuals, organizations, behaviors, equipment, and environment were extracted from subway construction accidents, and a network of subway construction accident causes was constructed based on this. The complex Network modeling and analysis tool Network X is used to measure network indicators such as node degree, shortest path, and intermediary centrality, combining the measurement results to conduct comparative analysis on key nodes and causal links in networks caused by serious accidents, major accidents, and catastrophic accidents. The results indicate that serious accidents are mainly caused by the superposition of management factors and equipment factors, as well as the superposition of management factors and illegal operation behavior factors. Major accidents are mainly caused by the combination of illegal operations and unfavorable geological conditions. Unfavorable geological conditions are the key factors for improving the level of accidents. Catastrophic accidents are mainly caused by the combination of organizational management factors and illegal operations, among which the lack of preparation of construction plans is the key factor for further escalation of accidents.
  • Du Mili, Li Xianglong, Wang Jianguo, Xu Jie
    Chinese Journal of Underground Space and Engineering. 2026, 22(3): 1115-1126. https://doi.org/10.20174/j.JUSE.2026.03.36
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To address the issue of noise interference in mining blasting vibration signals, a joint denoising algorithm combining the Crested Porcupine Optimizer (CPO), Improved Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (ICEEMDAN), and an enhanced wavelet thresholding method is proposed. First, the CPO algorithm was employed to globally optimize the key parameters of ICEEMDAN, adaptively decomposing the measured blasting vibration signals from an open-pit mine into a series of intrinsic mode functions (IMFs). A noise identification threshold was then constructed using multiscale permutation entropy (MPE) to screen out high-frequency noisy IMF components. These components were processed with the improved wavelet thresholding method and subsequently recombined with IMFs below the threshold to achieve denoising. Comparative experiments with CEEMDAN-MPE and ICEEMDAN-traditional wavelet thresholding methods demonstrate that the proposed method improves the average signal-to-noise ratio (SNR) by 39.33% and 19.93%, respectively, and reduces the root mean square error (RMSE) by 2~3 times across three sets of vibration signals. Furthermore, three-dimensional time-frequency energy analysis reveals that the main frequency energy distribution remains unchanged before and after denoising. These results indicate that the proposed method not only effectively eliminates noise interference but also fully preserves the main frequency energy characteristics of the original signal, demonstrating superior performance and engineering applicability.