CN112464992A - Water quantity distribution method and system based on ecological gate characteristics - Google Patents

Water quantity distribution method and system based on ecological gate characteristics Download PDF

Info

Publication number
CN112464992A
CN112464992A CN202011223107.4A CN202011223107A CN112464992A CN 112464992 A CN112464992 A CN 112464992A CN 202011223107 A CN202011223107 A CN 202011223107A CN 112464992 A CN112464992 A CN 112464992A
Authority
CN
China
Prior art keywords
ecological
water
gate
gates
benefit
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.)
Granted
Application number
CN202011223107.4A
Other languages
Chinese (zh)
Other versions
CN112464992B (en
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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202011223107.4A priority Critical patent/CN112464992B/en
Publication of CN112464992A publication Critical patent/CN112464992A/en
Application granted granted Critical
Publication of CN112464992B publication Critical patent/CN112464992B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply

Abstract

The invention discloses a water distribution method and a system based on ecological gate characteristics, wherein the method comprises the following steps: calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of each ecological gate water distribution; making total water quantity constraint on a reference section of a river reach to be decided; NDVI, the ecological sensitivity index and the water shortage degree are used as classification indexes of the ecological gates, and the ecological gates with the water shortage degree larger than 0 and the water shortage degree smaller than 0 are classified by using a systematic clustering method; setting water quantity permission coefficients of all the ecological gates according to the types of the ecological gates and the water incoming conditions; and solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method. Therefore, the method can carry out overall distribution on the limited water resources according to the water demand degree, the benefit return and the like of a plurality of gates so as to obtain the maximized ecological benefit, can provide decision support for relevant management organizations, and has important significance for improving the utilization rate of the water resources in the water-deficient areas.

Description

Water quantity distribution method and system based on ecological gate characteristics
Technical Field
The invention relates to the technical field of efficient water resource utilization and environmental protection, in particular to a water quantity distribution method and system based on ecological gate characteristics.
Background
In wide arid and semi-arid regions in China, river water is a main source of industrial, agricultural and ecological water, but the amount of the water reserved as ecological water supplement in drainage basins after the water meets the requirements of the industrial and agricultural water is very limited. When ecological water is supplemented in part of the watershed, besides supplementing water by depending on the main stream leakage of the river channel, a gate specially used for ecological water supplement is built along the river channel to introduce the main stream river water into the channel, and underground water is supplemented by the on-way leakage of the channel and the infiltration of water flowing out from the outlet and then supplied to the riparian vegetation for utilization.
In the existing research scheme, when the water demand of a river channel is calculated, mostly, hydrologic stations are simply used as calculation sections for calculation, and the water demand calculated by the method is too macroscopic and difficult to guide fine scheduling; in addition, in the method disclosed in the chinese invention patent (CN201810729786.9), the river channel is segmented according to the geographical position of the gate dam, and then the ecological water demand of the river channel is calculated, but no specific method for water distribution is mentioned; in the method disclosed in the Chinese invention patent (CN201910488384.9), a method for allocating water for production, living and ecology in inland river basin is provided, but the method does not relate to a method for distributing water among gates for the same purpose.
In view of the current situation, how to arrange the diversion amount of each gate in a large number of gates is an urgent problem to be solved.
Disclosure of Invention
Aiming at the defects and the improvement requirements of the prior art, the invention provides a water quantity distribution method and a water quantity distribution system based on the characteristics of an ecological gate, and aims to solve the technical problem that the water quantity calculated by the existing water quantity distribution method is too macroscopic and is difficult to guide fine scheduling.
In order to achieve the aim, the invention provides a water quantity distribution method based on the characteristics of an ecological gate, which comprises the following steps: s1: calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of water distribution of each ecological gate based on the vegetation type and area of each ecological gate water diversion benefit area and the water supply amount of different vegetation; s2: making total water quantity constraint on a reference section of a river reach to be decided so that the total water quantity converted to the reference section and available for the ecological gates to quote does not exceed the total water demand of all ecological gates; s3: NDVI, the ecological sensitivity index and the water shortage degree are used as classification indexes of the ecological gates, and a systematic clustering method is used for classifying the ecological gates with the water shortage degree greater than 0 and the water shortage degree less than 0 respectively; s4: setting water quantity permission coefficients of the ecological gates according to the types of the ecological gates and the incoming water conditions, and accordingly determining the upper and lower boundaries of the diversion quantity of the ecological gates; s5: and under the condition that the constraint condition of the step S4 is met, solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method, wherein the ecological benefit is the sum of the product of the water diversion amount of each ecological gate and the unit ecological benefit.
Further, in step S1, the water demand of each ecological gate water diversion benefit area is as follows:
Figure BDA0002762744450000021
wherein d isiRepresenting the water demand (m) of the ith ecological gate water diversion benefit area3) (ii) a iu represents the vegetation type in the range of the water diversion benefit area of the ecological gate; v. ofiNumber of representative vegetation types; e.g. of the typeiuRepresentative vegetation area (hm)2);wiuRepresenting annual water supply amount (m) of vegetation3/hm2);
The unit ecological benefits of water distribution of each ecological gate are as follows:
Figure BDA0002762744450000022
wherein p isiRepresents the unit ecological benefit (Yuan/m) of the water distribution of the ith ecological gate3);piuRepresenting the ecological value per unit area (Yuan/hm)2);θiuRepresenting the value correction factor.
Further, in the step S2, the total water amount constraint is expressed as:
Figure BDA0002762744450000031
wherein Q isiThe diversion quantity (m) of the ith ecological gate3) (ii) a n is the total gate number; alpha is alphaiRepresenting the water delivery efficiency from the reference section to the section of the ecological gate;
Figure BDA0002762744450000032
the total water quantity which is converted to the reference section and can be quoted by the ecological gate is represented; c is the ratio of the total water quantity which can be quoted by the ecological gate and converted to the reference section to the total water quantity required by all the ecological gates converted to the reference section, and c is more than or equal to 0 and less than or equal to 1.
Further, in step S3, the water shortage degree is:
Figure BDA0002762744450000033
wherein s isiThe water shortage degree of the ith ecological gate, wsiFor the actual water diversion volume (m) of the last year3)。
Further, in step S4, for the ecological gate with the water shortage degree greater than 0, the upper and lower limits of the water quantity permission coefficient are both greater than c and decrease with the decrease of the value of c; for the ecological gate with the water shortage degree less than 0, the upper limit and the lower limit of the water quantity permission coefficient are both less than c and are reduced along with the reduction of the value of c.
In another aspect, the present invention provides a water distribution system based on the characteristics of an ecological gate, comprising:
the calculation module is used for calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of water distribution of each ecological gate based on the vegetation type and area of each ecological gate water diversion benefit area and the water supply amount of different vegetation;
the first constraint module is used for making total water quantity constraint on a reference section of a river reach to be decided so that the total water quantity converted to the reference section and available for the ecological gates to quote does not exceed the total water demand of all the ecological gates;
the gate classification module is used for classifying the ecological gates with the water shortage degrees larger than 0 and the water shortage degrees smaller than 0 by using a systematic clustering method by taking the NDVI, the ecological sensitivity index and the water shortage degrees as classification indexes of the ecological gates;
the second constraint module is used for setting water quantity permission coefficients of the ecological gates according to the types of the ecological gates and the incoming water conditions so as to determine the upper and lower boundaries of the diversion quantity of the ecological gates;
and the output module is used for solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method under the condition that the constraint conditions set by the second constraint module are met, wherein the ecological benefit is the sum of the product of the water diversion amount of each ecological gate and the unit ecological benefit.
Further, the water demand of each ecological gate water diversion benefit area is as follows:
Figure BDA0002762744450000041
wherein d isiRepresenting the water demand (m) of the ith ecological gate water diversion benefit area3) (ii) a iu represents the vegetation type in the range of the water diversion benefit area of the ecological gate; v. ofiNumber of representative vegetation types; e.g. of the typeiuRepresentative vegetation area (hm)2);wiuRepresenting annual water supply amount (m) of vegetation3/hm2);
The unit ecological benefits of water distribution of each ecological gate are as follows:
Figure BDA0002762744450000042
wherein p isiRepresents the unit ecological benefit (Yuan/m) of the water distribution of the ith ecological gate3);piuRepresenting the ecological value per unit area (Yuan/hm)2);θiuRepresenting the value correction factor.
Further, the total water volume constraint is expressed as:
Figure BDA0002762744450000043
wherein Q isiThe diversion quantity (m) of the ith ecological gate3) (ii) a n is the total gate number; alpha is alphaiRepresenting the water delivery efficiency from the reference section to the section of the ecological gate;
Figure BDA0002762744450000044
the total water quantity which is converted to the reference section and can be quoted by the ecological gate is represented; c is the ratio of the total water quantity which can be quoted by the ecological gate and converted to the reference section to the total water quantity required by all the ecological gates converted to the reference section, and c is more than or equal to 0 and less than or equal to 1.
Further, the water shortage degree is as follows:
Figure BDA0002762744450000045
wherein s isiThe water shortage degree of the ith ecological gate, wsiFor the actual water diversion volume (m) of the last year3)。
Furthermore, for the ecological gate with the water shortage degree greater than 0, the upper limit and the lower limit of the water quantity permission coefficient are both greater than c and are reduced along with the reduction of the value of c; for the ecological gate with the water shortage degree less than 0, the upper limit and the lower limit of the water quantity permission coefficient are both less than c and are reduced along with the reduction of the value of c.
Generally, by the above technical solution conceived by the present invention, the following beneficial effects can be obtained:
the method takes NDVI, the ecological sensitivity index and the water shortage degree as indexes, adopts a systematic clustering mode to classify the ecological gates, and sets water quantity permission coefficients under different incoming water conditions for the gates according to the characteristics of different types of gates; and the specific diversion water distribution of each gate when the annual ecological benefit is maximum is obtained by using a linear programming mode by taking the annual ecological benefit as an objective function and taking water quantity permission coefficients under different water supply conditions as constraint conditions. Therefore, the method can carry out overall distribution on the limited water resources according to the water demand degree, the benefit return and the like of a plurality of gates so as to obtain the maximized ecological benefit, can provide decision support for relevant management organizations, and has important significance for improving the utilization rate of the water resources in the water-deficient areas.
Drawings
FIG. 1 is a flow chart of a water distribution method based on the characteristics of an ecological gate provided by the invention;
fig. 2 is a schematic tree diagram of classification of an ecological gate with a water shortage degree less than 0 according to an embodiment of the present invention, in which an abscissa represents a gate number and an ordinate represents a distance between different types;
fig. 3 is a schematic tree diagram of classification of an ecological gate with a water shortage degree greater than 0 according to an embodiment of the present invention, in which an abscissa represents a gate number and an ordinate represents a distance between different types;
fig. 4 is a schematic diagram of a value range of a water volume permission coefficient under different water supply conditions according to an embodiment of the present invention, where an abscissa represents the water supply condition and an ordinate is the water volume permission coefficient.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. The features, operations, or characteristics described in the specification may be combined in any suitable manner to form various embodiments. Also, the various steps or actions in the method descriptions may be transposed or transposed in order, as will be apparent to one of ordinary skill in the art. Thus, the various sequences in the specification and drawings are for the purpose of describing certain embodiments only and are not intended to imply a required sequence unless otherwise indicated where such sequence must be followed. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
Referring to fig. 1, the present invention provides a water distribution method based on the characteristics of an ecological gate, comprising the following steps:
s1: calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of water distribution of each ecological gate based on the vegetation type and area of each ecological gate water diversion benefit area and the water supply amount of different vegetation;
specifically, the water demand of each ecological gate water diversion benefit area is as follows:
Figure BDA0002762744450000061
wherein d isiRepresenting the water demand (m) of the ith ecological gate water diversion benefit area3) (ii) a iu represents the vegetation type in the range of the water diversion benefit area of the ecological gate; v. ofiNumber of representative vegetation types; e.g. of the typeiuRepresentative vegetation area (hm)2);wiuRepresenting annual water supply amount (m) of vegetation3/hm2);
The unit ecological benefits of water distribution of each ecological gate are as follows:
Figure BDA0002762744450000062
wherein p isiRepresents the unit ecological benefit (Yuan/m) of the water distribution of the ith ecological gate3);piuRepresenting the ecological value per unit area (Yuan/hm)2);θiuRepresenting the value correction factor.
When classifying vegetation, the vegetation may be classified into a type such as a wooded land, a sparsely-cut land, a high-coverage grassland, and a low-coverage grassland, based on vegetation coverage.
S2: making total water quantity constraint on a reference section of a river reach to be decided so that the total water quantity converted to the reference section and available for the ecological gates to quote does not exceed the total water demand of all ecological gates;
specifically, the total water volume constraint is expressed as:
Figure BDA0002762744450000063
wherein Q isiThe diversion quantity (m) of the ith ecological gate3) (ii) a n is the total gate number; alpha is alphaiRepresenting the water delivery efficiency from the reference section to the section of the ecological gate;
Figure BDA0002762744450000071
the total water quantity which is converted to the reference section and can be quoted by the ecological gate is represented; c is the ratio of the total water quantity which can be quoted by the ecological gate and converted to the reference section to the total water quantity required by all the ecological gates converted to the reference section, and c is more than or equal to 0 and less than or equal to 1.
S3: NDVI, the ecological sensitivity index and the water shortage degree are used as classification indexes of the ecological gates, and a systematic clustering method is used for classifying the ecological gates with the water shortage degree greater than 0 and the water shortage degree less than 0 respectively;
wherein, the degree index of lack of water is used for weighing each ecological gate drainage benefit district soil water storage level before the decision of diversion of present phase, and the degree of lack of water expresses as:
Figure BDA0002762744450000072
wherein s isiFor the ith ecological gateDegree of water deficit, wsiFor the actual water diversion volume (m) of the last year3)。
Specifically, from three aspects of ecological water supply effect, efficiency, early water supply level and the like, a vegetation coverage index NDVI, an ecological sensitivity index and a water shortage degree are respectively adopted as classification indexes of the ecological gates, a systematic clustering method is used for classifying the ecological gates with the water shortage degree being greater than 0 and the water shortage degree being less than 0, Euclidean distances are adopted among samples, and inner square distances are adopted among classes.
If the ecological gates with the ecological water shortage degree of more than 0 and less than 0 are respectively divided into 3 types, the types of the ecological gates with the ecological water shortage degree of more than 0 are sequentially defined as L1, L2 and L3, as shown in FIG. 2, the serial numbers corresponding to three vertical lines intersected with the horizontal dotted line in the figure are the three divided types of gates, wherein the water shortage degree of the L1 type gates is higher overall, the NDVI is higher overall, and the ecological sensitivity is lower overall; the water shortage degree of the L2 gate is integrally lower, and the ecological sensitivity is integrally higher; the water shortage degree of the L3 gate is midstream, and the ecological sensitivity is much higher than that of other types of gates. The ecological gate type with the water shortage degree smaller than 0 is defined as S1, S2 and S3 in sequence, as shown in FIG. 3, the serial numbers corresponding to three vertical lines intersected with the horizontal dotted line in the figure are three divided gates, wherein the absolute value of the water shortage degree of the S1 gate is lower, and the ecological sensitivity is higher as a whole; the absolute value of the water shortage degree of the S2 gate is higher, and the ecological sensitivity is lower as a whole; the absolute value of the water shortage degree of the S3 gate is extremely high, the ecological sensitivity is high, and the NDVI is low.
S4: setting water quantity permission coefficients of the ecological gates according to the types of the ecological gates and the incoming water conditions, and accordingly determining the upper and lower boundaries of the diversion quantity of the ecological gates;
specifically, the water quantity constraint conditions are set according to the type of the ecological gate:
αidi≤Qi≤bidi
in the formula, ai、biThe lowest and highest water volume allowable coefficients are respectively. The setting rule of the water quantity permission coefficient is as follows: (1) for the gate with water shortage degree more than 0, the water quantity allowance coefficientBoth the upper and lower limits of (c) are greater than c and decrease with decreasing value of c; for the gate with the water shortage degree less than 0, the upper limit and the lower limit of the water quantity permission coefficient are both less than c and are reduced along with the reduction of the value of c; (2) the speed of the water quantity permission coefficients of L1, L2 and L3 decreasing with the decreasing c value is increased, and the speed of the water quantity permission coefficients of S1, S2 and S3 decreasing with the decreasing c value is also increased; (3) the water volume allowance factor is established to ensure that the water volume to be distributed is not higher than the total water volume and that the acceptable range of water volume is completely distributed without waste. And setting water quantity allowable coefficients of different types of gates under different c values according to the principle, and fitting a corresponding function curve, as shown in FIG. 4.
S5: and under the condition that the constraint condition of the step S4 is met, solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method, wherein the ecological benefit is the sum of the product of the water diversion amount of each ecological gate and the unit ecological benefit.
Specifically, the ecological benefit is expressed as:
Figure BDA0002762744450000081
in the formula, i represents an ecological gate; n represents the number of ecological gates; qiAs decision variables, represent annual water distribution (m)3);piRepresenting the unit ecological benefit (Yuan/m) of water distribution of the gate3). Dividing the c value within the range of 0-1 according to a certain step length, taking the maximum annual ecological benefit as an objective function under the condition of different c values, taking the water quantity permission coefficients of the gates of different types prepared in the step S4 as limiting conditions, solving the actual water diversion quantity of each gate when the ecological benefit is maximum under different water supply conditions by a linear programming method, and calculating by using related functions in calculation tools such as MATLAB and the like in a specific solving process. In the practical application process, the optimal water distribution scheme of each ecological gate can be inquired according to the calculation result only by calculating the c value according to the actual water coming quantity of the reference section of the decision river reach.
In another aspect, the present invention provides a water distribution system based on the characteristics of an ecological gate, comprising:
the calculation module is used for calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of water distribution of each ecological gate based on the vegetation type and area of each ecological gate water diversion benefit area and the water supply amount of different vegetation;
the first constraint module is used for making total water quantity constraint on a reference section of a river reach to be decided so that the total water quantity converted to the reference section and available for the ecological gates to quote does not exceed the total water demand of all the ecological gates;
the gate classification module is used for classifying the ecological gates with the water shortage degrees larger than 0 and the water shortage degrees smaller than 0 by using a systematic clustering method by taking the NDVI, the ecological sensitivity index and the water shortage degrees as classification indexes of the ecological gates;
the second constraint module is used for setting water quantity permission coefficients of the ecological gates according to the types of the ecological gates and the incoming water conditions so as to determine the upper and lower boundaries of the diversion quantity of the ecological gates;
and the output module is used for solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method under the condition that the constraint conditions set by the second constraint module are met, wherein the ecological benefit is the sum of the product of the water diversion amount of each ecological gate and the unit ecological benefit.
The division of the modules in the water distribution system based on the ecological gate characteristics is only used for illustration, and in other embodiments, the water distribution system based on the ecological gate characteristics can be divided into different modules as required to complete all or part of the functions of the device.
It will be readily understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, so that various changes, modifications and substitutions may be made without departing from the spirit and scope of the invention.

