Fire Protection and Evacuation Analysis in Underground Interchange Tunnels by Integrating BIM and Numerical Simulation

: Rescue and evacuation of underground interchange tunnels after a ﬁre are challenging. Therefore, a method of integrating building information modeling (BIM) and a ﬁre dynamic simulator (FDS) was proposed to analyze ﬁre characteristics and personnel escapes in underground interchange tunnels. A BIM model of underground interchange tunnels was built, and then diﬀerent formats (DXF and CAD) were generated and imported into Pyrosim software and Pathﬁnder software. With an increase in ventilation velocity, the CO concentration and temperature downstream of the ﬁre source increased, and visibility decreased, according to simulation results. The critical ventilation velocity was 3.6 m/s at 30 MW. Evacuation simulation results suggested that the congestion of the transverse passage was very unfavorable for personnel escape: the escape time increased by 14.9% and 20% when the interior and entrance of the transverse passage were severely congested, while a 2.5 m wide transverse passage eﬀectively reduced the escape time. Visibility was the ﬁrst indicator that it did not meet the safety of the escape. After the tunnel’s personnel have been evacuated, the air supply or exhaust system should be started, and smoke should be expelled at a higher velocity. It is necessary to clear the passageway quickly or increase the automatic ﬁreﬁghting facilities when congestion is severe.


Introduction
The development and utilization of underground space interchanges are continuously improving, and sudden, uncertain, and uncontrollable safety risks are also increasing, which can easily lead to significant safety accidents [1,2].A large amount of flue gas released by combustion will rapidly reduce visibility, and the temperature will continue to rise.Meanwhile, the high temperature and toxic smoke generated by space restriction combustion cannot be discharged in time due to the long and narrow space of the underground interchange, which increases the risk of fire [3,4].Considering the particularities of the underground interchange and the severity of the fire, it is critical to analyze ventilation characteristics and fire rescue plans.
Researchers from industry and academia have performed much work studying fire smoke and have made some achievements in recent years [5][6][7].Currently, researchers focus on studying underground buildings in three aspects: physical experiments [8,9], model tests [10], and computer simulations [11,12].In terms of physical experiments, researchers have obtained actual physical quantitative data of the firing process, including the law of smoke diffusion, CO concentration change, and temperature field change, to Citation: Liu, Z.; Gu, X.; Hong, R. provide support for theoretical research [13].A smoke movement model was developed in underground space, and smoke model data were obtained by Matsushita [14].Then, a 1:20 test bed was built based on the underground space model to conduct an experimental simulation of the smoke.Experimental and model data were obtained to analyze the smoke situation in underground spaces.After that, many scholars carried out fire test analyses on rooms [15], roofs [16], shopping malls [17], and so on.However, large-scale physical experiments have the disadvantages of high cost, difficulty, and danger, and the participants are not afraid.Therefore, less effort is being made to apply this method.
Compared with physical experiments, model tests based on similarity theory and dimensional analysis are critical in studying fire development.Smoke flow characteristics in long and narrow underground spaces have been widely used [18,19].Li et al. used a reduced-scale model (1:6) to conduct fire experiments in tunnels, and there was a good agreement between the measured and simulated fire parameters (flame tilt angle and mean flame temperature downwind) [20].A study by Giache i et al. examined the characteristics of smoke in subway stations and investigated the relationship between the station's safe height and ventilation velocity.It is found that the lower the safety height of the station, the greater the ventilation volume [21].Based on full-scale experimental data and computational fluid dynamics theory, Wang et al. examined the fire characteristics of naturally ventilated tunnels [22].Computer fluid dynamics (CFD) software achieved greater accuracy than measured data in predicting flue gas temperatures, concentrations, and se lements.
Over the past few years, numerical simulations have become a more effective method of fire research due to the development of computer science technology.Numerical simulation studies [23][24][25] on narrow and long underground interchanges mainly focus on the plume characteristics near the fire source (maximum roof temperature, impact location, etc.), far-field temperature a enuation, smoke spread rate, and other issues [26,27].In particular, building information modeling (BIM) and artificial intelligence algorithms [28][29][30] have become more effective methods in the construction industry due to their three-dimensional (3D) capabilities [31][32][33][34].The combination of BIM technology and fire analysis has promoted fire simulation analysis to a new height.Chen et al. proposed a BIM fire rescue visual early warning system [35] to simulate various structural fire scenarios with the fire dynamic simulator (FDS) and BIM software.BIM models can display simulation results of temperature, carbon monoxide, and visibility for brief presentations before rescue operations begin.Li and Burcin have developed an emergency rescue system based on environmental awareness using the BIM model to improve the positioning accuracy at the room-level [36].Based on a BIM model and a cell-discrete event system (DEVS) evacuation algorithm, Wang et al. conducted cloud simulations and studied crowd evacuation behavior [37].Computer simulation has great application prospects.However, there are few simulation studies on the ventilation characteristics and fire transmission law of urban underground interchange, and even fewer studies combining BIM technology.In addition, few simulations consider both fires and escape.
Therefore, aiming at the research needs of the ventilation characteristics and fire rescue of underground interchanges, an integrated method of BIM and numerical simulation was developed to perform detailed fire characteristics and evacuation analysis of underground interchange tunnels in this paper.The extraction and rescue strategies were given according to the simulation results, which could provide valuable references and suggestions for the fire control design and safety management of similar underground interchange tunnels.This manuscript is arranged as follows.Section 2 describes the methodology of BIM modeling and fire numerical simulations in underground interchanges.Sections 3 and 4 discuss the simulation results of fire characteristics and evacuation, respectively.Finally, Section 5 concludes this study.

Methodology
The entire study design is illustrated in Figure 1, including the modeling of underground interchange, fire characteristics, and evacuation simulations.Details are described in the following sections.

BIM Models of Underground Interchange
BIM models can convey important information to fire simulation software.Therefore, an underground interchange model was built in the MicroStation V8i software of Bentley Company.This model adopted a typical longitudinal exhaust underground interchange, and the ventilation air flowed along the underground interchange.It was a full-jet underground interchange, and its ventilation was completed by the boost of the jet fan group.The cross section was uniformly set as a rectangle of 10 m × 6 m, and the length was 434 m based on the actual situation.Then, the BIM file was directly imported into the fire simulation software that supports the file format, which can save the development process of the software itself and is easy to operate.

Fire Characteristics Simulation
Pyrosim software is a simulation software used for the dynamic simulation of fire.As a preprocessing software of FDS, Pyrosim software can directly convert the information programmed by FDS into a visual simulation model, which can accurately predict the flow of fire smoke, the record of fire temperature field, the concentration of harmful gases, and their dynamic changes.Pyrosim can directly import the files of the BIM model without making any modifications.The mesh is the smallest calculation unit in the Pyrosim simulation system, and the mesh partition size is inversely proportional to the simulation accuracy.To verify mesh independence, the dimensionless formula D*/δx is usually used to verify the rationality of mesh division [35], and the solution expression of D* is shown in Equation (1).
where Q represents the heat release rate of fire source, ρ represents the air density, C represents the specific heat capacity of air, T represents the ambient temperature, and g represents the gravity acceleration.The D*/δx control range is between 5 and 15; thus, the grid size was fixed as 0.4 m × 0.4 m × 0.4 m according to the external conditions of this study.
With the help of the BIM 3D model and grid generation tools, the model presents the relevant parameters of fire simulation.It uses fire dynamics theory to simulate complex fire scenes, including fire direction, smoke spread, and temperature conditions.
After importing the Pyrosim model from the BIM model established above, the grid, vent, fire source, etc., were created, as shown in Figure 2.Then, the material and reaction were defined to obtain a complete numerical simulation model of fire.The t 2 fire model was used to describe the fire trend, and it was assumed that the fire development rate was very fast.After 400 s of fire exposure, the HRR reached a stable value of 30 MW, and the HRR remained unchanged until the end of the simulation calculation at 900 s.When the fire burning position is at the middle point of the simulation interval, the personnel can only choose the nearest cross passage or an underground vertical intersection for evacuation according to their positions.As a result, the evacuation distance of the personnel is the longest, making it the most unfavorable situation.
Then, a total of 11 measuring points was set at 202 m, 197 m, 150 m, 126 m (near the entrance of the left transverse passage), 100 m, 50 m, 0 m, and 50 m, 100 m, 126 m (near the entrance of the right transverse passage), 150 m, and 200 m downstream of the fire source.The CO concentration, visibility [38], and temperature during the fire simulation were measured at 1.8 m above ground level.According to the Fire Safety Manual of China (GB/T 31593.9-2015)[39], the CO concentration cannot exceed 0.0012 mol/L (1200 ppm), visibility cannot be lower than 10 m, and temperature cannot exceed 60 °C [40].

Determination and Verification of Critical Ventilation Velocity
Critical ventilation velocity correlates with underground interchange height, width, and fire heat release rate.It is the longitudinal ventilation velocity at which smoke does not regress during a fire that determines the critical ventilation velocity.A backflow layer occurs when longitudinal ventilation velocity is less than the critical ventilation velocity, causing flue gas to develop against the direction of ventilation.It is unfavorable for personnel to carry out rescue operations from upstream of the ventilation.Backflow is not caused by flue gases when longitudinal ventilation velocity exceeds critical ventilation velocity, but stratification is.The higher the ventilation velocity, the more noticeable this disorder phenomenon will be, which is not conducive to evacuating trapped people toward downstream ventilation.Therefore, determining the critical ventilation velocity of underground interchange is an essential part of fire rescue research.Figure 3  According to the design guidelines for tunnels of China (JTGT D702-2014), the longitudinal ventilation velocity was set as 0.5 m/s, 1.0 m/s, 2.0 m/s, 3.0 m/s, and 4.0 m/s.As shown in Figure 4, when the ventilation rate was 0 m/s, 0.5 m/s, 1.0 m/s, 2.0 m/s, and 3.0 m/s, the flue gas diffused asymmetrically to both sides of the underground interchange, and the stratification phenomenon occurred.It appears there was no stratification of smoke when ventilation velocity was 4.0 m/s.Then, the longitudinal ventilation velocity was further set as 3.1 m/s to 3.7 m/s.In the simulation, the flue gas did not experience flue gas regression when the ventilation velocity was 3.6 m/s.Consequently, 3.6 m/s was the critical ventilation velocity that met industry standards.From the front view, the flue gas did not produce regression at 3.6 m/s.Smoke can be seen entering the transverse passage downstream of the fire source from the top view.There was still smoke entering the transverse passage when the ventilation velocity was increased to 3 m/s.Therefore, it was recommended to close the transverse passage to prevent smoke diffusion through the transverse passage when 3.6 m/s was used as the ventilation velocity.

Evacuation Simulation
Pathfinder software is a simulator based on human evacuation and movement simulation, and is a commonly used software in tunnel and underground space human evacuation simulation.As described in Figure 1, the BIM model was converted into a CAD base map and imported into Pathfinder software, and then the elevation, evacuation exit, escape passage, and other components of the tunnel model were drawn.The steering mode was selected for simulation, and the calculation mechanism took real factors into consideration.The path planning, guidance mechanism, and collision processing were combined to control the movement of personnel so that the individual can act freely and avoid collision automatically, which can truly reflect the speed and evacuation time of the flow of people to be evacuated in the complex passage.
Evacuation standard: Trest refers to the time necessary for all personnel to be evacuated safely to another area after a fire; in fire situations, Taset refers to the interval between the time a fire occurs and when the entire environment cannot support safe evacuation.Then, the criterion of successful evacuation is Trest < Taset.Moreover, Trest can be obtained through simulations by Pathfinder software.
Evacuation scene: the movement mode was set to steer mode, the alarm time of the fire smoke alarm was set to 50 s, the walking speed was 1.2 m/s, the vehicle congestion density was 150 vehicles/km to simulate moderate congestion, and 200 vehicles/km to simulate heavy congestion in the underground interchange.The personnel evacuation conditions were set as follows:


