Research Article  Open Access
Qi Zhang, Yanhui Zhuang, Yuguang Wei, Hao Jiang, Hao Yang, "Railway Safety Risk Assessment and Control Optimization Method Based on FTAFPN: A Case Study of Chinese HighSpeed Railway Station", Journal of Advanced Transportation, vol. 2020, Article ID 3158468, 11 pages, 2020. https://doi.org/10.1155/2020/3158468
Railway Safety Risk Assessment and Control Optimization Method Based on FTAFPN: A Case Study of Chinese HighSpeed Railway Station
Abstract
In order to make safety risk assessment more accurately and more reasonably for highspeed railway station in China, this paper analyzes risk factors of fault tree and transfers the fault tree of risk accident into fuzzy petri net and then builds the FPNFTA model by combining the dynamic weighting fuzzy petri net (FPN) and fault tree analysis (FTA) based on the latter. This paper simulates the FTAFPN model with Stateflow of Matlab software. Then, it builds up a biobjective risk control model, making the minimum safety risk level and minimum necessary cost as the objectives, and it designs discrete particle swarm optimization algorithm to solve the risk control model. Finally, this paper selects stampede accident of Shijiazhuang highspeed railway station as an example in case study for assessing stampede risk level and gets the risk control schemes for this station. The results verify the feasibility and validity of the model and algorithm.
1. Introduction
The pressure on railway transport safety management has increased constantly with the rapid development of railway. Therefore, safety risk control is the key to secure railway transport. Safety risk control refers to the scientific management methods of realizing the maximum safety guarantee with the minimum cost by controlling risks effectively and optimizing safety risk control technologies based on the assessment of safety risks during the enterprise production activities. Safety Risk Assessment and Control Optimization can not only reduce probability of risk occurrence but also keep the risk control cost within reasonable boundaries, and many scholars have done researches on this topic [1–14]. Fault tree theory is widely applied to safety risk assessment at present [5–8], while fuzzy petri net is mostly used for safety diagnosis [9–14]. However, there is no existed literature that combines the fault tree and fuzzy petri net to make railway safety risk assessmentexcept our master thesis [15].
Although these methods or models have contributed to the safety risk control theory of railway transport in different aspects, there are still several problems need to be discussed.①In the previous safety risk assessment research, the indicator system is used to analyze the system as a whole, but the causal correlation of the development of security risk events is not described.②The traditional studies on the safety risk of transportation system consider the ambiguity and uncertainty of risk factors not enough. It is impossible to determine the risk status of the system accurately.③The weight value of each risk factor is relatively fixed, and it cannot directly reflect the influence degree of the development stage with different risk factors of the risk accident.④In previous studies, the control of the risk is lack of balanced analysis between benefit and cost since there is not enough consideration for the cost control.
In order to solve these problems, this paper takes the similar characteristics of dynamic weighting fuzzy petri net (FPN) and fault tree analysis (FTA) into consideration, introduces FTAFPN into railway safety risk assessment, analyzes risk factors of fault tree, maps the fault tree of risk accident into fuzzy petri net, and builds up railway safety risk assessment model combined with both FPN and FTA [4]. Then, this paper proposes to simulate FTAFPN model by Stateflow, considering the graphic characteristics of fuzzy petri net. Based on Target Cost Method, cost control is introduced into risk control optimization. A biobjective model with the minimum safety risk level and minimum necessary cost as the objectives is built up, and corresponding discrete particle swarm algorithm is designed for solving the model. The paper makes the risk control of passenger traffic congestion in highspeed rail passenger station as an example to carry out method verification. An empirical study is conducted at the Shijiazhuang Highspeed Railway Station to assess the safety risk level of the crowded stampede event and the control optimization plan.
2. Safety Risk Assessment and Control Model of Railway Based on FTAFPN
2.1. General Idea
In railway system, different control methods are taken to deal with different risk factors, then different costs and contributes to system general risk level are caused. Therefore, systematic risk control aims at reducing the system risk level effectively by controlling different risk factors selectively and saving the cost. This paper builds railway safety risk assessment model combining FPN and FTA and builds the double objective optimization model of risk control with the minimum safety risk level and minimum necessary cost of safety risk control.
2.2. Theoretical Research
2.2.1. Fault Tree Analysis
Fault tree can be used for analyzing the relevant cause events, including their interactions and logical relations with a specific selected risk accident as the target, building up the fault tree model and analyzing the importance of cause events and finally designing the scientific safety risk control measures contrapuntally. Fault tree consists of two parts: qualitative analysis and quantitative analysis. Qualitative analysis refers to comprehensive recognition of risk factors, mapping logical graph, and finding all cause events, including their combinations. Quantitative analysis refers to analysis of the structural importance of the cause events and determination of their influences on the top event, that is, the risk accident.
2.2.2. Dynamic Weighting Fuzzy Petri Net
Dynamic weighting fuzzy petri net is a transformation of the normal petri net; it consists of places, transitions, and the flow relations between them. However, the composition and element meaning of dynamic weighting fuzzy petri net are different from the normal petri net. Dynamic weighting fuzzy petri net can be defined by a tenelement group [8]:where refers to the nonempty set of places which represent cause events. refers to the nonempty set of transitions which represent the happening of the events. refers to the nonempty set of propositions. , . refers to the entry incidence matrix which represents the entry relations from places to transitions . . refers to the input arc. If there is input arc between place and transition , then , otherwise . refers to the export incidence matrix which represents the export relations from transitions to places . . refers to the export arc. If there is output arc between transition and place , then , otherwise . refers to the real number mapping of transitions. is the certainty factor (CF) of the inference rule of transition . , . refers to the real number mapping of places which represent the validity of propositions that correspond to places . because the places are mapped into propositions, so places and propositions are matched one by one. is the token of the petri net. , . refers to the mapping which represents the threshold values of the triggering of transitions . . refers to the dynamic weighting which represents the support of prerequisite proposition to result proposition; it varies with the changing of . Dynamic weighting reflects the degree of associations between a certain indicator and other indicators in the evaluation system. The dynamic performance of the petri nets means the dynamic weight will change dynamically with the assessed value of each evaluation of the cause events in the petri nets. Formula (2) shows the value of . , : refers to the identification of dynamic weighting fuzzy petri net.
The fuzzy rules of fuzzy petri net are as follows:(1), , are prerequisite propositions and is result proposition. represent the different supports of prerequisite propositions to result proposition. If and , then the transition is triggered. The validity of prerequisite proposition is not changed, and the validity of result proposition is .(2) . are prerequisite propositions; are validities of different rules, is result proposition, and are the threshold values of the triggering of different transitions. If , transition will be triggered. The validity of prerequisite proposition is not changed, and the validity of is .
2.3. The Safety Risk Assessment Model of Railway System Based on FTAFPN
2.3.1. Railway Safety Risk Assessment Model Based on FPN and FTA
Based on the features of railway system and safety risk assessment, this paper combined the advantages of fault tree method and fuzzy petri net theory and proposed the railway safety risk assessment method based on FTAFPN, as shown in Figure 1.
The specific rules for mapping fault tree model into fuzzy petri net model are as follows:①The events of fault tree and the places of petri net can be transformed into each other. The basic events of fault tree are corresponding to initial places of petri net; the top event of fault tree is corresponding to target place of petri net.②The causal relations between events of fault tree and the transitions of petri net can be transformed into each other.③The logical gates mainly consist of AND gates, OR gates, NOT gates. The negativity of the event can be set as a new event, so there will be no NOT gate existing in the fault tree. The logical gates of fault tree and symbols of petri net can be transformed into each other, as shown in Figure 2.④The special logical gates such as voting gates, exclusiveOR gates, and exclusion gates can be transformed into AND gates, OR gates, or their combinations. Therefore, the fault tree can only consist of AND gates and OR gates.
2.3.2. The Solving Method
There are two methods used for solving fuzzy petri net. The search algorithm based on the net model and the inference method based on the incidence matrix of petri net. The second method is popular at present for its convenience. However, it is the lack of efficiency while applied to largescale net. The first method can be visualized, but it is complicated and difficult for solving manually. According to the disadvantages of the first method, this paper proposes to use Stateflow for simulating petri net model and finally solve the fuzzy petri net. Stateflow is a kind of design tool with strong visualization and integrated into Simulink of Matlab software. It is widely used for the simulation of complicated eventdriven system.
2.4. The Safety Risk Control Model of Railway System
2.4.1. Assumptions
(1)The superior limit of system risk level and the limit of total cost of risk control is both known.(2)Cost spent on one risk source would not change the risk levels of other risk factors.(3)The risk levels of all risk factors are known before the risk control.(4)Leave force majeure events out of account.
2.4.2. Parameters
—general risk level of the whole system —cost for controlling the risk source —maximum permissible general risk level —limit of total cost for risk control —amount of risk factors —amount of controlled risk factors —decision variable of risk source —risk calculation function based on the FTAFPN model —risk level of risk source after control —risk level of risk source before control
2.4.3. The Railway Safety Risk Control Model
The safety risk control model with the minimum risk level and minimum necessary cost as the objectives is built up as follows:
Formula (3) refers to risk level of railway system, and formula (4) refers to the necessary cost for safety risk control of railway system.
Constraints:
Formula (5) represents that the risk level of railway system cannot exceed the maximum permissible general risk level. Formula (6) represents that the cost of risk control cannot exceed the limit of total cost of risk control. Formula (7) represents that the amount of controlled risk factors cannot exceed the amount of risk factors. Formula (8) represents the decision variable of risk source , if the risk source is controlled, then it is 1, otherwise, it is 0. Formula (9) represents that risk level of risk source after control is lower than that before control.
2.5. Discrete Particle Swarm Optimization Algorithm
This paper chooses Binary Particle Swarm Optimization (BPSO), also known as Discrete Particle Swarm Optimization Algorithm, to solve the safety risk control model. BPSO belongs to Particle Swarm Optimization (PSO). It has several features as follows:
2.5.1. Update Rule
where is the speed of particle in dimension in iteration ; is the position coordinates of particle in dimension in iteration ; is inertia weight of the current speed during state adjustment; Rand() is random coefficients evenly distributed between 0 and 1; and is acceleration limit factor. is the optimal position of particles in the population up to the iteration. is the optimal position experienced by all particles in the population until the iteration.
Different from basic PSO, the coding of BPSO is binary, the code of each particle is 0 or 1. Therefore, the update rule of basic PSO is changed by introducing the fuzzy function, as formula (10) shows.
For basic PSO, , which is definite, is the superior limit of particle speed, while for BPSO, the particle speed can be represented in the form of probability, and refers to the probability range of particle speed. Now, define the probability of particle position as , then . The transformation probability of particle will be
2.5.2. Inertia Weight
Inertia weight has influence on algorithm search capability. This article proposed to use nonlinear adaptive strategy for calculating the inertia weight , as formula (14) shows:where is the number of iterations at present; is the maximum number of iterations; is the initial inertia weight of particle at beginning; is the final inertia weight of particle at ending; is the number of iterations at present during calculation, initial value is 0; and is the inertia weight at present during calculation, initial value is ,
Specific calculation process: Step 1: Initialize parameters of model, algorithm, and particle initial status such as position and speed. Calculate the sufficiency for each particle by the following formula , and . Step 2: Compare the sufficiency with the current best sufficiency Pbest of each particle; update Pbest if it is better. Step 3: Compare the best sufficiency Pbest of each particle with the current best sufficiency Gbest of the whole group; update Gbest if it is better. Step 4: Calculate inertia weights by formula (12), then update particle positions by formula (10), correct unsuitable particles of constraints, and generate the new particle group. Step 5: Check the number of iterations, if it has reached the maximum number, end calculation, and turn to Step 6, otherwise, turn to Step 1. Step 6: Output result.
3. Case Study
This paper uses stampede accident in highspeed railway stations as the example of method validation. Comparing with the normal railway station, highspeed railway station covers a larger area, has more equipment, and more complicated passenger flows. Stampede accident may easily occur if there is passenger tumble or escalator failure. Therefore, it is vital to control the risks of stampede accidents accurately and efficiently in the highspeed railway station.
3.1. Stampede Risk Assessment Model of HighSpeed Railway Station Based on FTAFPN
This paper builds the fault tree model of stampede accident of highspeed railway station and maps it into fuzzy petri net based on the transformation methods. Then, the stampede risk assessment model of highspeed railway station based on FTAFPN is built. Figure 3 and Table 1 specify the contents of places. Place represents factor that may cause stampede accident, the corresponding proposition represents the status of the factor. The reliability of represents its assessed value under different statuses. Transition represents the development of event. Fuzzy factor of transition represents the possibility of development. The corresponding risk level of token value of place can represent the stampede risk level in highspeed railway station.