Claims (10)

1. A water quantity distribution method based on ecological gate characteristics is characterized by comprising the following steps:
s1: calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of water distribution of each ecological gate based on the vegetation type and area of each ecological gate water diversion benefit area and the water supply amount of different vegetation;
s2: making total water quantity constraint on a reference section of a river reach to be decided so that the total water quantity converted to the reference section and available for the ecological gates to quote does not exceed the total water demand of all ecological gates;
s3: NDVI, the ecological sensitivity index and the water shortage degree are used as classification indexes of the ecological gates, and a systematic clustering method is used for classifying the ecological gates with the water shortage degree greater than 0 and the water shortage degree less than 0 respectively;
s4: setting water quantity permission coefficients of the ecological gates according to the types of the ecological gates and the incoming water conditions, and accordingly determining the upper and lower boundaries of the diversion quantity of the ecological gates;
s5: and under the condition that the constraint condition of the step S4 is met, solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method, wherein the ecological benefit is the sum of the product of the water diversion amount of each ecological gate and the unit ecological benefit.
2. The ecological sluice feature-based water distribution method as recited in claim 1, wherein in step S1,
the water demand of each ecological gate water diversion benefit area is as follows:
Figure FDA0002762744440000011
wherein d isiRepresenting the water demand (m) of the ith ecological gate water diversion benefit area3) (ii) a iu represents the vegetation type in the range of the water diversion benefit area of the ecological gate; v. ofiNumber of representative vegetation types; e.g. of the typeiuRepresentative vegetation area(hm2);wiuRepresenting annual water supply amount (m) of vegetation3/hm2);
The unit ecological benefits of water distribution of each ecological gate are as follows:
Figure FDA0002762744440000021
wherein p isiRepresents the unit ecological benefit (Yuan/m) of the water distribution of the ith ecological gate3);piuRepresenting the ecological value per unit area (Yuan/hm)2);θiuRepresenting the value correction factor.
3. The ecological sluice feature-based water distribution method as claimed in claim 2, wherein in step S2, the total water quantity constraint is expressed as:
Figure FDA0002762744440000022
wherein Q isiThe diversion quantity (m) of the ith ecological gate3) (ii) a n is the total gate number; alpha is alphaiRepresenting the water delivery efficiency from the reference section to the section of the ecological gate;
Figure FDA0002762744440000023
the total water quantity which is converted to the reference section and can be quoted by the ecological gate is represented; c is the ratio of the total water quantity which can be quoted by the ecological gate and converted to the reference section to the total water quantity required by all the ecological gates converted to the reference section, and c is more than or equal to 0 and less than or equal to 1.
4. The ecological sluice feature-based water distribution method as recited in claim 3, wherein in step S3,
the water shortage degree is as follows:
Figure FDA0002762744440000024
wherein s isiThe water shortage degree of the ith ecological gate, wsiFor the actual water diversion volume (m) of the last year3)。
5. The ecological sluice feature-based water distribution method as recited in claim 4, wherein in step S4,
for the ecological gate with the water shortage degree greater than 0, the upper limit and the lower limit of the water quantity permission coefficient are both greater than c and are reduced along with the reduction of the value of c; for the ecological gate with the water shortage degree less than 0, the upper limit and the lower limit of the water quantity permission coefficient are both less than c and are reduced along with the reduction of the value of c.
6. A water distribution system based on ecological gate characteristics, comprising:
the calculation module is used for calculating the water demand of each ecological gate water diversion benefit area and the unit ecological benefit of water distribution of each ecological gate based on the vegetation type and area of each ecological gate water diversion benefit area and the water supply amount of different vegetation;
the first constraint module is used for making total water quantity constraint on a reference section of a river reach to be decided so that the total water quantity converted to the reference section and available for the ecological gates to quote does not exceed the total water demand of all the ecological gates;
the gate classification module is used for classifying the ecological gates with the water shortage degrees larger than 0 and the water shortage degrees smaller than 0 by using a systematic clustering method by taking the NDVI, the ecological sensitivity index and the water shortage degrees as classification indexes of the ecological gates;
the second constraint module is used for setting water quantity permission coefficients of the ecological gates according to the types of the ecological gates and the incoming water conditions so as to determine the upper and lower boundaries of the diversion quantity of the ecological gates;
and the output module is used for solving the actual water diversion amount of each ecological gate when the ecological benefit is maximum under different water supply conditions by using a linear programming method under the condition that the constraint conditions set by the second constraint module are met, wherein the ecological benefit is the sum of the product of the water diversion amount of each ecological gate and the unit ecological benefit.
7. The ecological sluice feature based water distribution system as recited in claim 6,
the water demand of each ecological gate water diversion benefit area is as follows:
Figure FDA0002762744440000031
wherein d isiRepresenting the water demand (m) of the ith ecological gate water diversion benefit area3) (ii) a iu represents the vegetation type in the range of the water diversion benefit area of the ecological gate; v. ofiNumber of representative vegetation types; e.g. of the typeiuRepresentative vegetation area (hm)2);wiuRepresenting annual water supply amount (m) of vegetation3/hm2);
The unit ecological benefits of water distribution of each ecological gate are as follows:
Figure FDA0002762744440000032
wherein p isiRepresents the unit ecological benefit (Yuan/m) of the water distribution of the ith ecological gate3);piuRepresenting the ecological value per unit area (Yuan/hm)2);θiuRepresenting the value correction factor.
8. The ecological gate feature based water distribution system of claim 7, wherein the total water volume constraint is expressed as:
Figure FDA0002762744440000041
wherein Q isiThe diversion quantity (m) of the ith ecological gate3) (ii) a n is the total gate number; alpha is alphaiRepresenting the water delivery efficiency from the reference section to the section of the ecological gate;
Figure FDA0002762744440000042
the total water quantity which is converted to the reference section and can be quoted by the ecological gate is represented; c is the ratio of the total water quantity which can be quoted by the ecological gate and converted to the reference section to the total water quantity required by all the ecological gates converted to the reference section, and c is more than or equal to 0 and less than or equal to 1.
9. The ecological sluice feature based water distribution system as recited in claim 8,
the water shortage degree is as follows:
Figure FDA0002762744440000043
wherein s isiThe water shortage degree of the ith ecological gate, wsiFor the actual water diversion volume (m) of the last year3)。
10. The ecological-gate-feature-based water distribution system of claim 9,
for the ecological gate with the water shortage degree greater than 0, the upper limit and the lower limit of the water quantity permission coefficient are both greater than c and are reduced along with the reduction of the value of c; for the ecological gate with the water shortage degree less than 0, the upper limit and the lower limit of the water quantity permission coefficient are both less than c and are reduced along with the reduction of the value of c.
CN202011223107.4A 2020-11-05 2020-11-05 Water quantity distribution method and system based on ecological gate characteristics Active CN112464992B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011223107.4A CN112464992B (en) 2020-11-05 2020-11-05 Water quantity distribution method and system based on ecological gate characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011223107.4A CN112464992B (en) 2020-11-05 2020-11-05 Water quantity distribution method and system based on ecological gate characteristics

