Successive approximation optimal solution maximum peak clipping criterion reservoir flood control optimal scheduling solving method

Abstract

The invention discloses a maximum peak clipping criterion reservoir flood control optimal scheduling solving method for successive approximation of an optimal solution, and the method comprises the steps: constructing a maximum peak clipping criterion reservoir flood control optimal scheduling model, and inputting the data needed by model solving; defining A = 0 and B = qA; setting the safe discharge q 'MA of the downstream flood control section to be equal to (A + B) / 2; according to the set safety discharge q 'MA of the downstream flood control section, Muskinggen inverse inversion is conducted, and the reservoir safety compensation discharge constraint process of the iteration is obtained through calculation; carrying out moment-by-moment rolling scheduling; checking whether the obtained scheduling process accords with the reservoir water level constraint and the reservoir discharge flow constraint, if so, enabling B to be equal to (A + B) / 2, otherwise, enabling A to be equal to (A + B) / 2; and when A-B is smaller than a threshold value epsilon, outputting a scheduling process. According to the method, the reservoir discharge upper limit constraint is adjusted through iteration, and the reservoir flood control optimization scheduling optimal solution which conforms to the constraint and meets the required precision is gradually approximated.

Description

Technical Field

[0001] The invention belongs to the technical field of reservoir flood control dispatching, and in particular relates to a maximum peak-cutting criterion reservoir flood control optimization dispatching solution method that successively approximates an optimal solution. Background Art

[0002] Reservoir flood control dispatching refers to the planned control and storage of flood water under the premise of ensuring the safety of the project and the flood control requirements of the downstream, so as to achieve the purpose of reducing the harm of floods. Conventional flood control dispatching generally uses flood control dispatching diagrams or dispatching rules for dispatching. Conventional dispatching is highly operational, but it is difficult to fully utilize the flood control storage capacity of the reservoir, and the flood control efficiency is low. In order to further reduce the harm of floods to the downstream, reservoir flood control optimization dispatching is required. Reservoir flood control optimization dispatching is a high-dimensional complex decision-making process with multiple constraints. Its solution algorithms usually include dynamic programming algorithms, step-by-step optimization algorithms, intelligent optimization algorithms, etc.

[0003] When the discrete precision is high, the dynamic programming algorithm takes a long time to solve, and the curse of dimensionality will occur when faced with high-dimensional complex decision-making problems; the solution effect of the step-by-step optimization algorithm depends on the initial solution. When a better initial solution cannot be obtained, it is difficult to gradually optimize to obtain the optimal solution; when processing simple models, the intelligent optimization algorithm has strong search capabilities and fast convergence speed, but its solution space expands exponentially with the number of scheduling times. When faced with high-dimensional complex decision-making problems, it is difficult to effectively search the solution space, and it is easy to fall into the local optimal solution, and the solution result is unstable.

[0004] Therefore, proposing a reservoir flood control optimization scheduling solution method with fast solution speed and the ability to obtain the optimal solution that meets the required accuracy is the key to formulating a reservoir flood control optimization scheduling plan. Summary of the invention

[0005] Technical problem solved: The present invention proposes a method for solving reservoir flood control optimization scheduling with the maximum peak-shaving criterion that approaches the optimal solution one after another. Based on the maximum peak-shaving criterion for reservoir flood control optimization scheduling, the contradictory relationship between the reservoir discharge flow and the upper limit constraint of the reservoir water level is considered, and a feasible interval contraction rule for the reservoir discharge flow is formulated. By iteratively adjusting the upper limit constraint of the reservoir discharge flow, the optimal solution for reservoir flood control optimization scheduling that meets the constraint and meets the required accuracy is approached one after another. The present invention not only has a fast calculation speed, but also can flexibly set the solution accuracy, providing a new method for solving the reservoir flood control optimization scheduling model.

[0006] Technical solution:

[0007] A maximum peak shaving criterion reservoir flood control optimization dispatching solution method that successively approximates the optimal solution, the maximum peak shaving criterion reservoir flood control optimization dispatching solution method comprising the following steps:

[0008] Step 1: Construct a reservoir flood control optimization dispatching model with maximum peak reduction criteria, and input the data required for model solution, including the characteristic water level of the reservoir, the relationship table between the reservoir water level and storage capacity, the relationship table between the reservoir water level and maximum discharge capacity, the starting water level at the beginning of the reservoir dispatching period, the actual safe discharge of the downstream flood control section, the inflow flood flow process, the interval flood flow process, the Muskingum parameter and the time period length;

[0009] Step 2: Define A = 0, B = q A ,q A The actual safe discharge volume of the downstream flood control section;

[0010] Step 3: Set the safe discharge volume q' of the downstream flood control section MA =(A+B) / 2;

[0011] Step 4: According to the set downstream flood control section safe discharge q' MA , Muskingum inverse model is performed, and the reservoir flood control optimization dispatching model based on the maximum peak shaving criterion is used to calculate the reservoir safety compensation discharge flow constraint value for this iteration;

[0012] Step 5: Perform rolling scheduling moment by moment; specifically:

[0013] Take the current moment reservoir discharge flow equal to the current moment reservoir safety compensation discharge flow constraint value, check whether it meets the reservoir maximum discharge capacity constraint, if not, take the current moment reservoir safety compensation discharge flow constraint value and the reservoir maximum discharge capacity constraint value as the current moment reservoir discharge flow, rolling moment by moment, to get the dispatching process;

[0014] Step 6: Based on the maximum peak shaving criterion reservoir flood control optimization dispatching model, check whether the dispatching process meets the constraints, including reservoir water level constraints and reservoir discharge flow constraints. If so, set B = (A + B) / 2, otherwise set A = (A + B) / 2;

[0015] Step 7: Determine whether |AB| is less than the threshold ε. If so, output the scheduling process, otherwise return to step 3; ε is related to the discharge accuracy.

[0016] Furthermore, in step 1, the objective function of the maximum peak shaving criterion reservoir flood control optimization scheduling model is:

[0017]

[0018] Where: f is the target of the reservoir flood control optimization scheduling model with maximum peak shaving criterion; q t is the discharge flow of the reservoir at time t; Q 区,t is the interval flow at time t; T is the scheduling time number;

[0019] The constraints of the maximum peak shaving criterion reservoir flood control optimization scheduling model include:

[0020] (1) Water balance constraints:

[0021] V t+1 =V t +(Q t+1 +Q t -q t+1 -q t )Δt / 2

[0022] Where: V t 、V t+1 is the reservoir capacity at time t and time t+1; Q t , Q t+1 is the inflow of the reservoir at time t and time t+1; q t ,q t+1 is the discharge flow of the reservoir at time t and time t+1; Δt is the time period length;

[0023] (2) Reservoir water level upper and lower limit constraints:

[0024] Z 死 ≤Z t ≤Z 防

[0025] Where: Z 死 is the dead water level; Z t is the water level of the reservoir at time t; Z 防 To prevent high water levels from flooding;

[0026] (3) Reservoir initial water level constraint:

[0027] Z0=Z start

[0028] Where: Z0 is the initial water level of the reservoir; Z start The starting water level at the beginning of the reservoir dispatch period;

[0029] (4) Constraints on the maximum discharge capacity of the reservoir:

[0030] q t ≤q 泄,t

[0031] Where: q 泄,t is the maximum discharge capacity of the reservoir at time t;

[0032] (5) Reservoir safety compensation discharge flow constraints:

[0033] q t ≤q 补,t

[0034]

[0035] Where: q 补,t ,q 补,t+1 The maximum compensatory discharge flow of the reservoir to ensure that the downstream flood control section does not exceed the safe discharge volume at time t and time t+1; q' MA The safe discharge of the downstream flood control section set for dispatching shall not be greater than the actual safe discharge of the downstream flood control section q A ;Q 区,t , Q 区,t+1 is the interval flow at time t and time t+1; C0, C1, C2 are Muskingum method parameters; q 补,t The boundary value is taken as q 补,T =q' MA -Q 区,T .

