Integrated Assessment Method of Emergency Plan for Sudden Water Pollution Accidents Based on Improved TOPSIS , Shannon Entropy and a Coordinated Development Degree Model

Water is the source of all things, so it can be said that without the sustainable development of water resources, there can be no sustainable development of human beings. In recent years, sudden water pollution accidents have occurred frequently. Emergency response plan optimization is the key to handling accidents. Nevertheless, the non-linear relationship between various indicators and emergency plans has greatly prevented researchers from making reasonable assessments. Thus, an integrated assessment method is proposed by incorporating an improved technique for order preference by similarity to ideal solution, Shannon entropy and a Coordinated development degree model to evaluate emergency plans. The Shannon entropy method was used to analyze different types of index values. TOPSIS is used to calculate the relative closeness to the ideal solution. The coordinated development degree model is applied to express the relationship between the relative closeness and inhomogeneity of the emergency plan. This method is tested in the decision support system of the Middle Route Construction and Administration Bureau, China. By considering the different nature of the indicators, the integrated assessment method is eventually proven as a highly realistic method for assessing emergency plans. The advantages of this method are more prominent when there are more indicators of the evaluation object and the nature of each indicator is quite different. In summary, this integrated assessment method can provide a targeted reference or guidance for emergency control decision makers.


Introduction
Sustainability has become an important topic in the political agenda of many countries [1][2][3].Water plays an extremely important role in human survival and regional socioeconomic development, particularly in arid and semi-arid areas [4].How to use the water resources in a sustainable way is an important problem for sustainable development, especially in areas where water resources are extremely scarce.[5][6][7][8].One of the important factors that restricts the sustainable development of water resources is uneven spatial temporal distribution [9].Therefore, a continuously growing number of water transfer projects have been built in recent years [10], such as the Central Arizona Project in the United States [11] and the South-North Water Transfer Project (SNWTP) in China [12].Open channels, complex operation condition, higher water quality requirement, diverse types and an increasing number of water pollution risk sources are the characteristics of these water transfer projects (WTPs).However, these buildings have greatly increased the pollution risk of water transfer projects [13].
Unlike other pollution accidents, sudden water pollution accidents have varied and complex pollution sources [14,15], while the pollution pathway and degree are unpredictable [16].Also, sudden water pollution accidents not only cause significant economic losses, but also seriously endanger the safety of industrial and domestic water supplies, thus causing social panic [17,18].Therefore, there is much research on how to deal with sudden water pollution accidents in WTPs.Lei et al. (2018) [19] developed and implemented a decision support system that can provide technical support in managing pollution accidents.Xu et al. (2016) [20] proposed a real-time, rapid emergency control mode.Wang et al. (2018) [21] proposed a method for identifying pollution sources in rivers that could quickly control the spread of pollution.Long et al. (2016) [22] proposed a more comprehensive emergency control system.Tang et al. (2015) [23] proposed a set of reasonable risk forecasting models based on an integrated Bayesian Network and water quality model.These systems offer valuable warning and response tools to reduce the impact of accidents.But for accidents that have already occurred, these models and systems are less effective.At this time, a more comprehensive plan library is needed, for which the key step is to select the optimal plan quickly and accurately from a plan library.
In the field of evaluation, decision-making models based on statistical data of historical events and subjective expert judgment have been widely used.In the decision-making model, the rationality of the emergency plan for sudden water pollution accidents is determined by analyzing the important indicators affecting sudden water pollution accidents through the technique for order preference by similarity to ideal solution (TOPSIS) [24][25][26], fuzzy fault tree analysis [27,28], the vulnerability model [29], multi-criteria analysis [30], and other decision methods.These methods are used in different situations.Compared with other multi-criteria decision analysis methods, TOPSIS has some advantages, such as the most adequate use of raw data information and no strict restrictions on the number of indicators [31].Due to its practicability and scientific nature, TOPSIS has been verified and applied by experts in many fields [32][33][34].Unfortunately, only a handful of experts have applied TOPSIS to assess emergency plans and generate decision support options.Therefore, we tried to apply TOPSIS to calculate the relative closeness of the ideal solution.
According to the wooden bucket theory [35,36], the merits of an organization depend not only on the overall conditions, but also on every element of the organization.Just like the two companies participating in an election, the overall score can be the same, but the first company is more average in all aspects, while another company has obvious disadvantages.In this case, the first company will win.Therefore, the unevenness between the indicators of the evaluation object will influence emergency plan assessment and warrants consideration.Accordingly, this research utilizes the Shannon entropy method (1948) to evaluate the inhomogeneity among the attributes of one evaluation object [37].
For the relationship between multiple systems, the coordinated development degree model (CDDM) has been investigated [22,38].In this paper, relative closeness and the inhomogeneity are regarded as two systems that influence the optimality of the emergency plan and the CDDM is used to assess their relationship.
The emergency plan assessments for global water transfer projects as well for the projects of China have been gradually increasing in quantity.Nevertheless, the existing emergency plan evaluation methods only focus on the final scores of different evaluation objects, and do not reflect the different quality between these evaluation objects.Therefore, this study develops an integrated assessment method of emergency plans for sudden water pollution accidents based on an improved TOPSIS and Shannon entropy method, which can effectively compensate for this drawback.
The method creates a combination of multi-indices system and the inhomogeneity of multiple indicators to assess the emergency plan.Our work aims to (1) establish a multi-indices system of risk, timelines, economy and feasibility; (2) determine the weight of each indicator based on AHP and entropy method; (3) develop an integrated emergency plan assessment method based on improved TOPSIS and Shannon entropy; (4) verify the rationality of the method with a case study; (5) compare the integrated method with other methods; (6) discuss the advantages and deficiencies of the proposed method.

Integrated Assessment Method of Emergency Plan
Our study proposes a multiple emergency plan assessment method, which divided into five phases.Firstly, the evaluation indicator system for emergency plan is established by literature research and statistics.Secondly, the weight of each indicator is determined based on AHP and the entropy method.Thirdly, the relative closeness to the ideal solution is calculated based on TOPSIS.Then, the inhomogeneity of the different indicators is assessed based on Shannon entropy.Finally, the emergency plan is comprehensively evaluated using the coordinated development degree model.The framework of the emergency plan assessment method is shown in Figure 1.
The method creates a combination of multi-indices system and the inhomogeneity of multiple indicators to assess the emergency plan.Our work aims to (1) establish a multi-indices system of risk, timelines, economy and feasibility; (2) determine the weight of each indicator based on AHP and entropy method; (3) develop an integrated emergency plan assessment method based on improved TOPSIS and Shannon entropy; (4) verify the rationality of the method with a case study; (5) compare the integrated method with other methods; (6) discuss the advantages and deficiencies of the proposed method.

Integrated Assessment Method of Emergency Plan
Our study proposes a multiple emergency plan assessment method, which divided into five phases.Firstly, the evaluation indicator system for emergency plan is established by literature research and statistics.Secondly, the weight of each indicator is determined based on AHP and the entropy method.Thirdly, the relative closeness to the ideal solution is calculated based on TOPSIS.Then, the inhomogeneity of the different indicators is assessed based on Shannon entropy.Finally, the emergency plan is comprehensively evaluated using the coordinated development degree model.The framework of the emergency plan assessment method is shown in Figure 1.

Establishment of Indicator System
Based on the literature review for emergency plan indictors [39,40] and the concept of decision making, which is one of the most important part of modern decision science, the indicator system for emergency plan assessment, is established.The hierarchy structure of the emergency plan assessment is shown in Figure 2, which includes indicators in resistance risk (B1), timeliness (B2), economy (B3) and feasibility (B4).Resistance risk mainly refers to the ability to respond to emergencies, and resistance risk indices here are focused on comprehensiveness of prevention measures (C11), accuracy of potential risk assessment (C12) and the operator's emergency response to

Establishment of Indicator System
Based on the literature review for emergency plan indictors [39,40] and the concept of decision making, which is one of the most important part of modern decision science, the indicator system for emergency plan assessment, is established.The hierarchy structure of the emergency plan assessment is shown in Figure 2, which includes indicators in resistance risk (B 1 ), timeliness (B 2 ), economy (B 3 ) and feasibility (B 4 ).Resistance risk mainly refers to the ability to respond to emergencies, and resistance risk indices here are focused on comprehensiveness of prevention measures (C 11 ), accuracy of potential risk assessment (C 12 ) and the operator's emergency response to danger (C 13 ).Timeliness mainly refers to the speed of the response to an emergency, and here are focused on the speed of the emergency starting (C 21 ), speed of arrival of personnel (C 22 ) and supplies and speed of restoring water delivery (C 23 ).Economy refers to the implementation effect and resource consumption.The main concerns are emergency funds (C 31 ), utilization of emergency resources (C 32 ) and implementation effects (C 33 ).Feasibility mainly refers to the difficulty degree of implementing emergency plans, and this paper mainly considers the rationality and scientific aspects of the plan (C 41 ), controllability of resource allocation (C 42 ), contingency response to accident evolution (C 43 ) and coordination between departments (C 44 ).The factors affecting each indicator are shown in Table 1.When assessing the risk of pollution accidents, we need to consider the types of pollutants, origin of pollution, levels of pollutants, extent of the pollution, location of accidents, economic losses, environmental impacts, and engineering damage.danger (C13).Timeliness mainly refers to the speed of the response to an emergency, and here are focused on the speed of the emergency starting (C21), speed of arrival of personnel (C22) and supplies and speed of restoring water delivery (C23).Economy refers to the implementation effect and resource consumption.The main concerns are emergency funds (C31), utilization of emergency resources (C32) and implementation effects (C33).Feasibility mainly refers to the difficulty degree of implementing emergency plans, and this paper mainly considers the rationality and scientific aspects of the plan (C41), controllability of resource allocation (C42), contingency response to accident evolution (C43) and coordination between departments (C44).The factors affecting each indicator are shown in Table 1.When assessing the risk of pollution accidents, we need to consider the types of pollutants, origin of pollution, levels of pollutants, extent of the pollution, location of accidents, economic losses, environmental impacts, and engineering damage.

Calculation of the Weights for the Indicator
At present, there are many methods to determine the indicator weight, which can be roughly divided into subjective weighting assessment, objective weighting assessment and combined integration weighting assessment.The subjective weighting method refers to the method of comparing, assigning and calculating the weight of each indicator according to the knowledge and experience or preferences of the decision makers, among which AHP [41,42] and Delphi [43,44] are commonly chosen.The objective weighting method is a method for determining the weight of each indicator based on the difference in objective data of each program evaluation indicator value.It mainly includes principal component analysis [45,46], the mean squared error method [47,48], and the entropy method [49,50], and so on.No matter which method is selected, the weight obtained will have a certain deviation, so in order to compensate for this defect, two weight methods need to be combined, which is the combined integration weighting method.To this end, the AHP and entropy weight method are combined to determine the indicator weights [25].The main steps of the combined integration weighting method are summarized as follows: (1) Constructing a raw data matrix of m objects and (2) Calculating the weight p ij of the jth indicator: (3) The entropy e j of the jth indicator is defined as follows: (4) Calculating the deviation degree g j of the jth indicator: (5) The indicator of entropy weight is defined as follows: (6) We calculate the subjective weight λ j by using the AHP method.The detailed steps of AHP method are listed in Long and Xu (2016) [22]; (7) We calculate the final weight ω j as follows: In order to make the distribution of objective weight and subjective weight more reasonable, the distance function is introduced to combine the subjective weight with the objective weight.The distance function between the subjective weight λ j and the objective weight µ j is: Sustainability 2019, 11, 510 6 of 18 The subjective weight and the objective weight distribution coefficient are α and β, respectively, and the calculation is performed using Equations ( 7) and (8).
The combined integration weighting ω j of the jth indicator is calculated as follows:

Calculate the Relative Closeness to the Ideal Solution Based on TOPSIS
Compared with other multi-criteria decision analysis methods, TOPSIS has some advantages, such as the most adequate use of raw data information and no strict restrictions on the number of indicators [31].TOPSIS is a ranking method that attempts to choose alternatives that simultaneously have the shortest distance from the positive ideal solution and the farthest distance from the negative ideal solution [51].TOPSIS makes full use of the attribute information and provides a cardinal ranking of the alternatives [52].The TOPSIS method was first proposed by C. L. Hwang and K. Yoon in 1981.The TOPSIS method is based on the closeness of a limited number of evaluation objects to the ideal solution [52].There are two ideal solutions, one is the positive ideal solution, and the other is the negative ideal solution.The best object should be the closest to the positive ideal solution and the farthest from the negative ideal solution [51].In this study, the relative closeness to the ideal solution is calculated based on TOPSIS.
The calculation process of the TOPSIS method is as follows [34]: Step 1: Constructing the original decision matrix: In Figure 2, the determination of each indicator is related to a number of factors, and the score C i of each indicator is determined in the form of expert consultation, as shown in Equation (10).
Where S k is the expert score of the k th factor which affect the ith indicator.
According to the calculated value, the number of evaluation objects is denoted by m.At the same time, the number of evaluation indicators is recorded as n.The performance matrix can be presented as Equation (1).
Step 2: Normalize the decision matrix: Because the evaluation indicator value types are different, the values can be accurate real values, interval numbers, triangular fuzzy numbers, etc.Therefore, the evaluation indicators in this paper are uniformly normalized by the following Equation.
Step 3: Calculate the weighted normalized decision matrix: where ω j is the combined integration weighting of the jth indicator.
Step 4: Calculate the distance from each alternative to the positive ideal solution and the negative ideal solution: where d i + denotes the distance between the ith alternative and the positive ideal solution, while d i

−
denotes the distance between the ith alternative and the negative ideal solution, where c j Step 5: Calculate the relative closeness to the ideal solution:

Inhomogeneity Evaluation Based on the Shannon Entropy Method
The optimal solution determined by TOPSIS is that it does not consider the inhomogeneity of evaluation indicator.The inhomogeneity of the evaluation indicators and the non-hierarchical value of the indicator weights will affect the determination of the optimal emergency plan.
In this study, the Shannon entropy method is employed to calculate the inhomogeneity of all indicators.In 1948, Shannon proposed the concept of "information entropy" to solve the problem of quantitative measurement of information.The amount of information in a piece of information is directly related to its uncertainty.For example, if we want to figure out a very uncertain thing, or something we don't know, we need to know a lot of information.On the contrary, if we have a lot of knowledge about something, we don't need much information to figure it out.So, from this perspective, we think the measure of the amount of information is equal to the amount of uncertainty.As proposed by Shannon, Shannon entropy is widely accepted [53,54].According to Shannon entropy, the larger the discrete degree of an indicator is, the larger influence it may have in the decision-making process, and thus, the smaller the information entropy, the reliability of the scheme, and vice versa.The calculation process is presented as follows: where H i is the entropy of the evaluation object, m is the number of evaluation objects, n is the number of attributes indicators.

Comprehensive Evaluation Based on Coordinated Development Degree Model
Based on the above relative closeness and inhomogeneity evaluation results, we try to apply a coordinated development degree model (CDDM) for emergency plan assessment, focusing on the relationship between relative closeness and inhomogeneity of emergency plans.The calculation process is presented as follows: Step 1: Data normalization to (0,1) Step 2: Coordinated development degree calculation In order to quantify the relationship between relative closeness and inhomogeneity, the coordinated development degree is used to characterize it.
where M i is the coordination degree of the ith emergency plan; T is the comprehensive evaluation indicator of the system, here T = 13; D is the coordinated development degree of the ith emergency plan.
Step 3: Determining the optimal emergency plan According to the above method, the D of each emergency plan is determined, and then the plan is sorted.The larger D is, the better the corresponding emergency plan is, that is:

Scenario Design
The SNWT Project, one of the most important and strategic projects in China, will be beneficial to suppress and reduce groundwater withdrawal.The decision support system for the Middle Route has been deployed by the Middle Route Construction and Administration Bureau since 2010 [19].Emergency response plans for 388 pollutants have been stored in this system.In order to test the rationality of the method, we extracted a scenario in the system for testing.A schematic diagram of the test area is shown in Figure 3.The detailed descriptions of the designed scenarios are shown in Table 2.According to the rapid prediction formulas proposed by Xu [20], the time for the pollutants to arrive at the Qi River Gate and the Shier-li River Gate was 20 and 180 min respectively.Given that the designed water conveyance flow of the Taocha Dam is 350 m 3 /s, four emergency plans were proposed, as shown in Table 3.
Sustainability 2019, 11, x 8 of 17 Step 3: Determining the optimal emergency plan According to the above method, the D of each emergency plan is determined, and then the plan is sorted.The larger D is, the better the corresponding emergency plan is, that is:

Scenario Design
The SNWT Project, one of the most important and strategic projects in China, will be beneficial to suppress and reduce groundwater withdrawal.The decision support system for the Middle Route has been deployed by the Middle Route Construction and Administration Bureau since 2010 [19].Emergency response plans for 388 pollutants have been stored in this system.In order to test the rationality of the method, we extracted a scenario in the system for testing.A schematic diagram of the test area is shown in Figure 3.The detailed descriptions of the designed scenarios are shown in Table 2.According to the rapid prediction formulas proposed by Xu [20], the time for the pollutants to arrive at the Qi River Gate and the Shier-li River Gate was 20 and 180 min respectively.Given that the designed water conveyance flow of the Taocha Dam is 350 m 3 /s, four emergency plans were proposed, as shown in Table 3.The water supply is normally supplied to each water inlet.The water level fluctuation before the sluice gates is 30 cm.

EP 3
The pollutants are confined to the following two channels.The order of closing the gates is: First, the Yanling gate was closed in 15 min.Secondly, after 30 min, the Qi River gate was closed and the closing time was 60 min.Finally, after 120 min, the Shier-li River gate was closed and the gate was closed for 60 min.At the same time, the Weihe River sluice gate is activated.
The water supply is normally supplied to each water inlet.Each check gate is closed synchronously for 15 min, and each backwater gate is not activated.

EP 4
The water supply is normally supplied to each water inlet.The water level fluctuation before the sluice gates is 30 cm.
Note: EP stands for emergency plan.

Integrated Assessment Method of Emergency Plan
According to the proposed integrated assessment method of emergency plan, the four kinds of emergency plans discussed previously are evaluated, and the optimal plan is selected.The calculation process and algorithm are shown in Figure 4 and the methodology steps together with the results are explained as follows.
Step 1: Determine the indicator value.
Based on Figure 1 and the principles of Table 1 and on experts' experience and knowledge, the expert score sheets for the four emergency plans are as shown in Tables S1-S4 (Electronic Supplementary Material).According to Equation (10), the data in Tables S1-S4 are calculated to obtain indicator values, as shown in Table S5 (Electronic Supplementary Material).The average value of the five expert scores was taken as the final value of the indicator, and the results are shown in Table 4.
Step 2: Calculation of objective weight µ j for the indicator The data in Table 4 was normalized and then the normalized matrix Y are obtained as shown in Equation (S1) (Electronic Supplementary Material).The entropy of each indicator is obtained according to Equations (2) to (3).The entropy e j and entropy weight of the indicators µ j are obtained according to Equations (4) to (5).The results are shown in Table 5.
Step 3: Calculation of subjective weight λ j for the indicator We calculate the subjective weight λ j by using the AHP method for taking the subjective attributes of the data into consideration.After expert consultation, the judgment matrices of four kinds of emergency plans are obtained.Take the emergency plan one as an example for detailed introduction.Based on the principles of Table 1 and on experts' experience and knowledge, the calculations of the local weights of events B 1 -B 4 , resistance risk event C 11 -C 13 , timeliness event C 21 -C 23 , economy events C 31 -C 33 , and feasibility events C 41 -C 44 were shown as the values in Tables S6-S10 (in Electronic Supplementary Material).The local weight and global weight of each evaluation indicator could be achieved according to Equations ( 1) to ( 4) in Long's paper (2016).
Step 4: Calculation of integration weighting ω j The water supply is normally supplied to each water inlet.The water level fluctuation before the sluice gates is 30 cm.

EP 3
The pollutants are confined to the following two channels.The order of closing the gates is: First, the Yanling gate was closed in 15 min.Secondly, after 30 min, the Qi River gate was closed and the closing time was 60 min.Finally, after 120 min, the Shier-li River gate was closed and the gate was closed for 60 min.At the same time, the Weihe River sluice gate is activated.
The water supply is normally supplied to each water inlet.
Each check gate is closed synchronously for 15 min, and each backwater gate is not activated.

EP 4
The water supply is normally supplied to each water inlet.
The water level fluctuation before the sluice gates is 30 cm.
Note: EP stands for emergency plan.

Integrated Assessment Method of Emergency Plan
According to the proposed integrated assessment method of emergency plan, the four kinds of emergency plans discussed previously are evaluated, and the optimal plan is selected.The calculation process and algorithm are shown in Figure 4 and the methodology steps together with the results are explained as follows.The integration weighting ω j of each indicator is calculated according to Equations ( 6) to ( 9).The results are shown in Table 6.Step 5: Building a weighted normalized decision matrix According to Equations ( 11) to ( 12), a weighted normalized decision matrix of the indicator is constructed as shown in Equation (S2) (in Electronic Supplementary Material).
Step 6: Calculate the relative closeness to the ideal solution According to Equations ( 13) to ( 14), the distance from each alternative to the positive ideal solution and the negative ideal solution are calculated as shown in Equations (S3) and (S4) (in Electronic Supplementary Material).Then, the relative closeness to the ideal solution of each emergency plan are calculated based on Equation (15), the result is shown in Equation (22).

V =
0.4305 0.3507 0.3756 0.3551 Step 7: Calculate the inhomogeneity The larger the discrete degree of an indicator is, the larger influence it may have in the decision-making process, and thus, the smaller the information entropy, the reliability of the scheme, and vice versa.The Shannon entropy of each emergency plan is calculated according to Equations ( 16) to (17), and the result is as shown in Equation (23).
Step 8: Determining the optimal emergency plan According to Equation (18), data of relative closeness and Shannon entropy are normalized as shown in Equations (S5) and (S6) (in Electronic Supplementary Material).And then the coordination degree and coordinated development degree are calculated as shown in Equation (S7) (in Electronic Supplementary Material) and Equation (24).

D =
3.5896 3.5771 3.6055 3.5772 (24) The larger D is, the better the corresponding emergency plan is.Therefore, according to Equation ( 21), the optimal emergency plan is the emergency plan 3.

Discussion
This work presents an approach for assessing the emergency plan by combining improved TOPSIS, Shannon entropy and Coordinated development degree model.The Shannon entropy method was used to analyze different types of index values.TOPSIS is used to calculate the relative closeness to the ideal solution.Coordinated development degree model is applied to express the relationship between the relative closeness and inhomogeneity of the emergency plan, but their potential in emergency plan assessment has rarely been investigated in the literature.
The assessment results of two kinds of emergency plan by different methods such as TOPSIS and the integrated emergency plan are presented in Table 7.We find out the differences between the two methods in the emergency plans and the results are shown in Table 8.EP 2 and EP 4 gain the same ranking in both assessment methods.At the same time, the rankings of EP 1 and EP 3 are reversed by two methods.The optimal plan will change when we take into account the internal inhomogeneity of the indicators in one evaluated object.The major driving factors of the water transfer project contingency plan can also be determined by the integrated assessment method.The weights of the indicators determined by the subjective weighting method and the objective weighting method are quite different, and the main driving factors cannot be determined.Therefore, we determine the main drivers based on the rankings of the indicators that affect the emergency plan.As shown in Figure 5, it can be seen that the ranking of index weights determined by subjective weighting method such as indicators C21, C32 and C33 are the main influencing factors; while the index weighting order determined by objective weighting method can be seen, indicators C23, C41 and C43 are the main influencing factors.The rankings of indicators determined by objective weighting method and subjective weighting method are not the same.The reasons why the discrepancies occurred may be the incompleteness of the data and the limitations of the evaluation experts.The combined integration weight will effectively make up for these deficiencies, and improve the weight of each index to be more reasonable, in order to show the status of this indicator in the evaluation of the whole event more accurately.In the combined integration weight ranking, indicators C33, C31, C21 and C43 are in a prominent position.These rankings indicate that indicator C33 (implementation effect) plays an important and decisive role in the evaluation of emergency plans.The main factors affecting indicator C33 are water quality level after disposal, the quantity of abandoned water and damage degree caused in the water transfer project.The level of water quality after disposal indicates the degree of pollution impact, and the quantity of abandoned water and damage degree caused in the water transfer project represent the economic losses.In the realistic practice, the main considerations for decision makers when they need to determine the emergency plan are the effect of planned disposal and the economic loss that may be caused.Therefore, the combined integration weight method proposed in this paper is not only reasonable, but also relatively scientific.
project.The level of water quality after disposal indicates the degree of pollution impact, and the quantity of abandoned water and damage degree caused in the water transfer project represent the economic losses.In the realistic practice, the main considerations for decision makers when they need to determine the emergency plan are the effect of planned disposal and the economic loss that may be caused.Therefore, the combined integration weight method proposed in this paper is not only reasonable, but also relatively scientific.In the process of evaluating the four emergency plans, we consulted five relevant experts.Their evaluation results for each plan are shown in Figure 6 by different colors of the line.It can be seen from the evaluation results that the index fluctuation is small, except for individual indicators.The small fluctuations in the evaluation results of the indicators are mainly due to the fact that the evaluation experts we chosen are both familiar with environmental and water transport engineering, as well as, mainly because the knowledge level and the engineering experience between them are not too much different.Therefore, it is necessary to further improve the expert database in order to deal with other more different situations.Experts are not limited to the related scientific majors, and we In the process of evaluating the four emergency plans, we consulted five relevant experts.Their evaluation results for each plan are shown in Figure 6 by different colors of the line.It can be seen from the evaluation results that the index fluctuation is small, except for individual indicators.The small fluctuations in the evaluation results of the indicators are mainly due to the fact that the evaluation experts we chosen are both familiar with environmental and water transport engineering, as well as, mainly because the knowledge level and the engineering experience between them are not too much different.Therefore, it is necessary to further improve the expert database in order to deal with other more different situations.Experts are not limited to the related scientific majors, and we should consult more staff member in many aspects, such as skilled technicians and managers in relevant projects.Comparing the four figures in Figure 6, it can be seen that the indicator C43 (contingency response to accident evolution) in the emergency plans 1 and 2 is relatively low.The main reason is that these two plans did not consider controlling pollutants in the accident pool as soon as possible, leading to an increase in the diffusion range of pollutants.Moreover, in relative terms, the variability and resilience of emergency plans 1 and 2 are both poor.That is why we did not take these two plans into the consideration of our actual project.
In the emergency plan 4, the value of indicator C23 (supplies and speed of restoring water delivery) is slightly lower.This is mainly because in emergency plan 4, in order to ensure the safety of the water transfer project, the water level before the sluice gates is required to be allowed to increase by 30 cm, resulting in a complicated control of the gates and a long time being required for the water level to stabilize.However, in actual projects, if a major pollution event occurs, the scope of pollution should be reduced at first, then secondly, the safety of the water delivery project should be considered.Therefore, comprehensive consideration shows that contingency plan 4 is not very reasonable.
In emergency plan 3, the evaluation shows that each index is relatively uniform, and there is no excessive fluctuation.There is no doubt that the experts have high recognition of emergency plan 3, so emergency plan 3 is considered to be reasonable.
In summary, although the emergency plan 1 and 2 are relatively simple to operate, (1) the gate closing time is too fast, causing great damage to the channel; (2) the pollutants are disposed in the channel, and the cost is relatively high; (3) the polluted water after the disposal treatment still below the standard level and we cannot abandon the water directly, (4) the flexibility of the program is weak.While emergency plan 3 looks rather cumbersome, the gates on the Middle Route of the South-to-North Water Transfer Project are all automated, so they're simple to operate in an actual project.The emergency plan 3 is the recommended solution in the decision support system [19].Therefore, the integrated assessment method proposed in this paper is relatively reasonable.channel, and the cost is relatively high; (3) the polluted water after the disposal treatment still below the standard level and we cannot abandon the water directly, (4) the flexibility of the program is weak.While emergency plan 3 looks rather cumbersome, the gates on the Middle Route of the Southto-North Water Transfer Project are all automated, so they're simple to operate in an actual project.The emergency plan 3 is the recommended solution in the decision support system [19].Therefore, the integrated assessment method proposed in this paper is relatively reasonable.In this study, we introduced an improved integrated assessment method that can evaluate the emergency plan more objectively and reasonably.The main advantages of the proposed method are as follows: (1) The relative closeness to the ideal solution using TOPSIS and the inhomogeneity of all indicators based on the Shannon entropy method are synthetically considered to make the evaluation results closer to the actual situation.(2) The entropy method that considers the objective attributes to the data and the AHP method that takes into account the subjective characteristics of the data are combined, and the final determined index weight is more reasonable.
However, the proposed method also has certain deficiencies.(1) The optimal emergency plan is still limited by the experts' experience and knowledge to some extent.Therefore, we will further improve the expert database to reduce the impact of professional knowledge and experience.
(2) Although the indicator system can reflect the relationship between human society and the environment, it was imperfectly scientific and reasonable.Therefore, we will explore the multi-faceted factors affecting the effectiveness of the emergency plan and establish a more perfect indicator system in further research.

Conclusions
The study proposed a novel method for assessing emergency plans based on improved TOPSIS, Shannon entropy and the Coordinated development degree model.The main conclusions are as follows: (1) Thirteen assessment indices were selected to establish the indicator system based on the principle of being scientific, systematic, comprehensive, hierarchical, regional and dynamic.The indicator system was analyzed layer by layer, which reflects the relationship between human beings, society, the environment and transfer projects.
(2) If only the subjective method or objective method is used to determine the index weight, the result is somewhat unreasonable.This study proposes a comprehensive weighting method based on the entropy method and AHP method.The key of the comprehensive weighting method is to

Figure 1 .
Figure 1.The framework of the emergency plan assessment method.

Figure 3 .
Figure 3. Schematic of the test area.

Figure 3 .
Figure 3. Schematic of the test area.

Figure 4 .
Figure 4. Calculation process and algorithm.Step 1: Determine the indicator value.Figure 4. Calculation process and algorithm.

Figure 4 .
Figure 4. Calculation process and algorithm.Step 1: Determine the indicator value.Figure 4. Calculation process and algorithm.

Figure 5 .
Figure 5. Analysis of the main driving factors of integrated emergency plan.

Figure 5 .
Figure 5. Analysis of the main driving factors of integrated emergency plan.

Figure 6 .
Figure 6.Evaluation results of various emergency plan indicators.Figure 6. Evaluation results of various emergency plan indicators.

Figure 6 .
Figure 6.Evaluation results of various emergency plan indicators.Figure 6. Evaluation results of various emergency plan indicators.

Table 1 .
Description of indicator.
Figure 2. Emergency plan evaluation indicator system.

Table 1 .
Description of indicator.
33 Water quality level after disposal, The number of abandoned waters, Damage degree of water transfer project C 41 Technical maturity, Operability, Rationality, Scientific basis, Degree of acceptance C 42 Number of materials and personnel, Allocation rights C 43 Variability, Resilience C 44 Co-movement, Coordination, Communication ability Sustainability 2019, 11, 510 5 of 18

Table 2 .
Scenario set for sudden water pollution accident in the study area.

Table 2 .
Scenario set for sudden water pollution accident in the study area.

Table 5 .
The entropy of each indicator.

Table 6 .
The integration weighting ω j of each indicator.

Table 7 .
Assessment results of different methods.

Table 8 .
Ranking differences of emergency plans between the two methods.