Rivers are drying up most frequently in West Liaohe River plain and the bare river beds present fine sand belts on land. These sand belts, which yield a dust heavily in windy days, stress the local environment deeply as the riverbeds are eroded by wind. The optimal operation of water resources, thus, is one of the most important methods for preventing the wind erosion of riverbeds. In this paper, optimal operation model for water resources based on riverbed wind erosion control has been established, which contains objective function, constraints, and solution method. The objective function considers factors which include water volume diverted into reservoirs, river length and lower threshold of flow rate, etc. On the basis of ensuring the water requirement of each reservoir, the destruction of the vegetation in the riverbed by the frequent river flow is avoided. The multi core parallel solving method for optimal water resources operation in the West Liaohe River Plain is proposed, which the optimal solution is found by DPSA method under the POA framework and the parallel computing program is designed in Fork/Join mode. Based on the optimal operation results, the basic rules of water resources operation in the West Liaohe River Plain are summarized. Calculation results show that, on the basis of meeting the requirement of water volume of every reservoir, the frequency of reach river flow which from Taihekou to Talagan Water Diversion Project in the Xinkai River is reduced effectively. The speedup and parallel efficiency of parallel algorithm are 1.51 and 0.76 respectively, and the computing time is significantly decreased. The research results show in this paper can provide technical support for the prevention and control of riverbed wind erosion in the West Liaohe River plain.
With the development and utilization of water resources,
partial rivers were dried up in arid or semi-arid area. After dried up, bare
riverbed becomes the new cradle land of sandstorm. Since 1972, the lower
reaches of the Yellow River has began to dry up (Gao, 1998). The reasons for
dry up in the watershed of the Yellow River can be concluded into three
aspects: (1) natural reasons, the Yellow River basin is located in an arid
and semi-arid region, at the same time, the river dry up may be aggravated by
global warming; (2) human activities, about 30 billion m
Hydrological station distribution schematic in West Liaohe River Plain.
The statistical results about runoff of each hydrological station in the West Liaohe River Plain.
Units:
Wind erosion caused soil loss and then impoverishes lands (Lal, 2001). The blown dust emission from eroding lands impairs air quality (Lee et al., 2003) and poses a threat to human health (Wilson and Sprengler, 1996). The research of wind erosion prevention and control is mainly focused on the following two aspects: (1) the nonerodible grains are placed on the erodible surface to reduce the degree of wind erosion (Bagnold, 1941); (2) the use of vegetation cover or crop residues to prevent wind erosion (Englehorn et al., 1952; Zhang et al., 2003; Jia et al., 2015). The wind erosion control of bare riverbed or lakebed is quite different from sand land. Due to the technical difficulties, it is hard to implement these measures. Bao et al. (2006) estimated the rational water area and inflow of the Ebinur Lake for controlling wind erosion in the dried up lake basin (Bao et al., 2006). For intermittent dried up reach, there is no better erosion control measures so far.
In this paper we analyze the hydrological characteristics West Liaohe River plain and establish the optimal operation model for water resources based on river desertification control.
The West Liaohe River Plain is formed by alluvial deposits of the West Liaohe
River and its main tributaries. West Liaohe River Basin area is
Comparison of large cross-section measured before and after flood
season at the Taihekou (Xi) station (
Sketch of water conservancy project in West Liaohe River Plain.
This paper collected nearly 8 years (2006
As shown in Table 1, the annual average runoff of Taihekou (Xi) station is
According to the analysis results of hydrological characteristics, the reach of West Liaohe River Plain is divided into intermittent river reach and long dry reach. This paper focuses on the prevention and control of riverbed wind erosion in intermittent river reach.
The river of West Liaohe River Plain has a large amount of sediment.
According to the Taihekou (Xi) station monitoring data, the average sediment
concentration is 9.07 kg m
In the summer flood season in 2011, the maximum flow rate is 585 m
The parameters of river flow evolution in West Liaohe River Plain.
The comparison of daily flow rate of Taihekou (Xin) station and Talagan station.
There are three large beside reservoirs (the total capacity of more than
In the intermittent river reach, a thick unsaturated zone formed beneath the riverbed. When the flow of water through the riverbed, seepage loss is very strong. There are two stages of seepage loss in the intermittent river reach, one is initial infiltration stage and the other is stable infiltration stage. Under the initial infiltration stage, the infiltration rate is large, but the decay rate is faster, and the cumulative infiltration of this stage accounted for a large proportion (Cheng et al., 2015).
Due to the difficulty of accurately dividing the seepage stage, the two stages of seepage loss are considered as a whole. From 2006 to 2013, only 3 years (summer flood season in 2011 and 2012, spring flood season in 2013) that the runoff diverted into Xinkai River by Taihekou Water Diversion Project can reached Talagan Hydrological Station. The daily average flow rate of Taihekou (Xin) station and Talagan station is compared (Fig. 4).
By statistics, the runoff of Talagan station accounted for only 20.38 % of the Taihekou (Xin) station in summer flood season in 2011, and the same ratio is 20.07 % in summer flood season in 2012. This indicates that the water seepage loss from the Taihekou (Xin) station to Talagan station reaches about 80 % of the total runoff. In spring flood season in 2013, the runoff of Talagan station accounted for 49.34 % of the Taihekou (Xin) station, and the seepage loss is relatively small, mainly because of the seasonal frozen soil in the riverbed has not been fully melted.
The Muskingum method is used to calculate river water flow evolution in the West Liaohe River plain (Eq. 1). The water flow evolution parameters are the results (Table 2) of the “West Liaohe River flood control planning report” (Northeast investigation and Design Institute of Ministry of water resources, 2005).
Considering the optimization calculation is too large, and the degree of concern for water quantity is higher than the flow process, simplified flow evolution and seepage loss coupling and used the comprehensive seepage loss coefficient. However, it is necessary to distinguish between spring flood season and summer flood season, expression is as flow:
To verify the comprehensive river seepage loss coefficient, the measured data of Daxingye station and Sanhetang station on the Xinkai River were selected from 1999 to 2000. The reach length is 36.7 km between the two hydrological stations. The comprehensive loss coefficients of spring and summer flood season were 0.78 and 0.65 respectively using Eq. (3). However, the measured values were 0.72 and 0.69 respectively, so the calculated result is close to the measured value. It is proved that the assumption of linear correlation between the loss coefficient and the length of the reach is reasonable. It is needed to be explained that the comprehensive seepage loss coefficient of reach was obtained under the condition that the river dried up for a long period previously. So the river water conditions will have a certain impact on the coefficient.
The distribution of river flow is controlled by water diversion project. Through the water diversion project, river flow can be diverted into beside reservoir or lower reach. The hydraulic connections which take water diversion project as a node are described below.
According to the analysis results of hydraulic system, the system is a mixed
type system and the hydraulic connection between upstream and downstream is
closely. The distribution of water resources were controlled by four water
diversion projects, so the dimension of the system space is four dimensional.
Based on the results of hydrological statistical analysis, the Water
Diversion Project of Talagan and Zongban had no outflow in recent 8 years
(from 2006 to 2013). So we use regular operation mode for the Water Diversion
Projects of Talagan and Zongban, and the two Water Diversion Projects no
longer participate in optimal operation. Specific to Zongban Water Diversion
Project as an example:
When When
After using regular operation mode for the Water Diversion Projects of
Talagan and Zongban, the four dimensional optimization problem can be
reduced to two dimensional optimization problem, and only the Water
Diversion Projects of Taihekou and Sujiapu can be seen as targets of the
optimal operation. The outflow and diversion flow are optional state
variables for every water diversion project. In order to meet the threshold
setting and the realization of the target, the diversion flow is chosen as
the state variable in Taihekou Water Diversion Project and the outflow is
chosen as state variable in Sujiapu Water Diversion Project.
The inflow of Taihekou Water Diversion Project is mainly concentrated in the
spring flood season (from 10 March to 20 April) and summer flood season
(from 25 June to 20 August). Other times the riverbed is in the drying
up state. The maximum inflow were 147 and 1074 m
Because the water diversion project has no capacity for runoff regulation,
the inflow either diverted into beside reservoir or discharged into lower
reach. Due to the extreme shortage of local water resources, the water
eventually was diverted into beside reservoir regardless of diversion flow or
outflow. Therefore, it is a prerequisite to optimize the operation of water
resources to meet the requirements of the local industry and agriculture.
Based on the prevention and control of wind erosion of riverbed, it is
beneficial that the river discharge or dry up for a long time. If the
upstream flow is less, the flow is firstly diverted into nearest beside
reservoir and reduce interference on the downstream of the river. If the
upstream flow is larger, the flow is diverted into long distance beside
reservoir preferentially, concentrated flow and reduce the frequency of water
through of riverbed. Therefore, as the objective function, it is needed to
consider diversion flow or outflow and river length. Specific expressions are
as follows:
The objective function includes two aspects: (1) the objection of water
volume that diverted into beside reservoir and (2) the objection of wind
erosion prevention and control of riverbed. The objective function is
essentially a function of the reach length (
The lower threshold flow is one of the key variables of the objective function. Through this variable, it can avoid the interference of high-frequency of small flow rate to the lower reach. To obtain the lower threshold flow, we can use two methods of optimization and simulation. If the optimization method is used to solve the problem, it is necessary to add two dimension on the basis of the original two dimensional optimization problem; and the use of simulation method to solve, we can learn from the operation experience, given a reasonable lower threshold flow. This paper uses simulation method to obtain the lower threshold of diversion flow of Taihekou Water Diversion Project and outflow of Sujiapu Water Diversion Project respectively. The lower threshold flow setting is directly related to the volume of diversion flow of Talagan and Zongban Water Diversion Project. The optimal value of the lower threshold of the diversion flow is based on the measured water diversion of the Talagan and Zongban Water Diversion Project.
In this paper, the problem of optimal flood control is not studied, and the flood discharge of each river is given in the form of constraint.
In the process of water resources optimal operation, the constrained
optimization problem is transformed into unconstrained optimization problem
by penalty function generally. In this paper, the discharge capacity and the
upper bound of reservoir water storage are all mandatory constraint. When the
variables violate the constraints, the penalty term is added on the basis of
the results of the original objective function. Specific expressions are as
follows:
Water diversion project itself does not have the ability to regulate runoff, only transform the inflow into diversion flow or outflow. Importantly there is a time effect of water flow evolution, that is, there is post-efficiency, which is essentially different from reservoirs. Therefore, the solution method of the optimal operation of the water diversion project group is quite different from the reservoir group.
In order to simplify the problem solving, multi-dimensional optimization problems are usually transformed into one-dimensional optimization problems, such as DPSA (Dynamic Programming Successive Approximation) in dynamic planning (Opan, 2010; Yi et al., 2003). Multi-stage optimization can also be translated into a single two-stage optimization problem, such as POA (Progressive Optimality Algorithm) in dynamic plan (Liu et al., 2011; Guo et al., 2011). If DPSA and POA are combined (DPSA-POA), the multi-dimensional and multi-stage optimization problems can be transformed into one-dimensional and two-stage optimization problems (Bai et al., 2015; Zhang et al., 2016), and the difficulty of solving optimization problems reduced greatly, at the same time, the problem of post-efficiency can be solved effectively.
To solve the problem of the water resources optimal operation based on riverbed wind erosion control in West Liaohe River plain, and improve computational efficiency, this paper presents a “multi-core parallel solution algorithm for optimal operation of water diversion project group”.
From the above analysis, we can see that the upper and lower reaches are
closely linked to each other of the optimal operation system in this paper,
and two independent state variables find the optimal combination of states
through the DPSA method under the framework of POA to get the optimal
operation scheme. According to the relationship between upstream and
downstream, Taihekou Water Diversion Project is the first level, Sujiapu
Water Diversion Project is the final level. Under the framework of POA, we
only change the state variable of the last stage
From the entire calculation process, the calculation is the largest in the final cycle. It is necessary to calculate the evolution of water flow under fixed stages of two water diversion project. It is necessary to calculate the amount of diversion flow and out flow of each water diversion project and calculate the entire objective function value. In this paper, the last stage is changed into parallel computing mode to balance the load of each thread, so as to avoid the cost of thread synchronization at the end of the cycle. The Fork/Join framework is a framework for performing parallel tasks. In the final cycle, each subtask is independent of each other, so the task can be distributed to multiple threads, and the results of multiple subtasks can be combined into a total calculation result.
The performance of parallel algorithm is evaluated by speedup ratio and
parallel efficiency (Peng et al., 2014), and the expression is as
follows:
Taking Taihekou Water Diversion Project daily measured flow from 2006 to 2013 as input condition, on the basis of guaranteeing the total water separation of various reservoirs, avoid the destruction of riverbed vegetation by the frequent overflow of the river through the optimal operation of Taihekou and Sujiapu Water Diversion Project, and achieve the goals of riverbed wind erosion control. The combination of the different thresholds of water diversion at the Taihekou and the Sujiapu Water Diversion Project are analyzed. Take Talagan and Zongban Water Diversion Project measured water diversion volume value as the verification. The results are shown in Table 3.
The flow chart of the modeling procedure.
Analysis of threshold of diversion flow at Taihekou Water Diversion Project and threshold of outflow at Sujiapu Water Diversion Project. Values in bold indicates that Scheme 2 is the best scheme.
Comparison of simulated and measured diversion flow into Xinkai River at Taihekou Water Diversion Project.
It can be seen from the calculation results in the table that the lower
threshold of water diversion for Taihekou Water Diversion Project is set to
150 m
Under scheme 2, Taihekou Water Diversion Project has diverted into Xinkai
River for 3 years, summer flood in 2008, summer flood in 2011 and summer
flood in 2012 respectively. The number of draining days at Talagan Water
Diversion Project is 29 days. In the measured value, Taihekou Water Diversion
Project has diverted into Xinkai River for 4 years, spring flood in 2006,
summer flood in 2011, summer flood in 2012 and summer flood in 2013
respectively. The number of draining days at Talagan Water Diversion Project
is 44 days. The results of the comparison between the calculated and the
measured are shown in Fig. 6. In the spring flood season of 2006, the maximum
inflow of Taihekou Water Diversion Project was 132 m
Under scheme 2, Zongban Water Diversion Project has inflow water for 7 years and the total number of inflow days is 176 days. In the measured values, Zongban Water Diversion Project also has inflow water for 7 years and the total number of inflow days is 200 days.
The results of comparison between the optimized operation and the measured data show that the number of years of water passing through the Zongban Water Diversion Project is the same with the measured value and the number of years that the Taihekou Water Diversion Project diverted into Xinkai River decreased by 1 year after the optimal operation is implemented. For the number of draining days, the number of draining days at Talagan and Zongban Water Diversion Project decreased to some extent contrast with the measured value, but the reduction in the number of draining days at Talagan Water Diversion Project was relatively large.
Through the above analysis we can see that the setting of the lower threshold of diversion flow and the amount of water diverted into the reservoir are conflicting. On the basis of guaranteeing the amount of water diverted into the reservoir, setting a reasonable lower threshold of flow diverted into the new Xinkai River of Taihekou Water Diversion Project and lower threshold of outflow of the Sujiapu Water Diversion Project can make the Xinkai River from Taihekou to Talagan Water Diversion Project (river length of 85 km) and the West Liaohe River from Sujiapu to Zongban Water Diversion Project (length of 65.7 km) through a larger water flow and control the lower flow in the reach. Reduce the frequency of river flow as much as possible, and reduce the frequent damage of water flow to the vegetation cover of the riverbed. The goal of water resources optimization of West Liaohe River Plain based on the prevention and control wind erosion of riverbed has been achieved.
By optimizing the operation, we can summarize the following rules: (1) When
the inflow of Tahekou Water Diversion Project is greater than or equal to
150 m
The optimized operation solution test environment is a notebook computer. CPU type is CORE i7, dual-core, 4G memory. After calculation, the serial algorithm takes 14.61 min and the parallel algorithm takes 9.66 min. The speed ratio of the parallel computation is 1.51, and the parallel efficiency is 0.76. It can be seen that the parallel computation of multi core shortens the time consuming and improves the computational efficiency.
According to the measured data and actual operation, the West Liaohe
River Plain water conservancy project system is simplified to be operation
system that is composed of the four water diversion project and three
reservoirs. On the basis of this, through the rules and regulations of
Talagan and Zongban Water Diversion Projection, the four-dimensional
optimization problem is further reduced to a two-dimensional optimization
problem. An optimal water resources operation model of West Liaohe River Plain
based on the prevention and control of wind erosion of riverbed was
established. The model objective function involves the factors such as the
amount of water diverted into the reservoir, the length of the river reach
and the lower threshold of the flow rate, and on the basis of ensuring the
requirements of water diversion in each reservoir, make the river flow
through the reach in a short period of time and avoid the destruction of
vegetation above the riverbed by the frequent overflow of the river. Proposed a multi-core parallel solution method of water resources
optimal operation in West Liaohe River Plain. Find the optimal combination
of states by DPSA method under the framework of POA and adopting the
Fork/Join mode design parallel computing algorithm in OpenPM programming
model, we proposed the calculation flow. The calculation results show that on the basis of ensuring the demand of
water diversion into Talagan and Molimiao reservoirs, the threshold of
diversion flow from Taihekou Water Diversion Project to Xinkai River is
150 m The research results of this paper show that it is feasible to achieve
the prevention and control of wind erosion of riverbed in the intermittent
over-water reach of West Liaohe River Plain through the optimized water
resources operation.
The hydrological data used in this study all come from the hydrologic Yearbook (Hydrology data of Liaohe River Basin). According to the People's Republic of China hydrological regulations, hydrological data are not allowed to be made public.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Innovative water resources management – understanding and balancing interactions between humankind and nature”. It is a result of the 8th International Water Resources Management Conference of ICWRS, Beijing, China, 13–15 June 2018.
The study is financially supported by the National Non-Profit Research Program of China (No.201401015). Edited by: Depeng Zuo Reviewed by: two anonymous referee