Next Article in Journal
Spatial Variability of Preferential Flow and Infiltration Redistribution along a Rocky-Mountain Hillslope, Northern China
Previous Article in Journal
Drought Trends and the Extreme Drought Frequency and Characteristics under Climate Change Based on SPI and HI in the Upper and Middle Reaches of the Huai River Basin, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Joint Optimal Allocation Study of Regional Total Water Consumption and Pollutant Carrying Capacity of Water Function Areas Based on Emergy Theory

1
School of Water Conservancy Engineering, Zhengzhou University, Zhengzhou 450001, China
2
School of Civil Engineering, Zhengzhou University, Zhengzhou 450001, China
3
Kaifeng Hydrology and Water Resources Survey Bureau, Kaifeng 475000, Henan, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(4), 1101; https://doi.org/10.3390/w12041101
Submission received: 19 February 2020 / Revised: 2 April 2020 / Accepted: 10 April 2020 / Published: 13 April 2020
(This article belongs to the Section Water Resources Management, Policy and Governance)

Abstract

:
How to allocate the limited water consumption and pollutant carrying capacity of water function areas to sub-areas and water sectors is key for the practice of the most stringent water resource management decision-making. In response to the joint optimal allocation mechanism of water consumption and pollutant-carrying capacity, joint optimal allocation and the unified quantification problem of multiform benefits, the whole process of joint optimal distribution is refined into the links “water supply—water use and consumption—drainage—pollutant carrying” and the unified quantification method based on the emergy theory to achieve various economic, social, and ecological benefits in different links is proposed. Then the joint optimal allocation joint optimal allocation model of regional water consumption and pollutant carrying capacity is constructed, which realizes the reverse constraint relationship between water quantity and pollution carrying capacity in the water quality and quantity conversion processes. Using Zhengzhou city as an example, the models and methods were applied and researched. Countermeasures and suggestions for the rational allocation of available water resources and water pollution capacity were put forward accordingly.

1. Introduction

According to the most stringent water resource management system, the total amount of control exercisable over water resources and the benefit of such control over water use and the pollutant carrying capacity limits of the water function area must be implemented simultaneously in the process of water resource utilization and protection [1,2,3]. Given the whole process of “water supply—water use and consumption—drainage—pollutant carrying,” the available water quantity is closely related to the pollutants in the river, and the total amount of pollutants in the river is interrelated with sewage discharge. At the same time, sewage discharge is also related to the allocation of domestic water, production water, and ecological water, as well as the benefit of using water. This requires that the water distribution, water function zone assimilative capacity distribution, and water use benefit improvements be unified across water resources in order to optimize allocation.
Therefore, the joint optimal allocation of water consumption and pollutant carrying capacity is essentially a manual regulation that realizes the reverse constraint relationship between water quantity and pollution carrying capacity in the water quality and quantity conversion processes in order to control the water consumption quantity from the source water supply link and the total pollutants amount into the river from terminal drainage link, and meet the requirements of available water supply of water sources and pollutant carrying capacity of the water functional area at the same time. Then, the water use efficiency is improved and the reasonable distribution of the limited water consumption and pollutant carrying capacity is realized.
Additionally, several imperfections exist in the current related research on the joint optimal allocation of water quality and water quantity. Yan et al. [4] proposed the technical framework for a dual cross control study of rivers based on the reasonable allocation of water resources from the water quantity control layer and the water quality control layer. However, this study only meets the total pollutant carrying capacity constraints of the water function area from the point of total amount control. It is difficult to guarantee the pollutant carrying capacity restrictions of water function zoning. Liu et al. [5] regarded both water withdrawal and pollutant discharge as control indexes, and constructed the multi-objective model for optimal allocation of water resources in based on dual-control of water quality and quantity. However, this water quality target is difficult to ensure that the minimum pollutants discharge into the river is within the pollutant carrying capacity range of the water function area. You et al. [6] summarized the regulation idea and research progress of unified regulation of water quantity and quality and pointed out that there are few research studies on unified regulation of water quantity and quality aiming at decision support at present, and put forward the total quantity regulation idea aiming at pollution control, water quantity allocation, and water ecological protection. Fan et al. [7] proposed the overall idea that the study water demand threshold through total pollution carrying capacity control, starting from the carrying capacity of water environment, which deduced the maximum water demand of the basin, according to pollution carrying capacity. However, as long as the water requirement threshold is met, it is regarded as the control index of total pollution carrying capacity. Zhang [8] took water resources utilization, water efficiency, and total pollution discharge into the river as a control target simultaneously, and constructed the optimization allocation of water quantity and quality based on life cycle theory of water resources. However, it cannot realize the allocation of pollutant carrying capacity in each sub-area as the same as the water quantity element. The ecological and economic impacts after pollutants into the river are not considered in the target selection. Therefore, the joint distribution benefits cannot be reflected fully.
In summary, the current state of integrated optimal allocation of water quality and water quantity still exists the following problems that need further study: (1) Configuration elements. The normal configuration mainly focuses on the water quantity allocation among different users, while water-quality factors have generally been simplified. However, it is difficult to ensure that pollutant emissions meet the pollutant carrying capacity limits of each water function area, and lacks the handling analysis for pollutant reduction and the mutual restriction mechanism study of water quantity and pollutant carrying capacity. (2) Configuration model. Configuration mainly uses the multi-objective model, but it is difficult to achieve the target of economic, social, and ecological environments—that is, for a uniform measurement—because they have different forms and dimensions. Thus, the weighted method is usually adopted to transfer the multiple objectives into a single objective, but its subjectivity is too strong to ensure that the configuration results accord with reality.
To sum up, the purpose of the work is to realize the joint distribution of total water consumption and pollutant carrying capacity of water functional areas from the perspective of eco-economics. Therefore, this paper builds the coupling relationship among water consumption, sewage discharge, and pollutants discharged into rivers, proposes the expression form and calculation method of various benefits in “water supply—water use and consumption—drainage—pollutant carrying” links, which simultaneously takes water consumption and pollutant carrying capacity in water functions as decision-making variables, and achieves the maximum net benefit in the process of joint distribution as an objective function in order to establish the joint optimal allocation model of the total water consumption and pollutant carrying capacity based on emergy theory.

2. Material and Methods

2.1. The Emergy Analysis of the Benefits of Joint Optimal Allocation of Water Consumption and Pollutant Carrying Capacity

The water resource eco-economic system transforms the water resources system, ecological environment system, and social economic system into an indispensable organic component of its internal structure and operation. Then, water resource utilization and protection, and the organic connection between water resources and other economic and natural resources, are emphatically focused on. However, as a multi-objective service for the natural and social economy, water resources not only support economic and social development, but also satisfy the stability of the ecological system [9]. The benefits of integrating the distribution of water consumption and pollutant carrying capacity are accompanied by the circulation process of “water supply—water consumption—drainage—pollutant carrying.” Therefore, in contrast to the target benefit of traditional water allocation methods, the benefits of joint optimal allocation not only needs to consider the economic, social, and ecological environment benefits of the water consumption link, but must also consider the environmental damage caused by sewage and the cost of maintaining ecological health in drainage and pollutant carrying links.

2.1.1. The Component of Joint Optimal Allocation Benefits

(1) Benefits of the water supply link
In the water supply link, dams and ponds are built to store water resources, reservoirs, and diversion projects that are built to develop water resources. Water plants and purified water plants are built to purify water bodies. In summary, the main purpose of the water supply link is to introduce natural water bodies into socio-economic systems, and to form engineering water bodies that are affected by construction works, which can be used by the water consumption link. Therefore, the benefits of the water supply link are shown in the supply of engineering water, which is used by various water use sectors, and belong to the category of ecological environmental benefits.
(2) Benefits of water use and the consumption link
In the water use and consumption links, as production materials, water resources provide power and raw materials for economic and social production, such as agricultural water, industrial water, hydroelectric power generation, shipping, and so on, and reflects the positive economic benefits of water resources utilization. As domestic materials, water resources provide drinking water for human beings and other animals in order to maintain domestic health, which brings the positive social benefit of water resources utilization [10]. As ecological materials, water resources are used for urban greening, environmental sanitation, river and lake hydration, and plays an irreplaceable role in regulating the climate, purifying the environment, fixing carbon, releasing oxygen, improving social welfare, and more. This creates the positive ecological and environmental benefit of water resource utilization. Therefore, due to the different sectors of water use, the corresponding forms of various benefits are different, and it is difficult to achieve a unified measurement.
(3) Benefits of the drainage link
In the drainage link, the sewage produced by agricultural water and domestic rural water are discharged into the natural ecosystem directly. It not only reduces the water quality and destroys the water environment, but also makes the originally available water resources lose their use value [9], which reflects the negative ecological-environment effects. The sewage from domestic urban water and industrial water are usually discharged into the natural water body or are reused after the treatment of sewage treatment projects. This process requires a lot of economic investment, such as energy consumption, pharmaceutical consumption, maintenance costs, personnel costs, and more, and brings corresponding negative economic effects.
(4) Benefits of the pollutant-carrying link
In the pollutant-carrying link, although sewage treatment in most regions of China have implemented the discharge standard, the amount of pollutants discharged into rivers after treatment still cannot meet the requirements of the restriction on the amount of pollutants in a regional water function area. On the one hand, we can add new sewage treatment systems in order to reduce the pollutants that exceed the pollutant carrying capacity of water function areas [11], so as to reduce the chemical energy of sewage. On the other hand, we can start from the source and control the water supply quantity of water sectors, which reduces the sewage amount directly. Lastly, it is also possible to improve water use efficiency, adjust the industrial structure, and reduce the pollutant discharge coefficient. The economic investment of these pollutant-carrying measures has negative economic effects. In summary, in the process of integrating the distribution of water consumption and the pollutant carrying capacity, the water cycle of different links in different water sectors show the different types and forms of benefits, which cover the positive and negative economic, social, and ecological environment effects. They are independent and self-contained, but interact with each other due to the circulation and flow of water.

2.1.2. Emergy Quantification Method for Joint Optimal Allocation Benefits

The benefits analysis of the integrated optimal distribution of regional water consumption and pollutant carrying capacity shows that the forms and dimensions of benefits are different. It is difficult for the traditional evaluation method to achieve a unified quantitative evaluation. Emergy analysis is a quantitative research method in the process of the ecological economic system, which can transform all kinds of materials with different levels of energy and money being put into the solar energy value, and the unified measurement of various materials and energies is realized [12]. Emergy is the energy amount of another category contained in the fluid or stored energy. There is a transformational relationship between energies of different qualities—emergy transformity, namely the contained amount of another energy that forms per unit of matter or energy [9]. Due to various energies being derived from solar energy directly or indirectly, the emergy analysis converts different types of energy that cannot be compared into the same standard. Solar emergy’s unit is solar emjoules. The abbreviation is sej.
Therefore, environmental factors, production factors, technological factors, and economic indicators can be transformed into the same “emergy”, and emergy analysis is applied to calculate all the integrated optimal distribution efficiencies, either directly or indirectly.

2.2. The Joint Optimal Allocation Model of Water Consumption and Pollutant Carrying Capacity

2.2.1. Decision Variables

In contrast to the traditional optimization allocation of water resources, the decision variables of macroscopic integrated optimal distribution of regional water consumption and pollutant carrying capacity includes two types: first is the quantity of water from different water sources that are allocated to different users. The other type is the pollutant-carrying capacity of different sub-regions by different users, which are allocated to different water function areas.
(1) The quantity of water from different water sources that are allocated to different sub-regions and different water sectors. X was used, respectively.
X = ( X 1 , X 2 , , X S L ) X k = ( x i j k )    ( i = 1 , 2 , , P W Y + W Y   ( k ) ;   j = 1 , 2 , , W S   ( k ) ;   k = 1 , 2 , , S L )
x i j k , represents the amount of water of user j in sub-region k, which was allocated from the water source i.
(2) The pollutant carrying capacity of different sectors in different sub-regions, which were allocated from different water function areas. Y was used, respectively.
Y = ( Y 1 p , Y 2 p , , Y S L p ) Y k p = ( Y q j k p )    ( q = 1 , 2 , , P W F + W F   ( k ) ;   j = 1 , 2 , , W S   ( k ) ;   k = 1 , 2 , , S L )
y q j k p , represents the controlling quantity of the pollutant into the river of pollutant p of user j in sub-region k , which was allocated from water function area q .

2.2.2. Objective Function

Based on the positive and negative effects generated by the links of “water supply—water use and consumption—drainage—pollutant carrying,” the final objective function of macroscopic integrated optimal distribution of regional water consumption and pollutant-carrying capacity is the maximum net benefit in the process of integrated optimal distribution. Among them, as the intermediate product of the system, the benefits in the water supply link are utilized by the next link. In order to avoid repetition, it is not included in the final net benefit [22].
Therefore, the objective function can be expressed as the positive benefits in water use and consumption link, by deducting the negative effects in the drainage link and pollutant-carrying link, which is the net benefit of joint optimal allocation, which reaches its maximum. Through the emergy calculation methods of all the benefits mentioned above, the objective function of joint optimal allocation can be expressed as follows.
F = max k = 1 S L j = 1 W S ( k ) f j k = max ( k = 1 S L j = 1 W S ( k ) E U j k ( x ) k = 1 S L j = 1 W S ( k ) E D j k ( x ) m = 1 p k = 1 S L j = 1 W S ( k ) E P j k p ( y ) ) .
In this formula, E U represents the positive economic, social, and ecological benefits of water use and the consumption link in various water sectors, including the industrial water use benefits, agricultural water use benefits, domestic water use benefits, and ecological–environmental water use benefits. E D represents the negative environmental benefits in drainage link, which consists of the sewage treatment loss caused by industrial water and urban domestic water, and the water quality pollution loss caused by the direct discharge of agricultural drainage water and rural domestic sewage. E P represents the negative economic effects of the pollutant-carrying link, including the pollutant reduction costs in industrial water and urban domestic water.
According to the emergy evaluation methods, which have varying benefits, the concrete calculation measures are as follows.
{ f 1 = E M p i E M i s E M i p        ( j = 1 , represents   the   water - use   department   in   industry ) f 2 = E M p a E M a d        ( j = 2 , represents   the   water - use   department   in   agriculture ) f 3 = E M l r E M d d E M d s E M d p     ( j = 3 , represents   the   water - use   department   in   life ) f 4 = E M d c + E M c h + E M c o + E M s d    ( j = 4   , represents   the   water - use   department   in   ecological ) }
In this formula, each symbol has the same meaning as above.

2.2.3. Constraint Condition

(1) Total water consumption constraint
Following the total quantity red line constraints from the development and utilization of water resources, regional water consumption should not exceed the control objective from the total regional water use.
k = 1 S L i = 1 P W Y + W Y ( k ) j = 1 W S ( k ) x i j k T W .
In this formula, x i j k represents the amount of water of user j in sub-region k, which was allocated from the water source i, and T W represents the control objective of the total regional water consumption.
(2) Available water resources constraint
a. Public water supplies: The sum of the total water resources of each sector in the relevant sub-area, supplied by the public water supplies, is not more than the quantity of the available water resources for allocation.
k = 1 S L j = 1 W S ( k ) x i j k G Y W i ,     i P W Y .
In this formula, P W Y means the collection of public water supplies, and G Y W i represents the available water resources in the public water supplies.
b. Special water supplies: The total water resources of the relevant sector, supplied by special water supplies to sub-areas, does not exceed the amount of water resources available for allocation.
Surface water sources
j = 1 W S ( k ) x i j k W S i k ,     i S W Y ( k ) , k = 1 , 2 , , S L .
Groundwater sources
j = 1 W S ( k ) x i j k W G i k + Δ W i k ,    i G W Y ( k ) , k = 1 , 2 , , S L .
In this formula, W S i k and W G i k represent the amount of water resources that are available for allocation in the i special surface water and groundwater sources of sub-region k, respectively. S W Y ( k ) and G W Y ( k ) represent the collection of special surface water and groundwater sources of sub-region k, and Δ W i k is the amount of over-exploited groundwater that is allowable in sub-region k.
(3) Water conveyance capacity constraint
The total amount of water from all the sources that are allocated for use by sector j in sub-region k cannot exceed the maximum carrying capacity of the water transfer project in use in sector j from relevant sources.
x i j k Q i k max ,    i P W Y W Y ( k ) , k = 1 , 2 , , S L .
In this formula, Q i k max is the maximum carrying capacity from source i to sub-region k in water use sector j.
(4) Water use efficiency constraint
Following the red line of controlling over water use efficiency, the quota of the water requirement of industry, agriculture, and national economy should be less than the control index examined for water use efficiency.
D i D i d , D a D a d , D e D e d .
In this formula, D i , D a , and D e represent the quota of the water requirement for industries, agriculture, and the national economy, respectively. D i d , D a d , and D e d are the control indexes of water use efficiency for industries, agriculture, and the national economy, according to the red line of controlling water use efficiency, respectively.
(5) Pollutant carrying capacity of water function areas constraint
Pollutant carrying capacity constraints of the water function area means the total pollutant carrying capacity of all the water use sectors cannot exceed the maximum amount of pollution that is allowed in the water function area.
k S L j G N W S ( q ) y q j k p W q p .
In this formula, y q j k p represents the pollutant-carrying capacity of pollutant p for user j in sub-region k, allocated from the water function area q. G N W S ( q ) represents all the water use sectors related to water function area q in this region. W q p represents the maximum amount of pollutant p that is allowed in the water functional area q, and other symbols have the same meaning as above.
(6) Mutual relation between water quantity and the pollution carrying capacity constraint
By using the pollutant-carrying capacity allocated to each water sector, the amount of pollutants discharged by pollutant discharging units and its allocated water amount are constrained. That is, the amount of pollutants discharge under a certain amount of allocated water conditions by each water sector in each sub-area, after the status quo treatment and subsequent reduction links, shows that the final pollutants discharge into the water body cannot exceed the controlled amount of pollutants into the river allocated to this water sector by the corresponding water functional area.
x j k y q j k p + λ q j k p c y q j k p λ q j k p α j [ P C j k p η j k ( P C j k p T P C j k p ) ] .
In this formula, P C j k p is the concentration of the pollutant p in the wastewater discharged by water use sector j in sub-region k. T P C j k p represents the concentration of the pollutant p by processing in the wastewater discharged by the water used by sector j in sub-area k. c y q j k p represents the quantity of pollutant p needed to be reduced in the sub-region k water use sector j before being discharged into water function area q, t. α j is the sewage discharge rate of user j. η j k is the treatment rate for wastewater of user j in sub-area k, and λ q j k p represents the coefficient in the river of pollutant p discharged into the water function area q by the sub-region k of water use sector j. Other symbols have the same meaning as above.
(7) Pollutant reduction funding constraint
The investment scale of regional pollution control is affected by economic and social development and is limited to a certain period. The funds put into pollutant reduction measures should meet the funding constraints.
S D C k S max   k
In this formula, S max   k represents the total amount of pollutant reduction funds in the corresponding period.
(8) Other constraints
Minimum water consumption constraints in various sectors (industrial, agricultural, domestic, and ecological water sectors) are specifically set according to regional characteristics.
(9) Variable non-negative constraint
The above objective functions and constraints constitute the joint optimal allocation model of regional water consumption and pollutant-carrying capacity.

3. Results

3.1. The Construction of the Joint Optimal Allocation Model of Water Consumption and Pollutant-Carrying Capacity in Zhengzhou City

3.1.1. The Sub-Area Division and Water Source and User Composition

(1) Sub-Area Division
Based on «The comprehensive planning of water resources in Zhengzhou», Zhengzhou city is divided into seven sub-districts (SL = 7), as administrative districts, including Zhengzhou district, Dengfeng city, Xinzheng city, Xinmi city, Xingyang city, Gongyi city, and Zhongmou county.
(2) Water functional zoning
The main water functional areas in Zhengzhou are shown in Table 2.
(3) Water source and user classification
Combined with the actual situation and planning results of water resources in Zhengzhou, the common water sources include the yellow river water diversion and south-to-north water diversion, which can supply water to any subarea within the region during water distribution. Special water sources include the water storage project, water diversion project (excluding the yellow river water diversion), and groundwater by reclaiming water and rainwater. The water use sectors in each sub-district are divided into industrial water use, agricultural water use, urban domestic water use, rural domestic water use, and ecological environmental water use.

3.1.2. Scene Settings of Joint Optimal Allocation

In terms of the water supply source, water distribution engineering, water use sector, wastewater treatment, water function areas, and the engineering and non-engineering measures, which have influence over water distribution, pollutant-carrying capacity distribution, and water use efficiency, are comprehensively analyzed. In accordance with the principle of feasible technology and possible economic benefit, the joint optimal allocation scenes of water consumption and pollutant-carrying capacity are constructed. The scene mainly refers to the combination of different control measures that affect joint optimal allocation results in order to analyze the mutual restriction relationship between water quantity and pollution carrying capacity intuitively, so as to obtain the joint optimal allocation situations of different combinations. The scene settings mainly consider the aspects of water quantity, pollutant-carrying capacity, and integrated regulation.
(1)
Water demand side: From the perspective of water saving measures, improving water use efficiency, and analyzing the influence of the amount of water in the joint optimal allocation.
(2)
Sewage treatment side: Starting from the sewage treatment level, analyze the influence of current and advanced treatments on the results of the joint optimal allocation.
(3)
Pollutant reduction funds side: From the perspective of pollutant reduction capital investment, analyze the limitation and influence of different capital investment levels on the water use structure of each sector.
On this basis, combined with the water conservation development planning and socio-economic development planning in Zhengzhou, the above regulation measures were combined in various ways. Then, the scene settings of the joint optimal allocation in Zhengzhou was constructed, and the configuration effects of different measures were compared, as shown in the Table 3.

3.2. The Results of Joint Optimal Allocation of Water Consumption and Pollutant Carrying Capacity in Zhengzhou City

3.2.1. The Prediction of Water Supply and Demand of Zhengzhou in 2030

According to the actual situation of Zhengzhou city, the moderate water-saving and intensified water-saving scenes are put forward in 2030. Referring to the water resources evaluation and other planning achievements of Zhengzhou, the water demand prediction results in 2030 are shown in Table 4 and Table 5, and the available water supply prediction result of Zhengzhou in 2030 is shown in Table 6.

3.2.2. Calculation of the Emergy Parameters of Water Resources

Combined with the basic data of water resources, economy, society, and ecology in ‹‹the statistical yearbook of Zhengzhou city›› and ‹‹The water resources bulletin of Zhengzhou city››, the emergy parameters involved in joint optimal allocation process are calculated by the emergy calculation method, which mainly include emergy transformity of different water bodies and emergy parameters of water sectors. With regard to the reference value for planning annual emergy parameters, the results are shown in Table 7, Table 8, Table 9, Table 10, Table 11 and Table 12.

3.2.3. Calculation of the Pollutant Carrying Capacity Parameters of the Water Function Area

According to the water quality objectives determined in ‹‹The comprehensive planning of water resources in Zhengzhou city››, and combined with the comprehensive hydrological conditions and attenuation coefficient of Zhengzhou, COD and NH3-N are determined to be the key factors for the pollutant carrying capacity calculation. The calculation results of Zhengzhou in 2030 are shown in Table 13.
On this basis, combining with the statistical data of sewage treatment quantity, treatment cost and total investment of the sewage treatment project in Zhengzhou city, the unit sewage treatment costs of Zhengzhou are determined: the industrial wastewater treatment cost is 0.914 yuan/ton, and the urban domestic sewage treatment cost is 0.528 yuan/ton. Then, using the emergy calculation method of water resources eco-economic system in literature [9], the emergy transformity of water bodies before and after pollution in Zhengzhou are 6.26 × 1011 sej/m and 8.05 × 1012 sej/m, respectively.

3.2.4. The Integrated Optimal Distribution Results of Zhengzhou in 2030

Using the solution methods established above, the joint optimal allocation models under different scenes are solved, and the results of joint optimal allocation of water consumption and pollutant carrying capacity in Zhengzhou city are obtained. The joint optimal allocation result under each scene is called a scheme. Then, the detailed results of the joint optimal allocation under the 50% guarantee rate of Zhengzhou city in 2030 are analyzed, as shown in Table 14.

4. Discussion

4.1. Analysis of Water Consumption Distribution Results

