Risk Evaluation Based on Variable Fuzzy Sets and Information Diffusion Method

Abstract

Risk evaluation is an effective way to reduce the impacts of natural hazards and it plays an increasingly important role in emergency management. Traditional methods of assessing risks mainly utilize Geographic Information System (GIS) to get risk map, and information diffusion method (IDM) to deal with incomplete data sets. However, there are few papers discuss the uncertainty of integrated hazards and consider dynamic risk under time dimension. The model proposed in this study combines the variable fuzzy set theory with information diffusion method (VFS-IDM) to solve the uncertainness of multiple hazards dynamic risk assessment when data sets are incomplete. This study employs fuzzy set theory (VFS) to calculate the relative membership degree and applies information entropy method (IEM) to obtain the weights of criteria indicators for multiple hazards evaluation. Then applies information diffusion method (IDM) to estimate condition probability distribution and vulnerability curve with the VFS-IEM model results, time data and multiple hazards losses. Then the expected value of multiple hazards dynamic risk can be calculated by using the normal information diffusion estimator so as to improve the accuracy of risk evaluation results.

Keywords

Share and Cite:

Gong, W. , Chen, K. , Huggins, T. and Yang, L. (2020) Risk Evaluation Based on Variable Fuzzy Sets and Information Diffusion Method. Journal of Applied Mathematics and Physics, 8, 821-830. doi: 10.4236/jamp.2020.85064.

1. Introduction

With the global climate change and rapid population growth, many cities suffered extreme natural hazards frequently [1]. The uncontrolled development and unplanned change of land use belong to the highly sensitive areas, where natural hazards cause devastating economic and social losses [2] [3]. These and many other impacts mean it is important to assess the risk of natural hazards and growing economy requires improved risk assessments methods. Risk evaluation is an effective way to reduce the negative impacts of natural hazards and it is an effective way to reduce the effects of hazards which plays an increasingly important role in emergency management and to reduce the losses caused by hazards.

Risk evaluation should consider the occurrence likelihood of a specific event and the severity of the outcome. The result is composed by the probability of all relevant hazards and the severity of losses scenarios. In order to evaluate risk of any hazards, it needs to build a reasonable and effective assessment system. Generally, risk evaluation can be divided into qualitative and quantitative assessment method. The qualitative assessment methods are by means of several software such as geographic information systems (GIS). For example, Chen et al. [4] provided useful detailed information for flood risk management by combining with analytical hierarchy process and GIS in flood risk assessment, and the method can be easily applied to most areas where required data sets are readily available. Aye et al. [5] provided a prototype of an interactive web-GIS tool for the risk evaluation and management evaluation in Central East Moldavian Region, considering the occurrence of floods and earthquakes.

However, there is an unavoidable fact that GIS adopted to risk evaluation can often not be quantified and results are limited by insufficient risk-related data. In order to quantify risk, T. J. Huggins et al. [6] analyzed the cascading disaster risk by using longitudinal data. Ribeiro et al. [7] proposed a probabilistic model, based on a bivariate copula approach using elliptical and Archimedean copulas, to estimate the probability of loss. And Huang [8] put forward the concept about probabilistic risk, which can be quantified as the expected value to predict future risk based on historical data. However, long time series of disaster-related data usually often does not exist, and the information contained in existing data sets is typically incomplete. Then in order to assess the risk of hazards quantitatively, especially when the recorded data sets are incomplete, there must be another method. Many studies have focused on IDM which belongs to fuzzy sets theory and is easier to quantify the probabilistic risk when the data is insufficient. Li [9] proposed a flood risk assessment model based on information diffusion method to deal with the small sample size and gave one example to demonstrate the model is successful. Xu et al. [10] developed another method for using IDM to quantitatively assess the risk of multiple hazards. Huang [11] gave the definition of integrated probability risk and assess the risk of annual loss by using information diffusion method. Although some papers have given the definition and assessed the integrated hazards risk [12] [13], the randomness and fuzziness which determine the reliable of evaluation result have been ignored.