[0036] Step 5 further includes:

[0037] Step 51, calculate the reservoir safety compensation flow constraint process q 补,t , t=1,2,3,…,T;

[0038] Step 52, giving an initial storage capacity V1;

[0039] Step 53, set t=1;

[0040] Step 54, let q t =q 补,t , calculate the reservoir water level at time t, and the corresponding maximum discharge capacity q at time t 泄,t ;

[0041] Step 55, determine q t Is it less than or equal to q? 泄,t , if yes, go to step 57, otherwise, go to step 56;

[0042] Step 56: Calculate the maximum discharge capacity q of the reservoir at time t using flood control calculations 泄,t , let q t =q 泄,t ; Go to step 57;

[0043] Step 57, let t=t+1, determine whether t is greater than T, if so, output the scheduling process and end the process; otherwise, go to step 54.

[0044] Furthermore, in step 6, the following conditions are judged to check whether the dispatching process complies with the reservoir water level constraint and the reservoir discharge flow constraint:

[0045] a. Determine whether the current discharge flow of the reservoir is less than or equal to the reservoir safety compensation discharge flow at that moment;

[0046] b. Determine whether the current discharge flow of the reservoir is less than or equal to the maximum discharge capacity of the reservoir at that moment;

[0047] c. Determine whether the current discharge flow of the reservoir is greater than or equal to the minimum discharge flow of the reservoir;

[0048] d. Determine whether the current water level of the reservoir is less than or equal to the flood control high water level;

[0049] e. Determine whether the current water level of the reservoir is greater than or equal to the dead water level;

[0050] If all conditions are met, the dispatching process is considered to comply with the reservoir water level constraints and reservoir discharge flow constraints.

[0051] Beneficial effects:

[0052] First, the maximum peak-cutting criterion reservoir flood control optimization scheduling solution method of the present invention, which successively approaches the optimal solution, continuously updates the values ​​of the left and right boundaries A and B of the safe discharge through the dichotomy method, thereby approaching the minimum q' that makes the scheduling process meet the constraints. MA The corresponding dispatching process satisfies the goal of minimizing the maximum flow rate of the downstream flood control section during the dispatching period and achieves the optimal solution. This method has a specified optimization direction and can stably obtain the theoretical optimal solution.

[0053] Secondly, the maximum peak-shaving criterion reservoir flood control optimization scheduling solution method of the present invention, which successively approximates the optimal solution, does not require a large number of discrete and repetitive solution space searches, and has the advantages of high efficiency and high speed.

[0054] Third, the maximum peak-shaving criterion reservoir flood control optimization scheduling solution method of the present invention, which successively approximates the optimal solution, has the advantages of stable solution and fast solution speed. It can not only be used for single reservoir optimization scheduling solution, but also is suitable for combination with system decomposition and coordination methods as a solution algorithm for the single reservoir subsystem model in a complex flood control system. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] Figure 1 This is a schematic block diagram of the implementation steps of the maximum peak-shaving criterion reservoir flood control optimization scheduling solution method that successively approximates the optimal solution proposed in the present invention.

[0056] Figure 2 The result diagram of applying the invented solution algorithm to the Zhexi single reservoir flood control optimization scheduling model. DETAILED DESCRIPTION

[0057] The following examples will enable those skilled in the art to more fully understand the present invention, but are not intended to limit the present invention in any way.

[0058] See also Figure 1The present invention discloses a method for solving the optimal dispatching of reservoir flood control based on the maximum peak shaving criterion that approaches the optimal solution one by one. The method comprises the following steps:

[0059] Step 1: Construct a reservoir flood control optimization dispatching model with maximum peak reduction criteria, and input the data required for model solution, including the characteristic water level of the reservoir, the relationship table between the reservoir water level and storage capacity, the relationship table between the reservoir water level and maximum discharge capacity, the starting water level at the beginning of the reservoir dispatching period, the actual safe discharge of the downstream flood control section, the inflow flood flow process, the interval flood flow process, the Muskingum parameter and the time period length;

[0060] Step 2: Define A = 0, B = q A ,q A The actual safe discharge volume of the downstream flood control section;

[0061] Step 3: Set the safe discharge volume q' of the downstream flood control section MA =(A+B) / 2;

[0062] Step 4: According to the set downstream flood control section safe discharge q' MA , Muskingum inverse model is performed, and the reservoir flood control optimization dispatching model based on the maximum peak shaving criterion is used to calculate the reservoir safety compensation discharge flow constraint value for this iteration;

[0063] Step 5: Perform rolling scheduling at each moment; specifically:

[0064] Take the current moment reservoir discharge flow equal to the current moment reservoir safety compensation discharge flow constraint value, check whether it meets the reservoir maximum discharge capacity constraint, if not, take the current moment reservoir safety compensation discharge flow constraint value and the reservoir maximum discharge capacity constraint value as the current moment reservoir discharge flow, rolling moment by moment, to get the dispatching process;

[0065] Step 6: Based on the reservoir flood control optimization scheduling model with maximum peak shaving criterion, check whether the obtained scheduling process meets the reservoir constraints, including water level constraints and reservoir discharge flow constraints. If so, set B = (A + B) / 2, otherwise set A = (A + B) / 2;

[0066] Step 7: Determine whether |AB| is less than the threshold ε. If so, output the scheduling process, otherwise return to step 3; ε is related to the discharge accuracy.

[0067] Under the condition of known inflow flood and downstream water flow, the maximum peak reduction target of reservoir flood control scheduling is achieved, that is, the maximum flow rate of downstream flood control section is minimized while satisfying the constraints. The objective function and constraints of the optimization scheduling model are as follows:

[0068] The objective function expression is as follows:

[0069]

[0070] Where: f is the target of the reservoir flood control optimization scheduling model with maximum peak shaving criterion; q t is the discharge flow of the reservoir at time t; Q 区,t is the interval flow at time t; T is the scheduling time number.

[0071] The constraints are as follows:

[0072] (1) Water balance constraints

[0073] V t+1 =V t +(Q t+1 +Q t -q t+1 -q t )Δt / 2

[0074] Where: V t 、V t+1 is the reservoir capacity at time t and time t+1; Q t , Q t+1 is the inflow of the reservoir at time t and time t+1; q t ,q t+1 is the discharge flow of the reservoir at time t and time t+1; Δt is the time period length.

[0075] (2) Reservoir water level upper and lower limit constraints

[0076] Z 死 ≤Z t ≤Z 防

[0077] Where: Z 死 is the dead water level; Z t is the water level of the reservoir at time t; Z 防 To prevent high water levels from flooding.

[0078] (3) Reservoir initial water level constraints

[0079] Z0=Z start

[0080] Where: Z0 is the initial water level of the reservoir; Z start The starting water level for the reservoir at the beginning of the regulation period is usually the flood limit water level.

[0081] (4) Constraints on the maximum discharge capacity of the reservoir

[0082] q t ≤q 泄,t

[0083] Where: q 泄,tis the maximum discharge capacity of the reservoir at time t.

[0084] (5) Reservoir safety compensation discharge flow constraints

[0085] In order to ensure the safety of the downstream flood control section, the interval flow is taken into account and the reservoir safety compensation discharge flow is calculated according to the Muskingum inverse method. The river section outflow process in the Muskingum formula is taken as the safety discharge of the downstream flood control section minus the interval flood flow process. The Muskingum method is used to perform the inverse modeling of the river flow. The calculated inflow process is the reservoir safety compensation discharge flow process. Replace the corresponding variables and transform the Muskingum formula. The expression is as follows:

[0086] q t ≤q 补,t

[0087]

