#### Abstract

Most cities in the world have extensive underground facilities, including public transport and commercial facilities, such as shopping malls, subways, and parks. Safety of underground space has been an important aspect of urban safety. Due to the influences of global climate change and human activities, waterlogging disasters in cities are becoming increasingly serious and underground facilities easily suffer from waterlogging disasters. Moreover, waterlogging disasters in cities can cause different degrees of damage to construction period and operating period of urban subway stations. By using the optimal combination weighting method combining subjective and objective weighting, this study assigned weights to evaluation factors. Based on this, a fuzzy synthetic evaluation model for waterlogging risk in the construction and operating periods of urban subway stations was established. Furthermore, the model was applied and verified to be effective in Chengdu Metro Line 4, Sichuan province, China, and the evaluation results coincided with the actual situation. The model in the study provides a new idea for evaluating waterlogging risks in urban subway stations, and these results can offer important information for local government to strengthen management on waterlogging risks.

#### 1. Introduction

In recent years, due to urbanization effect, climate warming, and sea-level rise, meteorological and hydrological disasters, such as typhoon, rainstorm, and extreme rainfall, frequently occur, which exerts great adverse impacts on construction and operation of urban subways. Accidents caused by floods in urban underground rail transit are also frequently reported. For example, in October 1996, heavy rainfall in Boston, USA, flooded the subway system between Kenmore Square and Hynes Convention Center so that Kenmore square was flooded with water of 5.3 × 10^{4} m^{3} at a depth of more than 7 m [1]. In July, 2001, Seoul in South Korea was hit by an extraordinary rainstorm, flooding the subways, which caused forty-nine deaths, including 21 deaths from electrocution. Moreover, on September 6, 2001, rainstorm and flood induced by Typhoon Nari flooded 18 subway stations in Taipei, China. Together with water accumulated in running tunnels, the total volume of accumulated water was about 3 × 10^{5} t so that the subway system was paralyzed, which brought profound lessons. On April 1, 2004, in the construction site of Lijiao Station of Guangzhou Metro Line 3 in China, the enclosure structure of underground continuous walls suddenly collapsed due to continuous heavy rain. A building only 30 m from the construction site leaned to the south for 60°, and the foundation pit by the side of the building collapsed, forming a hollow place with the area of more than 1,000 m^{2}, which caused great losses. From June 30 to July 6, 2016, weekly rainfall in Wuhan (Hubei province, China) hit a record high, and heavy rain flooded 206 places in urban areas and forced to shut down three subway stations [2]. The subway stations became flood discharge areas.

For the underlying reason, urban flooding is the result of urbanization. Owing to the fact that many rivers are reduced in size with city development, water level rises even though rainfall is same. The constant increase of impermeable area in cities relatively reduces time for flood reaching the downstream cities and results in quicker flooding. Under the circumstances, if there is heavy rain for a long time, drainage network in most cities cannot discharge the flood in time, so urban underground rail transit has potential hazards of being flooded by accumulation of rainstorm water. Groundwater can also seep into engineering linings surrounded by saturated soil for a long time even if the inlet is not flooded. In serious cases, the structure is damaged and ground subsidence happens, which affects the safety of surface buildings in the vicinity [3]. Moreover, most serious collapse accidents tend to happen in underground rail transit projects in the construction period, resulting in serious casualties, property loss, and unpredictable indirect loss. Therefore, evaluation of waterlogging risk and drainage in design phase, construction period, and operating period of urban subways is particularly important.

Although previous studies have made some achievements, there are still some deficiencies in the risk evaluation of urban subways [4–8]. The indexes of risk evaluation are determined qualitatively or based on semiquantitative experience. In fact, qualitative evaluation can effectively guide quantitative evaluation, and the latter can make the former more accurate. Therefore, waterlogging prevention in urban subways was used as research object. Based on this, a risk evaluation method based on optimal combination weighting of analytic hierarchy process (AHP) and entropy evaluation method and the fuzzy synthetic evaluation model was proposed by referring to risk evaluation methods in other fields, such as economics, aviation, and geology [9, 10]. Furthermore, the method was verified to be feasible and effective by applying it in Chengdu Metro Line 4, thus improving the current situations that the traditional methods cannot be effectively utilized in engineering practices.

#### 2. Materials and Methods