Considering that integrated risk evaluation contains fuzzy concept with multiple indicators and classes, the variable fuzzy sets theory introduced by Chen [14], which give membership degrees and relative membership function to evaluate the fuzzy concept, has been successfully used in risk evaluation. For example, Carreno et al. [15] applied the fuzzy sets theory into the seismic risk evaluation when the data required to assess risk are not available or are insufficient. Li et al. [16] used the fuzzy comprehensive evaluation method, which proved effective in solving problems of fuzzy boundaries and controlling the effect of monitoring errors on assessment results, to analysis flood risk. Wang et al. [17] proposed an integrated variable fuzzy evaluation model, to overcome the limitations of the traditional evaluations which only used a point value instead of an interval for grading standards, for the assessment of river water quality.

There are many approaches to assess risk including uncertainty theory and incomplete sets. Many papers have proved that these methods can be integrated to improve the accuracy of assessment result [18] [19]. Besides, there are many improved models, such as combined the projection pursuit method optimized by immune evolutionary algorithm with information diffusion method has been used to assess risk of drought [20]. However, there are still some obstacles for dynamically assessing integrated hazard risk and we say dynamically means that time dimensional should be considered in risk evaluation.

In our study, the definition of improved probability risk is identified first, then use VFS to eliminate both randomness and fuzziness in integrated risk evaluation. Based on the processed data, take time dimension into consideration and use the IDM to extract as much useful underlying information as possible from samples to estimate relationships behind the incomplete data sets, thus assess the integrated risk dynamically and improve the accuracy of evaluation result. This paper is organized as follows: In Section 2, we give some definitions of risk evaluation and describe some basic concepts and principles of the evaluation model. The process of how to build our model and how to evaluate the integrated risk dynamically is followed and discussed in detail. The results will be shown in Section 3. Finally, Section 4 discusses the importance of the new model and outlines implications for further work.

2. Basic Concepts

2.1. Improved Probabilistic Risk

According to what Huang have said, the risk could be classified into four categories: pseudo risk, probability risk, fuzzy risk, and uncertainty risk [11]. In the following contents, the definition of probability risk and fuzzy risk will be introduced, then the improved probability risk will also be defined.

The probability risk is a scenario in the future that is associated with some specified adverse incident that can be statistically predicted by using a probability model. There are limitations for determining this kind of risk, due to corresponding events with fuzzy boundaries and the predicative information is incomplete. Then we have the fuzzy risk, which is the scenario in the future associated with some specified adverse incident and that we are able to approximately infer by using fuzzy logic and incomplete information. IDM is used to solve this type of risk. Based on the two preceding definitions, we have the improved probability risk which is a scenario in the future associated with integrated adverse incident characterized by fuzzy sets of variables, and where VFS-IDM is used to dynamically assess integrated risk when the data sample is incomplete. Therefore, the equation for improved probabilistic risk can be defined as:

$Risk=probability\text{}distribution\cdot vulnerability\text{}curve$

In the case of that we can estimate the probability distribution $p\left(x\right)$ of the occurrence of hazards with respect to factor x and the relationship $f\left(x\right)$ between the factor and losses, the risk can be quantified as the expected value of hazards, just as show in Equation (1):

$Risk=\int p\left(x\right)\cdot f\left(x\right)dx.$ (1)

2.2. Variable Fuzzy Sets Theory

For a variable fuzzy set (VFS) U and a random element $u\in U$, we can find two relative membership degree functions (RMDF), denoted by ${\mu }_{A}\left(u\right)$ and ${\mu }_{A}^{c}\left(u\right)$, which express the extent of acceptability (A) and repellency ( ${A}^{c}$ ) respectively. That is, ${\mu }_{A}\left(u\right)+{\mu }_{A}^{c}\left(u\right)=1$ and when the RMDF ${\mu }_{A}\left(u\right)$ is larger than ${\mu }_{A}^{c}\left(u\right)$, the major property of is acceptability, and the minor property is repellency. This means a ratio can be used to represent this relationship in the interval and the expression of RMDFs can be derivate [14].

We define interval ${X}_{0}$ = [a, b] as the attracting sets of VFS U and X is a certain extended interval [c, d] including ${X}_{0}$ on the real axis. For any fixed u, we have ${\mu }_{A}\left(u\right)>{\mu }_{A}^{c}\left(u\right)$, where r stands for the assessment object set. If we have ${\mu }_{A}\left(u\right)>{\mu }_{A}^{c}\left(u\right)$ and M is the balance boundaries of interval ${X}_{0}$, then the following Figure 1 has shown the possible locations of x when M is fixed.

Thanks to the research of how to define the balance boundaries of VFS [21], we can get the important parameter ${M}_{rl}$ by following equation:

${M}_{rl}=\frac{L-l}{L-1}{a}_{rl}+\frac{l-1}{L-1}{b}_{rl}.$ (2)

where l denotes the assessment indicator set, $l=1,2,\dots ,L$. The position of $M=\left({M}_{rl}\right)$ decides the expression of RMDF and satisfies following suppositions: 1) when $l=1$, then ${M}_{r1}={a}_{r1}$ ; 2) when $l=L$, then ${M}_{rL}={b}_{rL}$ ; 3) when

$l=\frac{L+1}{2}$, then ${M}_{rl}=\frac{{a}_{rL}+{b}_{rL}}{2}$. To figure out the expression of RMDFs, define relative difference degree of u to A:

${D}_{A}\left(u\right)={\mu }_{A}\left(u\right)-{\mu }_{A}^{c}\left(u\right)\text{}{D}_{A}\left(u\right)\in \left[-1,1\right]$ (3)

Given the point $M=\left({M}_{rl}\right)$, if x is located at the right of M, the relative

Figure 1. The location between different point x and fix point M, [a, b], [c, d].

difference degree can be replaced by a ratio:

$\left\{\begin{array}{l}{D}_{A}\left(u\right)={\left(\frac{x-a}{M-a}\right)}^{p}\text{}x\in \left[a,M\right]\hfill \\ {D}_{A}\left(u\right)=-{\left(\frac{x-a}{c-a}\right)}^{p}\text{}x\in \left[c,a\right]\hfill \end{array}$ (4)

Combined with the property ${\mu }_{A}\left(u\right)+{\mu }_{A}^{c}\left(u\right)=1$ and Equations (3) and (4), the RMDF is denoted by Equation (5).

$\left\{\begin{array}{l}{\mu }_{A}\left(u\right)=0.5\left[1+{\left(\frac{x-a}{M-a}\right)}^{p}\text{}x\in \left[a,M\right]\right]\hfill \\ {\mu }_{A}\left(u\right)=0.5\left[1-{\left(\frac{x-a}{c-a}\right)}^{p}\text{}x\in \left[c,a\right]\right]\hfill \end{array}.$ (5)

From the above Equations, it can be determined that the RMDF is affected by hyper-parameter p and the position between random point x with parameters a, b, c, d, and M. The following conditions are RMDF should satisfy: 1) when $x\notin \left[c,d\right]$ or $x=c,x=d$, ${\mu }_{A}\left(u\right)=0$ ; 2) when $x=a,x=b$, ${\mu }_{A}\left(u\right)=0.5$ ; 3) when $x=M$, ${\mu }_{A}\left(u\right)=1$. In general, we take $p=1$ such that the Equation (5) become linear functions. For the different location of random point x with respect to the object interval, the significant point is that the calculation of relative membership degree (RMD) can be classified into two types: When the random point x located in the lowest or highest point of the standard interval, the RMD sum of this and adjacent level is equal to 1. When random point x located in the interval of a mid-level, the RMD sum of its adjacent levels is 0.5. The detailed process can be found in Fang [22] and we will apply this result to find out the RMD matrix in case study.

2.3. VFS-IEM to Evaluate the Comprehensive Degree Value

For each sample x, the measured values from rth indicator to lth class can be denoted as relative membership degree matrix by using RMDFs. Because the subjectivity of determining index weight and assessment standard entails fuzzy sets, results are often incompatible and even lead to unreliable conclusions. This study proposed the VFS-IEM model to make the result more reliable.

The degree value from rth object to lth indicator called measured value $U={\left(u\right)}_{rl}$ can be defined. The information entropy method (IEM) is used in the fuzzy comprehensive evaluation [23], to calculate the weight of indicators at each monitoring point over the whole assessment standard. For each element in U, the regularization membership vector ${f}_{rl}$ which represents entropy of the lth indicator is

${f}_{rl}={u}_{rl}/\underset{r=1}{\overset{R}{\sum }}\text{ }\text{ }{u}_{rl}\text{}{h}_{l}=-1/lnn\cdot \underset{r=1}{\overset{R}{\sum }}\left({f}_{rl}ln{f}_{rl}\right)$ (6)

Then the entropy coefficient can be defined as

${\omega }_{l}=\left(1-{h}_{l}\right)/\left(L-\underset{l=1}{\overset{L}{\sum }}\text{ }\text{ }{h}_{l}\right),{\omega }_{l}\in \left[0,1\right]\text{}and\text{}\underset{l=1}{\overset{L}{\sum }}\text{ }\text{ }{\omega }_{l}.$ (7)

where ${\omega }_{l}$ is the entropy coefficient of lth indicator.

To assess the variable fuzzy sets more reliable, we should find out the relative membership degree of u to each indicator in a reasonable way. Let ${\nu }_{A}{\left(u\right)}_{l}$ denote as the relative membership degree of u to lth indicator and it remains to be determined. By the study of Wang [21], we say that the fuzzy set can be assessed by VFS-IEM comprehensive evaluation model, denoted by Equation (8), and each sample had been converted from multiple-dimensional indicators into one-dimensional degree value H.

$\left\{\begin{array}{l}{\nu }_{A}{\left(u\right)}_{l}=1+\frac{{\sum }_{r=1}^{R}\left[{\omega }_{l}\left(1-{\mu }_{A}{\left(u\right)}_{rl}\right)\right]}{{\sum }_{r=1}^{R}\left[{\omega }_{l}\left(1-{\mu }_{A}^{c}{\left(u\right)}_{rl}\right)\right]}\hfill \\ {\omega }_{l}=\left(1-{h}_{l}\right)/\left(L-\underset{l=1}{\overset{L}{\sum }}\text{ }{h}_{l}\right)\hfill \\ {\nu }_{A}^{o}{\left(u\right)}_{l}=\frac{{\nu }_{A}{\left(u\right)}_{l}}{{\sum }_{l=1}^{L}{\nu }_{A}{\left(u\right)}_{l}}\hfill \\ H=\left(1\text{}2\dots L\right)\cdot {\left({\nu }_{A}^{o}{\left(u\right)}_{l}\right)}^{T}\hfill \end{array}.$ (8)

2.4. Information Diffusion Method

In order to assess the risk of hazards quantitatively, especially when the recorded data sets are incomplete, IDM which belongs to fuzzy sets theory is used to extract as much useful underlying information as possible from the samples to estimate relationships behind the incomplete data. We can make full use of the diffusion information given by samples to estimate the probability density of samples or the relationship between sample data without knowing which distribution the samples come from.

The research by Huang [24] [25] have given us many results about IDM which can transform data sample point into a fuzzy set with membership function so as to improve the precision of estimators. Let $X=\left\{{x}_{i}|i=1,2,\dots ,n\right\}$ be one dimensional random sample, $V=\left\{{v}_{j}|j=1,2,\dots ,J\right\}$ is the universal field. $\forall {x}_{i}\in X$, the result of the normal diffusion function for sample ${x}_{i}$ diffuse to monitoring point ${v}_{j}$ calculated by

${\mu }_{\left(1\right)}\left({x}_{i},{v}_{j}\right)=exp\left[-\frac{{\left({x}_{i}-{v}_{j}\right)}^{2}}{2{h}^{2}}\right],\text{}{x}_{i}\in X,{v}_{j}\in V$ (9)

And the diffusion coefficient is calculated by the following formula:

$h=\left\{\begin{array}{l}0.6841\left(b-a\right),\text{ }\text{ }\text{ }for\text{\hspace{0.17em}}\text{ }\text{ }n=5;\hfill \\ 0.5404\left(b-a\right),\text{ }\text{ }\text{ }for\text{\hspace{0.17em}}\text{ }\text{ }n=6;\hfill \\ 0.4482\left(b-a\right),\text{ }\text{ }\text{ }for\text{\hspace{0.17em}}\text{ }\text{ }n=7;\hfill \\ 0.3839\left(b-a\right),\text{ }\text{ }\text{ }for\text{\hspace{0.17em}}\text{ }\text{ }n=8;\hfill \\ 2.6581\left(b-a\right)/\left(n-1\right),\text{ }for\text{\hspace{0.17em}}\text{ }\text{ }n\ge 9.\hfill \end{array}\text{}where\text{\hspace{0.17em}}\text{\hspace{0.17em}}b=\underset{1\le i\le n}{\mathrm{max}}\left\{{x}_{i}\right\},\text{\hspace{0.17em}}a=\underset{1\le i\le n}{\mathrm{min}}\left\{{x}_{i}\right\}.$ (10)

Furthermore, we have the m-dimensional diffusion function of X on V,

$\mu \left({x}_{i},{v}_{j}\right)=\underset{k=1}{\overset{m}{\prod }}{\mu }_{\left(k\right)}\left({x}_{ki},{v}_{kj}\right).$ (11)

The Principle of Information Diffusion have proved that it is useful to estimate the probability density of samples or the relationship between sample data without knowing sample distribution by using normal diffusion function [24]. Such as the probability density can be estimated by the value of $\mu \left({x}_{i},{v}_{j}\right)$ which says that observation ${x}_{i}$ gives information to the monitoring point ${v}_{j}$. According to Equation (9), Let

${Q}_{j}=\underset{i=1}{\overset{n}{\sum }}\text{ }\mu \left({x}_{i},{v}_{j}\right)=\underset{i=1}{\overset{n}{\sum }}\text{ }{p}_{ij}\text{}and\text{}H=\underset{j=1}{\overset{J}{\sum }}\text{ }{Q}_{j}$ (12)

The probablity density function can be denoted by matrix q:

${q}_{j}={Q}_{j}/H,j=1,2,\cdots ,J;\text{}and\text{}q=\left({q}_{j}\right).$ (13)

For estimating the relationship between samples, let sample point with input ${x}_{i}$ and output ${y}_{i}$, the universal field $\left(U,V\right)=\left\{\left({u}_{j},{v}_{k}\right)|j=1,2,\dots ,J;k=1,2,\dots ,K\right\}$ and normal diffusion function ${\mu }_{\left(u,v\right)}$. For the value of ${\mu }_{\left({u}_{j},{v}_{k}\right)}\left({x}_{i},{y}_{i}\right)$, observation $\left({x}_{i},{y}_{i}\right)$ can diffuse information to the monitoring point $\left({u}_{j},{v}_{k}\right)$, and obtain a fuzzy relation ${R}_{f}$. Let ${Q}_{jk}={\sum }_{i=1}^{n}{\mu }_{\left({u}_{j},{v}_{k}\right)}\left({x}_{i},{y}_{i}\right)$, it is defined as:

${s}_{k}=\underset{1\le j\le J}{\text{max}}{Q}_{jk},\text{}k=1,2,\cdots ,K;\text{}{\mu }_{k}\left({u}_{j}\right)=\frac{{Q}_{jk}}{{s}_{k}},\text{}j=1,2,\cdots ,J.$ (14)

Then based on Equation (14), we have ${R}_{f}$ model

${r}_{jk}={\mu }_{k}\left({u}_{j}\right)\text{}and\text{}{R}_{f}={\left({r}_{jk}\right)}_{J×K},\text{}j=1,2,\cdots ,J;\text{}k=1,2,\cdots ,K.$ (15)

And it has shown that the max-min fuzzy composition rule, i.e.,

${\mu }_{y}\left({v}_{k}\right)=\underset{{u}_{j}\in U}{\text{max}}\left\{\text{min}\left\{{\mu }_{x}\left({u}_{j}\right),{r}_{jk}\right\}\right\}\text{}{v}_{k}\in V.$ (16)

can make more accurate inference when the samples are incomplete [25].

3. Dynamic Integrated Probabilistic Risk Evaluation Model

Risk evaluation is an effective way to reduce the impacts of hazards and it plays an increasingly important role in emergency management. However, the fuzziness related to integrated hazards and timeless always be ignored in risk evaluation. By the definition of improved probability risk, a combined model based on variable fuzzy sets and information diffusion method to assess integrated hazards risk dynamically is proposed in this study.

By defining the interval criterion matrix ${I}_{ab}=\left({\left[a,b\right]}_{rl}\right)$ and the variable interval matrix ${I}_{cd}=\left({\left[c,d\right]}_{rl}\right)$, the balance boundaries M and RMD matrix can be calculated. According to the VFS-IEM comprehensive evaluation model, the integrated hazards level value of every sample is denoted by Equation (8) and the result is

$H=\left(1\text{}2\dots L\right)\cdot {\left({\nu }_{A}^{o}{\left(u\right)}_{l}\right)}^{T}$

So, the fuzziness in sample data sets have been eliminated and the processed samples have reflected the level value of integrated hazards.

To evaluate the integrated risk dynamically, the time dimensional should be considered. This study used the processed samples to estimate conditional probability distribution and the vulnerability curve by applying the information diffusion method. The Equations (9) (13) show that the discrete probability distribution ${q}_{jk}$ denotes the hazards occurrence.

${q}_{jk}=\frac{{\sum }_{i=1}^{n}\mu \left({x}_{i};{u}_{j},{v}_{k}\right)}{{\sum }_{j=1}^{J}{\sum }_{k=1}^{K}{\sum }_{i=1}^{n}\mu \left({x}_{i};{u}_{j},{v}_{k}\right)},j=1,2,\cdots ,J\text{}and\text{}k=1,\cdots ,K$ (17)

Then the conditional probability distribution with respect to time dimensional denoted by Equation (18)

$p\left(x;{u}_{j}\right)={p}_{v|{u}_{j}}\left(v|{u}_{j}\right)=\frac{{q}_{jk}}{{\sum }_{k=1}^{K}{q}_{jk}},j=1,\cdots ,J.$ (18)

and the vulnerability curve obtained by Equation (16) and center of gravity method [26]

$f\left(x;{u}_{j}\right)=\frac{{\sum }_{k=1}^{K}{\mu }_{y}\left({v}_{k}\right)\cdot {v}_{k}}{{\sum }_{k=1}^{K}{\mu }_{y}\left({v}_{k}\right)},j=1,2,\cdots ,J.$ (19)

Then the integrated hazards dynamic risk can be quantified as the expected value of conditional probability distribution and vulnerability curve, just as shown in Equation (20).

$Risk=\int p\left(x\right)\cdot f\left(x\right)dx=\underset{j=1}{\overset{J}{\sum }}p\left(x;{u}_{j}\right)\cdot f\left(x;{u}_{j}\right).$ (20)

4. Conclusion and Discussion

Risk evaluation is a very important issue in emergency management, but there are few papers discuss the uncertainty of integrated hazards and consider dynamic risk under time dimension when the data set is incomplete. In our study, the definition of probabilistic risk has been modified and takes time dimension into consideration to introduce the concept of dynamic probabilistic risk. Then employs fuzzy set theory (VFS) to calculate the relative membership degree and applies information entropy method (IEM) to obtain the weights of criteria indicators for integrated hazards evaluation. Based on the results obtained by VFS-IEM model, this paper applies information diffusion method (IDM) to estimate condition probability distribution and vulnerability curve with the time data and integrated hazards losses. Then the dynamic risk can be calculated by using the normal information diffusion estimator so as to improve the accuracy of risk evaluation results. The proposed model highlights: the integrated hazards could be processed by VFS theory and combines with IEM to make the integrated hazards level more reliable; solve the problem of limited information in dynamic risk and improve the accuracy of sample data estimation by converting the sample points into fuzzy sets. In the further study, case study will be the focus and data sets will be collected to test out methods.

Acknowledgements

The authors gratefully acknowledge the funding from the National Natural Science Foundation of China, Project No. 71771113 and the National Key R&D Program of China, 2018YFC0807000.

Conflicts of Interest

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

 [1] Liao, X., Zheng, J.X., Huang, C. and Huang, G. (2018) Approach for Evaluating LID Measure Layout Scenarios Based on Random Forest: Case of Guangzhou China. Water, 10, 894. https://doi.org/10.3390/w10070894 [2] Gigovic, L, Pamucar, D, Bajic, Z., et al. (2017) Application of GIS-Interval Rough AHP Methodology for Flood Hazard Mapping in Urban Areas. Water, 9. https://doi.org/10.3390/w9060360 [3] Li, X., Yang, L., Zheng, H., et al. (2019) City-Level Water-Energy Nexus in Beijing-Tianjin-Hebei Region. Applied Energy, 235, 827-834. https://doi.org/10.1016/j.apenergy.2018.10.097 [4] Chen, Y.-R., Yeh, C.-H. and Yu, B.F. (2011) Integrated Application of the Analytic Hierarchy Process and the Geographic Information System for Flood Risk Assessment and Flood Plain Management in Taiwan. Natural Hazards, 59, 1261-1276. https://doi.org/10.1007/s11069-011-9831-7 [5] Aye, Z.C., Jaboyedoff, M., Derron, M.H., et al. (2016) An Interactive Web-GIS Tool for Risk Analysis: A Case Study in the Fella River Basin, Italy. Natural Hazards and Earth System Sciences, 16, 85-101. https://doi.org/10.5194/nhess-16-85-2016 [6] Huggins, T.J., Chen, K., Gong, W. and Yang, L. (2019) The Razor in the Waterfall: Using Longitudinal Data to Sharpen the Analysis of Cascading Disaster Risk. Proceedings of the ICPEE 2019, Zhuhai. https://doi.org/10.1088/1755-1315/432/1/012015 [7] Ribeiro, A.F.S., Russo, A., Gouveia, C.M., et al. (2019) Probabilistic Modelling of the Dependence between Rainfed Crops and Drought Hazard. Natural Hazards and Earth System Sciences, 19, 2795-809. https://doi.org/10.5194/nhess-19-2795-2019 [8] Huang, C. (2011) Discussion on Basic Methods of Risk Analysis. Journal of Natural Disasters, 20, 61-10. [9] Qiong, L. (2011) Flood Risk Assessment Based on the Information Diffusion Method. Communications in Computer and Information Science, 214, 111-117. https://doi.org/10.1007/978-3-642-23321-0_17 [10] Xu, W., Zhuo, L., Zheng, J., et al. (2016) Assessment of the Casualty Risk of Multiple Meteorological Hazards in China. Int J Environ Res Public Health, 13, 222. https://doi.org/10.3390/ijerph13020222 [11] Huang, C. and Huang, Y. (2018) An Information Diffusion Technique to Assess Integrated Hazard Risks. Environ Res, 161, 104-113. https://doi.org/10.1016/j.envres.2017.10.037 [12] Guo, E., Zhang, J., Ren, X.H., Zhang, Q. and Sun, Z.Y. (2014) Integrated Risk Assessment of Flood Disaster Based on Improved Set Pair Analysis and the Variable Fuzzy Set Theory in Central Liaoning Province, China. Natural Hazards, 74, 947-965. https://doi.org/10.1007/s11069-014-1238-9 [13] HAzarika, N., Barman, D., Das, A.K., et al. (2018) Assessing and Mapping Flood Hazard, Vulnerability and Risk in the Upper Brahmaputra River Valley Using Stakeholders' Knowledge AND Multicriteria Evaluation (MCE). J Flood Risk Manag, 11, S700-S716. https://doi.org/10.1111/jfr3.12237 [14] Shouyu, C. and Yu, G. (2006) Variable Fuzzy Sets and Its Application in Comprehensive Risk Evaluation for Flood-Control Engineering System. Fuzzy Optimization and Decision Making, 5, 153-162. https://doi.org/10.1007/s10700-006-7333-y [15] Carreno, M.L., Barbat, A.H., Cardona, O.D. (2014) Numerical Method for the Holistic Evaluation of the Seismic Risk Based on the Fuzzy Sets Theory. Rev Int Metod Numer Calc Dise, 30, 25-34. [16] Li, Q., Zhou, J., Liu, D., et al. (2012) Research on Flood Risk Analysis and Evaluation method Based on Variable Fuzzy Sets and Information Diffusion. Safety Science, 50, 1275-1283. https://doi.org/10.1016/j.ssci.2012.01.007 [17] Wang, W.-C., Xu, D.-M., Chau, K.-W., et al. (2014) Assessment of River Water Quality Based on Theory of Variable Fuzzy Sets and Fuzzy Binary Comparison Method. Water Resources Management, 28, 4183-4200. https://doi.org/10.1007/s11269-014-0738-4 [18] LI, Q. (2013) Fuzzy Approach to Analysis of Flood Risk Based on Variable Fuzzy Sets and Improved Information Diffusion Methods. Natural Hazards and Earth System Sciences, 13, 239-249. https://doi.org/10.5194/nhess-13-239-2013 [19] Chen, J.A., Deng, M., Xia, L. and Wang, H. (2017) Risk Assessment of Drought, Based on IDM-VFS in the Nanpan River Basin, Yunnan Province, China. Sustainability, 9. https://doi.org/10.3390/su9071124 [20] Deng, M., Chen, J., Liu, G., et al. (2016) Risk Assessment of Drought Based on IEAPP-IDM in Qujing, Yunnan Province, China. Advances in Meteorology, 2016, 1-10. https://doi.org/10.1155/2016/7082467 [21] Wang, Y., Sheng, D., Wang, D., et al. (2014) Variable Fuzzy Set Theory to Assess Water Quality of the Meiliang Bay in Taihu Lake Basin. Water Resources Management, 28, 867-880. https://doi.org/10.1007/s11269-014-0521-6 [22] Fang, Y., Zheng, X., Peng, H., et al. (2019) A New Method of the Relative Membership Degree Calculation in Variable Fuzzy Sets for Water Quality Assessment. Ecological Indicators, 98, 515-522. https://doi.org/10.1016/j.ecolind.2018.11.032 [23] Liu, L., Zhou, J., An, X., et al. (2010) Using Fuzzy Theory and Information Entropy for Water Quality Assessment in Three Gorges Region, China. Expert Systems with Applications, 37, 2517-2521. https://doi.org/10.1016/j.eswa.2009.08.004 [24] Chongfu, H. (1997) Principle of Information Diffusion. Fuzzy Sets and Systems, 91, 69-90. https://doi.org/10.1016/S0165-0114(96)00257-6 [25] Huang, C.F. and Shi, Y. (2002) Towards Efficient Fuzzy Information Processing: Using the Principle of Information Diffusion. Physical-Verlag, Heidelberg, German. https://doi.org/10.1007/978-3-7908-1785-0_5 [26] Huang, C. (2015) A Formal Model of Dynamic Risk Analysis of Natural Disasters. Journal of Catastrophology, 30, 1-9.