[0088] Where: q 补,t ,q 补,t+1 The maximum compensatory discharge flow of the reservoir to ensure that the downstream flood control section does not exceed the safe discharge volume at time t and time t+1; q' MA The safe discharge of the downstream flood control section set for dispatching must not be greater than the actual safe discharge of the downstream flood control section q A ;Q 区,t , Q 区,t+1 is the interval flow at time t and time t+1; C0, C1, C2 are Muskingum method parameters; q 补,t The boundary value is taken as q 补,T =q' MA -Q 区,T .

[0089] The decision variable of the optimal dispatching model is the reservoir discharge flow process, which is a T×1 matrix. Analyzing the model, there is the following relationship between the reservoir discharge flow and the maximum water level of the reservoir dispatching:

[0090] (1) Safe discharge volume q' of the downstream flood control section MA The smaller it is, the smaller the maximum flow rate of the downstream flood control section will be, the smaller the upper limit of the reservoir discharge flow will be, the higher the maximum water level of the reservoir during the scheduling process will be, and the easier it will be to break through the upper limit constraint of the reservoir water level.

[0091] (2) Safe discharge volume q' of the downstream flood control section MA The larger it is, the greater the maximum flow rate of the downstream flood control section is, and the greater the upper limit of the reservoir discharge flow is, and the easier it is to meet the upper limit constraint of the reservoir water level.

[0092] Therefore, to minimize the maximum flow rate of the downstream flood control section, it is only necessary to obtain the minimum q' that makes the model have a feasible solution.MA value.

[0093] Therefore, according to the invention, a method for solving the optimal dispatching of reservoir flood control with the maximum peak-cutting criterion that approaches the optimal solution one by one is characterized by comprising the following specific steps:

[0094] Step (1): Construct a reservoir flood control optimization scheduling model with maximum peak shaving criterion, and input the data required for model solution, including the characteristic water level of the reservoir, the relationship table between the reservoir water level and storage capacity, the relationship table between the reservoir water level and maximum discharge capacity, the starting water level at the beginning of the reservoir scheduling period, the actual safe discharge of the downstream flood control section, the inflow flood flow process, the interval flood flow process, the Muskingum parameter and the time period length.

[0095] Step (2): Define A = 0, B = q A .q A It is the actual safe discharge volume of the downstream flood control section.

[0096] Step (3): Set the safe discharge volume q' of the downstream flood control section MA =(A+B) / 2.

[0097] Step (4): According to the set downstream flood control section safe discharge q' MA , Muskingum inverse modeling is performed, and the reservoir flood control optimization scheduling model based on the maximum peak shaving criterion is used to calculate the reservoir safety compensation discharge flow constraint process for this iteration.

[0098] Step (5): Perform rolling scheduling moment by moment; specifically:

[0099] Take the current moment reservoir discharge flow equal to the current moment reservoir safety compensation discharge flow constraint value, check whether it meets the reservoir maximum discharge capacity constraint, if not, take the current moment reservoir safety compensation discharge flow constraint value and the reservoir maximum discharge capacity constraint value as the current moment reservoir discharge flow, rolling moment by moment, to get the dispatching process;

[0100] Step (6): Based on the maximum peak shaving criterion reservoir flood control optimization dispatching model, check whether the dispatching process meets the constraints, including reservoir water level constraints and reservoir discharge flow constraints. If it meets the constraints, set B = (A + B) / 2, otherwise set A = (A + B) / 2.

[0101] Step (7): Determine whether |AB| is less than a threshold value ε. If so, output the scheduling process, otherwise return to step (3); ε is related to the discharge accuracy.

[0102] A further preferred embodiment of the present invention is:

[0103] In step (6) of the present invention, checking whether the constraints are met requires rolling moment by moment to determine whether the following constraints are met:

[0104] a. The reservoir discharge flow at that moment must be less than or equal to the reservoir safety compensation discharge flow at that moment.

[0105] b. The outflow rate of the reservoir at that moment must be less than or equal to the maximum discharge capacity of the reservoir at that moment.

[0106] c. The outflow rate of the reservoir at that moment must be greater than or equal to the minimum outflow rate of the reservoir.

[0107] d. The water level of the reservoir at that moment must be less than or equal to the flood control high water level.