Working condition #1: the vehicle congestion density was 150 vehicles/km, and the total number of evacuees was 410.


Working condition #2: the vehicle congestion density was 200 vehicles/km, and the total number of evacuees was 540.
In both conditions, the evacuation terminus was two cross passageways and an underground interchange exit.
A simulation of evacuation was performed after se ing the relevant parameters.Frequent doors were set in key locations to facilitate the recording of critical control point data, as shown in Figure 7.The main line and the cross passage were set as separate rooms.Then, the main line of the underground interchange was recorded as the main hole, and the passageways on the left and right were recorded as PASS1 and PASS2, respectively.Door-PASS1 and Door-PASS2 were set at the entrances of the two passageways, and EXIT1 and EXIT2 were set at the ends of the two passageways.Moreover, the personnel was randomly distributed in the underground part of the underground interchange according to different working conditions.

CO Concentration
Figure 8 shows the changes in the longitudinal center section of CO concentration at different distributions.Furthermore, Figure 9 describes the 3D curves of CO concentrations at different distances.The small longitudinal and critical ventilation velocities were 0.5 m/s and 3.6 m/s, respectively.Under a 0.5 m/s longitudinal ventilation velocity, as the distance from the fire source increased, the CO concentration also changed.A longer underground interchange may have contributed to the higher concentrations upstream of the fire source than downstream.There was a lower CO concentration 50 m away from the fire and a higher concentration 50 m away.At 3.6 m/s longitudinal ventilation velocity, conditions were similar to 0.5 m/s.The distance of 50 m was also the dividing point of CO concentration, and the peak of CO concentration at both velocities was obtained directly above the ignition source.V changed with distance from the fire source when the longitudinal ventilation velocity was 0.5 m/s.After 200 s, the visibility in the whole tunnel was below the critical value.A fire with 3.6 m/s velocity had very different visibility distributions on the left and right sides.After 200 s, visibility was almost entirely below 10 m downstream and above 10 m upstream.Within 300 s after the fire, it was almost impossible to see downstream.From the figures, T variations were centered around the fire source at 0.5 m/s longitudinal ventilation velocity, and upstream and downstream T differed slightly.Inside 100 m of the fire source, T was higher than the critical value, even though the T outside 100 m was lower.A velocity of 3.6 m/s reduced the temperature upstream of the fire to less than 60 °C, while that downstream of the fire was higher than 60 °C after 300 s.In addition, the temperature difference between upstream and downstream was much larger than that at low ventilation velocity.

Ventilation Velocity of the Cross Passage
To measure CO concentrations, visibility, and temperature near the entrance of the transverse passage in response to ventilation velocity in the transverse passage, the working conditions were set as follows: In both conditions, the ventilation velocity of the mainline was 0.5 m/s.Moreover, the change curves of these conditions are shown in Figures 14 and 15  According to the simulation results, working conditions A and B were basically the same in CO concentration, visibility, and temperature.The possible reason was that the cross section and ventilation velocity of the cross passage were low, which had li le influence on the fire development of the main line.Reviewing the available literature [41,42], a transverse ventilation velocity of 0.5 m/s in the escape stage is thought to prevent smoke from flowing into the transverse passage and affect its ability to escape.In the meantime, it will not cause a smoke disorder because it has li le effect on smoke diffusion in the main line.

Analysis of Fire Characteristics under Different Fire Scales
For the convenience of comparison, the diagram of CO concentration, visibility, and temperature change is described in Figure 16 when the fire scale was 30 MW, 20 MW, and 5 MW, and the ventilation velocity of the cross channel was 0 m/s.Evidently, when the fire scale was 20 MW and 5 MW, compared to the fire scale of 30 MW, there were no significant differences in CO concentration, visibility, or temperature with distance from the fire source.Therefore, a critical ventilation velocity of 3.6 m/s can be used as the exhaust ventilation velocity in the rescue stage when a minor fire occurs in the local underground interchange.

Working Condition #1
Working condition #1 was simulated, as shown in Figure 17.People started to enter PASS1 and PASS2 at 21 s and 17 s, respectively.They entered all the above positions at 218 s and 211 s, respectively.Namely, no people were on the main line, and everyone was relatively safe.The personnel began to arrive at the escape exit at 60 s and 56 s, respectively, and all left the exit at 277 s and 270 s, respectively.That is, at 277 s, all the personnel successfully escaped.Working condition #2 was simulated, as shown in Figure 18; people started to enter PASS1 and PASS2 at 15 s and 16 s, respectively, and all entered the above two positions at 229 s and 216 s, respectively.Namely, no people were on the main line, and everyone was relatively safe at 229 s.The personnel, respectively, at 50 s and 57 s, began to arrive at the escape exit and left at 288 s and 275 s.Namely, all the personnel successfully escaped at 288 s.According to the simulation results of moderate congestion (walking speed was set as 0.7 m/s) described in Figure 19, the time it takes for everyone to enter the safe area remains the same, while it takes 16.6% longer for everyone to escape in the simulated moderate congestion situation.For heavy congestion simulations, these two times increased by 2.2% and 26.4%, respectively.Therefore, it is necessary to check and clean up the debris inside the cross channel regularly to ensure the normal use of the cross channel.(2) The patency of the transverse channel entrance A 3 m by 0.5 m obstacle was set up 0.5 m in front of the PASS1 and PASS2 entrances to simulate obstructions blocking the transverse passage entrance.According to the simulation results that are presented in Figure 20, compared with the situation without debris, the time for all people to enter the safe area increased by 19.3%, and the time for all people to escape increased by 19.6% under moderate congestion.Furthermore, in heavy congestion, these two times increased by 45.9% and 34.0%, respectively.Therefore, it is necessary to regularly check and clean up the obstacles outside the entrance of the transverse channel to ensure the normal use of the transverse channel.

Effect of the Fire Escape Width on the Evacuation Simulation
The transverse passage width was 2.5 m (the initial se ing was 2.0 m) to analyze the effect of the width of the transverse passage on personnel escape.According to the simulation results, the escape time of this situation was the same as the condition of the 2.0 m transverse passage, which illustrates that the width of the cross passage is not the main factor affecting the escape simulation of the underground interchange.

Safety Analysis of Personnel Escape
For personnel to escape, the fire heat release rate and the ventilation velocity were set at 30 MW and 0.5 m/s, respectively (a small ventilation velocity should be selected when personnel escape).Thus, fire simulation and personnel escape simulation results in each working condition are shown in Table 1.First, the congestion degree of the fire escape had a great impact on the start time and duration of the escape.The more severe the crowding, the later the escape began, and the longer it lasted.This was especially true when the transverse passage was externally congested.However, the above two times were basically consistent with working condition #1 when a 2.5 m wide transverse passage was set, indicating that this measure is very effective.Depending on the results in Table 1, working condition #1 and working condition #7 met the set requirements of fire simulation.In the evacuation process of other working conditions, the visibility reduction caused by smoke diffusion is the main factor affecting evacuation.Strayed personnel may be at risk in other working conditions due to poor visibility.Taking working condition #6 as an example, the 2D slice and 3D cloud views of the simulation process are displayed in Figure 21.It can be intuitively seen that the value of CO concentration and temperature rose slowly, which had li le influence on the escape of personnel in each working condition, while with the passage of time, the visibility was undeniable and was below the critical value after 200 s, which was the main factor affecting people's escape.

Extraction and Rescue Strategies
The influencing factors of underground interchange crowd evacuation mainly include building and human factors.Among them, the influence degree and adjustability of firefighting facility factors are relatively large, which can be improved through optimized design to improve evacuation efficiency.In addition, the extraction and rescue strategies were given from the following two aspects according to the above simulations.

(A) Ventilation and smoke evacuation:
(1) To ensure that the stratification of high-temperature smoke is not destroyed during the safe evacuation stage, the airflow velocity near the fire point should not be too large; (2) In the fire suppression stage, the ventilation velocity in the underground interchange must be greater than the critical ventilation velocity of fire to ensure that the firefighters can safely arrive at the fire site from the upwind side of the underground interchange; ( 3) The open fans should be increased when the smoke is upstream of the fire countercurrent.When the downstream smoke speed of the fire source is too fast or is obviously blown away, the number of fans should be reduced, the air supply or exhaust system should be started after the evacuation of personnel in the tunnel, and the smoke exhaust speed should be increased.
(B) Personnel guidance and evacuation: (1) In the case of fire, there is usually no guide for an underground interchange.Thus, it is necessary to set up eye-catching guide signs in the underground interchange; (2) In the early stage of tunnel fires, personnel evacuation should be given priority, and fire suppression and smoke exhaust measures should be carried out after the evacuation of personnel; (3) When a traffic jam inside the local underpass is serious, firefighting vehicles should drive into the adjacent tunnel and enter the transverse passage entrance near the fire-starting point for rescue and fire suppression.

Conclusions
A method of combining BIM and numerical simulation was performed to conduct the fire protection and evacuation analysis in large-scale underground interchange tunnels.The following are the main conclusions: (1) The peak values of CO concentration, visibility, and temperature all appeared above the fire source.The CO concentration and temperature downstream of the fire source increased, and the visibility decreased with increasing ventilation velocity.Therefore, taking a low ventilation speed is conducive to the evacuation of personnel from the fire downstream in the escape stage; the critical ventilation velocity during the rescue phase is conducive to the entry of fire personnel from upstream of the fire source; (2) When the critical ventilation velocity was 3.6 m/s, the ventilation volume in the tunnel met the requirement of no adverse flow of smoke in the longitudinal smoke evacuation mode.As well as meeting the requirements for the safe evacuation of vehicles and personnel, it can ensure the safety of fire sources upstream and downstream without causing tunnel congestion.After the fire, the longitudinal ventilation velocity of the non-fire side can be appropriately increased to prevent flue gas backflow in the main tunnel; (3) Visibility is the main factor affecting escape.When obstacles were placed at the crosspassage entrance, these two increased by 20% and 8.7%, respectively, which was mainly due to the reduced visibility caused by smoke diffusion.A 2.5-m-wide transverse passage was effective in reducing the escape time.The patency of the fire exits has an obvious effect on the escape time.
The combination of BIM technology and fire escape software provides a new idea for the fire evacuation research of underground interchanges in the future.However, there

Figure 1 .
Figure 1.Technology roadmap of BIM modeling and fire simulation in underground interchange tunnels.

Figure 2 .
Figure 2. Simplified model diagram of underground interchange in Pyrosim software.

Figure 3 .
Figure 3. Schematic diagram of flue gas development without longitudinal ventilation: (a) top view; (b) front view.

Figure 7 .
Figure 7. Simplified model diagram of underground interchange in Pathfinder software.

Figure 8 .
Figure 8. Changes in the longitudinal center section of the CO concentration: (a) 0.5 m/s; (b) 3.6 m/s.

3. 1
Figure12illustrates the temperature (T) changes in the longitudinal center section at different distances.In Figure13, T is plo ed at 0.5 m/s longitudinal ventilation velocity and 3.6 m/s critical ventilation velocity.


Working condition A: the ventilation velocity of the left and right transverse passageways was 0.5 m/s. Working condition B: the ventilation velocity of the left and right transverse passageways was 1.0 m/s. .

Figure 17 .
Figure 17.Results of evacuation simulation of working condition #1: (a) number of people passing through doors and exits; (b) relationship between the number of people and time in different areas.

Figure 18 .
Figure 18.Results of evacuation simulation of working condition #2: (a) number of people passing through doors and exits; (b) time and number of people in different areas.

4. 2 . 1 .
Effect of the Fire Escape Patency on the Evacuation Simulation (1) Internal patency of the transverse channel

Figure 21 .
Figure 21.Results of evacuation simulation of heavy congestion.

Table 1 .
Escape time under different working conditions.