(1) Analysis from the perspective of the water shortage
With the increase of various measures in each distribution scene, the total water shortage of Zhengzhou in 2030 decreased from 24,779 × 10 4 m3 (Scheme 1) to 1533× 10 4 m3 (Scheme 8), and the water shortage rate decreased from 10.9% to 0.74%, which alleviated the water shortage situation of Zhengzhou greatly. The change trend in the water shortage under the 50% guarantee rate of Zhengzhou in 2030 is shown in Figure 1.
Taking scheme 5 as an example, although the implementation of water-saving measures has alleviated the contradiction of the water supply and demand to a certain extent, the total water supply in Zhengzhou has not reached the maximum water supply capacity from the water supply source fully. This is because, under the certain circumstances of the sewage treatment level and pollutant treatment investment, the restraint of the pollutant carrying capacity in the water function area limits the water use quantity of each production and domestic sector directly. For scheme 7 and 8 (compared to scheme 5 and 6), due to the improvement in the sewage treatment level and pollutant reduction cost, the constraint in the pollutant reduction funding restrictions on the pollution discharge amount, by the high-polluting industrial water sector, is weakened to some extent. Water is preferentially distributed to the industrial water sector, with higher efficiency, for water consumption, which increased water shortages in the agricultural sector to some extent. However, with the increase in the scheme’s different measures, the water shortage of Zhengzhou city in 2030 presents a decreased trend integrally.
(2) Analysis from the perspective of the water use structure
As shown in Figure 2, domestic and ecological water is satisfied in different distribution schemes, which reflects the domestic water constraint and ecological objective proposal in the model’s settings and the calculation methods of the model’s solution are indirectly effective. In addition, if the pollutant reduction cost and pollutant carrying capacity limitation are not considered, industrial water use is given priority. However, for scheme 1 to scheme 8, due to the restrictions from the pollutant reduction costs and the pollutant-carrying capacity of the water functional area, water distribution decreases water consumption in industrial sectors in order to reduce pollutant discharge with a corresponding reduction in treatment costs. Thus, to achieve the goal of limiting the pollutant carrying factor in water functional areas, the reduced water amount is allocated to the agricultural water. Due to the increase in the allocation measures of each scheme, the ratio of the water content in each sector to the total water distribution amount also changes accordingly, which reflects the role of the model in targeting the setting and measure regulations.

4.2. Analysis of Pollutant Carrying Capacity Distribution Results

(1) Analysis from the perspective of pollutants into the river
The amount of pollutant in the river decreased from 23,997.2 t (COD) and 1479.4 t (NH3-N) in scheme 1 to 21,940.9 t (COD) and 530.8 t (NH3-N) in scheme 8, under the 50% guarantee rate of Zhengzhou in 2030, which shows a decrease trend on the whole, and meets the requirements of the pollutant-carrying capacity in the water functional area and each water use sector.
As the allocated water amount increases, the pollutant discharged amount of each water sector increases accordingly, but, with the increase in various sewage treatments and pollutant reduction measures in different schemes, the pollutant reduction capacity is improved significantly. Therefore, the final amount of pollutants into the river in each scheme also decreased significantly. The pollutants into the river in various sectors showed an overall decrease trend from schemes 1–8 (as shown in Figure 3, Figure 4 and Figure 5), which also verifies the mutual influence and restriction between water allocation measures and pollutant treatment measures in the model.
(2) Analysis from the perspective of the pollutant-carrying capacity allocation structure
Take scheme 1 and 2 as examples. Scheme 2 increases the pollutant reduction funding on the basis of scheme 1, which reduces the restrictions of water supply constraints for the high-polluting sectors. Industries, from pollutant reduction funding, means that water supply tends to be allocated to industries, whose unit water use efficiency is high. Therefore, the pollutant carrying capacity of industrial water sectors increased in comparison to agricultural and domestic water sectors, which reduced (as shown in Figure 6). To sum up, while increasing the water supply to solve water shortage, it is necessary to adjust the water use structure reasonably and to reduce water consumption in high-pollutant discharge sectors appropriately in order to control the pollutant discharge and the reduction amount directly. At the same time, it is necessary to increase the construction of sewage treatment facilities so that the pollutant to be treated can be treated effectively and the unity of water resource utilization and water environment protection can be achieved.

4.3. Analysis of Joint Optimal Allocation Results of Water Consumption and Pollutant-Carrying Capacity

According to the objective function and the calculation methods of various benefits in chapter 2, the joint optimal allocation benefits that are realized by different distribution schemes are calculated, as shown in Table 15. Since the dimensional level of the positive benefit and negative benefit are clearly different, the comparative analysis is carried out on them, respectively, as shown in Figure 7 and Figure 8.
As shown in Table 15 and Figure 7 and Figure 8, the joint optimal allocation benefits showed an overall increasing trend in Schemes 1–8. In order to compare the effects of various distribution measures on the joint optimal allocation results, the analysis is as follows.
(1) Analysis of the measures from the water demand side (schemes 5–8 vs. schemes 1–4)
Under the moderate water-saving scheme in 2030, the total water shortage of Zhengzhou is 19,552 × 104 m3. After the implementation of the strengthening water-saving scheme, the total water shortage in Zhengzhou was reduced to 1533 × 104 m3.
Under the circumstance that all other measures in schemes 1–4 remain unchanged, only strengthening water-saving is implemented in schemes 5–8. Then, on the one hand, the unit water efficiency is improved, the shortage in water supply is alleviated, and the water use guarantee rate of each sector is improved. On the other hand, due to the implementation of water-saving measures, the water supply amount of various sectors is reduced, which leads to the decrease in the sewage discharge amount and pollutant amount into the river in industrial and other high-discharge sectors directly, which results in a corresponding decrease in direct sewage discharge and sewage treatment costs. Therefore, the joint optimal allocation benefits in schemes 5–8 are significantly higher than that of schemes 1–4.
To sum up, the implementation of strengthened water-saving measures alleviate the contradiction between the supply and demand for water resources, while adjusting the water supply structure reduces the water supply of high-pollutant discharge sectors appropriately, reduces the final pollutant amount into the river indirectly, and decreases the allocation amount of the pollutant-carrying capacity. Therefore, water-saving measures are a powerful means to effectively alleviate the contradiction between supply and demand for water resources and the deterioration of the water environment.
(2) Analysis from the measures from the sewage treatment side (schemes 3–4 vs. schemes 1–2)
Scheme 3–4 improve the sewage treatment level on the basis of schemes 1–2 even though it increases the sewage treatment cost to some extent. On the one hand, it reduced the direct pollutant amount into rivers in industrial and domestic sectors greatly. On the other hand, it also reduced the pollutant reduction amount and had corresponding reduction costs. Therefore, the joint optimal allocation benefits still improved, correspondingly, from 1868.98 × 1019 sej (scheme 1) and 1873.14 × 1019 sej (scheme 2) to 1904.35 × 1019 sej and 1908.54 × 1019 sej in schemes 3 and 4, respectively. Thus, improving the sewage treatment level and paying attention to the construction of sewage treatment facilities can more effectively treat the pollutants discharged from production and living, realize more comprehensive water resources utilization benefits, and realize the coordination and unification of water resources utilization and water environmental protection.
(3) Analysis from the measures of the pollutant reduction funds side (scheme 7 vs. 8)
After the restriction on pollutant reduction investment is added, the pollutant reduction capacity is reduced indirectly, and the water allocation will reduce the quantity of water distribution in high-emission sectors correspondingly, in order to realize the constraint of pollutant-carrying capacity in water functional areas. Therefore, the quantity of water distribution in industrial water sectors, whose pollutant discharge amount per unit of water is the largest, has the clearest decrease. The reduced water amount is allocated to agricultural and other water sectors accordingly.
On the basis of scheme 7, the pollutant reduction funding level is improved in scheme 8, which reduces the restrictions on the water structure by reducing the funding limit, and makes the water supply more likely to be allocated to industrial and domestic water sectors, which leads to an increase in the pollution loss by direct discharge loss and sewage treatment cost, but the unit water use efficiency of the industrial and domestic water sector is significantly higher than the agricultural water sector. Therefore, the increasing trend of positive benefits in the water consumption link is clearer, and the joint optimal allocation benefit realized in scheme 8 is, in general, still higher than that in scheme 7.
In conclusion, with the model optimization objective as the criterion, scheme 8 is recommended as the optimal solution under the 50% guarantee rate in Zhengzhou in 2030.

4.4. Comparative Analysis of the Joint Optimal Allocation Results of Water Consumption and Pollutant-Carrying Capacity with the Traditional Allocation Method

