Next Article in Journal
Estimating Rainfall Erosivity from Daily Precipitation Using Generalized Additive Models
Previous Article in Journal
Development of an IoT System for the Generation of a Database of Residential Water End-Use Consumption Time Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Assessment of the Couple between the Historical Sample and the Theoretical Probability Distributions for Maximum flow Values Based on a Fuzzy Methodology †

by
Matthaios Saridakis
,
Mike Spiliotis
*,
Panagiotis Angelidis
and
Basil Papadopoulos
Section of Hydraulic Engineering, Department of Civil Engineering, Democritus University of Thrace, Kimmeria Campus, 67100 Xanthi, Greece
*
Author to whom correspondence should be addressed.
Presented at the 4thEWaS International Conference: Valuing the Water, Carbon, Ecological Footprints of Human Activities, Online, 24–27 June 2020.
Published: 13 August 2020

Abstract

:
In this article, an adjustment of the extreme theoretical probability distributions upon the sample is proposed, based on the conventional fuzzy linear regression model of Tanaka [1], where all the data must be included within the produced fuzzy band. This is achieved by using the quintile approach, which relates the observed return period with the theoretical cumulative probability. A new contribution of this work is the use of the fuzzified maximum likelihood, as a measure of goodness of fit. The model is applied for real data from the Strymonas River, regarding the annual maximum flow, and finally, useful conclusions are made.

1. Introduction

Fuzzy logic has proven to be a particularly useful tool in the hands of engineers, and its use in recent decades has been widespread in hydrology and hydraulics. Fuzzy linear regression provides a functional fuzzy relationship between dependent and independent variables, where uncertainty manifests itself in the coefficients of the independent variables.
Next, the fuzzy linear regression of Tanaka [1] is used. In general, the fuzzy regression analysis gives a fuzzy functional relationship between the dependent and independent variables [2,3]. In contrast to the statistical regression, the fuzzy regression model of Tanaka [1] has no error term, while the uncertainty is incorporated into the model with the use of fuzzy numbers [4,5]. The data of the fuzzy regression can be either fuzzy or crisp. Usually, the data are rather crisp numbers (observed data), and thus, the uncertainty arises from the used fuzzy model, that is, the fuzzy coefficients. The inclusion property of the produced fuzzy band, that is, the requirement that all the data must be included within the produced fuzzy band, creates the constraints.
Papadopoulos et al. and Spiliotis et al. [5,6] proposed a fuzzy hybrid frequency factor-based method, with the use of fuzzy regression, in order to improve the couple between the theoretical and the observed probability distributions. These articles deal with either annual cumulative streamflow or precipitation. In this work, the annual maximum flow is studied, and furthermore, the fuzzified likelihood is used as an additional measure of suitability.

2. Basic Notations and the Proposed Methodology

2.1. Fundamentals of Fuzzy Sets and Logic

A fuzzy set A on a universe set X is a mapping A : X [ 0 , 1 ] , assigning to each element x X a degree of membership 0 A ( x ) 1 . The membership function A(x) can be also presented as μ A ( x ) .
If A is a fuzzy set, and let any number a [ 0 , 1 ] by the α-cut, A [ α ] and the strong α-cut, A [ α ] + ,the crisp sets are defined, respectively as:
A [ α ] = { x X : A ( x ) α }
A [ α ] + = { x X : A ( x ) > α } ( strong   α - cut ) ,
The 0-cut can be defined as follows:
A [ 0 ] + = { x X : A ( x ) > 0 } ,
In order to have a closed interval containing the boundaries, Hanss [7] proposed the phrase worst-case interval W, which is the union of the 0-strongcut and the boundaries. It is worth noting that, by using the α-cut concept, we can move from the fuzzy sets to the conventional crisp mathematical methodologies.
A special kind of fuzzy sets is the fuzzy numbers. In this work, fuzzy symmetric triangular numbers are used, which are special kinds of fuzzy numbers. The fuzzy symmetric triangular numbers have the following membership function:
μ A ( x ) = { 1 | x a | w , i f a w x a + w 0 , o t h e r w i s e w > 0
in which a is the center and w the spread of the fuzzy number.
The operation of the usual crisp functions, if the inputs are fuzzy sets, can be extended based on the extension principle. In most cases, it is preferable to use a-cuts in the fuzzy analysis [8]. If g is a continuous function in the extension principle, the use of a-cuts can be made by determining the a-cuts of the function f, as follows [9,10]. Then, based on the min intersection, it holds:
{ g L ( x , A ˜ 1 , A ˜ 2 ) α = min { g ( y 1 , y 2 ) | x i s g i v e n , y 1 A 1 [ α ] , y 2 A 2 [ α ] } , g R ( x , A ˜ 1 , A ˜ 2 ) α = max { g ( y 1 , y 2 ) | x i s g i v e n , y 1 A 1 [ α ] , y 2 A 2 [ α ] } .
From the theorem of global existence for maxima and minima of functions with many variables, it is known that, if the domain of a real function is closed and bounded and the real function is continuous, then the function will have its absolute minimum and maximum values at some points in the domain [8,11]. Based on this theorem, it is evident that the α-cut for any real continuous function with real variables in this domain can be determined, given that the inputs are fuzzy triangular numbers [8].

2.2. Fuzzy Linear Regression

The fuzzy linear regression model proposed by Tanaka [1] has the following form:
Y ˜ j = A ˜ 0 + A ˜ 1 x 1 j + ... A ˜ i x i j + .... + A ˜ n x n j with j = , , m   , , n ,
where n is the number of independent variables xij (here only n =1 which is related with the observed return period), M is the number of data, Y ˜ j is the fuzzy predicted value of the dependent variable considering the jth data (here, the maximum annual flow).
According to Tanaka [1], all the data must be included within the produced fuzzy band (inclusion principle). Based on the extension principle and the concept of α-cuts, the inclusion principle is equivalent to (by using fuzzy symmetrical triangular numbers as coefficients):
i = 0 n a i x i j ( 1 h ) i = 0 n c i | x i j | y j i = 0 n a i x i j + ( 1 h ) i = 0 n c i | x i j | , j = 1 , ... , M ,
where ai is the center and ci the spread of the fuzzy number which represents the fuzzy coefficient.
Finally, the sum of the produced semi-widths for the produced dependent variable for all the data is proposed as an objective function:
J = m i n { M c 0 + j = 1 M i = 1 n c i | x i j | } ,
In other words, the problem of fuzzy linear regression is reduced to a linear programming problem [1,4].

2.3. Observed Probabilities

Let a historical sample. The rank order method involves ordering the data from the largest hydrological value to the smallest hydrological value, assigning a rank of 1 to the largest value and a rank of N to the smallest value. An empirical distribution is used to compute the plotting position probabilities as follows [12]:
P ( Q q ) = m Ν + 1 ,
Therefore, the cumulative probability or non-exceedance probability can be determined as follows:
P ( Q < q ) = 1 m Ν + 1 ,
The concept of the observed probability is critical, since the suitability of the used theoretical probability function is based on the comparison between the observed and the theoretical probability values.

2.4. Generalized Extreme Value (GEV) Distribution

The probability density function of the GEV distribution is of the form [13]:
f ( x ) = 1 α [ 1 k ( x u α ) ] 1 / k 1 e [ 1 k ( x u α ) ] 1 / k ,
The GEV distribution function is as follows [14]:
F ( x ) = exp { [ 1 k ( x u α ) ] 1 / k } ,
The distribution function of x given by Equation (12) can be written in the inverse form [13]:
x = u + α k [ 1 ( log F ) k ] ,
By substituting F = 1 1 / T where T is the return period, the T-year quantile estimate of the annual maximum flow, q ^ T , is obtained as follows:
q ^ T = u ^ + α ^ k ^ [ 1 { log ( 1 1 T ) } k ^ ] ,
Csis the skewness coefficient, which can be estimated as:
C S a s 3 ,
Furthermore, a’ is the asymmetry of a sample, of which the unbiased estimation is:
a = N ( N 1 ) ( N 2 ) j = 1 N ( x j x ¯ ) 3 ,
where N is the magnitude of the historical sample. Approximate relationships between the value of k and the skewness coefficient CS, obtained through regression analysis, are given in pages 225, 226 and 227 (7.1.12, 7.1.13 and 7.1.14) in [13]. This approach was adopted by the authors.
In this article, the term [ 1 { log ( 1 1 T j ) } k ^ ] is considered as a crisp number, while the terms u ^ and a ^ k ^ are considered as fuzzy numbers. Therefore, the following fuzzy regression model is modulated:
q ˜ T j = a ˜ 0 + a ˜ 1 x j w h e r e : x j = [ 1 { log ( 1 1 T j ) } k ^ ]
in which q ˜ T j , a ˜ 0 , a ˜ 1 are the fuzzified dependent variable, the constant term and the coefficient of the considered independent variable. It should be clarified that the independent variables, as well as all the observed values, are crisp numbers. The method does not require any transformation of the crisp data to fuzzy data.

2.5. The Extreme Value Type I EV1(2) Distribution

The probability density function of the EV1(2) distribution is given by the following equation [13]:
f ( x ) = 1 α exp [ ( x β α ) e ( x β α ) ] ,
The variable x takes values in the range < x < . The distribution function of x is given by the following equation:
F ( x ) = exp [ e ( x β α ) ] ,
The EV1(2) distribution is a special case of the GEV distribution, discussed in Section 2.4, in which the shape parameter k is equal to zero. The distribution function of EV1(2) (Equation (14)) can be obtained in the inverse form as follows [13]:
x = β α log ( log F ) ,
The T-year quantile is calculated by substituting F = 1 ( 1 / T ) , where T is the return period, to obtainthe following equation:
q ^ T = β ^ α ^ log [ log ( 1 1 / T ) ] ,
More specifically, the term log [ log ( 1 1 / T ) ] is considered as a crisp number, while the terms β ^ and α ^ are considered as fuzzy numbers. Inthe same way, by using the fuzzy linear regression, the following fuzzy relation is determined, with the aim of fuzzy linear regression model:
q ˜ T j = a ˜ 0 + a ˜ 1 x j w h e r e : x j = log [ log ( 1 1 / T j ) ]

2.6. Measures of Suitability

Since a fuzzy approach is dealt, the usual statistical test cannot be used to check the suitability of the proposed model. The first model of suitability, which can be used, is the produced uncertainty of the produced fuzzy band, J (Equation (8)):
M c 0 + c 1 j = 1 M | x j | ,
where M , c 0 , c 1 represent the number of data, the semi-width of the constant term and the semi- width of the independent variable coefficient, respectively.
In this article, the use of the likelihood is used as a measure of suitability (and not to determine the parameters). Since the parameters are fuzzy, then the fuzzified likelihood is constructed by determining the borders of the α-cut for the likelihood, for several representative values of the α-cut, as follows:
{ g L ( x j , a ˜ 0 , a ˜ 1 ) α = min { j = 1 M ln f ( y 1 , y 2 ) | x j i s g i v e n ( c r i s p ) , y 1 a 0 [ α ] , y 2 a 1 [ α ] } , g R ( x j , a ˜ 0 , a ˜ 1 , ) α = max { j = 1 M ln f ( y 1 , y 2 ) | x j i s g i v e n ( c r i s p ) , y 1 a 0 [ α ] , y 2 a 1 [ α ] } ,
To summarize the proposed methodology, the following steps are adopted: (1) modulating the independent and dependent variables based on the Weibull 1939 empirical distribution and the examined theoretical probability distribution. (2) Applying fuzzy regression; (3) the magnitude of the produced fuzzy band and the fuzzfied likelihoods are used as suitability measures.

3. Case Study

The case under investigation is a northern region (Marino Pole) of the Strymonas River (Figure 1). The Strymonas (Struma) River is one of the largest rivers in theBalkan Peninsula in terms of length, since it crosses the Bulgarian and Greek borders. It rises in the Vitosha Mountain in Bulgaria, and has its outlet in the Aegean Sea. Its drainage area is 17,330 km2, of which 10,797 km2 is in Bulgaria, 6295 km2 is in Greece, and the rest is in North Macedonia. In Greece, it is the main waterway feeding and exiting from Lake Kerkini, a significant center for migratory wildfowl. The river’s length is 415 km (of which 290 km is in Bulgaria), making it the country’s fifth-longest and one of the longest rivers that run solely in the interior of the Balkans. Τhe location of the analysis is the Marino Pole, which is located in Bulgaria.
Firstly, the GEV theoretical probability distribution is examined. Since the fuzzy regression enables the determination of two parameters, the third parameter is estimated at the beginning of the method, based on the moment method. Hence, based on the sample, the skewness coefficient is C s = 0.406632 and finally, k ^ = 0.157612 . Subsequently, the fuzzy linear regression of Tanaka with the worst-case interval W is applied, which leads to the following fuzzy curve (Figure 2):
q ˜ T j = ( 145.7 , 27.8 ) + ( 1128.6 , 43.3 ) [ 1 { log ( 1 1 T j ) } k ^ ]
Based on the solution, the sum of the produced semi-widths (J) is:
J = 971 . 7549 m 3 s
Another interesting point of view is that, as can be seen from Figure 2, all the (crisp)observations are included within the produced fuzzy band a property, which is not satisfied according to the conventional methodologies.
Secondly, the extreme value type I [EV1(2)] theoretical distribution is examined. The fuzzy linear regression was applied (Tanaka et al. 1987) based on the aforementioned quantile approach, which leads to the following fuzzy curve:
q ˜ T j = ( 158.62 , 27.01 ) + ( 323.18 , 20.1 ) { log [ log ( 1 1 / T j ) ] }
while the sum of the produced semi-widths is:
J = 1060.4 m 3 s
Both fuzzy curves are depicted in Figure 2 and Figure 3. Based on the magnitude of the produced fuzzy band (that is, the value of the objective function), the GEV theoretical probability distribution is preferred. An interesting point is that every fuzzy linear regression problem based on the Tanaka methodology has a solution. The critical point is the magnitude of the produced fuzzy band. From Figure 2 and Figure 3, it is evident that the fuzzy band can be characterized as functional. Another important point of view is that all the data must be included within the produced fuzzy band, and hence all the empirical probabilities are within the produced fuzzy band.
As mentionedbefore, the second measure of suitability is the produced fuzzy likelihood g, which in its crisp expression, is a continuous function. Based on the extension principle (Equation (24)), the fuzzified likelihood is built by determining several α-cuts. The fuzzified likelihoods for both examined theoretical probability distributions are presented in Figure 4. The two likelihoods have the shape of fuzzy numbers. Even if the comparison between two fuzzy numbers is an ill constructed problem, by applying several widely used measures, we conclude that the fuzzy likelihood of GEV is greater than the fuzzy likelihood of EV1(2) and hence the GEV is preferred according to this criterion. We highlight that the fuzzified likelihood is determined based on the fuzzy regression, or in other words, the fuzzified likelihood is used as a suitability measure, and not in order to determine the parameters of the probability distribution.

4. Conclusions

In order to achieve the couple between the examined extreme theoretical probability distribution and the sample, a hybrid fuzzy regression based approach is developed. The fuzzy regression model is formulated according to the quantile approach, which relates the observed return period with the theoretical cumulative probability. The problem of fuzzy linear regression for crisp data and fuzzy triangular numbers as coefficients concludes to an equivalent linear programming problem.
Two theoretical probability distributions were applied to study the annual maximum river in Strymonas River: the GEV and the EV1(2). To evaluate the proposed approach, two measures of suitability are proposed. The first one is the magnitude of the produced fuzzy bands, which is the objective function regarding the equivalent optimization problem. The second measure is the fuzzified likelihood, which is built according to the extension principle and the solution of the fuzzy regression. It is concluded that, in both cases, a functional uncertainty appears; therefore, the use of the fuzzified theoretical probability function is successful and can be further utilized. Secondly, by comparing both objective functions and the fuzzified likelihood, the use of the fuzzified GEV must be preferred. Hence, the proposed measure of suitability can be used, in order to select the fuzzified theoretical probability function.

Author Contributions

All authors have read and agree to the published version of the manuscript. All authors contributed equally.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tanaka, H. Fuzzy data analysis by possibilistic linear models. Fuzzy Sets Syst. 1987, 24, 363–375. [Google Scholar] [CrossRef]
  2. Tzimopoulos, C.H.; Papadopoulos, K.; Papadopoulos, B. Fuzzyregressionwithapplications in hydrology. Int. J. Eng. Innov. Technol. (IJEIT) 2016, 5, 69–75. [Google Scholar]
  3. Tsakiris, G.; Tigkas, D.; Spiliotis, M. Assessment of interconnection between two adjacent watersheds using deterministic and fuzzy approaches. Eur. Water 2006, 15, 15–22. [Google Scholar]
  4. Spiliotis, M.; Hrissanthou, V. Fuzzy and crisp regression analysis between sediment transport rates and stream discharge in the case of two basins in northeastern Greece. In Conventional and Fuzzy Regression: Theory and Engineering Applications, 1st ed.; Hrissanthou, V., Spiliotis, M., Eds.; Nova Science Publishers: New York, NY, USA, 2018; pp. 1–49. [Google Scholar]
  5. Papadopoulos, C.H.; Spiliotis, M.; Angelidis, P.; Papadopoulos, B. A hybrid fuzzy frequency factor based methodology for analyzing the hydrological drought. J. Desal. Water Treat. 2019, 167, 385–397. [Google Scholar] [CrossRef]
  6. Spiliotis, M.; Angelidis, P.; Papadopoulos, B. A hybrid probabilistic bi-sector fuzzy regression based methodology for normal distributed hydrological variable. Evolv. Syst. 2019, 1–14. [Google Scholar] [CrossRef]
  7. Hanss, M. Applied Fuzzy Arithmetic, an Introduction with Engineering Applications; Springer: Berlin, Germany, 2005; p. 256. [Google Scholar]
  8. Tsakiris, G.; Spiliotis, M. Uncertainty in the analysis of urban water supply and distribution systems. J. Hydroinformatics 2017, 19, 823–837. [Google Scholar] [CrossRef]
  9. Buckley, J.; Eslami, E. An Introduction to Fuzzy Logic and Fuzzy Sets; Springer: Berlin, Germany, 2002; p. 285. [Google Scholar]
  10. Buckley, J.; Eslami, E.; Feuring, T. Solvingfuzzyequations. Fuzzy Math. Econ. Eng. 2002, 91, 19–46. [Google Scholar]
  11. Marsden, J.; Tromba, A. Vector Calculus, 5th ed.; W.H. Freeman and Company: New York, NY, USA, 2003. [Google Scholar]
  12. Weibull, W. A Statistical Theory of the Strength of Materials; Generalstabens Litografiska Anstalts Förlag: Stockholm, Sweden, 1939. [Google Scholar]
  13. Ramachandra Rao, A.; Khaled Hamed, H. Flood Frequency Analysis; CRC-Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  14. Jenkinson, A.F. The frequency Distribution of the Annual Maximum (or Minimum) Values of Meteorological Elements. Q. J. R. Meteorol. Soc. 1955, 87, 158–171. [Google Scholar] [CrossRef]
Figure 1. The Basin of Strymonas River in Bulgaria and the location of the case study.
Figure 1. The Basin of Strymonas River in Bulgaria and the location of the case study.
Environsciproc 02 00022 g001
Figure 2. Graphical representation of the fuzzy linear regression for the Generalized Extreme Value (GEV) distribution.
Figure 2. Graphical representation of the fuzzy linear regression for the Generalized Extreme Value (GEV) distribution.
Environsciproc 02 00022 g002
Figure 3. Graphical representation of the fuzzy linear regression for the extreme value type I [EV1(2)] distribution.
Figure 3. Graphical representation of the fuzzy linear regression for the extreme value type I [EV1(2)] distribution.
Environsciproc 02 00022 g003
Figure 4. The fuzzified likelihoods for GEV and EV1(2) theoretical probability distributions.
Figure 4. The fuzzified likelihoods for GEV and EV1(2) theoretical probability distributions.
Environsciproc 02 00022 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saridakis, M.; Spiliotis, M.; Angelidis, P.; Papadopoulos, B. Assessment of the Couple between the Historical Sample and the Theoretical Probability Distributions for Maximum flow Values Based on a Fuzzy Methodology. Environ. Sci. Proc. 2020, 2, 22. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002022

AMA Style

Saridakis M, Spiliotis M, Angelidis P, Papadopoulos B. Assessment of the Couple between the Historical Sample and the Theoretical Probability Distributions for Maximum flow Values Based on a Fuzzy Methodology. Environmental Sciences Proceedings. 2020; 2(1):22. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002022

Chicago/Turabian Style

Saridakis, Matthaios, Mike Spiliotis, Panagiotis Angelidis, and Basil Papadopoulos. 2020. "Assessment of the Couple between the Historical Sample and the Theoretical Probability Distributions for Maximum flow Values Based on a Fuzzy Methodology" Environmental Sciences Proceedings 2, no. 1: 22. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002022

Article Metrics

Back to TopTop