The Group Evacuation Behavior Based on Fire Effect in the Complicated Three-Dimensional Space

In order to effectively depict the group evacuation behavior in the complicated three-dimensional space, a novel pedestrian flow model is proposed with three-dimensional cellular automata. In this model the calculation methods of floor field and fire gain are elaborated at first, and the transition gain of target position at the next moment is defined. Then, in consideration of pedestrian intimacy and velocity change, the group evacuation strategy and evolution rules are given. Finally, the experiments were conducted with the simulation platform to study the relationships of evacuation time, pedestrian density, average system velocity, and smoke spreading velocity. The results had shown that large-scale group evacuation should be avoided, and in case of large pedestrian density, the shortest route of evacuation strategy would extend system evacuation time.


Introduction
As the density of urban population increases rapidly, the sudden disaster in the complicated environment may bring about serious consequences, so the study focus of the pedestrian flow is how to provide effective evacuation strategy [1][2][3][4][5].At present, the evacuation methods of pedestrian flow are mainly divided into continuous model and discrete model.The typical continuous model is social force, which was used by Helbing et al. [6] to study evacuation process of pedestrians, discovering the formation mechanism of phenomena such as "faster is slower" and "arching" effect.The typical discrete model is lattice gas and cellular automata.The pedestrian evacuation efficiency in environment such as bottleneck channel, T-type channel, and hall was studied by Tajima et al. [7][8][9] with lattice gas model, analyzing the critical threshold value with phase transformation under different conditions.The evacuation process of pedestrian flow was studied by Guo and Huang [10] with dynamic lattice gas model.Through considering the repulsive force in cellular automata, Kirchner et al. [11,12] studied the phenomena of pedestrians avoiding each other when competing for the same position.The pedestrian evacuation model of twodimensional cellular automata was established by Burstedde et al. [13], which proposed "floor field" concept to study typical phenomenon.Nagai et al. [14] studied the evacuation characteristics of pedestrians whose view is blocked completely.
Moreover, the pedestrians' intimacy often leads to group evacuation.Based on A * algorithm, the evacuation efficiency of different-scale group was studied by Ji and Gao [15] through establishing the group evacuation model with leader who determined the direction of group evacuation.Qiu and Hu [16] established evacuation model and depicted the group behavior for evacuation efficiency, and the simulation results showed that the congestion degree was influenced by group scale.By combining conformity mentality, Low [17] elaborated the phenomena of some group following another group in evacuation.Furthermore, the factors such as pedestrians' panic [18,19], game [20], interaction between pedestrians [21], and fire spreading [22] were also introduced to cellular automata model to study evacuation efficiency.
Meanwhile, pedestrians' velocity of walking upward or downward the stairway is clearly different from their velocity 2 Mathematical Problems in Engineering (a) Cellular structure of walking horizontally, which has great effect on the evacuation efficiency.Kretz et al. [23] found that the distribution of horizontal walking velocity exhibited a maximum between 0.4 m/s and 0.5 m/s in the local density.And some pedestrians actually accelerate when walking upward a short stairway; however, such a phenomenon could neither be observed in downward motion on the short stairway nor in upward motion on the long stairway.However, the presented cellular automata model mainly focuses on the two-dimensional platform, which is different from the actual three-dimensional evacuation space containing stairways.But stairway factor has a great influence on the evacuation velocity and evacuation efficiency.Therefore, in this thesis a novel cellular automata model is proposed for the three-dimensional space in the complicated environment, and the group behavior is studied with different evacuation strategies when fire breaks out.

Methodology
As shown in Figure 1, the two-dimensional Moore structure is extended to the three-dimensional structure, every cellular of which is a basic unit.There are three statuses for every cellular: occupied by pedestrian, occupied by obstacle, or vacant.It is assumed that the pedestrian is located at the cellular (, , ); in each time step, the transition gain   of the neighboring twenty-six positions and his own position is calculated by (1), making the smallest transition gain as the target position at the next moment: In this formula, , , and  are the position coordinates,   is the transition gain of cellular (, , );   is the vacancy coefficient; if cellular (, , ) has been occupied by other pedestrians, then   = 1; otherwise   = 0;   is the obstacle coefficient, if cellular (, , ) is the obstacle, then   = 0; otherwise   = 1;   and  indicate the floor field and its influence coefficient; the bigger   is, the further it is from the exit;   and  indicate the fire gain and its influence coefficient, the bigger   is, the more serious it is influenced by fire.And  +  = 1,0 ≤  ≤ 1,0 ≤  ≤ 1; the bigger  is, the more possible it is for pedestrians to adopt the shortest route strategy; the bigger  is, the more possible it is for pedestrians to consider avoiding being influenced by fire.
Meanwhile, different from the two-dimensional cellular automata Moore structure, the parameter   is added to indicate whether the cellular (, , ) exists or not.If cellular (, , ) does not exist, then   = 1; otherwise   = 0.And the parameter   is used to indicate the transition gains for walking upstairs, downstairs, or horizontally.Usually, the pedestrians tend to walk horizontally more, which saves physical power best, and then the second best is walking downward the stairways; the last is walking upward the stairways.Therefore, the parameter   is set as 5, for walking horizontally, 1, for walking downstairs, 0.5, for walking upstairs. (2) 2.1.The Floor Field.In the two-dimensional platform, the calculation methods of "floor field" have been put forward in many studies [11][12][13][24][25][26].In this thesis, the calculation method of floor field is proposed in three-dimensional space.Considering the complicated three-dimensional space as shown in Figure 2, there exists an exit (the width of which is five cellular), and the white areas stand for wall or obstacle.The floor field   of the cellular (, , ) is determined based on the following principles: (1) for the wall or obstacle, the floor field   = 500, indicating that the pedestrians cannot pass; (2) for the exit, the floor field   = 0; (3) the floor field of the straight adjacent cellular at the same level as the exit is added by 1, the floor field of the diagonal adjacent cellular at the same level as the exit is added by √ 2, the floor field of the straight adjacent cellular at the adjacent level as the exit is added by √ 2, and the floor field of the diagonal adjacent cellular at the adjacent level as the exit is added by √ 3; (4) the floor field of other cellular is determined according to (3).Based on the principles above, the floor field of threedimensional space is shown in Figure 3: if in the different layer 0, if (, , ) exists 500, if (, , ) is wall or obstacle. (3) 2.2.The Fire Gain.Usually, after the fire breaks out, the spreading velocity of smoke has major impact on evacuation.
For the sake of depicting the problems clearly, the influence of smoke on the cellular (, , ) is divided into three grades: no effect, general effect, and serious effect.If the cellular (, , ) is the fire source or influenced seriously, then the fire gain   = 500, indicating that pedestrians cannot pass; if the cellular (, , ) is not influenced by the smoke, then   = 0.The fire gain   is shown in for general effect, 500, for serious effect or fire source.(4) In this formula, ( > 0) is the fire coefficient.It is assumed that the cellular (  ,   ,   ) is the fire source, which spreads to the neighboring cellular after 2m time steps.The influence of the smoke exerted on the cellular (, , ) is set to be general effect initially; then it is changed to serious effect after 14m time steps.The smoke continues to spread until the whole evacuation space is evacuated.
According to the findings of [23], pedestrians' velocity of walking horizontally on the stairways V = 0.42 m/s, the average velocity of walking upstairs V = 0.71 m/s, and the average velocity of walking downstairs V = 0.65 m/s.But in reality, when influenced by the smoke, the pedestrians' walking velocity would be slower.The pedestrians' walking velocity under influence of the smoke is set as V, for no effect, V  , for general effect, 0, for serious effect or fire source. (5)

Group Behavior.
As influenced by intimacy and conformity mentality, under emergency and in the complicated environment, pedestrians tend to evacuate in herding.Therefore, in this thesis, the effect of group behavior on the evacuation efficiency is studied.Assume group  is formed; the location of group becomes the key factor of characterization.As the location information of group has some range, for the convenience of depiction, the average location of all members is defined as the location of group GLocation (, , ).But we should exclude the fact that some members are far from the group, for it makes the central location of the group not the actual center.The location coordinates of all members in the group are sorted based on axes , , and , obtaining three sorted arrays ( 1 ,  2 , . . .,  NUM ), ( 1 ,  2 , . . .,  NUM ), and ( 1 ,  2 , . . .,  NUM ), and NUM is the number of group.Then % (generally, 5 <  < 20) data in front and back of every array are dislodged, the group center GLocation (, , ) is calculated according to (6) with the remaining data num (num = NUM(1 − %)).This is to avoid the impact of some far away members exerted on the group central location: In this formula, ⌊ ⌋ indicates the rounding down.Meanwhile, the biggest distance between the location of member  Location (  ,   ,   ) and the group center GLocation (, , ) is defined as , if Then, the member  should immediately adjust the current position by moving to group center GLocation (, , ) at once.

Evolution Rule.
The model adopts parallel update mechanism, and the pedestrians meet the following rules from moment of  to moment of  + 1 in evacuation.
(1) At the initial moment, the groups are initialized based on intimacy; the members of group calculate the transition gain   of the alternative position at moment of  + 1 according to (1).Meanwhile, the smallest transition gain   is set as the target  position at moment of  + 1.If there exist many identical smallest transition gains, one position is randomly chosen with equal probability.
(2) According to ( 6) and (7), if the distance of some group member with the GLocation (, , ) exceeds maximum , the group member should immediately adjust the current position by moving to group center GLocation (, , ).
(3) When many pedestrians compete for one target position at the same time, one of the pedestrians is chosen randomly by the system with equal probability to move to the target position, and the other pedestrians stay still.
(4) The pedestrians would evacuate out of the system at moment of  + 1 when moving to the exit.
(5) When all pedestrians finish evacuating, the simulation is completed.Meanwhile, pedestrian density  is defined as the ratio of the number of pedestrians to the number of all cellular in evacuation space, average system velocity  is defined as the average velocity of all pedestrians, and system evacuation time  is defined as the time for all pedestrians to leave the evacuation.It is set as  = 0.3,  =  = 0.5,  = 3,  = 2,  = 15, and  = 2.The fire source is set on the left side; the smoke is produced when  = 30 s.In Figure 4, the evacuation under the different time step is shown.When  = 5 s, the pedestrians are distributed evenly in the evacuation space.After the fire breaks out, the pedestrians choose the middle part or the right side for evacuation to avoid the fire source or smoke.Moreover, the "arching" effect at the exit is consistent with the reality