Publications (2)

Publication Number Publication Date
CN112464992A true CN112464992A (en) 2021-03-09
CN112464992B CN112464992B (en) 2022-05-24

Family

ID=74826716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011223107.4A Active CN112464992B (en) 2020-11-05 2020-11-05 Water quantity distribution method and system based on ecological gate characteristics

Country Status (1)

Country Link
CN (1) CN112464992B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113269438A (en) * 2021-05-25 2021-08-17 华中科技大学 Ecological water right space-time configuration method and system

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102034001A (en) * 2010-12-16 2011-04-27 南京大学 Design method for distributed hydrological model by using grid as analog unit
US20120101861A1 (en) * 2010-10-25 2012-04-26 Lindores Robert J Wide-area agricultural monitoring and prediction
US20120109614A1 (en) * 2010-10-25 2012-05-03 Lindores Robert J Crop characteristic estimation
CN103930919A (en) * 2011-10-24 2014-07-16 天宝导航有限公司 Agricultural and soil management
US20150294154A1 (en) * 2014-04-15 2015-10-15 Open Range Consulting System and method for assessing riparian habitats
CN107944603A (en) * 2017-11-09 2018-04-20 中山大学 Water resource optimal allocation newsboy's method based on water total amount control
CN109101687A (en) * 2018-07-05 2018-12-28 沈阳建筑大学 A kind of Water Requirement calculation method based on gate dam subregion
CN109685256A (en) * 2018-12-13 2019-04-26 武汉大学 A kind of basin water dynamic optimization distribution method towards overall control
CN110245854A (en) * 2019-06-05 2019-09-17 中国科学院新疆生态与地理研究所 A kind of configuration method of inland river basin in arid zone production, life and ecological water
CN110852518A (en) * 2019-11-15 2020-02-28 中国水利水电科学研究院 Ecological water supplementing method for river estuary area of sediment-laden river

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120101861A1 (en) * 2010-10-25 2012-04-26 Lindores Robert J Wide-area agricultural monitoring and prediction
US20120109614A1 (en) * 2010-10-25 2012-05-03 Lindores Robert J Crop characteristic estimation
CN102034001A (en) * 2010-12-16 2011-04-27 南京大学 Design method for distributed hydrological model by using grid as analog unit
CN103930919A (en) * 2011-10-24 2014-07-16 天宝导航有限公司 Agricultural and soil management
US20150294154A1 (en) * 2014-04-15 2015-10-15 Open Range Consulting System and method for assessing riparian habitats
CN107944603A (en) * 2017-11-09 2018-04-20 中山大学 Water resource optimal allocation newsboy's method based on water total amount control
CN109101687A (en) * 2018-07-05 2018-12-28 沈阳建筑大学 A kind of Water Requirement calculation method based on gate dam subregion
CN109685256A (en) * 2018-12-13 2019-04-26 武汉大学 A kind of basin water dynamic optimization distribution method towards overall control
CN110245854A (en) * 2019-06-05 2019-09-17 中国科学院新疆生态与地理研究所 A kind of configuration method of inland river basin in arid zone production, life and ecological water
CN110852518A (en) * 2019-11-15 2020-02-28 中国水利水电科学研究院 Ecological water supplementing method for river estuary area of sediment-laden river

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A.PARENT等: "Control of Coliform Growth in Drinking Water Distribution Systems", 《WATER AND ENVIRONMENT JOURNAL》 *
秦劲军: "崇阳溪生态低影响开发模式及通航措施研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113269438A (en) * 2021-05-25 2021-08-17 华中科技大学 Ecological water right space-time configuration method and system
CN113269438B (en) * 2021-05-25 2024-03-19 华中科技大学 Ecological water right space-time configuration method and system

Also Published As

Publication number Publication date
CN112464992B (en) 2022-05-24

Similar Documents

Publication Publication Date Title
CN108764573B (en) Inland arid region-oriented multidimensional equilibrium configuration system for water resources
Wang et al. Modelling water quality and quantity with the influence of inter-basin water diversion projects and cascade reservoirs in the Middle-lower Hanjiang River
CN108108838B (en) Season regulation reservoir optimal scheduling method with high water utilization rate
CN106202978B (en) Fully-loaded stream multiple purpose reservoir flood regulation by level and stages optimization method and its system
CN110728035A (en) Pollutant total amount control method based on control of section water quality reaching standard
CN111476676B (en) Classification, classification and classification intelligent scheduling method for small river basin in town
CN105297827A (en) Water resource allocation method taking multi-user water demand and multi-source water supply into consideration
CN103276686A (en) Method of determining downstream design flood of cascade reservoirs
CN112101818B (en) Sponge city flood optimal scheduling method suitable for complex hydraulic connection
Diao et al. Risk analysis of flood control operation mode with forecast information based on a combination of risk sources
CN112464992B (en) Water quantity distribution method and system based on ecological gate characteristics
CN115271304A (en) Cross-basin water transfer engineering optimal scheduling method based on climate change
CN109377107A (en) A kind of optimization method of industry water multi-water resources system
CN109002946B (en) Water resource optimization scheduling method for 'two-reservoir-two-station' system for replenishing water in rivers and lakes
Hu et al. Total control-based unified allocation model for allowable basin water withdrawal and sewage discharge
Yuan et al. Water quantity and quality joint-operation modeling of dams and floodgates in Huai River Basin, China
Sun et al. A framework for deriving dispatching rules of integrated urban drainage systems
CN109919397A (en) Drain flooded fields capacity verification optimization method for river network in plain areas
CN105225077A (en) A kind of lake type basin water ecological function three grades of partition methods
CN106292298B (en) Series-parallel cascade pumping station water-carriage system space-time hierarchy optimization operation method and system
CN112396306B (en) Hierarchical water distribution method and system based on space-time water quantity balance of ecological gate
CN116780508A (en) Multi-uncertainty-based gradient hydropower-photovoltaic complementary system medium-long term interval optimal scheduling method
Zhou Integrated allocation model of water quantity and quality and its application in the Yellow River
CN103473716B (en) Heterogeneous module adopts the step Navigation & Hydropower Junction Engineering Optimal Scheduling of Standardization Communication
Li et al. Effect of pollution load reduction on water quality in rural lakes in the shallow hill water network area

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant