Home Latest Articles
Latest Articles
  • Daolong Xiong, Shiqi Zhang, Yan Ouyang, Yan Tong, Zelin Liu, Kaiquan Chen, Peng Wang, Yinqing Zhou
    Underground Space. 2026, 27(0): 340-361.

    Shallow-buried urban road tunnels with shafts (URTS) have reduced traffic congestion in large cities. During fire scenarios, the backflows occur at those shafts far away from the fire source inhibiting smoke exhaust, but its rules have been unknown. A 400 m (length) × 12 m (width) × 5.5 m (height) physical model with 6-7 shafts over the ceiling is established using fire dynamic simulator software. The simulations are carried out after validation by both a small-scale experiment and a full-scale experiment. A total of 16 cases with 4 heat release rates (HRRs) and 4 spacing of fire source from the nearest unit shaft #1-1 (sf-us1), are designed. Results indicate that the smoke spreading length is nearly independent of HRR but increases with sf-us1. Ceiling smoke temperatures follow the power exponential laws, and the attenuation coefficients decrease with the increase of HRR and sf-us1. The farther away from the fire, the more likely the occurrence of shaft backflow. A good power exponential rule of the shaft negative mass flow rate is fitted out, and values of decay coefficient b2 range from 0.56 to 1.0. Based on dimensional analysis, a power exponential rule of the shaft dimensionless net mass flow rate is fitted for the exhaust shafts and a linear rule for the backflow shafts. The shaft neutral plane heights range from 1.4 m to 3.6 m for the exhaust shafts and 3.2 m to 5.4 m for the backflow shafts. Also, a linear rule is fitted. This study establishes the smoke backflow theory in URTS during fire scenarios and contributes to the tunnel fire protection engineering.

  • Zeyuan Sun, Cheng Zhao, Qinyuan Luo, Ran Huang, Ning Jia, Jinquan Xing, Jialun Niu, Chao Wei
    Underground Space. 2026, 27(0): 112-131.

    Burial depth and spacing are key factors influencing the stability of compressed air energy storage chamber groups. However, methods for determining safe burial depth and spacing in aligned multi-chamber systems remain underexplored. This study introduces a modified linear superposition method (M-LSM) for calculating elastic stress tensor fields in the surrounding rock around multiple chambers. The method analyzes the distribution of elastic stress under varying design parameters and surrounding rock conditions, identifying the most critical stress locations, which are the sidewalls affected by adjacent chambers, and providing deeper insights into stress concentrations caused by chamber interactions. The relationship between safe burial depth and spacing follows an approximately inverse proportionality, with the safety criterion of no plastic zones developing in the surrounding rock. A closed-form solution for safe burial depth and spacing is derived, which can be used to quickly determine the design parameters. The M-LSM is applicable to a wide range of internal pressurized chamber groups and borehole problems, capturing mutual interactions between adjacent chambers more realistically. Compared to finite element method simulations based on a practical engineering case, the results show errors that are typically negligible, validating the reliability of the analytical approach. This fully analytical method is mesh-free, iteration-free, and offers infinite resolution, making it highly efficient for both computations and practical applications. The closed-form solution derived from this method provides significant value for trend analysis, practical calculations, and engineering applications.

  • Zi-Yun Zhang, Fang-Le Peng, Chen-Xiao Ma, Yong-Kang Qiao
    Underground Space. 2026, 27(0): 58-71.

    Despite the thriving development of metro-led urban underground public space (UUPS) and its significant benefits and costs, there remains a critical research gap in understanding and evaluating its efficiency. This paper intends to improve the post-evaluation system of metro-led UUPS by proposing an efficiency evaluation framework based on data envelopment analysis. The public and the private sectors are taken as different coexisting decision-makers, and a pair of linear programming is built accordingly (with different assignments of discretionary and non-discretionary inputs) for each decision-making unit. The directional vector is calculated based on CRITIC weights to model the searching process of referential cases in terms of urban renewal. The empirical study of twenty metro-led UUPSs in central Shanghai reveals that (1) the proposed evaluation framework is feasible and discriminative, (2) the efficient form of metro-led UUPS in Shanghai is mainly limited to a compact pattern with a low proportion of pure public space, (3) the essential solution to promote efficiencies is closer cooperation between different parties, and (4) efficiency evaluation is crucial to avoiding the "the-more-the-better" type of development. The findings of this study are expected to shed light on the future planning and operation of metro-led UUPS.

  • Xing-Tao Lin, Deng Wang, Dong Su, Hui Zeng, Ruixiao Zhang, Xiangsheng Chen
    Underground Space. 2026, 27(0): 424-444.

    As shallow underground resources are depleted, urban development is extending to greater depths, necessitating a clear understanding of soil arching at various burial depth conditions. Laboratory trapdoor tests equipped with embedded soil-pressure cells and digital image correlation captured the ground-reaction curve and soil deformation. The results reveal pronounced discrepancies between shallow and deep burial. In shallow conditions, soil arching undergoes a "failure-reconstruction" process: soil pressure plunges, then rebounds to stability. In deep strata, the arching forms rapidly and attains stability almost immediately after the minimum pressure is reached. Shallow tests generate several horizontal displacement bands rising to 4.8B (B, trapdoor width); deep tests yield one stable band, with its influence height reduced to about 3.0B. Vertical displacement above the trapdoor evolves through 'triangular-tower-para bolic" stages to 4.3B in shallow tests, but follows a persistent parabolic profile limited to 2.7B in deep tests. Additionally, shear bands under deep conditions form at smaller angles and are more vertically oriented. These findings expose the fundamental differences in deformation mechanisms between shallow and deep burial and provide quantitative criteria for depth zoning in urban underground space development.

  • Oveis Farzay, Marilena Cardu
    Underground Space. 2026, 27(0): 72-91.

    Accurate TBM performance estimation is essential for effective tunnel design and planning. This study introduces a one-dimensional (1D) estimation model that estimates thrust, torque, power, cutterhead speed, and tool count using only excavation diameter. The model was developed across four TBM types—open, single shield (SS), double shield (DS), and earth pressure balance (EPB)—to isolate the influence of diameter from other variables. Validation against existing models and a 52-case independent dataset confirmed strong correlations: torque scales with the cube of the excavation diameter (R2 = 0.89 for EPB), power grows faster than linearly (R2 = 0.83 for EPB), thrust increases supra-linearly (R2 = 0.79 for EPB), and cutterhead speed decreases with diameter (R2 = 0.87 for open TBM). Tool count grows proportionally. A reliability matrix compares model accuracy and data support, aiding selection based on both fitness and robustness. This 1D model offers fast, consistent estimates for early-stage assessments. While it excludes detailed geological input, it is suited for feasibility studies and preliminary design. Future work will incorporate additional ground and machine parameters and extend validation across a broader range of tunneling conditions to enhance generalizability.

  • Haitao Yu, Yibo Wei, Yong Yuan, Guoliang Li, Hehua Zhu
    Underground Space. 2026, 27(0): 175-193.

    An earthquake with a magnitude of 6.6 (MW) occurred in Menyuan County, Qinghai Province, China, on January 8th, 2022. The Daliang Tunnel, which traverses the seismogenic fault, was severely damaged during this seismic event. The seismic damage investigation of the tunnel is introduced, and the damage characteristics along the tunnel are also presented. It is found that the damage severity of the tunnel is highly correlated with the distance to the fault. Damage modes for different tunnel zones (cross-fault zone, portal zone, and ordinary zone far from the fault and the portal) are quite distinct. Based on the understanding of seismic damage to the Daliang Tunnel in 2022 Menyuan Earthquake, as well as other damaged tunnels during the 2008 Wenchuan Earthquake, the critical influence factors for damage to cross-fault tunnels are discussed. From the seismic investigation, coupling effects of strong ground motions and fault dislocations are highlighted, which pose significant risks to cross-fault tunnels and should be considered in the seismic design for such tunnels. The deformation joints of the tunnel could help the tunnel adapt to the deformation caused by fault dislocation, and thus protect the main structure, but inversely an added local deformation may be expected around the deformation joints. The instability of the overlying slope may have contributed to the local amplification of the fault dislocation along the Daliang Tunnel. Finally, the requirements of seismic resilience for cross-fault tunnels are proposed, and detailed suggestions are provided to enhance the seismic performance of such cross-fault tunnels.

  • Rongmin Bai, Chuan He, Guowen Xu, Bo Wang, Xu Chen, Gaoyu Ma, You Zhou
    Underground Space. 2026, 27(0): 150-174.

    Reliable data sources are essential for intelligent tunnel construction, yet on-site data are often insufficient to meet sample require-ments. Previous numerical modeling studies have seldom considered the combined effects of different excavation methods' spatial effects and the distribution characteristics of joints. This paper develops a method to construct a stability database for jointed rock tunnels with primary support systems using a computational framework combining the finite difference method-discrete element method (FDM-DEM). The framework constructs a 2D model using the Mohr-Coulomb criterion and a 3D model with the Hoek-Brown failure criterion, enabling the stress release process to accurately replicate the influence of joint distribution features and excavation space effects in the 2D calculations by utilizing the longitudinal deformation profile parameters of the bench sections and fine-grained ground reaction curves. The computational circle is determined by grid research and data analysis, while the performance differences of various primary support components and their correlations with surrounding jointed rock are analyzed using the control variable method. The validity of the framework is initially confirmed by case comparisons and macroscopically validated using the Mantel test and Spearman analysis on the constructed simulation database—containing tunnel construction information, joint distribution, rock mechanics parameters, and stability indices—thereby establishing a reliable foundation for machine learning and transfer learning applications.

  • Ruikun Wang, Gang Zheng, Huayang Lei, Xuesong Cheng, Eng-Choon Leong, Yetao Ji
    Underground Space. 2026, 27(0): 24-44.

    This study investigated the long-term settlement behaviour of piled buildings induced by shield tunnelling in soft ground conditions within urban environments. By integrating a detailed case study with advanced numerical modelling techniques, this study provided a nuanced understanding of the interactions between tunnel construction and existing pile foundations. Central to the investigation is the role of soil consolidation, which significantly contributes to the settlement of piled buildings. To address this, this study emphasizes the critical need for the precise calibration of tunnelling parameters such as face pressure and grouting pressures. These parameters are meticulously controlled to mitigate the adverse effects on nearby piled buildings, ensuring their stability and integrity. It is established that an optimal face pressure, set at 90% of the lateral earth pressure, consistently minimizes the settlement of piled buildings, primarily due to the minimal reduction in the pile toe resistance observed near the tunnel. Similarly, the ideal grouting pressure was identified to be within the range of 120%-160% of the vertical earth pressure, with the smallest building settlement and decrease in pile toe resistance observed at a grouting pressure of 150% of the overburden pressure. This finding elucidates the load transfer mechanism within piled buildings. This study further demonstrated that the settlement induced by the second tunnel excavation is smaller than that caused by the first tunnel excavation owing to the sheltering effects of the adjacent first tunnel and pile foundations. During the consolidation phase following tunnel excavation, the settlement caused by the second tunnel is smaller than that caused by the first tunnel, which is attributed to the dissipation of the negative excess pore pressure around the first tunnel, leading to soil volume expansion. These insights not only validate the effectiveness of the numerical model but also contribute significantly to the field of geotechnical engineering by providing actionable guidelines for future tunnelling projects.

  • Yuyang Cao, Xiongyao Xie, Kun Zeng, Yangbin Zhang, Haiyang Tian, Jian Yao, Junliang Zong
    Underground Space. 2026, 27(0): 386-405.

    Super-large-diameter shield tunneling inevitably induces deformations in the surrounding soil and nearby existing tunnels due to ground-tunnel interactions. This study developed and validated a numerical model to simulate these interactions in typical soft soil strata in Shanghai, with a focus on stress and displacement responses during the undercrossing of an existing tunnel by a new super-large-diameter shield tunnel. The study identified an incomplete soil arching (ISA) effect and proposed methods to delineate the ISA, loosened, and compaction zones, categorizing the influenced areas into reinforced, stable, and safe zones. Parametric analyses examined the influence of tunnel spacing (S) and volume loss ratio (V) on ground deformation, loosened zone height, and existing tunnel deformation. Results indicate that greater volume loss ratios and smaller tunnel spacings amplify ground settlement, while the loosened zone height is affected by both the volume loss ratio and the stratigraphic boundary. Among the considered scenarios, a volume loss ratio of 0.2% minimizes the loosened zone height across various spacings. Changes in the convergence of the existing tunnel occur in two phases, characterized by rapid changes (S/D of 0.1-0.3, where D is the diameter of the newly constructed tunnel) and gradual changes (S/D of 0.3-0.7). To mitigate adverse effects on the ground and the existing tunnel, it is recommended to maintain the volume loss ratio below 0.2% and the tunnel spacing over 0.3D. Additionally, reinforcing the loosened zone is advised to enhance the stability of the existing tunnel.

  • Pengfei Chen, Quan Jiang, Yong Xia, Zhijun Liu, Long Li
    Underground Space. 2026, 27(0): 236-258.

    Actual seismic events typically exhibit mainshock-aftershock sequence characteristics, and source characteristics have a significant impact on cavern response. Currently, the influence of near-fault mainshock-aftershock sequences (NFMA) and far-field mainshock-aftershock sequences (FFMA) on underground caverns is generally ignored. This study aims to establish a framework for evaluating the dynamic response characteristics and seismic fragility of large-scale underground caverns under NFMA/FFMA. The response laws of residual displacement and rock fracture degree of cavern under NFMA/FFMA are comparatively studied, and the failure probability of different damage states is quantified by the fragility function. The results show that the surrounding rock of underground caverns exhibits significant cumulative damage effects and non-uniform failure characteristics under mainshock-aftershock sequences. Aftershock fragility is strongly related to the mainshock-damaged state for underground caverns. The collapse probability of underground caverns after 0.9g aftershocks in NFMA increased from 0.76% in slight damage to 21.12% in moderate damage and 53.51% in severe damage. This study can provide a probabilistic basis for seismic design, aftershock risk warning, and post-earthquake emergency assessment in underground engineering.