A reservoir water supply decision optimization method, device, medium and product

Abstract

The application discloses a reservoir water supply decision optimization method, device, medium and product, relates to the technical field of water resource optimization management, and the reservoir water supply decision optimization method comprises the following steps: obtaining reservoir basic data; determining a preliminary water supply decision of the reservoir based on the reservoir basic data; the preliminary water supply decision comprises a preliminary water supply amount and an initial final reservoir capacity of each time period in a scheduling period; determining whether the initial final reservoir capacity of the reservoir in each time period meets the reservoir capacity constraint; if yes, determining the preliminary water supply decision of the reservoir as the optimal water supply decision of the reservoir; if not, performing multi-round multi-section interval optimization based on the preliminary water supply decision by using the reservoir capacity constraint, obtaining the optimal water supply amount and the optimal final reservoir capacity of all time periods in the scheduling period, and determining the optimal water supply amount and the optimal final reservoir capacity of each time period in the scheduling period as the optimal water supply decision of the reservoir. The application realizes efficient and accurate optimization of long-term water supply decision of the reservoir.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of water resource optimal management, in particular to a reservoir water supply decision optimization method, device, medium and product. BACKGROUND

[0002] As the core infrastructure of water resource comprehensive regulation, reservoirs play a key role in water resource allocation. According to statistics, about 30% of the reservoirs in the world are built primarily for water supply. Due to the spatio-temporal heterogeneity of natural runoff, how to fully exert the regulation and storage capacity of reservoirs to maximize the efficiency of regional water supply has always been a core scientific problem in the field of water resource system optimization.

[0003] The traditional dynamic programming method needs to discretize the reservoir capacity and water supply when solving the long-term optimal scheduling problem of reservoirs. This method not only has the problem of exponential growth of computational complexity with the increase of discretization accuracy, but also has the limitation of not being able to accurately approximate the optimal process when the discretization accuracy is low. In contrast, the reservoir optimization scheduling method based on hedging rules can directly obtain the optimal decision through analytical derivation, which has significant advantages in computational efficiency and theoretical optimality. However, the related reservoir scheduling based on analytical optimization is mainly limited to single-stage or two-stage water supply decision, which is difficult to extend to continuous multi-stage long-term optimization scenarios. Secondly, in long-term water supply scheduling, the constraint conditions increase step by step. How to construct a long-term optimal water supply process based on the convex programming optimality condition, and how to determine the optimal correction order when the multi-period decision variable violates the constraint, these two key problems have not been effectively solved, which restricts the application of this method in actual scheduling engineering. SUMMARY

[0004] The purpose of the present application is to provide a reservoir water supply decision optimization method, device, medium and product to realize efficient and accurate optimization of long-term water supply decision of reservoirs.

[0005] To achieve the above purpose, the present application provides the following solutions.

[0006] In a first aspect, the present application provides a reservoir water supply decision optimization method, comprising:

[0007] obtaining reservoir basic data;

[0008] determining a preliminary water supply decision of the reservoir based on the reservoir basic data; the preliminary water supply decision includes preliminary water supply quantity and initial final reservoir capacity of each period in the scheduling period;

[0009] judging whether the initial final reservoir capacity of the reservoir in each period meets the reservoir capacity constraint to obtain a first judgment result;

[0010] if the first judgment result is yes, determining the preliminary water supply decision of the reservoir as the optimal water supply decision of the reservoir;

[0011] If the first determination result is no, then a multi-iteration multi-section interval optimization is performed based on the preliminary water supply decision using the reservoir capacity constraint to obtain the optimized water supply amount and the optimized end reservoir capacity of all time periods in the dispatch period, and the optimized water supply amount and the optimized end reservoir capacity of each time period in the dispatch period are determined as the optimal water supply decision of the reservoir.

[0012] In an embodiment, the reservoir basic data includes: an initial reservoir capacity in the dispatch period, upstream inflow amounts of each time period in the dispatch period, a number of time periods in the dispatch period, and an end-of-period reservoir capacity in the dispatch period.

[0013] In an embodiment, determining a preliminary water supply decision of the reservoir based on the reservoir basic data includes:

[0014] determining a cumulative available water supply amount of the reservoir in the dispatch period according to the initial reservoir capacity in the dispatch period, the upstream inflow amounts of each time period in the dispatch period, and the end-of-period reservoir capacity in the dispatch period;

[0015] determining a preliminary water supply amount of each time period in the dispatch period according to the cumulative available water supply amount of the reservoir in the dispatch period and the number of time periods in the dispatch period;

[0016] determining a to-be-calculated time period in the dispatch period;

[0017] determining an initial end reservoir capacity of the to-be-calculated time period according to the initial reservoir capacity in the dispatch period, upstream inflow amounts of all time periods from a starting time period in the dispatch period to the to-be-calculated time period, and preliminary water supply amounts of all time periods from the starting time period in the dispatch period to the to-be-calculated time period.

[0018] In an embodiment, the reservoir capacity constraint includes:

[0019] ;

[0020] ;

[0021] wherein, is an end reservoir capacity of time period t; is a minimum reservoir capacity; is a maximum reservoir capacity.

[0022] In an embodiment, the multi-iteration multi-section interval optimization process of any current iteration includes:

[0023] determining all violation time periods from the initial round to the last round as a violation time period set of the current round, and determining optimized end storages of all violation time periods from the initial round to the last round as optimized end storages of each violation time period in the violation time period set of the current round; when the current round is the initial round, the last round is a time period in which the initial end storage does not satisfy the storage constraint; and the optimized end storages of all violation time periods in the last round are obtained by adjusting the storages in the multiple segmented intervals of the round using the storage constraint;

[0024] dividing all time periods in the scheduling period based on the violation time period set of the current round to obtain multiple segmented intervals of the current round; the first time period and the last time period in each segmented interval of the current round correspond to a violation time period in the violation time period set of the current round, a starting time period and an ending time period in the scheduling period;

[0025] determining any segmented interval of the current round as a current segmented interval of the current round;

[0026] determining water supply amounts and end storages of each time period in the current segmented interval of the current round based on the optimized end storage of the last time period of the first time period in the current segmented interval of the current round, the optimized end storage of the last time period in the current segmented interval of the current round, and upstream inflow amounts of each time period in the current segmented interval of the current round; when the current segmented interval of the current round is the first segmented interval of the current round, the optimized end storage of the last time period of the first time period in the current segmented interval of the current round is the initial storage in the scheduling period; and when the current segmented interval of the current round is the last segmented interval of the current round, the optimized end storage of the last time period in the current segmented interval of the current round is the final storage in the scheduling period;

[0027] determining whether the end storages of all time periods in the scheduling period of the current round satisfy the storage constraint to obtain a second determination result;

[0028] if the second determination result is yes, determining the water supply amounts and the end storages of all time periods in the scheduling period of the current round as optimized water supply amounts and optimized end storages of all time periods in the scheduling period, and obtaining the optimal water supply decision of the reservoir;

[0029] if the second determination result is no, determining violation time periods of the current round based on all time periods in the current round that do not satisfy the storage constraint, adjusting the end storages of each violation time period of the current round to obtain optimized end storages of each violation time period of the current round, determining all violation time periods from the initial round to the current round as a violation time period set of the next round, and performing optimization of the next round until the end storages of all time periods in the scheduling period satisfy the storage constraint to obtain the optimal water supply decision of the reservoir.

[0030] In an embodiment, the water supply amount and the end storage capacity of each period in the current segment interval of the current round are determined based on the optimized end storage capacity of the previous period of the first period in the current segment interval of the current round, the optimized end storage capacity of the last period in the current segment interval of the current round, and the upstream inflow of each period in the current segment interval of the current round, including:

[0031] The cumulative water supply amount of the reservoir in the current segment interval of the current round is determined according to the optimized end storage capacity of the previous period of the first period in the current segment interval of the current round, the optimized end storage capacity of the last period in the current segment interval of the current round, and the upstream inflow of each period in the current segment interval of the current round.

[0032] The water supply amount of each period in the current segment interval of the current round is determined according to the cumulative water supply amount of the reservoir in the current segment interval of the current round and the number of periods in the current segment interval of the current round.

[0033] Any period in the current segment interval of the current round is determined as a target period in the current segment interval of the current round.

[0034] The end storage capacity of the target period in the current segment interval of the current round is determined according to the optimized end storage capacity of the previous period of the first period in the current segment interval of the current round, the upstream inflow of all periods from the first period to the target period in the current segment interval of the current round, and the water supply amount of all periods from the first period to the target period in the current segment interval of the current round, thereby obtaining the end storage capacity of each period in the current segment interval of the current round.

[0035] In an embodiment, the end storage capacity of each violation period in the current round is adjusted to obtain the optimized end storage capacity of each violation period in the current round, including:

[0036] When the end storage capacity of each violation period in the current round is less than the minimum storage capacity, the end storage capacity of each violation period in the current round is adjusted to the minimum storage capacity to obtain the optimized end storage capacity of each violation period in the current round.

[0037] When the end storage capacity of each violation period in the current round is greater than the maximum storage capacity, the end storage capacity of each violation period in the current round is adjusted to the maximum storage capacity to obtain the optimized end storage capacity of each violation period in the current round.

[0038] In a second aspect, the present application provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the reservoir water supply decision optimization method described above.

[0039] In a third aspect, the present application provides a computer readable storage medium, having stored thereon a computer program, which, when executed by a processor, implements the reservoir water supply decision optimization method described above.

[0040] In a fourth aspect, the present application provides a computer program product, comprising a computer program, which, when executed by a processor, implements the reservoir water supply decision optimization method described above.

[0041] According to the specific embodiments provided by the present application, the following technical effects are disclosed:

[0042] The present application discloses a reservoir water supply decision optimization method, device, medium and product. When there is a time period that does not meet the reservoir capacity constraint in the preliminary water supply decision, the reservoir capacity constraint is used for multi-round multi-section interval optimization, that is, the division and optimization of the section interval are performed based on the reservoir capacity constraint in each round. After multi-round multi-section interval optimization, the final reservoir capacity of each time period meets the reservoir capacity constraint. At this time, the optimized water supply quantity and the optimized final reservoir capacity of the corresponding all time periods are determined as the optimal water supply decision of the reservoir, and the reservoir capacity constraint is optimized through multi-round multi-section interval optimization to realize efficient and accurate optimization of the long-term water supply decision of the reservoir. BRIEF DESCRIPTION OF DRAWINGS

[0043] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the related art, the drawings needed in the embodiments will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor on the basis of these drawings.

[0044] Figure 1 The reservoir water supply decision optimization method flowchart provided by an embodiment of the present application.

[0045] Figure 2 The upstream inflow quantity diagram of each time period in the dispatching period.

[0046] Figure 3 The final reservoir capacity diagram of each time period in the dispatching period optimized in each round.

[0047] Figure 4 The water supply quantity diagram of each time period in the dispatching period optimized in each round.

[0048] Figure 5 The structure diagram of a computer device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0049] With reference to the accompanying drawings, the technical solutions in the embodiments of the present application will be described clearly and completely. Obviously, the described embodiments are only a part of the embodiments of the present application, but not all the embodiments of the present application. Based on the embodiments in the present application, all the other embodiments obtained by those ordinarily skilled in the art without creative effort belong to the scope of the present application.

[0050] The present application aims to provide a reservoir water supply decision optimization method, device, medium and product, and aims to realize efficient and accurate optimization of long-term reservoir water supply decision.

[0051] In order to make the above-mentioned purposes, features and advantages of the present application more obvious and easy to understand, the present application will be further described in detail below with reference to the drawings and specific embodiments.

[0052] In an exemplary embodiment, as shown in Figure 1 A reservoir water supply decision optimization method is provided, comprising the following steps.

[0053] Step 1: Obtain reservoir basic data.

[0054] As an optional implementation, the reservoir basic data includes initial reservoir capacity in the scheduling period, upstream inflow in each period in the scheduling period, number of periods in the scheduling period, and final reservoir capacity in the scheduling period.

[0055] Step 2: Determine the preliminary water supply decision of the reservoir based on the reservoir basic data; the preliminary water supply decision includes preliminary water supply in each period in the scheduling period and initial final reservoir capacity.

[0056] As an optional implementation, step 2 includes steps 21-24.

[0057] Step 21: Determine the cumulative available water of the reservoir in the scheduling period according to the initial reservoir capacity in the scheduling period, the upstream inflow in each period in the scheduling period, and the final reservoir capacity in the scheduling period.

[0058] Specifically, since the initial reservoir capacity in the scheduling period, the upstream inflow in each period in the scheduling period, the number of periods in the scheduling period, and the final reservoir capacity in the scheduling period are all known quantities, the calculation formula of the cumulative available water of the reservoir in the scheduling period is:

[0059] (1)

[0060] wherein, is the cumulative available water of the reservoir in the scheduling period; is the number of periods in the scheduling period; is the initial reservoir capacity in the scheduling period; is the upstream inflow in period t in the scheduling period; to the end of the dispatch period.

[0061] Step 22: determining the preliminary water supply amount of each time period in the dispatch period according to the cumulative water supply amount of the reservoir in the dispatch period and the number of time periods in the dispatch period.

[0062] Specifically, in general, in Step 22, when determining the preliminary water supply amount of each time period in the dispatch period according to the cumulative water supply amount of the reservoir in the dispatch period and the number of time periods in the dispatch period, the conditions to be met are shown in Equations (2) and (3). Equations (2) and (3) indicate that, when the reservoir is only subject to the water balance constraint, the optimal water supply analytical solution should satisfy the marginal benefit of water supply of each time period being equal.

[0063] (2)

[0064] (3)

[0065] wherein, is the water supply benefit of the reservoir in the dispatch period under the optimal water supply amount of the time period t, , ; is the preliminary water supply amount of the time period t in the dispatch period.

[0066] Particularly, if the water supply benefit function form does not change with time, and the water demand of each time period is the same, the calculation formula of the preliminary water supply amount is:

[0067] (4)

[0068] Step 23: determining any time period in the dispatch period as a to-be-calculated time period.

[0069] Step 24: determining the initial end storage capacity of the to-be-calculated time period according to the initial storage capacity in the dispatch period, the upstream inflow amount from the initial time period in the dispatch period to all time periods in the to-be-calculated time period, and the preliminary water supply amount from the initial time period in the dispatch period to all time periods in the to-be-calculated time period.

[0070] Specifically, the calculation formula of the initial end storage capacity of any time period in the dispatch period is:

[0071] (5)

[0072] wherein, is the initial end storage capacity of the time period m in the dispatch period.