[0108] e. The water level of the reservoir at that moment must be greater than or equal to the dead water level.

[0109] The specific implementation of the present invention is now further described by taking the application of the present invention to the Zhexi Reservoir Project in the Zishui River Basin as an example. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0110] 1. Basic Information

[0111] a. Characteristic water level of the reservoir:

[0112] Dead water level: Z 死 =144m; flood limit water level: Z 汛 =162m; Flood control high water level: Z 防 =170m.

[0113] b. Water level and storage capacity relationship table of Zhexi Reservoir, see Table 1 for details.

[0114] c. The relationship table between the water level and maximum discharge capacity of Zhexi Reservoir is shown in Table 2.

[0115] d. The starting water level Z at the beginning of the reservoir operation period start =162m

[0116] e. The downstream flood control section is Taojiang section, and the actual safe discharge volume is q A =12000m 3 / s.

[0117] f. Typical flood data (p=2%): Data table of flood flow process into Zhexi Reservoir and flood flow process in Zhetao area, see Table 3 for details.

[0118] g. Muskingum parameters, C0 = 0.187192, C1 = 0.674877, C2 = 0.137831.

[0119] h. Time period length Δt = 2h. Table 1: Relationship between water level and storage capacity of Zhexi Reservoir

[0120] Water level(m) <![CDATA[Storage capacity (100 million m 3 )]]> Water level(m) <![CDATA[Storage capacity (100 million m 3 )]]> Water level(m) <![CDATA[Storage capacity (100 million m 3 )]]> 144 7.97 154.5 16.19 165 24 144.5 8.34 155 16.61 165.5 25.45 145 8.7 155.5 17.03 166 26.33 145.5 9.07 156 17.32 166.5 26.65 146 9.44 156.5 17.58 167 26.98 146.5 9.8 157 17.85 167.5 27.3 147 10.17 157.5 18.12 168 28.02 147.5 10.54 158 18.39 168.5 28.75 148 10.9 158.5 18.65 169 29.48 148.5 11.27 159 18.92 169.5 30.2 149 11.64 159.5 19.19 170 31 149.5 12.01 160 19.45 170.5 31.55 150 12.41 160.5 19.72 171 32.1 150.5 12.83 161 19.99 171.5 32.9 151 13.25 161.5 20.25 172 33.78 151.5 13.67 162 20.52 172.5 34.66 152 14.09 162.5 21.1 173 36.4 152.5 14.51 163 21.68 173.5 37.45 153 14.93 163.5 22.26 174 38.49 153.5 15.35 164 22.84 174.21 38.8 154 15.77 164.5 23.42

[0121] Table 2: Relationship between water level and maximum discharge capacity of Zhexi Reservoir

[0122]

[0123]

[0124] Table 3: Data table of flood flow process into Zhexi Reservoir and flood flow process in Zhetao section

[0125]

[0126] 2. Implementation Process

[0127] Step 1: Construct a reservoir flood control optimization scheduling model based on the maximum peak shaving criterion, and input the data required for model solution, including the reservoir characteristic water level, the reservoir water level and storage capacity relationship table, the reservoir water level maximum discharge capacity relationship table, the starting water level at the beginning of the reservoir characteristic water level scheduling period, the actual safe discharge of the downstream flood control section, the flood inflow process, the interval flow process, the Muskingum parameter and the time period length.

[0128] Step 2, define A = 0, B = q A =12000m 3 / s.

[0129] Step 3: Set the safe discharge volume q' of the downstream flood control section MA =(A+B) / 2=6000m 3 / s.

[0130] Step 4: According to the set downstream flood control section safe discharge q' MA , Muskingum inverse modeling is performed, and the reservoir flood control optimization scheduling model based on the maximum peak shaving criterion is used to calculate the reservoir safety compensation discharge flow constraint process of this iteration, as shown in Table 4.

[0131] Table 4 Muskingum inverse solution of reservoir safety compensation discharge flow process table

[0132]

[0133]

[0134]

