CN113449993B - Urban water source water supply scheduling method - Google Patents
Urban water source water supply scheduling method Download PDFInfo
- Publication number
- CN113449993B CN113449993B CN202110728549.2A CN202110728549A CN113449993B CN 113449993 B CN113449993 B CN 113449993B CN 202110728549 A CN202110728549 A CN 202110728549A CN 113449993 B CN113449993 B CN 113449993B
- Authority
- CN
- China
- Prior art keywords
- water
- node
- water supply
- period
- supply
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 454
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 48
- 239000011159 matrix material Substances 0.000 claims abstract description 48
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 230000005540 biological transmission Effects 0.000 claims description 28
- 238000005086 pumping Methods 0.000 claims description 17
- 238000012546 transfer Methods 0.000 claims description 10
- 239000008400 supply water Substances 0.000 claims description 9
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 4
- 238000001704 evaporation Methods 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 3
- 238000001595 flow curve Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 2
- 241000196324 Embryophyta Species 0.000 description 49
- 238000010586 diagram Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 239000004576 sand Substances 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000012821 model calculation Methods 0.000 description 4
- 244000261422 Lysimachia clethroides Species 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000011144 upstream manufacturing Methods 0.000 description 3
- 241000123069 Ocyurus chrysurus Species 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003809 water extraction Methods 0.000 description 2
- 244000171263 Ribes grossularia Species 0.000 description 1
- 235000002357 Ribes grossularia Nutrition 0.000 description 1
- 238000011010 flushing procedure Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Feedback Control In General (AREA)
- Flow Control (AREA)
Abstract
The invention discloses a city water source water supply scheduling method, which belongs to the relevant technical field of city water source supply and demand scheduling and comprises the following specific steps: obtaining a target function; constructing an RR relation matrix; constructing a water supply to water supply model; constructing an RU relation matrix; calculating the engineering capacity of each node; constructing a global equilibrium scheduling model based on the lowest water shortage rate; the invention provides a reservoir/gate pump-reservoir/gate pump, reservoir-gate pump-water plant user topology matrix algorithm on the basis of water quantity balance, and can support the calculation of water quantity distribution of a gate, pump and reservoir multi-node urban water supply system; and the urban water fluctuation analysis is carried out by combining the water demand requirement data of the water plant, so that the urban water supply scheduling scheme is guaranteed.
Description
Technical Field
The invention relates to the technical field related to urban water source supply and demand scheduling, in particular to a method for scheduling urban water source supply.
Background
More and more urban water supply safety faces increasingly serious challenges, local water sources are difficult to meet increasing water demand, and water transfer across boundary has become one of important means for improving urban water supply guarantee, such as Shenzhen city, the average total water supply amount of many years is about 19.73 hundred million m3In the middle ofThe external water diversion amount is about 15.51 hundred million meters3The water transfer amount across the border is up to 78% of the total amount of urban water supply. However, as diversion engineering systems become larger and larger, diversion pump stations and gate control nodes are increased continuously, hydraulic connections of all projects in regions are tighter, how to scientifically and reasonably realize efficient water quantity allocation among water source-water plant users is provided, an effective decision basis is provided for water resource scheduling and management, the method is a basis for guaranteeing urban water safety, and the method has important significance in urban water supply scheduling research.
In the traditional research, the urban local reservoir is used as a water storage node to be dispatched with a water plant, so that the structure is simpler; along with the water diversion engineering, the water network engineering is increasingly complicated, more transfer supply nodes appear in the water supply network, simultaneously, different water storage nodes, the gate pump node, the water supply scheduling coordination between the pipe network water delivery is more and more complicated, how to overcome many objects between water source nodes and users in the urban water supply system, the scientific allocation of multi-target water quantity, realize the guarantee target of urban global balanced water supply, face the challenge, therefore, the problem that technical personnel need to solve urgently in the field can be solved by comprehensively scheduling by utilizing various hydraulic engineering in research and development.
Disclosure of Invention
In view of the above, the invention provides a method for scheduling water supply of an urban water source, which can utilize various hydraulic engineering to perform scheduling to meet the requirement of urban water use.
In order to achieve the above purpose, the invention provides the following technical scheme:
a city water source water supply scheduling method comprises the following specific steps:
an objective function: acquiring a target function from actual monitoring data and engineering parameters of a reservoir, a water plant, a pump, a gate and a pipeline;
constructing an RR relation matrix: establishing an RR relation matrix according to the connection relation between the sorted nodes;
constructing a model for converting urban water supply into water supply: establishing a water supply to water supply amount model according to the RR relation matrix;
constructing an RU relation matrix: establishing an RU relation matrix according to the node sequencing and the water supply relation between each node and each computing unit;
calculating the engineering capacity of each node: calculating engineering capacity according to the constraint conditions;
constructing a minimum water shortage rate model: constructing a minimum water shortage rate model according to the RU relation matrix, the objective function and the engineering capacity;
global water plant equalization: and acquiring global equilibrium target values of a plurality of computing units in the whole area according to the lowest water shortage rate model of the computing units.
Preferably, in the construction of the RR relation matrix, the plurality of nodes are ordered according to hydraulic connection and importance of trunk and branch lines in the water quantity allocation process.
Preferably, matrix elements of the RR relationship matrix represent a connection relationship between two adjacent nodes, and 1 represents that a connection relationship exists between two adjacent nodes; 0 represents that no connection relation exists between two adjacent nodes.
Preferably, the model for converting urban water supply into water supply is specifically as follows:
(1) pump brake-reservoir transfer algorithm
Wherein: wm,rRepresenting the t time interval of the pumping station water lifting node/gate transmission node to the r water storage node to convert the water supply amount; qmRepresenting the total water supply quantity in the period t of the pumping node/gate transmission node of the mth pump station; qm,sRepresenting the water supply amount of the mth pumping station water lifting node/gate transmission node in the t period of the s calculation unit;
(2) reservoir-reservoir transfer algorithm
Wherein: wn,rRepresenting the time t of the nth water storage node to the r water storage node to supply water; h isnRepresents the water level of the nth water storage node, hrRepresenting the water level of the r-th water storage node, q1、q2、q3、q4Curve parameter, p, representing water level and delivery flow of the nth node of water storage1、p2、p3、p4Representing the water level and delivery flow curve parameters of the r water storage node.
Preferably, the row number in the RU relation matrix represents a computing unit, the column number represents the sorted nodes, and the matrix elements represent the water supply relation between the nodes and the computing unit; 1 represents that the node supplies water to the computing unit; 0 represents that the node is not supplying water to the computing unit.
Preferably, the nodes comprise a water storage node, a pump station water lifting node and a gate transmission node.
Preferably, the constraint conditions are specifically:
1) water storage node water balance condition:
Vt1,n=Vt0,n+Qt,n-Qt,rn-Et,n-St,n; (3)
in the formula: vt1,nFor the nth water storage node storage capacity at the end of the t period, Qt,nThe water inflow amount of the nth water storage node in the t period, Qt,rnIs the outlet flow of the nth water storage node in the time period t, Vt0,nFor the initial nth storage node storage capacity at t time interval, Et,nFor the evaporation loss amount of the water storage node n in the time period of t, St,nLeakage loss amount of the water storage node n in a time period t;
2) and pump station water lifting node constraint conditions:
in the formula: b isn,tSupplying water for the water lifting node of the nth pump station in a time period t;the maximum water supply capacity of the nth pump station in the t period;
3) gate transmission node constraint conditions:
in the formula: zn,tFor the nth gate to transmit a node ZnAt the time of the t-th period of the excessive flow,for the t-th time period, the gate transmits a node ZnThe maximum flow rate of;
4) pipe network water volume transmission constraint conditions:
Gi,tnode pipeline G for pumping stationiWater supply amount at time t;node pipeline G for pumping stationiMaximum water delivery capacity at time t;
5) calculating the water balance condition of the unit:
in the formula: di,tTo calculate the water demand, Q, of unit i during the t-th time periodBji,tThe water quantity, Q, supplied to the computing unit i for the pumping node of the jth pump station during the period tZki,tSupplying the water amount of the computing unit i to the kth gate transmission node in the time period t; k is the total number of gate transmission nodes.
Preferably, the minimum water shortage model is specifically as follows:
in the formula: l (xt) is a water supply safety target, Di,tThe water demand of the ith calculation unit in the t period; qji,tNet water supply to the ith computing unit for the jth water supply node during the tth period; c is the total number of the water supply nodes; the water supply node comprises a water storage node, a pump station water lifting node and a gate transmission node; i is the total number of the calculation units; t is the total number of simulation periods.
Preferably, the global water plant balance is specifically as follows:
in the formula: f (xt) is the regional global equalization target, li,tCalculating the water shortage rate of the ith calculation unit in the t period;the average value of the water shortage rate of the unit in the t period is calculated.
Preferably, the method further comprises the step of checking the model,
Re=(Mt,out-Mt,real)/Mt,real; (10)
in the formula: reAs a relative error, Mt,outFor the calculation of the model over the time period t, Mt,realThe actual scheduling result is t period.
According to the technical scheme, compared with the prior art, the urban water source water supply scheduling method is provided, a reservoir/gate pump-reservoir/gate pump and reservoir-gate pump-water plant user topology matrix algorithm is provided on the basis of water quantity balance, and calculation of water quantity distribution of an urban water supply system with multiple nodes of gates, pumps and reservoirs can be supported; the urban water fluctuation analysis can be carried out by combining water demand prediction data of a water plant, and a guarantee is provided for an urban water supply scheduling scheme.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 is a general flow diagram of the present invention;
FIG. 2 is a partial schematic view of the water supply system of the present invention;
FIG. 3 is a diagrammatic view of Shenzhen water supply network;
FIG. 4(a) is a diagram showing a result of calculation of water quantity of a water source of a Meilin water plant; FIG. 4(b) is a diagram showing a result of water source calculation in a large flush plant;
FIG. 5(a) is a diagram showing the water content of the water divided before the sand-bay aqueduct; FIG. 5(b) is a diagram showing the water diversion amount from the Shawan aqueduct to the Xili reservoir; FIG. 5(c) is a water content chart of a water source part in the north; FIG. 5(d) is a water content diagram of the iron ore communicating section;
wherein, the water diversion trunk line project 1; a pump station 2; a first gate 3; a first water plant 4; a first reservoir 5; a second shutter 6; a second reservoir 7; a second water plant 8; water diversion branch line engineering 9.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the invention discloses a city water source water supply scheduling method, which provides reservoir/gate pump-reservoir/gate pump and reservoir-gate pump-water plant user topological algorithm based on the water quantity balance principle according to the characteristics of a city water supply network, and comprises the following steps as shown in figure 1 of the attached drawings, and the target function is as follows: estimating a target function according to historical data; constructing an RR relation matrix: establishing an RR relation matrix according to the connection relation between the sequencing nodes; constructing a water supply to water supply model: establishing a water supply to water supply amount model according to the RR relation matrix; constructing an RU relation matrix: establishing an RU relation matrix according to the node sequencing and the water supply relation among the nodes; calculating the engineering capacity of each node: calculating engineering capacity according to the constraint conditions; constructing a minimum water shortage rate model: constructing a water supply scheduling model with the lowest water shortage rate according to the RU relation matrix, the objective function and the engineering capacity; and checking the model precision according to the relative error index.
For urban water supply system of overseas diversion, the supply side: the water supply source water quantity is lifted by a water diversion source pump station and is conveyed to each branch project through a water diversion gate, and the surplus water quantity regulation and reservoir storage transfer supply are realized by combining each reservoir project along the way; a demand side: the system has the condition that the reservoir directly supplies water to the water plant and the water plant directly takes water from the water diversion engineering pipeline. The method is based on the water quantity allocation process from a water source to a city water supply plant, and is used for describing the water quantity transmission condition of each engineering node, and a system generalization method is adopted on the basis of the original model reservoir node to generalize the reservoir into water storage nodes; a water diversion project for water extraction through a pump station is generalized into a pump station water extraction node; the diversion project supplying water through the gate is generalized into gate transmission nodes, and each water plant is generalized into an independent computing unit, wherein the water plant is communicated with a corresponding water supply source through a pipe network.
The urban water supply is different from natural confluence, and the water quantity is gathered and transported according to the water diversion engineering capacity and can be directionally and quantitatively distributed according to different requirements. In order to depict the distribution process of urban water quantity, the method combines the theory of a topological structure, aims at the nodes of a reservoir, a pump station and a gate, sequences different nodes according to the principle that the hydraulic connection of different nodes in the water quantity distribution process and each diversion project follows from upstream to downstream, firstly, a main line project and then a branch line project, wherein the nodes of the reservoir are sequenced from a source to a tail end according to a water supply main line, the nodes of the gate and the pump station follow the condition that the node of the main line of the diversion project is in front and the node of the branch line of the diversion project is in back according to the pipeline of the diversion project, and establishes an RR (reservoir/gate pump-reservoir/gate pump) relation matrix. In the water supply system partial schematic diagram shown in the attached figure 2, a water diversion trunk line project 1 is connected with a pump station 2, the pump station is respectively connected with a first gate 3 and a first reservoir 5, the first gate 3 is connected with a first water plant 4, the first reservoir 5 is connected with a second gate 6, the second gate 6 is connected with a second reservoir 7, the second reservoir 7 is respectively connected with a second water plant 8 and a water diversion branch line project 9, different nodes in the attached figure 2 are sequenced according to the principle, the pump station 2, the first gate 3, the first reservoir 5, the second gate 6 and the second reservoir 8 are sequentially arranged, an RR relation matrix is constructed according to the node arrangement sequence, the row number and the column number of the matrix correspond to the node arrangement sequence, the matrix element represents the communication relationship between a calculation node (upstream node) in the sequence and a node in the sequence, and 1 represents that the communication relationship exists between two nodes (when the node has no downstream node, considered to be in communication with itself); 0 represents that no connection relation exists between two nodes, and the RR relation matrix specifically comprises:
note: i is a row number, j is a column number (the same below), different nodes are sequenced according to the upstream and the downstream, and the RR relation matrix generated correspondingly is a lower triangular matrix;
wherein rr isi,j1 represents that a communication relation exists between the j node and the i node; if the number is 0, the connection relation does not exist among the nodes, and C represents the number of the nodes.
When water quantity calculation is carried out, the function that a pump station and a gate node do not have water quantity regulation and storage is considered to be mainly restricted by engineering capacity, and the reservoir can regulate surplus water diversion quantity and supply water quantity again when daily scheduling is carried out, so that the calculation of converting urban water supply into water supply is carried out according to the constructed topological relation matrix:
(1) pump brake-reservoir transfer algorithm
Wherein: wm,rRepresenting the t time interval of the pumping station water lifting node/gate transmission node to the r water storage node to convert the water supply amount; qmRepresenting the total water supply quantity in the period t of the pumping node/gate transmission node of the mth pump station; qm,sRepresenting the water supply amount of the mth pumping station water lifting node/gate transmission node in the t period of the s calculation unit;
(2) reservoir-reservoir transfer algorithm
Wherein: wn,rRepresenting the time t of the nth water storage node to the r water storage node to supply water; h isnRepresents the water level of the nth water storage node, hrRepresenting the water level of the r-th water storage node, q1、q2、q3、q4Curve parameter, p, representing water level and delivery flow of the nth node of water storage1、p2、p3、p4Representing the water level and delivery flow curve parameters of the r water storage node.
In an urban water supply dispatching system, in order to improve the reliability of water supply of a water plant, the situation that a reservoir directly fetches water from the water plant and the water plant directly supplies water by a diversion project may exist, the model calculates to generalize the water plant into calculation units, the calculation units are sorted according to nodes, and an RU (reservoir-gate pump-water plant user) relation matrix is established according to the corresponding water supply relation of the water plant in the figure 2:
wherein, ru m,j1 represents that j nodes directly supply water like m water plants; ru is a Chinese characterm,jIf the number is 0, the water plant does not take water from the corresponding node, C represents the number of the nodes, and N represents the number of the computing units;
the RU relation matrix comprises row numbers, column numbers and matrix elements, wherein the row numbers represent different computing units, the column numbers represent nodes after sequencing, and the matrix elements represent water supply relations between the nodes and the computing units; 1 represents that the node supplies water to the computing unit; 0 represents that the node is not supplying water to the computing unit.
And determining the selection of the water taking path of the computing unit according to the RU topological matrix, and carrying out total water consumption constraint on the corresponding computing unit according to the actual engineering capacity of the node to realize water supply allocation calculation.
The engineering capacity is realized by depending on constraint conditions, and the constraint conditions are as follows:
(1) water storage node water balance condition:
Vt1,n=Vt0,n+Qt,n-Qt,rn-Et,n-St,n; (3)
in the formula: vt1,nFor the nth water storage node storage capacity at the end of the t period, Qt,nThe water inflow amount of the nth water storage node in the t period, Qt,rnIs the outlet flow of the nth water storage node in the time period t, Vt0,nFor the initial nth storage node storage capacity at t time interval, Et,nFor the evaporation loss amount of the water storage node n in the time period of t, St,nLeakage loss amount of the water storage node n in a time period t;
(2) and pump station water lifting node constraint conditions:
in the formula: b isn,tSupplying water for the water lifting node of the nth pump station in a time period t;the maximum water supply capacity of the nth pump station in the t period;
(3) gate transmission node constraint conditions:
in the formula: zn,tFor the nth gate to transmit a node ZnAt the time of the t-th period of the excessive flow,for the t-th time period, the gate transmits a node ZnThe maximum flow rate of;
(4) pipe network water volume transmission constraint conditions:
wherein,Gi,tnode pipeline G for pumping stationiWater supply amount at time t;node pipeline G for pumping stationiMaximum water delivery capacity at time t;
(5) calculating the water balance condition of the unit:
in the formula: di,tTo calculate the water demand, Q, of unit i during the t-th time periodBji,tThe water quantity, Q, supplied to the computing unit i for the pumping node of the jth pump station during the period tZki,tSupplying the water quantity of the computing unit i to the jth gate transmission node; k is the total number of gate transmission nodes.
Constructing a minimum water shortage rate model according to the conditions, specifically comprising the following steps:
in the formula: l (xt) is a water supply safety target, Di,tThe water demand of the ith calculation unit in the t period; qji,tNet water supply to the ith computing unit for the jth water supply node during the tth period; c is the total number of the water supply nodes; the water supply node comprises a water storage node, a pump station water lifting node and a gate transmission node; i is the total number of the calculation units; t is the total number of simulation periods.
The global water plant equilibrium model is constructed according to the lowest water shortage rate model, and specifically comprises the following steps:
in the formula: f (xt) is the regional global equalization target, li,tCalculating the water shortage rate of the ith calculation unit in the t period;the average value of the water shortage rate of the unit in the t period is calculated.
In the embodiment, the Shenzhen city is taken as an example, the urban water supply of the Shenzhen city is mainly domestic water for residents, the water supply source mainly depends on overseas water diversion due to the small regulation and storage capacity of the local reservoir, and the average overseas water diversion water quantity of many years accounts for 78% of the total urban water supply quantity. The water diversion project is used for communicating 29 main water supply reservoirs in the whole city to supply and regulate water quantity so as to meet daily water supply quantity requirements of 47 water plants in the whole city.
In order to improve the guarantee condition of water supply in a district, part of water plants adopt double water sources for water taking in daily scheduling, and simultaneously, in order to increase the reliability of a water supply source, the urban water supply network realizes the communication between the Dongjiang water source project and the Dongdeep water supply project through a sand-bay aqueduct. However, under the influence of engineering benefits, the aqueduct must run at full load, and the water of the east river water supply engineering is supplemented by the water of the east river under the condition that the residual water of the east river trunk line junction is insufficient. Therefore, the calculation of the urban diversion water quantity needs to be combined with the time interval and the end storage capacity of each reservoir during daily dispatching, the influence of engineering factors is considered, and the water quantity allocation among different diversion projects, between the diversion project and the reservoir and between the reservoir and the reservoir is realized through gate control so as to ensure the normal water supply of the water works in the whole city.
Aiming at the current situation of the Shenzhen city urban water supply network, the urban water supply scheduling model conforming to the Shenzhen water supply system is established by taking the improved model provided by the above as a framework. The model inputs include: basic engineering capacity information of reservoirs, gates, pumps, diversion engineering pipe networks and the like; an RR relation matrix; ③ RU relationship matrix; fourthly, daily water demand of each water plant user in the dispatching period; scheduling initial storage capacity of the day reservoir at the beginning of the time interval; sixthly, the total amount of water diversion outside the initial period of time; and the initial state information of each gate pump.
And constructing the model according to the principle, and verifying the calculation result of the water supply and water allocation of the whole market for 7 continuous days of the model by adopting an actual scheduling water regime report form of the whole market from 30 days to 6 days in month 7 in 2020. In order to determine the deviation of the calculation result of the model in each time period, the calculation accuracy of the model is checked by adopting a relative error index.
Re=(Mt,out-Mt,real)/Mt,real; (10)
In the formula: reAs a relative error, Mt,outFor the calculation of the model over the time period t, Mt,realThe actual scheduling result is t period.
And calculating the total water demand in each time period in the whole city. The actual water demand of the water plant in the simulation period is 3633 ten thousand meters3Model calculation total water supply amount of water lifting node of pump station is 3626 km3The relative error between the total water demand of the water plant and the total supply quantity calculated by the model in the time period is 0.19 percent, and the calculation error value of the model in each time period is within 0.5 percent. The model calculation result is better, and the requirement of water resource management precision is met.
TABLE 1 Total Water supply distribution in Water works
Time period | Water plant time period water demand total (ten thousand meters)3) | Model calculation Total supply (ten thousand m)3) | RE |
1 | 535 | 534 | -0.19% |
2 | 537 | 535 | -0.37% |
3 | 525 | 523 | -0.38% |
4 | 516 | 515 | -0.19% |
5 | 506 | 504 | -0.40% |
6 | 502 | 501 | -0.20% |
7 | 512 | 514 | 0.39% |
Total up to | 3633 | 3626 | -0.19% |
Water plant-water source water intake verification
Aiming at the situation that double water sources of water plants simultaneously take water in a Shenzhen city water supply scheduling system, a plum water plant and a Dachong water plant are taken as examples for verification of calculation results. Wherein, the Meilin water plant directly fetches water from the Dongjiang water source engineering main line, and takes the Dongjiang water supply engineering from the Shenzhen reservoir to supply water; the Dachong water works take Dongjiang water source engineering from the Xilishui reservoir to supply water, and take DongShen engineering from the Shenzhen reservoir to supply water. Sequencing the pump station water lifting nodes by adopting an RR relation matrix, calculating the water supply amount of the corresponding pump station water lifting nodes of the water plant, analyzing the water taking rule of the water plant by using a model because of lack of actual pipeline overflowing data, and finally setting the historical maximum water taking amount data of the water plant as pipeline overflowing capacity constraint; the results of the water diversion source calculation in the water plant are shown in fig. 4(a) and 4 (b).
In order to determine the specific water quantity deviation of double water sources water intake of a water plant, absolute errors are adopted for verification and calculation
The results are shown in Table 2. Wherein, the water yield calculation error of the Meilin water plant is the largest in the 2 nd period, and the water yield calculation result of the Dongdi water supply project is larger by 3.45 ten thousand meters3The water quantity calculation result of Dongjiang water source engineering is smaller by 3.07 ten thousand meters3(ii) a The maximum error of the calculation result of the large water flushing plant is the same as the 2 nd time period, and the east deep water volume is smaller than 1.96 ten thousand meters3The water quantity of Dongjiang is 1.78 ten thousand meters larger3. The model can be seen to have better calculation accuracy for the water intake quantity of the double water sources of the water plant.
TABLE 2 Water works diversion source calculation results
Note: dong Shen is water supply engineering; dongjiang river as Dongjiang river water source engineering
Reservoir dispatching end-of-term storage capacity verification
The model construction takes the storage capacity of a reservoir at the first day of a dispatching period as an initial storage capacity, water allocation calculation of the total city of 7 continuous days is carried out, the storage capacity calculation results of the main water supply reservoirs in the total city at the end of the dispatching period are shown in a table 3, the relative errors of the storage capacities at the end of the reservoir calculation period are all within 5%, wherein the largest errors are madder pit reservoir-3.94% (the minus sign represents that the calculation result is smaller than the actual value, the same below), gooseneck reservoir-3.13% and rocky reservoir 2.21%. The reason for the phenomenon is that the madder pit reservoir is a primary regulation reservoir of a water diversion project of a north line, the gooseberry gate is communicated with the gooseneck reservoir, after the gooseneck reservoir supplies water to a bright water plant, surplus water is transferred to the rock reservoir through the goosestone gate, and water supply of the water plant in a district is finally met, namely three reservoir nodes are in a series connection relation.
TABLE 3 calculation of reservoir capacity at end of time period
Engineering water distribution verification
In order to verify the feasibility of the model improvement method, 4 key water diversion projects are selected to verify the water distribution result among different projects: dongjiang water source engineering (dividing the sand bay aqueduct into a front sand bay aqueduct and a rear sand bay aqueduct by taking the sand bay aqueduct as a node); northern water source engineering (taking an upper Cambodia pump station as a starting point and a rock reservoir as a terminal point); model calculation is performed on the iron hillock and the rock reservoir, and the calculation results are shown in fig. 5(a) to 5 (d).
Compared with an actual scheduling result, the error of the water distribution calculation result after the sand bay aqueduct is found to be 7 percent, because the water plant unit behind the aqueduct has the condition of simultaneously taking water from double water sources, because the overflow data of a pipe network is lacked, the accumulated error occurs during the water distribution calculation in a continuous time period, and the engineering water distribution result in the section has larger deviation, besides, the model has higher calculation precision on other engineering water distribution results, and the deviation of the calculation result and the actual scheduling water is within 1 percent.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. The device disclosed by the embodiment corresponds to the method disclosed by the embodiment, so that the description is simple, and the relevant points can be referred to the method part for description.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (9)
1. A city water source water supply scheduling method is characterized by comprising the following specific steps:
an objective function: acquiring a target function from actual monitoring data and engineering parameters of a reservoir, a water plant, a pump, a gate and a pipeline;
constructing an RR relation matrix: establishing an RR relation matrix according to the connection relation between the sorted nodes;
constructing a model for converting urban water supply into water supply: establishing a water supply to water supply amount model according to the RR relation matrix;
constructing an RU relation matrix: establishing an RU relation matrix according to the node sequencing and the water supply relation between each node and each computing unit;
calculating the engineering capacity of each node: calculating engineering capacity according to the constraint conditions;
constructing a minimum water shortage rate model: constructing a minimum water shortage rate model according to the RU relation matrix, the objective function and the engineering capacity;
global water plant equalization: acquiring global equilibrium target values of a plurality of computing units in the whole area according to the lowest water shortage rate model of the computing units;
wherein, the city water supply changes the water supply model specifically to be:
(1) pump brake-reservoir transfer algorithm
Wherein: wm,rRepresenting that the pumping station water lifting node/gate transmission node of the mth pump station converts the water supply amount of the water storage node in the time t; qmRepresenting the total water supply of the pumping node/gate transmission node of the mth pump station in the period t; qm,sRepresenting the water lifting node/gate transmission node of the mth pump station to supply water to the mth computing unit in the time period t;
(2) reservoir-reservoir transfer algorithm
Wherein: wn,rRepresenting the time t of the nth water storage node to the r water storage node to supply water; h isnRepresents the water level of the nth water storage node, hrRepresenting the water level of the r-th water storage node, q1、q2、q3、q4Curve parameter, p, representing water level and delivery flow of the nth node of water storage1、p2、p3、p4Representing the water level and delivery flow curve parameters of the r water storage node.
2. The method of claim 1, wherein in constructing the RR relationship matrix, the nodes are ordered according to hydraulic connection and importance of trunk and branch lines in water quantity allocation.
3. The urban water source water supply scheduling method according to claim 2, wherein a matrix element of the RR relation matrix represents a connection relation between two adjacent nodes, and 1 represents that a connection relation exists between two adjacent nodes; 0 represents that no connection relation exists between two adjacent nodes.
4. The city water source water supply scheduling method of claim 3, wherein matrix elements in the RU relationship matrix represent water supply relationships between nodes and the computing units; 1 represents that the node supplies water to the computing unit; 0 represents that the node is not supplying water to the computing unit.
5. The urban water source water supply scheduling method of claim 4, wherein the nodes comprise a water storage node, a pump station water lifting node and a gate transmission node.
6. The urban water source water supply scheduling method according to claim 5, wherein the constraint conditions are specifically:
1) water storage node water balance condition:
Vt1,n=Vt0,n+Qt,n-Qt,rn-Et,n-St,n; (3)
in the formula: vt1,nFor the nth water storage node storage capacity at the end of the t period, Qt,nThe water inflow amount of the nth water storage node in the t period, Qt,rnIs the outlet flow of the nth water storage node in the time period t, Vt0,nFor the initial nth storage node storage capacity at t time interval, Et,nFor the evaporation loss amount of the water storage node n in the time period of t, St,nLeakage loss amount of the water storage node n in a time period t;
2) and pump station water lifting node constraint conditions:
in the formula: b isn,tSupplying water for the water lifting node of the nth pump station in a time period t;the maximum water supply capacity of the nth pump station in the t period;
3) gate transmission node constraint conditions:
in the formula: zn,tFor the nth gate to transmit a node ZnAt the time of the t-th period of the excessive flow,for the t-th time period, the gate transmits a node ZnThe maximum flow rate of;
4) pipe network water volume transmission constraint conditions:
in the formula, Gi,tNode pipeline G for pumping stationiWater supply amount at time t;node pipeline G for pumping stationiMaximum water delivery capacity at time t;
5) calculating the water balance condition of the unit:
in the formula: di,tTo calculate the water demand, Q, of unit i during the t-th time periodBji,tThe water quantity, Q, supplied to the computing unit i for the pumping node of the jth pump station during the period tZki,tSupplying the water amount of the computing unit i to the kth gate transmission node in the time period t; k is the total number of gate transmission nodes.
7. The urban water source water supply scheduling method according to claim 6, wherein the lowest water shortage rate model is specifically:
in the formula: l (xt) is a water supply safety target, Di,tThe water demand of the ith calculation unit in the t period; qji,tNet water supply to the ith computing unit for the jth water supply node during the tth period; c is the total number of the water supply nodes; the water supply node comprises a water storage node, a pump station water lifting node and a gate transmission node; i is the total number of the calculation units; t is the total number of simulation periods.
8. The urban water source water supply scheduling method according to claim 7, wherein the global water plant balance is specifically:
9. The municipal water supply scheduling method according to claim 8, further comprising a verification model,
Re=(Mt,out-Mt,real)/Mt,real; (10)
in the formula: reAs a relative error, Mt,outFor the calculation of the model over the time period t, Mt,realThe actual scheduling result is t period.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110728549.2A CN113449993B (en) | 2021-06-29 | 2021-06-29 | Urban water source water supply scheduling method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110728549.2A CN113449993B (en) | 2021-06-29 | 2021-06-29 | Urban water source water supply scheduling method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113449993A CN113449993A (en) | 2021-09-28 |
CN113449993B true CN113449993B (en) | 2022-02-08 |
Family
ID=77814028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110728549.2A Active CN113449993B (en) | 2021-06-29 | 2021-06-29 | Urban water source water supply scheduling method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113449993B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114707913B (en) * | 2022-06-02 | 2022-09-13 | 长江水利委员会长江科学院 | Method and device for determining water supply scheme of reservoir |
CN115271453B (en) * | 2022-07-29 | 2023-06-16 | 中国水利水电科学研究院 | Urban raw water supply allocation path identification method, system and storable medium |
CN116562600B (en) * | 2023-07-11 | 2023-09-19 | 中关村科学城城市大脑股份有限公司 | Water supply control method, device, electronic equipment and computer readable medium |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109658287A (en) * | 2018-12-27 | 2019-04-19 | 中国水利水电科学研究院 | A kind of basin water dispatching method evenly distributed based on water resource space-time |
CN110991687A (en) * | 2019-09-26 | 2020-04-10 | 深圳市东深电子股份有限公司 | Water resource scheduling optimization method based on empirical model |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140014763A (en) * | 2012-07-26 | 2014-02-06 | 한국전력공사 | Integrated supply apparatus for energy and fresh water used microgrid |
CN106204333B (en) * | 2016-07-20 | 2018-08-03 | 中国水利水电科学研究院 | A kind of water resource dispatching method based on Diabatic slow wave |
EP3450638A1 (en) * | 2017-08-30 | 2019-03-06 | Siemens Aktiengesellschaft | Control of a water distribution network |
CN108427459B (en) * | 2018-05-18 | 2020-12-04 | 成都市自来水有限责任公司 | Automatic control method for piston type flow regulating valve of multi-water-source annular water supply pipe network |
CN110174843B (en) * | 2019-05-30 | 2020-03-10 | 中国水利水电科学研究院 | Intelligent regulation and control method for water used in irrigation area |
-
2021
- 2021-06-29 CN CN202110728549.2A patent/CN113449993B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109658287A (en) * | 2018-12-27 | 2019-04-19 | 中国水利水电科学研究院 | A kind of basin water dispatching method evenly distributed based on water resource space-time |
CN110991687A (en) * | 2019-09-26 | 2020-04-10 | 深圳市东深电子股份有限公司 | Water resource scheduling optimization method based on empirical model |
Also Published As
Publication number | Publication date |
---|---|
CN113449993A (en) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113449993B (en) | Urban water source water supply scheduling method | |
CN108764573B (en) | Inland arid region-oriented multidimensional equilibrium configuration system for water resources | |
CN105913146B (en) | Southern wet area water resource optimal allocation system | |
Koutsoyiannis et al. | A decision support tool for the management OF multi‐reservoir systems 1 | |
Gomes et al. | Decision support system to divide a large network into suitable District Metered Areas | |
CN110689283B (en) | River system and canal system staggered water supply simulation method and device based on underground water model | |
CN109992909B (en) | Tree-shaped river network cascade reservoir hydrodynamic water quality sediment coupling simulation method and system | |
CN106502220B (en) | Cascade pumping station water-carriage system optimization operation control coupling coordination approach and system | |
CN112101818B (en) | Sponge city flood optimal scheduling method suitable for complex hydraulic connection | |
CN106503313B (en) | Series-parallel cascade pumping station water-carriage system operational efficiency calculation method and system | |
CN113256028B (en) | Emergency water supply allocation method for large-range long-duration drought | |
CN110633849B (en) | Water resource configuration method based on heuristic large system decomposition coordination idea | |
WO2015033269A1 (en) | A control system for operation of irrigation canals | |
CN107506894A (en) | A kind of water power group's dispatching method for considering non-constant coupling constraint | |
CN110533326A (en) | A kind of water-related Eco-environmental Carrying Capacity estimates regulator control system | |
CN114663254A (en) | Water resource-grain-energy-ecological cooperative regulation and control method | |
CN115271304A (en) | Cross-basin water transfer engineering optimal scheduling method based on climate change | |
CN109472463A (en) | A kind of more peak load regulation network dispatching methods of especially big basin water station group and system | |
CN109919397B (en) | Checking and optimizing method for drainage capacity of plain river network area | |
CN114357868A (en) | Multi-target cooperative scheduling method and device for complex flood control system | |
Belsnes et al. | Unit commitment in hydro power operation scheduling | |
CN112632733A (en) | Arrangement method of newly added flow meters of water supply network for maximizing monitoring information | |
Chang et al. | Water Dispatch Model for Middle Route of a South‐to‐North Water Transfer Project in China 1 | |
CN104047257B (en) | A kind of ecological characteristic storages of reservoir defining method of single service reservoir of Environment Oriented protection | |
Frijns et al. | Management tools for hydro energy interventions in water supply systems |
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 |