Regional water utilization right initial allocation method

Abstract

The invention discloses a regional water utilization right initial allocation method, which comprises the following steps of: collecting regional basic data, and determining the total amount of distributable water utilization right of a region; obtaining the distributable non-rigid water use total amount of the region; proposing a non-rigid water utilization right distribution principle to set a target function, and assigning a distribution index weight in the target function; a water utilization right initial distribution model based on an NSGA-II algorithm is constructed and solved, and the model is used as a regional non-rigid water utilization right result of model distribution; and adding the regional non-rigid water utilization right result allocated by the model to the regional rigid water utilization right result to obtain a final allocation result of initial allocation of the regional water utilization right. The method has the beneficial effects that regional water saving and efficient utilization of water resources can be effectively promoted through the distribution method, and a water utilization right distribution scheme can be better formulated.

Description

TECHNICAL FIELD

[0001] The present application relates to water resources optimal allocation technology, and particularly relates to a regional water right initial allocation method. BACKGROUND

[0002] The regional water right initial allocation is a research topic that needs to be solved urgently at present, is related to the economic and social development and benefit distribution of each region, involves water resource endowment, society and economy and other factors, and is a complex decision-making process involving multiple targets and mutual restraint between targets. In the existing related researches, the rigid water right demand is not considered, and the ecological protection is considered less, so that the water right allocation scheme cannot find a suitable balance point between targets.

[0003] Therefore, the present application provides a regional water right initial allocation method. SUMMARY

[0004] In view of the deficiencies in the prior art, the present application provides a regional water right initial allocation method, which is a regional water right initial allocation method considering rigid water right and based on NSGA- Under the premise of ensuring basic water, the non-rigid water right is allocated according to the allocation strategy, the fair rights and interests between water regions, water resource utilization efficiency and ecological protection can be considered in the allocation, and the allocation strategy can be adjusted as needed, so that the method can be applied to various scenes.

[0005] To achieve the above purpose, the technical scheme adopted by the present application is as follows: a regional water right initial allocation method, comprising the following steps: Step S1, collecting regional basic data to determine the total amount of the regional allocable water right; Step S2, providing a method for calculating the total amount of the regional allocable non-rigid water right, to calculate the total amount of the regional rigid water right and the total amount of the regional allocable non-rigid water right; Step S3, providing non-rigid water right allocation principles, including efficiency principle, ecological protection principle, fairness principle and water saving priority principle, setting target functions for the non-rigid water right allocation principles respectively, and assigning weights to the allocation indexes in the target functions; Step S4, constructing a non-rigid water right initial allocation model based on NSGA- Step S4, constructing a non-rigid water right initial allocation model based on NSGA- Step S5, adding the total amount of the regional rigid water right to the regional non-rigid water right result allocated by the non-rigid water right initial allocation model as the final allocation result of the regional water right initial allocation.

[0006] Furthermore, step S1 involves collecting basic regional data to determine the allocatable water rights data for the region; specifically: Step S11: Regional basic data includes the region's multi-year average total water resources, base year total water consumption, base year ecological and environmental water consumption, base year unconventional water source water consumption, base year industrial water consumption, river and lake pollution carrying capacity, base year forest coverage rate, base year population, base year irrigated area, base year agricultural irrigation water effective utilization coefficient, base year GDP value, base year industrial added value, current water rights allocation, base year domestic water consumption index, base year average irrigation water consumption per mu, and base year water consumption per 10,000 yuan of industrial added value. Step S12: Determine the water rights data in the regional basic data as the total amount of water rights that can be allocated in the region.

[0007] Furthermore, step S2 proposes a method for calculating the total amount of allocable non-rigid water rights in a region, which calculates the total amount of rigid water rights and the total amount of allocable non-rigid water rights in the region, specifically as follows: Step S21: The total amount of rigid water rights is calculated using the quota method. The total amount of rigid water rights includes rigid domestic water rights, rigid agricultural irrigation water rights, and rigid industrial water rights. The rigid domestic water use index, rigid irrigation water quota, and rigid industrial water use quota are calculated respectively. Step S22: Subtract the total amount of rigid water rights from the total amount of allocable water rights in the region to obtain the total amount of allocable non-rigid water rights. The calculation formula is as follows: ; ; ; In the formula: The total amount of rigid water rights in region i; Let be the population of region i. The rigid domestic water consumption quota for the region; Areai is the irrigated area of ​​the i-th region; For regional rigid irrigation water quotas; It is the industrial added value of region i; A rigid industrial water consumption quota for the regional baseline year; WU i Let be the total amount of non-rigid water rights that can be allocated in region i; Let be the total amount of water rights that can be allocated in region i; W t denoted as the total amount of non-rigid water rights that can be allocated in a region; n represents the number of regions to be allocated.

[0008] Furthermore, step S3 proposes non-rigid water rights allocation principles, including the principles of efficiency, ecological protection, fairness, and water conservation priority. Objective functions are set for each of these principles, and weights are assigned to the allocation indicators within the objective functions. Specifically: Step S31, the objective function of the efficiency principle; The principle of high efficiency is used to allocate water rights based on the GDP value per unit of water and the industrial added value per unit of water. Water use efficiency is quantified using the following formula: ; ; ; Where: WGR i Let be the unilateral water GDP value of region i; GDP i WP represents the baseline year GDP value for region i; i The baseline annual water consumption for region i; The added value of water industry per unit volume in region i; The baseline annual industrial water consumption for region i; To satisfy the principle of high efficiency, WD i The initial non-rigid water allocation for region i; , These are the weighting coefficients for the GDP per unit of water and the industrial added value per unit of water, respectively, in the principle of high efficiency. + =1; Efficiency principle satisfied The closer a value is to 1, the higher the level of satisfaction. Step S32, the objective function of ecological protection principles; In the principle of ecological protection, water rights are allocated based on water resource utilization rate, ecological environment water use rate, forest coverage rate, and river and lake pollution carrying capacity. The principle of ecological protection is quantified using the following formula: ; In the formula: To assess the degree of satisfaction with ecological protection principles, W i Let We be the multi-year average total water resources of region i; i Let be the baseline annual ecological water consumption for region i; For i Let be the forest cover rate of region i. The pollution carrying capacity of rivers and lakes in region i is calculated based on the chemical oxygen demand limit for total pollution discharge. , , , These are the weighting coefficients for water resource utilization rate, ecological environment water use rate, forest coverage rate, and river and lake pollution carrying capacity, respectively, in the principle of ecological protection. + + + =1; Satisfaction rate of ecological protection principles The closer a value is to 1, the higher the level of satisfaction. Step S33, the objective function of the fairness principle; The principle of fairness is primarily based on the current allocation of water rights, population size, and irrigated area within a region. It is quantified through a comparative analysis of the ratio of initial allocated water volume to current water rights, the initial allocated water volume per capita, and the initial allocated water volume per mu (unit of land area). ; In the formula: To satisfy the principle of fairness; Pop i Let be the population of region i; Area i It is the irrigated area of ​​region i; , , These are the weighting coefficients for water rights, population size, and irrigated area allocated according to the principle of fairness. + + =1; Satisfaction of the fairness principle The closer a value is to 1, the higher the level of satisfaction. Step S34, the objective function of the water conservation priority principle; The principle of prioritizing water conservation is mainly based on water consumption per 10,000 yuan of industrial added value, utilization rate of unconventional water sources, and effective utilization coefficient of agricultural irrigation water, and is quantified using the following formula: ; In the formula: To ensure compliance with the principle of prioritizing water conservation; Water consumption per 10,000 yuan of industrial added value in region i; The baseline year's unconventional water resource utilization for region i; The effective utilization coefficient of agricultural irrigation water in the i-th region is the baseline year. , , These are the weighting coefficients for water consumption per 10,000 yuan of industrial added value, utilization rate of unconventional water sources, and effective utilization coefficient of agricultural irrigation water, respectively, in the principle of prioritizing water conservation. + + =1; Satisfaction rate of the water conservation priority principle The closer a value is to 1, the higher the level of satisfaction. Step S35: Assign values ​​to all weight coefficients; By using subjective or objective assignment methods , , , , , , , , , , , Perform the assignment.

