# Sediment Deposition Risk Analysis and PLSR Model Research for Cascade Reservoirs Upstream of the Yellow River

Authors: Jie Yang; Jing Ma; De-xiu Hu; Lu Wang; Ji-na Yin; Jie Ren
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

It is difficult to effectively identify and eliminate the multiple correlation influence among the independent factors by least-squares regression. Focusing on this insufficiency, the sediment deposition risk of cascade reservoirs and fitting model of sediment flux into the reservoir are studied. The partial least-squares regression (PLSR) method is adopted for modeling analysis; the model fitting is organically combined with the non-model-style data content analysis, so as to realize the regression model, data structure simplification, and multiple correlations analysis among factors; meanwhile the accuracy of the model is ensured through cross validity check. The modeling analysis of sediment flux into the cascade reservoirs of Long-Liu section upstream of the Yellow River indicates that partial least-squares regression can effectively overcome the multiple correlation influence among factors, and the isolated factor variables have better ability to explain the physical cause of measured results.


---

## Body

## 1. Introduction

Sediment deposition and the balance problems [1, 2] of silt erosion and deposition are bound to be faced in the development and construction of the reservoir project, and they are also the key technical problems which have to be properly solved during project planning, design, construction, and operation management. The Yellow River is a heavy sediment-laden river. A rational prediction of sediment production and ensuring enough sediment transport water are the basic conditions of river flow with sand into sea; otherwise the sedimentation may occur, raise the riverbed, and aggravate flood threat. For this reason, researchers have attached great importance to the sedimentation analysis of the Yellow River and the fitting, forecasting, and other scientific research problems for sediment deposition, sediment production, and water demand for sediment transport. Compared with a single reservoir, basin cascade development mode, the river flow is under the influences of all levels of reservoirs storage and water discharges so that the risk of reservoir sedimentation becomes more serious and the relation between silt erosion and deposition is more complex. After the completion of a reservoir, starting with its initial impoundment, sediment is accumulation in the upriver reservoir area, especially reservoir area by dam. Bed load sediment blocked by dam is gradually accumulating in the reservoir area; meanwhile suspended load sediment is gradually deposited because of the flow velocity’s decrease [3]. Improperly desilting design may also cause severe sedimentation after impoundment [4]. Downstream of this reservoir, driven by economic returns, the discharged volume may be much lower than the demand for sediment transportation which causes the risk of sand sedimentation.Therefore, research of the risk of sediment deposition, the amount of sediment, and the balance of silt erosion and deposition in river basin cascade development mode can optimize sediment configuration, protect the built reservoir’s usable capacity, extend the life service of the reservoir, and operate cascade reservoirs correctly. This has far-reaching significance for the overall efficiency of the cascade development.In the research on prediction model of sediment erosion’s changes, the domestic and foreign researchers have achieved some achievement, and the prediction models can be roughly divided into three categories: first: conceptual model of sediment deposition, such as regression analysis, mathematical statistics method [5]; second: numerical sediment model, such as hydrodynamics, hydrology, and sediment mathematical model [6]; and third: black box sediment model, such as fuzzy mathematical method, the gray system method, and the ANN model [7]. Although each prediction model has a certain level of applicability, with the influence of mathematical method, optimization solution, and practical application, there are some inevitable limitations of each of them. In the modeling analysis of reservoir sediment storage, various multiple regression analyses adopted intensively nowadays were established based on the assumption that there is no close correlation between each model’s factor. While the actual situation is among impact factors of reservoir sediment meteorology, engineering, and social and economic factors, usually multiple correlation and varying degrees of uncertainty are existing. The most conventional least-squares regression analysis cannot effectively overcome the impact of the multiple correlation and relevant uncertainties among various model factors; even if multiple correlation coefficient R has a high accuracy index, it can only represent that fitting accuracy between the model and measured data is high, but not showing the model is really effective in the analysis of the influence on each variable factor. If the model results are directly used to operationally manage the sediment flux into the reservoir, it may take a great risk.Partial least-squares regression (PLSR) is a new multivariate data analysis method presented from application areas and mainly adopted for regression modeling analysis of single or multiple dependent variables on the multiple independent variables; it can effectively solve many problems which ordinary regression cannot solve. PLSR can also organically combine the basic functions of regression modeling, principal component analysis, and canonical correlation analysis. Within an algorithm, it can achieve regression modeling data structure simplification and correlation analysis between the two groups of variables. PLSR has been widely applied in many areas. Chun and Keleş used PLSR for high- dimensional genomic data analysis and they proposed an efficient implementation of sparse partial least-squares regression and compared it with well-known variable selection and dimension reduction approaches via simulation experiments and proved the practical utility of sparse partial least-squares regression in a joint analysis of gene expression and genome-wide binding data [8]. P. P. Roy and K. Roy explored the fact that the optimum variable selection strategy for PLSR used a model dataset of cytoprotection data [9]. Carrascal et al. used PLSR as an alternative to current regression method in ecological studies, which pointed out that PLSR was more reliable than other techniques in identifying relevant variables and their magnitudes of influence [10].Based on analysis of the complex causal relationships between the amount of sediment and their influencing factors, this paper proposed to adopt PLSR method to eliminate interference caused by multiple correlations and uncertainties of influencing factors of sediment flux into the reservoir in order to achieve model fitting and prediction of sediment amount of the reservoir in the cascade development mode, providing reliable technical support for joint dispatching, optimization sediment configuration, and extending the life of reservoir.The main contents of this paper are as follows. The second section is the risk analysis of sediment deposition influence under the mode of drainage basin cascade development. The third section introduces the modeling ideas and procedures of sediment flux into the reservoir PLSR model. The fourth section selects the most representative cascade reservoir upstream of the Yellow River as the modeling instantiated application and modeling analysis. Finally, the results obtained are discussed.


