CN110210710B - Water resource bearing capacity quantization method based on load balancing - Google Patents
Water resource bearing capacity quantization method based on load balancing Download PDFInfo
- Publication number
- CN110210710B CN110210710B CN201910370288.4A CN201910370288A CN110210710B CN 110210710 B CN110210710 B CN 110210710B CN 201910370288 A CN201910370288 A CN 201910370288A CN 110210710 B CN110210710 B CN 110210710B
- Authority
- CN
- China
- Prior art keywords
- water
- value
- kth
- partition
- water consumption
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 192
- 238000000034 method Methods 0.000 title claims abstract description 103
- 238000013139 quantization Methods 0.000 title claims description 14
- 238000011156 evaluation Methods 0.000 claims abstract description 17
- 238000011002 quantification Methods 0.000 claims abstract description 9
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 5
- 230000002068 genetic effect Effects 0.000 claims abstract description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 4
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 4
- 239000001301 oxygen Substances 0.000 claims abstract description 4
- 239000000126 substance Substances 0.000 claims abstract description 4
- 238000005192 partition Methods 0.000 claims description 72
- 238000011161 development Methods 0.000 claims description 42
- 230000008569 process Effects 0.000 claims description 37
- 238000004519 manufacturing process Methods 0.000 claims description 25
- 241001465754 Metazoa Species 0.000 claims description 11
- 230000007613 environmental effect Effects 0.000 claims description 10
- 239000008235 industrial water Substances 0.000 claims description 9
- 239000010865 sewage Substances 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000009776 industrial production Methods 0.000 claims description 6
- 239000002352 surface water Substances 0.000 claims description 6
- 230000008685 targeting Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 239000003673 groundwater Substances 0.000 claims description 3
- 244000144972 livestock Species 0.000 claims 1
- 238000012545 processing Methods 0.000 claims 1
- 238000005728 strengthening Methods 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a load balancing-based water resource bearing capacity quantification method, which belongs to the technical field of water resource evaluation and management, and establishes a connection among the available water resource amount and water consumption, the COD (chemical oxygen demand) pollution discharge amount and the total COD limit discharge amount and the load balancing of a socioeconomic subsystem and an ecological environment subsystem, so that a load balancing-based water resource bearing capacity quantification model is constructed, a genetic algorithm is adopted to convert multiple targets into single targets in a proper mode to solve the model, and a water resource bearing capacity quantification result is obtained.
Description
Technical Field
The invention belongs to the technical field of water resource evaluation and management, and particularly relates to a water resource bearing capacity quantization method based on load balancing.
Background
The research on the water resource bearing capacity goes through three main stages of conceptual proposal (around 1990), theoretical discussion (1990 to 2000) and method application (2000 to date).
To date, numerous mathematical methods have been applied to the research of water resource bearing capacity, but the requirement of the coordinated development of socioeconomic and ecological environment is not considered in the quantification process. In the field of water resources, the load refers to the pressure caused by the load-bearing object, namely the socioeconomic and ecological environment. On the premise of maintaining sustainable development of social economy and ecological systems, the bearing main body, namely the maximum load borne by the water resource system, is the water resource bearing capacity, and when the actual load of the water resource exceeds the water resource bearing capacity, the water resource system enters an overload state.
In the process of water resource development and utilization, direct contradiction exists between social economy and ecological environment water distribution, so that the development relationship between two systems is coordinated, and the water resource bearing capacity in an equilibrium state is quantized.
Disclosure of Invention
The invention aims to: the invention aims to provide a water resource bearing capacity quantization method based on load balancing, which can embody a sustainable development idea and an ecological idea in the quantization process so as to maximize the overall benefit of a system, so that the obtained water resource bearing capacity quantization result in an equilibrium state is more reasonable and reliable.
The technical scheme is as follows: in order to achieve the above purpose, the present invention provides the following technical solutions:
a water resource bearing capacity quantization method based on load balancing comprises the following steps:
1) Based on a load balancing theory, a water resource system, a social and economic system and an ecological environment system are built to be correlated and restricted with each other to form a water resource bearing power model, and the objective function of the model is respectively determined from three aspects of social scale, economic development and environmental protection;
2) Constraint conditions of the water resource bearing capacity quantization model comprise four aspects of water resource quantity constraint, socioeconomic constraint, water environment constraint and load balancing constraint;
3) After the multi-objective problem of the water resource bearing capacity quantification model is converted into a single-objective problem, a genetic algorithm is adopted to solve the problem, so that the population scale and the socioeconomic development scale which can be borne by the water resource system under different conditions are obtained.
Further, in step 1), the objective function includes:
2.1 Social goal: targeting regional loadable population (POP) max:
in the method, in the process of the invention,
population for the kth partition;
2.2 Economic goal: targeting regional domestic total production (GDP) max:
in the method, in the process of the invention,
a domestic production total value for the kth partition;
2.3 Environmental objective): the minimum COD discharge of the regional chemical oxygen demand is taken as the target:
in the method, in the process of the invention,
COD discharge amount of the kth zone.
Further, in step 2), the water resource amount constraint includes two aspects that the water consumption amount in each partition is not more than the available water resource amount of the unit, and the sum of the water consumption amounts of each user in each partition is not more than the available water resource amount in the area:
the water consumption in each partition does not exceed the available amount constraint of the unit water resource:
T ZHSK ≤T KLYK ；
wherein T is ZHSK For each partition water consumption, T KLYK The amount of water resource available for the unit.
The sum of the water consumption of each water user in each partition is not more than the water resource availability constraint of the area:
T ZHS ≤T KLY ；
wherein T is KLY The total available water resource amount; t (T) ZHS Is the total water consumption.
The method for calculating the total available water resource comprises the following steps:
T KLY ＝T DBS +T DXS -Δ RE ；
wherein T is DBS T is the available amount of surface water resources DXS Delta as the available amount of groundwater resources RE Repeating the amount of the surface water resource and the underground water resource;
the total water consumption calculating method comprises the following steps:
in the method, in the process of the invention,
rural domestic water consumption for kth zone, </i >>
Urban domestic water consumption for the kth partition, < >>
For the kth zone-water consumption, -/-, water production>
For the k-th partition, water consumption, < >>
For the k-th zone, water consumption, three-way production,/->
The ecological water consumption outside the river channel of the kth subarea; wherein->
Water consumption for the k-th partition planting industry,/->
For the k-th zone fishery water consumption, < >>
Water is consumed for the kth zone animal husbandry.
Further, in step 2), the socioeconomic constraints: including two aspects of industrial structure constraints and people average GDP constraints:
the industrial structure constraint is that the proportion of the industrial production value increment value of each partition to the total GDP value of the area cannot exceed the upper limit and the lower limit:
in the method, in the process of the invention,
is regional GDP total value,/->
Increasing the value of the ith industry production value for the kth partition,/->
Increasing the lower limit of the ratio of the value to the value of the i-th industrial production value for the kth partition,/->
Increasing the upper limit of the value proportion for the ith industry production value of the kth partition;
the regional GDP total value calculation method comprises the following steps:
in the method, in the process of the invention,
water consumption for ten thousand yuan per yield for the kth zone, < >>
Water consumption for ten thousand yuan production value of k-th partition>
The water consumption is ten thousand yuan for the third production of the kth partition;
the people average GDP is restricted, and the people average GDP of each partition is larger than or equal to the lower limit value of the people average GDP:
in the method, in the process of the invention,
for the total population of the kth partition, A perGDP Is the lower limit value of the average person GDP;
the total population calculation method of the kth partition comprises the following steps:
in the method, in the process of the invention,
for the k-th partition rural population, +.>
For the kth computational unit town population,
for the k-th subarea rural people daily water consumption, < >>
And daily water consumption is daily used for people in towns in the kth subarea.
Further, in step 2), the water environment is constrained: the total COD discharged by all pollution sources is smaller than the total COD limit discharge:
in the method, in the process of the invention,
point source emissions for the kth zone, +.>
For the k-th subarea non-point source pollution emission quantity, RECOD k The total amount is limited for the kth partition.
The k-th partition point source emission amount calculating method comprises the following steps:
wherein:
urban domestic water consumption rate for kth computing unit,/->
Third industrial water consumption rate for kth computing unit,/->
For the kth computing unit, the second industrial water consumption rate, A CLXS Is a sewage treatment coefficient A SHCOD Is COD concentration of town domestic water after sewage treatment, A GYCOD Is the COD concentration of the second industrial water after sewage treatment.
The k-th partition area source pollution emission amount calculating method comprises the following steps:
wherein:
the average water consumption per mu of the unit planting industry and the average water consumption per mu of the fishery industry are calculated for the kth unit planting industry. A is that ZZYCOD 、A YYCOD 、A XMYCOD 、A NCCOD COD generation coefficients of planting industry, fishery industry, animal husbandry and rural life are respectively A ZZYin 、A YYin 、A XMYin 、A NCSHin The river coefficients of the COD of the planting industry, the fishery industry and the rural life are respectively.
Further, in step 2), the load balancing constraint: the coordination development degree of the social and economic systems and the ecological environment system meets the requirement to ensure the stable development of the whole system.
Wherein: c is the degree of coordination, and the value is between 0 and 1; n is an adjustment coefficient, taken as 2 according to the number of the system of the invention; i 1 The method is a coordinated development evaluation index of a social and economic subsystem; i 2 The coordinated development evaluation index of the ecological environment subsystem; c (C) min Is the lowest value of the required coordination level.
Further, in the step 2), the load balancing constraint selects the average-human GDP as a coordinated development evaluation index of the social and economic subsystem of the region, and the average-human GDP is a forward index; and selecting the COD discharge amount of the ten thousand-yuan GDP as a coordinated development evaluation index of an ecological environment subsystem in the region, and taking the COD discharge amount as a negative index.
The method adopts a threshold method to carry out dimensionality removal treatment on the indexes, and respectively reflects the indexes of social economic development and ecological environment conditions, and different formulas are required to be adopted for index standardization; the formula is as follows:
positive indicators, i.e. when the greater the indicator value the better the efficacy,
I i ＝C i /S i ；
negative indicators, i.e. when the efficacy is better the smaller the indicator value,
I i ＝S i /C i ；
wherein I is i Is the index value of index i, S i Is the standard value of index i, C i Is the actual value of index i.
Further, the coordination in the load balancing constraint can be divided into 5 coordination levels of extremely uncoordinated, barely coordinated, basic coordinated and coordinated according to the magnitude of the value, wherein the coordination level in the extremely uncoordinated level is more than or equal to 0 and less than or equal to 0.49; c is more than or equal to 0.5 and less than or equal to 0.74 in the uncoordinated level; the coordination degree in the marginal coordination grade is more than or equal to 0.75 and less than or equal to 0.84; the coordination degree C in the basic coordination grade is more than or equal to 0.85 and less than or equal to 0.94; the coordination degree C in the coordination grade is more than or equal to 0.95 and less than or equal to 1.
Further, in the step 3), the specific process is as follows:
to achieve the optimal level f for the nth target nbest The corresponding standard value is set to 1, the worst level f nworst And setting the corresponding standard value to be 0.01, and calculating the corresponding standard value by each target according to an interpolation formula:
for social and economic goals, solve the maximization problem:
f' n (x)＝f n (x)/f nbest ；
for an environmental objective solution objective, solving a minimization problem:
f' n (x)＝1-f n (x)/f nbest ；
converting each target into a standard value between 0 and 1, and then converting the multiple targets into a single-target construction fitness function in a weighted summation mode:
in the method, in the process of the invention,
ρ n and determining the target weight coefficient by adopting an expert consultation method, an analytic hierarchy process and the like.
The beneficial effects are that: compared with the prior art, the water resource bearing capacity quantification method based on the load balancing increases the load balancing constraint when a model is built, introduces the load balancing into the water resource research field, considers the coordination and balanced development of a social and economic system and an ecological environment system, can embody the sustainable development concept and the ecological concept in the quantification process, is consistent with the primary principle based on the water resource bearing capacity research institute, namely requires the coordination development of economy, society, resource and environmental protection so as to maintain the load balancing, and has stronger practicability and wide applicability.
Detailed Description
The invention will be further described with reference to specific examples.
According to the water resource bearing capacity quantization method based on load balancing, provided by the invention, the maximum bearing population, the maximum GDP and the minimum COD discharge are selected as objective functions, model constraint conditions are established from four aspects of water resource quantity, social economy, water environment and load balancing, and a multi-target water resource bearing capacity quantization model based on load balancing is established by utilizing the relation of the water consumption of each water department and the available water resource quantity of a watershed, the pollution discharge quantity and the total limit discharge quantity of each water department and the load balancing constraint of a social economy subsystem and an ecological environment subsystem, and solving the model by adopting a genetic algorithm.
The load balancing constraint selects the average-human GDP as a coordinated development evaluation index of the social and economic subsystem of the region, and the average-human GDP is a forward index; and selecting the COD discharge amount of the ten thousand-yuan GDP as a coordinated development evaluation index of an ecological environment subsystem in the region, and taking the COD discharge amount as a negative index.
A water resource bearing capacity quantization method based on load balancing comprises the following steps:
1) Based on a load balancing theory, a water resource system, a social and economic system and an ecological environment system are built to be correlated and restricted with each other to form a water resource bearing power model, and the objective function of the model is respectively determined from three aspects of social scale, economic development and environmental protection;
2) Constraint conditions of the water resource bearing capacity quantization model comprise four aspects of water resource quantity constraint, socioeconomic constraint, water environment constraint and load balancing constraint;
3) After the multi-objective problem of the water resource bearing capacity quantification model is converted into a single-objective problem, a genetic algorithm is adopted to solve the problem, so that the population scale and the socioeconomic development scale which can be borne by the water resource system under different conditions are obtained.
Step 1) determining an objective function from three aspects, respectively, comprising:
2.1 Social goal: targeting regional loadable population (POP) max:
in the method, in the process of the invention,
population for the kth partition;
2.2 Economic goal: targeting regional domestic total production (GDP) max:
in the method, in the process of the invention,
a domestic production total value for the kth partition;
2.3 Environmental objective): the minimum COD discharge of the regional chemical oxygen demand is taken as the target:
in the method, in the process of the invention,
COD discharge amount of the kth zone.
In the step 2), the constraint of the water resource quantity comprises two aspects that the water consumption in each partition is not more than the available water resource quantity of the unit, and the sum of the water consumption of each water user in each partition is not more than the available water resource quantity in the area:
the water consumption in each partition does not exceed the water resource availability constraint of the unit:
T ZHSK ≤T KLYK ；
wherein T is ZHSK For each partition water consumption, T KLYK The amount of water resource available for the unit.
The sum of the water consumption of each water user in each partition is not more than the water resource availability constraint of the area:
T ZHS ≤T KLY ；
wherein T is KLY The total available water resource amount; t (T) ZHS Is the total water consumption.
The method for calculating the total available water resource comprises the following steps:
T KLY ＝T DBS +T DXS -Δ RE ；
wherein T is DBS T is the available amount of surface water resources DXS Delta as the available amount of groundwater resources RE Repeating the amount of the surface water resource and the underground water resource;
the total water consumption calculating method comprises the following steps:
in the method, in the process of the invention,
rural domestic water consumption for kth zone, </i >>
Urban domestic water consumption for the kth partition, < >>
For the kth zone-water consumption, -/-, water production>
For the k-th partition, water consumption, < >>
For the k-th zone, water consumption, three-way production,/->
The ecological water consumption outside the river channel of the kth subarea; wherein->
Water consumption for the k-th partition planting industry,/->
For the k-th zone fishery water consumption, < >>
Water is consumed for the kth zone animal husbandry.
In step 2), socioeconomic constraints: including two aspects of industrial structure constraints and people average GDP constraints:
the constraint of the industrial structure is that the proportion of the increased value of each industrial production value of each partition to the total value of GDP in the area cannot exceed the upper limit and the lower limit of the GDP:
in the method, in the process of the invention,
is regional GDP total value,/->
Increasing the value of the ith industry production value for the kth partition,/->
Increasing the lower limit of the ratio of the value to the value of the i-th industrial production value for the kth partition,/->
Increasing the upper limit of the value proportion for the ith industry production value of the kth partition;
the regional GDP total value calculating method comprises the following steps:
in the method, in the process of the invention,
water consumption for ten thousand yuan per yield for the kth zone, < >>
Water consumption for ten thousand yuan production value of k-th partition>
The water consumption is ten thousand yuan for the third production of the kth partition;
and the people average GDP is restricted, and the people average GDP of each partition is larger than or equal to the lower limit value of the people average GDP:
in the method, in the process of the invention,
for the total population of the kth partition, A perGDP Under GDP for people's averageA limit value;
the total population calculation method of the kth partition is as follows:
in the method, in the process of the invention,
for the k-th partition rural population, +.>
For the kth computational unit town population,
for the k-th subarea rural people daily water consumption, < >>
And daily water consumption is daily used for people in towns in the kth subarea.
In step 2), water environment constraint: the total COD discharged by all pollution sources is smaller than the total COD limit discharge:
in the method, in the process of the invention,
point source emissions for the kth zone, +.>
For the k-th subarea non-point source pollution emission quantity, RECOD k The total amount is limited for the kth partition.
The k-th partition point source emission amount calculating method comprises the following steps:
wherein:
urban domestic water consumption rate for kth computing unit,/->
Third industrial water consumption rate for kth computing unit,/->
For the kth computing unit, the second industrial water consumption rate, A CLXS Is a sewage treatment coefficient A SHCOD Is COD concentration of town domestic water after sewage treatment, A GYCOD Is the COD concentration of the second industrial water after sewage treatment.
The k-th partition area source pollution emission amount calculating method comprises the following steps:
wherein:
the average water consumption per mu of the unit planting industry and the average water consumption per mu of the fishery industry are calculated for the kth unit planting industry. A is that ZZYCOD 、A YYCOD 、A XMYCOD 、A NCCOD COD generation coefficients of planting industry, fishery industry, animal husbandry and rural life are respectively A ZZYin 、A YYin 、A XMYin 、A NCSHin The river coefficients of the COD of the planting industry, the fishery industry and the rural life are respectively.
In step 2), the load balancing constraint: the coordination development degree of the social and economic systems and the ecological environment system meets the requirement to ensure the stable development of the whole system.
Wherein: c is the degree of coordination, and the value is between 0 and 1; n is an adjustment coefficient, taken as 2 according to the number of the system of the invention; i 1 The method is a coordinated development evaluation index of a social and economic subsystem; i 2 The coordinated development evaluation index of the ecological environment subsystem; c (C) min Is the lowest value of the required coordination level.
In the step 2), the average-human GDP is selected from the load balancing constraint as a coordinated development evaluation index of the social and economic subsystem in the region, and the coordinated development evaluation index is a forward index; and selecting the COD discharge amount of the ten thousand-yuan GDP as a coordinated development evaluation index of an ecological environment subsystem in the region, and taking the COD discharge amount as a negative index.
The method adopts a threshold method to carry out dimensionality removal treatment on the indexes, and respectively reflects the indexes of social economic development and ecological environment conditions, and different formulas are required to be adopted for index standardization; the formula is as follows:
positive indicators, i.e. when the greater the indicator value the better the efficacy,
I i ＝C i /S i ；
negative indicators, i.e. when the efficacy is better the smaller the indicator value,
I i ＝S i /C i ；
wherein I is i Is the index value of index i, S i Is the standard value of index i, C i Is the actual value of index i.
The coordination in the load balancing constraint can be divided into 5 coordination levels of extremely uncoordinated, barely coordinated, basic coordinated and coordinated according to the magnitude of the value, wherein the coordination level in the extremely uncoordinated level is more than or equal to 0 and less than or equal to 0.49; c is more than or equal to 0.5 and less than or equal to 0.74 in the uncoordinated level; the coordination degree in the marginal coordination grade is more than or equal to 0.75 and less than or equal to 0.84; the coordination degree C in the basic coordination grade is more than or equal to 0.85 and less than or equal to 0.94; the coordination degree C in the coordination grade is more than or equal to 0.95 and less than or equal to 1.
In step 3), the specific process is as follows:
to achieve the optimal level f for the nth target nbest The corresponding standard value is set to 1, the worst level f nworst The corresponding standard value is set to 0.01, thenStandard values corresponding to the interpolation formula are calculated:
for social and economic goals, solve the maximization problem:
f' n (x)＝f n (x)/f nbest ；
for an environmental objective solution objective, solving a minimization problem:
f' n (x)＝1-f n (x)/f nbest ；
converting each target into a standard value between 0 and 1, and then converting the multiple targets into a single-target construction fitness function in a weighted summation mode:
in the method, in the process of the invention,
ρ n the target weight coefficient can be determined by adopting an expert consultation method, a hierarchical analysis method and the like.
Examples
Step 1, assume a region with 5 sub-regions, each of which is partition 1-5. In this embodiment, the current 2015 year, the recent planning year 2020 and the long-term planning year 2030 are selected to reflect different historical periods and characteristics of bearing capacity of water resources under socioeconomic development, and meanwhile, the influences of the available water resources on bearing capacity of water resources of each planning year are considered under different incoming water frequencies of 50% and 75%, so that 6 scenario schemes are set in total, as shown in table 1.
Table 1 scenario set table
And 2, calculating and determining constraint condition parameters such as regional water resource availability, daily life water consumption of towns and villages in each partition, water consumption of increment value of each industry ten thousand yuan, upper and lower limits of development proportion, water consumption rate of each water department, unit scale water consumption of each horizontal year planting industry, fishery and animal husbandry, relevant parameters of point and surface source pollution, total COD limit discharge of each partition and the like, and the constraint condition parameters are shown in tables 2-13.
TABLE 2 total amount of Water resource availability
TABLE 3 daily Water consumption for urban and rural people
Table 4 increases the water consumption by each of the industrial ten thousand yuan
TABLE 5 upper limit of the industrial development ratio
TABLE 6 lower limit of the development ratio of industries
TABLE 7 Water consumption rates of Water departments
Table 8 Water consumption per unit of animal husbandry in 2015
Table 9 Water consumption per unit of the plant, fishery and animal husbandry in 2020
Table 10 2030 represents the water consumption per unit scale of the plantation, fishery and animal husbandry
TABLE 11 Point Source pollution related parameters
TABLE 12 parameters related to surface source pollution
Table 13 COD limit discharge Total amount of each partition
And 3, determining parameters of the load balancing constraint as follows:
the standard value of the average human GDP in 2015 is determined as an average human GDP value 49680 yuan/person corresponding to the basic modernization level of social and economic development, and the standard values of the average human GDP in 2020 and 2030 are respectively determined as 80848 yuan/person and 115613 yuan/person. COD discharge amount of the ten thousand-membered GDP was determined to be 0.4kg, 0.36kg and 0.24kg as standard values in 2015, 2020 and 2030.
In addition, the present embodiment uses a "basic coordination" level achieved by 0.85 or more as a condition for satisfying load balancing, i.e., C min Take 0.85.
And 4, in the fitness function construction, each target weight coefficient in the embodiment is 1/3.
And 5, carrying out solution on the parameters in the model, wherein the calculation result is shown in a table 14.
Table 14 water resource load bearing capability calculation results table
The foregoing has shown and described the basic principles, principal features and advantages of the invention. It will be appreciated by persons skilled in the art that the above embodiments are not intended to limit the invention in any way, and that all technical solutions obtained by means of equivalent substitutions or equivalent transformations fall within the scope of the invention.
The present invention is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the technical scope of the present invention are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.
Claims (2)
1. A water resource bearing capacity quantization method based on load balancing is characterized by comprising the following steps: the method comprises the following steps:
1) Based on a load balancing theory, establishing a water resource system, a social and economic system and an ecological environment system interrelated and restrained water resource bearing strengthening model, and respectively determining an objective function of the model from three aspects of social scale, economic development and environmental protection;
2) Constraint conditions of the water resource bearing capacity quantization model comprise four aspects of water resource quantity constraint, socioeconomic constraint, water environment constraint and load balancing constraint;
3) After converting the multi-objective problem of the water resource bearing capacity quantification model into a single-objective problem, solving by adopting a genetic algorithm, thereby obtaining the population scale and the socioeconomic development scale which can be borne by the water resource system under different conditions;
step 1), the objective function includes:
1.1 Social goal: targeting regional loadable population values maximum:
in the method, in the process of the invention,
total population for the kth partition;
1.2 Economic goal: the maximum production total value in the country is as follows:
in the method, in the process of the invention,
a domestic production total value for the kth partition;
1.3 Environmental objective): targeting the minimum emission of regional chemical oxygen demand:
in the method, in the process of the invention,
COD discharge amount of the kth zone;
in the step 2), the water resource quantity constraint comprises two aspects that the water consumption in each partition is not more than the available unit water resource quantity and the sum of the water consumption of each water user in each partition is not more than the available area water resource quantity:
the internal water consumption of each partition is not more than the constraint of the available unit water resource:
T ZHSK ≤T KLYK ；
wherein T is ZHSK Inside each partitionConsumption of water, T KLYK The available amount of unit water resources;
the sum of the water consumption of each water user in each partition is not greater than the constraint of the available water resource in the area:
T ZHS ≤T KLY ；
wherein T is KLY The total available water resource amount; t (T) ZHS Is the total water consumption;
the method for calculating the total available water resource comprises the following steps:
T KLY ＝T DBS +T DXS -Δ RE ；
wherein T is DBS T is the available amount of surface water resources DXS Delta as the available amount of groundwater resources RE Repeating the amount of the surface water resource and the underground water resource;
the total water consumption calculating method comprises the following steps:
in the method, in the process of the invention,
rural domestic water consumption for kth zone, </i >>
Is town domestic water consumption of the kth subarea,
for the kth zone-water consumption, -/-, water production>
For the k-th partition, water consumption, < >>
For the k-th zone three-producing water consumption,
the ecological water consumption outside the river channel of the kth subarea; wherein->
Water consumption for the k-th partition planting industry,/->
For the k-th zone fishery water consumption, < >>
Water is consumed for the kth subarea animal husbandry;
in the step 2), the socioeconomic constraint comprises two aspects of industrial structure constraint and people average GDP constraint:
the industrial structure constraint is that the ratio of the increased value of each industrial production value of each partition to the total domestic production value of the area cannot exceed the upper limit and the lower limit of the increased value:
in the method, in the process of the invention,
total domestic production value for kth partition, < >>
Increasing the value of the ith industry production value for the kth partition,/->
Increasing the lower limit of the ratio of the value to the value of the i-th industrial production value for the kth partition,/->
Increasing the upper limit of the value proportion for the ith industry production value of the kth partition;
the method for calculating the domestic production total value of the kth partition comprises the following steps:
in the method, in the process of the invention,
water consumption for ten thousand yuan per yield for the kth zone, < >>
The water consumption is ten thousand yuan for the k-th partition,
the water consumption is ten thousand yuan for the third production of the kth partition;
the people average GDP is restricted, and the people average GDP of each partition is larger than or equal to the lower limit value of the people average GDP:
in the method, in the process of the invention,
for the total population of the kth partition, A perGDP Is the lower limit value of the average person GDP;
the total population calculation method of the kth partition comprises the following steps:
in the method, in the process of the invention,
for the k-th partition rural population, +.>
Town population for kth computing unit, < ->
For the k-th subarea rural people daily water consumption, < >>
Daily water consumption for people in towns in the kth subarea;
in step 2), the water environment constraint: the total COD discharged by all pollution sources is less than or equal to the total COD limit discharge:
in the method, in the process of the invention,
point source emissions for the kth zone, +.>
For the k-th subarea non-point source pollution emission quantity, RECOD k The total COD of the kth subarea is limited;
the k-th partition point source emission amount calculating method comprises the following steps:
wherein:
urban domestic water consumption for kth computing unitRate of->
A third industrial water consumption rate for the kth computing unit,
for the kth computing unit, the second industrial water consumption rate, A CLXS Is a sewage treatment coefficient A SHCOD Is COD concentration of town domestic water after sewage treatment, A GYCOD COD concentration of the second industrial water after sewage treatment;
the k-th partition area source pollution emission amount calculating method comprises the following steps:
wherein:
calculating the per mu water consumption of the unit planting industry, the per mu water consumption of the fishery industry and the water consumption of the unit livestock for the kth calculation unit planting industry; a is that ZZYCOD 、A YYCOD 、A XMYCOD 、A NCCOD COD generation coefficients of planting industry, fishery industry, animal husbandry and rural life are respectively A ZZYin 、A YYin 、A XMYin 、A NCSHin The river coefficients of the COD of the rural life are respectively the planting industry, the fishery industry, the animal husbandry and the planting industry;
in step 2), the load balancing constraint: the coordination degree of the social and economic system and the ecological environment system meets the requirement when the load balancing constraint is met, so that the stable development of the whole system is ensured;
wherein: c is the degree of coordination, and the value is between 0 and 1; n is an adjustment coefficient; i 1 Evaluation for coordinated development of socioeconomic subsystemsA price index; i 2 The coordinated development evaluation index of the ecological environment subsystem; c (C) min The lowest value of the coordination level is required to be reached;
the load balancing constraint selects the average-human GDP as a coordinated development evaluation index of the regional socioeconomic subsystem, and the average-human GDP is a forward index; the COD emission of the ten thousand-yuan GDP is selected as a coordinated development evaluation index of an ecological environment subsystem in the region and is a negative index; performing dimensionality removal treatment on the indexes by adopting a threshold value method, respectively reflecting the indexes of the social economic development and the ecological environment condition, and performing index standardization by adopting different formulas, wherein the formulas are as follows:
positive indicators, i.e. when the greater the indicator value the better the efficacy,
I i ＝C i /S i ；
negative indicators, i.e. when the efficacy is better the smaller the indicator value,
I i ＝S i /C i ；
wherein I is i Is the index value of index i, S i Is the standard value of index i, C i Is the actual value of index i;
the coordination degree in the load balancing constraint is divided into 5 coordination levels of extremely uncoordinated, barely coordinated, basic coordinated and coordinated according to the magnitude of the value, wherein the coordination degree in the extremely uncoordinated level is more than or equal to 0 and less than or equal to 0.49; c is more than or equal to 0.5 and less than or equal to 0.74 in the uncoordinated level; the coordination degree in the marginal coordination grade is more than or equal to 0.75 and less than or equal to 0.84; the coordination degree C in the basic coordination grade is more than or equal to 0.85 and less than or equal to 0.94; the coordination degree C in the coordination grade is more than or equal to 0.95 and less than or equal to 1.
2. The method for quantifying water resource load bearing capacity based on load balancing according to claim 1, wherein: in the step 3), the specific processing is as follows:
to achieve the optimal level f for the nth target nbest The corresponding standard value is set to 1, the worst level f nworst And setting the corresponding standard value to be 0.01, and calculating the corresponding standard value by each target according to an interpolation formula:
for social and economic goals, solve the maximization problem:
f′ n (x)＝f n (x)/f nbest ；
for an environmental objective solution objective, solving a minimization problem:
f′ n (x)＝1-f n (x)/f nbest ；
converting each target into a standard value between 0 and 1, and then converting the multiple targets into a single-target construction fitness function in a weighted summation mode:
in the method, in the process of the invention,
ρ n is the target weight coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910370288.4A CN110210710B (en) | 2019-05-06 | 2019-05-06 | Water resource bearing capacity quantization method based on load balancing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910370288.4A CN110210710B (en) | 2019-05-06 | 2019-05-06 | Water resource bearing capacity quantization method based on load balancing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110210710A CN110210710A (en) | 2019-09-06 |
| CN110210710B true CN110210710B (en) | 2023-06-27 |
Family
ID=67786854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910370288.4A Active CN110210710B (en) | 2019-05-06 | 2019-05-06 | Water resource bearing capacity quantization method based on load balancing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110210710B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111738554B (en) * | 2020-05-25 | 2024-01-19 | 华北水利水电大学 | Water consumption-based mountain water Lin Tianhu grass balanced water resource allocation method for each system |
| CN112465206B (en) * | 2020-11-20 | 2023-01-20 | 广东工业大学 | Water resource optimization configuration and optimization result evaluation feedback method |
| CN112966902B (en) * | 2021-02-04 | 2022-05-17 | 长江勘测规划设计研究有限责任公司 | Water pollution load distribution method considering region fairness and pollution source treatment difference |
| CN112966866A (en) * | 2021-03-05 | 2021-06-15 | 长春工程学院 | Water resource bearing capacity evaluation method based on optimized configuration |
| CN113343413B (en) * | 2021-04-22 | 2023-06-30 | 中国环境科学研究院 | Water environment bearing capacity evaluation method, device, equipment and medium |
| CN113344440B (en) * | 2021-06-30 | 2023-05-23 | 西北农林科技大学 | Construction method of underground water bearing capacity assessment model |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930954A (en) * | 2016-04-08 | 2016-09-07 | 中国水利水电科学研究院 | Three-layer risk assessment-based drought disaster risk coping method |
| CN109377089A (en) * | 2018-11-20 | 2019-02-22 | 河海大学 | A kind of evaluation method of the water resource loaded state based on coordinated development degree |
| CN109657933A (en) * | 2018-11-29 | 2019-04-19 | 河海大学 | A kind of regional water resources environment loaded state method for early warning |
| CN109670679A (en) * | 2018-11-29 | 2019-04-23 | 河海大学 | Six grades of water resource loaded state method for early warning of three classes based on Variable Fuzzy method |
-
2019
- 2019-05-06 CN CN201910370288.4A patent/CN110210710B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930954A (en) * | 2016-04-08 | 2016-09-07 | 中国水利水电科学研究院 | Three-layer risk assessment-based drought disaster risk coping method |
| CN109377089A (en) * | 2018-11-20 | 2019-02-22 | 河海大学 | A kind of evaluation method of the water resource loaded state based on coordinated development degree |
| CN109657933A (en) * | 2018-11-29 | 2019-04-19 | 河海大学 | A kind of regional water resources environment loaded state method for early warning |
| CN109670679A (en) * | 2018-11-29 | 2019-04-23 | 河海大学 | Six grades of water resource loaded state method for early warning of three classes based on Variable Fuzzy method |
Non-Patent Citations (3)
| Title |
|---|
| 基于多层次模糊综合评判法的南宁市水资源承载能力现状评价;熊健;《中国优秀硕士学位论文全文数据库工程科技II辑》;20190115;第42-50页 * |
| 水资源配置对水资源承载能力影响研究;俞祎波等;《人民黄河》;20180731;正文第19-23、26页 * |
| 熊健.基于多层次模糊综合评判法的南宁市水资源承载能力现状评价.《中国优秀硕士学位论文全文数据库工程科技II辑》.2019, * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110210710A (en) | 2019-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110210710B (en) | Water resource bearing capacity quantization method based on load balancing | |
| CN107301472A (en) | Distributed photovoltaic planing method based on scene analysis method and voltage-regulation strategy | |
| Jin-Yan et al. | Water resources allocation model based on ecological priority in the arid region | |
| CN106410781B (en) | A kind of power consumer demand response Potential Determining Method | |
| Yang et al. | Benefit evaluation of East Route Project of South to North Water Transfer based on trapezoid cloud model | |
| CN109214095A (en) | Electric car charge and discharge Multiobjective Optimal Operation method | |
| CN109598408B (en) | Annual water quantity scheduling plan compilation method considering water use fairness and importance | |
| CN111092451B (en) | Photovoltaic consumption early warning method based on power distribution network grid | |
| CN108649556A (en) | A kind of electric network emergency scene random optimization dispatching method | |
| Lin et al. | Comparative analysis of environmental carrying capacity of the Bohai Sea Rim area in China | |
| Yan et al. | Evaluation and prediction of water resources carrying capacity in Jiangsu Province, China | |
| CN108879659A (en) | A kind of evaluation method of elasticity power distribution network flexibility | |
| CN111428970A (en) | Large-scale hydropower station group trans-provincial delivery capacity analysis model and solving method | |
| CN104134102B (en) | Long-term electricity needs distribution forecasting method in power network based on LEAP models | |
| Song et al. | Optimal water allocation scheme in integrated water-ecosystem-economy system | |
| CN113240330A (en) | Multi-dimensional value evaluation method and scheduling strategy for demand side virtual power plant | |
| CN110310035A (en) | Active distribution network integrated evaluating method, device and storage medium | |
| CN114677064B (en) | Cascade reservoir scheduling decision support method coupling optimality and stability | |
| Bai et al. | Allocating total emission pollutant control based on water environmental carrying capacity: model establishment and case study | |
| CN115293644A (en) | Hybrid time scale hydrogen-electricity combined energy storage system planning method and system | |
| CN115239117A (en) | Classification and feature selection method based on multi-objective optimization | |
| CN114297913A (en) | Mining area water resource allocation method based on multi-objective optimization | |
| CN114493227A (en) | New energy scheduling evaluation method and medium based on colluding fuzzy-TOPSIS | |
| CN109255102B (en) | Photovoltaic power generation maximum grid-connected capacity evaluation method based on convex relaxation | |
| Lin et al. | Conception and evaluation methodology of water resources carrying capacity based on three-level analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |