Multi-water-source urban raw water combined dispatching scheme compilation method

Abstract

The invention relates to a multi-water-source urban raw water combined scheduling scheme compiling method, which is characterized in that a multi-target scheduling model considering water supply vacancy and comprehensive benefits is constructed, so that a compiled scheduling scheme can consider regional water supply vacancy and comprehensive benefits, and meanwhile, determined constraint conditions are corrected through a constraint correction strategy, so that the scheduling efficiency is improved. The accuracy of the scheduling scheme is improved, so that the utilization rate of water resources is improved, and the probability of flood and drought disasters is reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of water resource dispatching and management, and in particular to a method for compiling a joint dispatching plan for raw water in a city with multiple water sources. Background Art

[0002] Cross-regional water transfer projects and joint optimization of multiple water sources are important means to achieve optimal allocation of water resources, and are also one of the key measures to ensure the safety of water supply in cities with seasonal and regional water shortages. Urban raw water reservoir group water supply scheduling utilizes the differences in storage capacity and runoff characteristics of multiple reservoirs, and through joint optimization scheduling, it gives full play to the comprehensive role of reservoir groups in hydrological compensation and storage capacity compensation, thereby solving the uneven spatial distribution of water resources, improving the utilization rate of water resources in the basin, and achieving the sustainability of water resources. However, the scheduling methods in the existing technology are often relatively simple, can only meet a certain scheduling target, have low comprehensive performance, and cannot make full use of water resources. Summary of the invention

[0003] The technical problem to be solved by the present invention is to provide a method for compiling a joint scheduling plan for raw water in a multi-water source city, which is beneficial to improving the comprehensive performance of the scheduling plan and improving the utilization rate of water resources.

[0004] The technical solution adopted by the present invention to solve the above technical problems is: a method for compiling a joint scheduling plan for raw water in a multi-water source city, comprising:

[0005] S1, comprehensively consider the water supply shortage and comprehensive benefits of the study area to construct the objective function of the scheduling model;

[0006] S2, determining the constraints of the dispatching model by considering the parameters of regional reservoirs and water pipelines;

[0007] S3, the differential optimization algorithm and the stepwise optimization algorithm are combined to solve the scheduling model, obtain the optimal raw water supply flow process of the raw water supply pipeline in each time period, and then formulate a joint scheduling plan for raw water in multiple water source cities.

[0008] Furthermore, the scheduling model objective function constructed in step S1 includes a first objective function, which is as follows:

[0009]

[0010] In the formula, OBJ1 represents the first objective function value, V represents the overall water supply deficit of the region, n and m represent the number of reservoirs and water plants respectively, i and j represent reservoirs and water plants respectively, t represents the month, and S i,t represents the water supply of the i-th reservoir in the t-th month, D j,t Represents the water demand of the j-th water plant in the t-th month.

[0011] Furthermore, the scheduling model objective function constructed in step S1 includes a second objective function, which is as follows:

[0012]

[0013] In the formula, OBJ2 represents the second objective function value, that is, the maximum comprehensive benefit, B w,i,t represents the water supply benefit of the i-th reservoir in the t-th month, n represents the number of reservoirs, B g,i,t represents the power generation benefit of the i-th reservoir in the t-th month.

[0014] Furthermore, the water supply benefit of the i-th reservoir in the t-th month is B w,i,t The formula is as follows:

[0015] B w,i,t =(Q w,i,t ×T w,i,t ×24×3600×P w ) / 10000

[0016] In the formula, Q w,i,t represents the average water supply flow of the i-th reservoir in the t-th month, m 3 / s,T w,i,t represents the number of days the i-th reservoir supplies water in the t-th month, P w Indicates the raw water price.

[0017] Furthermore, the power generation benefit of the i-th reservoir in the t-th month is B g,i,t The formula is as follows:

[0018] B g,i,t =(N g,i,t ×T g,i,t ×24×P g ) / 10000

[0019] =(K g,i ×H g,i,t ×Q g,i,t ×T g,i,t ×24×P g ) / 10000

[0020] Where N g,i,t represents the average output of the i-th reservoir in the t-th month, K g,i is the output coefficient, H g,i,t represents the average head difference of the ith reservoir in the tth month, Q g,i,t represents the average power generation flow of the i-th reservoir in the t-th month, T g,i,t represents the number of days when the ith reservoir generates electricity in the tth month, P g Indicates the electricity price.

[0021] Furthermore, the scheduling model constraints determined in step S2 include:

[0022] Pipeline water delivery capacity constraints: SD i,j,t ≤SD i,j,t,max

[0023]

[0024] Reservoir water level constraint: Z i,t,min ≤Z i,t ≤Z i,t,max

[0025] Reservoir flow constraint: Q i,t,min ≤Q i,t ≤Q i,t,max

[0026] Reservoir output constraint: N i,t,min ≤N i,t ≤N i,t,max

[0027] Reservoir water level fluctuation constraint: ΔZ i,t,min ≤ΔZ i,t ≤ΔZ i,t,max

[0028] Variable non-negative constraint: SD i,j,t ≥0,D j,t ≥0,S i,t ≥0

[0029] In the above formula, SD i,j,t represents the upper limit of water flow from the i-th reservoir to the j-th water plant in the t-th month, SD i,t,max The upper limit of water flow from the i-th reservoir to all water plants in the t-th month, Z i,t,min and Z i,t,max They represent the lower and upper limits of the water level constraints of the i-th reservoir in the t-th month, Q i,t,min and Q i,t,max They represent the lower and upper limits of the outflow of the i-th reservoir in the t-th month, respectively. i,t,min and N i,t,max They represent the minimum and maximum guaranteed output of the i-th reservoir in the t-th month, ΔZ i,t,min and ΔZ i,t,max Represents the lower and upper limits of the water level variation of the ith reservoir in the tth month.

[0030] Furthermore, after executing the joint differential optimization algorithm and the stepwise optimization algorithm described in step S3 to solve the scheduling model, it also includes a step of correcting the minimum downstream flow of the reservoir, specifically including: when the water supply mode is water extraction from the reservoir, if the upstream water can meet the water supply needs, no correction is made; if the upstream water cannot meet the water supply needs, the minimum downstream flow of the reservoir is updated to the sum of the upstream water and the water supply flow; when the water supply mode is water extraction from the reservoir, if the upstream water can meet the water supply needs, no correction is made; if the upstream water cannot meet the water supply needs, the minimum downstream flow of the reservoir is updated to the water supply flow.

[0031] Furthermore, after executing the joint differential optimization algorithm and the stepwise optimization algorithm described in step S3 to solve the scheduling model, it also includes the step of correcting the water flow of the water supply pipeline, specifically including: when the water supply mode is water intake from the reservoir, determine whether the reservoir capacity meets the water intake conditions, if so, no correction is required, if not, the water supply of the water supply pipeline of this reservoir needs to be reduced according to the proportion of the project water supply capacity to balance the water volume, when the water supply mode is water intake from the reservoir, determine whether the water supply flow is greater than the downstream flow of the reservoir, if so, reduce the water flow of the water supply pipeline of this reservoir according to the proportion of the project water supply capacity to balance the water volume, if not, no correction is made.

[0032] Furthermore, before executing step S1, a step of generalizing the topological relationship of the water supply system is also included.

[0033] Compared with the prior art, the advantages of the present invention are:

[0034] The method of the present invention constructs a multi-objective scheduling model that considers water supply shortages and comprehensive benefits, so that the compiled scheduling plan can take into account regional water supply shortages and comprehensive benefits. At the same time, a constraint correction strategy is used to improve the accuracy of the scheduling plan, thereby improving the utilization rate of water resources and reducing the probability of water and drought disasters. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] Figure 1 This is a generalized topological diagram of the joint water supply from multiple water sources in a central urban area of ​​a certain city in an embodiment of the present invention.

[0036] Figure 2 Schematic diagram of the topological relationship traversal process of a central urban area in an embodiment of the present invention.

[0037] Figure 3 This is a comparison chart of the reservoir water levels of the Baixi Reservoir using the existing technology and the scheduling method of the present invention.

[0038] Figure 4 This is a comparison chart of the reservoir water levels of the Hengshan Reservoir using the prior art and the scheduling method of the present invention.

