CN116681262A - Water resource optimal configuration method based on remote sensing data - Google Patents
Water resource optimal configuration method based on remote sensing data Download PDFInfo
- Publication number
- CN116681262A CN116681262A CN202310934103.4A CN202310934103A CN116681262A CN 116681262 A CN116681262 A CN 116681262A CN 202310934103 A CN202310934103 A CN 202310934103A CN 116681262 A CN116681262 A CN 116681262A
- Authority
- CN
- China
- Prior art keywords
- water
- vegetation
- ecological
- data
- transpiration
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 212
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000005068 transpiration Effects 0.000 claims abstract description 48
- 238000011160 research Methods 0.000 claims abstract description 42
- 238000005457 optimization Methods 0.000 claims description 22
- 238000001704 evaporation Methods 0.000 claims description 9
- 230000008020 evaporation Effects 0.000 claims description 9
- 238000005192 partition Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000007812 deficiency Effects 0.000 claims description 4
- 239000002352 surface water Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 238000013468 resource allocation Methods 0.000 claims 1
- 239000000463 material Substances 0.000 abstract description 2
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- GVVPGTZRZFNKDS-JXMROGBWSA-N geranyl diphosphate Chemical compound CC(C)=CCC\C(C)=C\CO[P@](O)(=O)OP(O)(O)=O GVVPGTZRZFNKDS-JXMROGBWSA-N 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000002420 orchard Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 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/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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- 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
-
- 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/10—Services
- G06Q50/26—Government or public services
-
- 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Data Mining & Analysis (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Water Supply & Treatment (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a water resource optimal allocation method based on remote sensing data, and belongs to the technical field of water resource optimal allocation. Firstly, extracting ecological vegetation coverage conditions of a research area by using remote sensing image data; extracting vegetation transpiration data in a research area PML_V2, calculating a mean value of the vegetation transpiration data in the research area for many years, taking the obtained ecological vegetation spatial distribution condition as a frame, extracting a mean value of the vegetation transpiration data in a corresponding grid for many years, and obtaining a mean spatial distribution map of vegetation transpiration in the research area for many years; calculating the minimum ecological water requirement and the optimal ecological water requirement of the ecological vegetation in the research area according to the average spatial distribution map of the vegetation transpiration for years; constructing a water resource optimal configuration model: the method comprises the steps of solving an objective by an objective function and a condition constraint. The invention considers the space variability of ecological water demand more reasonably, not only can effectively reduce the investment of manpower and material resources, but also can realize the effective protection of ecological vegetation and improve the water resource management efficiency.
Description
Technical Field
The invention relates to a water resource optimal allocation method based on remote sensing data, and belongs to the technical field of water resource optimal allocation.
Background
The optimal allocation of water resources refers to scientific allocation of limited, different forms of water resources among users through engineering and non-engineering measures in a specific river basin or region according to effective, fair and sustainable principles. The optimal allocation of the water resources is the basis for realizing reasonable development and utilization of the water resources, and is the fundamental guarantee of sustainable utilization of the water resources. At present, a mathematical optimization model is an effective means for realizing water resource optimization configuration, and more mathematical algorithms are applied to the water resource optimization model, and uncertainty methods such as dynamic planning, interval planning, fuzzy mathematical planning and the like are widely applied. In the process of constructing the optimization model, some system parameters in the optimization model may be more complex, but most of parameters in the optimization model at present adopt a point-to-plane mode, and spatial variability of the parameters is not considered. As patent CN 112465206A discloses a water resource optimizing configuration and optimizing result evaluation feedback method, collecting monitoring data and regional population data; establishing a water resource bearing capacity index system; confirming index weight of the evaluation model, establishing a correlation function, and constructing a water resource bearing capacity evaluation model; the method comprises the steps of taking total water available for regional water resources as a core, taking maximized social and economic benefits and maximized unit system ecological benefits of each department as optimization targets, taking water available for each water source and water required by each department as constraint conditions, constructing a water resource optimal configuration model, and solving the model to obtain a regional water resource optimal configuration scheme. However, the patent does not consider the spatial variability of the ecological water demand, or use point values to calculate the model.
At present, with the development of remote sensing technology, space remote sensing technology is utilized to observe a large area with space-time continuity, continuous data can be acquired from remote sensing observation even if the number of observation sites is small, and the utilization of remote sensing technology and the development of a geographic information system enable more space information to be acquired.
Disclosure of Invention
The invention aims to provide a water resource optimal allocation method based on remote sensing data, aiming at areas with outstanding water resource supply and demand contradiction, the method can be combined with remote sensing basic monitoring data, and the scientific water resource optimal allocation can be realized by utilizing the actual evaporation of the ecological vegetation in an unmanned aerial vehicle or satellite remote sensing data observation research area and then calculating and acquiring the minimum and optimal ecological water demand in the research area.
The technical scheme adopted by the invention is as follows:
a water resource optimal configuration method based on remote sensing data comprises the following steps:
s1, extracting ecological vegetation coverage conditions of a research area by using remote sensing image data;
s2, extracting vegetation transpiration data in a research area PML_V2 by utilizing a PML_V2 land evapotranspiration and total primary productivity data set, calculating a multi-year vegetation transpiration data average value in the research area, and extracting a multi-year vegetation transpiration data average value in a corresponding grid by taking the ecological vegetation space distribution situation obtained in the step S1 as a frame to obtain a vegetation transpiration multi-year average spatial distribution map in the research area;
s3, calculating the minimum ecological water requirement and the optimal ecological water requirement of the ecological vegetation in the research area according to the average spatial distribution map of the vegetation transpiration for a plurality of years;
s4, constructing a water resource optimal configuration model: the method comprises the following steps of objective function and condition constraint, adopting 95% confidence interval numbers of annual water resource gazettes to represent the minimum and maximum water demand of agriculture, life and industry in a research area, and establishing a double-objective optimization model by combining the minimum ecological water demand and the optimum ecological water demand of each county in the research area estimated by remote sensing in the step S3:
1) The economic benefit is maximum
,
Wherein: f (F) 1 Is economic benefit, unit element; i is each county (i=1, 2, …) within the study area; j is water for each user, including industrial, agricultural, domestic and ecological water, unit m 3 ;B ij For i subareas j water use benefit coefficients, X ij The water supply amount of j water supply targets (users) in i subareas is represented by the number of intervals; wherein the lowest ecological water requirement and the most suitable ecological water requirement are the interval number (upper and lower limits) of the ecological water in the objective function;
2) Minimum water deficiency
,
Wherein: f (F) 2 For water shortage, the regional water shortage degree is used for expressing social benefit, W ijmax Maximum water demand for j water supply targets (users) in i subareas, X ij The water supply amount of j water supply targets (users) in i subareas is represented by the number of intervals;
converting the double-target model into a single-target model by adopting a weighted minimum deviation method to solve the double-target model;
constraints include a water availability limit, a water availability limit for each water user, and a non-negative constraint.
In the method, the ecological vegetation coverage in step S1 mainly extracts three ecological vegetation types including low coverage vegetation, high coverage vegetation and woodland, wherein the low coverage grassland mainly includes two types of attributes of 130 (grassland), 11 (herbaceous coverage), the high coverage grassland mainly includes four types of attributes of 120 (shrub), 121 (evergreen shrub), 122 (deciduous shrub), 12 (tree or shrub coverage (orchard)), and the woodland mainly includes 10 types of attributes of 51 (evergreen broadleaf forest), 52 (evergreen broadleaf forest), 61 (evergreen broadleaf forest), 62 (enclosed broadleaf forest), 71 (evergreen needle forest), 72 (evergreen needle forest), 81 (broadleaf needle forest), 82 (enclosed needle forest), 91 (evergreen mixed forest) and 92 (enclosed mixed forest). The attribute classification adopts the existing 30 m global fine earth surface coverage classification product in 2020 as a reference.
The pml_v2 land evaporation and total primary productivity dataset described in step S2 is downloaded from the national tibetan plateau data center. The calculation method of the average value of the annual vegetation transpiration data in the research area comprises the following steps: and extracting vegetation transpiration data in a study area PML_V2, overlapping Ec day by day in ArcGis software to obtain spatial distribution data of the vegetation transpiration of the study area for years, overlapping the spatial data of the vegetation transpiration of many years one by one, and then averaging to obtain an average value of the vegetation transpiration data of the study area for many years.
The minimum ecological water requirement and the optimum ecological water requirement in the step S3 are calculated according to the following formula:
,
,
i includes each county within the study area for each partition (i=1, 2, …); k is various ecological vegetation (k=1, 2, 3) including three ecological vegetation of low coverage vegetation and high coverage vegetation and woodland; ET (electric T) imin The lowest ecological water demand is i subareas; a is that ik Planting a quilt coverage area for the partition k of the i; ET (electric T) ikmin 10% quantiles of the transpiration quantity of k plants in the i partition are sequenced from small to large; ET (electric T) ikmax 90% quantile number for sorting transpiration amount of k plants in i subareas from small to large, ET imax Is the most suitable ecological water demand.
In the step S4, a weighted minimum deviation method is adopted to convert the double-target model into a single-target model for solving, and the specific formula is as follows:
,
wherein: f (F) 1 ,F 2 Respectively two objective functions in the optimization model; f (F) 1 max ,F 1 min ,F 2 max ,F 2 min Respectively the maximum value and the minimum value of the objective function; omega 1 ,ω 2 Respectively of two purposesThe weight of the scalar function, here taken as 0.5,0.5.
The constraint conditions in step S4 are specifically:
1) Limited water supply
,
Wherein: w (W) Public use Is commonly owned water, namely yellow water and river water; w (W) Independent and independent Is water unique to itself, namely surface water and underground water;
2) Water supply limit for each water user
,
Wherein: w (W) ijmin Minimum water demand for j water supply targets for i zones; w (W) ijmax The maximum water demand of j water supply targets is i subareas;
3) Non-negative constraints
。
The beneficial effects of the invention are as follows:
(1) And the remote sensing data is adopted to analyze the ecological water demand, so that the traditional artificial estimation method is replaced, and the accuracy of ecological water demand estimation is improved.
(2) The ecological water demand in the model is no longer in the form of point values, the minimum ecological water demand and the optimal ecological water demand of each subarea in the research area are no longer the same, and the space variability of the ecological water demand is more reasonably considered.
(3) Not only can the investment of manpower and material resources be effectively reduced, but also the effective protection of ecological vegetation can be realized, and the water resource management efficiency is improved.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
Further description will be provided below in connection with specific embodiments and exemplary examples.
A water resource optimal configuration method based on remote sensing data comprises the following steps:
s1, extracting ecological vegetation coverage conditions of a research area by using remote sensing image data:
the dynamic monitoring product is covered on the fine ground surface 30 m by using the space information innovation institute of China academy of sciences to obtain vegetation coverage conditions of a research area, and the data are downloaded from the space information innovation institute of China academy of sciences. The data are processed in ArcGIS software, and the ecological vegetation coverage conditions of a research area are extracted through tools such as mask extraction, attribute extraction and the like, so that three ecological vegetation types including low-coverage vegetation, high-coverage vegetation and woodland are mainly extracted, wherein the low-coverage grassland mainly comprises two types of attributes 130 and 11, the high-coverage grassland mainly comprises four types of attributes 120, 121, 122 and 12, and the woodland mainly comprises 10 types of attributes 51, 52, 61, 62, 71, 72, 81, 82, 91 and 92.
S2, extracting vegetation transpiration data in a research area PML_V2 by utilizing a PML_V2 land evapotranspiration and total primary productivity data set, calculating a multi-year vegetation transpiration data average value in the research area, and extracting a multi-year vegetation transpiration data average value in a corresponding grid by taking the ecological vegetation space distribution situation obtained in the step S1 as a frame to obtain a vegetation transpiration multi-year average spatial distribution map in the research area:
and selecting proper remote sensing data, and further acquiring the spatial distribution condition of vegetation actual evaporation. The water-carbon coupled land evaporation and total primary productivity dataset of the selected area pml_v2, pml_v2 (China), includes 5 elements of total primary productivity (gross primary product, GPP), vegetation transpiration (vegetation transpiration, ec), soil evaporation (Es), canopy interception evaporation (vaporization of intercepted rainfall, ei) and water and ice and snow evaporation (et_water). The data format is TIFF, and the space-time resolution is 1 day, 500 meters. Compared with the global version, the simulation precision of the PML_V2 (China) product in China is greatly improved, and the data are downloaded from the national Qinghai-Tibet plateau data center.
And extracting vegetation transpiration data in the PML_V2 of the research area, overlapping Ec day by day in ArcGis software to obtain vegetation transpiration spatial distribution data of the research area for years, and averaging the vegetation transpiration data of the research area for years. And (3) taking the space distribution condition of the ecological vegetation as a frame, extracting the average value of vegetation transpiration data in the corresponding grid, and obtaining an average space distribution map of vegetation transpiration in the research area for years.
S3, calculating the minimum ecological water requirement and the optimal ecological water requirement of the ecological vegetation in the research area according to the average spatial distribution map of the vegetation transpiration for years:
the ArcGis software is adopted to extract the transpiration values of all grids in three ecological vegetation distribution areas in each county in a research area respectively, the values are ordered from small to large in Excel respectively, and the minimum ecological water demand and the optimal ecological water demand of each county are calculated according to the following calculation formula:
,
,
i is each county within each sector (i=1, 2, …) pack study area; k is various ecological vegetation (k=1, 2, 3) including three ecological vegetation of low coverage vegetation and high coverage vegetation and woodland; ET (electric T) imin The lowest ecological water demand is i subareas; a is that ik Planting a quilt coverage area for the partition k of the i; ET (electric T) ikmin 10% quantiles of the transpiration quantity of k plants in the i partition are sequenced from small to large; ET (electric T) ikmax 90% quantile number for sorting transpiration amount of k plants in i subareas from small to large, ET imax Is the most suitable ecological water demand.
S4, constructing a water resource optimization configuration model, wherein the water resource optimization configuration model comprises an objective function and a condition constraint:
(1) Objective function
In order to accurately reflect the ecological water demand conditions of different counties in a research area and maximize economic benefit, a water resource optimization configuration based on remote sensing estimation of ecological water demand is established, the minimum and maximum water demand of agriculture, life and industry in the research area is represented by 95% confidence interval numbers of the past year water resource gazette, and a double-target optimization model is established by combining the minimum and maximum ecological water demand of each county in the research area estimated by remote sensing:
1) The economic benefit is maximum
,
Wherein: f (F) 1 Is economic benefit and primary; i is each county (i=1, 2, …) within the study area; j is water for each user, including industrial, agricultural, domestic and ecological water, m 3 ;B ij The water benefit coefficient of j water users in i subareas; x is X ij The water supply amount of j water supply targets (users) in i subareas is represented by the number of intervals; wherein the lowest ecological water requirement and the most suitable ecological water requirement are the interval number (upper and lower limits) of the ecological water in the objective function;
2) Minimum water deficiency
,
Wherein: f (F) 2 For water shortage, the regional water shortage degree is used for expressing social benefit, W ijmax The maximum water demand of j water supply targets in the i subareas.
(2) The target model is constrained by the following conditions:
1) Limited water supply
,
Wherein: w (W) Public use Is commonly owned water, namely yellow water and river water; w (W) Independent and independent Is water unique to itself, namely surface water and underground water;
2) Water supply limit for each water user
,
Wherein: w (W) ijmin Minimum water demand for j water supply targets for i zones; w (W) ijmax The most suitable for j water supply targets in i subareasLarge water demand;
3) Non-negative constraints
。
The model solving method comprises the following steps: and converting the double-target model into a single-target model by adopting a weighted minimum deviation method to solve, wherein the specific formula is as follows:
,
wherein: f (F) 1 ,F 2 Respectively two objective functions in the optimization model; f (F) 1 max ,F 1 min ,F 2 max ,F 2 min Respectively the maximum value and the minimum value of the objective function; omega 1 ,ω 2 The weights of the two objective functions are taken here as 0.5, respectively.
Examples: taking the optimal allocation of the city water resource of the Weifang as an example, the optimal allocation of the city water resource of the Weifang is carried out.
Firstly, the land classification data of the Weifang city is downloaded, and tools such as reclassification, mask extraction and the like are used in ArcGis to obtain the spatial distribution condition of three ecological vegetation including low coverage vegetation, high coverage vegetation and woodland in the Weifang city.
And downloading the water-carbon coupled land evaporation distribution and total primary productivity data set of the PML_V2 of the Weifang city by adopting ArcGis software, respectively overlapping Ec day by day in ArcGis according to the period of the year to obtain the vegetation transpiration space distribution condition of each year of the Weifang city, overlapping a plurality of years of data layers in ArcGis, and obtaining an average value to obtain the space distribution data of the average value of the vegetation transpiration years of the Weifang city.
Continuously processing in ArcGis, respectively extracting the transpiration values of three ecological vegetation in each county, respectively sequencing the values from small to large, and calculating the minimum ecological water demand and the optimal ecological water demand of each county, wherein the calculation formula is as follows:
,
,
i is each subarea (i=1, 2, …, 9) comprising 9 counties of prefecture, qingzhou city, all cities, shou city, anchu city, high density city, changyi city, lin_264; k is various ecological vegetation (k=1, 2, 3) including three ecological vegetation of low coverage vegetation and high coverage vegetation and woodland; ET (electric T) imin The lowest ecological water demand is i subareas; a is that ik Planting a quilt coverage area for the partition k of the i; ET (electric T) ikmin 10% quantiles of the transpiration quantity of k plants in the i partition are sequenced from small to large; ET (electric T) ikmax 90% quantile number for sorting transpiration amount of k plants in i subareas from small to large, ET imax Is the most suitable ecological water demand.
The calculated minimum ecological water requirement and the optimal ecological water requirement of each county are brought into an optimization model, and the optimization configuration model is as follows:
(1) Objective function
In order to accurately reflect the ecological water demand conditions of different counties in a research area and maximize economic benefit, a water resource optimization configuration based on remote sensing estimation of ecological water demand is established, the minimum and maximum water demand of agriculture, life and industry in the research area is represented by 95% confidence interval numbers of the past year water resource gazette, and a double-target optimization model is established by combining the minimum and maximum ecological water demand of each county in the research area estimated by remote sensing:
1) The economic benefit is maximum
(1),
Wherein: f (F) 1 Is economic benefit and primary; i is nine sectorial counties (i=1, 2, …) of the district of interest, qingzhou city, zhugu city, shou city, anchu city, high density city, changyi city, lin_264; j is water for each user, including industrial, agricultural, domestic and ecological water, m 3 ;B ij Partitioning j for iWater benefit coefficient, X of water user ij Water supply amounts for j kinds of water supply targets are i-zone.
2) Minimum water deficiency
(2) ,
Wherein: f (F) 2 For water shortage, the regional water shortage degree is used for expressing social benefit, W ijmax The maximum water demand of j water supply targets in the i subareas.
(2) Model targets are constrained by:
1) Limited water supply
(3),
Wherein: w (W) Public use Is commonly owned water, namely yellow water and river water; w (W) Independent and independent Is water unique to itself, namely surface water and underground water;
2) Water supply limit for each water user
(4),
Wherein: w (W) ijmin Minimum water demand for j water supply targets for i zones; w (W) ijmax The maximum water demand of j water supply targets is i subareas;
3) Non-negative constraints
(5)。
And converting the double-target model into a single-target model by adopting a weighted minimum deviation method to solve, wherein the specific formula is as follows:
(6),
wherein: f (F) 1 ,F 2 Respectively two objective functions in the optimization model; f (F) 1 max ,F 1 min ,F 2 max ,F 2 min Respectively the maximum value and the minimum value of the objective function; omega 1 ,ω 2 The weights of the two objective functions are taken here as 0.5, respectively.
The minimum ecological water demand and the maximum ecological water demand estimated by remote sensing are taken as the upper limit and the lower limit of the ecological water demand of the city of the Weifang in the model, the solution is carried out in Lingo and Excel, and formulas (1) (2) (3) (4) (5) (6) are written and operated in the Lingo in the form of codes, so that the results are shown in the following table:
(unit: ten thousand m) 3 )。
Claims (7)
1. The water resource optimizing configuration method based on remote sensing data is characterized by comprising the following steps:
s1, extracting ecological vegetation coverage conditions of a research area by using remote sensing image data;
s2, extracting vegetation transpiration data in a research area PML_V2 by utilizing a PML_V2 land evapotranspiration and total primary productivity data set, calculating a multi-year vegetation transpiration data average value in the research area, and extracting a multi-year vegetation transpiration data average value in a corresponding grid by taking the ecological vegetation space distribution situation obtained in the step S1 as a frame to obtain a vegetation transpiration multi-year average spatial distribution map in the research area;
s3, calculating the minimum ecological water requirement and the optimal ecological water requirement of the ecological vegetation in the research area according to the average spatial distribution map of the vegetation transpiration for a plurality of years;
s4, constructing a water resource optimal configuration model: the method comprises the following steps of objective function and condition constraint, adopting 95% confidence interval numbers of annual water resource gazettes to represent the minimum and maximum water demand of agriculture, life and industry in a research area, and establishing a double-objective optimization model by combining the minimum ecological water demand and the optimum ecological water demand of each county in the research area estimated by remote sensing in the step S3:
1) The economic benefit is maximum
,
Wherein: f (F) 1 Is economic benefit, unit element; i is each county (i=1, 2, …) within each study area; j is water for each user, including industrial, agricultural, domestic and ecological water, unit m 3 ;B ij For i subareas j water use benefit coefficients, X ij The water supply quantity of j water supply targets in i subareas is expressed by the number of intervals; wherein the minimum ecological water requirement and the most suitable ecological water requirement are the interval number of ecological water in the objective function;
2) Minimum water deficiency
,
Wherein: f (F) 2 For water shortage, the regional water shortage degree is used for expressing social benefit, W ijmax Maximum water demand for j water supply targets in i subareas, X ij The water supply quantity of j water supply targets in i subareas is expressed by the number of intervals;
converting the double-target model into a single-target model by adopting a weighted minimum deviation method to solve the double-target model;
constraints include a water availability limit, a water availability limit for each water user, and a non-negative constraint.
2. The water resource optimal allocation method based on remote sensing data according to claim 1, wherein the ecological vegetation coverage condition in the step S1 mainly extracts three ecological vegetation of low coverage vegetation, high coverage vegetation and woodland.
3. The method for optimizing water resource allocation based on remote sensing data according to claim 1, wherein the pml_v2 land evaporation and total primary productivity dataset of step S2 is downloaded from a national Qinghai-Tibet plateau data center.
4. The water resource optimal allocation method based on remote sensing data as claimed in claim 1, wherein the calculation method of the average value of the years vegetation transpiration data of the research area in step S2 is as follows: and extracting vegetation transpiration data in a study area PML-V2, overlapping Ec day by day in ArcGis software to obtain vegetation transpiration spatial distribution data of the study area for years, overlapping the spatial data for years one by one, and averaging to obtain an average value of vegetation transpiration data of the study area for years.
5. The water resource optimizing configuration method based on remote sensing data according to claim 1, wherein the minimum ecological water requirement and the optimum ecological water requirement in the step S3 are calculated as follows:
,
,
i is each county within each sector (i=1, 2, …) pack study area; k is various ecological vegetation (k=1, 2, 3) including three ecological vegetation of low coverage vegetation and high coverage vegetation and woodland; ET (electric T) imin The lowest ecological water demand is i subareas; a is that ik Planting a quilt coverage area for the partition k of the i; ET (electric T) ikmin 10% quantiles of the transpiration quantity of k plants in the i partition are sequenced from small to large; ET (electric T) ikmax 90% quantile number for sorting transpiration amount of k plants in i subareas from small to large, ET imax Is the most suitable ecological water demand.
6. The water resource optimization configuration method based on remote sensing data according to claim 1, wherein in step S4, a weighted minimum deviation method is adopted to convert a double-target model into a single-target model for solving, and the specific formula is as follows:
,
wherein: f (F) 1 ,F 2 Respectively two objective functions in the optimization model; f (F) 1 max ,F 1 min ,F 2 max ,F 2 min Respectively the maximum value and the minimum value of the objective function; omega 1 ,ω 2 The weights of the two objective functions are taken here as 0.5, respectively.
7. The water resource optimizing configuration method based on remote sensing data as claimed in claim 1, wherein the constraint condition in the step S4 is specifically:
1) Limited water supply
,
Wherein: w (W) Public use Is commonly owned water, namely yellow water and river water; w (W) Independent and independent Is water unique to itself, namely surface water and underground water;
2) Water supply limit for each water user
,
Wherein: w (W) ijmin Minimum water demand for j water supply targets for i zones; w (W) ijmax The maximum water demand of j water supply targets is i subareas;
3) Non-negative constraints
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310934103.4A CN116681262A (en) | 2023-07-28 | 2023-07-28 | Water resource optimal configuration method based on remote sensing data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310934103.4A CN116681262A (en) | 2023-07-28 | 2023-07-28 | Water resource optimal configuration method based on remote sensing data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116681262A true CN116681262A (en) | 2023-09-01 |
Family
ID=87787618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310934103.4A Pending CN116681262A (en) | 2023-07-28 | 2023-07-28 | Water resource optimal configuration method based on remote sensing data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116681262A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117236668A (en) * | 2023-11-15 | 2023-12-15 | 山东锋士信息技术有限公司 | Regional water resource allocation analysis decision method and system based on supply and consumption |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615285A (en) * | 2009-08-04 | 2009-12-30 | 北京师范大学 | A kind of computing method of ecological environment water demand threshold value |
CN107491868A (en) * | 2017-08-11 | 2017-12-19 | 南京信息工程大学 | A kind of forest land vegetation ecological needs the Quantizing Method of water |
CN113159560A (en) * | 2021-04-15 | 2021-07-23 | 中国科学院地理科学与资源研究所 | Optimized configuration method for psammophyte industry |
US20220092238A1 (en) * | 2020-07-13 | 2022-03-24 | Zhejiang University | Optimal allocation method of complex water resource system based on "zoning-grading" theory |
CN116341825A (en) * | 2023-02-08 | 2023-06-27 | 西安理工大学 | Water resource allocation method based on supply-demand balance and ecological benefit optimization |
-
2023
- 2023-07-28 CN CN202310934103.4A patent/CN116681262A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615285A (en) * | 2009-08-04 | 2009-12-30 | 北京师范大学 | A kind of computing method of ecological environment water demand threshold value |
CN107491868A (en) * | 2017-08-11 | 2017-12-19 | 南京信息工程大学 | A kind of forest land vegetation ecological needs the Quantizing Method of water |
US20220092238A1 (en) * | 2020-07-13 | 2022-03-24 | Zhejiang University | Optimal allocation method of complex water resource system based on "zoning-grading" theory |
CN113159560A (en) * | 2021-04-15 | 2021-07-23 | 中国科学院地理科学与资源研究所 | Optimized configuration method for psammophyte industry |
CN116341825A (en) * | 2023-02-08 | 2023-06-27 | 西安理工大学 | Water resource allocation method based on supply-demand balance and ecological benefit optimization |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117236668A (en) * | 2023-11-15 | 2023-12-15 | 山东锋士信息技术有限公司 | Regional water resource allocation analysis decision method and system based on supply and consumption |
CN117236668B (en) * | 2023-11-15 | 2024-03-08 | 山东锋士信息技术有限公司 | Regional water resource allocation analysis decision method and system based on supply and consumption |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101916337B (en) | Method for dynamically predicting potential productivity of paddy rice based on geographical information system | |
Worqlul et al. | Effect of climate change on land suitability for surface irrigation and irrigation potential of the shallow groundwater in Ghana | |
CN104486435B (en) | Low energy consumption ecological environmental monitoring node deployment method based on sensor network | |
Yan-Sui et al. | GIS-based assessment of land suitability for optimal allocation in the Qinling Mountains, China | |
CN109002621B (en) | Forest stand average height and breast diameter calculation method considering neighborhood and geographic difference | |
CN104794336B (en) | A kind of farmland stalk resources spatial distribution evaluation method | |
McConnell et al. | Precision conservation: A geospatial decision support tool for optimizing conservation and profitability in agricultural landscapes | |
Eweg et al. | Analysing degradation and rehabilitation for sustainable land management in the highlands of Ethiopia | |
CN102360458A (en) | Fuzzy evaluation method for selection of immigrant settlement areas | |
Fatih et al. | Multi criteria decision analysis to determine the suitability of agricultural crops for land consolidation areas | |
CN116681262A (en) | Water resource optimal configuration method based on remote sensing data | |
Handayani et al. | Estimation of built-up and green volume using geospatial techniques: A case study of Surabaya, Indonesia | |
CN113256101B (en) | Key driving force analysis method for water storage capacity change of urban and rural lakes | |
CN106803208A (en) | Agriculture block number is according to pair system and its matching method | |
Muleta et al. | GIS-based assessment of suitability area of rainwater harvesting in Daro Labu District, Oromia, Ethiopia | |
Jiang et al. | A geogrid-based framework of agricultural zoning for planning and management of water & land resources: A case study of northwest arid region of China | |
Yang et al. | Reconstruction of the spatial distribution of historical farmland in the Taiwan Province of China for 1659–1945 | |
Besbes et al. | Predictability of water resources with global climate models. Case of Northern Tunisia | |
Bandyopadhyay et al. | Strategies for backward-area development: A systems approach | |
Barchia et al. | Land Use Change Threat to Paddy Cultivation Sustainability on the Irrigated Rice Fields in Bengkulu Province, Indonesia | |
Lin et al. | Applying the grey assessment to the evaluation system of ecological green space on greening projects in Taiwan | |
Kai et al. | Evaluation on water source conservation capacity of West Liaohe River Basin based on InVEST model | |
Trigunasih et al. | Model of Numerical Spatial Classification for Sustainable Agriculture in Badung Regency and Denpasar City, Indonesia | |
CN116757357B (en) | Land ecological condition assessment method coupled with multisource remote sensing information | |
Mo et al. | Ecosystem Service Evaluation and Multi-Objective Management of Pinus massoniana Plantations in the North Tropics—the Case Study Area of Paiyangshan Forest Farm in Guangxi, China. Forests 2023, 14, 213 |
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 |