[0073] In fact, the preliminary water supply amount of each time period in the dispatch period is obtained by averaging the cumulative water supply amount in the dispatch period, and the theoretical process of derivation is as follows.

[0074] ​The optimization method of the reservoir water supply decision of the application optimizes the cumulative water supply benefit of the entire scheduling period to be maximum, and the objective function is:

[0075] (6)

[0076] wherein, is the water demand of the period t in the scheduling period, and the water demand of each period is the same.

[0077] Step 3: Determine whether the initial final reservoir capacity of the reservoir in each period meets the reservoir capacity constraint to obtain a first determination result.

[0078] As an optional implementation, the reservoir capacity constraint comprises:

[0079] (7)

[0080] (8)

[0081] wherein, is the final reservoir capacity of the period t; is the minimum reservoir capacity; is the maximum reservoir capacity. Formula (7) is the minimum reservoir capacity constraint, and formula (8) is the maximum reservoir capacity constraint.

[0082] Step 4: If the first determination result is yes, the preliminary water supply decision of the reservoir is determined as the optimal water supply decision of the reservoir.

[0083] Step 5: If the first determination result is no, the multi-division interval optimization of multiple rounds is performed based on the preliminary water supply decision by using the reservoir capacity constraint, the optimized water supply amount and the optimized final reservoir capacity of all periods in the scheduling period are obtained, and the optimized water supply amount and the optimized final reservoir capacity of each period in the scheduling period are determined as the optimal water supply decision of the reservoir.

[0084] As an optional implementation, in step 5, the multi-division interval optimization process of any current round comprises steps 51-57.

[0085] Step 51: Determine all violation periods from the initial round to the last round as the violation period set of the current round, and determine the optimized final reservoir capacity of all violation periods from the initial round to the last round as the optimized final reservoir capacity of each violation period in the violation period set of the current round; when the current round is the initial round, the violation period of the last round is the period whose initial final reservoir capacity does not meet the reservoir capacity constraint; the optimized final reservoir capacity of all violation periods of the last round is obtained by adjusting the multi-division interval optimization of the round by using the reservoir capacity constraint.

[0086] Step 52: dividing all time periods in the scheduling period based on the violation time period set of the current round to obtain a plurality of segmented intervals of the current round; the first time period and the last time period in each segmented interval of the current round correspond to the violation time period in the violation time period set of the current round, the starting time period and the ending time period in the scheduling period.

[0087] Step 53: determining any segmented interval of the current round as the current segmented interval of the current round.

[0088] Step 54: determining the water supply amount and the final reservoir capacity of each time period in the current segmented interval of the current round based on the optimized final reservoir capacity of the last time period of the first time period in the current segmented interval of the current round, the optimized final reservoir capacity of the first time period in the current segmented interval of the current round and the upstream inflow of each time period in the current segmented interval of the current round; when the current segmented interval of the current round is the first segmented interval of the current round, the optimized final reservoir capacity of the last time period of the first time period in the current segmented interval of the current round is the initial reservoir capacity in the scheduling period; when the current segmented interval of the current round is the last segmented interval of the current round, the optimized final reservoir capacity of the first time period in the current segmented interval of the current round is the final reservoir capacity in the scheduling period.

[0089] As an optional implementation, step 54 includes steps 541-544.

[0090] Step 541: determining the cumulative water supply amount of the reservoir in the current segmented interval of the current round according to the optimized final reservoir capacity of the last time period of the first time period in the current segmented interval of the current round, the optimized final reservoir capacity of the first time period in the current segmented interval of the current round and the upstream inflow of each time period in the current segmented interval of the current round.

[0091] Step 542: determining the water supply amount of each time period in the current segmented interval of the current round according to the cumulative water supply amount of the reservoir in the current segmented interval of the current round and the number of time periods in the current segmented interval of the current round.

[0092] Step 543: determining any time period in the current segmented interval of the current round as the target time period in the current segmented interval of the current round.

[0093] Step 544: determining the final reservoir capacity of the target time period in the current segmented interval of the current round according to the optimized final reservoir capacity of the last time period of the first time period in the current segmented interval of the current round, the upstream inflow of all time periods from the first time period to the target time period in the current segmented interval of the current round and the water supply amount of all time periods from the first time period to the target time period in the current segmented interval of the current round, thereby obtaining the final reservoir capacity of each time period in the current segmented interval of the current round.

[0094] Step 55: judging whether the end storage capacity of all time periods in the scheduling period of the current round meets the storage capacity constraint, to obtain a second judgment result.

[0095] Step 56: if the second judgment result is yes, determining the water supply amount and the end storage capacity of all time periods in the scheduling period of the current round as the optimized water supply amount and the optimized end storage capacity of all time periods in the scheduling period, to obtain the optimal water supply decision of the reservoir.