Results and Discussion
. By combining the above initial state and experimental results, the discussion of performance is made.In order to reduce the impact of the initial experiment status or waiting probability on the statistical indicators, the average value of thirty runtime value is taken for every indicator.In Figure 5, the relationship between system evacuation time  and the number of group members NUM under the different parameter  is given.From Figure 5 it can be seen that as NUM increases, system evacuation time  tends to increase as well; that is to say, the larger the scale of group, the more harmful it is to improve the evacuation efficiency.When NUM is small, the smaller the parameter , the shorter the system evacuation time.But when NUM is large, the curve suddenly changes: the bigger the parameter , the shorter the system evacuation time.This indicates that when other factors are unchanged, the small group adopting the shortest route strategy is good for improving evacuation efficiency, while large-scale group who choose to keep away from the fire is more beneficial for improving evacuation efficiency.It reflects that when the group scale is constant, the best evacuation efficiency may be achieved with the balanced evacuation method.Furthermore, in Figure 6 the relationship between system evacuation time  and pedestrian density  under the different parameter  is given.Similar to Figure 5, there occurs sudden change of curve in Figure 6.When pedestrian density is small, the smaller the parameter , the shorter the system evacuation time.But when pedestrian density is large, the bigger the parameter , the smaller the evacuation time.In case of large pedestrian density, if the shortest route strategy is preferred for evacuation, the congestion may be produced as the smoke would block some part of evacuation routes, which extends system evacuation time.
Meanwhile, the relationship between system evacuation time  and average system velocity  is shown in Figure 7. From Figure 7, as average system velocity  increases, system evacuation time  tends to decrease firstly and increase secondly, and the inflection of curve emerges when  = 1.8.This indicates that it cannot reduce the evacuation time by improving the evacuation velocity only.As the whole evacuation velocity of the pedestrians improves, serious congestion may be created at the exit, which is consistent with the "faster is slower" effect [19].
Finally, the relationship between system evacuation time  and the time step of smoke spreading velocity  under the different parameter  is given in Figure 8.The smaller the , the faster the smoke spreading.From the figure it can be seen that as  increases, system evacuation time  tends to decrease.The smaller the impact of smoke exerted on the evacuation space, the higher the evacuation efficiency.Meanwhile, when  is smaller, system evacuation time corresponding to the shortest route strategy (such as  = 0.1) shows longer evacuation time.Under such condition, that is, the smoke spreading velocity is big, it is bad for reducing system evacuation time to adopt the shortest route strategy.When  is big, the smoke spreading velocity is small; it can improve the evacuation efficiency to adopt the shortest route strategy (such as  = 0.1).

Conclusions
In order to depict the group evacuation behavior and evacuation efficiency when fire breaks out in the complicated threedimensional environment, a novel evacuation model is established based on the three-dimensional cellular automata.The calculation methods of floor field and fire gain are given in the model, in which the transition gain of the target position at the next moment is defined, and pedestrians will move to the position which is the minimum gain.Meanwhile, based on pedestrians' intimacy and velocity change, the group  evacuation strategy is elaborated.At last, the established experiment platform is used for simulation; the results show that (1) in case of large pedestrian density, the shortest route strategy would extend the system evacuation time; (2) the larger scale the group, the lower the evacuation efficiency; (3) it cannot effectively reduce the evacuation time by improving the evacuation velocity only.In our future studies, the panic mentality and the interaction of pedestrians would be taken into account in the three-dimensional automata model.to thank the reviewers for their valuable comments on the paper.

3. 1 .
Simulation.Based on the three-dimensional automata model above, the experiment platform is established to study the performance of group evacuation when the fire breaks out.The complicated three-dimensional evacuation space is considered as shown in Figure2.The volume of evacuation space is set as 14 m × 14 m × 1 m (1 m is the height from the first floor to the last floor), and there are five floors and (a)  = 5 (b)  = 58

Figure 5 :
Figure 5: The relationship between system evacuation time  and the number of group members NUM.

Figure 6 :
Figure 6: The relationship between system evacuation time  and pedestrian density .

Figure 7 :
Figure7: The relationship between system evacuation time  and average system velocity .

Figure 8 :
Figure 8: The relationship between system evacuation time  and the time step of smoke spreading .