[0009] Furthermore, in step S4, NSGA - The algorithm's calculation process is as follows: Step S41: Initialize the population. After determining the upper and lower limits of the ratio of the initial allocated water volume to the current water rights, generate random numbers to randomly generate the initial population with a population size of N, and calculate the objective function value for each individual. Step S42: Calculate the dominance level and distance for each individual, and sort the individuals according to their dominance level; Step S43: Select the top N / 2 individuals in the hierarchy to produce offspring. If the number of individuals in the hierarchy where the N / 2 individual is located plus the number of individuals in the preceding hierarchy exceeds N / 2, sort by distance and eliminate individuals with smaller distances. Step S44: Randomly select N / 2 individuals generated in step S43 according to certain rules to hybridize or mutate and produce offspring, and calculate the objective function value of the offspring; Step S45: Calculate the dominance level and distance of the offspring individuals, sort the parents and offspring together according to the dominance level, and select the top N individuals to keep the population size unchanged; within the same level, individuals with larger distances are better than individuals with smaller distances. Step S46: Repeat steps S43-S45 until the set number of generations is reached, then output and save the results.

[0010] The beneficial effects of this invention are: firstly, it allocates rigid water rights to each region, ensuring the basic water needs of each region; then, it redistributes non-rigid water rights, constructing a system based on NSGA. The algorithm's non-rigid water rights allocation method can consider multiple allocation principles in the initial allocation of water rights, and can also adjust the allocation strategy by assigning values ​​to the index weights in the objective function, which helps decision-makers to better formulate water rights allocation schemes; the above allocation method can effectively promote regional water conservation and efficient use of water resources. Attached Figure Description

[0011] Figure 1 This is a schematic diagram of the operation of the present invention. Detailed Implementation

[0012] Taking the initial allocation of water rights in various cities and prefectures of XX province as an example, this method is used for analysis.

[0013] Step S1: Collect basic regional data and determine the total amount of water rights that can be allocated in the region; Step S11: Regional basic data includes the region's multi-year average total water resources, base year total water consumption, base year ecological and environmental water consumption, base year unconventional water source water consumption, base year industrial water consumption, river and lake pollution carrying capacity, base year forest coverage rate, base year population, base year irrigated area, base year agricultural irrigation water effective utilization coefficient, base year GDP value, base year industrial added value, current water rights allocation, base year domestic water consumption index, base year average irrigation water consumption per mu, and base year water consumption per 10,000 yuan of industrial added value. Table 1. Basic Data of Various Cities in XX Province

[0014] Step S12: Determine the total amount of water rights that can be allocated in the region. According to relevant documents, the total water consumption control target for all cities and prefectures in XX Province in 2030 is 25.863 billion m³. 3 The current effective plans all use this value as the total water rights for all cities and prefectures in XX province, therefore this value is used as the total water rights allocated this time. Step S2 proposes a method for calculating the total amount of non-rigid water rights that can be allocated in a region, and calculates the total amount of rigid water rights and the total amount of non-rigid water rights that can be allocated in the region. The rigid domestic water consumption quota for various cities and prefectures in XX Province adopts the minimum per capita domestic water consumption of rural residents in various cities and prefectures of XX Province in 2022, which is 89L / d; the rigid irrigation water consumption quota adopts the minimum net irrigation quota per mu (unit of land area) in various cities and prefectures of XX Province in 2022, which is 360m³. 3 / mu; the rigid industrial water quota adopts the minimum water consumption per 10,000 yuan of industrial added value in various cities and prefectures of XX Province in 2022, which is 11m³. 3 / 10,000 yuan.

[0015] The total amount of rigid water rights and the total amount of allocable non-rigid water rights in various cities and prefectures of XX province were calculated.

[0016] Table 2. Calculation Results of Total Allocatable Non-Rigid Water Rights in Various Cities of XX Province

[0017] Step S3: Set target principles. In this embodiment, the target principles are selected simultaneously: efficiency principle, ecological protection principle, fairness principle, and water conservation priority principle. The weights of the indicators in the objective function are assigned. In this embodiment, the weights are assigned subjectively according to the importance of the indicators, resulting in the weight table 3.

[0018] Table 3. Results of Objective Function Weights under Each Allocation Principle

[0019] Step S4, build based on NSGA - The algorithm's initial water rights allocation model is derived and solved. Perform population initialization in NSGA - In the algorithm, each individual contains 11 variables. i = 1, 2, ..., 11 The ratio of initially allocated non-rigid water rights to allocable non-rigid water rights. A population size of 200 and 500 iterations were selected for simulation calculations; the objective functions for each calculation were: The objective function of the efficiency principle: ; Objective function of ecological protection principles: ; The objective function of the fairness principle: ; The objective function of the water conservation priority principle: ; After reaching the required number of iterations, the model calculates multiple Pareto optimal solutions. The average of the allocation schemes for all Pareto optimal solutions is then used as the final result of the model.

