# Study on Evaluation Model of Emergency Rescue Capability of Chemical Accidents Based on PCA-BP

Authors: Jianghong Liu; Junfeng Wu; Weisi Liu
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

The emergency management of chemical accidents plays an important role in preventing the expansion of chemical accidents. In recent years, the evaluation and research of emergency management of chemical accidents has attracted the attention of many scholars. However, as an important part of emergency management, the professional rescue team of chemicals has few evaluation models for their capabilities. In this study, an emergency rescue capability assessment model based on the PCA-BP neural network is proposed. Firstly, the construction status of 11 emergency rescue teams for chemical accidents in Shanghai is analyzed, and an index system for evaluating the capabilities of emergency rescue teams for chemicals is established. Secondly, the principal component analysis (PCA) is used to perform dimension reduction and indicators’ weight acquisition on the original index system to achieve an effective evaluation of the capabilities of 11 rescue teams. Finally, the indicators after dimensionality reduction are used as the input neurons of the backpropagation (BP) neural network, the characteristic data of eight rescue teams are used as the training set, and the comprehensive scores of three rescue teams are used for verifying the generalization ability of the evaluation model. The result shows that the proposed evaluation model based on the PCA-BP neural network can effectively evaluate the rescue capability of the emergency rescue teams for chemical accidents and provide a new idea for emergency rescue capability assessment.


---

## Body

## 1. Introduction

Due to the properties of hazardous chemicals, such as toxicity, corrosiveness, explosiveness, flammability, and combustion support, there are huge risks in their production, transportation, storage, sales, use, and disposal. Once a hazardous chemical accident occurs, it is easy to cause many casualties, huge property losses, and serious environmental pollution and bring catastrophic consequences to both enterprises and the society. For example, the explosion of a dangerous good warehouse in Tianjin Port on August 15, 2015, resulted in 165 deaths, 798 casualties, and 8 missing. The direct economic loss reached 6.866 billion yuan [1]. Therefore, the emergency treatment of chemical accidents must be timely and efficient to prevent accidents from expanding and causing even greater losses.In recent years, the assessment of emergency management capabilities for hazardous chemical accidents has attracted the attention of many scholars. Wang et al. [2] proposed a disaster management control capability assessment model based on the Capability Maturity Model (CMM). This model evaluates the capability of the organization from eight aspects and divides the capability assessment results into four levels, which provides general assessment guidelines for different types of emergency management organizations. Lin [3] analyzed the nature of emergency rescue capabilities from the perspective of the entire city and established an urban emergency rescue capability evaluation system based on AHP and Fuzzy Comprehensive Evaluation (FCE). Yang et al. [4] analyzed many factors that affect the emergency capacity of enterprises, established an emergency capacity evaluation index system, and determined the weight for each indicator through the Analytic Hierarchy Process (AHP). Yu and Guan [5] analyzed the current situation and difficulties of emergency treatment of hazardous chemical accidents, discussed the emergency training system, and provided a reference for improving the emergency capabilities of professional teams for fire and hazardous chemical rescue. Zhu et al. [6] used Bayesian networks to propose a framework for dynamically evaluating explosion accidents in chemical plants to support prevention, management, and real-time warning. He et al. [7] established a Petri net model of emergency process of chemical accidents in order to evaluate the emergency capability, which can dynamically evaluate the emergency capability of chemical accidents.In addition to the abovementioned traditional evaluation methods, the application of artificial neural networks to the evaluation of chemical accidents has made some progress. Yuan et al. [8] used back propagation neural networks, generalized regression neural networks and radial basis function neural networks to evaluate the safety production management of chemical companies and found that the prediction ability of the radial basis function neural network is more accurate. Aiming at the shortcomings of the current chemical production safety evaluation system and combining the knowledge of artificial neural networks, Yang established a new evaluation index system and proposed the advantages of applying neural networks to the chemical production safety evaluation system [9].These studies mainly focus on the establishment of an emergency management assessment system for hazardous chemical accidents and the application of emergency assessment methods. Or, consider the emergency management activity itself as a project management process and study the capacity assessment model for emergency management control. However, as an important part of emergency management, the professional rescue team of hazardous chemicals has few evaluation models for their capabilities. In addition, traditional evaluation methods, such as Analytic Hierarchy Process (AHP), are greatly affected by human factors in the implementation process, and it is difficult to obtain objective evaluation results. When there are many evaluation indicators, it will complicate the structure of the artificial neural network model and increase the computational complexity. Shanghai is an important petrochemical and fine chemical industry base in China, with a solid chemical industry foundation. Through the assessment of the emergency response capabilities of the 11 professional rescue teams for hazardous chemicals in Shanghai, the capabilities of the rescue teams can be strengthened using targeted countermeasures.In order to reasonably evaluate the capabilities of professional emergency rescue teams for hazardous chemical accidents, this study surveyed 11 professional rescue teams in Shanghai, analyzed the status of these teams, constructed a rescue capacity assessment index system, and built a rescue capability evaluation model combined with BP neural network. At the same time, in order to determine the indicators’ weight and reduce the number of neurons in the input layer of the backpropagation (BP) neural network, the principal component analysis (PCA) was used to reduce the dimension of the evaluation index system and obtain weight. The dimensionality-reduced feature factors were used as the input units of the BP neural network. This method can not only reduce the influence of human factors in the evaluation process but also simplify the structure of the artificial neural network and reduce the computational complexity of the evaluation model. The trained BP neural network evaluation model can well evaluate the capabilities of professional emergency rescue teams for hazardous chemical accidents, providing a new idea for emergency rescue capability assessment.


## 2. Methods

### 2.1. Construction of Rescue Capability Evaluation Index System

Shanghai has a total of 11 emergency rescue teams for production safety, as shown in Table1. At present, the 11 emergency rescue teams for production safety are managed by the company where they work. The Shanghai Emergency Management Bureau is responsible for providing business guidance. As a result of a thorough investigation, following construction problems were found with these teams:(1)
Inefficient cooperation mechanism: there is a lack of coordination between the emergency rescue team and other departments. The team’s responsibilities are unclear, and there is no unified command.(2)
Slow emergency response: the lack of classification and corresponding response plans based on the type and scale of hazardous chemical accidents makes the emergency rescue scene more chaotic.(3)
Inappropriate team building: the positioning of these rescue teams is unclear, and the rescue areas are not divided. There are no long-term full-time members in these rescue teams, and these team members have not received any specialized training in dealing with hazardous chemical accidents.(4)
Insufficient emergency equipment, materials, and maintenance funds. The maintenance of professional rescue equipment for hazardous chemical accidents lacks government financial support. The necessary equipment and materials cannot be timely supplemented.(5)
Noncompliant emergency handling. There is a dearth of norms and standards on the emergency handling of hazardous chemical accidents. There is no targeted emergency response plan for different hazardous chemical accidents. Due to the problem of team building, rescuers with nonprofessional characteristics, and temporary combination, it is difficult to conduct emergency response scientifically and quickly.Table 1
The 11 emergency rescue teams for hazardous chemicals in Shanghai.
Serial numberTeam nameTeam 1Emergency Rescue Team of Shanghai Jiemeng Chemical Co., Ltd.Team 2Chemical Accident Emergency Rescue Center of Shanghai Institute of Occupational Disease for Chemical IndustryTeam 3Emergency Rescue Team of Testing Center of Shanghai Research Institute of Chemical Industry Co., Ltd.Team 4Emergency Rescue Team of Shanghai Huayi Energy Chemical Co., Ltd.Team 5Fire Brigade of BayerTeam 6Emergency Rescue Team of Shanghai Gaoqiao Petrochemical Chemical Transportation Co., Ltd.Team 7Chemical Rescue Team of BASFTeam 8Rescue Team of Shanghai Zhongshi Chemical Logistics Co., Ltd.Team 9Emergency Rescue Team of Shanghai Chlor-Alkali Chemical Co., Ltd.Team 10Fire Brigade of Refining Department of Sinopec Shanghai Gaoqiao Petrochemical Co., Ltd.Team 11Ambulance Team of Jinshan Hospital of Fudan UniversityAccording to the above construction status, an emergency rescue team assessment index system was established. The emergency response capabilities of professional rescue teams for hazardous chemical accidents include the following elements:(1)
Emergency cooperation: the main consideration is whether the division of responsibilities within the rescue team is reasonable and clear and whether smooth information can be communicated between various departments; whether the management of human resources has considered a perfect reward and punishment system, employee benefits and incentives; and whether the communication is efficient enough to ensure the normal operation of the emergency mechanism.(2)
Emergency command: the emergency rescue of hazardous chemical accidents mainly includes two aspects: an emergency disposal plan and the emergency professional and technical personnel. Different schemes are needed to respond to different chemical accidents as the quantity and types of hazardous chemicals always vary between different plants and regions. Therefore, certain requirements are put forward for the pertinence and completeness of the emergency response plan and the allocation of emergency professional and technical personnel.(3)
Emergency foundation: personnel, materials, equipment, and funds are the basis for emergency rescue of hazardous chemical accidents. In this study, the factors that affect basic emergency support are divided into four parts: (1) the emergency team, considering whether the stability, quantity, and quality of emergency personnel; (2) emergency equipment (including personnel protection equipment), the functionality, safety, quality, and quantity of equipment should satisfy the emergency disposal requirements; (3) whether emergency materials could meet different types of hazardous chemical accidents; (4), emergency funding, whether governments and enterprises had been given economic support to ensure better operation of the emergency rescue teams.(4)
Training and education: consider the training of professional emergency rescue knowledge and skills. Assess the improvement of the emergency rescue ability of the corresponding emergency personnel.(5)
Emergency drills: consider the workload, such as whether the number and time-frequency of drills is reasonable to meet the demands. In addition, the factors that need to be considered are the effects of the emergency drills, whether the personnel is familiar with the emergency procedures and more effective in responding to special chemical accidents through emergency drills.Accordingly, this research proposed a rescue capability evaluation index system that includes 5 first-level indicators, 14 second-level indicators, and 28 third-level indicators, as shown in Table2:Table 2
Index system for evaluating the capabilities of emergency rescue teams for hazardous chemical accidents.
First-level indicatorsSecond-level indicatorsThird-level indicatorsEmergency cooperation B1Duties within the team C1Post setting D1Clarity of post responsibilities D2Duty fulfillment D3Information exchange between departments C2Smooth communication D4Timeliness of information transfer D5Integrity of information transfer D6Human resources C3Reward and punishment system D7Salary and benefits D8Career opportunities D9Communication network C4Basic intercom equipment D10Advanced communication system and platform D11Emergency command B2Emergency plan C5Whether the emergency response plan is scientific and accurate D12Whether the emergency response plan is targeted D13Whether the emergency treatment plan is operable D14Emergency expert C6Whether emergency experts are qualified for emergency work D15Emergency foundation B3Emergency team C7Whether there is sufficient staff D16Emergency equipment C8Completeness of emergency rescue equipment D17Completeness of personal protective equipment D18The ability of emergency rescue equipment to control and reduce the impact of an accident D19Protective capabilities of personal protective equipment D20Emergency materials C9Completeness of emergency materials D21Replenishment of emergency materials D22Emergency funding C10Emergency funding support in Shanghai D23Territory emergency funding support D24Training and education B4Emergency rescue expertise training C11Emergency rescue expertise training D25Emergency rescue professional skills training C12Emergency rescue professional skills training D26Emergency drills B5Workload of emergency drill C13Workload of emergency drill D27Effect of emergency drill C14Effect of emergency drill D28


### 2.2. Dimension Reduction of Evaluation Indicators and Weight Acquisition

The above index system is too complicated; using the original indicator as the input unit of the BP neural network will face problems of such as high data dimensions, poor fitting effects, and inaccurate prediction results. Therefore, the principal component analysis (PCA) was required to reduce the dimensions of the indicators to eliminate the correlation [10].Principal component analysis (PCA) is an important statistical method that uses the idea of dimensionality reduction to transform multiple indicators into a few comprehensive indicators. These comprehensive indicators are not explanatory but retain most of the original information [11]. The new comprehensive indicator is a linear combination of all the original indicators which remain independent of each other. The principal component analysis (PCA) can reduce the number of evaluation indicators, thereby reducing the number of neurons in the input layer to simplify the structure of the BP neural network.In a geometric sense, the principal component analysis (PCA) method is to project the original data onto a new coordinate axis, which is the principal component. In order to enable the principal component to contain more information about the original data, the variance of the principal component must be maximized. The process of principal component analysis (PCA) is to find linear combination coefficients. The coefficients must maximize the variance of the principal components, and the sum of the squares of the coefficients must be equal to one. In addition, starting from the second principal component, each principal component must be independent of the existing principal components.The specific steps of principal component analysis (PCA) are as follows:Step 1: standardize the raw data.Assuming the original data is ann × m matrix:(1)X=x11x12⋯x1mx21x22⋯x2m⋮⋮⋮⋮xn1xn2⋯xnm.The rows in the data matrix represent different samples, and the columns represent different evaluation indicators. It can be seen that there aren samples and m evaluation indicators. The matrix can be normalized by the following formula:(2)xij∗=xij−xj¯varxj,i=1,2,…,n;j=1,2,…,m,among them,(3)xj¯=1n∑i=1nxij.Variance:(4)varxj=1n−1∑i=1nxij−xj¯2,j=1,2,…,m.Step 2: calculate the correlation coefficient matrix of the sample indicators:(5)R=r11r12⋯r1mr21r22⋯r2m⋮⋮⋮⋮rm1rm2⋯rmm,among them,(6)rij=covxi,xjvarxivarxj=∑k=1nxki−xi¯xkj−xj¯∑k=1nxki−xi¯2∑k=1nxkj−xj¯2,n>1.Step 3: calculate the eigenvalueλii=1,2,…,m and eigenvector vi=vi1,vi2,vi3,…,vim,i=1,2,…,m of the correlation coefficient matrix.Step 4: select the principal componentPt.All the eigenvaluesλii=1,2,…,m are arranged in descending order. The larger the eigenvalue, the more the system information contained in the principal component. Calculate the contribution rate of each principal component by the following formula:(7)ωi=λi∑i=1mλi,i=1,2,…,m.It is generally considered that the cumulative contribution rate of the firstq principal components exceeds 85% is reasonable, indicating that the total amount of system information they contain exceeds 85%. At this time, the principal component Ptt=1,2,…,q is the characteristic index after dimensionality reduction.Step 5: find the unit orthogonal feature vectorvt=vt1,vt2,vt3,…,vtm,t=1,2,…,q of the first q feature vectors. Each principal component is a linear combination of all the original indicators, and the coefficient atj is the element of the unit orthogonal eigenvector:(8)atj=vtj,t=1,2,…,q;j=1,2,…,m.Then, the expression of the principal componentPt of the n-th sample is(9)P1=a11xn1+a12xn2+a13xn3+⋯+a1mxnmP2=a21xn1+a22xn2+a23xn3+⋯+a2mxnm⋮Pq=aq1xn1+aq2xn2+aq3xn3+⋯+aqmxnm.Step 6: the comprehensive evaluation functionFn of the n-th sample is shown in the following formula:(10)Fn=ω1P1+ω2P2+⋯+ωqPq.The weight of the principal component is the contribution rateωi,i=1,2,…,q.


### 2.3. Construction of Evaluation Model Based on BP Neural Network

BP neural network is a multilayer feedforward network model. Its network is mainly composed of three parts: input layer, hidden layer, and output layer. As shown in Figure1, it can map t-dimensional data to l-dimensional data. The neurons in each layer of the BP neural network are not connected, and the output of the neurons in each layer only affects the output of the next layer. At the same time, the network will backpropagate errors during operation to continuously adjust the weights and thresholds of the network to achieve self-adjustment [12]. BP neural network is a nonlinear adaptive system, so it is more suitable for dealing with fuzzy or nonlinear problems. This method can effectively reduce subjective factors in the evaluation process and reduce the evaluation time [13].Figure 1
Structure of BP neural network.The main parameter settings of the BP neural network are the number of network layers, the number of nodes in the hidden layer, the transfer function, and the training function, which will be introduced one by one as follows.


#### 2.3.1. Number of Network Layers

As the number of network layers increases, the structure of the BP neural network will become more and more complex. Correspondingly, the complex BP neural network will prolong the learning time and cause the phenomenon of “overfitting.” Through previous testing and research, the neural network is generally set up as a three-layer network, that is, the input layer, the hidden layer, and the output layer [14].


#### 2.3.2. Number of Hidden Layer Nodes

Too few hidden layer neurons may not train the desired network, or the trained network is not strong enough and has poor generalization ability. To contrary, too many hidden layer neurons will increase the learning time and the error may not be smaller. The number of neurons in the hidden layer can be obtained by the following empirical formula [15]:(11)b=t+o+a,where b is the number of neurons in the hidden layer, t is the number of neurons in the input layer, o is the number of neurons in the output layer, and a is a constant between [1, 10].


#### 2.3.3. Transfer Function

There are three main transfer functions: logsig function, tansig function, and purelin function. The choice of the transfer function of the hidden layer and the output layer has a greater impact on the prediction accuracy of the BP neural network. The logsig transfer function is an S-type logarithmic function, the tansig function is an S-type hyperbolic tangent function, and both are nonlinear functions. After that, the purelin function is a linear function. Generally, the transfer function of the hidden layer selects the logsig function or the tansig function, and the transfer function of the output layer selects the purelin function. As far as the nonlinear transfer function is concerned, if the output of the samples is greater than zero, the logsig function is mostly used; otherwise, the tansig function is used.


#### 2.3.4. Training Function

The common training functions are as follows:trainlm: LevenbergMarquardt methodtraingd: gradient descent methodtraingdm: gradient descent method with momentum factortraingda: gradient descent method with adaptive learning ratetraingdx: gradient descent method with adaptive learning rate and momentum factor


## 2.1. Construction of Rescue Capability Evaluation Index System

Shanghai has a total of 11 emergency rescue teams for production safety, as shown in Table1. At present, the 11 emergency rescue teams for production safety are managed by the company where they work. The Shanghai Emergency Management Bureau is responsible for providing business guidance. As a result of a thorough investigation, following construction problems were found with these teams:(1)
Inefficient cooperation mechanism: there is a lack of coordination between the emergency rescue team and other departments. The team’s responsibilities are unclear, and there is no unified command.(2)
Slow emergency response: the lack of classification and corresponding response plans based on the type and scale of hazardous chemical accidents makes the emergency rescue scene more chaotic.(3)
Inappropriate team building: the positioning of these rescue teams is unclear, and the rescue areas are not divided. There are no long-term full-time members in these rescue teams, and these team members have not received any specialized training in dealing with hazardous chemical accidents.(4)
Insufficient emergency equipment, materials, and maintenance funds. The maintenance of professional rescue equipment for hazardous chemical accidents lacks government financial support. The necessary equipment and materials cannot be timely supplemented.(5)
Noncompliant emergency handling. There is a dearth of norms and standards on the emergency handling of hazardous chemical accidents. There is no targeted emergency response plan for different hazardous chemical accidents. Due to the problem of team building, rescuers with nonprofessional characteristics, and temporary combination, it is difficult to conduct emergency response scientifically and quickly.Table 1
The 11 emergency rescue teams for hazardous chemicals in Shanghai.
Serial numberTeam nameTeam 1Emergency Rescue Team of Shanghai Jiemeng Chemical Co., Ltd.Team 2Chemical Accident Emergency Rescue Center of Shanghai Institute of Occupational Disease for Chemical IndustryTeam 3Emergency Rescue Team of Testing Center of Shanghai Research Institute of Chemical Industry Co., Ltd.Team 4Emergency Rescue Team of Shanghai Huayi Energy Chemical Co., Ltd.Team 5Fire Brigade of BayerTeam 6Emergency Rescue Team of Shanghai Gaoqiao Petrochemical Chemical Transportation Co., Ltd.Team 7Chemical Rescue Team of BASFTeam 8Rescue Team of Shanghai Zhongshi Chemical Logistics Co., Ltd.Team 9Emergency Rescue Team of Shanghai Chlor-Alkali Chemical Co., Ltd.Team 10Fire Brigade of Refining Department of Sinopec Shanghai Gaoqiao Petrochemical Co., Ltd.Team 11Ambulance Team of Jinshan Hospital of Fudan UniversityAccording to the above construction status, an emergency rescue team assessment index system was established. The emergency response capabilities of professional rescue teams for hazardous chemical accidents include the following elements:(1)
Emergency cooperation: the main consideration is whether the division of responsibilities within the rescue team is reasonable and clear and whether smooth information can be communicated between various departments; whether the management of human resources has considered a perfect reward and punishment system, employee benefits and incentives; and whether the communication is efficient enough to ensure the normal operation of the emergency mechanism.(2)
Emergency command: the emergency rescue of hazardous chemical accidents mainly includes two aspects: an emergency disposal plan and the emergency professional and technical personnel. Different schemes are needed to respond to different chemical accidents as the quantity and types of hazardous chemicals always vary between different plants and regions. Therefore, certain requirements are put forward for the pertinence and completeness of the emergency response plan and the allocation of emergency professional and technical personnel.(3)
Emergency foundation: personnel, materials, equipment, and funds are the basis for emergency rescue of hazardous chemical accidents. In this study, the factors that affect basic emergency support are divided into four parts: (1) the emergency team, considering whether the stability, quantity, and quality of emergency personnel; (2) emergency equipment (including personnel protection equipment), the functionality, safety, quality, and quantity of equipment should satisfy the emergency disposal requirements; (3) whether emergency materials could meet different types of hazardous chemical accidents; (4), emergency funding, whether governments and enterprises had been given economic support to ensure better operation of the emergency rescue teams.(4)
Training and education: consider the training of professional emergency rescue knowledge and skills. Assess the improvement of the emergency rescue ability of the corresponding emergency personnel.(5)
Emergency drills: consider the workload, such as whether the number and time-frequency of drills is reasonable to meet the demands. In addition, the factors that need to be considered are the effects of the emergency drills, whether the personnel is familiar with the emergency procedures and more effective in responding to special chemical accidents through emergency drills.Accordingly, this research proposed a rescue capability evaluation index system that includes 5 first-level indicators, 14 second-level indicators, and 28 third-level indicators, as shown in Table2:Table 2
Index system for evaluating the capabilities of emergency rescue teams for hazardous chemical accidents.
First-level indicatorsSecond-level indicatorsThird-level indicatorsEmergency cooperation B1Duties within the team C1Post setting D1Clarity of post responsibilities D2Duty fulfillment D3Information exchange between departments C2Smooth communication D4Timeliness of information transfer D5Integrity of information transfer D6Human resources C3Reward and punishment system D7Salary and benefits D8Career opportunities D9Communication network C4Basic intercom equipment D10Advanced communication system and platform D11Emergency command B2Emergency plan C5Whether the emergency response plan is scientific and accurate D12Whether the emergency response plan is targeted D13Whether the emergency treatment plan is operable D14Emergency expert C6Whether emergency experts are qualified for emergency work D15Emergency foundation B3Emergency team C7Whether there is sufficient staff D16Emergency equipment C8Completeness of emergency rescue equipment D17Completeness of personal protective equipment D18The ability of emergency rescue equipment to control and reduce the impact of an accident D19Protective capabilities of personal protective equipment D20Emergency materials C9Completeness of emergency materials D21Replenishment of emergency materials D22Emergency funding C10Emergency funding support in Shanghai D23Territory emergency funding support D24Training and education B4Emergency rescue expertise training C11Emergency rescue expertise training D25Emergency rescue professional skills training C12Emergency rescue professional skills training D26Emergency drills B5Workload of emergency drill C13Workload of emergency drill D27Effect of emergency drill C14Effect of emergency drill D28


## 2.2. Dimension Reduction of Evaluation Indicators and Weight Acquisition

The above index system is too complicated; using the original indicator as the input unit of the BP neural network will face problems of such as high data dimensions, poor fitting effects, and inaccurate prediction results. Therefore, the principal component analysis (PCA) was required to reduce the dimensions of the indicators to eliminate the correlation [10].Principal component analysis (PCA) is an important statistical method that uses the idea of dimensionality reduction to transform multiple indicators into a few comprehensive indicators. These comprehensive indicators are not explanatory but retain most of the original information [11]. The new comprehensive indicator is a linear combination of all the original indicators which remain independent of each other. The principal component analysis (PCA) can reduce the number of evaluation indicators, thereby reducing the number of neurons in the input layer to simplify the structure of the BP neural network.In a geometric sense, the principal component analysis (PCA) method is to project the original data onto a new coordinate axis, which is the principal component. In order to enable the principal component to contain more information about the original data, the variance of the principal component must be maximized. The process of principal component analysis (PCA) is to find linear combination coefficients. The coefficients must maximize the variance of the principal components, and the sum of the squares of the coefficients must be equal to one. In addition, starting from the second principal component, each principal component must be independent of the existing principal components.The specific steps of principal component analysis (PCA) are as follows:Step 1: standardize the raw data.Assuming the original data is ann × m matrix:(1)X=x11x12⋯x1mx21x22⋯x2m⋮⋮⋮⋮xn1xn2⋯xnm.The rows in the data matrix represent different samples, and the columns represent different evaluation indicators. It can be seen that there aren samples and m evaluation indicators. The matrix can be normalized by the following formula:(2)xij∗=xij−xj¯varxj,i=1,2,…,n;j=1,2,…,m,among them,(3)xj¯=1n∑i=1nxij.Variance:(4)varxj=1n−1∑i=1nxij−xj¯2,j=1,2,…,m.Step 2: calculate the correlation coefficient matrix of the sample indicators:(5)R=r11r12⋯r1mr21r22⋯r2m⋮⋮⋮⋮rm1rm2⋯rmm,among them,(6)rij=covxi,xjvarxivarxj=∑k=1nxki−xi¯xkj−xj¯∑k=1nxki−xi¯2∑k=1nxkj−xj¯2,n>1.Step 3: calculate the eigenvalueλii=1,2,…,m and eigenvector vi=vi1,vi2,vi3,…,vim,i=1,2,…,m of the correlation coefficient matrix.Step 4: select the principal componentPt.All the eigenvaluesλii=1,2,…,m are arranged in descending order. The larger the eigenvalue, the more the system information contained in the principal component. Calculate the contribution rate of each principal component by the following formula:(7)ωi=λi∑i=1mλi,i=1,2,…,m.It is generally considered that the cumulative contribution rate of the firstq principal components exceeds 85% is reasonable, indicating that the total amount of system information they contain exceeds 85%. At this time, the principal component Ptt=1,2,…,q is the characteristic index after dimensionality reduction.Step 5: find the unit orthogonal feature vectorvt=vt1,vt2,vt3,…,vtm,t=1,2,…,q of the first q feature vectors. Each principal component is a linear combination of all the original indicators, and the coefficient atj is the element of the unit orthogonal eigenvector:(8)atj=vtj,t=1,2,…,q;j=1,2,…,m.Then, the expression of the principal componentPt of the n-th sample is(9)P1=a11xn1+a12xn2+a13xn3+⋯+a1mxnmP2=a21xn1+a22xn2+a23xn3+⋯+a2mxnm⋮Pq=aq1xn1+aq2xn2+aq3xn3+⋯+aqmxnm.Step 6: the comprehensive evaluation functionFn of the n-th sample is shown in the following formula:(10)Fn=ω1P1+ω2P2+⋯+ωqPq.The weight of the principal component is the contribution rateωi,i=1,2,…,q.


## 2.3. Construction of Evaluation Model Based on BP Neural Network

BP neural network is a multilayer feedforward network model. Its network is mainly composed of three parts: input layer, hidden layer, and output layer. As shown in Figure1, it can map t-dimensional data to l-dimensional data. The neurons in each layer of the BP neural network are not connected, and the output of the neurons in each layer only affects the output of the next layer. At the same time, the network will backpropagate errors during operation to continuously adjust the weights and thresholds of the network to achieve self-adjustment [12]. BP neural network is a nonlinear adaptive system, so it is more suitable for dealing with fuzzy or nonlinear problems. This method can effectively reduce subjective factors in the evaluation process and reduce the evaluation time [13].Figure 1
Structure of BP neural network.The main parameter settings of the BP neural network are the number of network layers, the number of nodes in the hidden layer, the transfer function, and the training function, which will be introduced one by one as follows.


### 2.3.1. Number of Network Layers

As the number of network layers increases, the structure of the BP neural network will become more and more complex. Correspondingly, the complex BP neural network will prolong the learning time and cause the phenomenon of “overfitting.” Through previous testing and research, the neural network is generally set up as a three-layer network, that is, the input layer, the hidden layer, and the output layer [14].


### 2.3.2. Number of Hidden Layer Nodes

Too few hidden layer neurons may not train the desired network, or the trained network is not strong enough and has poor generalization ability. To contrary, too many hidden layer neurons will increase the learning time and the error may not be smaller. The number of neurons in the hidden layer can be obtained by the following empirical formula [15]:(11)b=t+o+a,where b is the number of neurons in the hidden layer, t is the number of neurons in the input layer, o is the number of neurons in the output layer, and a is a constant between [1, 10].


### 2.3.3. Transfer Function

There are three main transfer functions: logsig function, tansig function, and purelin function. The choice of the transfer function of the hidden layer and the output layer has a greater impact on the prediction accuracy of the BP neural network. The logsig transfer function is an S-type logarithmic function, the tansig function is an S-type hyperbolic tangent function, and both are nonlinear functions. After that, the purelin function is a linear function. Generally, the transfer function of the hidden layer selects the logsig function or the tansig function, and the transfer function of the output layer selects the purelin function. As far as the nonlinear transfer function is concerned, if the output of the samples is greater than zero, the logsig function is mostly used; otherwise, the tansig function is used.


### 2.3.4. Training Function

The common training functions are as follows:trainlm: LevenbergMarquardt methodtraingd: gradient descent methodtraingdm: gradient descent method with momentum factortraingda: gradient descent method with adaptive learning ratetraingdx: gradient descent method with adaptive learning rate and momentum factor


## 2.3.1. Number of Network Layers

As the number of network layers increases, the structure of the BP neural network will become more and more complex. Correspondingly, the complex BP neural network will prolong the learning time and cause the phenomenon of “overfitting.” Through previous testing and research, the neural network is generally set up as a three-layer network, that is, the input layer, the hidden layer, and the output layer [14].


## 2.3.2. Number of Hidden Layer Nodes

Too few hidden layer neurons may not train the desired network, or the trained network is not strong enough and has poor generalization ability. To contrary, too many hidden layer neurons will increase the learning time and the error may not be smaller. The number of neurons in the hidden layer can be obtained by the following empirical formula [15]:(11)b=t+o+a,where b is the number of neurons in the hidden layer, t is the number of neurons in the input layer, o is the number of neurons in the output layer, and a is a constant between [1, 10].


## 2.3.3. Transfer Function

There are three main transfer functions: logsig function, tansig function, and purelin function. The choice of the transfer function of the hidden layer and the output layer has a greater impact on the prediction accuracy of the BP neural network. The logsig transfer function is an S-type logarithmic function, the tansig function is an S-type hyperbolic tangent function, and both are nonlinear functions. After that, the purelin function is a linear function. Generally, the transfer function of the hidden layer selects the logsig function or the tansig function, and the transfer function of the output layer selects the purelin function. As far as the nonlinear transfer function is concerned, if the output of the samples is greater than zero, the logsig function is mostly used; otherwise, the tansig function is used.


## 2.3.4. Training Function

The common training functions are as follows:trainlm: LevenbergMarquardt methodtraingd: gradient descent methodtraingdm: gradient descent method with momentum factortraingda: gradient descent method with adaptive learning ratetraingdx: gradient descent method with adaptive learning rate and momentum factor


## 3. Data Processing

### 3.1. Data Sources

The abovementioned evaluation index system was determined by analyzing the construction status of 11 emergency rescue teams for hazardous chemical accidents in Shanghai. Then, this paper produced a score sheet for these professional rescue teams. To reduce the subjective arbitrariness when scoring, the results were divided into 5 levels, and each level was guaranteed to take a positive value, from the worst to the best, respectively, 1, 2, 3, 4, and 5 points. See Table3 for scoring criterion. The scorer could select the corresponding score according to the actual situation. Table 4 shows the scores of the 11 rescue teams:Table 3
The scoring criteria of the original evaluation index.
Evaluation indexPointsScoring criteriaPostsetting D15The emergency team has necessary positions, and the position setting is very reasonable, so that the team can operate most efficiently4The emergency team has necessary positions, and the position setting is reasonable, so that the team can operate effectively3The emergency team has necessary positions, and the position setting is not very reasonable, so that the team can operate2The emergency team does not fully satisfy the necessary positions, but it does not affect the operation of the team1The emergency team did not meet the necessary positions and affected the operation of the teamClarity of post responsibilities D25Very clear and no ambiguous parts4Clear, but a little bit ambiguous3More clear, but there are still some unclear parts2Partly clear, but there are still many unclear parts1UnclearDuty fulfillmentD35Personnel in all positions perform job duties4The majority of positions perform their duties3More than half of the staff perform job duties2Personnel in some positions perform job duties1A very small number of personnel perform job dutiesSmooth communication D45Very smooth, no need for improvement4Smooth, best if it can be improved3It is smoother and can be improved2Generally smooth, needs improvement1Not smoothTimeliness of information transfer D55Timely and effective4Not very timely but still effective3Timely but not very effective2Timely but almost ineffective1Not timely and ineffectiveIntegrity of information transfer D65Both accurate and complete4Accurate and relatively complete3Accurate but not very complete2Not very accurate but complete1Inaccurate and incompleteReward and punishment system D75Effectively mobilize the enthusiasm of the team members, and the effect is very good4Aroused the enthusiasm of the team members, and the effect was good3Aroused the enthusiasm of the team members, but the effect was average2Aroused the enthusiasm of the team members, but the effect was very small1Cannot mobilize the enthusiasm of team membersSalary and benefits D85Very satisfied4Satisfied3Generally satisfied2Not very satisfied1Not satisfiedCareer opportunities D95Have career development prospects, can maintain the stability of the team members, and attract new team members to join4Have career development prospects, can maintain the stability of the team members, but cannot attract new team members to join3Have certain career development prospects and can basically maintain the stability of the team members2There are very few career development prospects, and it is difficult to maintain the stability of the team members1There is almost no career development prospects, and the stability of team members cannot be maintainedBasic intercom equipment D105Very satisfied4Satisfied3Generally satisfied2Not very satisfied1Not satisfiedAdvanced communication system and platform D115Have1Have notWhether the emergency response plan is scientific and accurate D125The emergency response plan is highly scientific and correct and can handle emergency work very effectively4The emergency response plan has strong scientificity and accuracy and can handle emergency work3The emergency response plan is generally scientific and correct; although it can handle emergency work, there are some problems2The scientificity and accuracy of the emergency response plan are low, and there are problems and difficulties in handling emergency work1The emergency response plan is hardly scientific and correct and cannot handle emergency workWhether the emergency response plan is targeted D135The emergency response plan has strong pertinence4The emergency response plan has relatively strong pertinence3The emergency response plan is generally targeted2The emergency response plan is less targeted1The emergency response plan is hardly targetedWhether the emergency treatment plan is operable D145The emergency response plan has strong operability4The emergency response plan has relatively strong operability3The operability of the emergency response plan is general2The operability of the emergency response plan is weak1The emergency response plan is almost inoperableWhether emergency experts are qualified for emergency work D155Totally satisfied4Relatively satisfied3Generally satisfied2Not very satisfied1Not satisfiedWhether there is sufficient staff D165Totally satisfied4Relatively satisfied3Generally satisfied2Not very satisfied1Not satisfiedCompleteness of emergency rescue equipment D175Fully equipped and able to perform all emergency rescue activities4The equipment is relatively complete and can perform most of the emergency rescue activities3The equipment needs to be supplemented and perfected but can perform basic emergency rescue activities2The equipment is not perfect, and there are certain difficulties in performing emergency rescue activities1The equipment is not perfect, and there are difficulties in performing emergency rescue activitiesCompleteness of personal protective equipment D185The equipment is perfect, and the rescue team is hardly injured4The equipment is relatively complete, and the rescue team has safety guarantee3The equipment still needs to be supplemented and perfected, but the rescue team can perform basic emergency rescue tasks relatively safely2The equipment is not perfect, and there are certain difficulties for rescue team members to perform emergency rescue tasks safely1The equipment is imperfect, and the rescue team cannot get effective safety protectionThe ability of emergency rescue equipment to control and reduce the impact of an accident D195It has obvious effects on controlling the source of danger and reducing the consequences of accidents4It has a relatively obvious effect on controlling the source of danger and reducing the consequences of accidents3It has a certain obvious effect on controlling the source of danger and reducing the consequences of accidents2Has a small effect on controlling the source of danger and reducing the consequences of accidents1No effect on controlling hazards and reducing the consequences of accidentsProtective capabilities of personal protective equipment D205Strong protection4Relatively strong protection3General protection2Relatively weak protection1Weak protectionCompleteness of emergency materials D215A complete range of emergency supplies and can perform all emergency activities4A relatively complete range of emergency supplies and can perform most emergency activities3A complete range of emergency supplies and can perform basic emergency activities, but the types of emergency supplies still need to be optimized2The types of emergency supplies are not very complete, and there are difficulties in implementing emergency activities1The types of emergency supplies are not complete, and it is difficult to implement emergency activitiesReplenishment of emergency materials D225Emergency supplies can be supplemented in time4Emergency supplies can be replenished in a short time3Emergency supplies cannot be supplemented in the short term, but subsequent supplementation will not affect the normal emergency response2Emergency supplies will take a long time to be replenished1Emergency supplies cannot be replenishedEmergency funding support in Shanghai D235Give strong support and provide sufficient capital guarantee4Give relatively strong support and provide more financial protection3Give a certain level of support2Give very little support1Almost no supportTerritory emergency funding support D245Give strong support and provide sufficient capital guarantee4Give relatively strong support and provide more financial protection3Give a certain level of support2Give very little support1Almost no supportEmergency rescue expertise training D255There is continuous professional knowledge training, and the effect is very good4There is continuous professional knowledge training, and the effect is relatively good3There is continuous professional knowledge training, but the effect is average2There is continuous professional knowledge training, but the effect is very small1There is no continuous professional knowledge trainingEmergency rescue professional skills training D265There is continuous professional skills training, and the effect is very good4There is continuous professional skills training, and the effect is relatively good3There is continuous professional skills training, but the effect is average2There is continuous professional skills training, but the effect is very small1There is no continuous professional skill trainingWorkload of emergency drill D275Fully meet the requirements4More satisfying requirements3Generally meet the requirements2Rarely meet the requirements1Does not meet the requirementsEffect of emergency drill D285The effect is obvious4The effect is relatively obvious3General effect2Has little effect1Has no effectTable 4
Scoring table of professional rescue team for hazardous chemical accidents.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Team 9Team 10Team 11D113445453554D254545453553D343545453554D434434443333D535555453334D634444553333D753424541314D824424443333D923523433333D1044434443334D1111555151555D1234535353334D1333545453334D1433444553334D1534545341314D1632434443334D1733534343335D1832445343334D1932445353534D2033544453534D2133534343535D2223545555535D2321331313333D2422332313313D2533514251315D2633524151335D2742445453335D2832444453334