[0096] Step 57: if the second judgment result is no, determining the violation time period of the current round based on all time periods of the current round that do not meet the storage capacity constraint, and adjusting the end storage capacity of each violation time period of the current round to obtain the optimized end storage capacity of each violation time period of the current round, determining all violation time periods from the initial round to the current round as the violation time period set of the next round, and performing optimization of the next round until the end storage capacity of all time periods in the scheduling period meets the storage capacity constraint, to obtain the optimal water supply decision of the reservoir.

[0097] Specifically, in step 57, determining the violation time period of the current round based on all time periods of the current round that do not meet the storage capacity constraint includes steps 571-572.

[0098] Step 571: in all time periods of the current round that do not meet the storage capacity constraint, determining all continuous time periods that do not meet the minimum storage capacity constraint as the minimum continuous time period of the current round, and determining the time period with the minimum end storage capacity in each minimum continuous time period of the current round as the violation time period of the current round.

[0099] Step 572: in all time periods of the current round that do not meet the storage capacity constraint, determining all continuous time periods that do not meet the maximum storage capacity constraint as the maximum continuous time period of the current round, and determining the time period with the maximum end storage capacity in each maximum continuous time period of the current round as the violation time period of the current round.

[0100] As an optional implementation, in step 57, adjusting the end storage capacity of each violation time period of the current round to obtain the optimized end storage capacity of each violation time period of the current round includes steps 571-572.

[0101] Step 571: when the end storage capacity of each violation time period of the current round is less than the minimum storage capacity, adjusting the end storage capacity of each violation time period of the current round to the minimum storage capacity to obtain the optimized end storage capacity of each violation time period of the current round.

[0102] Step 572: when the end storage capacity of each violation time period of the current round is greater than the maximum storage capacity, adjusting the end storage capacity of each violation time period of the current round to the maximum storage capacity to obtain the optimized end storage capacity of each violation time period of the current round.

[0103] Furthermore, taking a medium-sized reservoir in a certain river basin as an example, the scheduling period is from 1990 to 2020, the scheduling step is monthly, and the minimum reservoir capacity is... 4 million m 3 Maximum storage capacity 41 million m 3 The reservoir inflow process, i.e., the upstream inflow volume at each time period during the scheduling period, is as follows: Figure 2 As shown, the average annual water inflow is 2.25 million cubic meters. 3 The reservoir supplies water to downstream cities, with a water demand of 22 million cubic meters per hour. 3 The initial reservoir capacity in January 1990. 28.8 million m 3 The ending storage capacity in December 2020 28.8 million m 3 The water supply decision optimization method of this application for the reservoir includes the following process.

[0104] S1: Obtain basic data on medium-sized reservoirs in a certain watershed.

[0105] Basic data for a medium-sized reservoir in a certain river basin includes: an initial reservoir capacity of 28.8 million m³ during the scheduling period. 3 Upstream water inflow during each time period within the scheduling period; Number of time periods within the scheduling period. and the end-of-period storage capacity during the scheduling period 28.8 million m 3 .

[0106] S2: Using formulas (1), (4), and (5), the preliminary water supply decision for the reservoir is determined based on the reservoir's basic data. The preliminary water supply decision includes the preliminary water supply volume and initial final reservoir capacity for each period during the scheduling period. The preliminary water supply volume for each period is 18,752,600 m³. 3 .

[0107] S3: Determine whether the initial and final storage capacities of the reservoir in each time period meet the storage capacity constraints, and obtain the first judgment result.

[0108] S4: If the first judgment result is yes, then the preliminary water supply decision of the reservoir is determined as the optimal water supply decision of the reservoir.

[0109] S5: If the first judgment result is negative, then using the reservoir capacity constraint, multi-round multi-segment interval optimization is performed based on the preliminary water supply decision to obtain the optimized water supply volume and optimized final reservoir capacity for all time periods within the scheduling period. The optimized water supply volume and optimized final reservoir capacity for each time period within the scheduling period are then determined as the optimal water supply decision for the reservoir. S5 specifically includes S501-S510.

[0110] S501: determine the time period corresponding to the initial end reservoir capacity that does not satisfy the minimum reservoir capacity constraint in the reservoir capacity constraint as the minimum time period, and the initial end reservoir capacities corresponding to all minimum time periods constitute a minimum violated end reservoir capacity set : , is the initial end reservoir capacity of the time period t in the scheduling period.

[0111] In this embodiment, the initial end reservoir capacities of 22 time periods do not satisfy the minimum reservoir capacity constraint in the reservoir capacity constraint, that is, .

[0112] S502: determine the minimum violated end reservoir capacity set In the step S502, a plurality of continuous minimum time periods are determined as minimum continuous time periods, and the time period with the minimum end reservoir capacity in each minimum continuous time period is determined as the violated time period of the first round, and a minimum violated time period set of the first round is constituted .

[0113] In this embodiment, In the step S502, the time period 2 to the time period 5 are the minimum continuous time period 1, the time period with the minimum end reservoir capacity in the minimum continuous time period 1 is the time period 4, the time period 16 to the time period 18 are the minimum continuous time period 2, the time period with the minimum end reservoir capacity in the minimum continuous time period 2 is the time period 17, and the like, and the minimum violated time period set of the first round is obtained .