[0020] Table 4 Initial Allocation Results of Non-Rigid Water Rights in Various Cities of XX Province

[0021] The results of the satisfaction with each objective principle are shown in Table 5.

[0022] Table 5. Satisfaction of Initial Allocation Results of Non-Rigid Water Rights in Various Cities of XX Province with Each Objective Principle

[0023] The results above show that the non-rigid water rights allocation does not satisfy the allocation principles well. This is because too many objectives are selected at the same time, and the mutual constraints between the objectives make it impossible for all objectives to be well satisfied.

[0024] Step S5: Add the results of the non-rigid water rights allocation in the region to the results of the rigid water rights allocation in the region as the initial and final allocation results of the regional water rights allocation.

[0025] The initial allocation results of water rights in various cities and prefectures of XX province are shown in Table 6.

[0026]

[0027] Table 6 Initial allocation results of water rights in various cities and prefectures of XX Province.

Claims

1. A method for initial allocation of regional water rights, characterized in that: Includes the following steps: Step S1: Collect basic regional data and determine the total amount of water rights that can be allocated in the region; Step S2 proposes a method for calculating the total amount of non-rigid water rights that can be allocated in a region, and calculates the total amount of rigid water rights and the total amount of non-rigid water rights that can be allocated in the region. Step S3 proposes non-rigid water rights allocation principles, including the principles of efficiency, ecological protection, fairness, and water conservation priority. Objective functions are set for each non-rigid water rights allocation principle, and the weights of the allocation indicators within the objective functions are assigned. Step S4, build based on NSGA - The algorithm is used to model and solve the initial allocation model of non-rigid water rights. The average of all Pareto optimal solutions is calculated, and the resulting solution is used as the result of the allocation of non-rigid water rights in the region by the initial allocation model of non-rigid water rights. Step S5: Add the total amount of rigid water rights in the region to the regional non-rigid water rights allocated by the initial allocation model to obtain the final allocation result of the initial allocation of regional water rights.

2. The method for initial allocation of regional water rights according to claim 1, characterized in that: Step S1 involves collecting basic regional data and determining the allocatable water rights data for the region; specifically: Step S11: Regional basic data includes the region's multi-year average total water resources, base year total water consumption, base year ecological and environmental water consumption, base year unconventional water source water consumption, base year industrial water consumption, river and lake pollution carrying capacity, base year forest coverage rate, base year population, base year irrigated area, base year agricultural irrigation water effective utilization coefficient, base year GDP value, base year industrial added value, current water rights allocation, base year domestic water consumption index, base year average irrigation water consumption per mu, and base year water consumption per 10,000 yuan of industrial added value. Step S12: Determine the water rights data in the regional basic data as the total amount of water rights that can be allocated in the region.

3. The method for initial allocation of regional water rights according to claim 2, characterized in that: Step S2 proposes a method for calculating the total amount of allocable non-rigid water rights in a region, and calculates the total amount of rigid water rights and the total amount of allocable non-rigid water rights in the region, specifically as follows: Step S21: The total amount of rigid water rights is calculated using the quota method. The total amount of rigid water rights includes rigid domestic water rights, rigid agricultural irrigation water rights, and rigid industrial water rights. The rigid domestic water use index, rigid irrigation water quota, and rigid industrial water use quota are calculated respectively. Step S22: Subtract the total amount of rigid water rights from the total amount of allocable water rights in the region to obtain the total amount of allocable non-rigid water rights. The calculation formula is as follows: ； ； ； In the formula: The total amount of rigid water rights in region i; Let be the population of region i. The rigid domestic water consumption quota for the region; Areai is the irrigated area of ​​the i-th region; For regional rigid irrigation water quotas; It is the industrial added value of region i; A rigid industrial water consumption quota for the regional baseline year; WU i Let be the total amount of non-rigid water rights that can be allocated in region i; Let be the total amount of water rights that can be allocated in region i; W t denoted as the total amount of non-rigid water rights that can be allocated in a region; n represents the number of regions to be allocated.