### 3.2. Data Dimensionality Reduction and Weight Determination

In this research, SPSS 25.0 software was used to perform the principal component analysis (PCA) of the original data to achieve data reduction and weight determination. The eigenvalues of the correlation coefficient matrix between 28 indicators, the contribution rate and the cumulative contribution rate of the principal components (Table5), and the factor load matrix (Table 6) could be automatically calculated through the factor analysis tool in SPSS software. As can be seen from Table 5, there are 7 eigenvalues greater than 1, and the cumulative contribution rate of the corresponding 7 principal components reaches 94.743%, which meets the requirement that the cumulative contribution rate is generally greater than 85%. However, as shown in Table 6, none of the last two of the 7 principal components exceeds 0.5, and the cumulative contribution rate of the first 5 principal components reaches 85.849%, which also meets the requirement of greater than 85%. Therefore, this study chose the first 5 principal components as input units of the final BP neural network. Although the unit orthogonal eigenvector could not be directly obtained by the factor analysis tool of SPSS software, it could be calculated by the relationship formula (12) between the unit orthogonal eigenvector and the factor load:(12)atj=Ltλt,t=1,2,…,q;j=1,2,…,m,where Lt is the load of principal component Pt. So the coefficients of the five principal component linear expressions were known (Table 7).Table 5
Eigenvalue of the correlation coefficient matrix, contribution rate, and cumulative contribution rate of principal components.
Principal componentsEigen valueContribution rate (%)Cumulative contribution rate (%)P112.43144.39544.395P24.24215.15059.544P32.84810.17269.716P42.3828.50578.222P52.1357.62785.849P61.4135.04690.894P71.0783.84994.743P80.7132.54697.288P90.4891.74799.036P100.2700.964100.000P119.714 × 10–163.469 × 10–15100.000P125.642 × 10–162.015 × 10–15100.000P134.890× 10–161.746× 10–15100.000P143.687× 10–161.317× 10–15100.000P153.357× 10–161.199× 10–15100.000P162.120× 10–167.570× 10–16100.000P171.229× 10–164.390× 10–16100.000P185.703× 10–172.037× 10–16100.000P19−4.248× 10–17−1.517× 10–16100.000P20−1.253× 10–16−4.475× 10–16100.000P21−1.701× 10–16−6.075× 10–16100.000P22−2.841× 10–16−1.015× 10–15100.000P23−3.203× 10–16−1.144× 10–15100.000P24−5.065× 10–16−1.809× 10–15100.000P25−5.375× 10–16−1.920× 10–15100.000P26−8.421× 10–16−3.008× 10–15100.000P27−8.752× 10–16−3.126× 10–15100.000P28−2.195× 10–15−7.839× 10–15100.000Table 6
Load matrix of principal components.
P1P2P3P4P5P6P7D10.5220.444−0.044−0.5790.267−0.136−0.078D20.3120.001−0.688−0.2860.0780.4660.249D30.6210.444−0.437−0.2480.0750.3680.021D40.615−0.6800.086−0.2960.2190.0340.039D50.661−0.4180.029−0.152−0.041−0.4910.236D60.561−0.5130.260−0.514−0.1750.0480.158D70.545−0.4240.0590.381−0.2260.4910.218D80.572−0.4690.149−0.3110.5030.023−0.198D90.478−0.1190.408−0.0720.6280.286−0.161D100.570−0.686−0.0080.373−0.0780.180−0.098D110.5920.646−0.331−0.1180.073−0.2090.023D120.849−0.279−0.2150.0690.190−0.253−0.140D130.944−0.0580.093−0.112−0.158−0.107−0.034D140.742−0.1510.393−0.277−0.3540.112−0.038D150.764−0.274−0.0680.250−0.026−0.2280.430D160.8120.2070.2890.035−0.2210.273−0.279D170.7860.1100.0610.4620.170−0.140−0.272D180.7610.347−0.0390.027−0.361−0.214−0.073D190.7110.541−0.173−0.087−0.115−0.0230.220D200.7850.3230.095−0.0930.1790.1480.374D210.6550.431−0.0230.4040.4240.0410.117D220.6090.3370.567−0.1550.199−0.124−0.028D23−0.2940.6520.5190.1050.1150.239−0.025D24−0.0540.2550.7560.3740.119−0.0630.414D250.809−0.158−0.1830.4580.1750.039−0.028D260.7270.052−0.4760.3660.181−0.096−0.151D270.7580.1710.0420.179−0.5660.056−0.174D280.8020.2180.194−0.161−0.4370.057−0.110Table 7
Coefficients of linear expressions of principal components.
P1P2P3P4P5D10.1480.215−0.026−0.3750.182D20.0880.001−0.408−0.1850.053D30.1760.216−0.259−0.1610.051D40.175−0.3300.051−0.1920.150D50.187−0.2030.017−0.099−0.028D60.159−0.2490.154−0.333−0.120D70.155−0.2060.0350.247−0.154D80.162−0.2280.088−0.2010.344D90.136−0.0580.242−0.0460.430D100.162−0.333−0.0050.242−0.054D110.1680.314−0.196−0.0760.050D120.241−0.135−0.1280.0450.130D130.268−0.0280.055−0.073−0.108D140.211−0.0730.233−0.180−0.242D150.217−0.133−0.0410.162−0.018D160.2300.1000.1710.023−0.151D170.2230.0530.0360.2990.116D180.2160.169−0.0230.017−0.247D190.2020.263−0.103−0.056−0.079D200.2230.1570.056−0.0600.123D210.1860.209−0.0130.2620.290D220.1730.1630.336−0.1000.136D23−0.0830.3160.3080.0680.079D24−0.0150.1240.4480.2420.082D250.229−0.077−0.1090.2970.119D260.2060.025−0.2820.2370.124D270.2150.0830.0250.116−0.388D280.2280.1060.115−0.105−0.299The contribution rates of the 5 principal components are weights, and the normalized weights are shown in Table8. The specific values of the 5 principal components in the 11 rescue teams were obtained through Table 7 and formula (9), and the comprehensive scores of the emergency capabilities of the 11 rescue teams were finally calculated through formula (10), as shown in Table 9:Table 8
Weight of principal components.
P1P2P3P4P5Weight (%)51.7117.6511.859.918.88Table 9
Scores of 11 professional rescue teams represented by 5 principal components.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Team 9Team 10Team 11P114.69214.9922.36316.77521.98617.47522.82413.05417.38614.64620.322P2−0.436−2.5351.9083.1161.631−0.251.361.8634.4593.4593.038P30.3681.2191.6582.2220.2955.0340.0784.1891.105−0.1132.022P41.305−0.8360.252−2.469−1.367−2.422−1.932−1.889−1.457−3.4451.742P5−0.4152.9663.591−0.1811.1870.7911.0751.7433.3252.3271.713Comprehensive score7.6577.63112.4419.22511.6629.41811.9567.54110.0588.03511.607As can be seen from Table9, among the 11 rescue teams, team 3 gets the highest score of 12.441 points, while team 8 gets the lowest score of only 7.541 points. The difference between the highest and lowest scores is 4.9 points. Compared with the original scoring Table 4, it can be found that team 3 scores 4 or 5 points except D23 and D24. The difference is that for team 8, except for D22, all other indicators are 1 or 3 points. Therefore, the comprehensive scores obtained by the principal component analysis method are consistent with actual situations and can effectively reflect the emergency response status of 11 rescue teams.


## 3.1. Data Sources

