Application of set-pair analysis and extension coupling model in health evaluation of the huangchuan river, China

Firstly, based on the river health evaluation theory and starting from the integrity of the river system, the river health evaluation index system is selected. Secondly, the index systems of river health assessment including river morphology, hydrological and environmental characteristics, habitat elements and human activity characteristics are determined. Thirdly, set-pair (SPA) evaluation theory and extension were dealt with, a set-pair analysis-extension coupling model is constructed. Finally, taking Huangchuan Rive as an example, the ecological health evaluation on the river channel is carried out by using set-pair analysis-extension coupling model. The results of the comprehensive evaluation using the set-pair analysis-extension coupling method is obtained through calculation. It is proved by practice that the method of river health evaluation using the coupling model of set-pair analysis and extension theory could make full use of the data transformation and extension calculation, which shows the superiority in calculating correlation degree, and could reflect the determination and uncertainty of practical problems. It is possible to comprehensively consider all relevant factors. The paper provides a new idea for solving the relevant problems and obtains more accurate conclusion.


Introduction
River ecological health is an important aspect of river system management. With the development of industrial and agricultural production and social economy, river health has been paid more and more attention by all levels of society and governments (Meyer 1997: Gibson et al. 2016. With the understanding of the river natural attributes and water shortage challenges, the goal of river management has been evolved from a single water resources function to the multiple objectives of river ecosystem management (Thorpjh et al. 2006;Razmkhah et al. 2010). Accordingly, as an assessment tool and technical means, the contents of river health assessment are gradually changed to the integrity and sustainability of river ecosystem in long-term evolution under the dual action of natural and human activities (Zhang et al. 2010;Wei et al. 2019). A river health status should be a decision basis for regional development planning and river managers (Lee et al. 2010;OMID et al. 2014). Because the river health and life involve many factors including natural and social attributes, it has the characteristics of diversity, variability and uncertainty, therefore, the river ecological health evaluation system belongs to the typical complex system (Wei et al. 2019;Yang et al. 2016). Although the previous methods, such as variable fuzzy evaluation, gray theory, extension and set-pair analysis, have made beneficial progress, there still exist some defectives (Zhai et al. 2016;Ladon et al. 1999). For example, extension theory can take all factors into account comprehensively, but the correlation degree is often calculated with the interval important points, and consider those points as the best scheme, which could miss some important constraints and result in differences between the results and actual situations (Gao et al. 2017;Gao et al. 2007). The set pair analysis method could solve the intermediate transformation problems, but it is difficult to determine the coefficient of difference in the evaluation (Yang et al. 2015;Zhang et al. 2014). The combination of set-pair analysis with extension theory could not only apply extension theory to calculate correlation degree, 1 3 198 Page 2 of 9 but also reflect the determination and uncertainty of practical problems, therefore, the accuracy of river health evaluation should be improved. By using this method, the grade of samples to be tested could be reasonably assessed. Therefore, based on the construction of comprehensive evaluation index system of river health, this paper systematically studies SPA evaluation theory, extension theory and set-pair analysis-extension coupling method. The improved set-pair analysis method and extension theory are introduced. The set pair analysis and identical inverse partition method are applied to extensional set domain (Ellmore et al. 2015;Erst et al. 2012). A set-pair analysis-extension coupled comprehensive evaluation model for river health evaluation is constructed. A new perspective to evaluate river health is explored. Finally, taking the Huangchuan River (which is a main tributary of the Huaihe River) as an example, the ecological health evaluation is carried out by using the studied methods. The results of comprehensive evaluation using set-pair analysis-extension coupling method were obtained. It is proved from the results that the evaluation results have an important reference value for the water resources rational exploitation, the improvement of water efficiency and ecological environment.

Construction of river health comprehensive evaluation standard index system.
The comprehensive evaluation index system for river health is constructed according to the following five principles (Chen et al. 2014;Whited et al. 2012); (1) to form a comprehensive and relatively independent indicator system;(2) easy operation, easy access to evaluation data and information; (3) qualitatively ref lect the degree of river health; (4) quantitative calculation; (5) real description of the river status. Based on those consideration, through literature investigation, field investigation, expert interview, combined with the character istics of ecological environment and social-economic characteristics in river basin, the standard index system of river health evaluation was finally determined, which includes the characteristics of river morphology, hydrology, environment, habitat elements and human activities, see Table1.

Improved set-pair analysis based on extension theory
(1) SPA evaluation theory Set pair analysis was proposed by Zhao Keqin in 1989 (Zhao 2000), which is based on the idea that divides the system's certainty into two aspects:" concordant "and" opposite". Uncertainty could be defined by the "difference". The concordant, opposite and difference (three aspects) are interconnected, influenced and transformed mutually (Wohle 2017;Xia and Yan 2004).Here,we introduce the correlation degree to express the uncertainty of the system, to quantify and identify the uncertainty.
Let set pair H = (A,B) is composed of set A and B. The set pair H could be separated to obtain N characteristics, in which, A and B share the opposite number of characteristics S、P, the remaining characteristics are F, which can be expressed using the following formula (1) In which, μ AB is correlation degree; i is difference uncertainty coefficient; i ∈ (-1, 1); j is opposite degree coefficient; a、b、c is the concordant degree coefficient, difference degree and opposite degree coefficient under a certain characteristic respectively, which satisfies the condition a + b + c = 1.
Such kind of method can only make a rough analysis on the membership grade according to the a,b,c. It neglects the degree of comprehensive membership, and also ignores the importance of indicators and the coefficient of difference.
(2) Theory of extenics Born in 1983, Extenics is a new subject overlapped by mathematics, engineering sciences and philosophy (Liu. 2015). The extension theory takes physical element as the logical cell, and takes physical element theory and set-extension theory as the basis(Mingwu, Juliang, 2014). The physical element takes the ordered three dimension R = (thing, characteristic, quantity) = (N,C,V) as the basic elements to describe an event. Different objects can be represented by the elements with the same feature. Physical elements are extensible. Correlation degree k(x 0 ) is the degree which an element belongs to an evaluation level, which can be calculated with the following Eq. (2).
In which, ρ (x 0 , X) and ρ (x 0 , X 0 ) are physical moment of point x 0 to domain X and X 0, respectively ( Fig. 1), (x 0 , X)-(x 0 , X 0 ) is the location of the point x 0 to the relevant two domains.

Improved set pair analysis-extension coupling method
(1) The basic principle and evaluation steps The flow chart of improved set pair analysis and extension coupling method is shown in Fig. 2 (2) Model construction For an evaluation problem, let N = N 1 , N 2 , …, N k , …, N K (k = 1,2, …, K) be entire domain, C j is the jth (j = 1,2, …, m)evaluated indicator, v jk represents quantitative value of element N k relevant to indicator C j , then, the evaluated value for indicator k could be defined as the concordant R k by using physical element. i,e.
in which, V k = (a jk , b j ) is the value domain of Nk relevant to indicator Ci, then, the entire domain R k is composed of the concordant elements from evaluated indicators, i.e. (3) The section domain of the relevant evaluation indicators can be defined as; in which, N represents the total grades of evaluation subjects, (a md , b md ) represents the value domain of indicator C m under certain conditions. Then, R p (p = 1,2, …, P) (the concordant physical elements relevant to total evaluation samples) is as follows; It can be resulted from the relations between extension set and set pair (concordant, differences, opposite) that if sample j falls in the scope of standard k (k = 1,2, …, K), then, the relation between sample and k is concordant, i.e., it locates in the positive domain of the standard extension set, X 0 = (F j,k , F j,k-1 ), the calculation model of relevant correlation membership degree is as follows; in which, F j,k and F j,k+1 refer to limitation values for different grades, respectively, μ k (x ij ) is correlation membership degree of standard positive domain constructed from indicator j in the grades k relevant to evaluation sample i.
If the ith sample j falls in the scope of adjacent grades of k, i.e. k-1 (k > 2) or k + 1, and still in the corresponding transition positive domain, X 1 = (F j,k-1 , F j,k )or X 2 = (F j,k+1 , F j,k+2 ), and the relation of j and k could be defined as difference, the relevant calculation equation is as follows; in which, F j,k-1 and F j,k+2 are limitation value of grades respectively; ρ(x ij , X) and ρ(x ij , X 0 )are the distance of the measured value of sample i to the indicator j of the grade k in the extension positive and standard domain respectively. In formula (10), when k = 1, then F j,k-1 = F j,k , when k = K, then, F j,k+2 = F j,k. When the sample j is out of the extension positive domain,i.e. X = (F j,k-1 , F j,k+2 ), and falls in the adjacent grade k -2 or k + 2, or out of any standard grades, then, in the extension negative set domain Y, then, the relation of sample j and k is opposite, i.e. μ k (x ij ) = -1. As mentioned above, the relation of set pair and extension domain constructed by measured data and standard grade is shown in Fig. 3 Comprehensive indicator's weight and correspondent membership degree μ k (p) are calculated respectively by using the following equations.
in which, w i is indicator weight. The corresponding assessment criteria are defined as, in which, K refers to determined evaluation grade.

River background and data sources
The Huangchuan river is a tributary on the south bank of the Huaihe River, which is originated from the northern foot of Dabie Mountain and flows through Xinxian, Guangshan and Huangchuan county, empties into the Huaihe River at the west of Sihetai village. The total length of the main channel is . The upper reach of the Huangchuan river flows through the steep Dabie Mountains, the middle reach pass through hilly area with gentle slope. The region below Huangchuan is plain area and belongs to agricultural and industrial region. In the Huangchuan river basin, the mountain area amounts to 75%. The overall terrain is high in south and low in north. The trend of terrain is inclined from southwest to northeast. The land surface elevation varies from 35 to 1010 m. The average annual runoff of the river is 1200 million m 3 . The annual average suspended sediment transport is 216,000 t. Annual sediment transport modulus is 132 t/(km 2 ). A schematic diagram of the river system in the Huangchuan basin is shown in Fig. 4. By the end of 2016, several water projects have been built in the region, there were 1 large reservoir,4 medium-sized reservoirs,10 small reservoirs and many sluice dams in the river basin.
In this paper,14 control sections were used to evaluate the health grades of Huangchuan River. From August to October in 2017, the riparian zone, water quality and biological conditions at 14 control section points were investigated and monitored by the Project Team. Systematical data including water quality, ecosystem status, water resources and flow environment were obtained. We also conducted a questionnaire survey for the public satisfaction on river systematic management. On the basis of investigation and monitoring, the information provided by Xinxian, Guangshan, Huangchuan and Xinyang Water Conservancy Bureau was collected, which could be meet the requirement for basic data. Based on reviewing results using standard methods, the data reliability meets evaluation requirements.

Determination of physical elements model in the entire and section domain
The physical element matrix (R 1 、R 2 、R 3 、R 4 ) and section matrix of river health evaluation both in entire and section domain could be constructed respectively, as follows; in which, R 1 、R 2 、R 3 、R 4 and R m are the matrix in entire and section domain respectively, N is the grade of evaluation objectives; C is an evaluation characteristic factor.

Ecological health evaluation on Huangchuan River.
In the physical element extension matrix of health evaluation of Huangchuan river, there are 16 evaluation factors. The evaluation concordant physical element is R 2017. The methods of set pair and extension coupling are based on Eq.
(3)-(10). The evaluation results for different indicators are seen in Table 2. Comparing with Table 1, it can be seen from Table 2 that water quality category c 7, Sewage treatment rate c 8 are illness, which means the Huangchuan river water pollution is serious, at same time, it shows low sewage treatment ability. The riparian integrity index c 3 is low, which shows that there is no natural buffer zone, lack of ecological barriers between rivers and land, the ecological function is weak. Low selfregulation ability leads to weak sewage interception and low sewage load capacity as well as low self purification. The evaluation result of satisfaction index of ecological flow c 4 shows that the water resources in the Huangchuan river basin could meet the ecological water demand in the relevant river reaches without water transfer from outside basin, but, due to serious water pollution, water quality could not meet the requirement, thus, the river health is affected seriously. River channel corridor connectivity index c 11 shows that the river channel corridor system connectivity in Huangchuan basin is low. The main reason is that during recent years many reservoirs, sluice gates and water projects have been built. Ecological improvement projects are carried out in the downstream urban areas and many rubber dams have been constructed. Although such projects have played an active role in ensuring river channel ecological water use, urban water use, landscape construction and environmental improvement, it inevitably cuts off the longitudinal connectivity of the river channel, resulting in low connectivity of channel ecological corridors, which in turn leads to the reduction of the effective transfer of flow materials. The river energy and information transition would be affected. The diversity of species would be reduced. Eventually, the level of river health grade of the biodiversity index c 10 was reduced.
Indicators c 2 、c 14 fall in sub-health domain. The reason may be related to the unreasonable development of water use and sandstone exploitation, which increases the soil erosion in the riparian zone, results in river channel silt deposition and river channel block up. The flood control standards would be lowered, and the ecological security of the river is threatened.
Inputting the correlation degree in Table 2 and the index weights into formula (11), the integrated correlation degree of the health state of the Huangchuan river can be obtained, see Table 3. The results show that the health degree of the Huangchuan river is sub-health, the main reasons are water resources over-development, low river channel corridor longitudinal connectivity, lack of riparian natural buffer zone and water pollution. The evaluation results are consistent with the present actual situation of the river management. It also shows that the method of set pair and extension coupling model is suitable to carry out river ecological health assessment.

Discussion and conclusions
The paper firstly put forward the standard index system of river health evaluation, which is qualitatively reflection of river health degree, and real description of the river status. The standard index system of river health comprehensively makes consideration of the relevant characteristics such as river morphology, hydrology, environment, habitat elements and human activities.
Then, set pair analysis-extension coupling method was studied. The relationship between set pair and extension could be classified as concordant, different, opposite in different domains. The coupling model of set pair and extension was constructed. Comprehensive indicator's weight and correspondent membership degree are calculated respectively by using the relevant models. Taking Huangchuan river as an example, the physical elements model in the entire and section domain were determined. Ecological health evaluation on Huangchuan river was carried out. Through calculation, the evaluation results for different indicators were achieved.
The method of river health assessment based on the coupling evaluation theory of set-pair analysis and extension makes full use of the extension transformation and calculation, which shows the superiority in calculating correlation degree, and can reflect the determination and uncertainty of practical problems. The model comprehensively considers different factors. The present study provides a new idea in solving the problem of determining the different degree of set pair analysis, and obtains more accurate conclusions.
It is shown by the evaluation results that the health status ot the Huangchuan river is unhealthy. It is necessary to take relevant technical measures such as river regulation, pollution control and ecological restoration as soon as possible to improve the river ecological health so that it can be gradually restored to a healthy environment.
In view of illegal encroachment on river channels, low quality of water environment, insufficient environmental flow and serious damage to water ecology, further relevant measures are needed, such as strengthening management of livestock and poultry breeding, improving control of agricultural non-point source and municipal solid waste treatment in Guangshan and Huangchuan counties along the middle and lower reaches of the river, and rising treatment rate of sewage discharge. Therefore, it is necessary to carry out natural resources control and utilization with the aim of improving Huangchuan river ecological health. Through the rational development of water resources, rising water efficiency, realizing water pollution and other measures, the health level of the Huangchuan river would be raised with a view to sustainable water resources development.