1. Introduction
Considerable research has been carried out over the last few decades to better understand the factors that affect the occurrence of traffic accidents in road tunnels. These studies have led to the development of several statistical models to estimate the expected value of accident frequency (i.e., number of crashes per year) [
1,
2,
3,
4,
5,
6] as a function of geometry, traffic flow, and certain types of equipment implemented in tunnels.
The occurrence of traffic accidents in a road tunnel might also temporarily compromise its functionality over time due to the partial or complete closure of the tunnel (e.g., only one lane or all the lanes of the tunnel are closed), as well as this potentially causing negative effects on the service level of the traffic flow of the road network containing the tunnel. The duration of the disruptive event in the tunnel may depend on how serious it is. When there is a severe accident (i.e., involving fatalities and/or injuries), it might also take up to several hours before the tunnel is available again for traffic flow. This may be due to the time necessary for the rescue team to arrive, the road police to investigate the causes on site, the local emergency services to help the injured or take away fatalities, and the tunnel management staff to remove damaged vehicles and clean the road of any debris.
Given the essential role that road tunnels play in a road network in supporting the mobility of people and goods, following a traffic accident, the non-availability of a tunnel for the passage of traffic flow through it is a relevant issue for the social and economic interests of the community.
Road tunnels that present an adaptive response to a disruptive event (e.g., a traffic accident) are generally described as resilient. There is no shared concept of resilience in the current literature. However, with respect to infrastructure networks, resilience may be described as “the ability of a system to rapidly recover its functionality after a disturbing event” [
7].
Resilience of a transportation system has been investigated from different perspectives. Some researchers, for example, have discussed the most appropriate metrics for assessing the resilience of an infrastructure network [
7,
8,
9,
10]. Knoop et al. [
11], and Liao et al. [
12] focused, instead, on traffic flow simulation models to investigate both the robustness of a road system as well as the effectiveness of different management strategies. Omer et al. [
13], and Antoniou et al. [
14] proposed frameworks to evaluate the resilience of a road network. Reggiani et al. [
15] concentrated on the resilience and vulnerability of transportation networks. Kaviani et al. [
16] found the optimal location of road guidance tools for re-routing vehicular flow after a disturbing event. Amini et al. [
17] developed a methodology to evaluate the effectiveness of diversion on the resilience of a road network. Xu et al. [
18] investigated the redundancy of a transportation network, while Sohouenou and Neves [
19] concentrated on both the robustness of a network and the rapidity of the recovery process. Abudayyeh et al. [
20] proposed a method for setting traffic signals in order to improve resilience. Khetwal et al. [
21] set up a stochastic simulation model to evaluate tunnel resilience, while Zhao [
22] developed a genetic algorithm for optimizing a road network. PIARC [
23] presented the current state of knowledge of road tunnel resilience, whereas Borghetti et al. [
24] investigated the resilience and emergency management of road tunnels. Finally, Caliendo et al. [
25] set up a traffic simulation model for evaluating the resilience of a twin-tube motorway tunnel in the case of a traffic accident in a tube.
However, it is worth highlighting that, among the studies referred to, only a few have focused specifically on the resilience of road tunnels, such as Liao et al. [
12], Khetwal et al. [
21], PIARC [
23], Borghetti et al. [
24], and Caliendo et al. [
25].
Moreover, as far as the authors of this paper are aware, there are no studies in the current literature in which the restoration of the functionality of a tunnel (i.e., resilience) in the event of a traffic accident is also associated with user safety that might be compromised, for example, by a fire of vehicles transiting on the remaining undisrupted lanes of the partially closed tunnel. Furthermore, in the case of a twin-tube road tunnel in which a tube is completely closed due to a traffic accident, the resilience of the tunnel system has not been evaluated in combination with occupant safety that may be affected by a fire affecting vehicles that, having changed carriageway, use one lane of the adjacent tube to travel (i.e., where bi-directional traffic conditions are generated in the parallel tube). This represents a lacuna of knowledge that this paper seeks to fill.
With reference to our previous studies, it is also worth highlighting that we did not investigate user safety in the earlier paper [
25], while in other studies, we evaluated only the safety level of tunnel users in the event of a fire (e.g., Caliendo et al. [
26,
27,
28,
29,
30,
31,
32], Caliendo and De Guglielmo [
33,
34,
35,
36], and Caliendo and Genovese [
37]), but not resilience. In other words, tunnel resilience and user safety were not coupled to each other; consequently, the effects of a strategy to recover the functionality of a tunnel that might influence user safety or vice versa were not considered. A simultaneous analysis of user safety and resilience now represents an important topic of discussion in the field of tunnel engineering, in particular for tunnel management agencies (TMAs). This represents a primary justification for the present research.
This study complements our previous research studies by coupling user safety and tunnel resilience, to provide additional perspectives on relevant factors, as well as to suggest certain recovery strategies in order to enhance the resilience of tunnels without neglecting user safety. Moreover, it facilitates the dissemination of additional information for road tunnel management.
A twin-tube motorway tunnel was investigated. We developed two computational fluid dynamics (CFD) models for scenarios of the partial or complete closure of a tube; then, the corresponding evacuation processes were simulated to evaluate exposure risk to users due to different types of fire. The results obtained were applied in a quantitative risk analysis (QRA), based on a probabilistic approach, to build F-N curves for each scenario investigated. Comparisons between the F-N curves and tunnel resilience were carried out to optimize the choice of recovery strategy in relation to user safety.
Given that it was not obvious how the implementation of certain recovery measures for improving the tunnel resilience might simultaneously combine with user safety in the event of a fire, we reviewed our assumptions with a view to obtaining a deeper understanding of both safety and resilience at the same time.
In the light of the above considerations, it needs to be stressed that the main objective of the present study was to couple a resilience analysis of twin-tube road tunnels with a risk assessment for users in the event of a fire that might occur during the partial or complete closure of a tube. The findings might be useful for understanding whether the recovery strategies implemented to restore the functionality of the system in the event of a traffic accident in a road tunnel are suitable also from the point of view of user safety.
The scientific significance of the present study lies in the improvement of knowledge in the field of both fire safety engineering and tunnel management, by providing additional insights on both user safety and tunnel resilience, and in the dissemination of further information for the operating conditions of tunnels in the event of emergency.
The paper is organized as follows: the next section contains a description of the tunnel investigated. Subsequently, the computational fluid dynamics (CFD) modeling and people egress process are presented and implemented. Then, the quantitative risk analysis (QRA) is described, and the frequency (F) of potential fatalities (N) is quantified. The results are then reported and commented on, and comparisons are made with tunnel resilience for the different scenarios investigated. Finally, some conclusions are drawn and recommendations for practical applications are made, and future research possibilities are discussed.
2. Materials and Methods
2.1. Tunnel Description
A twin-tube motorway tunnel situated along an Italian road network is examined. Each tube has a length of 850 m and is used for one-way traffic under ordinary conditions of service. The tubes are flat, straight, and have no emergency lanes. Each tube is characterized by a horseshoe-shaped cross-section of 55 m2, with a total width of 9.5 m (i.e., two traffic lanes of 3.5 m, and two sidewalks of 1 m) and a maximum height of 6.8 m. An emergency exit, which consists of a cross-connection between the two tubes, is located in the middle of the tunnel length.
2.2. Scenarios Investigated
In our previous study [
25], we assumed that the occurrence of a traffic accident in the north tube would lead to different potential scenarios involving its partial or complete closure for 1 h, 2 h, or 3 h (however, in this investigation, we considered the worst-case scenario, which was a closure of 3 h). Specifically, the functionality of the investigated twin-tube motorway tunnel, after the occurrence of a traffic accident in the north tube, was assumed to be restored by utilizing the remaining available lane of the north tube (i.e., only one lane is closed) or rearranging the traffic flow using the south tube for two-way traffic (i.e., both lanes of the north tube are closed). The effects of the presence of an alternative route, along which to divert heavy goods vehicles (HGVs), were also investigated.
In this study, we investigated the same scenarios assuming that, during the partial or complete closure of the north tube due to the occurrence of a traffic accident in it, a fire had simultaneously occurred on the undisrupted remaining lane of the north tube, or in the parallel tube used for two-way traffic, in the case of the partial or complete closure of the north tube.
Figure 1 shows a cross-section of the twin-tube motorway tunnel examined.
Figure 2 shows the scenarios investigated: (i) Scenario 0 represents the case in which the tunnel system is used under its ordinary conditions of functionality (i.e., without any traffic accident and with both lanes of each tube used for one driving direction) when a fire occurs in the north tube; (ii) Scenario Ia corresponds to the case in which, assuming that the north tube is partially closed due to a traffic accident on the right lane, a fire occurs on the undisrupted remaining lane of the north tube (i.e., the left lane); (iii) Scenario Ib is similar to Scenario Ia, but also includes the use of an alternative route along which, by using variable message signs (VMSs), only HGVs traveling towards the north are diverted; (iv) Scenario IIa corresponds to the case of the complete closure of the north tube and the use of the parallel tube for two-way traffic; here, the fire is assumed to occur on the lane used by the traffic flow traveling towards the north; (v) Scenario IIb is similar to Scenario IIa, but it also includes the use of an alternative route along which, using VMSs, only HGVs traveling towards the north are redirected. With reference to Scenarios IIa and IIb, the two traffic by-passes located at the tunnel portals, which allow for a change of carriageway, are assumed to be opened by the emergency service team in a very short time (i.e., 10 min) from the start of the traffic accident. The nodes E and I represent the two motorway junctions which are before the portals of the twin-tube tunnel that was investigated. Specifically, node E is used by HGVs when the alternative itinerary is recommended for these vehicles only.
2.3. Hourly Traffic Volume
The traffic flow, expressed in terms of hourly traffic volume (VHP), transiting through the north tube during its ordinary conditions of functionality (i.e., Scenario 0) was found to be—according to the traffic database of the road management agency—equal to 2200 vehicles/h (i.e., 1100 vehicles/h per lane). However, it should be noted that the partial or complete closure of the north tube due to a traffic accident might lead to traffic congestion on its remaining undisrupted lane or in the adjacent tube used for two-way traffic, respectively. On this basis, the VHP was assumed to be equal to the capacity of a one-way lane (i.e., 1700 vehicles/h) or to that of a bi-directional carriageway having two lanes (i.e., 3200 vehicles/h) [
38] for the partial (i.e., Scenarios Ia and Ib) or complete (i.e., Scenarios IIa and IIb) closure of the north tube, respectively. The heavy vehicle percentage (i.e., buses and HGVs), which was obtained from the same traffic database, was always equal to 25%.
2.4. Types of Burning Vehicles
The effects of five different categories of burning vehicles were investigated. Each burning vehicle was assumed to be placed in five different positions along the right lane of the north tube in Scenario 0 (i.e., 145, 280, 420, 570, and 710 m from portal A). The burning vehicles were assumed to be two cars, a van, a bus, and two different types of HGVs capable of developing a maximum heat release rate (HRRmax) of 8, 15, 30, 50, and 100 MW, respectively. The same five different types of burning vehicles and locations were considered on the remaining undisrupted lane (i.e., the left lane) of the north tube, partially closed due to a traffic accident in Scenario Ia; while in Scenario Ib only the occurrence of a fire involving light vehicles was simulated (i.e., characterized by an HRRmax of 8 and 15 MW) since the HGVs were redirected to the alternative route.
In the case of the complete closure of the north tube and the use of the parallel tube for two-way traffic, the aforementioned five different types of burning vehicles were assumed to be located on the lane used by traffic traveling towards the north and at the aforementioned positions (i.e., 145, 280, 420, 570, and 710 m from portal D) in Scenario IIa; while in Scenario IIb, only light vehicles were considered (i.e., HRRmax of 8 and 15 MW) given that the HGVs were diverted to the alternative itinerary.
2.5. Natural Ventilation
The two tubes are not equipped with a mechanical ventilation system. Therefore, the longitudinal ventilation along each tube is natural and due to the piston effect of moving vehicles. In this respect, the natural ventilation was simulated by setting an appropriate pressure difference between the tunnel portals. Specifically, a positive pressure difference of 5 Pa [
29] was applied between the entrance portal (portal A) and the exit portal (portal B) of the north tube when used for unidirectional traffic towards the north.
In the case of the complete closure of the north tube and the use of the parallel tube for two-way traffic, due to the two opposite traveling directions in the tube, the piston effect is expected to be almost negligible. However, since the speed of the vehicles that changed tube is expected to be slightly lower than that of the vehicles traveling in the opposite direction, a negative pressure difference of −0.5 Pa [
30] was set for traffic flow traveling towards the north (i.e., from portal D to portal C of the adjacent tube).
2.6. Concrete Characteristics
The tunnel walls and ceiling were assumed to be made of ordinary concrete with a thickness of 0.4 m. The properties of this material, which were taken from the literature, are as follows: thermal conductivity of 1.67 W/m/K, specific heat of 0.94 kJ/kg/K, density of 2585 kg/m
3, and emissivity of 0.9 [
39].
2.7. Results of Resilience Analysis
Resilience analysis of a transportation system is generally performed using traffic simulation models [
40,
41,
42,
43,
44]. In our previous study [
25], the level of functionality of the twin-tube motorway tunnel, during the period in which the north tube was partially or completely closed due to a traffic accident, was calculated as the ratio between the average travel time before and after the traffic accident. We calculated the reduction in the functionality F* (%) for all the scenarios examined (the functionality level F* was assumed to be equal to 100% in absence of incidents). Then, we built the functionality curves vs. time (i.e., F* (%) − t) and for all the scenarios computed the resilience loss:
, in which t
o is the time instant when the traffic accident occurs in the north tube and t
h is the time instant when the twin-tube tunnel is again fully functional. Smaller resilience losses were found with the partial closure of the north tube rather than the complete one. However, when an alternative route was activated for HGVs only, both in the event of the partial and of the complete closure of the north tube, minor reductions of resilience losses were found compared to those without considering HGV rerouting.
In order to describe the resilience through a parameter, whose value ranges between 0 and 1, we computed the resilience index ( ) that is widely used in the literature for transportation networks. A low resilience index denotes a high resilience loss (i.e., a resilience index equal to one corresponds to a zero-resilience loss). The resilience index confirmed the abovementioned findings of the resilience loss in all the scenarios investigated.
Given these findings, the intent of this paper was to assess the impacts on user safety related to each of the mentioned fire scenarios in order to help tunnel management agencies (TMAs) decide which arrangement is more appropriate to adopt in the case of a traffic accident in a tube of a twin-tube road tunnel.
2.8. Research Framework
This study is set in the context of research on the simultaneous analysis of user safety and resilience of road tunnels, but expands the state-of-the-art by setting up a computational fluid dynamics (CFD) model to assess the risk exposure of users during their escape from the tunnel in the event of a fire occurring after a traffic accident had occurred, as well as estimating the corresponding potential number of victims with a probabilistic-based quantitative risk analysis (QRA). Therefore, given the lack of studies in which resilience is associated with user safety in road tunnels, this paper can serve as a possible reference for tunnel management agencies (TMAs), and a means of improving our knowledge in the field of the interaction between the recovery of the functionality of tunnels and user safety.
The methodology applied is briefly represented in
Figure 3.
5. CFD Results
Considering the worst-case combination of HRRmax (i.e., 100 MW for Scenarios 0, Ia, and IIa, and 15 MW for Scenarios Ib and IIb) and the location of the fire (i.e., in the middle of the tunnel tube length) for each scenario investigated, the longitudinal profiles of temperature, radiant heat flux, visibility distance, and toxic gases concentration predicted at breathing height (2 m) along the escape route (i.e., the nearest sidewalk) are reported in this section and compared with the acceptable tenability limits in order to verify if the environmental conditions were acceptable for human health.
5.1. Temperature
Figure 5 shows the longitudinal temperature profiles predicted at 2 m height along the escape path after 10 min from the start of the 100 MW (i.e., Scenarios 0, Ia, and IIa) or 15 MW (i.e., Scenarios Ib and IIb) fire located in the middle of the tunnel tube length, as well as the positions occupied by the last escaping user after 10 min of fire exposure both upstream (for all scenarios) and downstream (for Scenarios IIa and IIb only) of the burning vehicle. It is worth noting how the aforementioned evacuee, according to his/her effective walking speed and pre-movement time, reaches a distance upstream of the fire of about 371 m (i.e., walking speed of 0.7 m/s, the red user) for Scenarios 0, Ia, and Ib (escape route only towards portal A of the north tube), as well as for Scenario IIb both upstream and downstream of the fire (i.e., escape route in the direction of both portal D and portal C of the adjacent tube). In contrast, for Scenario IIa, due to a reduced walking speed (i.e., about 0.55 m/s), he/she (the green user) reaches, after 10 min from the start of the fire, a shorter distance from the burning vehicle (i.e., 295 m), which is due to the poor visibility and high CO concentration (see
Section 5.3 and
Section 5.4, respectively). However, the mentioned user is always exposed to a temperature ≤60 °C (i.e., dotted line), which is the tenability limit according to UPTUN [
50].
5.2. Radiant Heat Flux
Figure 6 shows the longitudinal profiles radiant heat flux predicted at 2 m height along the escape path after 10 min from the start of the 100 MW (i.e., Scenarios 0, Ia, and IIa) or 15 MW (i.e., Scenarios Ib and IIb) fire located in the middle of the tunnel length.
Regardless of the scenario considered,
Figure 6 shows how the last escaping user is always subjected to a radiant heat flux lower than the tenability limit of 2 kW/m
2 (i.e., dotted line) [
50].
5.3. Visibility Distance
Figure 7 shows the longitudinal profiles of visibility distance predicted at 2 m height along the escape path after 10 min from the start of the 100 MW (i.e., Scenarios 0, Ia, and IIa) or 15 MW (i.e., Scenarios Ib and IIb) fire located in the middle of the tunnel length.
Figure 7 shows how the last escaping user always has a visibility distance greater than the tenability limit of 10 m (i.e., dotted line) [
50], except for the 100 MW fire related to Scenario IIa in which he/she (the green user) might be exposed to unsafe environmental conditions due to his/her inability to discern evacuation signs located along the escape route.
5.4. CO Concentration
Figure 8 shows the longitudinal CO concentration profiles predicted at 2 m height along the escape path after 10 min from the start of the 100 MW (i.e., Scenarios 0, Ia, and IIa) or 15 MW (i.e., Scenarios Ib and IIb) fire located in the middle of the tunnel length. Based on these profiles, the last escaping user is exposed to a CO concentration lower than the tenability limit of 1200 ppm (i.e., dotted line) [
51] in Scenarios 0, Ia, Ib, and IIb, while he/she (the green user) might be exposed to unsafe environmental conditions in Scenario IIa.
5.5. CO2 Concentration
Figure 9 shows the longitudinal profiles of CO
2 concentration predicted at 2 m height along the escape path after 10 min from the start of the 100 MW (i.e., Scenarios 0, Ia, and IIa) or 15 MW (i.e., Scenarios Ib and IIb) fire located in the middle of the tunnel length.
Figure 9 shows how the last escaping user is always exposed to a CO
2 concentration below the tenability limit of 40,000 ppm (i.e., dotted line) [
51].
Summing up the results above, with reference to the last evacuating user: (i) better environmental conditions were observed to occur with partial closure of the tube rather than complete closure (i.e., Scenarios Ia and Ib against Scenarios IIa and IIb); (ii) additional benefits can be found by activating VMSs that suggest an alternative itinerary for HGVs only (i.e., Scenario Ib versus Scenario Ia, and Scenario IIb against Scenario IIa); this is due to the fact that when the HGVs are diverted to the alternative route, the most potentially dangerous event is represented only by a van fire (i.e., HRRmax of 15 MW); (iii) safety issues for human health may arise only in the event of 100 MW fire occurring during the complete closure of the north tube and the use of the parallel one for two-way traffic (i.e., Scenario IIa) due to exceeding the acceptable limit of 10 m and 1200 ppm for the visibility distance and CO concentration, respectively.
It is worth highlighting that similar evacuation conditions were also found for the remaining locations of the burning vehicle (i.e., 145, 280, 570, 710 m from portal A of the north tube for Scenarios 0, Ia, and Ib, or from portal D of the adjacent tube used for bi-directional traffic for Scenarios IIa and IIb). With reference to the remaining HRRmax investigated (i.e., 8, 15, 30, and 50 MW fires for Scenarios 0, Ia, and IIa; or 8 MW fire for Scenarios Ib and IIb) and locations of the fire in the tunnel tube, we did not find any unsafe conditions for the users evacuating the tube. However, it needs to be stressed that, in order to save space, in the present paper, we have reported only the results of the worst-case scenarios that were found to be the burning vehicles with HRRmax = 100 MW located in the middle of the tunnel tube length for Scenarios 0, Ia, and IIa, and for the fire of 15 MW in Scenarios Ib and IIb.
5.6. Back-Layering Phenomenon and Smoke Stratification
Figure 10 shows the back-layering phenomenon and smoke stratification predicted along the escape route after 10 min from the start of the fire of 100 MW (i.e., Scenarios 0, Ia, and IIa) or 15 MW (i.e., Scenarios Ib and IIb) located in the middle of the tunnel length. It shows that the last escaping user is always exposed to safe evacuation conditions, except for the 100 MW fire of Scenario IIa in which the aforementioned user (i.e., the green user) might be a victim of serious risks to human health due to smoke spread.
Figure 10 shows that: (i) the natural ventilation due to the piston effect of moving vehicles is not able to prevent the back-layering phenomenon (i.e., the smoke propagation in the opposite direction of the natural ventilation); (ii) with increase in HRR
max (i.e., the 100 MW for Scenario Ia in contrast to the 15 MW for Scenario Ib, and the 100 MW for Scenario IIa versus the 15 MW for Scenario IIb), worse smoke conditions are found on the left sidewalk of the north tube and both sidewalks of the parallel tube used for two-way traffic, respectively; (iii) the smoke layering length is not symmetrical for Scenarios 0, Ia, and Ib due to the natural ventilation that pushes the smoke towards portal B (i.e., pressure difference of 5 Pa from portal A to portal B); (iv) the smoke distribution in the adjacent tube (i.e., Scenarios IIa and IIb) is almost symmetrical in the two directions upstream and downstream of the fire because of the reduced pressure difference between the tunnel portals (i.e., 0.5 Pa from portal C to portal D, or equivalently −0.5 Pa from portal D to portal C, that very slowly pushes the air flow towards portal D).
Therefore,
Figure 10 also confirms, with respect to smoke layering, that better environmental conditions are found with the partial closure of the tube rather than complete closure, and that additional advantages can be found through the activation of VMSs that suggest an alternative route for HGVs only.
6. Quantitative Risk Analysis
6.1. Methodological Approach
A quantitative risk analysis (QRA) based on a probabilistic approach was carried out to estimate the consequences in terms of the number of potential victims associated with the scenarios simulated by means of the FDS + Evac modeling. The main steps for implementing a probabilistic approach, which involves an event tree, fault tree, and models for estimating the consequences, entail the choice of objective, specification of the tunnel, classification of each hazard and evaluation of its probability (or frequency) of occurrence, evaluation of the consequences on user health, estimation of the risk as the sum of the probabilities multiplied by the consequences, choice of the risk acceptability criterion, and quantification of the risk level.
The so-called societal risk is generally the main outcome of a QRA and is commonly expressed graphically in terms of F-N curves (F is the annual cumulative frequency of potential fatalities N) or numerically in terms of the expected value (EV) of risk (). The aforementioned F-N curves might be used to compare the risk level of a given tunnel under different scenarios and/or with the unacceptable and acceptable risk limits of the ALARP (As Low As Reasonably Practicable) region.
6.2. Event Tree
The identification of each hazard (i.e., fires characterized by an HRR
max of 8, 15, 30, 50, or 100 MW) and the estimation of its probability of occurrence are made using the event tree. It should be noted that the event tree is designed in a graphic form to depict several chronological series of successive events starting from an initial event. Specifically, it reports the probability (p) of the occurrence of each intermediate event based on the probability of the initial event. In this paper, we used the event tree proposed by ANAS [
52]. The initial event consists of the annual frequency of traffic accidents in road tunnels, which is followed by: traffic collisions (p = 94.9%) and fires (p = 5 %). Fires might involve light vehicles (p = 70%) or heavy vehicles (p = 30%) and are classified as relevant (p = 70% for light vehicles, or p = 20% for heavy vehicles) or non-relevant (p = 30% for light vehicles, or p = 80% for heavy vehicles). The relevant fires involving light vehicles might entail an HRR
max of 8 MW (p = 2.5%), while those involving heavy vehicles might be characterized by an HRR
max of 15, 30, 50, and 100 MW (p = 81.5, 14.5, 2.5, and 1.5%, respectively).
6.3. Annual Frequency of Traffic Accidents
To estimate the annual frequency of traffic accidents, which represents the initial event in the event tree, the average accident rate (i.e., the average number of accidents per 100 million vehicle-kilometers) of 43.5 accidents/108 vehicle-km provided by the tunnel management agency (TMA) for the investigated tunnel was considered. It is worth noting that the risk analysis was carried out considering the partial or complete closure of the north tube for three hours due to the occurrence of a traffic accident on the right lane of this tube, since this was the worst-case scenario in our previous study on tunnel resilience.
Therefore, by considering the three hours of peak-traffic volumes (see
Section 2.3), we computed the corresponding equivalent annual average daily traffic (EAADT) for each scenario investigated: 6600 vehicles/day (i.e., 2200 vehicles/h for the two unidirectional lanes) for Scenario 0; 5100 vehicles/day (i.e., 1700 vehicles/h for one lane) for Scenarios Ia and Ib; and 9600 vehicles/day (i.e., 3200 vehicles/h for the two bi-directional lanes) for Scenarios IIa and IIb.
Based on the EAADTs and accident rate, the annual frequency of traffic accidents (i.e., the average number of accidents per year) was found to be: 0.891 for Scenario 0; 0.688 for Scenarios Ia and Ib; and 1.3 for Scenarios IIa and IIb. These values were implemented as the initial event in the event tree.
6.4. Annual Fire Frequency
The annual frequency of each fire scenario reproduced through the FDS + Evac codes was estimated as the product of the annual frequency of traffic accidents and the probabilities of all the intermediate events included between the initial and final events in the event tree (e.g., with reference to the 100 MW fire for Scenario IIa, we calculated this as: 1.3 × 0.05 × 0.3 × 0.2 × 0.015 = 5.85 × 10
−5).
Table 2 reports the annual frequency of occurrence of each scenario investigated.
Moreover, each of the five positions at which the burning vehicle was assumed to be located along the tunnel length (i.e., 145, 280, 420, 570, and 710 m from portal A of the north tube for Scenarios 0, Ia, and Ib, or portal D of the parallel tube used for two-way traffic for Scenarios IIa and IIb) was given the same probability of occurrence, thus equal to 0.2. Consequently, the values reported in
Table 2 were multiplied by 0.2 to estimate the annual frequency of occurrence of a certain fire in each mentioned location of the burning vehicle in the tunnel.
6.5. F-N Curves
In this paper, the fractional effective dose (FED) was used to calculate the number of people at risk in the event of a fire in the tunnel tube investigated. Based on the FED concept, an evacuee is considered as a potential victim when the FED corresponding to his/her exposure to temperatures and radiant heat fluxes (FED
heat) and/or toxic substances (FED
toxic gases) exceeds unity (for more details see Caliendo et al. [
29]). It should be noted that the QRA was performed by considering the FED of all the evacuees and their entire egress process from the tunnel, which occurs using: (i) both the right and left sidewalk upstream of the burning vehicle for Scenario 0 (i.e., users escaping towards portal A); (ii) the left sidewalk upstream of the fire source for Scenarios Ia and Ib (i.e. users escaping towards portal A); (iii) the right sidewalk upstream of the fire for the users escaping towards portal D, and the right sidewalk downstream of the fire for the users escaping towards portal C for Scenarios IIa and IIb.
This was carried out for each HRRmax (i.e., 8, 15, 30, 50, and 100 MW for Scenarios 0, Ia, and IIa, and 8 and 15 MW for Scenarios Ib and IIb) and location (i.e., 145, 280, 420, 570, and 710 m from portal A of the north tube for Scenarios 0, Ia, and Ib, or portal D of the adjacent tube for Scenarios IIa and IIb) of the fire. For some of the combinations of HRRmax and location of the burning vehicle, the FEDtoxic gases was computed to be greater than 1, while the FEDheat was always found to be much less than 1. As a result, the estimated number of people at risk derives only from their exposure to toxic gases.
Then, the number of potential victims, referred to each HRRmax and location of the fire, was associated with the corresponding annual frequency of fire occurrence in order to compute the individual cumulative frequency of each fire. Subsequently, by combining the individual cumulative frequencies, we calculated the final cumulative frequency with respect to the escaping users from one of the tunnel portals and/or the emergency exit.
Figure 11 shows, in a bi-logarithmic diagram, the F-N curves for each scenario investigated, as well as the unacceptable and acceptable limits of the Italian ALARP region.
Figure 11 shows how all the F-N curves lie below the unacceptable limit of the ALARP region. The partial closure of the north tube presents a lower risk level than that associated with the tunnel in its ordinary conditions of functionality (i.e., Scenarios Ia and Ib against Scenario 0). This is because Scenarios Ia and Ib are characterized by a lower traffic volume (1700 vehicles/h for Scenarios Ia and Ib against 2200 vehicles/h for Scenario 0) and consequently by a lower number of users exposed to risk (only one lane is available to traffic before the start of the fire) than Scenario 0. In contrast, the F-N curves corresponding to the scenarios involving the complete closure of the north tube and the use of the parallel one for two-way traffic (i.e., Scenarios IIa and IIb) always lie above that of Scenario 0. This might be explained by the higher traffic volume (3200 vehicles/h for Scenarios IIa and IIb against 2200 vehicles/h for Scenario 0) and thus the greater number of potential victims expected in the adjacent tube when it is employed for two-way traffic.
Figure 11 also shows that, both in the case of the partial and of the complete closure of the north tube, certain benefits in terms of users’ safety, might be found by activating VMSs, situated at a reasonable distance from the motorway junction before the entrance portal of the north tube, for diverting only HGVs towards an alternative itinerary (i.e., Scenario Ib versus Scenario Ia, and Scenario IIb against Scenario IIa). These benefits are because the fires, involving only cars and vans traveling towards the north, are characterized by a reduced HRR
max compared to those of HGVs. Moreover, Scenario Ib (i.e., partial closure of the north tube with HGVs redirected to the alternative itinerary) presents the lowest level of risk among the scenarios investigated.
6.6. Resilience Index and Expected Value of Risk for Each Scenario
Figure 12 summarizes the results of the simultaneous analysis of user safety and resilience reported in Caliendo et al. [
25] for each scenario investigated expressed in terms of the expected value (EV) of risk and resilience index, respectively. It should be noted that the EV of risk represents the area subtended by the F-N curve (i.e.,
).
Figure 12 shows how the lowest value of risk, expressed by the EV of risk, is found in the case of the partial closure of the north tube when the alternative itinerary is activated for HGVs only (i.e., Scenario Ib), and that the associated resilience index is only a little lower than 1 (i.e., 0.89), which indicates a small resilience loss. The EV of risk increases when the north tube is partially closed with its left lane used by all the vehicles, while the corresponding resilience index decreases from 0.89 to 0.77 (i.e., Scenario Ia presents a slightly higher resilience loss than Scenario Ib). However, in both cases, the value of the EV of risk is less than that corresponding to the ordinary conditions of functionality of the north tube (i.e., Scenario 0).
In contrast, with reference to the complete closure of the north tube, the values of EV of risk are found to be higher than that of Scenario 0, and the resilience index decreases up to 0.66 and 0.50 for Scenarios IIb and IIa, respectively, which indicates a more significant resilience loss.
7. Summary and Conclusions
The present study was mainly motivated by the need to simultaneously couple user safety with resilience in order to identify the most appropriate arrangement to be adopted in the event of a traffic accident in a tube of a twin-tube road tunnel. In this respect, we investigated different scenarios by assuming that during the partial or complete closure of a tube, because of the occurrence of a traffic accident, a fire might simultaneously occur on the undisrupted remaining lane, or in the parallel tube used for two-way traffic in the case of the partial or complete closure of the disrupted tube, respectively.
For this purpose, we developed a computational fluid dynamics (CFD) model, and simulated the corresponding egress process, to evaluate the risk level of tunnel users in the event of a fire. The Fire Dynamics Simulator (FDS) code, as well as its evacuation module (Evac), were implemented, and their output were used to carry out a quantitative risk analysis (QRA) based on a probabilistic approach.
The effects of five different types of burning vehicles were investigated. Each burning vehicle was assumed to be located at five different positions along the remaining undisrupted lane, or in the parallel tube used for two-way traffic in the case of the partial or complete closure of the tube affected by the traffic accident. The burning vehicles were assumed to be two cars, a van, a bus, and two different types of HGVs capable of developing a maximum heat release rate (HRRmax) of 8, 15, 30, 50, and 100 MW, respectively.
With reference to the last user escaping from the tunnel tube in the event of a fire, the findings of the CFD simulations, in terms of temperature, radiant heat flux, visibility distance, and toxic gases concentration at breathing height along the escape route, showed that: (i) better environmental conditions were found with partial closure of the tube rather than complete closure; (ii) additional benefits can be found by activating variable message signs (VMSs) that suggest an alternative itinerary for HGVs only; (iii) safety issues for human health may arise only in the case of a 100 MW fire, occurring during the complete closure of the north tube and the use of the parallel one for bi-directional traffic, due to exceeding the acceptable tenability limits of both visibility distance and CO concentration.
The results of the QRA were found to be consistent with those obtained by CFD modeling, as follows: (i) the partial closure of the tube presents a lower risk level than that associated with the tube in its ordinary conditions of functionality; (ii) the complete closure of a tube and the use of the parallel tube for two-way traffic leads to a higher risk level compared to that of the ordinary conditions of functionality of the tunnel; (iii) additional reductions in risk can be found by activating VMSs to divert only HGVs towards an alternative itinerary.
The lowest value of risk (expressed by the expected value (EV) of risk), which corresponded to the case of the partial closure of the north tube and the activation of an alternative route for HGVs only, was found to be associated with a resilience index slightly less than one, which showed a small resilience loss. The EV of risk increases when the north tube is partially closed with its left lane used by all the vehicles, and the corresponding resilience index was found to be slightly decreased. However, in both cases the value of the EV is less than that corresponding to the ordinary conditions of functionality of the tunnel tube. In contrast, with reference to the complete closure of the north tube, the values of the EV of risk were found to be higher than that of the tunnel in its ordinary conditions of functionality, while the resilience index decreased, which indicated a more significant resilience loss.
The scientific significance of the present paper lies in increasing our knowledge in the research area of both fire safety engineering and tunnel management, and in the dissemination of additional information for the operating conditions of tunnels.
This study provides valuable perspectives on both relevant aspects of user safety and the resilience of road tunnels at the same time. It can serve as a possible reference for tunnel management agencies (TMAs) in the choice of the most appropriate strategy for the recovery process of system functionality in the event of an incident.
Although the present paper represents an advancement in the context of research on the simultaneous analysis of user safety and tunnel resilience, the randomness of certain parameters might create some uncertainties. Moreover, disruption in a road tunnel might also propagate to other transport systems. Our modeling is likely to be generally valid for tunnels having a length ≤ 1 km, which are by far the most frequent in Italy. However, additional studies should also involve much longer tunnels which use innovative equipment that may be automatically activated when a fire starts, as an alternative or supplementary measure to hydrants, in order to contrast fire growth phases in road tunnels.
Therefore, further studies, based on analysis of uncertainties and intermodal transport, need to be completed.