4. The method for initial allocation of regional water rights according to claim 3, characterized in that: Step S3 proposes non-rigid water rights allocation principles, including the principles of efficiency, ecological protection, fairness, and water conservation priority. Objective functions are set for each of these principles, and weights are assigned to the allocation indicators within the objective functions. Specifically: Step S31, the objective function of the efficiency principle; The principle of high efficiency is used to allocate water rights based on the GDP value per unit of water and the industrial added value per unit of water. Water use efficiency is quantified using the following formula: ； ； ; Where: WGR i Let be the unilateral water GDP value of region i; GDP i WP represents the baseline year GDP value for region i; i The baseline annual water consumption for region i; The added value of water industry per unit volume in region i; The baseline annual industrial water consumption for region i; To satisfy the principle of high efficiency, WD i The initial non-rigid water allocation for region i; , These are the weighting coefficients for the GDP per unit of water and the industrial added value per unit of water, respectively, in the principle of high efficiency. + =1; Efficiency principle satisfied The closer a value is to 1, the higher the level of satisfaction. Step S32, the objective function of ecological protection principles; In the principle of ecological protection, water rights are allocated based on water resource utilization rate, ecological environment water use rate, forest coverage rate, and river and lake pollution carrying capacity. The principle of ecological protection is quantified using the following formula: ； In the formula: To assess the degree of satisfaction with ecological protection principles, W i Let We be the multi-year average total water resources of region i; i Let be the baseline annual ecological water consumption for region i; For i Let be the forest cover rate of region i. The pollution carrying capacity of rivers and lakes in region i is calculated based on the chemical oxygen demand limit for total pollution discharge. , , , These are the weighting coefficients for water resource utilization rate, ecological environment water use rate, forest coverage rate, and river and lake pollution carrying capacity, respectively, in the principle of ecological protection. + + + =1; Satisfaction rate of ecological protection principles The closer a value is to 1, the higher the level of satisfaction. Step S33, the objective function of the fairness principle; The principle of fairness is primarily based on the current allocation of water rights, population size, and irrigated area within a region. It is quantified through a comparative analysis of the ratio of initial allocated water volume to current water rights, the initial allocated water volume per capita, and the initial allocated water volume per mu (unit of land area). ； In the formula: To satisfy the principle of fairness; Pop i Let be the population of region i; Area i It is the irrigated area of ​​region i; , , These are the weighting coefficients for water rights, population size, and irrigated area allocated according to the principle of fairness. + + =1; Satisfaction of the fairness principle The closer a value is to 1, the higher the level of satisfaction. Step S34, the objective function of the water conservation priority principle; The principle of prioritizing water conservation is mainly based on water consumption per 10,000 yuan of industrial added value, utilization rate of unconventional water sources, and effective utilization coefficient of agricultural irrigation water, and is quantified using the following formula: ; In the formula: To ensure compliance with the principle of prioritizing water conservation; Water consumption per 10,000 yuan of industrial added value in region i; The baseline year's unconventional water resource utilization for region i; The effective utilization coefficient of agricultural irrigation water in the i-th region is the baseline year. , , These are the weighting coefficients for water consumption per 10,000 yuan of industrial added value, utilization rate of unconventional water sources, and effective utilization coefficient of agricultural irrigation water, respectively, in the principle of prioritizing water conservation. + + =1; Satisfaction rate of the water conservation priority principle The closer a value is to 1, the higher the level of satisfaction. Step S35: Assign values ​​to all weight coefficients; By using subjective or objective assignment methods , , , , , , , , , , , Perform the assignment.

5. The method for initial allocation of regional water rights according to claim 4, characterized in that: NSGA in step S4 - The algorithm's calculation process is as follows: Step S41: Initialize the population. After determining the upper and lower limits of the ratio of the initial allocated water volume to the current water rights, generate random numbers to randomly generate the initial population with a population size of N, and calculate the objective function value for each individual. Step S42: Calculate the dominance level and distance for each individual, and sort the individuals according to their dominance level; Step S43: Select the top N / 2 individuals in the hierarchy to produce offspring. If the number of individuals in the hierarchy where the N / 2 individual is located plus the number of individuals in the preceding hierarchy exceeds N / 2, sort by distance and eliminate individuals with smaller distances. Step S44: Randomly select N / 2 individuals generated in step S43 according to certain rules to hybridize or mutate and produce offspring, and calculate the objective function value of the offspring; Step S45: Calculate the dominance level and distance of the offspring individuals, sort the parents and offspring together according to the dominance level, and select the top N individuals to keep the population size unchanged; within the same level, individuals with larger distances are better than individuals with smaller distances. Step S46: Repeat steps S43-S45 until the set number of generations is reached, then output and save the results.

Images