CN115907320A - Irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence - Google Patents
Irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence Download PDFInfo
- Publication number
- CN115907320A CN115907320A CN202210586171.1A CN202210586171A CN115907320A CN 115907320 A CN115907320 A CN 115907320A CN 202210586171 A CN202210586171 A CN 202210586171A CN 115907320 A CN115907320 A CN 115907320A
- Authority
- CN
- China
- Prior art keywords
- water
- water distribution
- channel
- distribution
- irrigation
- 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
Images
Abstract
The invention discloses an irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence, which is implemented as follows: step 1, collecting data required by the text; 2. distributing the wool water distribution amount in different water consumption departments in different areas; step 3, distributing water for each water department to each secondary channel; step 4, establishing a canal system optimization water distribution model, wherein decision variables are canal water distribution flow and water distribution starting time, and the objective is that the canal leakage loss is minimum and the water distribution duration is shortest; step 5, calculating water leakage loss of each water department in each region; step 6, calculating the net water distribution amount; step 7, configuring a model for water resources, wherein a decision variable is net water distribution amount, and a target function is water supply benefit and fairness; and 8, summarizing the indexes to obtain a water distribution scheme. The invention solves the problem that only water resource allocation is considered to influence a canal system water distribution list in the prior water distribution technology, and provides a water distribution method considering the water resource allocation and the canal system optimization water distribution mutual feedback influence.
Description
Technical Field
The invention belongs to the technical field of water resource optimization configuration of an irrigation district, and particularly relates to a water distribution method of the irrigation district, which considers the influence of water resource configuration and canal system water distribution bidirectional mutual feedback.
Background
The optimal allocation of water resources in the irrigation area refers to reasonably allocating limited water resources in each water-using department of each area within the range of the irrigation area so as to achieve the aims of maximizing the benefit of the irrigation area, considering water distribution fairness and the like. The channel optimization water distribution refers to that under the condition of certain channel overflow capacity, in order to meet the water distribution requirement of an irrigation area, a certain method and technology are adopted to optimize water distribution flow and water distribution time.
For a certain water distribution in an irrigation area, different allocation schemes of water resources among water using departments can cause different water distribution amounts of each channel, so that different influences are caused on channel system water distribution leakage loss, water distribution total time and the like. From the view of channel operation scheduling, different channel system water distribution schemes can generate different water distribution leakage losses, the net water distribution amount of water resources is influenced, and further the overall water distribution benefit and the water distribution fairness are influenced. The current research results are that an irrigation district water resource optimization configuration model and a canal system water distribution model are respectively researched, only the one-way influence of water resource configuration on canal system water distribution is considered, and the influence of canal system water distribution on water resource configuration is not considered.
Disclosure of Invention
The invention aims to provide an irrigation district water distribution method considering the influence of water resource allocation and canal system water distribution bidirectional mutual feedback, and solves the problem that the existing irrigation district water distribution technology only considers the unidirectional influence of water resource allocation on canal system water distribution.
The invention adopts the technical scheme that an irrigation area water distribution method considering the influence of water resource allocation and canal system water distribution bidirectional mutual feedback is implemented according to the following steps:
assuming that a irrigation area is formed by two levels of channels, wherein the irrigation area comprises m water supply areas, one primary channel in the irrigation area and z secondary channels for taking water from the primary channel, wherein the number of the secondary channels for supplying water to agriculture is z1, the number of the secondary channels for supplying water to industry is z2, the number of the secondary channels for supplying water to life is z3, and the number of the secondary channels for supplying water to ecology is z 4.
step 2, distributing the wool water distribution amount of different water departments in different areas;
step 3, the water distribution amount in the step 2 is redistributed to each secondary channel;
step 4, establishing a canal system optimization water distribution model, wherein decision variables are channel water distribution flow and water distribution starting time, and an objective function is the time when the channel leakage loss is minimum and the water distribution duration is shortest;
step 5, calculating the water leakage loss of each water consumption department in each region when water is distributed based on the channel leakage loss calculated in the step 4;
step 6, calculating the net water distribution amount, namely subtracting the loss water amount calculated in the step 5 from the gross water distribution amount distributed in the step 1 to obtain the net water distribution amount of each water use department of each area;
step 7, configuring a model for water resources, wherein a decision variable is the net water distribution amount calculated in the step 6, and a target function is the water supply benefit and fairness of the irrigation district;
and 8, summarizing the water distribution amount of each water use department in each area, the water distribution flow and the water distribution starting time of each channel, the water transmission and distribution loss, the water distribution duration, the fairness, the water distribution income of the irrigation area and other indexes in the steps to obtain a water distribution scheme.
The present invention is also characterized in that,
the data to be collected in step 1 are specifically as follows:
(1) Length of primary canal l u The unit is km;
(2) Design flow q of primary channel du In the unit of m 3 /s；
(3) Length of secondary canal l k The unit is km, wherein k is more than or equal to 1 and less than or equal to z, and z is the number of secondary channels;
(4) Design flow q of secondary channel dk Unit is m 3 S, wherein k is more than or equal to 1 and less than or equal to z;
(5) The total control area S of the irrigation area is hm 2 ；
(6) Control area S of each secondary channel k In units of hm 2’ Wherein k is more than or equal to 1 and less than or equal to z;
(7) The total water quantity W of the irrigation area in the water distribution time period is m 3 ；
(8) A water distribution wheel period T with the unit of d;
(9) Water demand W of different water using departments in different areas related to the irrigation area in the water distribution time period Need ij In the unit of m 3 I represents different zones, j represents different water use departments;
(10) Minimum water demand minW of different water using departments in different areas related to irrigation areas in water distribution time period Need j to In the unit of m 3 ；
(11) Water price p of different water consumption departments in different areas related to irrigation area in water distribution time period ij The unit is element;
(12) Soil permeability coefficient A of first-level channel bed u ；
(13) Leakage water yield reduction coefficient beta after seepage-proofing measures are taken in primary channel u ；
(14) Canal bed soil permeability index m of first-level channel u ；
(15) Second-level channel bed soil permeability coefficient A k ；
(16) Leakage water yield reduction coefficient beta after seepage prevention measures are taken for secondary channels k ；
(17) Canal bed soil permeability index m of second-level channel k 。
The step 2 is as follows:
adopting a getpy tool box to realize a genetic algorithm, generating an initial chromosome matrix, and specifically operating as follows: initialization: setting an evolution algebra counter h =1, namely initially setting an individual generation 1; setting the maximum evolution algebra H =50, namely, calculating the algorithm to the 50 th generation to terminate the calculation; generating a1 st generation population chromosome matrix P (h), h =1 by a crtpc function in a getpy tool box in combination with constraint conditions and decision variables to obtain the wool water distribution W Hair ij The crtpc function creates a population chromosome matrix;
(1) Decision variables:water purification W for distribution to different water departments in different regions Net ij ；
(2) Constraint conditions are as follows:
(1) the total water supply amount of the irrigation area is restricted, and the total water distribution amount meets the requirement of the water supply amount:
in the formula: m represents m areas in the irrigation area; n represents n water departments in the irrigation area;
(2) and (3) each water consumption department in each region is restricted by the minimum water distribution quantity, and the water distribution quantity of each water consumption department in each region is above the minimum water demand quantity:
W hair ij ≥minW Need ij (2)
The step 3 is specifically as follows:
distributing the agricultural irrigation water to each secondary channel, setting the secondary channel of the agricultural irrigation water in the area i as a, and distributing the rough water amount W to the channel Root of hairy Asia 1 Comprises the following steps:
in the formula: s represents the total control area of the area, with the unit of hm 2 ；W Wool i1 Represents the amount of agricultural water allocated to the area i in m 3 ；S a The control area of the secondary channel a is expressed in hm 2 ；
Distributing industrial water to each secondary channel, setting the secondary channel of industrial water in the area i as d, and distributing the crude water amount W to the channels Wool id2 ：
In the formula: w is a group of Need id2 The unit of the industrial water demand corresponding to the secondary channel d is m 3 ；W Need i2 Represents the total industrial water demand in the area and has the unit of m 3 ；
And calculating the water distribution amount when a certain channel is a reservoir, a living and ecological water supply channel in the same way.
Step 4, constructing an optimized water distribution model of the channel system as follows:
(1) Decision variables: water distribution flow q of secondary channel k And water distribution start time t sk ；
(2) An objective function:
establishing a model by taking the minimum channel leakage loss and the shortest water distribution duration in an irrigation area as targets in the water distribution process:
in the formula: s represents the leakage loss of the whole irrigation area and the unit is m 3 ；
S u Represents the leakage loss of the primary channel in m 3 ；
S k Represents the leakage loss of the secondary channel in m 3 ；
q ru The unit of the water distribution flow of the primary channel in the r-th time period is m 3 /s；
q k The unit of the flow of the water distribution of the secondary channel is m 3 /s；
t u The unit of the water distribution time of the primary channel is s;
t k representing the water distribution time of the secondary channel, wherein the unit is s;
r represents a water distribution period;
z represents the number of secondary channels;
k represents a secondary channel serial number;
objective function 2, irrigation area water distribution duration is shortest: the water distribution duration is the maximum value of the water distribution end time of the secondary channel:
minT＝max(t sk +t k ) (6)
in the formula: t is t sk Representing the starting time of the kth secondary channel, wherein the unit is h;
t k representing the water distribution duration of the kth secondary channel, wherein the unit is h;
(3) Constraint conditions are as follows:
(1) and secondary channel flow capacity constraint: the water distribution flow of any secondary channel is within 0.6-1.0 time of the design flow:
0.6q dk ≤q k ≤q dk (7)
in the formula: q. q of dk Represents the design flow of the secondary channel k, m 3 /s；
(2) And (3) cycle constraint: the water distribution starting time and the water distribution ending time of each secondary channel are within a period T:
0≤t sk ≤t sk +t k ≤T (8)
(3) water balance constraint 1: the product of the water distribution flow and the water distribution duration of the secondary channel is equal to the water distribution W of the channel Wool ijk Calculating by step 3:
W wool ijk ＝q k ·t k (9)
(4) And (3) water balance constraint 2: actual distribution flow q of primary channel at any time interval ru The sum of the water distribution flow of each secondary channel for water distribution in the time interval is equal to:
in the formula: x rk Showing the water distribution state of the channel k when t sk ≤r≤t sk +t k When, X rk =1, representing r period j channel water distribution; x rk =0 represents that the channel j of r period is not water;
(5) and primary channel overflow capacity constraint: at any time intervalWater distribution flow q of first-level channel ru Should be within 0.6-1.0 times of its design flow:
0.6q du ≤q ru ≤q du (11)
in the formula: q. q.s du Design flow m representing primary channel 3 /s；
The genetic algorithm solving model is realized by adopting a getpy tool box, and the solving process is as follows:
s1, adopting a getpy tool box to realize a genetic algorithm, and initializing: setting an evolution algebra counter g =1, namely, initially setting a generation 1 individual; setting the maximum evolutionary algebra G =100, namely, calculating the algorithm to the 100 th generation to terminate the calculation; generating a population chromosome matrix P (g) (g = 1) of the 1 st generation by a crtpc function in a getpy tool box, wherein the crtpc function combines the constraint conditions and the decision variables of the model built in the step, and the P (g) (g = 1) is a decision variable matrix meeting the constraint conditions, and the crtpc function creates a population chromosome matrix;
s2, individual evaluation: calculating the fitness of each individual in the group P (g) by combining a ranking function in a getpy toolbox and an objective function of the model built in the step;
s3, selection operation: applying a selection operator to the population P (g) through a tour function in the getpy toolkit;
s4, intersection operation: continuing to act the crossover operator on population P (g) in S3 through the xovdp function in the getpy toolkit;
s5, mutation operation: continuing to act on population P (g) in S4 with the mutation operator via the muturni function in the getpy toolbox; after the group P (x) is subjected to selection, intersection and mutation operation, g = g +1 to obtain a next generation group P (g);
s6, judging termination conditions: if G = G, outputting the individual with the maximum fitness obtained in the evolution process as an optimal solution to obtain the water distribution flow q of the secondary channel k And water distribution start time t sk And each target value, terminating the calculation; and if G is less than G, returning the group P (G) to the step S2, and calculating according to the steps in sequence.
The step 5 is specifically as follows:
setting and calculating the water loss caused by industrial water distribution in an area i, wherein the loss is equal to the sum of the water loss of the primary channel and the water loss of the secondary channel:
in the formula: w is a group of Loss of i2 Represents the water leakage loss caused by industrial water distribution of the area i and has the unit of m 3 ；
W To damage The unit of the leakage loss of a primary channel in the water distribution process is m 3 ；
d represents the d-th industrial water supply canal in the area 1;
p represents the number of industrial water canals in region 1;
v d represents the actual water supply flow rate of the channel d in m 3 /s；
v General assembly The sum of the water supply flow of all water supply channels is expressed in m 3 /s；
W Loss of i2d The leakage loss of the water delivery of the secondary channel d in the area i is expressed in the unit of m 3 ；
And the leakage loss including agriculture, industry, reservoir, life and ecological water distribution in different areas is obtained in the same way.
Step 6, the net water distribution W of each regional water department Net ij The calculation is as follows:
W net ij ＝W Hair ij -W Loss ij (13)
Step 7 is specifically as follows:
(1) Decision variables: the amount of purified water W distributed to different water consumption departments in different regions Net ij ；
(2) An objective function: establishing a model by taking the maximum benefit of water supply in an irrigation area and the fair distribution of water in different water departments in different areas as targets;
the objective function 1 is that the irrigation district water supply benefit is the maximum, and the irrigation district water supply benefit is the sum of the product of water cost and net water distribution amount of different water use departments in different regions:
in the formula: b represents the water supply benefit of the irrigation area, and the unit is Yuan;
i represents different areas related to the irrigation area;
m represents the number of regions;
j represents different water use departments;
n represents the number of water using departments;
p ij the unit of water charge is yuan/m 3 ；
w Net ij The unit of the purified water quantity is m, which represents the actual distributed purified water quantity of different water using departments in different areas 3 ；
And (2) an objective function 2 is that water distribution is fair for different water consumption departments in different areas of an irrigation area, and fairness is expressed by a standard deviation of a ratio value of actual water distribution and water demand:
in the formula: f represents the standard deviation of the ratio of the actual water distribution amount to the water demand amount;
ω ij the ratio of actual water distribution and water demand of different water consumption departments in different areas is represented;
(3) Constraint conditions are as follows:
(1) the water supply amount of the irrigation area is restricted, and the total water distribution amount meets the water supply amount requirement:
in the formula: m represents m areas in the irrigation area; n represents n water departments in the irrigation area;
(2) and (3) each water consumption department in each region is restricted by the minimum water distribution quantity, and the water distribution quantity of each water consumption department in each region is above the minimum water demand quantity:
W hair ij ≥minW Need ij (18)
And (4) solving the model established in the step (7) through a genetic algorithm on the basis of the net water distribution amount calculated in the step (6), wherein the model solving process is as follows:
b1, adopting a getpy tool box to realize a genetic algorithm, and initializing: step 1 is completed;
b2, individual evaluation: calculating the fitness of each individual in the population P (h) by combining a ranking function in a getpy toolbox and an objective function of the model built in the step;
b3, selection operation: applying a selection operator to the population P (h) through a tour function in the getpy toolkit;
b4, intersection operation: continuing to act the crossover operator on the population P (h) in B3 through the xovdp function in the getpy toolbox;
b5, mutation calculation: continuing to act on population P (h) in B4 with the mutation operator via the muturni function in the getpy toolbox; after the group P (x) is subjected to selection, intersection and mutation operation, h = h +1 to obtain a next generation group P (h);
b6, judging termination conditions: if H = H, the individual with the maximum fitness obtained in the evolution process is used as the optimal solution output to obtain the net water volume W of different water consumption departments in different regions Net ij And each target value, terminating the calculation; and if H is less than H, returning the population P (H) to the step B2, and calculating according to the steps in sequence.
The irrigation district water distribution method considering the bidirectional mutual feedback influence of water resource allocation and canal system water distribution has the advantages that the water resource allocation model is coupled with the canal system optimized water distribution model, after the water resource allocation model allocates the capillary water quantity, the canal system optimized water distribution model takes the capillary water quantity as the basis, the water leakage loss and other targets are obtained by adjusting the water distribution flow and the water distribution starting time, and the net water distribution quantity is obtained by subtracting the water leakage loss from the capillary water distribution quantity allocated by the water resource allocation model, so that the target of the water resource allocation model can be calculated more accurately, and a more reasonable water distribution scheme is obtained.
Drawings
FIG. 1 is a schematic diagram of an irrigation area in a water distribution method of the invention, considering the influence of water resource allocation and canal system water distribution bidirectional mutual feedback;
FIG. 2 is a flow chart of the present invention;
FIG. 3 is a diagram of the coupling model steps of the present invention;
FIG. 4 is a diagram illustrating a coupling model according to the present invention;
FIG. 5 is a distribution diagram of the irrigation canal system in the embodiment of the present invention.
Detailed Description
The invention is described in detail below with reference to the drawings and the detailed description.
The invention relates to an irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence, a flow chart is shown in figure 2, and the method is implemented according to the following steps:
assuming that a irrigation area is composed of two levels of channels, as shown in fig. 1, the irrigation area relates to m water supply areas, one primary channel in the irrigation area and z secondary channels take water from the primary channel, wherein z1 secondary channels for supplying water to agriculture, z2 secondary channels for supplying water to industry, z3 secondary channels for supplying water to life and z4 secondary channels for supplying water to ecology are provided;
the data to be collected in step 1 are specifically as follows:
(1) Length of primary canal l u The unit is km;
(2) Design flow q of primary channel du In the unit of m 3 /s；
(3) Length of secondary canal l k The unit is km, wherein k is more than or equal to 1 and less than or equal to z, and z is the number of secondary channels;
(4) Design flow q of secondary channel dk Unit is m 3 S, wherein k is more than or equal to 1 and less than or equal to z;
(5) The total control area S of the irrigation area is hm 2 ；
(6) Control area S of each secondary channel k In units of hm 2 ' wherein 1≤k≤z；
(7) The total water quantity W of the irrigation area in the water distribution time interval is m 3 ；
(8) A water distribution wheel period T with the unit of d;
(9) Water demand W of different water consumption departments in different areas related to the irrigation area in the water distribution time period Need ij Unit is m 3 I represents different zones, j represents different water use departments;
(10) Minimum water demand minW of different water consumption departments in different areas related to an irrigation area in a water distribution period Need ij In the unit of m 3 ；
(11) Water price p of different water consumption departments in different areas related to irrigation area in water distribution time period ij The unit is element; (12) Soil permeability coefficient A of first-level channel bed u ；
(13) Leakage water yield reduction coefficient beta after seepage-proofing measures are taken in primary channel u ；
(14) Bed soil permeability index m of first-level channel u ；
(15) Second-level channel bed soil permeability coefficient A k ；
(16) Leakage water yield reduction coefficient beta after seepage prevention measures are taken for secondary channels k ；
(17) Canal bed soil permeability index m of second-level channel k 。
Step 2, distributing the wool water distribution amount of different water consumption departments in different areas;
the step 2 is specifically as follows:
the genetic algorithm is realized by adopting a getpy tool box, only the first step of the genetic algorithm is involved in the step, and an initial chromosome matrix is generated, and the specific operation is as follows: initialization: setting an evolution algebra counter h =1, namely, initially setting a generation 1 individual; setting the maximum evolution algebra H =50, namely, stopping calculation when the algorithm is calculated to the 50 th generation; generating a1 st generation population chromosome matrix P (h), h =1 by a crtpc function in a getpy toolbox in combination with constraint conditions and decision variables (see below), and obtaining the wool water distribution W Hair ij The crtpc function creates a population chromosome matrix;
(1) Decision variables: is assigned toWater purification W for different water use departments in different regions Net ij ；
(2) Constraint conditions are as follows:
(1) the water supply amount of the irrigation area is restricted, and the total water distribution amount meets the requirement of the water supply amount:
in the formula: m represents m areas in the irrigation area; n represents n water departments in the irrigation area;
(2) and (3) each water consumption department in each region is restricted by the minimum water distribution quantity, and the water distribution quantity of each water consumption department in each region is above the minimum water demand quantity:
W hair ij ≥minW Need ij (2)
Step 3, distributing water to each secondary channel;
the step 3 is as follows:
the agricultural irrigation water is distributed to each secondary channel, and the secondary channel of the agricultural irrigation water in the area i is set as a (shown in figure 1), and the quantity W of the rough water distributed by the channel Root of hairy Asia 1 Comprises the following steps:
in the formula: s represents the total control area of the area, with the unit of hm 2 ；W Wool i1 Represents the amount of agricultural water allocated to the area i in m 3 ；S a The control area of the secondary channel a is expressed in hm 2 ；
Distributing industrial water to each secondary channel, and setting the industrial water secondary channel of the area i as d (shown in figure 1), wherein the channel distributes the crude water amount W Wool id2 ：
In the formula: w Need id2 The industrial water demand corresponding to the secondary channel d is representedBit is m 3 ；W Need i2 Represents the total industrial water demand in the area and has the unit of m 3 ；
And calculating the water distribution amount when a certain channel is a reservoir, a living and ecological water supply channel in the same way.
Step 4, establishing a canal system optimization water distribution model;
step 4, constructing an optimized water distribution model of the channel system as follows:
(1) Decision variables: water distribution flow q of secondary channel k And water distribution start time t sk ；
(2) An objective function:
establishing a model by taking the minimum channel leakage loss and the shortest water distribution duration in an irrigation area as targets in the water distribution process:
in the formula: s represents the leakage loss of the whole irrigation area in m 3 ；
S u Represents the leakage loss of the primary channel in m 3 ；
S k Represents the leakage loss of the secondary channel in m 3 ；
q ru The unit of the water distribution flow of the primary channel in the r-th period is m 3 /s；
q k The unit of the water distribution flow of the secondary channel is m 3 /s；
t u Representing the water distribution time of the primary channel, wherein the unit is s;
t k show the water distribution time of the second level channelThe bit is s;
r represents a water distribution period;
z represents the number of secondary channels;
k represents a secondary channel serial number;
objective function 2, irrigation area water distribution duration is shortest: the water distribution duration is the maximum value of the water distribution end time of the secondary channel:
minT＝max(t sk +t k ) (6)
in the formula: t is t sk Representing the starting time of the kth secondary channel, wherein the unit is h;
t k representing the water distribution duration of the kth secondary channel, wherein the unit is h;
(3) Constraint conditions are as follows:
(1) and secondary channel flow capacity constraint: the water distribution flow of any secondary channel is within 0.6-1.0 time of the design flow:
0.6q dk ≤q k ≤q dk (7)
in the formula: q. q.s dk Represents the design flow of the secondary channel k, m 3 /s；
(2) And (3) cycle restriction: the water distribution starting time and the water distribution ending time of each secondary channel are within a period T:
0≤t sk ≤t sk +t k ≤T (8)
(3) water balance constraint 1: the product of the water distribution flow and the water distribution duration of the secondary channel is equal to the water distribution W of the channel Wool ijk Calculating by step 3:
W wool ijk ＝q k ·t k (9)
(4) And (3) water balance constraint 2: actual distribution flow q of primary channel at any time interval ru The sum of the water distribution flow of each secondary channel for water distribution in the time period is equal to:
in the formula: x rk Represents the water distribution state of the channel k when t sk ≤r≤t sk +t k When the temperature of the water is higher than the set temperature,X rk =1, representing water distribution in the j channel in the r period; x rk =0 denotes that channel j is not water in r period;
(5) and primary channel overflowing capacity constraint: water distribution flow q of primary channel in any time period ru Should be within 0.6-1.0 times of its design flow:
0.6q du ≤q ru ≤q du (11)
in the formula: q. q.s du Design flow m representing primary channel 3 /s；
And (3) realizing a genetic algorithm solving model by adopting a getpy tool box, wherein the solving process is as follows (see a figure 4):
s1, realizing a genetic algorithm by adopting a getpy toolbox, and initializing: setting an evolution algebra counter g =1, namely initially setting an individual of the 1 st generation; setting the maximum evolutionary algebra G =100, namely, calculating the algorithm to the 100 th generation to terminate the calculation; generating a1 st generation population chromosome matrix P (g) (g = 1) by combining the constraint conditions and decision variables of the model built in the step through a crtpc function in a getpy tool box, wherein the P (g) (g = 1) is a decision variable matrix meeting the constraint conditions, and the crtpc function creates a population chromosome matrix;
s2, individual evaluation: calculating the fitness of each individual in the group P (g) by combining the objective function of the model built in the step through a ranking function in a getpy toolbox (the function can calculate the fitness based on the objective function value);
s3, selection operation: applying a selection operator to the population P (g) via the tour function (i.e., tournament selection operator) in the getpy toolbox;
s4, intersection operation: continuing to act on population P (g) in S3 by the xovdp function in the getpy toolkit (two-point intersection);
s5, mutation operation: continuing the mutation operator on population P (g) in S4 by the muturni function (uniform mutation) in the getpy toolbox; after the group P (x) is subjected to selection, intersection and mutation operation, g = g +1 to obtain a next generation group P (g);
s6, judging termination conditions: if G = G, the individual with the greatest fitness obtained in the evolution process is used as the candidate for the next generationObtaining the water distribution flow q of the secondary channel for the optimal solution output k And a water distribution start time t sk And each target value, terminating the calculation; and if G is less than G, returning the group P (G) to the step S2, and calculating according to the steps in sequence.
Step 5, calculating water leakage loss when water departments in various regions distribute water;
the step 5 is specifically as follows:
setting and calculating the water loss caused by industrial water distribution in an area i, wherein the loss is equal to the sum of the water loss of the primary channel and the water loss of the secondary channel:
in the formula: w is a group of Loss of i2 Represents the water leakage loss caused by industrial water distribution of the area i and has the unit of m 3 ；
W To damage The unit of the leakage loss of a primary channel in the water distribution process is m 3 ；
d represents the d-th industrial water supply canal in the area 1;
p represents the number of industrial water canals in region 1;
v d represents the actual water supply flow rate of the channel d in m 3 /s；
v General (1) The sum of the water supply flow of all water supply channels is expressed in m 3 /s；
W Loss of i2d Representing the loss of water delivery leakage of the secondary channel d in the area i in units of m 3 ；
And the leakage loss including agriculture, industry, reservoir, life and ecological water distribution in different areas is obtained in the same way.
Step 6, calculating the net water distribution amount, and subtracting the loss water amount calculated in the step 5 from the wool water distribution amount distributed in the step 1 to obtain the net water distribution amount of each regional water department;
step 6, the net water distribution W of each regional water department Net ij The calculation is as follows:
W net ij ＝W Hair ij -W Loss ij (13)
Step 7, a water resource configuration model, namely calculating two targets of water supply benefit and fairness of the irrigation district by taking the net water distribution amount calculated in the step 5 as a decision variable;
step 7 is specifically as follows:
(1) Decision variables: water purification W for distribution to different water departments in different regions Net ij ；
(2) An objective function: establishing a model by taking the maximum benefit of water supply in an irrigation area and the fair distribution of water in different water departments in different areas as targets;
the objective function 1 is that the irrigation district water supply benefit is the maximum, and the irrigation district water supply benefit is the sum of the product of water cost and net water distribution amount of different water use departments in different regions:
in the formula: b represents the water supply benefit of the irrigation area, and the unit is Yuan;
i represents different areas related to the irrigation area;
m represents the number of regions;
j represents different water use departments;
n represents the number of water using departments;
p ij the unit of water charge is yuan/m 3 ；
W Net ij The unit of the water quantity is m, which represents the actual distributed water quantity of different water departments in different areas 3 ；
And (2) an objective function 2 is that water distribution is fair for different water consumption departments in different areas of an irrigation area, and fairness is expressed by a standard deviation of a ratio value of actual water distribution and water demand:
in the formula: f represents the standard deviation of the ratio of the actual water distribution amount to the water demand amount;
ω ij the ratio of actual water distribution and water demand of different water consumption departments in different areas is represented;
(3) Constraint conditions are as follows:
(1) the water supply amount of the irrigation area is restricted, and the total water distribution amount meets the water supply amount requirement:
in the formula: m represents m areas in the irrigation area; n represents n water departments in the irrigation area;
(2) the minimum water distribution quantity of each water consumption department in each region is restricted, and the water distribution quantity of each water consumption department in each region is above the minimum water demand quantity:
W hair ij ≥minW Need ij (18)
And (3) solving the model established in the step (7) through a genetic algorithm on the basis of the net water distribution amount calculated in the step (6) (the solving process is shown in figure 4), wherein the model solving process is as follows:
b1, adopting a getpy tool box to realize a genetic algorithm, and initializing: step 1 is completed;
b2, individual evaluation: calculating the fitness of each individual in the group P (h) by combining the objective function of the model built in the step through a ranking function in a getpy toolbox (the function can calculate the fitness based on the objective function value);
b3, selection operation: applying a selection operator to the population P (h) through a tour function (i.e., a tournament selection operator) in the getpy toolbox;
b4, cross operation: continuing to act on the swarm P (h) in B3 by the xovdp function (two-point crossover) in the getpy toolbox;
b5, mutation operation: continuing the mutation operator on population P (h) in B4 by the muturni function (uniform mutation) in the getpy toolbox; after the group P (x) is subjected to selection, intersection and mutation operation, h = h +1 to obtain a next generation group P (h);
b6, judging termination conditions: if H = H, the individual with the maximum fitness obtained in the evolution process is used as the optimal solution output to obtain the net water volume W of different water consumption departments in different regions Net ij And each target value, terminating the calculation; if h<And H, returning the group P (H) to the step B2, and calculating according to the steps in sequence.
Step 8, summarizing the water distribution amount of each water use department in each area, the water distribution flow and water distribution starting time of each channel, water transmission and distribution loss, water distribution duration, fairness and irrigation district water distribution income indexes in the above steps to obtain a water distribution scheme, which is shown in the following table:
TABLE 1 Water distribution Scale for different water departments in different areas
Meter 2 channel water distribution information meter
Water target meter 3
Example application:
in the patent, model calculation is carried out by taking a dustpan plum irrigation area 2020 spring irrigation as an example. The dustpan plum irrigation area is located in coastal cities of Shandong province, and relates to three counties including Huimin county, yangxin county and non-county, and the irrigation area is designed to be 90 ten thousand mu. The water supply objects of the irrigation area are agricultural water of three counties and water for a reservoir, water is supplied by two main channels of one main channel, the total length of the main channel is 36.063km, the total length of one main channel is 46.384km, the total length of two main channels is 65.699km, 124 branch channels are arranged on the main channels in total, 67 branch channels are arranged in the counties of Huimin, and 1 branch channel is a reservoir water supply channel; 27 Yangxin county, 2 reservoir water supply channels; there are no 30 in the county, 2 are reservoir supply canals. The distribution diagram of the irrigation canal is shown in FIG. 5.
Step 1: data collection
(1) The total length of the main channel is 148.15km;
(2) Designing flow rate of a main channel: total main channel 75m 3 S, a main canal of 20m 3 Second main canal 40 m/s 3 /s；
(3) The branch canal length is shown in table 3;
TABLE 3 channel information Table
(4) The design flow of the branch channel is shown in a table 3;
(5) Total control area 78720hm of irrigation area 2 ；
(6) The control area of each lateral channel is shown in table 3;
(7) The water supply amount of the irrigation area in the water distribution time period is 27741.3 ten thousand meters 3 ；
(8) The water preparing period is 87 days;
(9) The water requirements of reservoirs in Huimin county, yangxin county and non-junkerer county and the agricultural water requirements of reservoirs in the third county are 610 ten thousand meters 3 1430 nm 3 2200 km 3 Wherein the water demand of two reservoirs in Yangxin county is 405 ten thousand meters 3 1025 ten thousand meters 3 The water demand of two reservoirs in the region without kerchien is 660m 3 、1540m 3 (ii) a The agricultural water requirement is 16337.88 ten thousand meters 3 4017.096 km 3 7984.08 km 3 ；
(10) The water requirement lower limit of a reservoir in the third county is 90% of the reported water quantity, and the agricultural water requirement lower limit is 70% of the reported water quantity;
(11) The water charges of reservoirs in three counties are all 0.31 yuan/m 3 The agricultural water cost is 0.04 yuan/m 3 ；
(12) Permeability coefficient of soil in canal bed of main canal and canal u ＝1.3；
(13) The leakage water yield reduction coefficient is respectively beta after the main canal adopts the anti-seepage measure u ＝0.1；
(14) Soil permeability index m of dry channel and canal bed u ＝0.35；
(15) Canal bed soil permeability coefficient A of branch canal channel u ＝1.3；
(16) The reduction coefficient of the leakage water amount is beta respectively after the branch channel adopts the anti-seepage measure u ＝1；
(17) Branch canal bed soil permeability index m u ＝0.35。
Step 2: and (3) distributing the wool water distribution amount in different water use departments in different areas: 27741.3m 3 The available water supply is distributed to two water using departments of three counties to obtain W Hair ij I =1,2,3 respectively represents Huimin county, yanxin county and Amelanchier county, and j =1,2 respectively represents agricultural water and reservoir water. The water distribution should satisfy the following constraints:
(1) the irrigated area available water supply constraint, which can be derived from equation (1):
(2) the minimum water distribution quantity constraint of each water consumption department in each region can be obtained by a formula (2):
W hair 11 ≥11436.52；W Wool 12 ≥549；
W Hair 21 ≥2811.97；W Hair 22 ≥1287：
W Hair 31 ≥5588.56：W Hair 32 ≥1980
And step 3: channel water distribution: take the branch channel 1 of Huimin county as an example, huimin branch channel 1 is the agricultural water supply channel, and the control area is 2333hm 2 . The total irrigation area of the dustpan plum irrigation area in Huimin county is 45383hm 2 The amount of the capillary water W distributed by the Huimin branch 1 can be obtained by the formula (3) Wool 111 Comprises the following steps:
taking the water supply channel 2 of Yangxin county as an example, the total water demand of the reservoir of Yangxin county is 1430 ten thousand meters 3 The water demand of the reservoir 2 is 405 ten thousand meters 3 The formula (4) can obtain the wool distributed by the 2 channels of the water inlet of the Yangxin county reservoirWater quantity W Hair 222 Comprises the following steps:
the water distribution amount of other agricultural and reservoir water supply channels can be calculated in the same way.
And 4, step 4: establishing an canal system optimization water distribution model:
(1) Decision variables: flow q of branch water distribution k And water distribution start time t sk 。
(2) An objective function:
the objective function 1 is that the channel leakage loss is minimum, and the water delivery loss of the trunk and branch channels in the water distribution process is obtained by the formula (5):
and (3) obtaining an objective function 2, wherein the total duration of water distribution of the irrigation area can be obtained by a formula (6):
minT＝max(t sk +t k )
(3) Constraint conditions are as follows:
(1) and (3) restricting the flow capacity of the branch channel: taking the branch canal 1 in Huimin county as an example, the design flow rate is 10m 3 And/s, the water distribution flow constraint can be obtained by the formula (7) as follows:
0.6×10≤q 1 ≤10
6≤q 1 ≤10
the same can be said for the flow constraints of other branches.
(2) And (3) cycle constraint: the irrigation period is 87 days, and can be obtained by the formula (8):
0≤t sk ≤t sk +t k ≤87
(3) water balance constraint 1: taking the branch canal 1 in Huimin county as an example, the product of water distribution flow and water distribution duration should be equal to W Wool 111 ，
From step 3 and equation (9):
q 1 ·t 1 ＝W wool 111 ＝0.05141W Hair 11
(4) And (3) water balance constraint 2: from equation (10):
(5) main canal flow capacity constraint: the design flow of the main canal is 75m 3 S, a design flow rate of 20m for a main canal 3 Design flow rate of 40m for the second main canal/s 3 And/s, as can be derived from equation (11):
45≤q r total ≤75
12≤q r is ≤20
24≤q e two ≤40
In the formula: q. q of r total of M represents the water distribution flow of the main canal in the period r 3 /s；
q r is M represents the water distribution flow of a main canal in the period r 3 /s；
q r two Shows the water distribution flow of the two main canals in the period r, m 3 /s。
And 5: and (4) calculating the water leakage loss when water consumption departments in various regions distribute water.
Taking the leakage loss of water caused by reservoir water distribution in Yangxin county as an example, the loss is equal to the sum of the water loss of the main channel and the water loss of the branch channel, and the formula (12) can be obtained:
in the formula: w is a group of 22 loss (c) M represents the water leakage loss caused by reservoir water distribution in Yangxin county 3 ；
v Library 2 M represents the actual water distribution flow of the reservoir water supply 2 channel 3 /s；
v Library 3 Shows the actual water distribution flow m of the reservoir water supply 3 channels 3 /s；
W Lossy 22 library 2 M represents the leakage water loss in the process of supplying water to the reservoir 2 channels 3 ；
W Lossy 22 library 3 M represents the leakage water loss in the process of 3 channel water distribution of reservoir water supply 3 。
And leakage lost water caused by agriculture in three counties and reservoir water distribution can be obtained in the same way.
Step 6: calculating the net water distribution quantity, and obtaining the formula (13):
W net ij ＝W Hair ij -W Loss ij
And 7: water resource allocation model
(1) Decision variables: water purification W for distribution to different water departments in different regions Net ij 。
(2) An objective function:
the objective function 1, which is the maximum benefit of water supply in the irrigation area, can be obtained from the formula (14):
maxB＝0.04·(W net 11 +W Medicine 21 +W Cleaning agent 31 )+0.31·(W Jing 12 +W Net 22 +W Medicine 32 )
Objective function 2, fair, given by equations (15) and (16):
and step 8: summarizing the indexes of the water distribution amount of each water use department in each area, the water distribution flow and the water distribution starting time of each channel, the water transmission and distribution loss, the water distribution duration, the fairness, the irrigation area water distribution income and the like in the steps, and obtaining a water distribution scheme as follows:
TABLE 4 Water distribution quantity of dual-purpose water intake of three counties
Meter 5 channel water distribution flow and water distribution time
The following table reflects the comparison of the optimized water distribution method with the traditional water distribution method, and it can be seen from table 6 that the four targets of the optimized water distribution method are superior to the traditional optimization method, which shows that the water distribution result of the patent is more efficient and reasonable.
Meter 6 optimized Water distribution method in contrast to traditional Water distribution methods
Claims (8)
1. The irrigation district water distribution method considering the water resource allocation and canal system water distribution bidirectional mutual feedback influence is characterized by being implemented according to the following steps:
assuming that one irrigation area is formed by two levels of channels, wherein the irrigation area relates to m water supply areas, one primary channel in the irrigation area and z secondary channels take water from the primary channel, wherein the number of the secondary channels for supplying water to agriculture is z1, the number of the secondary channels for supplying water to industry is z2, the number of the secondary channels for supplying water to life is z3, and the number of the secondary channels for supplying water to ecology is z 4;
step 1, collecting data in an irrigation area;
step 2, distributing the wool water distribution amount of different water consumption departments in different areas;
step 3, the water distribution amount in the step 2 is redistributed to each secondary channel;
step 4, establishing a canal system optimization water distribution model, wherein decision variables are channel water distribution flow and water distribution starting time, and an objective function is the time when the channel leakage loss is minimum and the water distribution duration is shortest;
step 5, calculating the water leakage loss of each water consumption department in each region when water is distributed based on the channel leakage loss calculated in the step 4;
step 6, calculating the net water distribution amount, and subtracting the loss water amount calculated in the step 5 from the wool water distribution amount distributed in the step 1 to obtain the net water distribution amount of each regional water department;
step 7, configuring a model for water resources, wherein a decision variable is the net water distribution amount calculated in the step 6, and a target function is the water supply benefit and fairness of the irrigation district;
and 8, summarizing the water distribution amount of each water use department in each area, the water distribution flow and the water distribution starting time of each channel, the water transmission and distribution loss, the water distribution duration, the fairness, the water distribution income of the irrigation area and other indexes in the steps to obtain a water distribution scheme.
2. The irrigation district water distribution method considering the influence of water resource allocation and canal system water distribution bidirectional mutual feedback as claimed in claim 1, wherein the data setting in step 1 is as follows:
(1) Length of primary canal l u The unit is km;
(2) Design flow q of primary channel du Unit is m 3 /s；
(3) Length of secondary canal l k The unit is km, wherein k is more than or equal to 1 and less than or equal to z, and z is the number of secondary channels;
(4) Design flow q of secondary channel dk Unit is m 3 S, wherein k is more than or equal to 1 and less than or equal to z;
(5) The total control area S of the irrigation area is hm 2 ；
(6) Control area S of each secondary channel k In units of hm 2 Wherein k is more than or equal to 1 and less than or equal to z;
(7) The total water quantity W of the irrigation area in the water distribution time interval is m 3 ；
(8) A water distribution wheel period T with the unit of d;
(9) Water demand W of different water using departments in different areas related to the irrigation area in the water distribution time period Need ij Unit is m 3 I represents different areas, j represents different water departments;
(10) Minimum water demand minW of different water consumption departments in different areas related to an irrigation area in a water distribution period Need ij In the unit of m 3 ；
(11) Water price p of different water consumption departments in different areas related to irrigation area in water distribution time period ij The unit is element;
(12) Soil permeability coefficient A of first-level channel bed u ；
(13) Leakage water yield reduction coefficient beta after seepage prevention measures are taken in primary channel u ；
(14) Canal bed soil of primary channelWater permeability index m of soil u ；
(15) Second-level channel bed soil permeability coefficient A k ；
(16) Leakage water yield reduction coefficient beta after seepage prevention measures are taken for secondary channels k ；
(17) Canal bed soil permeability index m of second-level channel k 。
3. The irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence as claimed in claim 2, wherein the step 2 is as follows:
adopting a getpy tool box to realize a genetic algorithm, generating an initial chromosome matrix, and specifically operating as follows: initialization: setting an evolution algebra counter h =1, namely, initially setting a generation 1 individual; setting the maximum evolution algebra H =50, namely, stopping calculation when the algorithm is calculated to the 50 th generation; generating a1 st generation population chromosome matrix P (h), h =1 by a crtpc function in a getpy tool box in combination with constraint conditions and decision variables to obtain the wool water distribution W Hair ij The crtpc function creates a population chromosome matrix;
(1) Decision variables: water purification W for distribution to different water departments in different regions Net ij ；
(2) Constraint conditions are as follows:
(1) the water supply amount of the irrigation area is restricted, and the total water distribution amount meets the water supply amount requirement:
in the formula: m represents m areas in the irrigation area; n represents n water departments in the irrigation area;
(2) the minimum water distribution quantity of each water consumption department in each region is restricted, and the water distribution quantity of each water consumption department in each region is above the minimum water demand quantity:
W hair ij ≥minW Need ij (2)。
4. The irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence as claimed in claim 3, wherein the step 3 is as follows:
distributing the agricultural irrigation water to each secondary channel, setting the secondary channel of the agricultural irrigation water in the area i as a, and distributing the rough water amount W to the channel Root of hairy Asia 1 Comprises the following steps:
in the formula: s represents the total control area of the area, and the unit is hm2; w Wool i1 Represents the amount of agricultural water allocated to the area i in m 3 ；S a The control area of the secondary channel a is expressed in hm 2 ；
Distributing industrial water to each secondary channel, setting the secondary channel of industrial water in the area i as d, and distributing the capillary water amount W to the channels Wool id2 ：
In the formula: w Need id2 The unit of the industrial water demand corresponding to the secondary channel d is m 3 ；W Need i2 Represents the total industrial water demand in the area and has the unit of m 3 ；
And calculating the water distribution amount when a certain channel is a reservoir, a living and ecological water supply channel in the same way.
5. The irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence according to claim 4, wherein the step 4 canal system optimized water distribution model is constructed as follows:
(1) Decision variables: water distribution flow q of secondary channel k And water distribution start time t sk ；
(2) An objective function:
establishing a model by taking the minimum channel leakage loss and the shortest water distribution duration in an irrigation area as targets in the water distribution process:
objective function 1, channel leakage loss is minimal: in the water distribution process, the flow of the primary channel at different time intervals is the sum of the flow of the secondary channel for water distribution at the time interval, the flow is constantly changed at any time interval, and the leakage loss of the primary channel needs to be calculated at different time intervals; the flow of each secondary channel is unchanged in the water distribution process, the leakage loss of the whole secondary channel is the sum of the losses of the secondary channels, and the leakage loss formula is as follows:
in the formula: s represents the leakage loss of the whole irrigation area and the unit is m 3 ；
S u Represents the leakage loss of the primary channel in m 3 ；
S k Represents the leakage loss of the secondary channel in m 3 ；
q ru The unit of the water distribution flow of the primary channel in the r-th period is m 3 /s；
q k The unit of the flow of the water distribution of the secondary channel is m 3 /s；
t u Representing the water distribution time of the primary channel, wherein the unit is s;
t k representing the water distribution time of the secondary channel, wherein the unit is s;
r represents a water distribution period;
z represents the number of secondary channels;
k represents a secondary channel serial number;
objective function 2, irrigation area water distribution duration is shortest: the water distribution duration is the maximum value of the water distribution end time of the secondary channel:
minT＝max(t sk +t k ) (6)
in the formula: t is t sk Representing the starting time of the kth secondary channel, wherein the unit is h;
t k representing the water distribution duration of the kth secondary channel, wherein the unit is h;
(3) Constraint conditions are as follows:
(1) and secondary channel overflowing capacity constraint: the water distribution flow of any secondary channel is within 0.6-1.0 time of the design flow:
0.6q dk ≤q k ≤q dk (7)
in the formula: q. q of dk Represents the design flow of the secondary channel k, m 3 /s；
(2) And (3) cycle restriction: the water distribution starting time and the water distribution ending time of each secondary channel are within a period T:
0≤t sk ≤t sk +t k ≤T (8)
(3) water balance constraint 1: the product of the water distribution flow and the water distribution duration of the secondary channel is equal to the water distribution W of the channel Wool ijk Calculating by step 3:
W wool ijk ＝q k ·t k (9)
(4) And (3) water balance constraint 2: actual distribution flow q of primary channel in any time interval ru The sum of the water distribution flow of each secondary channel for water distribution in the time period is equal to:
in the formula: x rk Showing the water distribution state of the channel k when t sk ≤r≤t sk +t k When, X rk =1, representing water distribution in the j channel in the r period; x rk =0 denotes that channel j is not water in r period;
(5) and primary channel overflow capacity constraint: water distribution flow q of primary channel at any time interval ru Should be within 0.6-1.0 times of its design flow:
0.6q du ≤q ru ≤q du (11)
in the formula: q. q.s du Design flow m representing primary channel 3 /s；
The genetic algorithm solving model is realized by adopting a getpy tool box, and the solving process is as follows:
s1, adopting a getpy tool box to realize a genetic algorithm, and initializing: setting an evolution algebra counter g =1, namely initially setting an individual of the 1 st generation; setting the maximum evolutionary algebra G =100, namely, stopping calculation when the algorithm is calculated to the 100 th generation; generating a population chromosome matrix P (g) (g = 1) of the 1 st generation by a crtpc function in a getpy tool box, wherein the crtpc function combines the constraint conditions and the decision variables of the model built in the step, and the P (g) (g = 1) is a decision variable matrix meeting the constraint conditions, and the crtpc function creates a population chromosome matrix;
s2, individual evaluation: calculating the fitness of each individual in the group P (g) by combining the target function of the model built in the step through a ranking function in a getpy toolbox;
s3, selection operation: applying a selection operator to the population P (g) through a tour function in the getpy toolkit;
s4, intersection operation: continuing to act the crossover operator on the swarm P (g) in S3 through the xovdp function in the getpy toolbox;
s5, mutation operation: continuing to act on population P (g) in S4 with the mutation operator via the muturni function in the getpy toolbox; after the group P (x) is subjected to selection, intersection and mutation operation, g = g +1 to obtain a next generation group P (g);
s6, judging termination conditions: if G = G, outputting the individual with the maximum fitness obtained in the evolution process as an optimal solution to obtain the water distribution flow q of the secondary channel k And a water distribution start time t sk And each target value, terminating the calculation; and if G is less than G, returning the population P (G) to the step S2, and calculating according to the steps in sequence.
6. The irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence as claimed in claim 5, wherein the step 5 is as follows:
and (3) setting and calculating the water loss caused by industrial water distribution in the area i, wherein the loss is equal to the sum of the water loss of the primary channel and the water loss of the secondary channel:
in the formula: w is a group of Loss of i2 Indicating area i industryWater leakage loss in m due to water distribution 3 ；
W To damage The unit of the leakage loss of a primary channel in the water distribution process is m 3 ；
d represents the d-th industrial water supply channel in the area 1;
p represents the number of industrial water supply channels in region 1;
v d represents the actual water supply flow rate of the channel d in m 3 /s；
v General assembly The sum of the water supply flow of all water supply channels is expressed in m 3 /s；
W Loss of i2d The leakage loss of the water delivery of the secondary channel d in the area i is expressed in the unit of m 3 ；
And leakage losses including agriculture, industry, reservoir, life and ecological water distribution in different areas are obtained in the same way.
7. The method of claim 6, wherein the net water distribution W of each regional water department in step 6 is determined by considering the mutual feedback effect of water resource allocation and canal system water distribution Net ij The calculation is as follows:
W net ij ＝W Hair ij -W Loss ij (13)。
8. The method for water distribution in irrigation areas considering the influence of water resource allocation and canal system water distribution mutual feedback according to claim 7, wherein the step 7 is as follows:
(1) Decision variables: water purification W for distribution to different water departments in different regions Net ij ；
(2) An objective function: establishing a model by taking the maximum benefit of water supply in an irrigation area and the water distribution fairness of different water consumption departments in different areas as targets;
the objective function 1 is that the irrigation district water supply benefit is the maximum, and the irrigation district water supply benefit is the sum of the products of water fees and net water distribution amounts of different water departments in different regions:
in the formula: b represents the water supply benefit of the irrigation area, and the unit is Yuan;
i represents different areas related to the irrigation area;
m represents the number of regions;
j represents different water use departments;
n represents the number of water using departments;
p ij the unit of water charge is yuan/m 3 ；
W Net ij The unit of the purified water quantity is m, which represents the actual distributed purified water quantity of different water using departments in different areas 3 ；
And (2) an objective function 2 is that water distribution is fair for different water consumption departments in different areas of an irrigation area, and fairness is expressed by a standard deviation of a ratio value of actual water distribution and water demand:
in the formula: f represents the standard deviation of the ratio of the actual water distribution amount to the water demand amount;
ω ij the ratio of actual water distribution and water demand of different water consumption departments in different areas is represented;
(3) Constraint conditions are as follows:
(1) the water supply amount of the irrigation area is restricted, and the total water distribution amount meets the requirement of the water supply amount:
in the formula: m represents m areas in the irrigation area; n represents n water departments in the irrigation area;
(2) and (3) each water consumption department in each region is restricted by the minimum water distribution quantity, and the water distribution quantity of each water consumption department in each region is above the minimum water demand quantity:
W hair ij ≥minW Need ij (18)
And (3) solving the model established in the step (7) by a genetic algorithm on the basis of the net water distribution calculated in the step (6), wherein the model solving process is as follows:
b1, adopting a getpy tool box to realize a genetic algorithm, and initializing: step 1 is completed;
b2, individual evaluation: calculating the fitness of each individual in the group P (h) by combining the target function of the model built in the step through a ranking function in the getpy toolbox;
b3, selection operation: applying a selection operator to the population P (h) through a tour function in the getpy toolkit;
b4, cross operation: continuing to act the crossover operator on the population P (h) in B3 through the xovdp function in the getpy toolbox;
b5, mutation calculation: continuing to act on population P (h) in B4 with the mutation operator via the muturni function in the getpy toolbox; after the group P (x) is subjected to selection, intersection and mutation operation, h = h +1 to obtain a next generation group P (h);
b6, judging termination conditions: if H = H, the individual with the maximum fitness obtained in the evolution process is used as the optimal solution output to obtain the net water volume W of different water consumption departments in different regions Net ij And each target value, terminating the calculation; and if H is less than H, returning the population P (H) to the step B2, and calculating according to the steps in sequence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210586171.1A CN115907320A (en) | 2022-05-27 | 2022-05-27 | Irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210586171.1A CN115907320A (en) | 2022-05-27 | 2022-05-27 | Irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115907320A true CN115907320A (en) | 2023-04-04 |
Family
ID=86469852
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210586171.1A Pending CN115907320A (en) | 2022-05-27 | 2022-05-27 | Irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115907320A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116562133A (en) * | 2023-04-26 | 2023-08-08 | 水利部交通运输部国家能源局南京水利科学研究院 | Canal system multi-target water distribution real-time optimization method for self-flow irrigation areas |
| CN116796997A (en) * | 2023-08-01 | 2023-09-22 | 东北农业大学 | Three-level canal system water distribution optimization method based on third-generation non-dominant sorting genetic algorithm |
-
2022
- 2022-05-27 CN CN202210586171.1A patent/CN115907320A/en active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116562133A (en) * | 2023-04-26 | 2023-08-08 | 水利部交通运输部国家能源局南京水利科学研究院 | Canal system multi-target water distribution real-time optimization method for self-flow irrigation areas |
| CN116562133B (en) * | 2023-04-26 | 2024-02-23 | 水利部交通运输部国家能源局南京水利科学研究院 | Canal system multi-target water distribution real-time optimization method for self-flow irrigation areas |
| CN116796997A (en) * | 2023-08-01 | 2023-09-22 | 东北农业大学 | Three-level canal system water distribution optimization method based on third-generation non-dominant sorting genetic algorithm |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115907320A (en) | Irrigation district water distribution method considering water resource allocation and canal system water distribution bidirectional mutual feedback influence | |
| CN109685256B (en) | Basin water quantity dynamic optimization distribution method facing total amount control | |
| CN109658287B (en) | Basin water quantity scheduling method based on water resource space-time uniform distribution | |
| CN110689283B (en) | River system and canal system staggered water supply simulation method and device based on underground water model | |
| Yang et al. | Optimization of the irrigation water resources for Shijin irrigation district in north China | |
| CN114202167B (en) | Multi-water-source irrigation area irrigation water consumption metering statistical method based on water circulation simulation | |
| CN110889553B (en) | Real-time water distribution method based on optimization strategy | |
| CN101356882A (en) | Determination method of irrigation district water delivery project | |
| CN105096004A (en) | Real-time scheduling method for reservoir group water supply and transfer system | |
| CN106327065A (en) | Water resource optimization configuration method for single pumping station - single reservoir system for direct canal supplement under full irrigation condition | |
| CN108805449A (en) | Cooperative game method towards integrated energy system cost sharing and distribution of income | |
| CN105913146A (en) | Southern moist region water resource optimal allocation system | |
| CN115796381B (en) | Actual runoff forecasting method based on improved Xinanjiang model | |
| CN113657811A (en) | Well and canal combined irrigation area water saving potential analysis method based on zero excess mining of underground water | |
| CN106651011A (en) | Particle swarm algorithm-based canal system optimization water distribution method | |
| CN112005679B (en) | Water and fertilizer integrated distribution method for precision agriculture | |
| CN112214885B (en) | Irrigation area crop irrigation water quantity optimal distribution method under insufficient irrigation condition | |
| Fan et al. | A Model Coupling Water Resource Allocation and Canal Optimization for Water Distribution | |
| CN110310022B (en) | Power plant monthly power generation plan decomposition method based on electric power spot market price | |
| CN113887882A (en) | Reservoir optimal scheduling method considering virtual water based on improved whale optimization algorithm | |
| CN104594281B (en) | Running dispatching diagram drawing method capable of enabling comprehensive benefits of reservoir to be optimal | |
| CN108960616A (en) | Strong anthropic zone's Water resources security diagnosis space atlas calculation based on county domain | |
| Zhao et al. | Rationality Analysis of the Heihe River'97'Water Diversion Scheme under Changing Environment | |
| Laufer et al. | Optimal weekly releases from a seasonal reservoir: 1. Deterministic future | |
| CN114118809A (en) | Method and device for distributing carbon emission, electronic device, and storage medium |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Fan Yu Inventor after: Gao Zhanyi Inventor after: Chen Haorui Inventor after: Liu Jing Inventor after: Gao Jian Inventor before: Fan Yu Inventor before: Gao Zhanyi Inventor before: Chen Haorui Inventor before: Liu Jing Inventor before: Gao Jian |
|
| CB03 | Change of inventor or designer information |