Chinese  |  English
Home Table of Contents

15 March 2025, Volume 65 Issue 3
    

  • Select all
    |
    Editorial
  • Journal of Tsinghua University(Science and Technology). 2025, 65(3): 413-413.
    Abstract ( ) Download PDF ( )   Knowledge map   Save
  • Review
  • ZHOU Wendong, CUI Yanwei, WANG Hetang, REN Gehui, CUI Xinyue, WANG Hao, SHEN Aojie
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 414-432. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.017
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Significance] Coal remains the most essential fossil energy source in China, with production from underground coal mines accounting for more than 80%. In underground coal mining operations, dust is a pervasive and hazardous material that can significantly compromise the safety and health of miners. High dust concentrations are associated with an increased incidence of pneumoconiosis and a heightened risk of catastrophic events, such as coal dust or gas explosions. The primary source of dust is the cutting processes of excavators and roadheaders during the extraction of coal and rock. Therefore, it is essential to comprehensively elucidate the mechanism of dust generation during cutting and heading operations for effective dust control in underground coal mines. [Progress] This paper presents an overview of the latest research developments in the field of dust generation, with a particular emphasis on three key areas: the behavior of dust generation during cutting, research methodology, and influencing mechanisms. First, regarding the behavior of dust generation during cutting by excavators and roadheaders, researchers have proposed several theoretical models to elucidate both the fragmentation processes of coal and rock bodies and the generation of dust under the influence of cutting. These models are based not only on the principle of energy conversion but also on the influence of the geometry of the cutting pick during the dust generation process. This provides a solid theoretical basis for understanding the physical nature of dust generation. Regarding the research tools employed, researchers have simulated the dust generation phenomena when mining machinery cuts coal and rock bodies through physical experiments conducted on self-designed experimental platforms. Researchers have also conducted numerical simulations using finite element or discrete element methods. These advanced experimental techniques elucidate the actual cutting conditions and offer a robust analytical tool for investigating the dust generation mechanism in depth. Additionally, this study provides a comprehensive analysis of the mechanisms influencing dust generation during cutting, examining both internal and external factors. These factors include the physicochemical properties of coal and rock, such as coal rank, moisture content, pore characteristics, strength, and brittleness, as well as the parameters of the cutting conditions, such as cutting depth, advance speed, drum rotation speed, and the morphology and arrangement of picks. [Conclusions and Prospects] Current research on the dust generation mechanism during cutting reveals several contradictions. Existing models often rely on simplified assumptions, neglecting the anisotropy of coal and the actual cutting conditions in the field, leading to discrepancies between the calculated results and experimental observations. Moreover, existing experimental platforms struggle to accurately replicate the motion of the cutting pick during actual operations. Although many studies have focused on the properties of coal and rock and dust characteristics, some of the conclusions are conflicting. Future research should prioritize the construction of full-scale experimental platforms and the development of high-precision monitoring technologies. Comprehensively investigating the dust generation characteristics of complex coal seams and quantifying the energy conversion mechanisms during the cutting process are crucial. These efforts are essential for improving the efficiency of cutting operations and achieving more effective dust control.
  • Research Article
  • XIANG Yunfei, LUO Yiming, NING Zeyu, LIU Yuanguang, YANG Zuobin, LI Zichang, LIN Peng
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 433-445. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.014
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Object] Hydropower underground engineering encounters significant safety management challenges owing to overlapping construction activities, diverse process stages, and dynamic resource flows. This involves multidisciplinary safety tasks, such as safety hazard identification and rectification, emergency response, and regulatory compliance checks, which require specialized domain knowledge. In this context, safety management knowledge is intricate, such as expert experience, patterns and characteristics, and management codes, and is dispersed across multimodal data formats, including text, tables, and images. Efficient extraction of these multimodal data sources can significantly enhance data utility and support intelligent safety management. However, owing to the diverse nature of data formats, the complexity of the knowledge system, and the various management scenarios, current research struggles with limited knowledge sources, acquisition difficulties, and poor generalization. [Methods] This study proposes a method of constructing a multimodal knowledge graph (KG) for safety management in hydropower underground engineering. (1) A large-scale, high-quality, multisource heterogeneous dataset is built from safety hazard identification and rectification records, regulations, and images. (2) Knowledge modeling employs top-down and bottom-up approaches to define entities, relationships, attributes, and events pertinent to safety management in hydropower underground engineering. (3) The entity and relationship information from text data is obtained using a knowledge extraction method that uses a large language model (LLM) tuned with domain knowledge, enriched by specific examples for each entity type to handle small sample sizes. This approach uses demonstrations to provide the model with prior knowledge. (4) Instance segmentation is used to annotate safety hazard images. The entities identified in the images are then converted into vectors. Image and text data are linked based on semantic similarity. Image data are integrated into the textual KG, enabling the transformation from multimodal data to multimodal knowledge. (5) The multimodal KG is stored in Neo4j, an open-source graph database management system. (6) A scenario-specific knowledge acquisition method addresses the specific needs of safety management scenarios, integrating KG with LLMs to enable retrieval-augmented generation and interpretable knowledge reasoning. [Results] (1) This paper collected more than 120 000 safety hazard records, 30 regulatory documents, and 300 000 images of safety hazards. Leveraging these comprehensive data, this paper constructed a large-scale, high-quality, multisource heterogeneous dataset specifically designed for managing safety in hydropower underground engineering projects. (2) Taking a hydropower underground engineering project as an example, the constructed multimodal KG was applied to intelligent recommendations for safety hazard rectification and compliance checks. (3) The workflow for generating intelligent recommendations for safety hazard rectification measures involved the following steps. After users input safety hazard information, the scene-KG was extracted from the multimodal KG and fed into an LLM to generate appropriate rectification measures. (4) Based on the scene-KG, an inference retrieval method extended neighboring nodes and constructed inference-KG for compliance checks. By integrating inference-KG with an LLM, the system retrieved relevant content from regulatory documents based on user input. [Conclusions] The proposed method effectively extracts and applies domain knowledge from multimodal data in the context of safety management in hydropower underground engineering. It also successfully applies domain knowledge for safety management. The results serve as a reference for transitioning infrastructure construction safety management from a data-driven approach to a knowledge-driven approach.
  • AN Ruinan, LIN Peng, WANG Xin, LI Zichang, LIU Yuanguang, HE Pinjie
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 446-454. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.016
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Intelligent ventilation on demand is crucial for ensuring environmental safety in underground caving groups and for the high-quality construction and development of hydropower projects. Ventilation systems for large underground caving groups during construction frequently exhibit complex three-dimensional layouts, different air loads across regions, and dynamic demand under varying regulation conditions. [Methods] To achieve spatial node extraction, branch correlation decoupling, and stable joint adjustment of complex flow fields, this paper examines the development characteristics of fluids under construction ventilation in extensive spatial structures. It demonstrates the necessity of constructing a graph structure based on the ventilation flow characteristics for analyzing and adjusting ventilation system parameters. The regional modeling theory is discussed, detailing the principles and methods of node extraction for one-dimensional tube bundle fluids (network) and three-dimensional spatial flow field elements (field). Among these, the area where fluid parameter information changes along the main airflow direction employs network node extraction, while the regions with multi-directional complex flow paths utilize the three-dimensional field node extraction method. Virtual branches address the network-field coupling problem, utilizing the nodal pressure approach. This method treats the nodal pressure as the unknown variable and airflow deviation as the assessment criterion. Nodes with known pressure values serve as reference nodes for solving the pressure at all network nodes, and are further assigned to field simulation boundaries. By numerically simulating the three-dimensional spatial flow field, the virtual branch air flow rates are iteratively fed back into the air network calculation for a coupled solution. This paper also introduces the node-property-edge triplet, which effectively reflects the structure, performance, and behavioral characteristics of nodes. Furthermore, to optimize the ventilation coordination efficiency, a hypergraph structure for joint adjustment, with edges as the analysis object, displays the coupling interactions between the ventilation branches and loops. Considering the joint adjustment sensitivity, an optimal resistance control method is proposed, which involves constructing target and response node sets, setting response efficiency constraints, and optimizing to form a ventilation adjustment plan. An intelligent ventilation coordination platform integrates the resistance control model of coupling interactions, including modules for network design, ventilation design, field-network integration, loop generation, and optimization analysis. Within this framework, the network design module is dedicated to reconstructing the physical model of the ventilation system, while the ventilation design and field-network integration modules are used to assign basic fluid characteristic parameters of ventilation to the established model. The loop generation and optimization analysis modules are employed for solving the overall wind network parameters, including air volume, air pressure, and wind resistance. [Results] The field-network coupling method using nodal pressure eliminated the need for loop identification and effectively addressed the interdependent coupling between network nodes and flow field boundaries. The intelligent ventilation coordination platform was integrated with online environmental monitoring devices to automatically gather critical ventilation environment parameters, thereby enabling real-time calculations of the ventilation system based on environmental monitoring data and providing 3D visualization and early warning capabilities. [Conclusions] The ventilation design parameters of an engineering project are used to implement targeted air volume control deployment. The integrated control system exhibits high responsiveness. On the premise that the air volume of each unit meets the threshold requirements, the air volume adjustment efficiency of the target unit and the overall stability of the air distribution network can always fulfill the specified requirements. The results indicate a timely and stable system response and can provide a reference for similar projects.
  • CHENG Huihang, MIAO Ruifeng, YANG Xiujun, PAN Ronglian, CHEN Junfeng, ZHONG Maohua
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 455-468. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.002
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] With the continuous development of social and economic levels and the increasing demand for a higher quality of life, the scale and quantity of transportation tunnel construction in China continue to expand. The risk and harm of tunnel fires are increasing. Previous research has focused on single or bifurcated tunnels, lacking experimental research on fires for long-distance tunnels with double-hole tunnels. [Methods] This study focuses on a typical fire scenario of a double-hole long-distance highway tunnel, conducting full-scale experiments to evaluate the diffusion characteristics and temperature distribution of smoke under natural ventilation conditions and obtains basic data on double-hole tunnel fires. From the full-scale results, a computational fluid dynamics model was built, and further numerical simulation analysis was conducted to discuss the ventilation linkage mode of double-hole tunnels under fire conditions. [Results] Smoke diffusion under different fire conditions was characterized by analyzing key parameters such as airflow velocity, smoke temperature distribution, and smoke diffusion time. The smoke control effects under different mechanical ventilation modes were compared using computational fluid dynamics tools. The results showed that: (1) Under natural ventilation conditions, when a smaller fire source power (eight oil pans) was used, the highest temperature point upstream of the fire source appeared at a height of 3 m instead of at the ceiling, and the temperature in the area between 3 and 4 m was higher. As the power of the fire source increased to 12 oil pans, the increase in thermal buoyancy increased the temperature to the highest point, approaching 3 m. (2) Mechanical ventilation reduced the doping effect of natural wind, stabilizing the distribution of the smoke layer upstream of the fire source, and the temperature upstream of the fire source was vertically distributed with the height gradient. Because of the opposite direction between mechanical and natural ventilation, the reduction in fresh air doping weakened the cooling effect of ventilation, resulting in a higher temperature under mechanical ventilation than under natural ventilation and a maximum temperature increase of 5-10℃. (3) For the flame inclination angle, as the combustion intensified, the thermal buoyancy gradually increased, and a larger plume buoyancy led to a smaller flame inclination angle. For flame length, as the heat release rate increased, the buoyancy of the plume increased, resulting in increased flame volume and length. (4) Based on the numerical simulation, the smoke control effects of single tunnel ventilation and left and right line linkage ventilation modes were compared. Under the set fire source power and position, the mode of smoke exhaust at end A and air supply at end B of the left tunnel while using the right tunnel for natural ventilation achieved the greatest benefits. [Conclusions] Smoke diffusion, temperature distribution, and fire source morphology in tunnel fires are discussed, and the ventilation mode for smoke control in tunnel fires is presented. The optimal ventilation mode under the set operating conditions is obtained from numerical simulation. The experimental results can provide data support and a technical reference for the smoke control design of tunnel projects with similar structures.
  • WU Le, YUE Shunyu, CHENG Huihang, ZHONG Maohua
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 469-478. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.001
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] In fires originating from the basement or ground floor of high-rise buildings, the path of toxic smoke often intersects with the escape routes of the building occupants. This especially poses a significant risk when the smoke spreads vertically at an accelerated pace, threatening escapees and emergency rescue personnel. This study uses comprehensive field experimental data to reflect the progression of a fire to its fullest extent. In addition, the patterns of smoke diffusion and sedimentation were studied by observing experimental phenomena. [Methods] A fire experiment was conducted on-site in the basement of a high-rise building. Four different fire scenarios were created at various locations: a room corner, the front room of the stairwell, and the stairwell leading from the basement to the ground floor. The horizontal and vertical smoke temperatures in the fire area, adjacent rooms, and stairwells were measured and dispersed in real time. In addition, fluctuations in the air velocity during the experiment were recorded. [Results] The results showed that in the fire area, there was a significant accumulation of smoke in the room after the fire. The smoke temperature maintained a relatively stable vertical gradient. However, the stairwell and its front room, which are connected to multiple areas, allowed a large amount of smoke to diffuse in neighboring zones, which resulted in reduced smoke accumulation and stable smoke stratification. Under experimental fire powers of 0.125 and 0.250 MW, the vertical temperature distribution in the fire area above a height of 3.50 m fluctuated significantly. Conversely, the smoke temperature below a height of 3.50 m remained consistent. In the room adjacent to the fire site, smoke initially spread across multiple interconnected areas. Despite this, smoke accumulation and settlement effects could still form further in the adjacent room and the stairwell's front room. In the stairwell's front room, which has a lower spatial limitation, smoke stratification was evident. Here, the spread of low-temperature smoke and air to distant rooms was more serious, resulting in a noticeable increase in the flue gas temperature in the lower space. In the stairwell, the smoke spread rapidly vertically. In the 0.125 MW fire scenario, smoke could reach the fourth or fifth floor area within 300 s, reducing the smoke layer height at the stairwell corner to 2.00 m, which poses a significant threat to personnel evacuation. [Conclusions] Through this analysis of smoke settlement characteristics in the fire area, the law of smoke spread and settlement characteristics in adjacent rooms and stairwells during 0.125 and 0.250 MW field fire experiments in the underground space of high-rise buildings is better understood. This knowledge of fire occurrence conditions, spread patterns, and smoke flow characteristics provides robust data support for smoke control design in the underground spaces of high-rise buildings. This approach improves the efficiency and effectiveness of fire response measures. This study successfully achieves its goal of analyzing smoke diffusion ranges and smoke layer heights.
  • HUANG Ran, ZHU Shiyou, HE Mengchen, LI Ruoyu, GE Xinru, WANG Qiao, CHEN Juan, LO Jacqueline T. Y., MA Jian
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 479-494. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.013
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] In case of catching fire, a train moving in a tunnel section might lose power and come to an emergency halt while approaching the next station for emergency rescue. Predicting the distribution of smoke and temperature in tunnels under such fire scenarios is difficult because of the influence of moving trains. This difficulty in prediction might seriously threaten the safety of passengers and impede metro evacuation management. [Methods] This study considers influencing factors, including train speed, fire source location, and fire heat release rate, to solve this problem, and 75 different fire scenarios are designed and simulated. Train movement in the simulated scenarios is realized using the equivalent piston wind method. The simulation results of the smoke and temperature distributions collected using sensors near the tunnel ceiling are used to construct a dataset for deep learning. Accordingly, a deep learning model comprising long short-term memory networks, a convolution (Conv) module, and a deconvolution (DeConv) module is then proposed for rapid prediction of temperature distribution in tunnels under moving train fire conditions. The train speed, train braking time, and temperature time-series information from the sensors together are fed as inputs to the model. [Results] The results indicated that: (1) Under various train movement states, the model was able to predict the temperature distribution of the lateral evacuation platform in the tunnel 30 s in advance using the current sensor data, with a mean absolute error (MAE) of only 2.2 ℃ and a mean absolute percentage error (MAPE) of 4%, indicating high accuracy. (2) In a stark contrast with the week-long time taken to obtain temperature distribution in a fire dynamics simulator (FDS), this deep learning model could make prediction within only 0.08 s, hence representing a computational efficiency improvement of four orders of magnitude versus the computational fluid dynamics method. (3) Validation with fire scenarios in none of the training, validation, and test datasets resulted in model MAE and MAPE values of 3.1 ℃ and 5%, respectively, indicating a strong generalization ability. (4) Considering the possibility of sensor failure within tunnels, this study investigated the influence of simulated sensor failures on the model's prediction accuracy by varying sensor spacing. The model continued to exhibit an effective predictive ability even at a sensor spacing of 4.00 m. At 8.00 m sensor spacing, the model's errors were larger albeit at very few time frames (with a maximum MAE lower than 10.0 ℃ and a maximum MAPE lower than 15%). However, for other sensor spacing cases, the model's MAE and MAPE were less than 5.0 ℃and 10%, respectively. Hence, it could be concluded that the model has strong robustness. [Conclusions] This study constructs a comprehensive dataset for a tunnel with moving train fire conditions using an FDS and leverages advanced deep neural networks to completely exploit the extensive information within the dataset, ultimately resulting in a high-precision, robust model for rapid prediction of temperature distribution in tunnels under moving train fire conditions. These advancements are highly important for effective emergency management and response planning in tunnels under these challenging conditions.
  • DENG Zhiyun, LIN Peng, MA Xiang, XIA Yong, LI Zichang
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 495-508. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.015
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] To address the issue of cracking in the lining structures of deep-buried tunnels, this paper proposes the use of basalt fiber-reinforced concrete (BFRC) to improve the load-bearing capacity of lining structures. [Methods] Tension and compression tests were conducted on BFRC specimens with varying volume basalt fiber fractions ranging from 0 to 0.5%. The optimal fiber content was determined, and the concrete damage plasticity (CDP) model parameters for plain and fiber-reinforced concrete with the optimal fiber content were validated. Then, numerical simulations were employed to create an integrated bearing model of surrounding rock-initial support and a secondary lining support model. The use of solid and structural elements and a secondary lining support characteristic curve (SCC) for deep-buried tunnels were obtained, revealing the crack propagation characteristics of the concrete. A quantitative analysis was conducted on the effects of the reinforcement ratio, secondary lining thickness, and fiber-reinforced concrete on the normal and ultimate state load-bearing capacity of the secondary lining. [Results] (1) Compared with those of plain concrete (B0), the maximum increases in the axial tensile strength and splitting tensile strength of BFRC with a fiber volume fraction of 0.2% were 12.81% and 14.79%, respectively. Furthermore, a maximum enhancement of 31.68% in the flexural strength of BFRC was noted when the fiber volume fraction was increased to 0.5%. The optimal fiber content was 0.2%. (2) The stress-strain curve of the BFRC could be fitted using peak compressive strength, peak compressive strain, and compressive shape parameters. The compressive shape parameter values for B0 and B0.2 were 6.50 and 3.00, respectively. The tensile stress-strain curve could be fitted using the peak tensile strength, peak tensile strain, and tensile shape parameter, with tensile shape parameter values for B0 and B0.2 being 3.00 and 1.86, respectively. (3) The CDP parameters for plain concrete and fiber-reinforced concrete accurately simulated the peak tensile and compressive strengths as well as the shapes of the tensile and compressive stress-strain curves. For compressive stress-strain curves, the error between numerical simulation and experimental fitting values at 0.50% compressive strain was 2.48% (2.01%) for B0 (B0.2). For tensile stress-strain curves, the error at 0.04% tensile strain was 4.08% (1.68%) for B0 (B0.2). (4) The SCC curve of the secondary lining exhibited rapid linear growth initially, slow growth in the middle, and a nearly horizontal trend in the later stages with increasing displacement. For class V surrounding rock, the secondary lining crack width showed slow linear growth in the initial stage and rapid linear growth after reaching approximately 0.10 mm. Higher reinforcement ratios effectively delayed crack propagation in the early stage, although increasing the reinforcement ratio beyond 0.6% or 0.8% was not economically reasonable. (5) Increases in reinforcement ratio and lining thickness resulted in almost linear increases in the 0.30 mm crack load and ultimate state load-bearing capacity. For every 0.1% increase in the reinforcement ratio, 0.30 mm crack load increased by an average of 5.40%. In addition, for every 0.10 m increase in the secondary lining thickness, 0.30 mm crack load increased by an average of 11.18%. Fiber addition considerably enhanced concrete resistance to crack propagation, especially in the early stages, increasing the initial cracking load by 38.64% and the 0.30 mm crack load by 5.54%. [Conclusions] This study provides theoretical and practical guidance for designing deep-buried tunnel lining structures and serves as a reference for applying fiber-reinforced concrete in secondary lining structures of deep-buried tunnels.
  • PENG Yuqi, LI Chao, YANG Ruihang, WANG Dachuan, ZHOU Tiejun
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 509-520. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.011
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Deep underground spaces present unique challenges owing to their strong closure, long evacuation distances, and high evacuation difficulty. In disaster situations, evacuees must navigate long, confined horizontal channels, increasing the risk of overcrowding and trampling, which can severely compromise safety. Existing national standards, industry norms, and research results involving underground space safe evacuation primarily target shallow underground spaces, rail transit, and civil defense areas, which are difficult to adapt to the evacuation requirements of deep underground spaces. This highlights the urgent need for specific research into the evacuation design of deep underground spaces. Given the scarcity of relevant guidelines and studies and the potential for congestion during large-scale evacuation in these horizontal channels of deep underground spaces, there is a critical need to optimize the design of horizontal channel widths. [Methods] A lattice-like deep underground space model was developed based on the underground space planning at Nanyang Technological University in Singapore. This model connects chambers in series through horizontal channels. In the first layer of the model, chamber exits and channels are identified, with internal chamber channels, evacuation channel intersections, and vertical evacuation facilities being abstracted as source, intersecting, and terminating nodes, respectively. The shortest evacuation paths from the source to the terminating nodes are considered as directed sides in the network. Then, the topology network model is established. The betweenness centrality index from complex network analysis is used to assess the importance of nodes. Based on these calculations, horizontal channels are divided into three levels: first, second, and third. This classification helps construct a graded road network model. Current specifications guide the assignment of basic widths for these channels at each level. The study uses GoAhead, a self-developed pedestrian evacuation dynamics simulation software, to simulate evacuation scenarios. By setting different evacuation widths for various working conditions, the simulation evaluates evacuation density and time as people move through key node areas in real time, providing insights into evacuation effectiveness. [Results] The evacuation simulation results showed the following: (1) When the width of the third-level horizontal channel was 6.0 m or less, increasing its width effectively reduced evacuee density, alleviated congestion, and shortened the evacuation time. However, beyond 6.0 m, further width increases did not affect the evacuation time. (2) When the width of the second-level horizontal channel was 9.0 m or less, the maximum evacuation density and time were reduced as width increases, thus effectively improving the evacuation efficiency. Between 9.0 and 11.0 m, the density tended to rise as wider channels increased the horizontal walking distance, resulting in longer evacuation times. Beyond 11.0 m, the density decreased again. (3) When the width of the first-level horizontal channel was 15.0 m or less, the maximum evacuation density and time decreased as width increased, thus effectively improving evacuation efficiency. Beyond 15.0 m, the maximum evacuation density continued to decrease, whereas the evacuation time remained stable. A width of 17.0 m was optimal for minimizing crowding and maintaining safety, though exceeding this brought unnecessary economic costs. (4) The suggested horizontal channel widths in deep underground spaces were 6.0 m for third-level channels, 9.0 m (left) and 13.0 m (right) for second-level channels, and 17.0 m for first-level channels. [Conclusion] By comparing the results of evacuation time and density across different widths, this study establishes reasonable horizontal channel widths at each level, providing both theoretical and technical support for the safe evacuation design of horizontal channels in deep underground spaces and helping establish a complete underground safety evacuation system. However, this study focuses on the preliminary design of safety evacuation in deep underground space. Future studies should incorporate the psychological and behavioral factors of pedestrians in deep underground spaces. Testing and calibrating theoretical results with practical engineering cases and actual evacuation data is crucial to improving the safety evacuation model for these deep underground spaces.
  • WU Xiangfei, NIU Jianlong, WANG Jingwu, ZHAO Xiaolong, TANG Fei
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 521-531. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.009
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Currently, the service life of buried gas pipelines in most Chinese cities has exceeded 20 a. Because of the combined effects of third-party damage, corrosion, and human error, the operational safety of these pipelines has decreased, leading to a high probability of accidents. Because gas pipelines are often buried underground, the transported medium is flammable and explosive, and they are located in densely populated areas, they possess characteristics of concealment and high risk. A leaking gas pipeline can easily trigger secondary fire and explosion incidents, causing severe casualties, property damage, and gas supply interruptions. [Methods] To prevent buried gas pipeline leakages and associated accidents in urban areas, a method for assessing the cascading risk evolution of pipeline leakages is proposed. Based on the operational characteristics of urban buried gas pipelines, risk factors, such as human factors, physical elements, environmental factors, and pipeline conditions are comprehensively integrated. These risk factors are further refined into a three-level indicator system. The logical relationships among the development process of leakage accidents and these risk factors are examined, and an integrated risk evolution and assessment model is developed through the Grey decision-making trial and evaluation laboratory (DEMATEL) method, interpretive structural modeling (ISM), and the dynamic Bayesian network (DBN) approach. The Grey-DEMATEL-ISM method is utilized to analyze causal hierarchical relationships among risk factors and identify key risk factors. The introduction of Grey numbers compensates for the subjective uncertainty of expert evaluations in the DEMATEL-ISM method, making the assessment results more aligned with actual conditions. Different types of risk factors are quantified through methods such as expert opinions, historical data, and probability distributions, and a pipeline leakage risk evolution DBN model is established based on the logical relationships among the risk factors. Furthermore, the probability of accident consequences is calculated and updated in real time to predict the potential development paths of accidents, thereby assisting in the safety operation and maintenance decision-making of buried gas pipelines. [Results] The results indicated that the corrosion and human factors exhibited high levels of causality, centrality, and sensitivity, making them critical risk factors for preventing pipeline leakages. A dynamic risk analysis model was conducted the use of the gas pipeline explosion accident in Shiyan, Hubei Province, as a case study. The results showed that the posterior probabilities of corrosion-related nodes were significantly higher than those of other nodes. Based on the assessment results, targeted measures were proposed to prevent initial incidents and cut off the accident propagation paths throughout the entire life cycle during the design and manufacturing phase, laying phase, and operational phase of the buried gas pipeline. The case study verified the feasibility and effectiveness of the proposed method. [Conclusions] The Grey-DEMATEL-ISM-DBN model can comprehensively consider the characteristics and quantitative representation of risk factors relate to humans, machines, environment, and management, even in situations with limited data. Based on historical accident information, it establishes a multi-level hierarchical structure model for the causes of buried gas pipeline leakage accidents. The model quantifies the degree of influence and sensitivity among various risk factors and provides an intuitive display of accident evolution scenarios. The proposed method enables the dynamic analysis of the risk evolution process of buried gas pipeline leakages, providing support for decision-making in safety operations, maintenance, and accident investigations involving urban buried gas pipelines.
  • WANG Chaozheng, GUO Dapeng, XU Bingzhou, LIU Chang, MIAO Zheng
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 532-546. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.021
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] To better understand how the structural characteristics of urban high-voltage cable tunnel networks affect fire propagation, this study aimed to provide a scientific basis for fire prevention, control strategies, and tunnel design optimization. Specifically, fire propagation behavior under different structural configurations and ventilation systems in cable tunnels was investigated. Key fire dynamics parameters analyzed included temperature distribution, smoke propagation, flame spread speed, and toxic gas concentrations for both single-layer and multilayer cable arrangements. [Methods] The investigation relied on three-dimensional fire dynamics simulations using Fire Dynamics Simulator software to explore fire propagation in urban high-voltage cable tunnels. Various tunnel structural configurations were analyzed, including single-layer and multilayer cable arrangements, with and without ventilation systems. Critical fire parameters, such as heat release rate (HRR), ceiling temperatures above the fire source, smoke flow patterns, and toxic gas concentrations, were examined under different fire scenarios. Numerical modeling provided detailed insights into how fire dynamics interact with tunnel structural features, emphasizing the significance of cable layout and ventilation. The study also explored how ventilation affects smoke behavior, assessing its influence on fire spread, temperature, and gas emissions. [Results] The results revealed distinct differences in fire behavior depending on tunnel structure and ventilation. Multilayer cable configurations caused ceiling temperatures above the fire source to rise significantly faster and reach higher peaks within 400 s compared to single-layer arrangements. This rapid increase in temperature indicated that the denser cable arrangement boosted the HRR, resulting in greater thermal effects. Smoke propagation was highly dynamic. During the early stages, it initially spread rapidly to one side of the fire source. Ventilation systems, however, altered this behavior by reversing the smoke flow direction over time. This reversal created localized temperature increases and more complex smoke distribution patterns. It also introduced cooler air into the system, influencing flame propagation and heat transfer dynamics. Toxic gas analysis showed that carbon monoxide levels were significantly greater for multilayer cables between 300 s and 500 s, indicating more incomplete combustion and increased hazardous gas emissions in these configurations. Flame propagation was faster, and heat transfer effects were more pronounced in multilayer configurations, highlighting the critical role of structural design in fire dynamics. These results underscore the heightened fire hazards posed by multilayer cable arrangements, including faster flame spread, greater thermal effects, and elevated toxic gas concentrations. [Conclusions] This study emphasizes the critical need for optimizing tunnel designs and ventilation systems to mitigate fire risk effectively in urban high-voltage cable tunnels. Multilayer cable arrangements notably increased heat release rates, toxic gas emissions, and flame propagation intensity, exacerbating fire hazards. While ventilation systems can positively influence smoke propagation and control localized temperatures, improper ventilation strategies may introduce risks, such as smoke flow reversal and uneven heat distribution. These findings provide essential technical insights for developing fire safety measures, highlighting the need for prevention and control strategies tailored to the structural and operational characteristics of urban high-voltage cable tunnels. By addressing these challenges, this research helps improve fire resilience, enhance infrastructure safety, and protect personnel during fire emergencies. This study provides a valuable scientific foundation for improving fire safety management in these tunnels, with practical implications for designing safer infrastructure and implementing effective fire prevention strategies.
  • HAN Hao, JIANG Xue, LIU Xudong, ZHANG Peihong
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 547-554. https://doi.org/10.16511/j.cnki.qhdxxb.2024.26.050
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] With the wide application of lithium-ion battery (LIB) in electronic devices, new energy vehicles, and energy storage power stations, the risk of fire and explosion accidents caused by the thermal runaway (TR) of LIB modules has increased, significantly hindering the development of the new energy industry of LIB. The influence of ventilation on TR and its propagation in LIB modules is complex and requires further investigation. [Methods] This study constructed a longitudinal ventilation environment and used an LIB module composed of five 5 Ah ternary lithium-ion batteries as the experimental object. By varying the wind speed, parameters such as the battery surface temperature, flue gas concentration, TR propagation time, and mass loss rate of LIB modules were compared and analyzed. This approach aims to clarify the dual action mechanism of the oxygen supply for combustion and heat dissipation from longitudinal ventilation on TR and its propagation in the LIB module. [Results] The experiments showed that under wind speeds of 2.0, 3.0, and 4.5 m/s, the combustion of the LIB module intensified during the TR process. The average temperature and temperature increased rate of the LIB module, the TR propagation time decreased, and the volume fraction of oxygen and carbon dioxide production increased. The radiant heat flux and mass loss rate of the LIB module were higher compared with nonwind conditions. At wind speeds of 6.0, 7.5, and 9.0 m/s, the TR phenomenon of the LIB module gradually weakened with increasing wind speed. At 7.5 and 9.0 m/s, the TR phenomenon did not appear in the initial phase of fire and combustion. The average temperature and temperature rise rate of the LIB module decreased, and the propagation time of TR was effectively extended. During TR, the oxygen volume fraction increased significantly, carbon dioxide production decreased, and the radiant heat flux and mass loss rate of the LIB module were correspondingly lower than those under no-wind conditions. [Conclusions] The results show that the influence of longitudinal ventilation on TR and TR propagation in the LIB module is determined by the synergistic effects of the oxygen supply, which promote combustion and heat dissipation. The dominant mechanism varies significantly with different wind speeds. At wind speeds of 2.0, 3.0, and 4.5 m/s, the oxygen supply promoting the combustion effect is dominant, enhancing TR and its propagation in the LIB module, with the strongest effect occurring at 2.0 m/s, followed by 3.0 m/s, and the weakest effect occurring at 4.5 m/s. At wind speeds of 6.0, 7.5, and 9.0 m/s, the heat dissipation effect is dominant, inhibiting TR and its propagation, with inhibitory effects increasing along with wind speed. The research results provide theoretical and technical support for the thermal safety of LIB and their application in the new energy sector.
  • ZHANG Jiangshi, LI Yongtun, WU Jingru, REN Xiaofeng, PAN Yu, ZHANG Qi
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 555-568. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.006
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] An effective causal analysis of accidents is essential for learning from and preventing coal mine accidents. Manual analysis of accidents is strongly influenced by the subjectivity of the personnel involved and becomes inefficient for analyses involving large volumes of accident and risk text data. Although considerable research has been conducted in the area of accident text mining, most studies directly apply data mining techniques to extract accident information and factors from texts without considering accident causation theories. This approach leads to results that lack systematic and logical coherence. [Methods] To address the aforementioned issues, this paper proposes a method for the intelligent identification of accident causes in the coal mining sector. This method integrates entity recognition, semantic dependency analysis, text classification, and the accident causation “2-4” model (24Model). Specific implementation steps for this method are also provided. Accident causation theory is crucial for ensuring the effectiveness and scientific validity of accident analysis. This paper introduces the 24Model as a theoretical basis for accident cause identification, and the advantages of the model in the intelligent analysis of accident causes are highlighted. Entity recognition technology is employed to identify key entity information in accident texts, including information on personnel, organizational structures, accidents, abnormal characteristics, values, safety management, building facilities, environments, equipment and materials, safety policies, procedural documents, and operational processes. To effectively identify this information, this paper integrates the bidirectional encoder representations from transformers(BERT)-bidirectional long short-term memory(BiLSTM)-conditional random fields(CRF) model and trains the combined model using 660 accident texts. This paper utilizes semantic dependency analysis technology to identify the semantic relationships among entity information. Text representation patterns were extracted according to the definitions of unsafe individual actions and organizational factors by the 24Model, and these definitions were used to determine the types of accident causes. This paper utilizes a text classification method to develop a model for identifying individual capability-related causes, and the focus is on five aspects: knowledge, awareness, habits, and psychological and physiological factors. The text classification model was based on BERT. This paper evaluates the accuracy of the proposed method in identifying accident causes by comparing both the entity recognition model and the text classification model with similar models. Test cases were selected, and the results of accident cause analysis via the proposed method were compared with those from manual analysis. Additionally, this paper develops an application based on the proposed method to facilitate the analysis and learning of onsite accident cases by employees of coal mining enterprises. [Results] This research results showed that the precision rates of the trained entity recognition model and the text classification model reached 95.42% and 96.11%, respectively. Additionally, the accuracy of the accident cause identification method, when combined with semantic dependency analysis, reached 73.09%. [Conclusions] The contributions of this paper are as follows: (1) Integration of the definition and concept of the 24Model and the automatic identification of unsafe behaviors according to the 24Model. This approach helps avoid the strong subjectivity and inconsistency often present in accident analysis conducted by different personnel. (2) Further identification of the actors of actions, operational procedures, materials, and equipment. (3) Fusion of multimodel algorithms to identify the causes of accidents, allowing for the rapid analysis of a large number of accidents. (4) Facilitating the application of accident causation theory in coal mining enterprises, enhancing the effectiveness of accident case analysis and learning, thereby achieving the objective of preventing related accidents.
  • LI Gang, XU Xiuping, LIU Jianguo, JIN Longzhe
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 569-579. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.008
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Blasting operations in mines generate significant amounts of impact dust, which poses a significant risk of pneumoconiosis and threatens workers' health. This also results in substantial direct economic losses annually, severely impacting China's pursuit of high-quality economic development. Existing ventilation and dust removal technologies have proven inadequate. Among various developed dust reduction methods, wet dust suppression technology has improved continuously. A key advancement is the use of ultrasonic atomizing nozzles, known for their low water consumption, effective atomization, and superior dust capture efficiency. Consequently, a dry mist dust suppression technology has been proposed to efficiently manage dust from mining blasting operations, improve working conditions in return airways, and safeguard the physical and mental health of workers. This study investigates the application of dry mist dust suppression technology in roadway-type mining faces of metal mines to achieve these objectives. [Methods] This study focused on a roadway-type mining face in a certain iron ore mine. The initial investigation involved analyzing the physical and chemical characteristics of dust generated during blasting operations, particularly examining mechanisms that influence its wetting properties. Dust samples were classified through flotation into hydrophilic and hydrophobic types. Further analysis was conducted on their wetting properties, surface morphology, particle size distribution, and surface pore structure to explore the physicochemical characteristics affecting dust wetting. Three types of ultrasonic atomizing nozzles were selected for testing their atomization characteristics under different pressure parameters using an atomization test platform. This study analyzes the influence of different air-water parameters on atomization characteristics and identified the optimal nozzle for dust reduction applications. Furthermore, a dry mist dust suppression device was designed and developed for use in mining. Field experiments evaluated dust distribution in return airways before and after blasting operations, with and without the application of mist spraying for dust suppression. [Results] This research indicated that the impact dust generated during blasting operations was predominantly hydrophilic. The dust wetting properties were primarily influenced by factors such as particle size and surface porosity. Critical to the efficiency of dry mist dust suppression were the droplet size and quantity. The median droplet size D50 showed an inverse relationship with the ratio of air pressure to water flow rate. Among the tested nozzles, SK-508, SV-980, and SV-882, the SK-508 ultrasonic atomizer exhibited the smallest average droplet size and consumed the least amount of water, thus conserving water resources effectively. Under conditions of 0.7 MPa air pressure and a water flow rate of 0.1 kg/s, the SK-508 demonstrated significant atomization effects, making it the optimal nozzle for dust suppression among those tested. Leveraging the advantages of dry mist dust suppression technology and the atomization characteristics of ultrasonic nozzles, a dry mist dust suppression device was developed. Field tests of the prototype demonstrated its notable effectiveness in reducing both total dust and respirable dust, achieving a high level of dust suppression efficiency. [Conclusions] A dry mist dust suppression device was developed to address the issue of dust in return airways, effectively managing both total dust and respirable dust in mining tunnels. It achieves efficient control of impact dust generated during blasting operations. This innovation provides a solid theoretical foundation for advancing the national green mining initiative and contributes to establishing a comprehensive technical system for dust control in mines.
  • TAO Zhenxiang, ZHOU Haozhi, LIU Xiaohan, HU Peifeng, LUO Ning, ZHOU Biao, YANG Rui, LIU Quanyi
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 580-588. https://doi.org/10.16511/j.cnki.qhdxxb.2024.26.047
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Unlike fire accidents caused by direct combustion, which are easy to detect and control, numerous accident investigations have shown that ignition caused by high-temperature molten metal droplets is a key link in the spread and escalation of fire events. Recently, a number of studies in China and other countries have investigated the ignition mechanism of fire spreading owing to high-temperature molten metal droplets. However, most of these studies have investigated the ignition mechanism alone, and only a few studies have explored the dynamic characteristics of molten metal droplets impacting the wall surface. Consequently, in this experimental study, the dynamic of molten aluminum particles with different particle diameters impacting noncombustible stainless-steel plates and combustible expanded polystyrene (EPS) foam were explored. [Methods] To address the problems of adhesion of molten metal particles to the heating vessel and the low upper limit of the heating temperature in traditional research, this study adopted a non-contact method to heat and control the metal particles and recorded the spreading and retraction of molten metal droplets and the pyrolysis and combustion phenomena of the substrate material in the collision process using a high-speed camera. The morphology of the molten metal can be used to visualize droplet motion and assess the heating conditions of the impacted substrate and surrounding objects. The droplet spreading diameter is one of the most important parameters for studying the dynamic characteristics of molten metal droplets impacting the wall surface; however, after impact, molten metal droplets do not have standard circular surfaces, and their contours cannot be directly obtained. Using Image-Pro Plus software to draw the outer contour of the thin layer of metallic aluminum particles, the number of pixels in the coil is calculated and multiplied by the scale; furthermore, the actual spreading area of the thin layer of molten metal S0 is calculated, and the equivalent spreading length D is subsequently obtained. In the experiment on molten droplets impacting on a stainless-steel metal surface, the abovementioned factors were analyzed by varying the particle diameter of the molten aluminum droplets, the distance of impingement, and the particle temperature to analyze the effects of the above factors on the spreading pattern, spreading distance, and motion changes of the droplets. In addition to considering the effects of pyrolysis of the stainless-steel wall, the effects of pyrolysis of EPS foam on the droplet spreading process should also be analyzed in an experiment on the impact of molten metal droplets on the EPS foam surface. [Results] The results of the experiments of molten droplets impacting stainless-steel metal surfaces and EPS foam surfaces showed that: (1) Aluminum droplets oxidize on the surface of the substrate particles during the impact process; this effect was more pronounced for particles with larger diameters, and larger particles took a longer time to reach the maximum spreading diameter. (2) The spreading characteristics of particles with different initial temperatures were roughly the same, and the maximum spreading length achieved by the particles was approximately 25 mm. (3) The collision of aluminum particles with the EPS foam reached a maximum spreading at t=20.0 ms during the aluminuny particles motion, followed by a decrease in the spreading length of the thin layer of metal particles. [Conclusions] Larger initial diameters and impact distances imply that the peak of the spreading diameter is also larger. The initial temperature has little effect on the spreading motion of molten droplets in the study range. When the base plate material is EPS foam, the spreading diameter generally increases and subsequently decreases with time, and there is no stabilization phase.
  • WANG Hetang, ZHANG Qi, WANG Yuxuan, LI Xiaojuan, YANG Panpan, XU Yifei
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 589-600. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.020
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Dust is easily generated in all parts of the mining process, which can cause dust explosions and lead to occupational pneumoconiosis in workers. Therefore, dust reduction is crucial in mining operations and plays an important role in protecting the environment and worker health. While traditional chemical dust suppressants provide short-term effectiveness, they pose considerable challenges. These include poor resistance to natural degradation, low environmental performance, and the risk of secondary pollution of soil and water sources, creating a growing demand for sustainable alternatives. To address these challenges, the authors proposed the idea of using microbial fermentation to synthesize biological dust suppressants. [Methods] To enhance the scalability and efficiency of microbial fermentation dust suppressant (MFDS) production, this study used a self-developed experimental device for fermentation and synthesis of biological dust suppressants. Six straight-blade disc paddles (6S-DR), six semicircular-blade disc paddles (6S-SDR) and four inclined straight-blade paddles (4-IR) were selected to design eight mixing combinations. Numerical simulations were performed to analyze key parameters, such as the matrix flow field velocity, turbulent kinetic energy, gas holdup, and stirring power. In addition, experimental tests were conducted to validate gas holdup and MFDS yield under different conditions. Furthermore, MFDS solutions of three purities, acid precipitation, single-stage ultrafiltration, and two-stage ultrafiltration, were characterized using liquid chromatography-mass spectrometry (LC-MS). MFDS was conducted using an LC-MS system. Tests for interfacial performance and wettability of MFDS were also performed to assess its dust suppression capabilities at different purity levels. [Results] The results indicated that the stirring combination, using four oblique straight paddles, formed an obvious liquid circulation area. Adding an extra layer of paddles further expanded this area, enhancing heat and mass transfer. Stirring combination H demonstrated smaller liquid surface fluctuations and a more uniform distribution range of turbulent kinetic energy. These features increased dissolved oxygen levels and improved gas-liquid mass transfer, fostering microorganism growth. Consequently, stirring combination H achieved better substrate homogenization and delivered better fermentation performance under identical conditions.Interfacial performance tests revealed that when two-stage ultrafiltration was used, the critical micellar concentration of the bioreactor solution decreased to 22.65 mg/L. The higher the purity of the MFDS solution was, the smaller the critical micelle concentration of MFDS, and the viscoelastic modulus stabilized with increasing concentration. The ultrafiltration technique had a greater effect on the MFDS viscoelastic modulus, but the number of ultrafiltration cycles had a smaller effect on the viscoelastic modulus. The viscoelastic modulus stabilized as the concentration increased. Moreover, wettability tests indicated that the two-stage ultrafiltration purification improved dust suppression efficiency, achieving the shortest dust settling time of 67 s and remarkable wettability performance. Considering the results of the surface interface performance and wetting performance tests, the dust suppression performance of the two-stage ultrafiltration MFDS was superior. [Conclusions] This study investigated large-scale synthesis methods for biodust suppressants, focusing on matrix homogenization to increase production efficiency and scalability. This approach addresses key shortcomings of traditional chemical dust suppressants, including poor degradability, low surface activity, and environmental harm. By overcoming these issues, this study offers a meaningful solution to reduce dust pollution, protect the environment, and improve the occupational health of miners. Meanwhile, this study provides valuable theoretical guidance for the efficient development of green and environmentally friendly biological dust suppressants.
  • QIN Lei, WANG Hui, LI Shugang, LIU Pengfei, LI Jiawei
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 601-613. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.005
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Coalbed methane (CBM) outbursts pose a threat to coal mining safety, particularly in China, where many coal seams exhibit low permeability, making CBM extraction inefficient. Liquid nitrogen fracturing technology is an important anhydrous method for enhancing coal seam permeability. The phase transition of coal pore water under low temperature conditions profoundly affects the coal pore structure. Studying the pore-water phase transition and pore changes during low-temperature liquid nitrogen freeze-thaw cycles is crucial for evaluating the effectiveness of this technology and improving CBM extraction efficiency. [Methods] This experiment involves four groups of water-saturated bituminous coal samples, each subjected to different numbers of liquid nitrogen freeze-thaw cycles. The transverse relaxation time T2 curve, cumulative porosity, and cumulative pore throat distribution of coal samples in the initial, frozen and thawed states were evaluated by using low-field nuclear magnetic resonance (NMR). [Results] These findings showed the following: (1) Even at -196 ℃, trace amounts of unfrozen water remained in the coal pores. As the temperature increased, the ice in the smaller pores first melted, producing micropore unfrozen water. Smaller pore diameters corresponded to lower pore melting pints of pore ice. (2) As the pore ice gradually melted, the cumulative porosity increased in three stages: exponential function type rapid growth, primary function type uniform growth, and quadratic function type growth. (3) By integrating the T2 curve, it is observed that the unfrozen water content changed exponentially with increasing temperature in the negative temperature intervals. Closer to 0 ℃, the unfrozen water content growth rate increased. At temperatures between -196 ℃ and -34 ℃, micropore unfrozen water constituted nearly 100%, whereas temperatures between -20 ℃ and 0 ℃ saw mesopore and macropore unfrozen water dominating. (4) NMR tests on raw and fully melted coal samples revealed that the permeability of coal samples increased with the number of freeze-thaw cycles. The effective porosity and permeability growth rates increased from 20.68% to 24.15% and 109.73% to 122.20% respectively, as the number of cycles increased from 5 to 30. However, the marginal utility of permeability enhancement of coal samples became increasingly pronounced. [Conclusions] Liquid nitrogen cyclic freezing of coal effectively increases the effective porosity and permeability of water-saturated coal. Liquid nitrogen fracturing can be used as an anhydrous fracturing technology to increase the coal seam permeability in water-scarce coal mining areas, thus improving the CBM extraction efficiency. The results of this study can provide valuable insights into the evolution of seepage pore structure during low-temperature medium cyclic freezing and thawing of coal and guide field operations related to liquid nitrogen fracturing and permeability enhancement.
  • WANG Xinyu, WANG Enyuan, YUE Jianhua, ZHU Guoqing, CHENG Deqiang, LIU Xiaofei, LI Dexing, ZHENG Sijian
    Journal of Tsinghua University(Science and Technology). 2025, 65(3): 614-624. https://doi.org/10.16511/j.cnki.qhdxxb.2025.26.010
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    [Objective] Coal and gas outbursts are significant disasters impeding the safe and efficient production of coal mines. The stress and structural anomaly zone ahead of the coal heading face is a high-risk incidence area for these outburst. Timely and precise detection of spatial characteristics in front of the heading face for stress anomalies is crucial for disaster prediction. However, relevant research and methodologies are currently inadequate. [Methods] This study explores the feasibility of using the direct current (DC) method for the advanced detection of stress distributions and anomalous zones in coal seams. Experiments were conducted on coal samples under uniaxial graded loading to test apparent resistivity and the temporal and spatial evolution laws of apparent resistivity during this process were investigated. Subsequently, continuous tracking and advanced detection were performed on the heading face of a mine using the DC method. The regional correspondence between the stress distribution and apparent resistivity was analyzed in the coal rock body in front of the heading face. This analysis was compared with drilling chip indexes and actual excavation results. Results revealed the characteristics of the geos-spatial response in the stress anomaly. [Results] The temporal and spatial evolution of the apparent resistivity in loaded coal samples was closely related to stress changes. High-resistance regions and mean apparent resistivity initially decreased and then increased with increasing stress, which was attributed primarily to structural evolution in the coal samples. The apparent resistivity of the heading face coincided with the stress distribution of the coal seam, showing an initial decrease followed by a gradual increase and stabilization from the headland to the inner part of the coal seam. This pattern refleced the distribution of pressure relief zones, stress concentration zones, and original stress zones in the coal body. When irregular high and low resistances were observed in the original stress region, they indicated stress or structural anomalies. [Conclusions] The spatiotemporal evolution law of apparent resistivity in loaded coal samples accurately reflects the stress levels and damage evolution. The apparent resistivity response characteristics at the coal seam front aligned well with the spatial stress distribution. Irregular high and low resistance anomalies in apparent resistivity, deviating from the typical “stress three zones” distribution, are key indicators of stress or structural anomalies in the coal rock body. The detection results are verified using the drilling chip index and on-site excavation. The findings indicate a good correspondence between the detection results and verification indices. Notably, all classified stress anomaly areas are identified as having a high risk of coal and gas outbursts. When detecte the stress distribution state of coal body by DC method, the ideal detection accuracy can be obtained by narrowing the detection range and increasing the density of electrode arrangement within a reasonable range. This suggests that the DC method can effectively detect and identify stress distributions and anomaly areas ahead of the heading face. The research results provide a solid theoretical basis and technical support for using the DC method to assess coal and gas outburst risks. Furthermore, this approach introduces possibilities for the advanced detection and early warning of coal rock dynamic disasters.