CN113887073A - Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium - Google Patents

Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium Download PDF

Info

Publication number
CN113887073A
CN113887073A CN202111231979.XA CN202111231979A CN113887073A CN 113887073 A CN113887073 A CN 113887073A CN 202111231979 A CN202111231979 A CN 202111231979A CN 113887073 A CN113887073 A CN 113887073A
Authority
CN
China
Prior art keywords
water
soil
basic configuration
configuration unit
resources
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.)
Pending
Application number
CN202111231979.XA
Other languages
Chinese (zh)
Inventor
秦天玲
王建伟
吕锡芝
刘姗姗
李晨昊
冯贱明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202111231979.XA priority Critical patent/CN113887073A/en
Publication of CN113887073A publication Critical patent/CN113887073A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Abstract

The invention provides a method and a system for optimal allocation of water and soil resources in a drainage basin, electronic equipment and a storage medium, and belongs to the technical field of optimal allocation of water and soil resources. The method comprises the following steps: determining a basic configuration unit for optimal configuration of water and soil resources of a drainage basin; determining an objective function of optimal allocation of watershed water and soil resources; determining constraint conditions for optimal configuration of water and soil resources of the drainage basin; and calculating by using a multi-objective optimization algorithm according to the objective function and the constraint condition to obtain an optimal allocation result of the watershed water and soil, and completing the optimal allocation of resources of the watershed water and soil. The method brings the land resources into the traditional water resource allocation, realizes the optimal allocation of water and soil resources under the mutual influence of the land resources, and simultaneously spreads the optimal allocation result of the land resources to the whole drainage basin, thereby providing technical support for the ecological civilization construction of the drainage basin and the construction of the 'mountain water forest field, lake, grass and sand' life community.

Description

Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium
Technical Field
The invention belongs to the technical field of water and soil resource optimal allocation, and particularly relates to a method and a system for optimal allocation of water and soil resources in a drainage basin, electronic equipment and a storage medium.
Background
At present, social economy develops rapidly, requirements of people on ecological environment are higher and higher, contradiction between social development water land and ecological water land is aggravated day by day, and sustainable development of society is influenced. In order to guarantee the sustainable development of human society and economy and build a living environment, water and soil resources need to be reasonably and optimally configured to meet the requirements of the sustainable development of human society.
At present, a great deal of research is carried out on water resource optimal allocation and land resource optimal allocation at home and abroad, the water resource optimal allocation is as large as a drainage basin and as small as a farmland, the land resource optimal allocation is mainly based on the farmland, and the research on the optimal allocation of water and soil resources from the drainage basin angle is relatively less. Aiming at the defects, a multi-objective water and soil resource optimal allocation method of the watershed scale is constructed, and a technical support system for the watershed ecological civilization construction is supplemented and expanded.
Disclosure of Invention
Aiming at the defects in the prior art, the method, the system, the electronic equipment and the storage medium for optimal allocation of the water and soil resources in the drainage basin provided by the invention can perform optimal allocation of the water and soil resources in the drainage basin from the perspective of the drainage basin, and provide technical support for ecological civilization construction of the drainage basin.
In order to achieve the above purpose, the invention adopts the technical scheme that:
the scheme provides a method for optimizing and configuring resources of water and soil in a drainage basin, which comprises the following steps:
s1, determining a basic configuration unit for optimal configuration of the water and soil resources of the drainage basin;
s2, determining a target function of the optimal allocation of the watershed water and soil resources according to the basic allocation unit of the optimal allocation of the watershed water and soil resources;
s3, determining constraint conditions of optimal allocation of the watershed water and soil resources;
and S4, calculating by using a multi-objective optimization algorithm according to the objective function and the constraint condition to obtain an optimal allocation result of the water and soil of the drainage basin, and completing the optimal allocation of the resources of the water and soil of the drainage basin.
The invention has the beneficial effects that: the invention relates to a 'quality-effect-raw' oriented water and soil resource optimal allocation method, which is characterized in that a land utilization prediction model and a multi-objective optimal allocation model are utilized, an optimization algorithm is combined for solving, the areas and the water consumption of different land utilization types of sub-basins in a basin range are finally obtained, and spatial distribution is carried out through a land utilization model, so that a spatial distribution result of water and soil resource optimal allocation is obtained. The method brings the land resources into the traditional water resource allocation, realizes the optimal allocation of water and soil resources under the mutual influence of the land resources, and simultaneously spreads the optimal allocation result of the land resources to the whole drainage basin, thereby providing technical support for the ecological civilization construction of the drainage basin and the construction of the 'mountain water forest field, lake, grass and sand' life community.
Further, the step S1 includes the following steps:
s101, hydrologic analysis is carried out on the basis of a basin digital elevation model and a basin exit point, and sub-basins of a research area are drawn;
s102, performing intersection processing on the sub-river basin and the administrative region to obtain a sub-river basin set administrative region basic unit;
s103, the areas of different land types in the sub-flow field administrative region covering basic units are counted to obtain basic configuration units for optimal configuration of water and soil resources in the river basin.
The beneficial effects of the further scheme are as follows: acquiring a sub-basin of a research area through hydrological analysis, and facilitating configuration research from the aspect of the topography and landform of the basin; through the intersection processing of the sub-watersheds and the administrative regions, the landform and the landform of the watershed and the administrative regions are considered in the configuration process; the areas of different land utilization types in the sub-river basin set administrative region are counted, the basic configuration unit is further refined to the land utilization type scale, and the basic configuration unit (the sub-river basin-administrative region-land utilization type) is finally obtained, so that the configuration from three angles of river basin, administrative region and land utilization type is facilitated, the three aspects are comprehensively considered, the configuration result considers both water resources and land resources, and the joint configuration of water and soil resources is realized.
Still further, the expression of the objective function of the watershed water and soil resource optimization configuration in step S2 is as follows:
Figure BDA0003316243070000031
Figure BDA0003316243070000032
Figure BDA0003316243070000033
Figure BDA0003316243070000034
Figure BDA0003316243070000035
Figure BDA0003316243070000036
Figure BDA0003316243070000037
Hij(X)＝-Pij×LnPij
Figure BDA0003316243070000038
Figure BDA0003316243070000039
Figure BDA00033162430700000310
HHij(Y)＝-PPij×LnPPij
Figure BDA0003316243070000041
wherein, minF1(w,L)、minF2(w,L)、minF3(w,L)、minF4(L) and maxF5(w, L) respectively representing a minimum objective function of total water shortage of a drainage basin, a minimum objective function of total pollutant discharge, a minimum objective function of ten-thousand-yuan water consumption, a maximum objective function of net primary productivity and a maximum function of water and soil resource allocation damping coefficient, wijkIndicating the water demand of the jth land utilization type to the kth water source in the ith basic configuration unit, LijArea G representing the jth type of land use in the ith basic configuration unitikThe water supply amount of the kth water source in the ith basic configuration unit is represented, beta represents a point source surface source pollution discrimination coefficient, the value of the point source surface source pollution discrimination coefficient is 0 or 1, 1 represents point source pollution, 0 represents surface source pollution, eijRepresents the concentration of pollutants in the discharge wastewater of the jth soil utilization type in the ith basic configuration unit, pijA sewage discharge coefficient, A, representing the jth soil utilization type in the ith basic configuration unitijRepresents the pollutant load of the jth land utilization type in the ith basic configuration unit, giRepresenting GDP output value of ith basic configuration unit, n representing total number of configuration units, p representing total number of land utilization types, m representing total number of water sources, NPPijRepresenting the net primary productivity of the jth type of land utilization in the ith basic configuration unit,
Figure BDA0003316243070000042
shows the water and soil resource matching kini coefficient, x, of the kth water source in the flow fieldijRepresents the cumulative proportion of the j-th land utilization type balance degree in the i-th basic configuration unit, yijRepresents the cumulative proportion, x, of the water consumption balance of the jth land utilization type in the ith basic configuration unitijAnd yijForm a space Lorentzian curve, u represents the cumulative number of the configuration units, JijCumulative structural balance, H, representing the jth type of land utilization in the ith basic configuration cellijEntropy of structural information, H, representing the jth type of land use in the ith basic configuration unitijmaxMaximum structure information entropy representing all land use types in all configuration units, Ln representing natural logarithm, Hij(X) Structure information entropy, P, representing the jth type of land utilization in the ith basic configuration UnitijRepresents the total area proportion of the jth land utilization type in the ith basic configuration unit, JJujRepresents the cumulative water consumption balance of the jth land utilization type in the ith basic configuration unit, JJijWater consumption balance, HH, representing the jth land utilization type in the ith basic configuration unitijEntropy of water quantity information, HH, representing the jth type of land utilization in the ith basic configuration unitijmaxEntropy of maximum water usage information, HH, representing all land use types in all configured unitsij(Y) entropy of water consumption information indicating the jth type of land use in the ith basic configuration unit, PPijThe water consumption of the jth land utilization type in the ith basic configuration unit is represented as the proportion of the total water consumption of the jth land utilization type in the whole flow field.
The beneficial effects of the further scheme are as follows: the multi-target function is established from five aspects of water quantity, water quality, water use efficiency, ecology and water and soil resource space matching degree, so that the consideration factors of the configuration result can more comprehensively meet the actual development requirement, and the problem of incomplete consideration factors of the configuration result caused by the optimization of a single target is avoided.
Still further, the expression of the constraint condition in step S3 is as follows:
Figure BDA0003316243070000051
Figure BDA0003316243070000052
Figure BDA0003316243070000053
Figure BDA0003316243070000054
Figure BDA0003316243070000055
wherein n represents the total number of configuration units, m represents the total number of land use types, GikIndicates the water supply available amount, W, of the kth water source of the ith basic configuration unitk(p) represents the total amount of available water resources of the kth water source when the incoming water frequency is p, SikIndicates the actual water supply amount, w, of the kth water source in the ith basic configuration unitijkThe water demand of the jth land utilization type to the kth water source in the ith basic configuration unit is represented, beta represents a point source and surface source pollution discrimination coefficient, the value is 0 or 1, 1 represents point source pollution, 0 represents surface source pollution, and e represents the surface source pollutionijRepresents the concentration of pollutants in the j soil utilization type discharged wastewater in the ith basic configuration unit, pijIn the ith basic configuration unitSewage discharge coefficient of jth land utilization type, wijkRepresenting the water demand of the jth land utilization type to the kth water source in the ith basic configuration unit, AijIndicating the pollutant load of the jth land utilization type of the ith basic configuration unit, LijAn area T representing the jth land utilization type in the ith basic configuration unitiRepresents the total amount of contamination consumed in the ith basic configuration unit, AiRepresents the total area, Ldown, of the ith basic configuration unitjIndicating the lower boundary, Lup, of the jth land utilization type red line on the ith basic configuration unitjAnd the upper boundary of the jth land utilization type in the drainage basin in the ith basic configuration unit is shown.
The beneficial effects of the further scheme are as follows: the constraint conditions are established through three aspects of water supply amount, pollutant carrying capacity and land utilization red line, so that the rapid solution of optimization calculation and the optimization result are more reasonable, the actual development requirements are met, and the method has actual guiding significance on the territory planning development.
Still further, the step S4 includes the steps of:
s401, unifying the objective functions into the minimum optimum;
s402, uniformly naming the optimized variables in a coding naming mode;
s403, calculating the water consumption of each configuration unit and the area of each land utilization type by using a multi-objective optimization algorithm according to the constraint conditions;
s404, inputting the water consumption of each configuration unit and the area of each land utilization type into the FLUS model to obtain the optimal configuration result of the water and soil in the drainage basin, and completing the optimal configuration of the resources of the water and soil in the drainage basin.
The beneficial effects of the further scheme are as follows: the target function is unified into the minimum optimum, so that the optimization, solution and calculation are facilitated; the optimization change is named uniformly by adopting a coding naming mode, so that the input, the output and the statistics of model data are facilitated in the optimization solving process; through optimization solution calculation, the water consumption in each configuration unit and the area of each land utilization type are obtained, and data support is provided for the space spread of land utilization; and the land use space distribution of the river basin optimal configuration result is favorably obtained by performing simulation calculation through the FLUS model.
The invention provides a system for optimizing and configuring resources of water and soil in a drainage basin, which comprises:
the basic configuration unit determining module is used for determining a basic configuration unit for optimal configuration of water and soil resources of the drainage basin;
the objective function determining module is used for determining an objective function of the optimal allocation of the watershed water and soil resources according to the basic allocation unit of the optimal allocation of the watershed water and soil resources;
the constraint condition determining module is used for determining constraint conditions of optimal configuration of water and soil resources of the drainage basin;
and the optimal configuration result calculation module is used for calculating to obtain an optimal configuration result of the watershed water and soil by using a multi-objective optimization algorithm according to the objective function and the constraint condition, and completing the optimal configuration of the resources of the watershed water and soil.
The invention has the beneficial effects that: the invention relates to a 'quality-effect-raw' oriented water and soil resource optimal allocation method, which is characterized in that a land utilization prediction model and a multi-objective optimal allocation model are utilized, an optimization algorithm is combined for solving, the areas and the water consumption of different land utilization types of sub-basins in a basin range are finally obtained, and spatial distribution is carried out through a land utilization model, so that a spatial distribution result of water and soil resource optimal allocation is obtained. The method brings the land resources into the traditional water resource allocation, realizes the optimal allocation of water and soil resources under the mutual influence of the land resources, and simultaneously spreads the optimal allocation result of the land resources to the whole drainage basin, thereby providing technical support for the ecological civilization construction of the drainage basin and the construction of the 'mountain water forest field, lake, grass and sand' life community.
The invention provides electronic equipment which comprises a memory, a processor and a computer program stored on the memory and running on the processor, wherein the processor executes the program to realize the optimal resource allocation method of the watershed water and soil.
The invention provides a computer-readable storage medium, which stores a computer program, wherein the computer program is executed by a processor to realize the method for optimally configuring the resources of the water and the soil in the drainage basin.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a network diagram of a water and soil resource system of the Sihefu river basin in this embodiment.
Fig. 3 is a schematic diagram of the optimal configuration result of the water and soil resources in the river basin in this embodiment.
Fig. 4 is a spatial layout view of the optimal arrangement result of the Sichuan river basin in this embodiment.
FIG. 5 is a schematic diagram of the system of the present invention.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
Example 1
The invention takes historical land utilization, landform, socioeconomic data and water resource utilization as the basis to predict water and soil resources and is used as constraint conditions, takes social, economic, environmental, ecological and spatial patterns as an objective function of optimal configuration, and carries out optimal solution by means of an NSGA-II optimization algorithm to obtain the water and soil resource distribution condition in each configuration unit in a basin range, as shown in figure 1, the invention provides a method for optimal configuration of water and soil resources in a basin, which comprises the following steps:
s1, determining a basic configuration unit for optimal configuration of water and soil resources in the drainage basin, wherein the implementation method comprises the following steps:
s101, hydrologic analysis is carried out on the basis of a basin digital elevation model and a basin exit point, and sub-basins of a research area are drawn;
s102, performing intersection processing on the sub-river basin and the administrative region to obtain a sub-river basin set administrative region basic unit;
s103, the areas of different land types in the sub-flow field administrative region covering basic units are counted to obtain basic configuration units for optimal configuration of water and soil resources in the river basin.
In the embodiment, hydrologic analysis is carried out by means of a tool (Hydrology) in ArcGIS software on the basis of a basin Digital Elevation Model (DEM) and a basin exit point (outlet), and sub-basins of a research area are drawn; performing intersection (intersector) treatment on the obtained sub-river basin and the administrative region to obtain a sub-river basin set administrative region basic unit; the areas of different land utilization types in the unit are counted by using the obtained drainage basin administrative region covering basic unit, and then a basic configuration unit (sub-drainage basin-administrative region-land utilization type) for optimal configuration of water and soil resources is obtained, as shown in fig. 2.
S2, determining an objective function of the optimal configuration of the watershed water and soil resources according to the basic configuration unit of the optimal configuration of the watershed water and soil resources, wherein the objective function comprises: the minimum total water shortage (quantity), the minimum total pollutant discharge (quality), the minimum ten-thousand-yuan water consumption (efficiency), the maximum net primary productivity (yield) and the maximum water and soil resource allocation damping coefficient of the basin are calculated as follows:
Figure BDA0003316243070000091
Figure BDA0003316243070000092
Figure BDA0003316243070000093
Figure BDA0003316243070000094
Figure BDA0003316243070000095
Figure BDA0003316243070000096
Figure BDA0003316243070000097
Hij(X)＝-Pij×LnPij
Figure BDA0003316243070000101
Figure BDA0003316243070000102
Figure BDA0003316243070000103
HHij(Y)＝-PPij×LnPPij
Figure BDA0003316243070000104
wherein, minF1(w,L)、minF2(w,L)、minF3(w,L)、minF4(L) and maxF5(w, L) respectively representing a minimum objective function of total water shortage of a drainage basin, a minimum objective function of total pollutant discharge, a minimum objective function of ten-thousand-yuan water consumption, a maximum objective function of net primary productivity and a maximum function of water and soil resource allocation damping coefficient, wijkIndicating the water demand of the jth land utilization type to the kth water source in the ith basic configuration unit, LijArea G representing the jth type of land use in the ith basic configuration unitikThe water supply amount of the kth water source in the ith basic configuration unit is shown, and beta represents the point source non-point source pollution judgmentCoefficient of difference, its value is 0 or 1, 1 represents point source pollution, 0 represents surface source pollution, eijRepresents the concentration of pollutants in the discharge wastewater of the jth soil utilization type in the ith basic configuration unit, pijA sewage discharge coefficient, A, representing the jth soil utilization type in the ith basic configuration unitijRepresents the pollutant load of the jth land utilization type in the ith basic configuration unit, giRepresenting GDP output value of ith basic configuration unit, n representing total number of configuration units, p representing total number of land utilization types, m representing total number of water sources, NPPijRepresenting the net primary productivity of the jth type of land utilization in the ith basic configuration unit,
Figure BDA0003316243070000105
shows the water and soil resource matching kini coefficient, x, of the kth water source in the flow fieldijRepresents the cumulative proportion of the j-th land utilization type balance degree in the i-th basic configuration unit, yijRepresents the cumulative proportion, x, of the water consumption balance of the jth land utilization type in the ith basic configuration unitijAnd yijForm a space Lorentzian curve, u represents the cumulative number of the configuration units, JijCumulative structural balance, H, representing the jth type of land utilization in the ith basic configuration cellijEntropy of structural information, H, representing the jth type of land use in the ith basic configuration unitijmaxMaximum structure information entropy representing all land use types in all configuration units, Ln representing natural logarithm, Hij(X) Structure information entropy, P, representing the jth type of land utilization in the ith basic configuration UnitijRepresents the total area proportion of the jth land utilization type in the ith basic configuration unit, JJujRepresents the cumulative water consumption balance of the jth land utilization type in the ith basic configuration unit, JJijWater consumption balance, HH, representing the jth land utilization type in the ith basic configuration unitijEntropy of water quantity information, HH, representing the jth type of land utilization in the ith basic configuration unitijmaxEntropy of maximum water usage information, HH, representing all land use types in all configured unitsij(Y)Entropy of water consumption information, PP, representing the jth type of land use in the ith basic configuration unitijThe water consumption of the jth land utilization type in the ith basic configuration unit is represented as the proportion of the total water consumption of the jth land utilization type in the whole flow field.
S3, determining constraint conditions of optimal configuration of water and soil resources in the drainage basin, wherein the constraint conditions comprise: the water balance constraint (basin total water consumption constraint and configuration unit water balance constraint), the pollutant discharge total amount constraint and the land resource balance constraint (basin land area total amount constraint and configuration unit land area balance constraint) are calculated according to the following formula:
Figure BDA0003316243070000111
Figure BDA0003316243070000112
Figure BDA0003316243070000121
Figure BDA0003316243070000122
Figure BDA0003316243070000123
wherein n represents the total number of configuration units, m represents the total number of land use types, GikIndicates the water supply available amount, W, of the kth water source of the ith basic configuration unitk(p) represents the total amount of available water resources of the kth water source when the incoming water frequency is p, SikIndicates the actual water supply amount, w, of the kth water source in the ith basic configuration unitijkThe water demand of the jth land utilization type to the kth water source in the ith basic configuration unit is represented, beta represents a point source surface source pollution discrimination coefficient, the value is 0 or 1, 1 is shown in the tablePoint source pollution, 0 surface source pollution, eijRepresents the concentration of pollutants in the j soil utilization type discharged wastewater in the ith basic configuration unit, pijSewage discharge coefficient, w, representing the jth type of land use in the ith basic configuration unitijkRepresenting the water demand of the jth land utilization type to the kth water source in the ith basic configuration unit, AijIndicating the pollutant load of the jth land utilization type of the ith basic configuration unit, LijAn area T representing the jth land utilization type in the ith basic configuration unitiRepresents the total amount of contamination consumed in the ith basic configuration unit, AiRepresents the total area, Ldown, of the ith basic configuration unitjIndicating the lower boundary, Lup, of the jth land utilization type red line on the ith basic configuration unitjAnd the upper boundary of the jth land utilization type in the drainage basin in the ith basic configuration unit is shown.
In this embodiment, first, a water resource constraint is determined. The total water supply of the Sihefu river region 2015 is 6.24 hundred million m3, wherein the surface water supply is 3.36 hundred million m32.45 hundred million meters of underground water supply30.4 hundred million m of other water source3. 5.25 hundred million meters of agricultural water according to water user statistics30.46 hundred million m of industrial water30.48 hundred million m of domestic water3Ecological water 0.4 hundred million m3. The water consumption of the water consumers is corresponding to the land utilization type (agricultural water consumption corresponds to cultivated land, domestic and industrial water consumption corresponds to construction sites, ecological water consumption corresponds to forest lands and grasslands), the water consumption of the water consumers is distributed to all the configuration units according to the land utilization type, the water consumption of the water areas and the unused land is not distributed, and the water consumption level of all the configuration units is further obtained. The water quantity is used as the current water quantity constraint of water resources.
Second, land resource constraints are determined. The total area of the Sihefu river basin is 2613.7km2Counting according to the land utilization type of 2014, and obtaining the cultivated land area 1726.3km in the drainage basin2Forest land area 119.5km2275.7km of grassland area295.7km of water area2And the area of the residential site is 384.8km211.6km of unused land area2. Counting various types of land in each configuration unit by means of ArcGIS toolAnd utilizing the area of the type, calculating the proportion of various land utilization types in the configuration unit to the areas of various land utilization types in the drainage basin, and taking the proportion as the basis for planning the land utilization area distribution of the year to be used as the constraint of land resources.
Finally, a contamination constraint is determined. The pollutants in the Siheyuen area are mainly polluted by a few sources (industry and life) and a surface source (farmland), and the load of the pollutants in the Siheyuen area is calculated as the constraint condition of planned annual pollutants according to the social and economic development and agricultural development level in the Jining city. The pollutant is mainly COD.
S4, calculating by using a multi-objective optimization algorithm according to the objective function and the constraint conditions to obtain an optimal allocation result of the water and soil in the drainage basin, and completing the optimal allocation of the resources of the water and soil in the drainage basin, wherein the implementation method comprises the following steps:
s401, unifying the objective functions into the minimum optimum;
in this embodiment, the objective functions include the maximum optimum and the minimum optimum, and the objective functions are unified into the minimum optimum to be solved when the solving calculation is performed.
S402, uniformly naming the optimized variables in a coding naming mode;
in the embodiment, optimization variables are determined, and the variables involved in the optimization process are divided into two types, namely the area of the land utilization type and the water consumption of different water sources by various land utilization types. Therefore, the variable is named uniformly by adopting a coding naming mode, so that the variable has universality. The encoding principle is as follows: each attribute takes two digits. The naming modes and the schematic diagrams of the variables of land use type and the variables of water resource type are shown in table 1 and table 2, respectively.
TABLE 1
Numbering Configuration Unit numbering Land use type numbering
101 1 01
6506 65 06
TABLE 2
Numbering Configuration Unit numbering Land use type numbering
10101 1 01
650602 65 06
Note: the configuration unit number refers to the total number of sub-flow areas formed by covering administrative areas, the land utilization types are six types (01-06), and the water source types are two types (01 is surface water and 02 is underground water).
S403, calculating the water consumption of each configuration unit and the area of each land utilization type by using a multi-objective optimization algorithm according to the constraint conditions;
s404, inputting the water consumption of each configuration unit and the area of each land utilization type into the FLUS model to obtain the optimal configuration result of the water and soil in the drainage basin, and completing the optimal configuration of the resources of the water and soil in the drainage basin.
In the embodiment, firstly, objective functions are unified, wherein the objective functions have maximum optimization and minimum optimization, and the objective functions are unified into the minimum optimization for solving when solving and calculating; secondly, determining optimization variables, wherein the variables involved in the optimization process are divided into two types, one is the area of the land utilization type, and the other is the water consumption of different water sources by various land utilization types. 822 optimization variables (274 land utilization variables and 548 water resource variables) are optimized, the number of the populations is set to be 100, the iteration times are set to be 100, the exchange probability is set to be 0.6, and the variation probability is set to be 0.05 by means of an NSGA-II algorithm. The calculation results are shown in fig. 3. And finally, performing space distribution, namely taking the various land utilization areas obtained by optimization solution as the input of the FLUS model, and finally obtaining the space distribution maps of various land utilizations. As shown in fig. 4.
Example 2
As shown in fig. 5, the present invention provides a system for optimizing and configuring resources of water and soil in a drainage basin, including:
the basic configuration unit determining module is used for determining a basic configuration unit for optimal configuration of water and soil resources of the drainage basin;
the objective function determining module is used for determining an objective function of the optimal allocation of the watershed water and soil resources according to the basic allocation unit of the optimal allocation of the watershed water and soil resources;
the constraint condition determining module is used for determining constraint conditions of optimal configuration of water and soil resources of the drainage basin;
and the optimal configuration result calculation module is used for calculating to obtain an optimal configuration result of the watershed water and soil by using a multi-objective optimization algorithm according to the objective function and the constraint condition, and completing the optimal configuration of the resources of the watershed water and soil.
The system for optimizing and configuring resources of water and soil in a drainage basin, which is provided in the embodiment shown in fig. 5, may implement the technical solutions shown in the above method embodiments, and the implementation principles and beneficial effects thereof are similar, and are not described herein again.
Example 3
The invention provides electronic equipment, which comprises a memory, a processor and a computer program, wherein the computer program is stored on the memory and runs on the processor, and the processor executes the program to realize the method for optimally configuring the resources of the water and the soil of the drainage basin in embodiment 1.
Example 4
The invention provides a computer-readable storage medium, which stores a computer program, wherein the computer program is executed by a processor to realize the method for optimally configuring the resources of the water and the soil in the drainage basin, which is described in the embodiment 1.
The invention provides a computer-readable storage medium, which stores a computer program, wherein the computer program is executed by a processor to realize the method for optimally configuring the resources of the water and the soil in the drainage basin, which is described in the embodiment 1.
The computer-readable storage medium described above may be implemented in any type or combination of volatile or non-volatile memory devices such as Static Random Access Memory (SRAM), erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disks, and may be any available medium that can be accessed by a general purpose or special purpose computer. A readable storage medium is coupled to the processor such that the processor can read information from, and write information to, the readable storage medium, which may also be a component of the processor, and the processor and the readable storage medium may reside in an Application Specific Integrated Circuit (ASIC), and the processor and the readable storage medium may also reside as discrete components in a resource optimization configuration system for a watershed water and soil.
Embodiments of the present application may be provided as a method, apparatus, or computer program product, and as such, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein. While the methods, apparatus (devices), and computer program products according to embodiments of the invention have been described with reference to flowchart illustrations and/or block diagrams, it is to be understood that each flowchart illustration and/or block diagram block or blocks, and combinations of flowchart illustrations and/or block diagrams, can be implemented by computer program instructions which are provided to a computer-readable memory of a general purpose computer, special purpose computer, embedded computer, or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart illustration of one or more flow diagrams and/or block diagrams block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks and/or flowchart block or blocks.

Claims (8)

1. A method for optimizing and configuring resources of water and soil in a drainage basin is characterized by comprising the following steps:
s1, determining a basic configuration unit for optimal configuration of the water and soil resources of the drainage basin;
s2, determining a target function of the optimal allocation of the watershed water and soil resources according to the basic allocation unit of the optimal allocation of the watershed water and soil resources;
s3, determining constraint conditions of optimal allocation of the watershed water and soil resources;
and S4, calculating by using a multi-objective optimization algorithm according to the objective function and the constraint condition to obtain an optimal allocation result of the water and soil of the drainage basin, and completing the optimal allocation of the resources of the water and soil of the drainage basin.
2. The method for resource optimal allocation of watershed water and soil according to claim 1, wherein the step S1 comprises the steps of:
s101, hydrologic analysis is carried out on the basis of a basin digital elevation model and a basin exit point, and sub-basins of a research area are drawn;
s102, performing intersection processing on the sub-river basin and the administrative region to obtain a sub-river basin set administrative region basic unit;
s103, the areas of different land types in the sub-flow field administrative region covering basic units are counted to obtain basic configuration units for optimal configuration of water and soil resources in the river basin.
3. The method for optimal allocation of resources in water and soil of a watershed as claimed in claim 1, wherein the expression of the objective function for optimal allocation of water and soil resources in the watershed in step S2 is as follows:
Figure FDA0003316243060000011
Figure FDA0003316243060000012
Figure FDA0003316243060000013
Figure FDA0003316243060000021
Figure FDA0003316243060000022
Figure FDA0003316243060000023
Figure FDA0003316243060000024
Hij(X)＝-Pij×LnPij
Figure FDA0003316243060000025
Figure FDA0003316243060000026
Figure FDA0003316243060000027
HHij(Y)＝-PPij×LnPPij
Figure FDA0003316243060000028
wherein, min F1(w,L)、min F2(w,L)、min F3(w,L)、min F4(L) and max F5(w, L) respectively representing a minimum objective function of total water shortage of a drainage basin, a minimum objective function of total pollutant discharge, a minimum objective function of ten-thousand-yuan water consumption, a maximum objective function of net primary productivity and a maximum function of water and soil resource allocation damping coefficient, wijkIndicating the water demand of the jth land utilization type to the kth water source in the ith basic configuration unit, LijArea G representing the jth type of land use in the ith basic configuration unitikThe water supply amount of the kth water source in the ith basic configuration unit is represented, beta represents a point source surface source pollution discrimination coefficient, and the value of the discrimination coefficient is 0 or 11 represents point source pollution, 0 represents area source pollution, eijRepresents the concentration of pollutants in the discharge wastewater of the jth soil utilization type in the ith basic configuration unit, pijA sewage discharge coefficient, A, representing the jth soil utilization type in the ith basic configuration unitijRepresents the pollutant load of the jth land utilization type in the ith basic configuration unit, giRepresenting GDP output value of ith basic configuration unit, n representing total number of configuration units, p representing total number of land utilization types, m representing total number of water sources, NPPijRepresenting the net primary productivity of the jth type of land utilization in the ith basic configuration unit,
Figure FDA0003316243060000031
shows the water and soil resource matching kini coefficient, x, of the kth water source in the flow fieldijRepresents the cumulative proportion of the j-th land utilization type balance degree in the i-th basic configuration unit, yijRepresents the cumulative proportion, x, of the water consumption balance of the jth land utilization type in the ith basic configuration unitijAnd yijForm a space Lorentzian curve, u represents the cumulative number of the configuration units, JijCumulative structural balance, H, representing the jth type of land utilization in the ith basic configuration cellijEntropy of structural information, H, representing the jth type of land use in the ith basic configuration unitijmaxMaximum structure information entropy representing all land use types in all configuration units, Ln representing natural logarithm, Hij(X) Structure information entropy, P, representing the jth type of land utilization in the ith basic configuration UnitijRepresents the total area proportion of the jth land utilization type in the ith basic configuration unit, JJujRepresents the cumulative water consumption balance of the jth land utilization type in the ith basic configuration unit, JJijWater consumption balance, HH, representing the jth land utilization type in the ith basic configuration unitijEntropy of water quantity information, HH, representing the jth type of land utilization in the ith basic configuration unitijmaxEntropy of maximum water usage information, HH, representing all land use types in all configured unitsij(Y) represents the i-th groupThe j-th land utilization type water consumption information entropy in the configuration unit, PPijThe water consumption of the jth land utilization type in the ith basic configuration unit is represented as the proportion of the total water consumption of the jth land utilization type in the whole flow field.
4. The method for optimal allocation of resources in water and soil of a drainage basin according to claim 1, wherein the expression of the constraint conditions in the step S3 is as follows:
Figure FDA0003316243060000041
Figure FDA0003316243060000042
Figure FDA0003316243060000043
Figure FDA0003316243060000044
Figure FDA0003316243060000045
wherein n represents the total number of configuration units, m represents the total number of land use types, GikIndicates the water supply available amount, W, of the kth water source of the ith basic configuration unitk(p) represents the total amount of available water resources of the kth water source when the incoming water frequency is p, SikIndicates the actual water supply amount, w, of the kth water source in the ith basic configuration unitijkThe water demand of the jth land utilization type to the kth water source in the ith basic configuration unit is represented, beta represents a point source and surface source pollution discrimination coefficient, the value is 0 or 1, 1 represents point source pollution, 0 represents surface source pollution, and e represents the surface source pollutionijIs shown asConcentration of pollutants in j' th land utilization type discharged wastewater, p, of i basic configuration unitsijSewage discharge coefficient, w, representing the jth type of land use in the ith basic configuration unitijkRepresenting the water demand of the jth land utilization type to the kth water source in the ith basic configuration unit, AijIndicating the pollutant load of the jth land utilization type of the ith basic configuration unit, LijAn area T representing the jth land utilization type in the ith basic configuration unitiRepresents the total amount of contamination consumed in the ith basic configuration unit, AiRepresents the total area, Ldown, of the ith basic configuration unitjIndicating the lower boundary, Lup, of the jth land utilization type red line on the ith basic configuration unitjAnd the upper boundary of the jth land utilization type in the drainage basin in the ith basic configuration unit is shown.
5. The method for resource optimal allocation of watershed water and soil according to claim 1, wherein the step S4 comprises the steps of:
s401, unifying the objective functions into the minimum optimum;
s402, uniformly naming the optimized variables in a coding naming mode;
s403, calculating the water consumption of each configuration unit and the area of each land utilization type by using a multi-objective optimization algorithm according to the constraint conditions;
s404, inputting the water consumption of each configuration unit and the area of each land utilization type into the FLUS model to obtain the optimal configuration result of the water and soil in the drainage basin, and completing the optimal configuration of the resources of the water and soil in the drainage basin.
6. A system for optimizing and configuring resources of water and soil in a drainage basin is characterized by comprising:
the basic configuration unit determining module is used for determining a basic configuration unit for optimal configuration of water and soil resources of the drainage basin;
the objective function determining module is used for determining an objective function of the optimal allocation of the watershed water and soil resources according to the basic allocation unit of the optimal allocation of the watershed water and soil resources;
the constraint condition determining module is used for determining constraint conditions of optimal configuration of water and soil resources of the drainage basin;
and the optimal configuration result calculation module is used for calculating to obtain an optimal configuration result of the watershed water and soil by using a multi-objective optimization algorithm according to the objective function and the constraint condition, and completing the optimal configuration of the resources of the watershed water and soil.
7. An electronic device comprising a memory, a processor and a computer program stored on the memory and running on the processor, wherein the processor executes the program to implement the method for resource-optimized allocation of watershed water and soil according to any one of claims 1 to 5.
8. A computer-readable storage medium storing a computer program, wherein the computer program is executed by a processor to implement the method for resource optimal allocation of watershed water and soil according to any one of claims 1 to 5.
CN202111231979.XA 2021-10-22 2021-10-22 Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium Pending CN113887073A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111231979.XA CN113887073A (en) 2021-10-22 2021-10-22 Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111231979.XA CN113887073A (en) 2021-10-22 2021-10-22 Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium

Publications (1)

Publication Number Publication Date
CN113887073A true CN113887073A (en) 2022-01-04

Family

ID=79004176

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111231979.XA Pending CN113887073A (en) 2021-10-22 2021-10-22 Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN113887073A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115081737A (en) * 2022-07-19 2022-09-20 中国水利水电科学研究院 Water-based land-fixing configuration method for regional full-aperture land type differential water

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115081737A (en) * 2022-07-19 2022-09-20 中国水利水电科学研究院 Water-based land-fixing configuration method for regional full-aperture land type differential water

Similar Documents

Publication Publication Date Title
Zarghami et al. System dynamics modeling for complex urban water systems: Application to the city of Tabriz, Iran
CN112765800B (en) Design method of distributed water resource configuration model
Simonovic et al. A new modeling approach for water resources policy analysis
CN105297827A (en) Water resource allocation method taking multi-user water demand and multi-source water supply into consideration
Wu et al. Prediction and analysis of water resources demand in Taiyuan City based on principal component analysis and BP neural network
Li et al. Multi-water resources optimal allocation based on multi-objective uncertain chance-constrained programming model
Song et al. Dynamic modeling application for simulating optimal policies on water conservation in Zhangjiakou City, China
Gao et al. Insights into water-energy cobenefits and trade-offs in water resource management
Liu et al. Comprehensive evaluation on water resources carrying capacity based on water-economy-ecology concept framework and EFAST-cloud model: A case study of Henan Province, China
Li et al. Dynamic and intelligent modeling methods for joint operation of a flood control system
CN113887073A (en) Method and system for optimal allocation of resources of water and soil in drainage basin, electronic equipment and storage medium
Li et al. Comprehensive evaluation and scenario simulation of water resources carrying capacity: a case study in Xiong’an New Area, China
Dong et al. Optimization and trade-off framework for coupled green-grey infrastructure considering environmental performance
Heydari Kushalshah et al. Hybrid modelling for urban water supply system management based on a bi-objective mathematical model and system dynamics: A case study in Guilan province.
CN116484647B (en) Distributed water resource allocation method and system for overall coordination
Chen et al. Towards Sustainable Cities: Studying Evaluation Index of Water Environment Carrying Capacity
Azizifard et al. Sustainable surface water management and wastewater treatment plant location: A case study of Urmia lake
Huang et al. Water-energy nexus analysis in an urban water supply system based on a water evaluation and planning model
Chen et al. Optimization of water resources utilization by GA–PSO in the Pinshuo open pit combined mining area, China
Zhao et al. Coupling research on spatio-temporal differentiation of water resources carrying capacity in mature mining cities
Zhang A study on urban water reuse management modeling
Lei et al. Study on ecological allocation of mine water in mining area based on long-term rainfall forecast
Qin et al. Risk assessment and configuration of water and land resources system network in the Huang-Huai-Hai watershed
Ge et al. Water resource optimization bi-level coupling model and carrying capacity of a typical plateau basin based on interval uncertainty stochastic programming
Chen et al. Optimization of river environmental management based on reinforcement learning algorithm: a case study of the Yellow River in China

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