According to the procedures, i.e., risk identification, risk assessment, risk evaluation, and risk countermeasures, the risk research was carried out following the above steps, along with case verification [11].

##### 2.1. Risk Identification

Risk identification is the first step of risk evaluation. Through certain methods, the step reveals influences of various uncertain factors on the whole system, and the identified influence factors should be able to comprehensively reflect all characteristics in the current situation of waterlogging prevention in urban subways. These characteristics should not only show engineering factors, such as safety and reliability that are concerned in terms of waterlogging prevention, but also reflect nonengineering factors, such as supposed effects on surrounding buildings, transportation, and human society after waterlogging disasters occur. In view of the presence of many factors influencing safety of preventing flood and draining waterlogging in the subways, this study determined risk factors for preventing flood and draining waterlogging in the urban subways through methods, including investigation, consultation, discussion, negotiation, and reference to other literatures. In addition, main factors affecting flood and waterlogging prevention in the urban subway were analyzed and summarized, and secondary factors with slight influences were eliminated. In this way, the list of influence factors for risk evaluation of the urban subway was formed, as shown in Table 1 [12, 13].

##### 2.2. Risk Assessment

After determining the factors influencing waterlogging risks, it needs to assess influence degrees of each factor. At present, scholars in the world have conducted a few research studies on evaluation of waterlogging risk. By referring to the existing research results in other fields, risk assessment methods mainly include AHP, matter-element model, and evidence theory [14–19]. The reliability and correctness of the evaluation results obtained through all methods are directly related to whether factor weights can be scientifically and reasonably determined, while most methods adopt the single subjective or objective weighting method in determining weights of evaluation indexes. The subjective weighting method is mainly based on the knowledge and experience of estimators for judgment, so it can make full use of the knowledge, experience, and prediction ability of risk estimators, but it is subjective and cannot make accurate judgment. These kinds of methods mainly include the AHP, line-relative method, and Delphi method. The objective weighting methods mainly refer to the method of determining corresponding index weights through mathematical methods by fully exploring the information contained in the original data of attributes. The estimated value is not affected by subjective factors, while sometimes the weight value is not consistent with the reality. The methods include the entropy evaluation method, principal component analysis, and gray correlation method. Therefore, the two types of methods for determining weights of assessment indexes have their own limitations.

This study attempted to estimate weights by using the optimal combination weighting method based on the subjective AHP and objective entropy weight method. The optimal combination weighting method is a method that fully uses objective information of indexes to determine weights of evaluation indexes based on the subjective opinions of estimators. By combining the subjective weighting method with objective one, it takes preferences of decision makers to attributes into account and reduces subjectivity of subjective weighting. In this way, unification of objectivity and subjectivity is reached when determining weights, thus reducing errors between the determined weight and the possible weight of the indexes so that evaluation results are more realistic and reasonable.

###### 2.2.1. Analytic Hierarchy Process (AHP)

AHP, which was firstly proposed by Professor Saaty in 1973, a famous mathematician in the United States, is a method combining qualitative weighting and quantitative weighting of risk factors [20]. Being a multicriteria decision-making method, AHP helps solving complex problems with multiple conflicting and subjective criteria, such as projects ranking and location or investment selection [21]. The basic idea of AHP is to decompose a complex problem into different levels and factors and construct judgment matrix by comparing the importance of each hierarchy and factor. Moreover, the maximum eigenvalue *λ*_{max} and corresponding eigenvector *W* are obtained through matrix calculation. After normalization processing, weights of each hierarchy of indexes can be obtained. The specific steps are shown as follows.

*(1) Establishing the Structure Model of AHP*. In general, the structural model of AHP is divided into three hierarchies. The top hierarchy is the intended objective to be reached through systematical analysis; the middle hierarchy includes some criteria and subcriteria involved in realizing the general objective on the top hierarchy. Moreover, the bottom hierarchy consists of various measures, decisions, and schemes selected for achieving the objective.

*(2) Constructing Judgment Matrix*. The relative importance of *n* evaluation factors is judged, and the judgment is quantified through 1∼9 scaling method. The meaning of the scaling method is demonstrated in Table 2 [22]. After that, the judgment matrix *A* with dimension of *n* × *n* is built through pairwise comparison, as shown in the following formula:where indicates the importance of factor *i* to factor *j.*

*(3) Calculating the Maximum Eigenvalue λ _{max} and Weight Vector W*. The weight vector is calculated by using the arithmetic mean method:where

*λ*

_{max}and represent the maximum eigenvalue of the judgment matrix

*A*and the weight vector corresponding to

*λ*

_{max}, respectively; , ; and

*n*indicates the number of evaluation factors.

*(4) Consistency Check*. Due to complexity of objective things and fuzzy judgment and comparison of things of human, in the above judgment matrix is likely to have deviation with the actual ratio. Therefore, it is necessary to carry out consistency check on the weights obtained through judgment matrix, and the values that pass the consistency test are reliable. The test method is shown as follows:(1)Calculating consistency index CI:(2)Calculating random consistency rate CR: where indicates the random consistency index and can be obtained in Table 3 [23].(3)Judgment: If , it means complete consistency and the weights obtained through judgment matrix are completely reliable. If , it needs to further test . If , the consistency of the matrix is generally accepted and the influence weights of the factors determined by the judgment matrix on the objective are set reasonably. If , it is generally considered that consistency of the matrix is unacceptable and the judgment matrix needs to be constructed until consistency meets the requirements.

*(5) Total Ranking of Each Hierarchy and Consistency Check*(1)Total ranking of levels: Total ranking of levels is to calculate total weights of each level, that is, to assign a weight to the importance of all factors at a given level relative to those at the highest level. It needs to be carried out and synthesized from top to bottom layer by layer. The steps are presented as follows. It is assumed that the weight vector of *m* elements in the *k*−1th level in relative to the total objective is It is supposed that the weight vector of single ranking of *n* elements in the *k*th hierarchy in relative to the *j*th element in the upper hierarchy (the *k*−1th hierarchy) is Weight of the elements that are not subjected to control of *j* is 0. represents the ordering of element in the *k*−1th hierarchy by elements in the *k*th hierarchy, so the total ranking of the elements in the *k*th hierarchy to the total objective is obtained as follows:(2)Consistency check on the total ranking results: Assuming that and has been calculated by using the *j*th element in the *k*−1th hierarchy as criterion, the comprehensive test index of the *k*th hierarchy is The consistency check is conducted hierarchy by hierarchy. When the consistency rate of the *k*th hierarchy , it is considered that the satisfactory consistency of recursion order hierarchy is shown in the *k*th hierarchy.

###### 2.2.2. Entropy Evaluation Method

In the theory of information, information entropy is a method to measure uncertainty and reflects disorder degree of information. In the formula, the smaller the information entropy, the smaller the disorder degree of a system and the larger the utility value of information; otherwise, contrary results are obtained.

Entropy evaluation method can objectively reveal disorder of information of data. It determines weight of indexes based on the judgment matrix constructed by values of evaluation indexes, which can eliminate subjectivity of each factor weight as far as possible so that evaluation results are practical. The calculation steps are as follows [24, 25].(1)Constructing the judgment matrix of *m* evaluation indexes in *n* schemes:(2)Normalizing judgment matrix to obtain the normalized judgment matrix : where and indicate the most or least satisfied results in different schemes under the same index (the smaller or larger the index is, the more satisfied the obtained results are).(3)According to definition of entropy, for *m* evaluation indexes in *n* schemes, the entropy of the evaluation indexes can be determined as follows: In order to make meaningful, it is generally assumed that when . However, when , is also equal to zero, which is obviously inconsistent with the reality and contradicts with the meaning of entropy. Therefore, needs to be corrected and defined as(4)Calculating entropy weight of evaluation indexes:

###### 2.2.3. Optimal Combination Weighting Method

The weights of evaluation indexes of waterlogging risk were finally determined by using the optimal combination weighting method. The weights determined based on the above steps combine the subjective weighting method with objective one so that the determined weight coefficient can simultaneously reflect subjective and objective information. The specific formula is shown as follows:where and separately represent the subjective and objective weights; *a* and *b* are corresponding undetermined coefficients, respectively; and denotes the combination weight. It is obvious that the key to determining the combination weight is to determine the undetermined coefficient. According to references, the following formulas are used for calculation [26]:

In order to meet normalization conditions, *a*_{1} and *b*_{1} are substituted into formulas (20) and (21) for normalization processing and *a* and *b* are finally obtained and substituted into formula (17), thus obtaining the combination weight.

##### 2.3. Fuzzy Evaluation Model

Based on the determination of weights of evaluation factors by using the optimal combination weighting method, it is considered that safety of waterlogging prevention of the urban subway in the construction and operating periods is affected by many factors. However, some of these factors are certain and some are uncertain. Moreover, some factors are difficult to be measured and index values are hard to be ascertained accurately. Therefore, safety risks of waterlogging prevention of the urban subway are quantitatively evaluated by combining the optimal combination weighting method with the risk evaluation model with fuzzy attributes. The purpose is to make evaluation results more accurate and scientific and provide references for early warning, prevention, and control of waterlogging disasters in the urban subway. The fuzzy comprehensive evaluation method mainly aims at the evaluation system whose index values are difficult to be ascertained accurately [27]. The specific modeling steps are shown in the following sections.

###### 2.3.1. Determining Evaluation Factor Set and Judgment Set

The evaluation index set is defined according to the specific evaluation system, where represents the index of the evaluation system. In risk analysis of waterlogging prevention, it is defined that (*i* indicates the number of indexes).

Possible judgment set of factors is . and *j* denote the level of risks and the number of evaluation levels, respectively.

###### 2.3.2. Determining Membership Function

The determination of membership function can quantitatively describe fuzzy concepts, so it is a key step of fuzzy evaluation. The specific modeling formula refers to the experience of relevant scholars, and the normal membership function is utilized [28].where *x* indicates the actual value and *x*_{i} and *x*_{i+1} separately represent the upper and lower limits of an evaluation factor corresponding to a certain risk level.

###### 2.3.3. Establishing Fuzzy Evaluation Matrix

Fuzzy evaluation matrix *R* shows the mapping from the evaluation factor set to the judgment set . r_{ij} represents the actual value of the *i*th factor to the *j*th rating grade in the evaluation system and can be obtained through formula (22):

###### 2.3.4. Fuzzy Comprehensive Evaluation

Through fuzzy evaluation and calculation on weight index set *W* and fuzzy evaluation matrix *R* obtained through the combination weighting method, fuzzy evaluation vector *B* can be obtained, as shown in the following formula:

Based on the calculated fuzzy evaluation vector *B*, the maximum value of the corresponding grade of membership and the risk judgment level can be judged, thus determining grades of fuzzy comprehensive evaluation on flood prevention in the urban subway station. In other words, based on grade of membership , the grades of risk levels can be determined.

#### 3. Case Study of the Model to Wenjiachang Metro Depot in Phase I of Chengdu Metro Line 4

##### 3.1. Analysis of General Situation of the Researched Region

By taking Wenjiachang Metro Depot in phase I of Chengdu Metro Line 4 as an example, the evaluation of flood risk in the subway was studied by using the established model. For the engineering line in phase I of Chengdu Metro Line 4 with the overall length of about 22.4 km, there are 16 stations in total, with average interval of about 1.32 km. All belong to underground lines and stations. The total investment is estimated to be 1.2562 × 10^{10} yuan. There is a metro depot, namely, Wenjiachang Metro Depot in the phase I project of the line 4 and two main transformer substations.

Wenjiachang Metro Depot, basically located within Lvzhou Park in Chengdu city, Qingyang district, has a large area of grassland and woods as well as densely distributed houses for residents (peasants). Planting soil that is plastic in local areas is mainly found on the surface in the field, and local areas are low lying and wet. In particular, for Supo branch canal and its influence scope, due to long-term immersion, clayey soil is mainly soft plastic to plastic, and some belong to soft to loose and soft soil, showing large porosity ratio and high water content and compressibility. Characteristic measured area with poor mechanical properties is located in the I-level terrace of Minjiang river system, on which it is clayey soil with the thickness of 2∼4 m of Quaternary Holocene alluvial deposit (), and the soil is mainly hard plastic, with plastic soil in local areas. Planting soil is mainly distributed at the surface of 1 m thick, and filling soil is seen in local areas. Below the surface, slightly dense to moderately dense lenticular sand layers sandwiched with boulder and pebble soil with a large particle size are mainly distributed, and the thickness is 7∼11.8 m. The upper Pleistocene series () containing boulder and pebble soil with a large particle size lies in the deeper layers, and the exposed part shows the thickness of 5.5∼11.8 m, and thin layers of lenticular sand with the thickness of about 0.4 m are found in local areas. In addition, ground surface of the metro depot is located in a plain area where there is no municipal drainage system at present, so rain drainage and flood discharge mainly rely on Supo branch canal.

##### 3.2. Determination of Factor Set and Evaluation Set

###### 3.2.1. Evaluation Factor Set

In analysis of waterlogging risk in the subway station, because of different influence degrees of some factors in the construction and operating periods, weight coefficients of each corresponding factor are distinct in the construction and operating periods. In this evaluation, six factors, i.e., terrain *u*_{1}, river *u*_{2}, rain facilities *u*_{3}, feed-water pipe *u*_{4}, groundwater *u*_{5}, and regional importance *u*_{6}, were mainly selected for analysis in the construction period. A reliable structural seepage-proofing system has been formed in the station during the operating period, so influences of groundwater are basically eliminated. For the above reasons, the five factors, i.e., terrain *u*_{1}, river *u*_{2}, rain facilities *u*_{3}, feed-water pipe *u*_{4}, and regional importance *u*_{6}, are used for analysis and evaluation during the operating period. Therefore, evaluation factor sets separately are and in the construction and operating periods.

###### 3.2.2. Establishment of Evaluation Set

The evaluation criteria for each factor are shown in Table 4. According to explanations of each factor in the table, each evaluation grade of each factor is specifically quantified. The quantifying idea is that 0 or 1 is taken in ultimate case under the most unfavorable conditions. Of them, the larger the terrain, groundwater scale, and regional importance are, the more unfavorable it is; the larger the river, rain facilities, and groundwater scale are, the more favorable the obtained results are. Therefore, the evaluation set is , and specific rating grades are demonstrated in Table 5.

##### 3.3. Determination of Weights of Evaluation Factors

By using the AHP and entropy evaluation method, the subjective and objective weights were calculated, and the combination weight was finally obtained by utilizing the optimal combination weighting method.

###### 3.3.1. Determining Subjective Weight Using AHP

In order to determine reasonable subjective weights, questionnaires were issued to experts in relevant fields in many provinces, including metro construction organizations, metro operation companies, metro design institutes, and professors in universities, and 15 effective questionnaires were collected. By using 1∼9 scaling method, the judgment made by experts and scholars on the relative importance of the above evaluation factors is summarized. By comparing and analyzing the contributions of influence factors to waterlogging risks in the station, the judgment matrix using AHP is constructed, as shown in Tables 6 and 7. The subjective weight vectors and in the construction and operating periods are calculated by formulas (3) and (4). Moreover, the corresponding maximum eigenvalues and are separately obtained. In the determination of consistency indexes, indexes of the construction period are CI = 0.0119 and CR *=* 0.0096 < 0.1, while those in the operating period are CI = 0.0010 and CR *=* 0.0009 < 0.1. Therefore, the constructed judgment matrix is reasonable.

###### 3.3.2. Determining Objective Weight Using Entropy Evaluation Method

Based on the calculation steps using the entropy evaluation method, the normalized judgment matrix is firstly determined and weights = (0.1477, 0.3047, 0.2974, 0.0820, 0.1680) and = (0.1007, 0.2077, 0.2027, 0.0559, 0.3186, 0.1145) are calculated by utilizing formulas (13), (15), and (16).

###### 3.3.3. Calculating Weights through Combination Weighting

By substituting the above calculated subjective and objective weights into formulas (18)–(20) of the optimal combination weighting method, the partial coefficients of the final weight of evaluation indexes for waterlogging risk in the subway station are calculated. They are *a* = 0.5385 and *b* = 0.4615 in the construction period as well as *a* = 0.6015 and *b* = 0.3985 in the operating period. Finally, by substituting them into formula (17), the combination weights *W*_{construction period} = (0.1366, 0.1860, 0.1423, 0.0559, 0.3777, 0.1016) and*W*_{operating period} = (0.1132, 0.1757, 0.2756, 0.0870, 0.3486) can be calculated.

##### 3.4. Fuzzy Comprehensive Evaluation

According to the method for building the fuzzy evaluation model, fuzzy evaluation matrix *R* was constructed so as to perform fuzzy evaluation and analysis on waterlogging risks in Wenjiachang Metro Depot. The fuzzy evaluation matrix *R* is shown as follows:

By calculating the fuzzy evaluation matrix *R* by formula (26), the fuzzy evaluation vector *B* is obtained, that is, *B*_{construction period} = (0.00045, 0.00000, 0.23590, 0.22460, 0.40860) and*B*_{operating period} = (0.00070, 0.00000, 0.44340, 0.34290, 0.33100). The specific corresponding relationship of fuzzy evaluation vectors and grades of membership is presented in Table 8.

Based on the maximum membership principle, the grades of waterlogging risks in Wenjiachang Metro Depot in the construction period and operating period are separately V and III.

#### 4. Results and Discussion

By utilizing the optimal combination weighting method combining subjective and objective weighting based on the AHP and entropy weight method to assign weights to evaluation factors, this study established the fuzzy synthetic evaluation model for waterlogging prevention in the construction and operating periods of the subway station in the city. Moreover, the model was applied and validated to be effective in Wenjiachang Metro Depot in phase I of Chengdu Metro Line 4, and the Wenjiachang Metro Depot showed V and III grades of flooding risks in the construction and operating periods. The evaluation results coincided with the actual situation.

In view of evaluation results, the following risk countermeasures are suggested to be used.(1)Engineering measures:(a)During the construction period, the combination of prevention and drainage should be adopted for safety of waterlogging prevention in the construction period in the metro depot. The measures are taken to ensure that the construction site has adequate draining and discharging capacities. For water accumulated on peripheral roads, the corresponding antiscour walls should be set according to grades of waterlogging prevention in the metro depot and temporary drainage ditches with enough water discharge capacity should be arranged in front of the antiscour walls. In the construction, special attention should be paid to protect pipes for sewage and running water in the station and its surrounding so as to prevent pipe-burst accident. In addition, reliable dewatering measures should be taken to reduce groundwater level, thus eliminating threats to the foundation pit.(b)During operating period, a policy dominated by prevention and assisted with drainage should be in line with municipal rainwater facilities in the surrounding of metro depot, and drainage standard for rainwater needs to be improved. Structures like steps in the entrances and exits of the metro depot should be high enough for ensuring safety, and flood gates should be set if necessary. Based on daily maintenance, it is necessary to increase the inspection frequency on flood gates and waterlogging prevention materials so as to ensure normal operation of equipment and allocation of waterlogging prevention materials. In the metro depot, the drainage pump station should have a certain scale to drain internal accumulated water. Furthermore, it needs to notice the treatment on nodes of pump station outlet with municipal pipe network or rivers. In terms of elevation, water sealing, and operation designs for the pump station, effective measures should be taken to ensure that external water (water from rivers and rainwater in pipes) is not poured back into the metro depot through the pump station and internal pipe network.(2)Nonengineering treatment measures:(a)The complete system for preventing waterlogging should be established, and special funds for flood control, disaster reduction, and flood prevention and waterlogging control need to be set up. In addition, social channels, such as flood control insurance should be used as much as possible to avoid disaster risks.

#### 5. Conclusions

This study determined risk factors for preventing flood and draining waterlogging in the urban subways through methods, including investigation, consultation, discussion, negotiation, and reference to other literatures. In addition, the list of influence factors for risk evaluation of the urban subway was formed. These can be applied to more subway constructions, and we need to make appropriate adjustments according to the difference of each subway. We attempted to estimate weights by using the optimal combination weighting method based on the subjective AHP and objective entropy weight method. In this way, we can reduce errors between the determined weight and the possible weight of the indexes so that evaluation results are more realistic and reasonable. This method can also be applied to other research fields.

This study established the fuzzy synthetic evaluation model for waterlogging prevention in the construction and operating periods of the subway station in the city. The model was applied and validated to be effective in Wenjiachang Metro Depot in phase I of Chengdu Metro Line 4. The evaluation results coincided with the actual situation. The model in the study provides a new idea for evaluating waterlogging risks in urban subway stations, and these results can offer important information for local government to strengthen management on waterlogging risks.

We try to provide a new risk evaluation method for security experts, architecture design engineer, construction risk management engineer, government, or other relevant person, but for the constraints, this method cannot be verified comprehensively through more actual projects. Next, further research will be carried out through the whole life cycle of subway projects, not just the construction and operation periods. We hope this model is widely and practically applied, and we can achieve further research results and provide more reference for more risk assessment participants.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This study was supported by the scientific research project funded by the Provincial Department of Education of Sichuan Province, China (18ZA0411) and the Science and Technology Planning Project of Sichuan Province, China (2017JY0514).