[0114] S503: determine the time period corresponding to the initial end reservoir capacity that does not satisfy the maximum reservoir capacity constraint in the reservoir capacity constraint as the maximum time period, and the initial end reservoir capacities corresponding to all maximum time periods constitute a maximum violated end reservoir capacity set : .

[0115] In this embodiment, the initial end reservoir capacities of 286 time periods do not satisfy the maximum reservoir capacity constraint in the reservoir capacity constraint, that is, .

[0116] S504: determine the maximum violated end reservoir capacity set In the step S504, a plurality of continuous maximum time periods are determined as maximum continuous time periods, and the time period with the maximum end reservoir capacity in each maximum continuous time period is determined as the violated time period of the first round, and a maximum violated time period set of the first round is constituted .

[0117] In this embodiment, In the step S504, the time period 8 to the time period 12 are the maximum continuous time period 1, the time period with the maximum end reservoir capacity in the maximum continuous time period 1 is the time period 9, the time period 21 to the time period 25 are the maximum continuous time period 2, the time period with the maximum end reservoir capacity in the maximum continuous time period 2 is the time period 22, and the like, and the maximum violated time period set of the first round is obtained .

[0118] S505: obtaining a minimum violation period set of the first round and a maximum violation period set of the first round The violation period set of the first round is composed of the minimum violation period set and the maximum violation period set of the first round, and the corresponding end storage capacity set of the first round is denoted as: .

[0119] In this embodiment, .

[0120] S506: adjusting the end storage capacity of each period in the violation period set of the first round to obtain the optimized end storage capacity of each period in the violation period set of the first round Specifically, if , let , if , let The end storage capacity of period 0 is the initial storage capacity , the end storage capacity of period T is the final storage capacity , and the end storage capacity of other periods is a variable to be optimized and solved, that is: .

[0121] S507: dividing all periods in the scheduling period based on each period in the violation period set of the first round to obtain a plurality of segmented intervals of the first round.

[0122] In this embodiment, the violation period set of the first round has 14 violation periods, and the optimized end storage capacity of the 14 violation periods is shown in Table 1. The first round has 15 segmented intervals, which are period 1~period 4, period 5~period 9, period 10~period 17, period 18~period 22, period 23~period 40, period 41~period 46, period 47~period 64, period 65~period 70, period 71~period 76, period 77~period 165, period 166~period 220, period 221~period 249, period 250~period 364, period 365~period 370 and period 371~period 372.

[0123] Table 1: Optimized end storage capacity table of 14 violation periods in the violation period set of the first round

[0124]

[0125] S508: determining any segmented interval of the first round as the current segmented interval of the first round.

[0126] S509: Based on the optimized final reservoir capacity of the previous time period within the first segment interval of the first round, the optimized final reservoir capacity of the last time period within the current segment interval of the first round, and the upstream inflow of each time period within the current segment interval of the first round, determine the water supply and final reservoir capacity for each time period within the current segment interval of the first round. For any segment interval of the first round... The formulas for calculating the water supply and final reservoir capacity for each time period include:

[0127] ;

[0128] ;

[0129] .

[0130] in, For the reservoir in the first round of segmented sections The cumulative available water volume within; The segmented interval for the first round Inner First Period The optimized storage capacity at the end of the previous time period; The last time period within the current segment interval of the first round Optimize the final storage capacity; The segmented interval for the first round Water supply during time period t; The segmented interval for the first round The number of time periods within a single period; The segmented interval for the first round The final storage capacity within time period t.

[0131] S510: Repeatedly execute similar operations from S501 to S509 until the final reservoir capacity of all time periods within the scheduling cycle meets the reservoir capacity constraint, thus obtaining the optimal water supply decision for the reservoir.

[0132] In this embodiment, after 6 iterations, the process of determining the final reservoir capacity and water supply for each time period within the optimized scheduling period in each round is as follows: Figure 3 and Figure 4 As shown. The optimal water supply process is ultimately obtained where water supply decisions in all time periods satisfy the reservoir capacity constraint. Figure 4 (Black line).

[0133] To further verify the implementation effect of this application, Table 2 compares its optimization results with those of the traditional dynamic programming method. The discrete storage capacity K set for dynamic programming is 100,000 m³. 3 and 10,000 m 3Two scenarios. As shown in Table 2, the method of the present application shows significant advantages in calculation efficiency and optimization effect. In terms of calculation efficiency, the calculation of water supply decision for 372 periods only takes 0.3s, and the calculation time of the dynamic programming method increases exponentially with the increase of the discrete storage capacity accuracy, reaching 100 to 10000 times of the method of the present application. In terms of water supply efficiency, the method of the present application obtains the theoretical optimal solution of the long-term water supply process through analytical derivation, and the average water supply efficiency is the largest, which is 0.8575; the dynamic programming method has similar results when the high discrete accuracy (K=10000 m 3 ) is used, but when the low discrete accuracy (K=100000 m 3 ) is used, the water supply efficiency is still 4% lower than the optimal value.

[0134] Table 2 Comparison table of calculation results of reservoir long-term optimization by different optimization methods

[0135]

[0136] From the above analysis, it can be seen that the reservoir water supply decision optimization method of the present application not only ensures the global optimality of the water supply process, but also significantly reduces the calculation complexity to the linear level compared with the traditional dynamic programming algorithm, greatly improving the efficiency of reservoir optimization solution.

[0137] In an exemplary embodiment, a computer device is provided, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor executing the computer program to implement the reservoir water supply decision optimization method.

[0138] In an exemplary embodiment, a computer readable storage medium is provided, which stores a computer program, and the computer program is executed by a processor to implement the reservoir water supply decision optimization method.

[0139] In an exemplary embodiment, the present application provides a computer program product comprising a computer program, which is executed by a processor to implement the above-mentioned reservoir water supply decision optimization method.

[0140] In an exemplary embodiment, a computer device is provided, which can be a server or a terminal, and its internal structure diagram can be as shown in Figure 5As shown in the figure. The computer device includes a processor, a memory, an Input / Output (I / O) interface and a communication interface. Among them, the processor, the memory and the input / output interface are connected through the system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capability. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and the external device. The communication interface of the computer device is used to communicate with the terminal outside through the network connection. The computer program is executed by the processor to implement a reservoir water supply decision optimization method.

[0141] Those skilled in the art can understand that, Figure 5 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0142] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of each method can be included. Any reference to memory, database or other medium used in each embodiment provided by the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (Read-Only Memory, ROM), magnetic tape, floppy disk, flash memory, optical storage, high-density embedded non-volatile memory, resistive memory (ReRAM), magnetoresistive random access memory (Magnetoresistive Random Access Memory, MRAM), ferroelectric memory (Ferroelectric Random Access Memory, FRAM), phase change memory (Phase Change Memory, PCM), graphene memory, etc. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, etc. As an illustration but not limitation, RAM can be in various forms, such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), etc.

[0143] The database involved in each embodiment provided by the present application can include at least one of a relational database and a non-relational database. The non-relational database can include a distributed database based on a blockchain, etc., without being limited thereto. The processor involved in each embodiment provided by the present application can be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without being limited thereto.

[0144] It should be noted that the user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of related data need to comply with relevant regulations.

[0145] Each technical feature of the above embodiments can be combined arbitrarily. In order to make the description simple, not all possible combinations of each technical feature in the above embodiments are described, but as long as the combination of these technical features does not exist contradictory, it should be considered as the scope of the present application.

[0146] The principles and implementations of the present application are described in detail with specific examples in this paper, and the above examples are only used to help understand the method of the present application and its core idea; at the same time, for those skilled in the art, according to the idea of the present application, the specific implementation and application range will be changed. Therefore, the content of the specification should not be understood as a limitation of the present application.

Claims

1. A method for optimizing reservoir water supply decisions, characterized in that, The reservoir water supply decision optimization method includes: Obtain basic data on the reservoir; The preliminary water supply decision for the reservoir is determined based on the basic data of the reservoir; the preliminary water supply decision includes the preliminary water supply volume and the initial final reservoir capacity for each period during the scheduling period; Determine whether the initial and final storage capacities of the reservoir in each time period meet the storage capacity constraints, and obtain the first judgment result; If the first judgment result is yes, then the preliminary water supply decision of the reservoir is determined as the optimal water supply decision of the reservoir; If the first judgment result is negative, then the reservoir capacity constraint is used to perform multi-round multi-segment interval optimization based on the preliminary water supply decision to obtain the optimized water supply volume and optimized final reservoir capacity for all time periods within the scheduling period, and the optimized water supply volume and optimized final reservoir capacity for each time period within the scheduling period are determined as the optimal water supply decision for the reservoir. The storage capacity constraint includes: ； ； in, The final storage capacity at time period t; Minimum storage capacity; Maximum storage capacity; The multi-segment interval optimization process in any current round includes: The set of all violation periods from the initial round to the previous round is defined as the set of violation periods for the current round. The optimized final storage capacity of all violation periods from the initial round to the previous round is defined as the optimized final storage capacity of each violation period in the set of violation periods for the current round. When the current round is the initial round, the violation periods of the previous round are the periods when the initial final storage capacity does not meet the storage capacity constraint. The optimized final storage capacity of all violation periods in the previous round is obtained by adjusting the storage capacity constraint during the multi-segment interval optimization of the round. Based on the set of violation periods in the current round, all time periods within the scheduling period are divided to obtain multiple segmented intervals for the current round; the first and last time periods in each segmented interval of the current round correspond to the violation periods in the set of violation periods in the current round, the start and end time periods within the scheduling period; Define any segment interval of the current round as the current segment interval of the current round; Based on the optimized final reservoir capacity of the previous time period within the current segment interval of the current cycle, the optimized final reservoir capacity of the last time period within the current segment interval of the current cycle, and the upstream inflow of each time period within the current segment interval of the current cycle, the water supply and final reservoir capacity of each time period within the current segment interval of the current cycle are determined. When the current segment interval of the current cycle is the first segment interval of the current cycle, the optimized final reservoir capacity of the previous time period within the first time period of the current segment interval of the current cycle is the initial reservoir capacity within the scheduling period. When the current segment interval of the current cycle is the last segment interval of the current cycle, the optimized final reservoir capacity of the last time period within the current segment interval of the current cycle is the final reservoir capacity within the scheduling period. Determine whether the final storage capacity of all time periods within the current round of scheduling meets the storage capacity constraint, and obtain the second judgment result; If the second judgment result is yes, then the water supply and final reservoir capacity of all time periods in the current round of scheduling are determined as the optimal water supply and optimal final reservoir capacity of all time periods in the scheduling period, thus obtaining the optimal water supply decision of the reservoir. If the second judgment result is negative, then the violation period of the current round is determined based on all the time periods that do not meet the reservoir capacity constraint in the current round, and the final reservoir capacity of each violation period in the current round is adjusted to obtain the optimized final reservoir capacity of each violation period in the current round. All violation periods from the initial round to the current round are determined as the set of violation periods for the next round, and the optimization for the next round is carried out until the final reservoir capacity of all time periods in the scheduling meets the reservoir capacity constraint, and the optimal water supply decision of the reservoir is obtained.