[0039] Figure 5This is a comparison chart of the reservoir water levels of the Hengxi Reservoir using the existing technology and the scheduling method of the present invention.

[0040] Figure 6 This is a comparison chart of the reservoir water levels of the Sanxipu Reservoir using the prior art and the scheduling method of the present invention.

[0041] Figure 7 This is a comparison chart of the reservoir water levels of the Tingxia Reservoir using the prior art and the scheduling method of the present invention.

[0042] Figure 8 This is a comparison chart of the reservoir water levels of Zhougongzhai Reservoir dispatched by the prior art and the dispatching method of the present invention.

[0043] Fig. 9 This is a comparison chart of the reservoir water levels of Jiaokou Reservoir dispatched by the prior art and the dispatching method of the present invention. DETAILED DESCRIPTION

[0044] The present invention is further described in detail below in conjunction with the embodiments of the drawings, examples of which are shown in the drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the drawings are exemplary and are intended to be used to explain the present invention, and should not be construed as limiting the present invention.

[0045] In the description of the present application, it should be noted that for directional words, such as the terms "center", "lateral", "longitudinal", "length", "width", "thickness", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", "clockwise", "counterclockwise" and the like indicating directions and positional relationships are based on the directions or positional relationships shown in the drawings, which are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific direction, be constructed and operated in a specific direction, and cannot be understood as limiting the specific protection scope of the present application. The terms "first" and "second" are only used for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. In the description of the present invention, the meaning of "multiple" is two or more, unless otherwise clearly and specifically defined.

[0046] The present invention relates to a method for compiling a joint scheduling plan for raw water in a multi-water source city, the method comprising the following steps:

[0047] S1, comprehensively consider the water supply shortage and comprehensive benefits of the study area to construct the objective function of the scheduling model;

[0048] S2, determining the constraints of the dispatching model by considering the parameters of regional reservoirs and water pipelines;

[0049] S3, the differential optimization algorithm and the stepwise optimization algorithm are combined to solve the scheduling model, obtain the optimal raw water supply flow process of the raw water supply pipeline in each time period, and then formulate a joint scheduling plan for raw water in multiple water source cities.

[0050] Specifically, the scheduling model objective function constructed in the above step S1 includes a first objective function, which is as follows:

[0051]

[0052] In the formula, OBJ1 represents the first objective function value, V represents the overall water supply deficit of the region, n and m represent the number of reservoirs and water plants respectively, i and j represent reservoirs and water plants respectively, t represents the month, and S i,t represents the water supply of the i-th reservoir in the t-th month, D j,t Represents the water demand of the j-th water plant in the t-th month.

[0053] Furthermore, the scheduling model objective function constructed in the above step S1 also includes a second objective function, which is as follows:

[0054]

[0055] In the formula, OBJ2 represents the second objective function value, that is, the maximum comprehensive benefit, B w,i,t represents the water supply benefit of the i-th reservoir in the t-th month, n represents the number of reservoirs, B g,i,t represents the power generation benefit of the i-th reservoir in the t-th month.

[0056] In this embodiment, the water supply benefit B of the above-mentioned reservoir in month t is w,i,t The formula is as follows:

[0057] B w,i,t =(Q w,i,t ×T w,i,t ×24×3600×P w ) / 10000

[0058] In the formula, Q w,i,t represents the average water supply flow of the i-th reservoir in the t-th month, m 3 / s,T w,i,t represents the number of days the i-th reservoir supplies water in the t-th month, P w Indicates the price of raw water;

[0059] The power generation benefit of the i-th reservoir in the t-th month is B g,i,t The formula is as follows:

[0060] B g,i,t =(Ng,i,t ×T g,i,t ×24×P g ) / 10000

[0061] =(K g,i ×H g,i,t ×Q g,i,t ×T g,i,t ×24×P g ) / 10000

[0062] Where N g,i,t represents the average output of the i-th reservoir in the t-th month, K g,i is the output coefficient, H g,i,t represents the average head difference of the ith reservoir in the tth month, Q g,i,t represents the average power generation flow of the i-th reservoir in the t-th month, T g,i,t represents the number of days when the ith reservoir generates electricity in the tth month, P g Indicates the electricity price.

[0063] In this embodiment, the scheduling model constraints determined in the above step S2 include:

[0064] Pipeline water delivery capacity constraints: SD i,j,t ≤SD i,j,t,max

[0065]

[0066] Reservoir water level constraint: Z i,t,min ≤Z i,t ≤Z i,t,max

[0067] Reservoir flow constraint: Q i,t,min ≤Q i,t ≤Q i,t,max

[0068] Reservoir output constraint: N i,t,min ≤N i,t ≤N i,t,max

[0069] Reservoir water level fluctuation constraint: ΔZ i,t,min ≤ΔZ i,t ≤ΔZ i,t,max

[0070] Variable non-negative constraint: SD i,j,t ≥0,D j,t ≥0,S i,t ≥0

[0071] In the above formula, SD i,j,t represents the upper limit of water flow from the i-th reservoir to the j-th water plant in the t-th month, SD i,t,maxThe upper limit of water flow from the i-th reservoir to all water plants in the t-th month, Z i,t,min and Z i,t,max They represent the lower and upper limits of the water level constraints of the i-th reservoir in the t-th month, Q i,t,min and Q i,t,max They represent the lower and upper limits of the outflow of the i-th reservoir in the t-th month, respectively. i,t,min and N i,t,max They represent the minimum and maximum guaranteed output of the i-th reservoir in the t-th month, ΔZ i,t,min and ΔZ i,t,max Represents the lower and upper limits of the water level variation of the ith reservoir in the tth month.

[0072] Since the optimal scheduling of raw water in cities with multiple water sources is a multi-dimensional and multi-constrained nonlinear optimization problem, it may be impossible to obtain the optimal solution of the optimization model through the algorithm. Therefore, it is necessary to propose a reasonable constraint correction strategy to convert the infeasible solution into a feasible solution to ensure the normal calculation of the algorithm. Specifically, after executing the combined differential optimization algorithm and the stepwise optimization algorithm described in step S3 to solve the scheduling model, it also includes the minimum downstream flow Q of the reservoir. i,t,min The steps for correction specifically include: when the water supply mode is water extraction from the reservoir, if the upstream water can meet the water supply needs, no correction is made; if the upstream water cannot meet the water supply needs, the minimum downstream flow of the reservoir is updated to the sum of the upstream water and the water supply flow; when the water supply mode is water extraction from the reservoir, if the upstream water can meet the water supply needs, no correction is made; if the upstream water cannot meet the water supply needs, the minimum downstream flow of the reservoir is updated to the water supply flow.

[0073] After executing the combined differential optimization algorithm and the stepwise optimization algorithm in step S3 to solve the scheduling model, it also includes calculating the water flow rate SD of the water supply pipeline. i,j,t,max The steps for making corrections specifically include: when the water supply mode is water extraction from the reservoir, determine whether the reservoir capacity meets the water extraction conditions. If so, no correction is required. If not, the water supply of the water supply pipeline of the reservoir needs to be reduced in proportion to the water supply capacity of the project to balance the water volume. When the water supply mode is water extraction from the reservoir, determine whether the water supply flow is greater than the downstream flow of the reservoir. If so, the water delivery flow of the water supply pipeline of the reservoir needs to be reduced in proportion to the water supply capacity of the project to balance the water volume. If not, no correction is required.

[0074] In this embodiment, before executing step S1, a step of generalizing the topological relationship of the water supply system is also included. First, it is necessary to establish a multi-water source-multi-water plant relationship matrix based on the interconnection of multiple water sources and multiple water plants; secondly, based on this relationship matrix, a network configuration relationship table of the complex water supply system is compiled. In addition, it is necessary to formulate reasonable traversal rules to traverse all entities of the urban raw water joint optimization scheduling model. Based on the above-mentioned water supply system configuration relationship table and scheduling model entity traversal rules, a multi-water source urban raw water joint optimization scheduling model that considers the comprehensive benefits of the region is constructed.

[0075] The urban raw water system is a complex system composed of multiple entities such as reservoirs, rivers, pipelines, and water-using units. Among them, the reservoir is a dispatching entity with a storage function composed of inherent attributes such as characteristic curves (e.g., water level-reservoir capacity curves) and characteristic parameters (e.g., normal water storage level); rivers and pipelines are key carriers for raw water delivery, which follow the principle of water balance (outflow is the difference between inflow and water diversion); water-using units are all water users in the region; key control sections include: the inflow section of the catchment area, the section where the reservoir is located, the water-using unit water intake section, the head and tail sections of the pipeline water delivery, etc. According to the above rules, a topological relationship diagram of the regional multi-water source joint optimization dispatching model is drawn.

[0076] It is very important to clarify the one-to-one correspondence between reservoir water supply and water plant water intake, and use a 0-1 matrix to indicate whether reservoir A can supply water to water plant B (if A supplies water to B, the value is 1, otherwise it is 0), and the network configuration relationship table of the complex water supply system based on the one-to-one correspondence matrix of multiple water sources and multiple water plants is very important, see the table below. The relationship matrix composed of reservoirs and water plants is as follows:

[0077]

[0078] In the matrix, R1~R m represents the water supply reservoir, m is the number of reservoirs; U1~U n Water intake plant, n is the number of water plants.

[0079] Table 1 Network configuration relationship table of complex water supply system

[0080]

[0081] In addition, the steps for traversing the entities of the urban raw water joint optimization scheduling model are as follows:

[0082] 1) Traverse all entity types contained in the regional optimization scheduling model to obtain all river / pipeline entities and store them one by one;

[0083] 2) Analyze and judge each river channel / pipeline, determine the traversal order based on the upstream and downstream relationship of the river channels corresponding to the front and rear sections corresponding to its location, and store them at the same time; during the traversal process, the sections need to be judged and follow the following rules: ① If the current section is the head and tail section or other important control section, the current section needs to be stored and automatically jump to the next section; ② If the current section is the inflow section of the catchment area, the reservoir entity section or the water intake section of the water unit, it is necessary to associate the entity type corresponding to the above position, store the current section and its corresponding entity type and jump to the next section.

[0084] 3) Determine whether all river / pipeline entities have been traversed. If all rivers / pipelines have been traversed, end the compilation.

[0085] In order to further illustrate the method of the present invention, the raw water dispatching in the central urban area of ​​a certain city is taken as an example for detailed description.

[0086] The city is a city with relatively scarce water resources, with per capita water resources of only 1285m 3 , which is about 60% of the average level in China. In this case, five large reservoirs are considered as the main water supply sources for the central urban area of ​​the city: Baixi, Zhougongzhai, Tingxia, Hengshan, Jiaokou Reservoirs, three medium-sized reservoirs: Hengxi, Sanxipu and Xixia Reservoirs, and one out-of-region reservoir: Qincun Reservoir.

[0087] In view of the complex network of the city's raw water system and the interconnection of multiple water sources and water plants, a multi-reservoir-multi-water plant relationship matrix is ​​constructed, among which the city's West Line Project reservoir-water plant relationship matrix is ​​shown below:

[0088]

[0089] The relationship matrix between the reservoir and water plant of the city's East Line Project is as follows:

[0090]

[0091] A network configuration table is constructed to describe the complex water supply system with multiple water sources, multiple pipelines, and multiple water plants interconnected, as shown in Table 2.

[0092] Table 2 Network configuration of a city's water supply system

[0093]

[0094]

[0095] In addition, if attached Figure 1As shown in the figure, the topological relationship of the eastern part of the city is taken as an example to compile the traversal process of the topological relationship, where R represents the reservoir entity, S represents the river or pipeline entity, U represents the water unit entity; P represents the cross-section element (belonging to the river or pipeline entity). The city is mainly composed of two rivers, and the order of the traversal domain topological structure is divided into the following two cases, as shown in the attached figure. Figure 2 As shown:

[0096] (1) Starting from the river channel of S1, traverse: S1→R1→S2→S3→U1→S4→S5→R2→S6→S7→S9→U2→S10→(P1)→S11→U3→S12→S13→U4→S14→(P2);

[0097] (2) Start traversing from the river channel belonging to S5: S5→R2→S6→S7→S9→U2→S10→S1→R1→S2→S3→U1→S4→(P1)→S11→U3→S12→S13→U4→S14→(P2).

[0098] Based on the above-mentioned water supply system configuration relationship table and scheduling model entity traversal rules, the implementation steps are as follows: Taking the water inflow process of the reservoir and the water demand process of the water plant in 2018 as an example, the scheduling process of the actual plan and the optimized plan is compared, and the set goals are compared to test the effectiveness of the optimization model. The actual water supply pattern of the city in 2018 is: 9 reservoirs jointly supply water to 8 water plants. The 9 reservoirs are 7 reservoirs under the water supply group (the "eight major reservoirs" do not include the Xixia Reservoir) and the Xixi (Huangtan) cascade. The 8 water plants are 4 water plants in the central urban area (the "five major water plants" in the central urban area do not include the Taoyuan Water Plant) and the local township water plants, Fenghua Water Plant, Xikou Water Plant and Zhenhai Refining and Chemical Water Plant. During the 2018 scheduling cycle, based on the water inflow process of each reservoir and the water demand process of each water plant, a regional multi-water source joint multi-objective urban raw water optimization scheduling model is established. Compare the actual scheduling and optimized scheduling process, see the attached Figures 3 to 9 . In the study, it is assumed that the reservoir's water level at the end of the month is the water level at the beginning of the next month. Therefore, the one-year study period includes 13 dispatching points. In addition, it should be noted that Zhougongzhai Reservoir and Jiaokou Reservoir, as cascade reservoirs, jointly supply water to Maojiaping Water Plant, and other reservoirs do not participate in the water supply dispatching of Maojiaping Water Plant. The results show that except for the optimized year-end water level of Baixi Reservoir, which is lower than the actual end-of-year water level, the optimized year-end water levels of the other reservoirs are higher than the actual year-end water level, and the monthly dispatching process is basically similar. Among them, Hengshan Reservoir, Hengxi Reservoir, Sanxipu Reservoir, Tingxia Reservoir, and Zhougongzhai Reservoir all operate at a relatively high water level. The comparison of the actual dispatching results and the dispatching results of the optimization model in 2018 is shown in Table 3.

[0099] Table 3 Summary and comparison of actual scheduling and optimization model scheduling effects in 2018

[0100] Scheduling Target Assessment criteria Practical solution Optimization plan Water supply shortage Water plant water demand - reservoir water supply 0 0 Comprehensive Benefits Water supply benefit + power generation benefit (100 million yuan) 2.88+0.39＝3.28 3.18+0.40＝3.58

[0101] The results show that compared with the actual dispatch (water supply guarantee rate 100%) in 2018, a normal year (the spatial distribution of rainfall in the city was uneven, with more rainfall in the north and less rainfall in the south), the regional multi-water source combined multi-objective annual optimization dispatch model increased the year-end water level of multiple reservoirs, with a total storage capacity of 68.11 million cubic meters, and increased power generation by 2 million kWh while ensuring water supply. It has a good comprehensive effect.

[0102] Although embodiments of the present invention have been shown and described, those skilled in the art will appreciate that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the claims and their equivalents.

Claims

1. A method for compiling a joint scheduling plan for raw water in a multi-water source city, characterized by: include, S1, comprehensively consider the water supply shortage and comprehensive benefits of the study area to construct the objective function of the scheduling model; S2, determining the constraints of the dispatching model by considering the parameters of regional reservoirs and water pipelines; S3, the differential optimization algorithm and the stepwise optimization algorithm are combined to solve the scheduling model, obtain the optimal raw water supply flow process of the raw water supply pipeline in each time period, and then formulate a joint scheduling plan for raw water in multiple water source cities.

2. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 1, characterized in that: The scheduling model objective function constructed in step S1 includes a first objective function, which is as follows: In the formula, OBJ1 represents the first objective function value, V represents the overall water supply deficit of the region, n and m represent the number of reservoirs and water plants respectively, i and j represent reservoirs and water plants respectively, t represents the month, and S i,t represents the water supply of the i-th reservoir in the t-th month, D j,t Represents the water demand of the j-th water plant in the t-th month.

3. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 2, characterized in that: The scheduling model objective function constructed in step S1 includes a second objective function, which is as follows: In the formula, OBJ2 represents the second objective function value, that is, the maximum comprehensive benefit, B w,i,t represents the water supply benefit of the i-th reservoir in the t-th month, n represents the number of reservoirs, B g,i,t represents the power generation benefit of the i-th reservoir in the t-th month.

4. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 3 is characterized by: The water supply benefit B of the i-th reservoir in month t w,i,t The formula is as follows: B w,i,t =(Q w,i,t ×T w,i,t ×24×3600×P w ) / 10000 In the formula, Q w,i,t represents the average water supply flow of the i-th reservoir in the t-th month, m 3 / s,T w,i,t represents the number of days the i-th reservoir supplies water in the t-th month, P w Indicates the raw water price.

5. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 4, characterized in that: The power generation benefit of the i-th reservoir in the t-th month is B g,i,t The formula is as follows: B g,i,t =(N g,i,t ×T g,i,t ×24×P g ) / 10000=(K g,i ×H g,i,t ×Q g,i,t ×T g,i,t ×24×P g ) / 10000 Where N g,i,t represents the average output of the i-th reservoir in the t-th month, K g,i is the output coefficient, H g,i,t represents the average head difference of the ith reservoir in the tth month, Q g,i,t represents the average power generation flow of the i-th reservoir in the t-th month, T g,i,t represents the number of days when the ith reservoir generates electricity in the tth month, P g Indicates the electricity price.

6. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 1, characterized in that: The scheduling model constraints determined in step S2 include: Pipeline water delivery capacity constraints: SD i,j,t ≤SD i,j,t,max Reservoir water level constraint: Z i,t,min ≤Z i,t ≤Z i,t,max Reservoir flow constraint: Q i,t,min ≤Q i,t ≤Q i,t,max Reservoir output constraint: N i,t,min ≤N i,t ≤N i,t,max Reservoir water level fluctuation constraint: ΔZ i,t,min ≤ΔZ i,t ≤ΔZ i,t,max Variable non-negative constraint: SD i,j,t ≥0,D j,t ≥0,S i,t ≥0 In the above formula, SD i,j,t represents the upper limit of water flow from the i-th reservoir to the j-th water plant in the t-th month, SD i,t,max The upper limit of water flow from the i-th reservoir to all water plants in the t-th month, Z i,t,min and Z i,t,max They represent the lower and upper limits of the water level constraints of the i-th reservoir in the t-th month, Q i,t,min and Q i,t,max They represent the lower and upper limits of the outflow of the i-th reservoir in the t-th month, respectively. i,t,min and N i,t,max They represent the minimum and maximum guaranteed output of the i-th reservoir in the t-th month, ΔZ i,t,min and ΔZ i,t,max Represents the lower and upper limits of the water level variation of the ith reservoir in the tth month.

7. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 1, characterized in that: After executing the joint differential optimization algorithm and the stepwise optimization algorithm described in step S3 to solve the scheduling model, it also includes the step of correcting the minimum downstream flow of the reservoir, specifically including: when the water supply mode is water intake from the reservoir, if the upstream water can meet the water supply needs, no correction is made; if the upstream water cannot meet the water supply needs, the minimum downstream flow of the reservoir is updated to the sum of the upstream water and the water supply flow; when the water supply mode is water intake from the reservoir, if the upstream water can meet the water supply needs, no correction is made; if the upstream water cannot meet the water supply needs, the minimum downstream flow of the reservoir is updated to the water supply flow.

8. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 7, characterized in that: After executing the joint differential optimization algorithm and the stepwise optimization algorithm described in step S3 to solve the scheduling model, it also includes the step of correcting the water flow of the water supply pipeline, specifically including: when the water supply mode is water intake from the reservoir, determine whether the reservoir capacity meets the water intake conditions, if so, no correction is required, if not, the water supply of the water supply pipeline of this reservoir needs to be reduced according to the proportion of the project water supply capacity to balance the water volume, when the water supply mode is water intake from the reservoir, determine whether the water supply flow is greater than the reservoir discharge flow, if so, reduce the water flow of the water supply pipeline of this reservoir according to the proportion of the project water supply capacity to balance the water volume, if not, no correction is made.

9. The method for compiling a joint scheduling plan for raw water in a multi-water source city according to claim 1, characterized in that: Before executing step S1, a step of generalizing the topological relationship of the water supply system is also included.

Images