## 2. Sediment Deposition Risk Analysis of River Basin Cascade Development

### 2.1. Analysis of Influencing Factors

The important means to prevent siltation of the Yellow River is strengthening its scientific use and management to ensure the ecological base flow and sediment transport. The main factors influencing the Yellow River sediment deposition are as follows [11, 12].(
1
) There are natural factors, including weather condition, vegetation condition, terrain condition, geological condition, and the inflow water and sediment condition.(
2
) There are engineering factors, including the impact of dam impoundment, the project construction, and the water storage of cascade reservoir.(
3
) Social and economic factors: during social economic development, irrational use of the river water resources may lead to a lack of base flow, discontinuous flow, and the river flow less than sediment transport demand, which are important causes of sediment deposition.


### 2.2. Risk Analysis of Sediment Deposition

The comprehensive analysis showed the risk of sedimentation basin cascade development mainly in the following aspects [13, 14].(
1
) There are environmental risks, reduction of arable land, soil desertification, fertility deterioration, and air quality deterioration.(
2
) Project benefit and risk: the constructions of cascade reservoirs block the original natural river, which lead to the river’s flow pattern changes and the multilevel noncontinuous river morphology. As the effect of reservoir sediment affects the rivers erosion and sediment, the balance of the original river sediment is destroyed. Sediment deposition will decrease flood control capacity and limited storage capacity and shorten the operation life of reservoir [15].(
3
) Aggravating flood risks: during the flood, sedimentation will raise the flood level, which significantly increases the possibility of overrun of water height and river capacity under the same flow rate. Meanwhile, flood losses will be aggregated by the sedimentation in the flood storage region.(
4
) Other risks: clogged or silted channel, shipping difficulty, aggravation of water turbine wear, impact on the safe operation of the hub, siltation range extended, increased upstream flooding losses, downstream river’s erosion and deformation influenced by reservoir discharged water, and pollutants attached to sediment affect water quality in the reservoir and so on.


## 2.1. Analysis of Influencing Factors

The important means to prevent siltation of the Yellow River is strengthening its scientific use and management to ensure the ecological base flow and sediment transport. The main factors influencing the Yellow River sediment deposition are as follows [11, 12].(
1
) There are natural factors, including weather condition, vegetation condition, terrain condition, geological condition, and the inflow water and sediment condition.(
2
) There are engineering factors, including the impact of dam impoundment, the project construction, and the water storage of cascade reservoir.(
3
) Social and economic factors: during social economic development, irrational use of the river water resources may lead to a lack of base flow, discontinuous flow, and the river flow less than sediment transport demand, which are important causes of sediment deposition.


## 2.2. Risk Analysis of Sediment Deposition

The comprehensive analysis showed the risk of sedimentation basin cascade development mainly in the following aspects [13, 14].(
1
) There are environmental risks, reduction of arable land, soil desertification, fertility deterioration, and air quality deterioration.(
2
) Project benefit and risk: the constructions of cascade reservoirs block the original natural river, which lead to the river’s flow pattern changes and the multilevel noncontinuous river morphology. As the effect of reservoir sediment affects the rivers erosion and sediment, the balance of the original river sediment is destroyed. Sediment deposition will decrease flood control capacity and limited storage capacity and shorten the operation life of reservoir [15].(
3
) Aggravating flood risks: during the flood, sedimentation will raise the flood level, which significantly increases the possibility of overrun of water height and river capacity under the same flow rate. Meanwhile, flood losses will be aggregated by the sedimentation in the flood storage region.(
4
) Other risks: clogged or silted channel, shipping difficulty, aggravation of water turbine wear, impact on the safe operation of the hub, siltation range extended, increased upstream flooding losses, downstream river’s erosion and deformation influenced by reservoir discharged water, and pollutants attached to sediment affect water quality in the reservoir and so on.


## 3. Partial Least-Squared Regression Model of the Sediment Storage Volume in Cascade Development Reservoirs

### 3.1. The Basic Idea of Partial Least-Squares Regression

Suppose that there arep independent variables {x1,…,xp} and q dependent variables {y1,…,yq}. In order to study the statistical relationship between dependent variables and independent variables, n sample points are collected and constitute data table of the independent variables and dependent variables: X={x1,…,xp}n×p and dependent variables Y=y1,…,yqn×q. Extract the principal components t1 and u1 from X and Y, respectively, where t1 and u1 are the linear combination of {x1,…,xp} and {y1,…,yq}. t1 and u1 should carry their variation information of their own data table as much as possible to maximize their relevance while extracting the two components [16, 17].After extracting the first principal componentst1 and u1, make the regression of X on t1 and Y on u1, respectively; then use residual information to extract the main ingredients for the second round. Continue to implement the return of X on t1 and Y on u1 until a satisfactory accuracy is achieved [18]. Ultimately extract the m components t1,t2,…,tm from X; the regression of yk on t1,t2,…,tm can then be transformed into the regression equation of yk (k=1,2,…,q) on the original variables x1,…,xp; this completes the partial least-squared regression modeling.


### 3.2. PLSR Model of Sediment Storage Volume in Cascade Development Reservoirs

Denoting dependent variables (sediment storage volume) byy∈Rn, then the independent variable set is X=[x1,x2,…,xm], xj∈Rn. The standardized independent variables matrix of X is denoted as E0, and the standardized dependent variables matrix of y is denoted as F0:(1)F0i=yi-y¯sy,i=1,2,…,n,where y¯ is the mean value of y and sy is the standard deviation error of y.In the first step, the unit vectorsE0 and F0 are known; for the single dependent variabley, its component u1 equals F0; then(2)w1=E0F0E0F0=1∑j=1mr2xj,yrx1,y⋮rxm,y,t1=E0w1=1∑j=1mr2xj,yrx1,yE01+⋯+rxm,yE0m,where r(xj,y) is the simple correlation coefficient.If so, the optimization value of objective function is(3)θ1=E0′F0=∑j=1mr2xj,y.The regression values of E0 and F0 on t1 are calculated, respectively;(4)E0=t1p1+E1,(5)F0=t1r1+F1,where E1 and F1 are residual matrices, p1 and r1 are regression coefficients, and r1 is invariant and p1=E0t1/t12.In the second step replaceE0, F0 with E1, F1, respectively. Repeat the first step to implement the regression of E1 and F1 on t2 and obtain the regression coefficients p2 and r2. Conduct partial least-squares regression and so on.In theh step, according to the cross-validation test result of the model’s accuracy, eventually supposing extracted components t1,t2,…,th, then(6)wh=Eh-1F0Eh-1F0,th=Eh-1wh,ph=Eh-1thth2,Eh=Eh-1-thph.Conduct the regression ofF0 on t1,t2,…,th; then will get(7)F^0=r1t1+r2t2+⋯+rhth.Becauset1,t2,…,th is a linear combination of E0, th=Eh-1wh=E0wh∗; therefore(8)F^0=r1E01w1∗+⋯+rhE0hwh∗=E0k∑k=1hrkwk∗.Ordery∗=F0, xj∗=E0j, and aj=∑k=1hrkwkj∗ (j=1,2,…,m) is regression coefficient of xj∗, then:(9)y^∗=a1x1∗+a2x2∗+⋯+amxm∗,where y^∗ is the fitted values of the sediment flux into the reservoir (the dependent variable), xj∗ (j=1,2,…,m) is the various influencing factors of the volume of sand storage (the independent variables), the formula (9) is just the PLSR model expression for the sediment storage of a reservoir in cascade development.


### 3.3. Factors Multicorrelation Analysis

According to principal component analysis, under the premise of the sediment storage data with minimal loss, make dimension reduction for high-dimensional data systems; when the originalp-dimensional data is reduced to two-dimensional system, the location of all the sample points can be drawn out on the two-dimensional plan, which can be directly observed as similar structure of p-dimensional space. After establishing the PLSR model, use the extracted components t1 and u1 directly, do a chart view of t1/u1, make (t1(i),u1(i)) as coordinate points of the sample points, and then plot the location of all sample points to analyze distribution and structural similarity of sample points in high-dimensional space.According to the theory and the specific steps of PLSR modeling, we can use the MATLAB language to compile the PLSR analysis program of sediment storage in cascade development.


## 3.1. The Basic Idea of Partial Least-Squares Regression

Suppose that there arep independent variables {x1,…,xp} and q dependent variables {y1,…,yq}. In order to study the statistical relationship between dependent variables and independent variables, n sample points are collected and constitute data table of the independent variables and dependent variables: X={x1,…,xp}n×p and dependent variables Y=y1,…,yqn×q. Extract the principal components t1 and u1 from X and Y, respectively, where t1 and u1 are the linear combination of {x1,…,xp} and {y1,…,yq}. t1 and u1 should carry their variation information of their own data table as much as possible to maximize their relevance while extracting the two components [16, 17].After extracting the first principal componentst1 and u1, make the regression of X on t1 and Y on u1, respectively; then use residual information to extract the main ingredients for the second round. Continue to implement the return of X on t1 and Y on u1 until a satisfactory accuracy is achieved [18]. Ultimately extract the m components t1,t2,…,tm from X; the regression of yk on t1,t2,…,tm can then be transformed into the regression equation of yk (k=1,2,…,q) on the original variables x1,…,xp; this completes the partial least-squared regression modeling.


## 3.2. PLSR Model of Sediment Storage Volume in Cascade Development Reservoirs

Denoting dependent variables (sediment storage volume) byy∈Rn, then the independent variable set is X=[x1,x2,…,xm], xj∈Rn. The standardized independent variables matrix of X is denoted as E0, and the standardized dependent variables matrix of y is denoted as F0:(1)F0i=yi-y¯sy,i=1,2,…,n,where y¯ is the mean value of y and sy is the standard deviation error of y.In the first step, the unit vectorsE0 and F0 are known; for the single dependent variabley, its component u1 equals F0; then(2)w1=E0F0E0F0=1∑j=1mr2xj,yrx1,y⋮rxm,y,t1=E0w1=1∑j=1mr2xj,yrx1,yE01+⋯+rxm,yE0m,where r(xj,y) is the simple correlation coefficient.If so, the optimization value of objective function is(3)θ1=E0′F0=∑j=1mr2xj,y.The regression values of E0 and F0 on t1 are calculated, respectively;(4)E0=t1p1+E1,(5)F0=t1r1+F1,where E1 and F1 are residual matrices, p1 and r1 are regression coefficients, and r1 is invariant and p1=E0t1/t12.In the second step replaceE0, F0 with E1, F1, respectively. Repeat the first step to implement the regression of E1 and F1 on t2 and obtain the regression coefficients p2 and r2. Conduct partial least-squares regression and so on.In theh step, according to the cross-validation test result of the model’s accuracy, eventually supposing extracted components t1,t2,…,th, then(6)wh=Eh-1F0Eh-1F0,th=Eh-1wh,ph=Eh-1thth2,Eh=Eh-1-thph.Conduct the regression ofF0 on t1,t2,…,th; then will get(7)F^0=r1t1+r2t2+⋯+rhth.Becauset1,t2,…,th is a linear combination of E0, th=Eh-1wh=E0wh∗; therefore(8)F^0=r1E01w1∗+⋯+rhE0hwh∗=E0k∑k=1hrkwk∗.Ordery∗=F0, xj∗=E0j, and aj=∑k=1hrkwkj∗ (j=1,2,…,m) is regression coefficient of xj∗, then:(9)y^∗=a1x1∗+a2x2∗+⋯+amxm∗,where y^∗ is the fitted values of the sediment flux into the reservoir (the dependent variable), xj∗ (j=1,2,…,m) is the various influencing factors of the volume of sand storage (the independent variables), the formula (9) is just the PLSR model expression for the sediment storage of a reservoir in cascade development.


## 3.3. Factors Multicorrelation Analysis

According to principal component analysis, under the premise of the sediment storage data with minimal loss, make dimension reduction for high-dimensional data systems; when the originalp-dimensional data is reduced to two-dimensional system, the location of all the sample points can be drawn out on the two-dimensional plan, which can be directly observed as similar structure of p-dimensional space. After establishing the PLSR model, use the extracted components t1 and u1 directly, do a chart view of t1/u1, make (t1(i),u1(i)) as coordinate points of the sample points, and then plot the location of all sample points to analyze distribution and structural similarity of sample points in high-dimensional space.According to the theory and the specific steps of PLSR modeling, we can use the MATLAB language to compile the PLSR analysis program of sediment storage in cascade development.