[0135] Step 5, perform rolling dispatch at each moment; specifically: take the reservoir discharge flow at the current moment equal to the reservoir safety compensation discharge flow constraint value at the current moment, check whether it meets the reservoir maximum discharge capacity constraint, if not, take the minimum value of the reservoir safety compensation discharge flow constraint value and the reservoir maximum discharge capacity constraint value at the current moment as the reservoir discharge flow at the current moment, and roll at each moment to obtain the dispatch process. As shown in Table 5.

[0136] Table 5 Results of the first iteration of reservoir flood control optimization scheduling

[0137]

[0138]

[0139]

[0140] Step 6: Based on the maximum peak-cutting criterion reservoir flood control optimization dispatching model, check whether the dispatching process meets the constraints, including reservoir water level constraints and reservoir discharge flow constraints. If it meets the constraints, set B = (A + B) / 2, otherwise set A = (A + B) / 2.

[0141] It is easy to observe that the reservoir water level in Table 2 (the first iteration result table) exceeds the flood control high water level by 170m, which does not meet the constraint. Therefore, let A = (A + B) / 2 = (0 + 12000) / 2 = 6000m 3 / s.

[0142] Step 7, determine whether |AB| is less than the threshold ε. If so, output the scheduling process, otherwise return to step 3; ε is related to the discharge accuracy and is taken as 0.1.

[0143] Continue to iterate until the result is output. The iterative process is shown in Table 6.

[0144] Table 6 Constraint parameters of reservoir safety compensation discharge flow Iteration process table

[0145]

[0146]

[0147] After 17 iterations, q' MA =9388.27m 3 The scheduling results corresponding to / s are shown in Table 7.

[0148] Table 7 Final scheduling results

[0149]

[0150]

[0151]

[0152] The invented solution algorithm was applied to solve the Zhexi single reservoir flood control optimization dispatching model. When encountering a flood that occurs once in 50 years, the peak flow of the downstream flood control section can be reduced to 9388m 3 / s.

[0153] The above are only preferred embodiments of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions under the concept of the present invention belong to the protection scope of the present invention. It should be pointed out that for ordinary technicians in this technical field, some improvements and modifications without departing from the principle of the present invention should be regarded as the protection scope of the present invention.

Claims

1. A method for solving reservoir flood control optimization scheduling with maximum peak shaving criterion that approaches the optimal solution successively, characterized in that: The maximum peak shaving criterion reservoir flood control optimization dispatch solution method comprises the following steps: Step 1: Construct a reservoir flood control optimization dispatching model with maximum peak reduction criteria, and input the data required for model solution, including the characteristic water level of the reservoir, the relationship table between the reservoir water level and storage capacity, the relationship table between the reservoir water level and maximum discharge capacity, the starting water level at the beginning of the reservoir dispatching period, the actual safe discharge of the downstream flood control section, the inflow flood flow process, the interval flood flow process, the Muskingum parameter and the time period length; Step 2: Define A = 0, B = q A ,q A The actual safe discharge volume of the downstream flood control section; Step 3: Set the safe discharge volume q' of the downstream flood control section MA =(A+B) / 2; Step 4: According to the set downstream flood control section safe discharge q' MA , Muskingum inverse model is performed, and the reservoir flood control optimization dispatching model based on the maximum peak shaving criterion is used to calculate the reservoir safety compensation discharge flow constraint value for this iteration; Step 5: Perform rolling scheduling moment by moment; specifically: Take the current moment reservoir discharge flow equal to the current moment reservoir safety compensation discharge flow constraint value, check whether it meets the reservoir maximum discharge capacity constraint, if not, take the current moment reservoir safety compensation discharge flow constraint value and the reservoir maximum discharge capacity constraint value as the current moment reservoir discharge flow, rolling moment by moment, to get the dispatching process; Step 6: Based on the maximum peak shaving criterion reservoir flood control optimization dispatching model, check whether the dispatching process meets the constraints, including reservoir water level constraints and reservoir discharge flow constraints. If so, set B = (A + B) / 2, otherwise set A = (A + B) / 2; Step 7: Determine whether |AB| is less than the threshold ε. If so, output the scheduling process, otherwise return to step 3; ε is related to the discharge accuracy.

2. The maximum peak-cutting criterion reservoir flood control optimization dispatching solution method according to claim 1 is characterized in that: In step 1, the objective function of the maximum peak shaving criterion reservoir flood control optimization scheduling model is: Where: f is the target of the reservoir flood control optimization scheduling model with maximum peak shaving criterion; q t is the discharge flow of the reservoir at time t; Q 区,t is the interval flow at time t; T is the scheduling time number; The constraints of the maximum peak shaving criterion reservoir flood control optimization scheduling model include: (1) Water balance constraints: 5 t+1 =V t +(Q t+1 +Q t -q t+1 -q t )Δt / 2 Where: V t 、V t+1 is the reservoir capacity at time t and time t+1; Q t , Q t+1 is the inflow of the reservoir at time t and time t+1; q t ,q t+1 is the discharge flow of the reservoir at time t and time t+1; Δt is the time period length; (2) Reservoir water level upper and lower limit constraints: WITH 死 ≤Z t ≤Z 防 Where: Z 死 is the dead water level; Z t is the water level of the reservoir at time t; Z 防 To prevent high water levels from flooding; (3) Reservoir initial water level constraint: Z0=Z start Where: Z0 is the initial water level of the reservoir; Z start The starting water level at the beginning of the reservoir dispatch period; (4) Constraints on the maximum discharge capacity of the reservoir: q t ≤q 泄,t Where: q 泄,t is the maximum discharge capacity of the reservoir at time t; (5) Reservoir safety compensation discharge flow constraints: q t ≤q 补,t Where: q 补,t ,q 补,t+1 The maximum compensatory discharge flow of the reservoir to ensure that the downstream flood control section does not exceed the safe discharge volume at time t and time t+1; q' MA The safe discharge of the downstream flood control section set for dispatching shall not be greater than the actual safe discharge of the downstream flood control section q A ;Q 区,t , Q 区,t+1 is the interval flow at time t and time t+1; C0, C1, C2 are Muskingum method parameters; q 补,t The boundary value is taken as q 补,T =q' MA -Q 区,T .

3. The maximum peak-cutting criterion reservoir flood control optimization dispatching solution method according to claim 1 is characterized in that: Step 5 further includes: Step 51, calculate the reservoir safety compensation flow constraint process q 补,t , t=1,2,3,…,T; Step 52, giving an initial storage capacity V1; Step 53, set t=1; Step 54, let q t =q 补,t , calculate the reservoir water level at time t, and the corresponding maximum discharge capacity q at time t 泄,t ; Step 55, determine q t Is it less than or equal to q? 泄,t , if yes, go to step 57, otherwise, go to step 56; Step 56: Calculate the maximum discharge capacity q of the reservoir at time t using flood control calculations 泄,t , let q t =q 泄,t ; Go to step 57; Step 57, let t=t+1, determine whether t is greater than T, if so, output the scheduling process and end the process; otherwise, go to step 54.

4. The maximum peak-cutting criterion reservoir flood control optimization dispatching solution method according to claim 1 is characterized in that: In step 6, the following conditions are judged to check whether the dispatching process meets the reservoir water level constraint and the reservoir discharge flow constraint: a. Determine whether the current discharge flow of the reservoir is less than or equal to the reservoir safety compensation discharge flow at that moment; b. Determine whether the current discharge flow of the reservoir is less than or equal to the maximum discharge capacity of the reservoir at that moment; c. Determine whether the current discharge flow of the reservoir is greater than or equal to the minimum discharge flow of the reservoir; d. Determine whether the current water level of the reservoir is less than or equal to the flood control high water level; e. Determine whether the current water level of the reservoir is greater than or equal to the dead water level; If all conditions are met, the dispatching process is considered to comply with the reservoir water level constraints and reservoir discharge flow constraints.

Images