3.2. Stampede Risk Control Model of HighSpeed Railway Station Based on FTAFPN
The stampede risk control model of highspeed railway station based on the FPN model of stampede risk is built, and its object function and constraints are shown as follows:
3.3. Assessment Criteria Quantization of HighSpeed Railway Station Based on Cloud Model
According to the characteristic of assessment of stampede, this paper classifies assessment criteria into 5 classifications, including Good (H), Fair (JH), Normal (Z), Substandard (JC), and Poor (C). It quantifies the assessment criteria based on cloud model. 5 cloud models that are can be generated with golden section method in the universe of . 5 cloud models are corresponding to 5 classifications, and numeric characteristics of each cloud model is shown in Table 2.

S experts, recorded as , are invited to assess stampede risk factors, and expert weights are . The assessment criteria of stampede accident are recorded as , and the assessed values are . All criteria are calculated based on formula (15). More than 99.73% clouds shall be situated in the range of according to the rule. Considering characteristics of Hyperentropy calculation, transform cloud model is used into interval number . The assessed value of interval number can be replaced by expectation of interval. In this paper, the cloud model is normal distribution. Therefore, the assessed values of criteria are .
3.4. Stampede Risk Assessment of Shijiazhuang HighSpeed Railway Station
Shijiazhuang highspeed railway station is a large transportation hub with multiple transport modes, including subway, highspeed railway, normalspeed railway, and road transport.
In this paper, 15 experts are invited to participate in the evaluation work. However, considering the daily safety assessment of the highspeed railway station, the number of experts participating in the evaluation is limited. In order to make the evaluation results be in line with the actual work, the authors selected 3 typical experts and their evaluation results for Shijiazhuang highspeed railway station. One of them is an expert on the field of passenger safety supervision for highspeed railway passenger stations with rich supervision experience. One is a senior management with high management experience in highspeed railway passenger transportation. The other is a professor of highspeed railway passenger safety management with rich theoretical knowledge. For the opinions of the other 12 experts, we have summarized their opinions with the 3 typical expert opinions, which are selected to a certain degree, and finally form an evaluation index system for the thesis. Taking into account the equality of experts from different fields, the paper sets their weights as .
After calculation, the assessed values of criteria can be received, as shown in Table 3.

Based on the relations between criteria and fuzzy petri net, Stateflow is used to simulate the fuzzy petri net model of stampede risk in Shijiazhuang highspeed railway station; The Stateflow model is shown in Figure 4. After substituting the assessed values of criteria, the risk levels of some risk factors can be received, as shown in Figure 5.
The risk status can be divided into 4 segments according to relevant Chinese national standards, including safe, alert, dangerous, and extreme dangerous, and the corresponding ranges of assessed value are shown in Table 4 based on the literature [4, 9]. Calculation results show that risk status of stampede accident is alert. Main causes include crowded station entrance, escalator, or passage and two escalator failures during the assessment.

After simulation, the result is shown in Figure 4. The red square frames represent that the risk levels of the corresponding risk factors have reached the alert level, which is 0.6, and the control methods are necessary.
The risk levels of risk factors of stampede accident in the Shijiazhuang highspeed railway station can be received, as shown in Figure 6. The risk statuses of risk factors P_{18}, P_{19}, P_{16}, P_{1}, P_{4}, and P_{8} are alert. The risk statuses of risk factors P_{2} and P_{3} are nearly dangerous.
The control methods for risk factors are shown in Table 5. Fire explosion P_{7} and terrorist attack P_{9} are left out of consideration for their assessed values are very low.

Different control methods need different amount of staff or materials, which leads to different costs. Therefore, this paper determines the costs of different control methods for different risk factors based on investigation and survey, which is shown in Table 6.

are the predicted risk levels after control, assessed by experts. The model is solved by discrete particle swarm optimization algorithm; initialization is necessary at the beginning. Population size , study factor , , , . Matlab is used for programming.
This paper sets 3 different cost constraints, including 700,000 Yuan/year, 900,000 Yuan/year, and 1150,000 Yuan/year for comparing different control results. The control results under different cost constraints are shown in Table 7.

Figure 7 shows the control results of different risk factors under different cost constraints. Although the cost is 900,000 Yuan/year, the result is obvious, which is 0.4143, and all risk factors remain in safe status. If safety redundancy is needed, that is, risk levels of all risk factors should remain under 0.5, and then, the control cost should be 1150,000 Yuan/year. This paper concludes that Scheme 2 is preferable under overall consideration.
4. Conclusions
This paper analyzes the relevant cause events of risk accidents, including their interactions and logical relations with fault tree model and then maps fault tree model into fuzzy petri net model. Risk control model with the minimum safety risk level and minimum necessary cost as the objectives is also built up with the risk control costs as constraints. Discrete particle swarms optimization algorithm is designed accordingly. Meanwhile, Stateflow of Matlab software is used for simulating fuzzy petri net and visualized simulation for risk assessment model based on FTAFPN is also achieved.
This paper selects stampede accident of the Shijiazhuang highspeed railway station as an example in case study for assessing stampede risk level and receiving risk control schemes. Results show that the risk level of Shijiazhuang highspeed railway station is 0.601 which is alert. Then, research for safety risk control was carried out, and results show that while the limit of total risk control cost is 900,000 Yuan/year, risk source P_{1}, P_{2}, P_{3}, P_{4}, and P_{8} should be controlled, and the risk level shall be reduced to 0.4143. The improvement is obvious and risk levels of all risk factors remain safe.
Data Availability
Major part of the data used to support the findings of this study are included within the article. Further information can be obtained from the author Mr. Zhuang.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is supported by the Chinese National Railway Administration Project, titled “Research on Safety Monitoring of Highspeed Railway Passenger Transport” (T17DJ00071), and 2016 China National Key Research and Strategic International Science and Technology Innovation Cooperation Project, titled “‘the Belt and Road’ Land Transport International Intermodal Research and Exchange Center” (2016YFE0201700).
References
 H. Luo, Z. Yu, and J. Ning, “Research on railway risk control method and its application,” Journal of Computer Applications, vol. 26, no. 4, pp. 40–47, 2017. View at: Google Scholar
 J. Zhang and L. Liu, “Research on highspeed railway station safety management based on risk control theory,” Journal of Transportation Engineering and Information, vol. 11, no. 4, pp. 99–104, 2013. View at: Google Scholar
 Z. An, X. Dai, and J. Liu, “Study on safety risk control of highspeed railway,” Railway Transport and Economy, vol. 36, no. 1, pp. 22–27, 2014. View at: Google Scholar
 Y. Zhuang, Study on Safety Risk Assessment and Control Optimization of Passenger Transport in HighSpeed Railway Passenger Station Based on Fuzzy Petri Net Theory, Beijing Jiaotong University, Beijing, China, 2018.
 J. Liu, “Synthesized risk assessment method for railway safety based on FTAAHP analysis,” China Railway Science, vol. 38, no. 2, pp. 138–144, 2017. View at: Google Scholar
 G. Zhan, X. Shang, and L. Hai, “AHPbased safety assessment model for rail transit system,” China Railway Science, vol. 32, no. 3, pp. 123–125, 2011. View at: Google Scholar
 J. Liu, X. Dai, G. Zhan et al., “Quantitative safety assessment method based on risk in railway system,” China Railway Science, vol. 30, no. 5, pp. 123–127, 2009. View at: Google Scholar
 E. Jafarian and M. A. Rezvani, “Application of fuzzy fault tree analysis for evaluation of railway safety risks: an evaluation of root causes for passenger train derailment,” Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, vol. 226, no. 1, pp. 14–25, 2012. View at: Publisher Site  Google Scholar
 P. Ghu, L. Liu, and J. Yin, “On the subway fire risk assessment based on the dynamically variable weight petri net,” Journal of Safety and Environment, vol. 16, no. 6, pp. 39–44, 2016. View at: Google Scholar
 P. Ghu and L. Liu, “Safety assessment for integrated transportation hub based on variable weight fuzzy petri net,” Journal of XIhua University (Natural Science), vol. 35, no. 5, pp. 13–17, 2016. View at: Google Scholar
 A. Giua and C. Seatzu, “Modeling and supervisory control of railway networks using petri nets,” IEEE Transactions on Automation Science and Engineering, vol. 5, no. 3, pp. 431–445, 2008. View at: Publisher Site  Google Scholar
 J. He, D. Zhuang, X. Zhang et al., “Study on safety of fire accident in urban mass transit system based on FTApetri net,” China Safety Science Journal, vol. 19, no. 10, pp. 77–177, 2009. View at: Google Scholar
 D. Zhang and J. He, “Constructing security model of vehicle in urban mass transit system based on petri net,” Journal of Transportation Engineering and Information, vol. 7, no. 4, pp. 109–113, 2009. View at: Google Scholar
 P. Lei, The Method Study of Fuzzy Petri Nets Based Accident Causation analysis, Guangdong University of Technology, Guangdong, China, 2015.
 J. Zhang, Research on Risk Management and Evaluation Model of Grid Enterprise Informatization Project, North China Electric Power University, Beijing, China, 2015.
Copyright
Copyright © 2020 Qi Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.