## 4. Application and Analysis of Engineering Case

### 4.1. Basic Information

The upper reaches of the Yellow River between Longyangxia and Liujiaxia section mainly have six large-sized cascade hydropower stations, which are Longyangxia, Laxiwa, Lijiaxia, Gongboxia, Jishixia, and Liujiaxia from upstream to downstream, as shown in Figure1.Figure 1
Sketch map of cascade hydropower stations layout between Longyangxia and Liujiaxia section in the Yellow River.Longyangxia hydropower station is the leading hydropower station in the main stream of the Yellow River, with a total reservoir capacity of 274 × 108 m3 and less sediment concentration, at which control basin is one of the few large-scale runoff regional stabilities and is located in 1687 km from the source of the Yellow River. Jishixia hydropower station is fifth large-scale cascade hydropower station in the Long-Liu segment of the Yellow River, with the dam located at the outlet of Jishixia canyon in the Yellow River, up from Gongboxia hydropower station about 55 km, down from Liujiaxia hydropower station about 93 km. It is a daily regulating reservoir with the total capacity of 2.94 × 108 m3 [19]. The basic statistics of other cascade hydropower stations in Long-Liu section of the Yellow River are showed in Table 1.Table 1
The basic statistics of other cascade hydropower stations between Longyangxia and Liujiaxia section.
Hydropower stations
Dam height (m)
Dam types
Basin size of controlling (km2)
Normal storage level/dead water level (m)
Total capacity/regulation capacity (108 m3)
Water storage date (year)
Longyangxia
178
VA
131420
2600/2530
247/193.5
1986
Laxiwa
250
VA
132160
2452/2440
10.56/1.5
2009
Lijiaxia
165
VA
136747
2180/2178
16.5/0.6
1996
Gongboxia
127
CFRD
143619
2005/2002
6.2/0.75
2004
Jishixia
100
CFRD
146749
1856/1852
2.94/0.45
2010
Liujiaxia
147
PG
173000
1735/1694
57/41.5
1968
Note: the dam type code in Table1 refers to the ICOLD code.Each cascade hydropower station in the Long-Liu section is distributed between Tangnaihai and Xunhua Hydrometrical Stations. Tangnaihai Hydrometrical Station is located along the eastern edge of the Qinghai-Tibet Plateau, the boundary between the natural runoff and manually adjusted section, which is an important control section on the upper reaches of the Yellow River. It is also a storage station of the Longyangxia Reservoir, with 121,972 km2 control drainage area, accounting for 92.8% of Longyangxia Gorge Reservoir control area; the annual average sediment concentration is 0.61 kg/m3 [20]. Xunhua Hydrometrical Station is the control station in Qinghai Province, as storage station of Jishixia reservoir, with 145,459 km2 control area. Before the year of 1985 the annual runoff was not affected by human activities, but, after 1986, with the building of the reservoirs such as Longyangxia Reservoir, the measured data has included and reflected the impact of the upstream cascade regulating reservoirs. The measured water and sediment statistics of Tangnaihai and Xunhua Hydrometrical Stations are shown in Tables 2~4.Table 2
The measured water and sediment statistics of Tangnaihai Hydrometrical Station (1956 to 2004).
Years
Average annual runoff(108 m3)
Average annual sedimentload (104 t)
Flood season
Full year
Flood season
Full year
1956–1985
127.8
213.1
950.3
1301.8
1986–1996
111.7
187.5
889.7
1289.1
1997–2004
98.9
167.9
839.5
1106.4Table 3
The average sediment load value monthly measured and year allocation of Tangnaihai Hydrometrical Station (1982).
Month
1
2
3
4
5
6
7
8
9
10
11
12
Full year
Sediment load (104 t)
2.2
2.2
11.2
29
145
675.7
461.6
303.3
405.9
165
24.5
4.4
2230
%
0.1
0.1
0.5
1.3
6.5
30.3
20.7
13.6
18.2
7.4
1.1
0.2
100Table 4
The measured water and sediment statistics of Xunhua Hydrometrical Station (1956 to 2004).
Years
Average annual runoff(108 m3)
Average annual sedimentload (104 t)
Flood season
Full year
Flood season
Full year
1956–1985
137.9
232.8
2840.2
3623.6
1986–1996
75.9
192
1490.2
1909
1997–2004
55.2
159.2
894.6
1169.8


### 4.2. Factor Selection and PLSR Model Establishment

Making an exact fitting about the amount of sediment in the reservoir is significant for reasonable prediction of reservoir sedimentation and scientifically determines water demand for sediment. Due to multi-correlations among each factor (such as rainfall amount, natural runoff, sediment transport rate, the amount of water discharged from upper reservoir, water level of the upper reservoir, water level of the downstream reservoir, water temperature and other factors), which influence sediment flux into each cascade reservoir, so adopting the PLSR model established before to fitting.Main factors affecting the sediment storage are rainfall amount within the drainage area, runoff, upstream water level and the discharging water, the backup water level of downstream reservoir and discharged volume of this level’s reservoir, natural river runoff, water temperature, sediment transport rate, and so on. For this reason, in the mode of cascade development, the PLSR model of sediment quantity is composed of runoff factor componentPJ, sediment discharge factor component PS, rainfall factor component PY, water level factor component PH, water temperature factor component PT, discharge water factor component PX, the constant a0, and so forth. That is,(10)P=PJ+PS+PY+PH+PT+PX+a0.According to (8) to reform the above equation into formula (9), we can get the PLSR model expression of sediment storage in the Long-Liu section of the Yellow River under cascade development mode:(11)y^∗=a1x1∗+a2x2∗+⋯+ajxj∗+a0,where y^∗ is fitted values of sediment storage, xj∗ is various influencing factors affecting sediment storage, aj is the partial least-squares regression coefficients corresponding to each xj∗, a0 is a constant, and j represents the amount of sediment storage factor. It is necessary to consider actual state to select factors. For different time periods, cascade reservoirs in different locations, the selected model factors may be different.


### 4.3. Fitting Results and Comparative Analysis of PLSR Model

Considering the influence on the downstream channel and the reservoir sediment transportation caused by water storage of Longyangxia, this paper takes the Longyangxia and Liujiaxia reservoirs as example of built PLSR fitting model of sediment storage. Selecting the monthly average sediment storage capacity from January 1982 to December 1987 as the fitting time, including important periods from the Longyangxia Reservoir, water level gradually increased to the normal water level during initial impoundment so that it can better analyze the influence that water storage of Longyangxia Reservoir for sediment flux into downstream reservoir. A PLSR fitting model for sediment storage of Longyangxia and Liujiaxia Reservoirs is shown in Figures2 and 3.Figure 2
Fitting chart of PLSR model for sediment storage of Longyangxia.Figure 3
Fitting chart of PLSR model for sediment storage of Liujiaxia.Multicorrelation analysis of factors: take rainfall factor setY and transport rate factor set S as example, extract the first components t1(i) and u1(i), and draw plan of t1-u1, as shown in Figure 4. As can be seen, compositions t1(i) and u1(i) appear in roughly linear correlation changes, which suggests that between rainfall factor and transport rate factor of influence reservoir sediment exists in a complex correlation. In this case, the conventional least-squares regression is difficult to separate pressure, temperature, rainfall, and aging components; thus using PLSR to do modeling analysis is more rational and effective.Figure 4
The first principal component distribution of rainfall and sediment transport rate factort1(i), u1(i).According to Figures2 and 3, the following can be found.(
1
) Focusing on the PLSR model of sediment storage of the cascade reservoirs in the Long-Liu section, it largely overcame the multicorrelation between the model factors and can effectively separate the various influencing factors. Meanwhile, model fitting has high precise results that can better reflect the actual changes of sediment storage and explain the impact of each factor on the amount of storage reservoir sediment.(
2
) The fitting sequence is made of monthly average sediment storage of Longyangxia and Liujiaxia Reservoirs, whose sediment storage in the previous 60 months is large while, due to the influence of water impoundment by Longyangxia, (corresponding to those in 1987), the incoming-reservoir sediment production decreases during later 60 months, which indicates that the leading reservoir’s impoundment has an important influence on the cascade development of river sediment transport.(
3
) The sediment storage is larger in the flood season every year, while it decreases obviously in dry seasons. It is characterized by cyclical changes significantly and is in line with the general rules of sediment yield and transport.(
4
) The main factors affecting the amount of sediment storage are the changes of transport rate, and its model factor PS had a large impact on maximum amount of sediment storage; additionally, the influence of runoff factor PJ, water factor PH, and discharge water factor PX is relatively large, while the impact of rainfall and temperature factors PY and PT is relatively small.


## 4.1. Basic Information

The upper reaches of the Yellow River between Longyangxia and Liujiaxia section mainly have six large-sized cascade hydropower stations, which are Longyangxia, Laxiwa, Lijiaxia, Gongboxia, Jishixia, and Liujiaxia from upstream to downstream, as shown in Figure1.Figure 1
Sketch map of cascade hydropower stations layout between Longyangxia and Liujiaxia section in the Yellow River.Longyangxia hydropower station is the leading hydropower station in the main stream of the Yellow River, with a total reservoir capacity of 274 × 108 m3 and less sediment concentration, at which control basin is one of the few large-scale runoff regional stabilities and is located in 1687 km from the source of the Yellow River. Jishixia hydropower station is fifth large-scale cascade hydropower station in the Long-Liu segment of the Yellow River, with the dam located at the outlet of Jishixia canyon in the Yellow River, up from Gongboxia hydropower station about 55 km, down from Liujiaxia hydropower station about 93 km. It is a daily regulating reservoir with the total capacity of 2.94 × 108 m3 [19]. The basic statistics of other cascade hydropower stations in Long-Liu section of the Yellow River are showed in Table 1.Table 1
The basic statistics of other cascade hydropower stations between Longyangxia and Liujiaxia section.
Hydropower stations
Dam height (m)
Dam types
Basin size of controlling (km2)
Normal storage level/dead water level (m)
Total capacity/regulation capacity (108 m3)
Water storage date (year)
Longyangxia
178
VA
131420
2600/2530
247/193.5
1986
Laxiwa
250
VA
132160
2452/2440
10.56/1.5
2009
Lijiaxia
165
VA
136747
2180/2178
16.5/0.6
1996
Gongboxia
127
CFRD
143619
2005/2002
6.2/0.75
2004
Jishixia
100
CFRD
146749
1856/1852
2.94/0.45
2010
Liujiaxia
147
PG
173000
1735/1694
57/41.5
1968
Note: the dam type code in Table1 refers to the ICOLD code.Each cascade hydropower station in the Long-Liu section is distributed between Tangnaihai and Xunhua Hydrometrical Stations. Tangnaihai Hydrometrical Station is located along the eastern edge of the Qinghai-Tibet Plateau, the boundary between the natural runoff and manually adjusted section, which is an important control section on the upper reaches of the Yellow River. It is also a storage station of the Longyangxia Reservoir, with 121,972 km2 control drainage area, accounting for 92.8% of Longyangxia Gorge Reservoir control area; the annual average sediment concentration is 0.61 kg/m3 [20]. Xunhua Hydrometrical Station is the control station in Qinghai Province, as storage station of Jishixia reservoir, with 145,459 km2 control area. Before the year of 1985 the annual runoff was not affected by human activities, but, after 1986, with the building of the reservoirs such as Longyangxia Reservoir, the measured data has included and reflected the impact of the upstream cascade regulating reservoirs. The measured water and sediment statistics of Tangnaihai and Xunhua Hydrometrical Stations are shown in Tables 2~4.Table 2
The measured water and sediment statistics of Tangnaihai Hydrometrical Station (1956 to 2004).
Years
Average annual runoff(108 m3)
Average annual sedimentload (104 t)
Flood season
Full year
Flood season
Full year
1956–1985
127.8
213.1
950.3
1301.8
1986–1996
111.7
187.5
889.7
1289.1
1997–2004
98.9
167.9
839.5
1106.4Table 3
The average sediment load value monthly measured and year allocation of Tangnaihai Hydrometrical Station (1982).
Month
1
2
3
4
5
6
7
8
9
10
11
12
Full year
Sediment load (104 t)
2.2
2.2
11.2
29
145
675.7
461.6
303.3
405.9
165
24.5
4.4
2230
%
0.1
0.1
0.5
1.3
6.5
30.3
20.7
13.6
18.2
7.4
1.1
0.2
100Table 4
The measured water and sediment statistics of Xunhua Hydrometrical Station (1956 to 2004).
Years
Average annual runoff(108 m3)
Average annual sedimentload (104 t)
Flood season
Full year
Flood season
Full year
1956–1985
137.9
232.8
2840.2
3623.6
1986–1996
75.9
192
1490.2
1909
1997–2004
55.2
159.2
894.6
1169.8


## 4.2. Factor Selection and PLSR Model Establishment

Making an exact fitting about the amount of sediment in the reservoir is significant for reasonable prediction of reservoir sedimentation and scientifically determines water demand for sediment. Due to multi-correlations among each factor (such as rainfall amount, natural runoff, sediment transport rate, the amount of water discharged from upper reservoir, water level of the upper reservoir, water level of the downstream reservoir, water temperature and other factors), which influence sediment flux into each cascade reservoir, so adopting the PLSR model established before to fitting.Main factors affecting the sediment storage are rainfall amount within the drainage area, runoff, upstream water level and the discharging water, the backup water level of downstream reservoir and discharged volume of this level’s reservoir, natural river runoff, water temperature, sediment transport rate, and so on. For this reason, in the mode of cascade development, the PLSR model of sediment quantity is composed of runoff factor componentPJ, sediment discharge factor component PS, rainfall factor component PY, water level factor component PH, water temperature factor component PT, discharge water factor component PX, the constant a0, and so forth. That is,(10)P=PJ+PS+PY+PH+PT+PX+a0.According to (8) to reform the above equation into formula (9), we can get the PLSR model expression of sediment storage in the Long-Liu section of the Yellow River under cascade development mode:(11)y^∗=a1x1∗+a2x2∗+⋯+ajxj∗+a0,where y^∗ is fitted values of sediment storage, xj∗ is various influencing factors affecting sediment storage, aj is the partial least-squares regression coefficients corresponding to each xj∗, a0 is a constant, and j represents the amount of sediment storage factor. It is necessary to consider actual state to select factors. For different time periods, cascade reservoirs in different locations, the selected model factors may be different.


## 4.3. Fitting Results and Comparative Analysis of PLSR Model

Considering the influence on the downstream channel and the reservoir sediment transportation caused by water storage of Longyangxia, this paper takes the Longyangxia and Liujiaxia reservoirs as example of built PLSR fitting model of sediment storage. Selecting the monthly average sediment storage capacity from January 1982 to December 1987 as the fitting time, including important periods from the Longyangxia Reservoir, water level gradually increased to the normal water level during initial impoundment so that it can better analyze the influence that water storage of Longyangxia Reservoir for sediment flux into downstream reservoir. A PLSR fitting model for sediment storage of Longyangxia and Liujiaxia Reservoirs is shown in Figures2 and 3.Figure 2
Fitting chart of PLSR model for sediment storage of Longyangxia.Figure 3
Fitting chart of PLSR model for sediment storage of Liujiaxia.Multicorrelation analysis of factors: take rainfall factor setY and transport rate factor set S as example, extract the first components t1(i) and u1(i), and draw plan of t1-u1, as shown in Figure 4. As can be seen, compositions t1(i) and u1(i) appear in roughly linear correlation changes, which suggests that between rainfall factor and transport rate factor of influence reservoir sediment exists in a complex correlation. In this case, the conventional least-squares regression is difficult to separate pressure, temperature, rainfall, and aging components; thus using PLSR to do modeling analysis is more rational and effective.Figure 4
The first principal component distribution of rainfall and sediment transport rate factort1(i), u1(i).According to Figures2 and 3, the following can be found.(
1
) Focusing on the PLSR model of sediment storage of the cascade reservoirs in the Long-Liu section, it largely overcame the multicorrelation between the model factors and can effectively separate the various influencing factors. Meanwhile, model fitting has high precise results that can better reflect the actual changes of sediment storage and explain the impact of each factor on the amount of storage reservoir sediment.(
2
) The fitting sequence is made of monthly average sediment storage of Longyangxia and Liujiaxia Reservoirs, whose sediment storage in the previous 60 months is large while, due to the influence of water impoundment by Longyangxia, (corresponding to those in 1987), the incoming-reservoir sediment production decreases during later 60 months, which indicates that the leading reservoir’s impoundment has an important influence on the cascade development of river sediment transport.(
3
) The sediment storage is larger in the flood season every year, while it decreases obviously in dry seasons. It is characterized by cyclical changes significantly and is in line with the general rules of sediment yield and transport.(
4
) The main factors affecting the amount of sediment storage are the changes of transport rate, and its model factor PS had a large impact on maximum amount of sediment storage; additionally, the influence of runoff factor PJ, water factor PH, and discharge water factor PX is relatively large, while the impact of rainfall and temperature factors PY and PT is relatively small.


## 5. Conclusions

Under the cascade development mode, river flow is influenced by the reservoir impoundment of each cascade and discharge flow, reservoir sedimentation encounters greater risks; the relations between the erosion and deposition become more complicated. For this reason, on the basis of systematic analysis for sediment deposition risk of cascade development in the upstream of the Yellow River’s, mainly studied the model and fitting method of sediment flux into the reservoir under cascade development mode.Considering the difficulty to effectively identify and eliminate the multiple correlation influence among the independent factors by least-squares regression, the partial least-squares regression (PLSR) method is adopted for modeling analysis; the model fitting is organically combined with the non-model-style data content analysis, so as to realize the regression model, data structure simplification, and multiple correlations analysis among factors; meanwhile the accuracy of the model is ensured through cross validity check.The modeling analysis of sediment flux into the cascade reservoirs of Long-Liu section upstream of the Yellow River indicates that partial least-squares regression can effectively overcome the multiple correlation influence among factors, and the isolated factor variables have better ability to explain the physical cause of measured results.



---

*Source: 696015-2015-10-11.xml*