The abovementioned evaluation index system was determined by analyzing the construction status of 11 emergency rescue teams for hazardous chemical accidents in Shanghai. Then, this paper produced a score sheet for these professional rescue teams. To reduce the subjective arbitrariness when scoring, the results were divided into 5 levels, and each level was guaranteed to take a positive value, from the worst to the best, respectively, 1, 2, 3, 4, and 5 points. See Table3 for scoring criterion. The scorer could select the corresponding score according to the actual situation. Table 4 shows the scores of the 11 rescue teams:Table 3
The scoring criteria of the original evaluation index.
Evaluation indexPointsScoring criteriaPostsetting D15The emergency team has necessary positions, and the position setting is very reasonable, so that the team can operate most efficiently4The emergency team has necessary positions, and the position setting is reasonable, so that the team can operate effectively3The emergency team has necessary positions, and the position setting is not very reasonable, so that the team can operate2The emergency team does not fully satisfy the necessary positions, but it does not affect the operation of the team1The emergency team did not meet the necessary positions and affected the operation of the teamClarity of post responsibilities D25Very clear and no ambiguous parts4Clear, but a little bit ambiguous3More clear, but there are still some unclear parts2Partly clear, but there are still many unclear parts1UnclearDuty fulfillmentD35Personnel in all positions perform job duties4The majority of positions perform their duties3More than half of the staff perform job duties2Personnel in some positions perform job duties1A very small number of personnel perform job dutiesSmooth communication D45Very smooth, no need for improvement4Smooth, best if it can be improved3It is smoother and can be improved2Generally smooth, needs improvement1Not smoothTimeliness of information transfer D55Timely and effective4Not very timely but still effective3Timely but not very effective2Timely but almost ineffective1Not timely and ineffectiveIntegrity of information transfer D65Both accurate and complete4Accurate and relatively complete3Accurate but not very complete2Not very accurate but complete1Inaccurate and incompleteReward and punishment system D75Effectively mobilize the enthusiasm of the team members, and the effect is very good4Aroused the enthusiasm of the team members, and the effect was good3Aroused the enthusiasm of the team members, but the effect was average2Aroused the enthusiasm of the team members, but the effect was very small1Cannot mobilize the enthusiasm of team membersSalary and benefits D85Very satisfied4Satisfied3Generally satisfied2Not very satisfied1Not satisfiedCareer opportunities D95Have career development prospects, can maintain the stability of the team members, and attract new team members to join4Have career development prospects, can maintain the stability of the team members, but cannot attract new team members to join3Have certain career development prospects and can basically maintain the stability of the team members2There are very few career development prospects, and it is difficult to maintain the stability of the team members1There is almost no career development prospects, and the stability of team members cannot be maintainedBasic intercom equipment D105Very satisfied4Satisfied3Generally satisfied2Not very satisfied1Not satisfiedAdvanced communication system and platform D115Have1Have notWhether the emergency response plan is scientific and accurate D125The emergency response plan is highly scientific and correct and can handle emergency work very effectively4The emergency response plan has strong scientificity and accuracy and can handle emergency work3The emergency response plan is generally scientific and correct; although it can handle emergency work, there are some problems2The scientificity and accuracy of the emergency response plan are low, and there are problems and difficulties in handling emergency work1The emergency response plan is hardly scientific and correct and cannot handle emergency workWhether the emergency response plan is targeted D135The emergency response plan has strong pertinence4The emergency response plan has relatively strong pertinence3The emergency response plan is generally targeted2The emergency response plan is less targeted1The emergency response plan is hardly targetedWhether the emergency treatment plan is operable D145The emergency response plan has strong operability4The emergency response plan has relatively strong operability3The operability of the emergency response plan is general2The operability of the emergency response plan is weak1The emergency response plan is almost inoperableWhether emergency experts are qualified for emergency work D155Totally satisfied4Relatively satisfied3Generally satisfied2Not very satisfied1Not satisfiedWhether there is sufficient staff D165Totally satisfied4Relatively satisfied3Generally satisfied2Not very satisfied1Not satisfiedCompleteness of emergency rescue equipment D175Fully equipped and able to perform all emergency rescue activities4The equipment is relatively complete and can perform most of the emergency rescue activities3The equipment needs to be supplemented and perfected but can perform basic emergency rescue activities2The equipment is not perfect, and there are certain difficulties in performing emergency rescue activities1The equipment is not perfect, and there are difficulties in performing emergency rescue activitiesCompleteness of personal protective equipment D185The equipment is perfect, and the rescue team is hardly injured4The equipment is relatively complete, and the rescue team has safety guarantee3The equipment still needs to be supplemented and perfected, but the rescue team can perform basic emergency rescue tasks relatively safely2The equipment is not perfect, and there are certain difficulties for rescue team members to perform emergency rescue tasks safely1The equipment is imperfect, and the rescue team cannot get effective safety protectionThe ability of emergency rescue equipment to control and reduce the impact of an accident D195It has obvious effects on controlling the source of danger and reducing the consequences of accidents4It has a relatively obvious effect on controlling the source of danger and reducing the consequences of accidents3It has a certain obvious effect on controlling the source of danger and reducing the consequences of accidents2Has a small effect on controlling the source of danger and reducing the consequences of accidents1No effect on controlling hazards and reducing the consequences of accidentsProtective capabilities of personal protective equipment D205Strong protection4Relatively strong protection3General protection2Relatively weak protection1Weak protectionCompleteness of emergency materials D215A complete range of emergency supplies and can perform all emergency activities4A relatively complete range of emergency supplies and can perform most emergency activities3A complete range of emergency supplies and can perform basic emergency activities, but the types of emergency supplies still need to be optimized2The types of emergency supplies are not very complete, and there are difficulties in implementing emergency activities1The types of emergency supplies are not complete, and it is difficult to implement emergency activitiesReplenishment of emergency materials D225Emergency supplies can be supplemented in time4Emergency supplies can be replenished in a short time3Emergency supplies cannot be supplemented in the short term, but subsequent supplementation will not affect the normal emergency response2Emergency supplies will take a long time to be replenished1Emergency supplies cannot be replenishedEmergency funding support in Shanghai D235Give strong support and provide sufficient capital guarantee4Give relatively strong support and provide more financial protection3Give a certain level of support2Give very little support1Almost no supportTerritory emergency funding support D245Give strong support and provide sufficient capital guarantee4Give relatively strong support and provide more financial protection3Give a certain level of support2Give very little support1Almost no supportEmergency rescue expertise training D255There is continuous professional knowledge training, and the effect is very good4There is continuous professional knowledge training, and the effect is relatively good3There is continuous professional knowledge training, but the effect is average2There is continuous professional knowledge training, but the effect is very small1There is no continuous professional knowledge trainingEmergency rescue professional skills training D265There is continuous professional skills training, and the effect is very good4There is continuous professional skills training, and the effect is relatively good3There is continuous professional skills training, but the effect is average2There is continuous professional skills training, but the effect is very small1There is no continuous professional skill trainingWorkload of emergency drill D275Fully meet the requirements4More satisfying requirements3Generally meet the requirements2Rarely meet the requirements1Does not meet the requirementsEffect of emergency drill D285The effect is obvious4The effect is relatively obvious3General effect2Has little effect1Has no effectTable 4
Scoring table of professional rescue team for hazardous chemical accidents.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Team 9Team 10Team 11D113445453554D254545453553D343545453554D434434443333D535555453334D634444553333D753424541314D824424443333D923523433333D1044434443334D1111555151555D1234535353334D1333545453334D1433444553334D1534545341314D1632434443334D1733534343335D1832445343334D1932445353534D2033544453534D2133534343535D2223545555535D2321331313333D2422332313313D2533514251315D2633524151335D2742445453335D2832444453334


## 3.2. Data Dimensionality Reduction and Weight Determination

In this research, SPSS 25.0 software was used to perform the principal component analysis (PCA) of the original data to achieve data reduction and weight determination. The eigenvalues of the correlation coefficient matrix between 28 indicators, the contribution rate and the cumulative contribution rate of the principal components (Table5), and the factor load matrix (Table 6) could be automatically calculated through the factor analysis tool in SPSS software. As can be seen from Table 5, there are 7 eigenvalues greater than 1, and the cumulative contribution rate of the corresponding 7 principal components reaches 94.743%, which meets the requirement that the cumulative contribution rate is generally greater than 85%. However, as shown in Table 6, none of the last two of the 7 principal components exceeds 0.5, and the cumulative contribution rate of the first 5 principal components reaches 85.849%, which also meets the requirement of greater than 85%. Therefore, this study chose the first 5 principal components as input units of the final BP neural network. Although the unit orthogonal eigenvector could not be directly obtained by the factor analysis tool of SPSS software, it could be calculated by the relationship formula (12) between the unit orthogonal eigenvector and the factor load:(12)atj=Ltλt,t=1,2,…,q;j=1,2,…,m,where Lt is the load of principal component Pt. So the coefficients of the five principal component linear expressions were known (Table 7).Table 5
Eigenvalue of the correlation coefficient matrix, contribution rate, and cumulative contribution rate of principal components.
Principal componentsEigen valueContribution rate (%)Cumulative contribution rate (%)P112.43144.39544.395P24.24215.15059.544P32.84810.17269.716P42.3828.50578.222P52.1357.62785.849P61.4135.04690.894P71.0783.84994.743P80.7132.54697.288P90.4891.74799.036P100.2700.964100.000P119.714 × 10–163.469 × 10–15100.000P125.642 × 10–162.015 × 10–15100.000P134.890× 10–161.746× 10–15100.000P143.687× 10–161.317× 10–15100.000P153.357× 10–161.199× 10–15100.000P162.120× 10–167.570× 10–16100.000P171.229× 10–164.390× 10–16100.000P185.703× 10–172.037× 10–16100.000P19−4.248× 10–17−1.517× 10–16100.000P20−1.253× 10–16−4.475× 10–16100.000P21−1.701× 10–16−6.075× 10–16100.000P22−2.841× 10–16−1.015× 10–15100.000P23−3.203× 10–16−1.144× 10–15100.000P24−5.065× 10–16−1.809× 10–15100.000P25−5.375× 10–16−1.920× 10–15100.000P26−8.421× 10–16−3.008× 10–15100.000P27−8.752× 10–16−3.126× 10–15100.000P28−2.195× 10–15−7.839× 10–15100.000Table 6
Load matrix of principal components.
P1P2P3P4P5P6P7D10.5220.444−0.044−0.5790.267−0.136−0.078D20.3120.001−0.688−0.2860.0780.4660.249D30.6210.444−0.437−0.2480.0750.3680.021D40.615−0.6800.086−0.2960.2190.0340.039D50.661−0.4180.029−0.152−0.041−0.4910.236D60.561−0.5130.260−0.514−0.1750.0480.158D70.545−0.4240.0590.381−0.2260.4910.218D80.572−0.4690.149−0.3110.5030.023−0.198D90.478−0.1190.408−0.0720.6280.286−0.161D100.570−0.686−0.0080.373−0.0780.180−0.098D110.5920.646−0.331−0.1180.073−0.2090.023D120.849−0.279−0.2150.0690.190−0.253−0.140D130.944−0.0580.093−0.112−0.158−0.107−0.034D140.742−0.1510.393−0.277−0.3540.112−0.038D150.764−0.274−0.0680.250−0.026−0.2280.430D160.8120.2070.2890.035−0.2210.273−0.279D170.7860.1100.0610.4620.170−0.140−0.272D180.7610.347−0.0390.027−0.361−0.214−0.073D190.7110.541−0.173−0.087−0.115−0.0230.220D200.7850.3230.095−0.0930.1790.1480.374D210.6550.431−0.0230.4040.4240.0410.117D220.6090.3370.567−0.1550.199−0.124−0.028D23−0.2940.6520.5190.1050.1150.239−0.025D24−0.0540.2550.7560.3740.119−0.0630.414D250.809−0.158−0.1830.4580.1750.039−0.028D260.7270.052−0.4760.3660.181−0.096−0.151D270.7580.1710.0420.179−0.5660.056−0.174D280.8020.2180.194−0.161−0.4370.057−0.110Table 7
Coefficients of linear expressions of principal components.
P1P2P3P4P5D10.1480.215−0.026−0.3750.182D20.0880.001−0.408−0.1850.053D30.1760.216−0.259−0.1610.051D40.175−0.3300.051−0.1920.150D50.187−0.2030.017−0.099−0.028D60.159−0.2490.154−0.333−0.120D70.155−0.2060.0350.247−0.154D80.162−0.2280.088−0.2010.344D90.136−0.0580.242−0.0460.430D100.162−0.333−0.0050.242−0.054D110.1680.314−0.196−0.0760.050D120.241−0.135−0.1280.0450.130D130.268−0.0280.055−0.073−0.108D140.211−0.0730.233−0.180−0.242D150.217−0.133−0.0410.162−0.018D160.2300.1000.1710.023−0.151D170.2230.0530.0360.2990.116D180.2160.169−0.0230.017−0.247D190.2020.263−0.103−0.056−0.079D200.2230.1570.056−0.0600.123D210.1860.209−0.0130.2620.290D220.1730.1630.336−0.1000.136D23−0.0830.3160.3080.0680.079D24−0.0150.1240.4480.2420.082D250.229−0.077−0.1090.2970.119D260.2060.025−0.2820.2370.124D270.2150.0830.0250.116−0.388D280.2280.1060.115−0.105−0.299The contribution rates of the 5 principal components are weights, and the normalized weights are shown in Table8. The specific values of the 5 principal components in the 11 rescue teams were obtained through Table 7 and formula (9), and the comprehensive scores of the emergency capabilities of the 11 rescue teams were finally calculated through formula (10), as shown in Table 9:Table 8
Weight of principal components.
P1P2P3P4P5Weight (%)51.7117.6511.859.918.88Table 9
Scores of 11 professional rescue teams represented by 5 principal components.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Team 9Team 10Team 11P114.69214.9922.36316.77521.98617.47522.82413.05417.38614.64620.322P2−0.436−2.5351.9083.1161.631−0.251.361.8634.4593.4593.038P30.3681.2191.6582.2220.2955.0340.0784.1891.105−0.1132.022P41.305−0.8360.252−2.469−1.367−2.422−1.932−1.889−1.457−3.4451.742P5−0.4152.9663.591−0.1811.1870.7911.0751.7433.3252.3271.713Comprehensive score7.6577.63112.4419.22511.6629.41811.9567.54110.0588.03511.607As can be seen from Table9, among the 11 rescue teams, team 3 gets the highest score of 12.441 points, while team 8 gets the lowest score of only 7.541 points. The difference between the highest and lowest scores is 4.9 points. Compared with the original scoring Table 4, it can be found that team 3 scores 4 or 5 points except D23 and D24. The difference is that for team 8, except for D22, all other indicators are 1 or 3 points. Therefore, the comprehensive scores obtained by the principal component analysis method are consistent with actual situations and can effectively reflect the emergency response status of 11 rescue teams.


## 4. Implementation of BP Neural Network

### 4.1. Sample Data Normalization

Normalization refers to limiting the input and output data of the network to [0, 1] or [−1, 1] through the processing of variables, which can improve the efficiency of the transfer function and the accuracy of the output of the neural network. There is a maximum-minimum method to limit the data between [0, 1], and the function form is as follows:(13)xi=xi−xminxmax−xmin.The formula that limits the data to [−1, 1] is(14)xi=2×xi−xminxmax−xmin−1.In this paper, the mapminmax function provided by MATLAB was used to obtain the normalized input and output data between [−1, 1]. The normalized data are shown in Table10:Table 10
Normalized statistics of 11 professional rescue teams for hazardous chemical accidents.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Team 9Team 10Team 11P1−0.665−0.6040.906−0.2380.828−0.0951.000−1.000−0.113−0.6740.488P2−0.400−1.0000.2710.6160.191−0.3470.1140.2581.0000.7140.594P3−0.813−0.482−0.312−0.093−0.8411.000−0.9260.672−0.527−1.000−0.170P40.8320.0060.425−0.624−0.199−0.606−0.417−0.400−0.233−1.0001.000P5−1.0000.6881.000−0.883−0.200−0.398−0.2560.0770.8670.3690.062Comprehensive score−0.953−0.9631.000−0.3130.682−0.2340.802−1.0000.027−0.7980.660Teams 1 to 8 were used as training samples for the BP neural network, and teams 9 to 11 were used as prediction samples for the BP neural network.


### 4.2. Determination of the Transfer Function

In this study, the data were limited to [−1, 1] during normalization. It could be known from the above description of the transfer function that the tansig function should be used. Therefore, it was determined that the transfer function of the hidden layer was the tansig function, and the transfer function of the output layer was selected as the purelin function.


### 4.3. Selection of Training Function

In order to determine the fast and accurate training function, this research used BP neural network toolbox of MATLAB software to experiment the above five training functions and then compared the training results to choose. According to the previous principal component analysis (PCA), five input neurons have been identified, namedP1, P2, P3, P4, and P5, and only one output neuron which was “comprehensive score.” From the above empirical formula, it could be known that the number of hidden layer neurons should be selected between [4, 13], and it was temporarily determined to be 9. The number of iterations and convergence accuracy were used as the evaluation indicators for the training function selection. Before using the above five training functions to perform prediction fitting on 8 training samples, the maximum number of iterations was set to 2000 and the target convergence accuracy was set to 0. The results are shown in Table 11:Table 11
Results of five training functions on training samples.
Training functionAlgorithmNumber of iterationsConvergence accuracytrainlmLevenbergMarquardt method54.4339× 10–24traingdGradient descent method20000.0022997traingdmGradient descent method with momentum factor20000.0021524traingdaGradient descent method with adaptive learning rate20005.854× 10–08traingdxGradient descent method with adaptive learning rate and momentum factor3138.3968× 10–10From the above training results, it can be seen that the trainlm training function achieves high accuracy in only 5 steps. The trainlm training function has the fastest convergence speed, but it is easy to fall into a local minimum. The traingd function and trackingdm function converge slowly in practical applications, and the convergence accuracy is not as high as the other three training functions. The traingda function and traingdx function have greatly improved the convergence accuracy of training, but the traingdx function converges faster than the traingda function and the traingdx function can avoid falling into local minima due to the additional momentum term and adaptive learning rate. At the same time, the training precision of the traingdx function reaches 8.3968× 10−10, which is consistent with the convergence accuracy in general cases. Therefore, the BP neural network model established in this paper used the traingdx function as the training function.


### 4.4. Setting of Training Parameters

In this study, two parameters of the BP neural network model were set. The maximum allowable error was set to 0.00001, the maximum number of learning times was set to 1000, and the remaining parameters adopted default values.


### 4.5. Determination of the Number of Hidden Layer Neurons

The value range of hidden layer neurons has been obtained through the empirical formula in the previous article. However, if the number of hidden layer neurons is too small, the ability of the neural network to obtain information from the sample is poor, and it is impossible to generalize and reflect the sample law. At the same time, if there are too many neurons, the irregular content in the sample may be learned and the phenomenon of “overfitting” may occur. Therefore, it is necessary to determine an optimal number of hidden layer neurons.The number of hidden layer neurons can be determined one by one [4, 13] through experiments. The number of training sessions was set to 1000, and the target accuracy was set to 0.00001. The number of iterations and the convergence accuracy of 10 trainings were compared to determine the appropriate number of hidden layer neurons. The training results are shown in Table 12:Table 12
Training results of different numbers of neurons in the hidden layer.
Number of hidden neuronsNumber of iterationsConvergence accuracy42689.7354× 10–653379.4521× 10–661878.8167× 10–671939.9709× 10–681449.7153× 10–691258.6169× 10–6102429.7497× 10–6112408.8618× 10–6121388.7017× 10–6131205.8329× 10–6As shown in Table12, the convergence accuracy is on an order of magnitude, but there is some gap in the number of iterations. When the number of hidden layer neurons is 8, 9, 12, and 13, the number of iterations is small. From the perspective of simplifying the structure of the BP neural network, this study determined that the number of hidden layer neurons was 8. The resulting structure of BP neural network is shown in Figure 2:Figure 2
The structure of the BP neural network model.


### 4.6. Computational Complexity

The number of neurons in a neural network has an important impact on the computational complexity. When the number of neurons increases, the network calculation becomes more complex. The complexity of calculating the gradient of a certain layer isOD3, assuming that the number of neurons in this layer is D. This study uses principal component analysis to reduce the number of neurons in the input layer of the neural network and obtained the final structure of neural network: 5 neurons in the input layer, 8 neurons in the hidden layer, and 1 neuron in the output layer. Therefore, the final calculation of the entire neural network is O53+83+13. In contrast, the computational complexity of the input layer of the neural network that directly uses all the evaluation indicators as the input layer neurons is O283, which is far greater than the computational complexity of the optimized neural network.


## 4.1. Sample Data Normalization

Normalization refers to limiting the input and output data of the network to [0, 1] or [−1, 1] through the processing of variables, which can improve the efficiency of the transfer function and the accuracy of the output of the neural network. There is a maximum-minimum method to limit the data between [0, 1], and the function form is as follows:(13)xi=xi−xminxmax−xmin.The formula that limits the data to [−1, 1] is(14)xi=2×xi−xminxmax−xmin−1.In this paper, the mapminmax function provided by MATLAB was used to obtain the normalized input and output data between [−1, 1]. The normalized data are shown in Table10:Table 10
Normalized statistics of 11 professional rescue teams for hazardous chemical accidents.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Team 9Team 10Team 11P1−0.665−0.6040.906−0.2380.828−0.0951.000−1.000−0.113−0.6740.488P2−0.400−1.0000.2710.6160.191−0.3470.1140.2581.0000.7140.594P3−0.813−0.482−0.312−0.093−0.8411.000−0.9260.672−0.527−1.000−0.170P40.8320.0060.425−0.624−0.199−0.606−0.417−0.400−0.233−1.0001.000P5−1.0000.6881.000−0.883−0.200−0.398−0.2560.0770.8670.3690.062Comprehensive score−0.953−0.9631.000−0.3130.682−0.2340.802−1.0000.027−0.7980.660Teams 1 to 8 were used as training samples for the BP neural network, and teams 9 to 11 were used as prediction samples for the BP neural network.


## 4.2. Determination of the Transfer Function

In this study, the data were limited to [−1, 1] during normalization. It could be known from the above description of the transfer function that the tansig function should be used. Therefore, it was determined that the transfer function of the hidden layer was the tansig function, and the transfer function of the output layer was selected as the purelin function.


## 4.3. Selection of Training Function

In order to determine the fast and accurate training function, this research used BP neural network toolbox of MATLAB software to experiment the above five training functions and then compared the training results to choose. According to the previous principal component analysis (PCA), five input neurons have been identified, namedP1, P2, P3, P4, and P5, and only one output neuron which was “comprehensive score.” From the above empirical formula, it could be known that the number of hidden layer neurons should be selected between [4, 13], and it was temporarily determined to be 9. The number of iterations and convergence accuracy were used as the evaluation indicators for the training function selection. Before using the above five training functions to perform prediction fitting on 8 training samples, the maximum number of iterations was set to 2000 and the target convergence accuracy was set to 0. The results are shown in Table 11:Table 11
Results of five training functions on training samples.
Training functionAlgorithmNumber of iterationsConvergence accuracytrainlmLevenbergMarquardt method54.4339× 10–24traingdGradient descent method20000.0022997traingdmGradient descent method with momentum factor20000.0021524traingdaGradient descent method with adaptive learning rate20005.854× 10–08traingdxGradient descent method with adaptive learning rate and momentum factor3138.3968× 10–10From the above training results, it can be seen that the trainlm training function achieves high accuracy in only 5 steps. The trainlm training function has the fastest convergence speed, but it is easy to fall into a local minimum. The traingd function and trackingdm function converge slowly in practical applications, and the convergence accuracy is not as high as the other three training functions. The traingda function and traingdx function have greatly improved the convergence accuracy of training, but the traingdx function converges faster than the traingda function and the traingdx function can avoid falling into local minima due to the additional momentum term and adaptive learning rate. At the same time, the training precision of the traingdx function reaches 8.3968× 10−10, which is consistent with the convergence accuracy in general cases. Therefore, the BP neural network model established in this paper used the traingdx function as the training function.


## 4.4. Setting of Training Parameters

In this study, two parameters of the BP neural network model were set. The maximum allowable error was set to 0.00001, the maximum number of learning times was set to 1000, and the remaining parameters adopted default values.


## 4.5. Determination of the Number of Hidden Layer Neurons

The value range of hidden layer neurons has been obtained through the empirical formula in the previous article. However, if the number of hidden layer neurons is too small, the ability of the neural network to obtain information from the sample is poor, and it is impossible to generalize and reflect the sample law. At the same time, if there are too many neurons, the irregular content in the sample may be learned and the phenomenon of “overfitting” may occur. Therefore, it is necessary to determine an optimal number of hidden layer neurons.The number of hidden layer neurons can be determined one by one [4, 13] through experiments. The number of training sessions was set to 1000, and the target accuracy was set to 0.00001. The number of iterations and the convergence accuracy of 10 trainings were compared to determine the appropriate number of hidden layer neurons. The training results are shown in Table 12:Table 12
Training results of different numbers of neurons in the hidden layer.
Number of hidden neuronsNumber of iterationsConvergence accuracy42689.7354× 10–653379.4521× 10–661878.8167× 10–671939.9709× 10–681449.7153× 10–691258.6169× 10–6102429.7497× 10–6112408.8618× 10–6121388.7017× 10–6131205.8329× 10–6As shown in Table12, the convergence accuracy is on an order of magnitude, but there is some gap in the number of iterations. When the number of hidden layer neurons is 8, 9, 12, and 13, the number of iterations is small. From the perspective of simplifying the structure of the BP neural network, this study determined that the number of hidden layer neurons was 8. The resulting structure of BP neural network is shown in Figure 2:Figure 2
The structure of the BP neural network model.


## 4.6. Computational Complexity

The number of neurons in a neural network has an important impact on the computational complexity. When the number of neurons increases, the network calculation becomes more complex. The complexity of calculating the gradient of a certain layer isOD3, assuming that the number of neurons in this layer is D. This study uses principal component analysis to reduce the number of neurons in the input layer of the neural network and obtained the final structure of neural network: 5 neurons in the input layer, 8 neurons in the hidden layer, and 1 neuron in the output layer. Therefore, the final calculation of the entire neural network is O53+83+13. In contrast, the computational complexity of the input layer of the neural network that directly uses all the evaluation indicators as the input layer neurons is O283, which is far greater than the computational complexity of the optimized neural network.


## 5. Simulation and Result Analysis of BP Neural Network Evaluation Model

### 5.1. Training of BP Neural Network

After the above discussion, the parameters of the BP neural network evaluation model had been determined. The first 8 teams were now used as training samples to train the BP neural network. The training result is shown in Figure3. After 150 iterations, the mean-squared error of the neural network reaches 9.688e − 6. As shown in Table 13, the maximum relative error between the training results and the comprehensive score is only 0.146%, indicating that the BP neural network has reached the training requirements.Figure 3
Mean square error trace.Table 13
Relative error between training result and theoretical comprehensive score.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Comprehensive score7.6577.63112.4419.22511.6629.41811.9567.541Predicted result7.6557.63212.4389.22511.6799.41511.9437.543Relative error (%)0.0260.0130.0240.0000.1460.0320.1090.027


### 5.2. Simulation of BP Neural Network

The trained BP neural network evaluation model was used to predict the comprehensive scores of the remaining three rescue teams. The relative errors between the predicted results and the theoretical values are shown in Table14. It can be known from Table 14 that the maximum relative error is 6.658%, which indicates that the generalized ability of the trained BP neural network evaluation model can meet the emergency rescue capability assessment needs of rescue teams.Table 14
Relative error between prediction result and theoretical comprehensive score.
Team 9Team 10Team 11Comprehensive score10.0588.03511.607Predicted result10.6617.511.699Relative error (%)5.9956.6580.793


## 5.1. Training of BP Neural Network

After the above discussion, the parameters of the BP neural network evaluation model had been determined. The first 8 teams were now used as training samples to train the BP neural network. The training result is shown in Figure3. After 150 iterations, the mean-squared error of the neural network reaches 9.688e − 6. As shown in Table 13, the maximum relative error between the training results and the comprehensive score is only 0.146%, indicating that the BP neural network has reached the training requirements.Figure 3
Mean square error trace.Table 13
Relative error between training result and theoretical comprehensive score.
Team 1Team 2Team 3Team 4Team 5Team 6Team 7Team 8Comprehensive score7.6577.63112.4419.22511.6629.41811.9567.541Predicted result7.6557.63212.4389.22511.6799.41511.9437.543Relative error (%)0.0260.0130.0240.0000.1460.0320.1090.027


## 5.2. Simulation of BP Neural Network

The trained BP neural network evaluation model was used to predict the comprehensive scores of the remaining three rescue teams. The relative errors between the predicted results and the theoretical values are shown in Table14. It can be known from Table 14 that the maximum relative error is 6.658%, which indicates that the generalized ability of the trained BP neural network evaluation model can meet the emergency rescue capability assessment needs of rescue teams.Table 14
Relative error between prediction result and theoretical comprehensive score.
Team 9Team 10Team 11Comprehensive score10.0588.03511.607Predicted result10.6617.511.699Relative error (%)5.9956.6580.793


## 6. Conclusion

(1)
The capacity evaluation index system for the emergency rescue teams for hazardous chemical accidents proposed in this study was considered from five aspects: emergency cooperation, emergency command, emergency foundation, training and education, and emergency drills. A total of 5 first-level indicators, 14 second-level indicators, and 28 third-level evaluation indicators were proposed, which fully considered the factors affecting emergency rescue capabilities and were relatively comprehensive and objective.(2)
The dimensionality reduction of the evaluation index was achieved by the principal component analysis (PCA) method, and the weights and comprehensive scores of 11 emergency rescue teams were obtained. According to the original index scoring table, it is found that the comprehensive scores obtained by the principal component analysis are in line with actual situations and can well reflect the emergency rescue capabilities of the 11 rescue teams.(3)
After the dimensionality reduction by the principal component analysis (PCA) method, the number of neurons in the input layer of the BP neural network was greatly reduced and the structure of the BP neural network was simplified. At the same time, the computational complexity of the neural network evaluation model has also been reduced. Teams 1 to 8 were used as training samples for the BP neural network evaluation model, and teams 9 to 11 were used as test samples for verifying the model. The trained BP neural network evaluation model showed good generalization ability, and the highest relative error with the theoretical comprehensive score of the test sample was 6.658%, which could meet the needs of emergency rescue capability assessment.The aim of this paper is to assess the capability of the emergency rescue team for hazardous chemical accidents in order to better understand the current situation of these rescue teams. This understanding will support the improvement the capability of these teams in a targeted way. In addition, competent authorities could use such assessments to improve their management level. The above research results showed that it is feasible to use the PCA-BP neural network-based evaluation model to evaluate the capability of emergency rescue teams for hazardous chemical accidents, which provides a new idea for emergency rescue capability assessment.



---

*Source: 8869608-2021-01-15.xml*