2. The reservoir water supply decision optimization method according to claim 1, characterized in that, The basic data of the reservoir includes: the initial reservoir capacity during the scheduling period, the upstream inflow during each period during the scheduling period, the number of periods during the scheduling period, and the reservoir capacity at the end of the scheduling period.

3. The reservoir water supply decision optimization method according to claim 2, characterized in that, Based on the aforementioned basic data of the reservoir, a preliminary water supply decision is made regarding the reservoir, including: The cumulative water supply of the reservoir during the scheduling period is determined based on the initial reservoir capacity during the scheduling period, the upstream inflow during each period during the scheduling period, and the reservoir capacity at the end of the scheduling period. Based on the cumulative water supply of the reservoir during the scheduling period and the number of time periods during the scheduling period, the initial water supply for each time period during the scheduling period is determined. Any time period within the scheduling period is designated as the time period to be calculated. The initial and final reservoir capacity for the period to be calculated is determined based on the initial reservoir capacity during the scheduling period, the upstream inflow from the start of the scheduling period to all periods in the period to be calculated, and the preliminary water supply from the start of the scheduling period to all periods in the period to be calculated.

4. The reservoir water supply decision optimization method according to claim 1, characterized in that, Based on the optimized final reservoir capacity of the previous time period within the first time period of the current segment interval in the current cycle, the optimized final reservoir capacity of the last time period within the current segment interval in the current cycle, and the upstream inflow of each time period within the current segment interval in the current cycle, the water supply and final reservoir capacity of each time period within the current segment interval in the current cycle are determined, including: Based on the optimized final reservoir capacity of the previous period in the first period of the current segment interval in the current round, the optimized final reservoir capacity of the last period in the current segment interval in the current round, and the upstream inflow of each period in the current segment interval in the current round, determine the cumulative water supply of the reservoir in the current segment interval of the current round. Based on the cumulative water supply of the reservoir in the current segment of the current cycle and the number of time periods in the current segment of the current cycle, determine the water supply for each time period in the current segment of the current cycle. Define any time period within the current segment interval of the current round as the target time period within the current segment interval of the current round; Based on the optimized final reservoir capacity of the previous time period within the current segment interval of the current cycle, the upstream inflow of water from the first time period to the target time period within the current segment interval of the current cycle, and the water supply of water from the first time period to the target time period within the current segment interval of the current cycle, the final reservoir capacity of the target time period within the current segment interval of the current cycle is determined, thereby obtaining the final reservoir capacity of each time period within the current segment interval of the current cycle.

5. The reservoir water supply decision optimization method according to claim 4, characterized in that, Adjust the final storage capacity for each violation period in the current round to obtain the optimized final storage capacity for each violation period in the current round, including: When the final storage capacity of each violation period in the current round is less than the minimum storage capacity, the final storage capacity of each violation period in the current round is adjusted to the minimum storage capacity to obtain the optimized final storage capacity of each violation period in the current round. When the final storage capacity of each violation period in the current round is greater than the maximum storage capacity, the final storage capacity of each violation period in the current round is adjusted to the maximum storage capacity to obtain the optimized final storage capacity of each violation period in the current round.

6. A computer device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the reservoir water supply decision optimization method according to any one of claims 1-5.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the reservoir water supply decision optimization method according to any one of claims 1-5.

8. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the reservoir water supply decision optimization method according to any one of claims 1-5.

Images