To identify the feasibility and advantages of joint optimal allocation methods of water consumption and pollutant carrying capacity based on emergy theory, the above joint distribution results are compared with those obtained by traditional water resources allocation. The reverse restriction of pollutant carrying capacity limits to the water amount allocation and benefits evaluation is expounded. The necessity of joint optimal allocation is verified. The results obtained by traditional methods consulted the results in “Study on the water resources allocation based on secondary distribution in Zhengzhou metropolitan area,” [23] and the comparison and analysis between the joint optimal allocation results based on emergy theory and traditional results in 2030 are shown in Table 16.
Compared with the traditional water resources allocation, the joint optimal allocation based on emergy theory in this paper focus on the mutual restriction relations between water quantity and pollutants discharge into the river, realizes the reverse regulation of water consumption and pollutant carrying capacity in the joint distribution process, and ensure to obtain the most optimal joint allocation scheme on the premise of meeting the pollutant carrying capacity of water function areas. On the contrary, the traditional allocation of water resources only aims at the maximum benefit by ignoring the limitation of the pollution carrying capacity, which may not meet the needs of an ecological environment. Therefore, the satisfaction degree of pollution carrying capacity under the optimal water allocation of two methods is compared and analyzed. The results are shown in Table 16.
As shown in Table 16, the joint optimal allocation method in this paper limits from water supply and pollution carrying. The ultimate limit of pollution carrying capacity reverses the initial water supply, and then achieve maximum benefits of joint allocation in ensuring to meet the pollution carrying capacity constraints, and satisfy the efficient utilization of water resources and water environmental protection targets simultaneously. However, the traditional method takes no account of the limitation and regulation function of pollutant carrying capacity. It only regards the final biggest benefit as a target, so its allocation water quantity is significantly greater than the emergy method. Then the corresponding pollutant discharge of COD and NH3-N calculated according to the allocation water quantity are higher, far beyond the pollution carrying capacity of the water function area. This is not conducive for the sustainable development of water resources’ ecological economic system. Therefore, it is necessary to carry out the research on the pollution carrying capacity constraint and its joint allocation in the optimal allocation of water resources.

5. Conclusions

Based on the emergy theory of ecological economics, the thought frame, model method, and overall application of the joint optimal allocation of regional water consumption and pollutant carrying capacity are studied, and enriches existing research ideas and technical methods of the optimal allocation of water resources. The main conclusions are as follows.
(1) Ecological economics emphasizes the coupling of the natural ecological system and social economy system. This study is based on the emergy theory of ecological economics with regard to the joint optimal allocation process of regional water consumption and pollutant-carrying capacity as the material circulation and energy conversion process, and takes the water consumption and pollutant-carrying capacity as the direct configuration elements at the same time. This can reflect the restriction relationship and reverse the regulation mechanism between them profoundly, as well as the theoretical basis of sustainable utilization and management research of water resources.
(2) By analyzing the water quality and water quantity conversion process of “water supply—water use and consumption—drainage—pollutant carrying” links, the positive and negative effects, and their manifestations of water on the economic, social, and ecological environment, can be described comprehensively. The results show that the components of the joint optimal allocation benefits include: ① benefits of the water use and consumption link (industrial and agricultural products benefits, labor recovery benefit, dilution purification benefits, cooling and humidifying benefits, carbon fixation and oxygen release benefits, spraying and dust removal benefits), ② benefits of the drainage link (water quality pollution loss caused by direct discharge and sewage treatment loss), and ③ benefits of the pollutant carrying link (pollutant reduction costs). Then, the emergy quantitative analysis technology can realize the unified quantitative evaluation of the economic, social, and ecological benefits in the process of the joint optimal allocation process of regional water consumption and pollutant-carrying capacity.
(3) In the process of establishing the joint optimal allocation model of regional water consumption and pollutant-carrying capacity, one of the key links is to construct the interaction relationship between water consumption and pollutants into the river, according to the physical contact between sewage discharge and pollutant discharge into the river and the water consumption allocated to each department in order to realize the coupling between water quantity and water quality. Then, taking the maximum emergy net benefit as the objective function cannot only realize the unified quantification of the economic, social, and ecological benefits of different sectors in the process of joint optimal allocation, but can also change the traditional multi-objective model of water resource optimal allocation into a single-objective model to avoid the difficulty of determining the multi-objective weight.
(4) The example verification of Zhengzhou city shows that, with the continuous improvement of various control measures and inputs, the contradiction between water supply and demand in Zhengzhou has been alleviated. The constraint of the pollutant-carrying capacity distribution in the water function areas was satisfied to the maximum extent, and the maximum joint optimal allocation benefit was achieved. This shows that the pollutant carrying capacity limitation restricts water supply from the source reversely, and there is a clear mutual influence and restriction relationship between the water distribution and pollutant carrying capacity distribution, which verifies the joint optimal allocation model of regional water consumption and pollutant carrying capacity. Its solution method in this paper is effective and feasible.

Author Contributions

This paper represents a result of collaborative teamwork. Z.W. developed the concept and supervised the work; X.G. conceived of the study, performed model validation, and drafted the manuscript; C.L. and M.G. participated in the design of the study and coordination and revised the manuscript; C.G. and Y.L. collected the data and revised the manuscript; and X.W. provided constructive suggestions and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The National Natural Science Foundation of China (No. NSCF-51909240 and No. NSCF-51739009) funded this research.

Acknowledgments

The authors are grateful to colleagues and friends who shared their meteorological and hydrological data with us. We also thank two reviewers for insightful comments that improved an earlier version of this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, Y. Basic function of the reasonable allocation of water resources management in the implementation of the most rigorous system. China Water Resour. 2010, 20, 26–28. (In Chinese) [Google Scholar]
  2. Wang, H.; You, J. Advancements and development course of research on water resources deployment. SHUILIXUEBAO 2008, 10, 1168–1175. (In Chinese) [Google Scholar]
  3. Jiang, Z.; Yang, J.; Ren, X. Optimal allocation of water resources based on ‘three red lines’ and se-d EA model. Water Sav. Irrig 2016, 11, 81–84. (In Chinese) [Google Scholar]
  4. Yan, D.; Luo, X.; Wang, H.; Qin, D.; She, J. The Dual Gross Control of Rivers Based on the Rational Water Resources Allocation—A Case Study on Tangshan. J. Nat. Resour. 2007, 22, 3–10. (In Chinese) [Google Scholar]
  5. Liu, B.; Chen, X.; Jiang, T. Water resources allocation model based on the dual-control of water quantity and quality. Adv. Water Sci. 2009, 20, 55–59. (In Chinese) [Google Scholar]
  6. You, J.; Xue, X.; Nu, C. Concept and study progress of joint regulation on water quantity and quality. Water Resour. Hydropower Eng. 2010, 41, 7–9, 18. (In Chinese) [Google Scholar]
  7. Fan, Q.; Wang, L.; Li, X.; Qin, J. Discussion on water demand management method based on pollutant discharge limits. Water Sav. Irrig. 2012, 1, 37–38, 43. (In Chinese) [Google Scholar]
  8. Zhang, W. Study on the overall optimal allocation and countermeasures of regional water quantity and quality. Ph.D. Thesis, China University of Mining and Technology, Beijing, China, 2016. (In Chinese). [Google Scholar]
  9. Wu, Z.; Lv, C.; Hu, C.; Guan, X. Research on the Analysis Method and Application of Ecological Economic Value Based on Emergy Theory; Science Press: Beijing, China, 2013. (In Chinese) [Google Scholar]
  10. Wu, Z.; Di, D.; Lv, C.; Guo, X.; Wang, H. Defining and evaluating the social value of regional water re-sources in terms of emergy. Water Policy 2019, 21, 73–90. [Google Scholar] [CrossRef]
  11. Dong, M. Study on Integrated Optimal Macro-Distribution Model of Regional Water Consumption and Pollutant Carrying Capacity. Master’s Thesis, Zhengzhou University, Zhengzhou, China, 2014. (In Chinese). [Google Scholar]
  12. Lan, S.; Qin, P.; Lu, H. Emergy in Ecosystems; Chemical industry press: Beijing, China, 2002. (In Chinese) [Google Scholar]
  13. Guo, X.; Wu, Z.; Wang, H. Emergy evaluation of positive and negative effects of water utilization in regional agricultural production system. Adv. Sci. Technol. Water Resour. 2016, 36, 1–5. (In Chinese) [Google Scholar]
  14. Tian, G. Study on Emergy Evaluation Method of the Value of Ecological Water Use Based on Material Cycle. Master’s Thesis, Zhengzhou University, Zhengzhou, China, 2016. (In Chinese). [Google Scholar]
  15. Li, C.; Yuan, Z.; Ren, H.; Fu, Z.; Cai, Z. Effects of water spray control on fugitive dusts originated from raw material piles in steel plants. China Power Sci. Technol. 2009, 15, 78–83. (In Chinese) [Google Scholar]
  16. Zhao, P.; Feng, Y.; Jing, J.; Han, B.; Bi, X.; Zhu, T.; Zhang, X. Characteristics and control indicators of fugitive dust from building construction sites. Acta Sci. Circumstantiate 2009, 29, 1618–1623. (In Chinese) [Google Scholar]
  17. Zhao, G. Ecological Economics; Chemical industry press: Beijing, China, 2008. (In Chinese) [Google Scholar]
  18. Wu, Z.; Guo, X.; Lv, C.; Wang, H.; Di, D. Study on the quantification method of water pollution ecological compensation standard based on emergy theory. Ecol. Indic. 2018, 92, 189–194. [Google Scholar] [CrossRef]
  19. Lv, C.; Wu, Z. Emergy analysis of regional water ecological—Economic system. Ecol. Eng. 2009, 35, 703–710. [Google Scholar] [CrossRef]
  20. Odum, H.T.; Odum, E.C. Energy Analysis Overview of Nations; Working Paper WP-83-82; International Institute for Applied Systems Analysis: Laxenburg, Austria, 1983; p. 421. [Google Scholar]
  21. Wu, Z.; Guo, X.; Di, D.; Lv, C.; Wang, H. Water efficiency evaluation of a regional water scheme—Zhengzhou, China, using a water ecological–economic system (WEES) and based on emergy theory. Water Sci. Technol. Water Supply 2016, 17, 674–687. [Google Scholar] [CrossRef]
  22. Guo, X. The Integrated Distribution Study of Regional Water Consumption and Pollutant Carrying Capacity Based on Emergy Theory. Ph.D. Thesis, Zhengzhou University, Zhengzhou, China, 2018. (In Chinese). [Google Scholar]
  23. Mu, L. Study on the Water Resources Allocation Based on Secondary Distribution in Zhengzhou Metropolitan Area. Master’s Thesis, Zhengzhou University, Zhengzhou, China, 2017. (In Chinese). [Google Scholar]
Figure 1. Water shortage change trend of each scheme and water use sector under the 50% guarantee rate of Zhengzhou in 2030.
Figure 1. Water shortage change trend of each scheme and water use sector under the 50% guarantee rate of Zhengzhou in 2030.
Water 12 01101 g001
Figure 2. Water structure diagram of each scheme and water use sector under the 50% guarantee rate of Zhengzhou in 2030.
Figure 2. Water structure diagram of each scheme and water use sector under the 50% guarantee rate of Zhengzhou in 2030.
Water 12 01101 g002
Figure 3. Changing trends of pollutants into the river by industrial water use sectors of Zhengzhou in 2030.
Figure 3. Changing trends of pollutants into the river by industrial water use sectors of Zhengzhou in 2030.
Water 12 01101 g003
Figure 4. Changing trends of pollutants into the river from agricultural water use sectors of Zhengzhou in 2030.
Figure 4. Changing trends of pollutants into the river from agricultural water use sectors of Zhengzhou in 2030.
Water 12 01101 g004
Figure 5. Changing trends of pollutants into the river from domestic water use sectors of Zhengzhou in 2030.
Figure 5. Changing trends of pollutants into the river from domestic water use sectors of Zhengzhou in 2030.
Water 12 01101 g005
Figure 6. Pollutant capacity distribution structure diagram of each scheme and water use sector of Zhengzhou in 2030.
Figure 6. Pollutant capacity distribution structure diagram of each scheme and water use sector of Zhengzhou in 2030.
Water 12 01101 g006
Figure 7. Positive benefits comparison of each scheme under the 50% guarantee rate of Zhengzhou in 2030.
Figure 7. Positive benefits comparison of each scheme under the 50% guarantee rate of Zhengzhou in 2030.
Water 12 01101 g007
Figure 8. Negative benefits comparison of each scheme under the 50% guarantee rate of Zhengzhou in 2030.
Figure 8. Negative benefits comparison of each scheme under the 50% guarantee rate of Zhengzhou in 2030.
Water 12 01101 g008
Table 1. The benefits composition of joint optimal allocation of water consumption and pollutant carrying capacity.
Table 1. The benefits composition of joint optimal allocation of water consumption and pollutant carrying capacity.
LinksWater SectorBenefitsEmergy Quantification Method
Water supply linkIndustrial, agricultural, domestic, and ecological water supplyEngineering water supply benefits (+) EM s = ( τ gw τ mw ) × W s
Water use and the consumption linkIndustrial water use and the consumptionIndustrial product benefits (+) EM pi = б i   ×   EM oi
Agricultural water use and the consumptionAgricultural product benefit (+) EM pa = б a   ×   EM oa
Domestic water use and the consumptionLabor recovery benefit (+) EM lr = б d   × EM od × E
Ecological water use and the consumptionDilution of the purification benefits (+) EM dc = W r   ×   τ r
Cooling and humidifying benefits (+) EM ch = W g   ×   L   ×   τ s
Carbon fixation and oxygen release benefits (+) EM co = ( G co 2 × τ co 2   G o 2 × τ o 2 )   ×   W g W p
Spraying and dust removal benefits (+) EM sd = A   ×   h   ×   Δ PM 10   ×   τ n   ×   48.62 %
Drainage linkIndustrial drainageloss of industrial wastewater treatment (-) EM is = W is   ×   UTC   ×   EDR
Agricultural drainagewater quality pollution loss caused by the direct discharge of agricultural drainage (-) EM ad = ( τ wa     τ wb )   ×   W ad
Domestic drainagesewage treatment loss of urban domestic (-) EM ds = W ds   ×   UTC   ×   EDR
water quality pollution loss of rural domestic (-) EM dd = ( τ wa   τ wb )   ×   W dd
Pollutant carrying linkIndustrial water pollutant-carryingPollutant reduction costs in the industrial water sector (-) EM ip = Q p   ×   UCR p   ×   EDR .
Urban domestic water pollutant carryingPollutant reduction costs in the urban domestic water sector (-) EM dp = Q p   ×   UCR p   ×   EDR
Table 2. The water function zoning table in Zhengzhou.
Table 2. The water function zoning table in Zhengzhou.
Name of Water Function AreaSub-AreaRiverSource of Pollutant
The sewage control area of Zhengzhou in Jialu RiverZhengzhou districtJialu RiverIndustry, Agriculture
The sewage control area of Zhengzhou in Qili RiverQili RiverUrban domestic
The sewage control area of Zhongmou in Jialu RiverZhongmou countyJialu RiverIndustry, Urban domestic
The agricultural water area of Zhongmou in Jialu RiverJialu RiverRural domestic, Agriculture
The agricultural water area of Gongyi in Yiluo RiverGongyi cityYiluo RiverRural domestic, Agriculture
The sewage control area of Gongyi in Housi RiverHousi RiverIndustry, Domestic
The sewage control area of Gongyi in Sishui RiverSishui RiverIndustry, Agriculture
The sewage control area of Xingyang in Suoxu RiverXingyang citySuoxu RiverIndustry, Urban domestic
The sewage control area of Xinmi in Shuangji RiverXinmi cityShuangji RiverIndustry, Agriculture, Domestic
The sewage control area of Xinzheng in Shuangji RiverXinzheng cityShuangji RiverIndustry, Agriculture, Domestic
The industrial water area of Dengfeng in Ying RiverDengfeng cityYing RiverIndustry
The sewage control area of Dengfeng in Ying riverYing RiverDomestic, Agriculture
Table 3. Scene settings of the joint optimal allocation in Zhengzhou city.
Table 3. Scene settings of the joint optimal allocation in Zhengzhou city.
Scene SettingsPollutant Carrying Capacity ConstraintsWater Demand SideSewage Treatment SidePollutant Reduction Funds Side
Moderate Water SavingStrengthening Water SavingModerate TreatmentHigh TreatmentCurrent LevelIncreasing Investment
Scene 1
Scene 2
Scene 3
Scene 4
Scene 5
Scene 6
Scene 7
Scene 8
Table 4. The water demand prediction results under a moderate water-saving scene. (Unit: 104 m3).
Table 4. The water demand prediction results under a moderate water-saving scene. (Unit: 104 m3).
Sub-AreaIndustrial Water DemandAgricultural Water DemandDomestic Water DemandEcological Water DemandTotal
P = 50%P = 75%P = 95%UrbanRuralP = 50%P = 75%P = 95%P = 50%P = 75%P = 95%
Gongyi city16,5313389370740293843158674818425,42325,74826,073
Dengfeng city74861894197320322267130510110611013,05313,13713,200
Xingyang city870112,05912,84113,1372223159210010711024,67525,46425,763
Xinmi city90222742286330473267157086939616,68716,81517,002
Zhengzhou district14,30013,54513,96814,35442,870417331,27732,21533,826106,165107,526109,523
Xinzheng city8293377040744254273278032133235815,89616,21116,417
Zhongmou county434731,99534,98435,8922650167743446948841,10344,12745,054
The whole city69,30169,39474,41076,74559,85212,68332,39333,40335,072243,623249,649253,653
Table 5. The water demand prediction results under an intensified water-saving scene. (Unit: 104 m3).
Table 5. The water demand prediction results under an intensified water-saving scene. (Unit: 104 m3).
Sub-AreaIndustrial Water DemandAgricultural Water DemandDomestic Water DemandEcological Water DemandTotal
P = 50%P = 75%P = 95%UrbanRuralP = 50%P = 75%P = 95%P = 50%P = 75%P = 95%
Gongyi city13,6142520280929643733158674818421,52721,82321,981
Dengfeng city68051564162416832236130510110611012,01112,07612,139
Xingyang city683710,79911,32911,7772155159210010711021,48322,02022,471
Xinmi city65162247236825523075157086939613,49413,62213,809
Zhengzhou district12,12410,59310,82211,05340,726417331,27732,21533,82698,893100,060101,902
Xinzheng city6984325433433438268978032133235814,02814,12814,249
Zhongmou county263925,38727,29028,6512618167743446948832,75534,69336,073
The whole city55,51956,36459,58562,11857,232758932,39333,40335,072209,097213,328217,530
Table 6. The available water supply prediction results of Zhengzhou in 2030. (Unit: 104 m3).
Table 6. The available water supply prediction results of Zhengzhou in 2030. (Unit: 104 m3).
Sub-AreaSpecial Water SourceCommon Water Source
Surface WaterGroundwaterUnconventional Water SourceTotalYellow River Water DiversionSouth-to-North Water Diversion
50%75%95%Rainwater CollectionReclaimed Water50%75%95%
50%75%95%
Gongyi city10,030984895507204292271245699024,51624,31323,989
Dengfeng city3823320922695801252238221435714,23313,60512,648
Xingyang city12,09411,83211,5065880260250238466122,89522,62322,285
Xinmi city2988264120796759193177158405913,99913,63613,055
Zhengzhou district111,760111,36911,099283339317795918096129,120128,577128,012
Xinzheng city11,47411,21610,71655451059480423421,35821,08920,575
Zhongmou county21,41921,41921,41916,87334272019,51057,83657,82957,822
The whole city173,588171,534168,53156,39420661836155251,907283,955281,671278,38442,00054,000
Table 7. The emergy parameters of industrial water sector in Zhengzhou.
Table 7. The emergy parameters of industrial water sector in Zhengzhou.
ProjectZhengzhou District Zhongmou County Gongyi City Xingyang City Xinmi CityXinzheng CityDengfeng City
The emergy of industrial water (1020 sej)7.492.425.084.864.5946.213.71
Total emergy input of industrial production (1020 sej)268.254.9209.9101.6143.6101.4225.6
Total emergy output of industrial production (1020 sej)685.8104.8558.4247.4372.2160.5610.9
Emergy contribution rate of industrial water (%)2.794.402.424.783.204.561.65
Emergy contribution amount of industrial water (1020 sej)19.14.6113.511.83 11.917.3110.08
Monetary value of industrial water (108 ¥)66.7516.0847.1441.2541.5525.5235.16
Emergy benefit of per m3 industrial water (1012 sej/m3)15.1311.1615.0414.3015.219.3115.35
Monetary benefit of per m3 industrial water (¥/m3)52.8038.9352.4749.8953.0732.4953.56
Table 8. The emergy parameters of agricultural water sector in Zhengzhou.
Table 8. The emergy parameters of agricultural water sector in Zhengzhou.
ProjectZhengzhou District Zhongmou County Gongyi City Xingyang CityXinmi CityXinzheng CityDengfeng City
The emergy of agricultural water (1020 sej)0.082.880.200.560.250.610.16
Total emergy input of agricultural production (1020 sej)10.416.613.59.499.0110.638.50
Total emergy output of agricultural production (1020 sej)26.678.139.8187.162.361.915.5
Emergy contribution rate of agricultural water (%)8.4417.381.505.952.765.711.93
Emergy contribution amount of agricultural water (1020 sej)2.2513.60.6011.131.723.540.30
Monetary value of agricultural water (108 ¥)7.8447.362.0838.836.0012.341.05
Emergy benefit per m3 agricultural water (1012 sej/m3)3.606.213.2126.519.277.921.99
Monetary benefit per m3 agricultural water (¥/m3)12.5521.6511.1892.4732.3627.616.97
Table 9. The emergy parameters of the urban domestic water sector in Zhengzhou.
Table 9. The emergy parameters of the urban domestic water sector in Zhengzhou.
ProjectZhengzhou District Zhongmou County Gongyi CityXingyang CityXinmi CityXinzheng CityDengfeng City
Emergy contribution rate of urban domestic water (%)42.24 49.44 63.72 47.43 46.40 48.26 57.42
Engel coefficient (%)37.40 29.90 33.10 31.70 25.90 28.20 30.50
Disposable income of urban residents (¥/p)22,477 16,588 18,186 18,188 18,194 18,204 17,565
Emergy value of urban domestic water (1014 sej/p)10.18 7.03 11.00 7.84 6.27 7.10 8.82
Monetary value of urban domestic water (¥/p)3550.86 2452.13 3835.67 2734.62 2186.48 2477.44 3076.18
Urban domestic water amount per person per year (m3/p)42.07 40.67 34.20 33.03 35.39 48.04 21.01
Emergy benefit per m3 urban domestic water (1012 sej/m3)24.20 17.29 32.16 23.74 17.72 14.78 41.98
Monetary benefit per m3 urban domestic water (¥/m3)84.40 60.29 112.15 82.79 61.78 51.57 146.42
Table 10. The emergy parameters of the rural domestic water sector in Zhengzhou.
Table 10. The emergy parameters of the rural domestic water sector in Zhengzhou.
ProjectZhengzhou District Zhongmou County Gongyi City Xingyang City Xinmi City Xinzheng CityDengfeng City
Emergy contribution rate of rural domestic water (%)31.83 26.10 28.20 30.07 19.81 21.30 25.88
Engel coefficient (%)37.40 29.90 33.10 31.70 25.90 28.20 30.50
Disposable income of rural residents (¥/p)13,522 10,216 11,392 10,846 10,835 11,344 9788
Emergy value of rural domestic water (1014 sej/p)4.62 2.29 3.05 2.96 1.59 1.95 2.21
Monetary value of rural domestic water (¥/p)1609.81 797.25 1063.39 1033.89 555.92 681.39 772.62
Rural domestic water amount per person per year (m3/p)42.07 40.67 34.20 33.03 35.39 48.04 21.01
Emergy benefit per m3 rural domestic water (1012 sej/ m3)10.98 5.63 8.92 8.96 4.49 4.06 10.52
Monetary benefit per m3 rural domestic water (¥/m3)38.27 19.60 31.09 31.30 15.71 14.18 36.77
Table 11. The emergy parameters of the ecological water sector in Zhengzhou.
Table 11. The emergy parameters of the ecological water sector in Zhengzhou.
ProjectZhengzhou District Zhongmou County Gongyi City Xingyang City Xinmi City Xinzheng CityDengfeng City
Average annual temperature (°C) 13.713.713.713.713.713.713.7
Latent heat of evaporation (J/g) 2475.02475.02475.02475.02475.02475.02475.0
Evaporation deep (mm) 620.00621.00622.00623.00624.00625.00626.00
Annual evaporation (1015 m3) 0.620.870.650.590.610.550.76
Evaporation energy (1019 J) 0.150.220.160.150.150.140.19
Emergy transformity of steam (sej/J) 13.2014.2015.2016.2017.2018.2019.20
Evaporation emergy (1020 sej) 0.200.310.240.240.260.250.36
Ecological water consumption (108 m3) 0.710.821.400.610.690.591.72
Emergy benefit per m3 ecological water (1012 sej/m3) 0.280.380.170.390.370.420.21
Monetary benefit per m3 ecological water (¥/m3)0.10 0.13 0.06 0.14 0.13 0.15 0.07
Table 12. The emergy transformity of different water bodies in Zhengzhou.
Table 12. The emergy transformity of different water bodies in Zhengzhou.
ProjectZhengzhou District Zhongmou County Gongyi CityXingyang CityXinmi CityXinzheng CityDengfeng City
Emergy transformity of surface water (1011 sej/m3)8.108.104.078.798.238.593.91
Emergy transformity of groundwater (1012 sej/m3)1.441.261.201.251.301.291.23
water diversion cost (¥/m3)0.290.290.290.290.290.290.29
Reservoir construction cost (¥/m3)1.141.141.141.141.141.141.14
Tube Wells pumping station construction cost (¥/m3)0.240.240.240.240.240.240.24
Water conservancy project cost of yellow river (¥/m3)0.590.590.590.590.590.590.59
Emergy transformity of surface water after engineering consideration (1012 sej/m3)1.221.220.821.291.231.270.80
Emergy transformity of groundwater after engineering consideration (1012 sej/m3)1.601.421.361.401.451.451.38
Emergy transformity of yellow river after engineering consideration (1011 sej/m3)1.691.691.691.691.691.691.69
Table 13. The pollutant carrying capacity parameters of the main water function area. (Unit: t/a).
Table 13. The pollutant carrying capacity parameters of the main water function area. (Unit: t/a).
Sub-AreaWater Function AreaCODNH3-N
Zhengzhou districtThe sewage control area of Zhengzhou in Jialu River4052.1344.4
The sewage control area of Zhengzhou in Qili River5287.1228.3
Zhongmou countyThe sewage control area of Zhongmou in Jialu River2376.1117.9
The agricultural water area of Zhongmou in Jialu River512.913.1
GongyicityThe agricultural water area of Gongyi in Yiluo River686.339
The sewage control area of Gongyi in Housi River1893.4186.3
Xingyang cityThe sewage control area of Gongyi in Sishui River1431.175.4
The sewage control area of Xingyang in Suoxu River2235.8139.2
Xinmi cityThe sewage control area of Xinmi in Shuangji River2613158.9
Xinzheng cityThe sewage control area of Xinzheng in Shuangji River1406.666.3
Dengfeng cityThe industrial water area of Dengfeng in Ying River625.934.3
The sewage control area of Dengfeng in Ying River1072.4111.8
Total 23,997.21479.4
Table 14. The integrated optimal distribution results of Zhengzhou in 2030 (p = 50%).
Table 14. The integrated optimal distribution results of Zhengzhou in 2030 (p = 50%).
SchemeWater Use SectorQuantity of Water DemandQuantity of Water Distribution Quantity of Water ShortagePollutant Discharge before Pollutant Carrying Link Pollutant Reduction in Pollutant Carrying LinkPollutant in the River after a Pollution Carrying LinkDistribution of Pollution Carrying Capacity
(104 m3)(104 m3)(104 m3)(t) (t) (t) (t)
CODNH3-NCODNH3-NCODNH3-NCODNH3-N
Scheme 1Industrial90,10978,23511,87427,351.35702.210,247.64680.117,103.71022.117,103.71022.1
Urban domestic50,24550,24503901.13051461.6250.32439.554.72439.554.7
Rural domestic809380930803.1181.4300.9148.9502.232.5502.232.5
Agricultural48,48937,90810,5816319.52064.92367.71694.83951.8370.13951.8370.1
Ecological environment31,35731,357000000000
Total228,294205,83722,45638,374.97799.314,377.76319.923,997.21479.423,997.21479.4
Scheme 2Industrial90,10979,47010,63927,792.15558.411,304.64799.916,487.5758.517,109.71022.7
Urban domestic50,24550,24503964.3297.31612.4256.72351.940.62440.554.7
Rural domestic809380930820.3177.7333.7153.4486.624.350532.7
Agricultural48,48938,98395066403.12007.12604.41733.23798.7273.93941.9369.3
Ecological environment31,35731,357000000000
Total228,294208,14820,14638,9808040.615855.16943.223124.91097.423,997.21479.4
Scheme 3Industrial90,10980,04310,06626,7145342.89827.94411.516,886.1931.317,068.71013.3
Urban domestic50,24550,24503610.6270.81328.3223.62282.347.2230751.4
Rural domestic809380930829.9179.8305.4148.5524.531.3530.334.1
Agricultural48,48939,17493156434.52016.92355.71657.34078.8359.64091.2380.7
Ecological environment31,35731,357000000000
Total228,29420,891219,38237,5897810.313,817.26440.823771.81369.523,997.21479.4
Scheme 4Industrial90,10980,195991426,738.85347.810,824.6485115914.2496.817,067.41013.6
Urban domestic50,24550,24503620.4270.81465.6245.72154.825.12310.951.3
Rural domestic809380930833.2179.8337.3163.1495.916.7531.834.1
Agricultural48,48938,85296376403.12007.12592.21820.73810.9186.44087.1380.4
Ecological environment31,35731,357000000000
Total228,294208,74219,55237,5957805.515,219.77080.422,375.3725.123,997.21479.4
Scheme 5Industrial77,73376,120161326,973.45462.210,001.94450.216,971.5101216,971.51012
Urban domestic48,81548,81503889.6295.31442.3240.62447.354.72447.354.7
Rural domestic809380930836.3183.5310.1149.5526.234526.234
Agricultural43,74543,5002456440.320442388.11665.34052.2378.74052.2378.7
Ecological environment31,32531,325000000000
Total209,711207,853185838,139.8785014,142.66370.623,997.21479.423,997.21479.4
Scheme 6Industrial77,73376,657107627,184.45436.911,047.34789.916,137.164716,9811012.9
Urban domestic48,81548,81503920.52941593.2259.12327.334.9244954.8
Rural domestic809380930846.1183.3343.9161.5502.221.8528.534.2
Agricultural43,74543,2804656465.42026.62627.41785.43838241.24038.7377.6
Ecological environment31,32531,325000000000
Total209,711208,170154138,416.47940.815,611.86995.922,804.6944.923,997.21479.4
Scheme 7Industrial77,73377,52820527184.45436.99582.3440716,550.2819.516,950.91003.8
Urban domestic48,81548,81503920.52941301.9224.52248.641.8230351.1
Rural domestic809380930846.1183.3310.6154.8536.528.7549.535.3
Agricultural43,74542,41213336432.82016.42370.71708.94094.7317.74193.8389.2
Ecological environment31,32531,325000 0 00
Total209,711208,173153838,383.87930.613,565.66495.123,429.91207.823,997.21479.4
Scheme 8Industrial77,73377,733026,2265245.210,5834843.915,643401.316,961.11004.6
Urban domestic48,81548,81503555.2266.61434.6246.32120.620.32299.251.1
Rural domestic809380930859186.1346.7171.9512.314.2555.535.6
Agricultural43,74542,212153362741966.62609.11871.63665954181.4388.1
Ecological environment31,32531,325000000000
Total209,711208,178153336,914.27664.514,973.37133.721940.9530.823,997.21479.4
Table 15. Joint optimal allocation benefits comparison of different schemes under the 50% guarantee rate of Zhengzhou in 2030 (unit: 1019 sej).
Table 15. Joint optimal allocation benefits comparison of different schemes under the 50% guarantee rate of Zhengzhou in 2030 (unit: 1019 sej).
BenefitScheme 1Scheme 2Scheme 3Scheme 4Scheme 5Scheme 6Scheme 7Scheme 8
Positive benefits of water use and consumption link2256.832264.092274.12280.492513.792519.682541.952548.16
Negative benefits of drainage link−150.81−152.31−121.66−122.81−140.12−141.63−122−123.29
Negative benefits of pollutant carrying link−17.22−18.66−17.22−18.66−17.22−18.66−17.22−18.66
Joint optimal allocation benefits2088.82093.122135.222139.022356.452359.392402.732406.21
Table 16. Allocation results’ comparison of an optimal scheme of Zhengzhou in 2030.
Table 16. Allocation results’ comparison of an optimal scheme of Zhengzhou in 2030.
Allocation MethodsAllocation Water Quantity (104 m3)Corresponding Pollution Discharge (t)Pollutant Carrying Capacity (t)Satisfaction Degree of Pollution Carrying Capacity
CODNH3-NCODNH3-NCODNH3-N
Emergy method205,77921,940.9 530.8 23,997.21479.4
Traditional method262,79646,888.6 9761.1 ××
Note: √ means satisfied, × means not satisfied.

Share and Cite

MDPI and ACS Style

Guo, X.; Wu, Z.; Wang, X.; Lv, C.; Gu, C.; Li, Y.; Gao, M. The Joint Optimal Allocation Study of Regional Total Water Consumption and Pollutant Carrying Capacity of Water Function Areas Based on Emergy Theory. Water 2020, 12, 1101. https://doi.org/10.3390/w12041101

AMA Style

Guo X, Wu Z, Wang X, Lv C, Gu C, Li Y, Gao M. The Joint Optimal Allocation Study of Regional Total Water Consumption and Pollutant Carrying Capacity of Water Function Areas Based on Emergy Theory. Water. 2020; 12(4):1101. https://doi.org/10.3390/w12041101

Chicago/Turabian Style

Guo, Xi, Zening Wu, Xinling Wang, Cuimei Lv, Changkuan Gu, Yang Li, and Meiyan Gao. 2020. "The Joint Optimal Allocation Study of Regional Total Water Consumption and Pollutant Carrying Capacity of Water Function Areas Based on Emergy Theory" Water 12, no. 4: 1101. https://doi.org/10.3390/w12041101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop