CN103646294B - A kind of power station Flood Dispatching Optimization method considering operating gate rule - Google Patents

Abstract

The invention discloses a kind of power station Flood Dispatching Optimization method considering operating gate rule, including five steps: (1) operating gate rule structuralisation describes, solve the statement problem of complicated operating gate rule, (2) operating gate rule instantiation, solves how to apply operating gate regular texture body example to the problem stating specific gate of hydropower station operation rules；(3) gatage combination sequence optimisation generates, and solves the coupled problem of operating gate rule and tradition power station Flood Dispatching Optimization method；(4) gatage combination sequence optimisation adjusts, and during solving Dispatching Flood, indivedual period water levels are out-of-limit and dispatch the problem that end of term water level exists deviation with target water level；(5) operating gate plan generates, the problem solving dynamically to generate operating gate plan according to gatage combination sequence.The present invention improves the accuracy that later stage flood discharge discharge process is predicted by power station, improves the capability of fast response of power station flood decision.

Description

A kind of power station Flood Dispatching Optimization method considering operating gate rule

Technical field

The present invention relates to a kind of Flood Dispatching Optimization method, be specifically related to a kind of water power considering operating gate rule Stand Flood Dispatching Optimization method.

Background technology

All gates are considered as a scheduler object by traditional power station Flood Dispatching Optimization method, only consider all The overall discharge capacity of gate, does not consider the gates such as the permission aperture of each gate, opening sequence, aperture combination Operation rules, is only capable of calculating optimization total flood discharge discharge curve that gate of hydropower station group is overall, and cannot be automatic Make the optimization operation planning of each gate.In actual motion link, need to be by dispatcher according to this optimization Total flood discharge discharge curve, the artificial operation planning formulating each gate.This mode is difficult to meet power station essence Refinement scheduling and the business demand of real-time rolling scheduling, also had a strong impact on the efficiency of flood decision simultaneously.Additionally, Owing to actual motion must comply with the operation rules of each gate, it is impossible to regulating gate it is anticipated that optimization total Flood discharge discharge curve runs, and causes gate actual total between letdown flow and the total letdown flow of intended optimization Deviation, and then cause the deviation of later stage actual library water level and expection reservoir level, it is further exacerbated by actual flood discharge flow And the extent of deviation between expection flood discharge flow, have a strong impact on the reality of traditional power station Flood Dispatching Optimization method By value.

Power station Flood Dispatching Optimization mainly solves the optimization storage of flood into reservoir and lets out problem, is guaranteeing power station self On the premise of safety, alleviate the flood damaging effect to downstream as far as possible, and the resource profit of with due regard to flood With.All gates are considered as an object by tradition power station Flood Dispatching Optimization method, are calculated by Optimized Operation Obtain the optimization flood discharge discharge curve in whole power station, in the process it is not intended that the operation of each gate is advised Then.In practical implementation, need the experience by dispatcher and continuous tentative calculation, with reference to the expection obtained Optimize flood discharge discharge curve and actual reservoir level, determine that meeting each gate allows aperture, opening sequence, opens The operation planning of the operation rules such as degree combination, and the artificial gatage regulating release floodwatering facility in real time so that actual Flood discharge flow with expection flood discharge flow as close possible to.Owing to have ignored the operation rules of each gate, tradition water power The optimization flood discharge discharge curve overflow performance true with power station that Flood Dispatching Optimization method of standing obtains is not inconsistent, from And cause perform link and plan link substantial deviation, significantly reduce conventional method result of calculation value and Referring to property.

The final purpose of Flood Dispatching Optimization is by the rational management to flood so that national economy loss reduction Change.But, the national economy loss that flood is caused itself is difficult to quantification and calculates, and therefore power station flood is excellent Changing scheduling generally uses the mode simplified to consider flood losses.Generally, its people's warp that flood is caused Ji loss and maximum flood discharge being proportionate property of flow, can be maximum by national economy minimization of loss targeted transformation Flood discharge flow minimizes target.Utilize mathematical derivation may certify that, in the Dispatching Flood phase, use the flood control storehouse determined Rong Shi, the scheme that maximum flood discharge flow minimizes is actually the most uniform scheme of flood discharge flow.Due to flood Initial water level and the end water level of scheduling determine the most in advance, and the water yield of flood warehouse-in determines the most in advance, therefore basis The principle of water balance can determine the total flood discharge water yield in schedule periods in advance, it means that day part in schedule periods Flood discharge flow sum be a fixed value.According to mathematical principle, when the timing of multiple independent variable sums one, permissible Its uniformity coefficient is stated by the quadratic sum of each independent variable.Therefore, it can set up following object function:

$\mathrm{Min}{\∫}_{t_0}^{t_d}{q}^2\left(t\right)\mathrm{dt},t\∈[t_0,t_d]$

For the ease of using computer solving, above-mentioned object function can be carried out discretization:

$\mathrm{Min}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{q}_i^2\·\mathrm{\Δt},i=\mathrm{1,2},L,N;N=\frac{t_d-t_0}{\mathrm{\Δt}};$

Constraints:

Storage capacity retrains $\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(Q_i-q_i)\·\mathrm{\Δt}={\mathrm{\ΔV}}_f,i=\mathrm{1,2},L,N$

Traffic constraints q during without letting out in advance_i≤Q_max, i=1,2, L, N

Maximum discharge capacity constraint q_i≤q_max(z_i), i=1,2, L, N

Reservoir water yield Constraints of Equilibrium V_i+1=V_i+(Q_i-q_i) Δ t, i=1,2, L, N

In formula, t₀For schedule periods initial time, t_dFor schedule periods finish time, q (t) is flood discharge discharge process, q_iFor the flood discharge flow of i-th period, Q_maxFor the maximum reservoir inflow (crest discharge) of flood, when Δ t is Segment length, hop count when N is total, Δ V_fThe storage capacity employed for plan, z_iFor the reservoir level of i-th period, q_max(z_i) it is that all release floodwatering facilities are at z_iMaximum discharge capacity under reservoir level, can be bent according to maximum discharge capacity Line interpolation obtains, V_iFor the storage capacity at the beginning of the i-th period.

Intend storage capacity size delta V employed_fShould be less than the total storage capacity equal to reservoir.ΔV_fSpecific size Should be according to later stage inflow forecast situation and combine artificial experience and determine.When later stage water is more, should set less Δ V_f, make schedule periods interior with flood discharge water under bigger peak clipping flow, in order to reservoir is most during flood water-break It is reduced to below flood season limit level soon.When later stage water is less, bigger Δ V can be set_f, make schedule periods interior with relatively Flood discharge water under little peak clipping flow, reduces downstream flood loss.

According to above-mentioned Mathematical Modeling, use traditional dynamic programming algorithm can be calculated the excellent of whole power station Change flood discharge discharge curve.

From above-mentioned Mathematical Modeling it can be seen that prior art only accounts for all release floodwatering facilities in power station Big discharge capacity constraint, does not consider the gate behaviour such as the permission aperture of each gate, opening sequence, aperture combination Make the impact that flood discharge process is produced by rule, be only capable of obtaining theoretic optimization flood discharge discharge curve, and cannot Obtain each operating gate plan corresponding to this optimization flood discharge discharge curve.During actual motion, need to rely on people Work monitors that reservoir level changes in real time, decides whether out according to current reservoir level and intended optimization flood discharge flow Opening or closed shutter, and calculate corresponding gatage, workload is the biggest.Additionally, due to release floodwatering facility lock The aperture of door is all confined to several discrete fixing opening values, and actual gate flood discharge flow is let out with expection optimization Bigger deviation can be there is in mighty torrent amount.This deviation all exists in whole schedule periods, and the constantly accumulation of deviation makes The optimization flood discharge graph of schedule periods posterior part loses prediction and directive significance.Traditional rolling scheduling method is also Cannot thoroughly solve this problem so that dispatcher cannot reservoir level and letdown flow in accurately anticipation goes out schedule periods Development trend, result even in the dispatcher lacked experience and judge flood-proofing terrain by accident, bring relatively to flood control Big potential safety hazard.

Summary of the invention

The deficiency existed for prior art, it is an object of the present invention to provide a kind of water power considering operating gate rule Stand Flood Dispatching Optimization method, solve conventional method actual execution link and there is severe deviations with plan link Problem, improves power station and optimizes the practice guiding action of flood discharge discharge curve, reduce the work of dispatcher Making intensity and error probability, meanwhile, auxiliary dispatching personnel study and judge the change of later stage flood-proofing terrain more accurately.

To achieve these goals, the present invention is to realize by the following technical solutions:

The present invention considers the power station Flood Dispatching Optimization method of operating gate rule, includes following step successively Rapid:

(1) operating gate rule structuralisation describes, and uses structurized description method to carry out complete expression gate and allows The complex informations such as aperture, opening sequence and aperture combination, for operating gate rule instantiation and computer analyzing Operating gate rule provides basis；

(2) operating gate rule instantiation is carried out for specific power station, according to the true gate in actual power station Allow the information such as aperture, opening sequence and aperture combination, dynamically generate according to certain rule and arrange all kinds of lock Multiple examples of door operation regular texture body, and set up the incidence relation between different instances；

(3) generate the gatage combination sequence optimized, push away according to tradition power station Flood Dispatching Optimization method The theoretical flood discharge discharge curve asked, calculates optimum aperture combination sequence number by the period, is formed in whole schedule periods Optimize gatage combination sequence, for solving the coupling of operating gate rule and tradition Flood Dispatching Optimization method Conjunction problem；

(4) gatage combination sequence is optimized adjustment, for reservoir levels indivedual during Dispatching Flood The out-of-limit period, adjust the previous operating gate moment by designation method, until eliminating the out-of-limit problem of reservoir level.Pin To dispatching the deviation situation of end of term reservoir level and target water level and previous operating gate situation, adjust by designation method The last operating gate moment, until scheduling end of term water level is closest；

(5) operating gate plan is generated, by optimizing resolving by the period of gatage combination sequence, dynamic State generates and meets the operating gate plan that all operating gate rules require.

In step (1), described operating gate rule includes the permission aperture of gate, opening sequence and aperture combination Information, its description method is as follows:

Be respectively adopted aperture structure (hereinafter referred to as OM structure), (hereinafter referred to as GT ties gate structure body Structure body) and aperture combining structure body (hereinafter referred to as CB structure) characterize the permission aperture of gate, unlatching Order and aperture combined information.

Described OM structure is for storing the information of the specific aperture of gate, and its content includes aperture sequence number, aperture The stage discharge relation that value, aperture type, aperture describe and be made up of several water levels, flow binary information group Curve；

Described GT structure is for storing the information of the specific gate in power station, and its content includes gate sequence number, lock The aperture sequence number that door describes and allows, GT structure is associated with OM structure by aperture sequence number；

Described CB structure is used for storing specific gatage combined information, its content include number of regulation and The aperture combined information being made up of several gate sequence numbers, aperture sequence number binary information group, CB structure passes through Gate sequence number associates with GT structure with aperture sequence number.

In step (2), operating gate rule is instantiated, refer to allow for the gate in actual power station Aperture, opening sequence and operation rules information, respectively instantiation OM structure, GT structure and CB structure Body, and form the aperture table for this power station, gate table and combination table.Specifically include following steps:

First, allowing aperture according to every kind of gate in actual power station, the OM structure building correspondence respectively is real Example, states this gate and allows the relevant information of aperture.All of OM structure example collectively forms complete opening Degree table.

Secondly, according to each floodgate in actual power station, build the GT structure example of correspondence respectively, State the relevant information of this gate.All of GT structure example collectively forms complete gate table.

Finally, according to every kind of gatage combination in actual power station, the CB structure building correspondence respectively is real Example, states the relevant information of this aperture combination.All of CB structure collectively forms complete combination table.

Rating curve in described each OM structure example should be arranged in strict accordance with water level progressive law Row, in order to calculate, according to reservoir level Fast Interpolation, the flood discharge flow that this aperture is corresponding in solution procedure.

Each OM structure example in described aperture table can arrange from small to large according to opening value；Described gate Each GT structure example in table can arrange according to the type of gate and physical location；In described combination table Each CB structure example then should carry out ordered arrangement in strict accordance with gate opening sequencing and combined situation.

Between described combination table each CB structure example, order is the most correct to take following methods to check:

(a1) between flood starting-point detection and the highest reservoir level of permission, several water is selected by equidistant principle Position；

(a2) one of them water level selected, calculates different CB corresponding the letting out of structure example under this water level successively Mighty torrent amount；

(a3) verify whether whether all CB structure examples meet flood discharge flow with a combination thereof sequence number monotonic increase Rule, if met; by checking, if being unsatisfactory for, the letter that wherein must have a CB structure example is described There is mistake in breath, need to verify CB structure sample order or the correctness of CB structure example information, for asking Topic amendment CB structure sample order or CB structure example information are till passing through checking；

(a4) if not traveled through all water levels, then jumping to step (a2), if having traveled through all water levels, then tying Bundle checking.

In step (3), the gatage combination sequence of optimization uses following step to generate:

(b1) according to the reservoir level at the beginning of the period, various combination sequence number CB structure example is calculated under this reservoir level Actual flood discharge flow, and record actual flood discharge flow and the expection period optimizes the absolute difference between flood discharge flow Or difference square；

(b2) using square minimum CB structure example of absolute difference or difference as period optimizing decision, will A combination thereof sequence number is as the optimum combination sequence number of this period；

(b3) according to reservoir level at the beginning of the period and this optimum combination sequence number, the optimum flood discharge flow of period is calculated；

(b4) water balance equation is utilized to calculate the reservoir level at the beginning of the next period；

(b5) whether having traveled through all periods, if not traveled through all periods, then having gone to step (b1), if Travel through all periods, then terminate ergodic process；

(b6) the multiple optimum gatage combination number obtained is kept to sequence, the lock of optimization can be formed Door aperture combination sequence.

In step (4), gatage combination sequence optimisation adjusts and includes the adjustment that scheduling process reservoir level is out-of-limit Adjustment with scheduling end of term reservoir level；

The out-of-limit method of adjustment of described scheduling process reservoir level is as follows:

(c1) required precision and the permission requirement of gate regulation frequency of Dispatching Flood process control are taken into account, really Segment length when fixed rational gatage combination sequence optimisation adjusts, and using the first period as present period；

(c2) judge the highest limit of reservoir level of period Mo, if the highest limit, be directly entered the next stage；

(c3) if the highest limit of period end reservoir level, recent operating gate is found the most forward；

(c4) judge that this gate adjusts to open a sluice gate or barrier gate, if gatage combination sequence number and front period phase Than increasing, then illustrate that this operating gate is for opening a sluice gate operation；If gatage combination sequence number subtracted compared with the front period Little, then illustrate that this operating gate is barrier gate operation；

(c5) if opening a sluice gate operation, then the moment will should be opened a sluice gate in advance.If gatage combination sequence optimisation is adjusted A length of Δ T ' of whole period, the former operating gate moment is t_r, then suppose that the operating gate moment is t_r-Δ T ', and from t_rProceed by the beginning of period residing for the-Δ T ' moment and calculate by period water level, if this period water level is the most out-of-limit, then Assuming that the operating gate moment is t_r-2 Δ T ', the like till present period water level is the most out-of-limit.If closing Lock operates, then should be postponed in the barrier gate moment.Assuming that the operating gate moment is t_r+ Δ T ', and from t_r+ Δ T ' the moment Proceed by the beginning of the residing period and calculate by period water level, if this period water level is the most out-of-limit, then suppose operating gate Moment is t_r+ 2 Δ T ', the like till present period water level is the most out-of-limit；

(c6) if in step (c5) the operating gate moment moved forward to an operating gate moment earlier or after Move to a slower operating gate moment, then should by earlier or slower gatage combination sequence number delete, If the operating gate moment moved forward to schedule periods just or after move to present period end, then terminate the water level of present period Out-of-limit adjustment, enters subsequent period and goes to step (c2)；

(c7) if not traveled through all periods, then enter subsequent period and go to step (c2), if having traveled through institute There is the period, then exit the out-of-limit adjustment of scheduling process water level.

The method of adjustment general principle adjustment side out-of-limit with scheduling process reservoir level of described scheduling end of term reservoir level Method is similar to, and specifically comprises the following steps that

(d1) according to the deviation requirement of scheduling end of term reservoir level with target water level, to gatage combination sequence Optimize and revise period length and carry out suitable encryption, the time segment length Δ T ' ' adjusted as scheduling end of term reservoir level；

(d2) the comprehensive descision scheduling out-of-limit situation of end of term reservoir level and last operating gate situation, just determining True gate adjustable strategies.If scheduling end of term water level is barrier gate behaviour higher than target water level and last operating gate Make, or the water level in the scheduling end of term less than target water level and last operating gate for opening a sluice gate to operate, then should Postpone the last operating gate moment.If scheduling end of term water level is higher than target water level and last operating gate For opening a sluice gate operation, or the water level in the scheduling end of term is barrier gate operation less than target water level and last operating gate, Then should shift to an earlier date the last operating gate moment；

(d3) suppose that the last operating gate moment is t_r.If the last operating gate moment should be postponed, First find the period at last operating gate moment place, then from the beginning of at the beginning of this period, suppose respectively A rear operating gate moment is t_r+ΔT′′,t_r+ 2 Δ T ' ' ..., until the scheduling end of term, and utilize the water yield to put down Weighing apparatus principle carries out calculating by period water level, obtains the scheduling end of term water level corresponding to above-mentioned each operating gate moment, Find wherein minimum with target water level deviation scheduling end of term water level and corresponding operating gate moment, be excellent Change the operating gate moment after adjusting；

(d4) if should be shifted to an earlier date, when first finding second from the bottom operating gate the last operating gate moment Carve place period, then from the beginning of at the beginning of this period (if the most for once operating gate, then from first time Start at the beginning of Duan), suppose that the last operating gate moment is t respectively_r-ΔT′′,t_r-2 Δ T ' ' ..., until should Till at the beginning of period, and principle of water balance is utilized to carry out, by the calculating of period water level, obtaining corresponding to above-mentioned each lock The scheduling end of term water level in door operation moment, find wherein minimum with target water level deviation scheduling end of term water level and In the corresponding operating gate moment, it is the operating gate moment after optimizing and revising.

In step (5), operating gate plan uses following methods to generate:

The gatage combination sequence number of adjacent two periods in gatage combination sequence is optimized by period contrast, And by the contrast of two gatage combined informations is resolved, draw the operating gate plan of this period.Traversal The operating gate plan of each period is also all saved in a sequence by all periods, forms the lock in schedule periods Door operation plan.

If the gatage combination sequence number of adjacent two periods is identical, then this period does not carry out operating gate；After if One period aperture combination sequence number combines sequence number more than previous period aperture, and this period opens a sluice gate operation；After if One period aperture combination sequence number combines sequence number less than previous period aperture, and this period is carried out barrier gate operation.

The present invention considers permission aperture, the opening sequence of gate of hydropower station in Optimum Flood Dispatching Model, opens The operation rules such as degree combination, it is possible to directly output each gate optimization operation planning during flood, significantly Reduce workload and the work difficulty of dispatcher；Compared with conventional solution, the present invention is calculated Power station optimizes flood discharge discharge curve and more conforms to the actual discharge capacity characteristic in power station, it is possible to more accurate Reservoir level and the Evolution of downstream flood discharge flow during ground prediction Dispatching Flood, auxiliary dispatching personnel are to flood Risk during scheduling is preferably studied and judged.

Accompanying drawing explanation

Fig. 1 is the power station Flood Dispatching Optimization method overall flow figure of the present invention；

Fig. 2 is the gate flood discharge flow of stepped change；

Fig. 3 is tri-kinds of structures of OM, GT, CB and interrelated graph of a relation；

Fig. 4 is that OM structure, GT structure and CB structure all use text structure block to represent；

Fig. 5 is operating gate rule example table；

Fig. 6 is the classification of correctness verification；

Fig. 7 is aperture combination number change example；

Fig. 8 is that operating gate plan generates.

Detailed description of the invention

For the technological means making the present invention realize, creation characteristic, reach purpose and be easy to understand with effect, under Face combines detailed description of the invention, and the present invention is expanded on further.

Seeing Fig. 1, the Flood Dispatching Optimization method of the present invention includes the description of operating gate rule structuralisation, gate Operation rules instantiation, gatage combination sequence optimisation generate, gatage combination sequence optimisation adjusts, Operating gate plan generates five steps.

Wherein, operating gate rule structuralisation describe have employed structurized description method come complete expression gate permit Permitted the complex informations such as aperture, opening sequence and aperture combination, for operating gate rule instantiation and machine solution Analysis operating gate rule provides basis.

Operating gate rule instantiation allows aperture, opening sequence and aperture according to the true gate in actual power station The information such as combination, dynamically generate and arrange multiple realities of all kinds of operating gate regular texture body according to certain rule Example, and set up the incidence relation between different instances.

Gatage combination sequence optimisation generates the theoretical flood discharge inquired into according to tradition Flood Dispatching Optimization method Discharge curve, calculates optimum aperture combination sequence number by the period, forms the optimization gatage in whole schedule periods Combination sequence, for solving the coupled problem of operating gate rule and tradition Flood Dispatching Optimization method.

Gatage combination sequence optimisation adjusts the period out-of-limit for reservoir level during Dispatching Flood, by referring to Method of determining adjusts the previous operating gate moment, until eliminating the out-of-limit problem of reservoir level；For scheduling end of term reservoir level With the deviation situation of target water level and previous operating gate situation, adjust last operating gate by designation method Moment, until scheduling end of term water level is closest.

Operating gate plan generates by optimizing resolving by the period of gatage combination sequence, dynamically generating Meet the operating gate plan that all operations rule requires.

As can be seen here, only gatage combination sequence optimisation generating portion belongs to tradition Flood Dispatching Optimization side The improvement of method, other four parts are by present invention firstly provides, excellent for supplementing perfect tradition flood further Change dispatching method system.

Operating gate rule structuralisation is described as remaining four part provides basic specification to support, and remaining four part Then constitute the entire protocol of power station Flood Dispatching Optimization method.

(1) operating gate rule structuralisation describes

Gate of hydropower station operation rules mainly includes the permission aperture of gate, opening sequence and aperture combination.Allow Aperture refers to the permission open position of each gate, can use the height (absolute aperture) of the actual unlatching of gate or real Border Lift represents with the ratio (relative opening degree) of maximum Lift.Actual during Dispatching Flood choose Gatage be necessary for one of several discrete apertures of allowing.Therefore, gatage can only strengthen step by step, Corresponding flood discharge flow is step-like change, it is impossible to accurately follow the tracks of expection flood discharge discharge process.Gate allows aperture It it is the principal element (see figure 2) of tradition Flood Dispatching Optimization method execution link flood discharge flow deviation.Owing to permitting Perhaps aperture is all for specific gate, and therefore its structuring expression is fairly simple, it is only necessary to by aperture type (absolutely To aperture, relative opening degree) and the set of permission aperture, as 0,0.25,0.5,0.75,1}, it is saved in In corresponding gate information.

Part gate of hydropower station scheduling rule has simply arranged the priority opening sequence between multiple gate, and advises Could open, after determining that gate standard-sized sheet that priority is the highest, the gate that priority is taken second place.When needs starting gate is let out Hong Shi, that gate that first operator precedence level is the highest, it is gradually increased along with flood into reservoir flow and strengthens this gate Aperture.When this gate reach fully open position and flood into reservoir continue to increase time, then choose that priority is taken second place one Individual or multiple gates, and in time increase or reduce gatage.In the flood water-break stage, then use contrary Regulation order, the gate opened the most afterwards first closes, and closes after the gate first opened.For this situation, it is only necessary to adopt Structured storage requirement can be met with " priority sequence ".

But, also there is part gate of hydropower station scheduling rule to require to regulate in turn between multiple gates, to optimize The HYDRODYNAMIC CONDITION RELATING TO in downstream, power station.Such as, some hydropower station has on the left of A(), in the middle of B(), on the right side of C() Three gates, it is desirable to along with being gradually increased of flood into reservoir flow, first open B gate to 0.25 aperture and 0.5 Aperture, then keeps B gate in 0.5 aperture, opens two gates of A, C to 0.25 aperture and 0.5 simultaneously Aperture, when flood discharge flow increases further, increasing B gate is to 0.75 aperture and standard-sized sheet again, the most again Strengthen A, C gate successively to 0.75 aperture and standard-sized sheet.In this case, can be due to low priority gate Open the most in advance before high priority gate standard-sized sheet, in the case of above-mentioned " priority sequence " method cannot meet this Gate opening sequential organization memory requirement.

For solving the problems referred to above, the present invention stores different gatage combined situation, and is opened by different gates The ordered arrangement of degree combined information, implies the complicated gate opening order information of storage, it is achieved complicated gate opening The indirect statement of order.The method is not only able to solve well the complicated opening sequence of above-mentioned regulating gate in turn Statement problem, and also be able to the most compatible and substitute simple " priority sequence " expression method.Additionally, The method also achieves the integration storage of different gate opening order and aperture combined information, and the information of reducing is deposited The redundancy of storage.

According to above-mentioned analysis, see Fig. 3, be respectively adopted tri-kinds of structures of OM, GT, CB and carry out characterized gate Permission aperture, opening sequence and aperture combined information.OM structure is for storing the letter of the specific aperture of gate Breath, content includes that aperture sequence number, opening value, aperture type, aperture describe, and by several water levels, stream The overflow performance curve that amount binary information group is constituted.GT structure is used for storing the information of the specific gate in power station, Content includes that gate sequence number, gate describe and allow the sequence number of aperture, by aperture sequence number and OM structure Association.CB structure is used for storing specific gatage combined information, and content includes number of regulation, and The aperture combined information being made up of several gate sequence numbers, aperture sequence number binary information group, by gate sequence number and Aperture sequence number associates with GT structure.This mode by opening information with gate unpack storage is permissible The multiple gates avoiding pattern the most identical with Installation Elevation re-define identical opening information so that redundancy Minimize.

(2) operating gate rule instantiation

Operating gate rule structuralisation describes the statement problem only solving different types of information.For specific water power For standing, need according to its concrete gate quantity, the permission aperture of each gate and aperture combined situation, Different structure body is instantiated, in order to utilize structural representation method to describe all information in specific power station.

Due to different power stations gate quantity, allow the combination of aperture, aperture the most different, the most all kinds of knots The quantity of structure body instantiation is the most different.The quantity of OM structure example depends on the quantity of gate of hydropower station With the quantity allowing aperture.In principle, it should all permission apertures of each gate are all set up an OM knot Structure body example.But, when there is multiple gate in power station, and it allows the overflow performance under aperture and corresponding aperture The while that curve being homogeneous, the OM structure example of identical aperture can not be re-defined.With above-mentioned A, B, C tri- As a example by individual gate, if three gate patterns, Installation Elevation are the most identical, it is allowed to aperture be 0.25,0.5, 0.75,1 }, and the overflow performance curve under arbitrary aperture is the most identical, and actual having only to defines 4 OM structure example, opening value is respectively 0.25,0.5,0.75,1.The quantity of GT structure example and water The quantity of power station gate is equal, is respectively directed to each gate and instantiates.The quantity of CB structure example is then The quantity of gatage each with power station combination is equal, is respectively directed to the combination of each aperture and instantiates.

In order to eliminate the dependence to commercial data library management system (DBMS), it is achieved light-weighted information configuration, The present invention uses the block structure of textual form to represent operating gate Rule Information, and computer program can resolve these Block structure, and instantiate corresponding structure.It should be noted that the stage discharge relation in OM example is bent Line should arrange in strict accordance with water level progressive law, in order to according to reservoir level Fast Interpolation meter in solution procedure Calculate the flood discharge flow that this aperture is corresponding.Still as a example by tri-gates of above-mentioned A, B, C, its operating gate rule list State as shown in Figure 4.

The structure example of three types is all loaded in calculator memory by the mode using object array, respectively Setting up aperture table, gate table and combination table as shown in Figure 5, each member in each table is a structure Body example.Aperture table and gate table to the order no requirement (NR) between each example, but in order to the most anti- Mirror aperture Changing Pattern and gate situation, preferably deposit in order according to certain rule.Such as, the knot in aperture table Structure body example can arrange from small to large according to opening value.Structure example in gate table can be according to gate Type and physical location arrangement.Each example in combination table then should in strict accordance with gate opening sequencing and Combined situation carries out ordered arrangement, in order to follow the rational operating gate plan of operating gate Rulemaking.

As given reservoir level z, the flood discharge flow that CB example is corresponding can uniquely determine.Method is that traversal CB is real All aperture sequence numbers in example, find the OM example that this aperture sequence number is corresponding respectively, then travel through this OM example In water level, flow information group, find respectively less than or equal to and more than or equal to two adjacent water level values of z, note For z ' and z ' ', corresponding flood discharge flow is respectively q ' and q ' ', then when reservoir level is z, the flow of this OM example can It is expressed as q_om=q′+(q′′-q′)·(z-z′)/(z′′-z′).If not finding z ' and z ' ', then should when reservoir level is z The flood discharge flow of OM example is zero.Add up the flood discharge flow of each OM example successively in CB example, obtains total Flood discharge flow, is the flood discharge flow of this CB example.Therefore, flood discharge flow is reservoir level and aperture combination sequence Number function, be designated as q_cb=f_rq(z, c), wherein z is reservoir level, and c is that aperture combines sequence number, q_cbFor this lock The total flood discharge in power station flow under door aperture combination.

After all structure examples are loaded into internal memory, reply example information carries out correctness verification, is divided into sequence number school Test, associate verification and data verification three kinds, see Fig. 6.Wherein, the sequence number verification aperture sequence to each OM example Number, the combination sequence number of the gate sequence number of GT example and CB example verifies.Require that each sequence number must be with certainly So Number Sequence is incremented by, and must not duplicate or jump.Association verification is for verifying that the relevance between each example is No correctly, including the aperture sequence number in gate table whether exist in aperture table correspondence OM example, combination table In gate sequence number in gate table, whether there is the GT example of correspondence, and whether the aperture sequence number of correspondence deposit It is in this GT example.Whether data verification exists mistake for verifying data data self, and this mistake can Firsthand information can be come from, it is also possible to come from the human error of Data Input Process.Verification content mainly includes Whether the rating curve of OM example reasonable, and combination table each CB example and order the most reasonable. According to the principle verification OM example that reservoir level the highest flood discharge flow under same aperture is the biggest, contrast OM example one by one Each water level, flow binary information group, check whether water level and flow are monotonic increase, are unsatisfactory for requiring, sentence There is mistake in this OM example fixed.Combination table stores and increases along with flood into reservoir and strengthen starting gate step by step Aperture combined information, under thus like reservoir level, smaller group is closed the letdown flow of sequence number CB example and should be less than The letdown flow of bigger combination sequence number CB example.Based on this principle, following steps can be taked to check, and combination table is each CB example and order are the most reasonable:

(a1) between flood starting-point detection and the highest reservoir level of permission, several water is selected by equidistant principle Position；

(a2) one of them water level selected, calculates the flood discharge flow that under this water level, different CB examples are corresponding successively；

(a3) verify all CB structure examples whether meet flood discharge flow with combination sequence number monotonically increasing rule, If meeting, by checking, if being unsatisfactory for, illustrate wherein to need with the presence of the information mistake of a CB example Verification CB structure sample order or the correctness of CB structure example information, revise CB structure for problem Body sample order or CB structure example information are till passing through checking；

(a4) if not traveled through all water levels, then step (a2) is jumped to, if having traveled through all water levels, then Terminate checking.

(3) gatage combination sequence optimisation generates

Gate of hydropower station needs to follow certain opening sequence and aperture combination requirement, and this two category information is all protected Exist in combination table.Single CB example represents a kind of specific gatage combination, and between CB example Ordered arrangement, it can be ensured that the operating gate scheme that Optimization Solution obtains meets opening sequence requirement all the time.? In the case of correct structure aperture table, gate table and combination table, it is only necessary to search in the combination table of ordered arrangement The optimum CB sequence number of different periods, can consider the permission aperture of gate, opening sequence and aperture combination comprehensively Constraint, efficiently solves the coupled problem that operating gate rule calculates with Flood Dispatching Optimization.

Increase very additionally, power station allows to cross over regulation, such as flood between different gatage combinations Time fast, the gate being in closed positions can be crossed over 0.25 aperture and directly open to 0.5 aperture, so arbitrarily Redirecting between two combination sequence numbers, all can meet gate opening order naturally and aperture combines requirement.Therefore, There is not interconnection constraint between the combination sequence number of adjacent time interval, the optimizing space of gatage combination sequence number is whole Individual combination table so that the optimizing operation of aperture combination sequence number has obtained great simplification.

As it was previously stated, power station discharge capacity is constrained to q in tradition Flood Dispatching Optimization method_i≤q_max(z_i), It is not intended that the operation rules of gate during this so that model fails truly to reflect the operating gate rule in power station Then retrain.By the structural representation of operating gate rule, power station flood discharge flow is represented by reservoir level and opens The function of degree combination sequence number, i.e. q_cb=f_rq(z,c).Owing to each CB example in combination table is opened according to gate Degree situation ordered arrangement, therefore available maximum combination sequence number tries to achieve the power station earial drainage energy under particular bin water level Power.So, the power station discharge capacity constraint in tradition Flood Dispatching Optimization method can be revised as q_i≤f_rq(z_i,c_max), wherein c_maxMaximum opening for power station combines sequence number.Utilize this amended mathematics Model and tradition method for solving, try to achieve intended optimization flood discharge discharge curve.

Using above-mentioned optimization flood discharge discharge curve as benchmark, travel through each period successively, find actual flood discharge Flow and expection optimize flood discharge flow immediate aperture combination sequence number.The Mathematical Modeling of the i-th period is as follows:

Object function:

$\mathrm{Min}\left[{(q_i^j-q_i)}^2\right],j=\mathrm{0,1,2},L,R$

Constraints:

Actual letdown flow retrains $q_i^j=f_{\mathrm{rq}}(z_i,r_i^j),j=\mathrm{1,2},L,R;j=\mathrm{1,2},L,N$

Traffic constraints during without letting out in advance $q_i^j\≤Q_{\max },j=\mathrm{1,2},L,R;j=\mathrm{1,2},L,N$

Reservoir water yield Constraints of Equilibrium $V_{i+1}^j=V_i+(Q_i-q_i^j)\·\mathrm{\Δt},i=\mathrm{1,2},L,N$

Water level storage capacity transformational relation $z_{i+1}=f_{\mathrm{zv}}\left(V_{i+1}^j\right)$

Reservoir peak level retrains $z_{i+1}^j\≤f_{\mathrm{vz}}(V_0+V_f),j=\mathrm{1,2},L,R;i=\mathrm{1,2},L,N$

In formula, q_iFor the optimization flood discharge flow of i-th period,Combine for i-th choosing period of time jth aperture Number flood discharge flow,For the combination number of jth gatage, R is gatage combination quantity, f_rqIt is Functional relation between flood discharge flow and reservoir level, gatage combination number, f_vzLetter for reservoir level Yu storage capacity Number relation,For the period Mo (at the beginning of the i+1 period) in the case of the combination number of i-th choosing period of time jth aperture Storage capacity,For period Mo (at the beginning of the i+1 period) storehouse in the case of the combination number of i-th choosing period of time jth aperture Water level.If object function is minimum during j=0, then this period should be by all gate Close Alls.Due under reality Vent flow constraint has reflected the gatage restriction effect to letdown flow, and therefore this model is without the concern for Big discharge capacity constraint.

During the period optimizes, there will be traffic constraints and peak level retrains situation about cannot simultaneously meet.If Keep storage outflow less than reservoir inflow, then the reservoir level of period Mo can exceed high limit.If keeping period end water level Below height limits, then must assure that storage outflow is more than reservoir inflow.Now, it should preferentially guarantee storage outflow Less than reservoir inflow, it is to avoid cause artificial flood, and by follow-up operating gate planning optimization method of adjustment solution The certainly out-of-limit problem of water level.

The concrete solution procedure using above-mentioned Mathematical Modeling solving-optimizing gatage combination sequence is as follows:

(b1) according to the reservoir level at the beginning of the period, various combination sequence number CB structure example is calculated under this reservoir level Actual flood discharge flow, and record actual flood discharge flow and the expection period optimizes the absolute difference between flood discharge flow Or difference square；

(b2) using square minimum CB structure example of absolute difference or difference as period optimizing decision, will A combination thereof sequence number is as the optimum combination sequence number of this period；

(b3) according to reservoir level at the beginning of the period and this optimum combination sequence number, the optimum flood discharge flow of period is calculated；

(b4) water balance equation is utilized to calculate the reservoir level at the beginning of the next period；

(b5) whether having traveled through all periods, if not traveled through all periods, then having gone to step (b1), if Travel through all periods, then terminate ergodic process.

(b6) the N number of optimum gatage combination number obtained is kept to sequence, the lock of optimization can be formed Door aperture combination sequence.

(4) gatage combination sequence optimisation adjusts

As it was previously stated, optimize gatage combination sequence only pursue the gate flood discharge flow in each stage closest to Intended optimization flood discharge flow, it is impossible to the problem that reservoir level transfinites occurs during avoiding Dispatching Flood.Additionally, The impact on reservoir risk of the later stage water situation, and the need of flood resource Optimum utilization are considered in part power station Ask, may require that Dispatching Flood end of term reservoir level returns back to designated water level.Gatage combination sequence optimisation adjusts The above-mentioned out-of-limit problem of scheduling process reservoir level and scheduling end of term water level adjustment problem is solved for simultaneously.

Seeing table 1, the out-of-limit method of adjustment of scheduling process water level is as follows:

(c1) required precision and the permission requirement of gate regulation frequency of Dispatching Flood process control are taken into account, really When fixed rational gatage combination sequence optimisation adjusts, segment length Δ T ', Δ T ' preferably take 15 minutes or 10 minutes, Using the first period as present period；

(c2) judge the highest limit of reservoir level of period Mo, if the highest limit, be directly entered the next stage；

(c3) if the highest limit of period end reservoir level, recent operating gate (gatage is found the most forward The period that combination sequence number changes)；

(c4) judge that this gate adjusts to open a sluice gate or barrier gate.If gatage combination sequence number and front period phase Than increasing, then illustrate that this operating gate is for opening a sluice gate operation；If gatage combination sequence number subtracted compared with the front period Little, then illustrate that this gate is barrier gate operation；

(c5) if opening a sluice gate operation, then the moment will should be opened a sluice gate in advance.If the former operating gate moment is t_r, then suppose The operating gate moment is t_r-Δ T ', and from t_rProceed by by period water-level gauge at the beginning of period residing for the-Δ T ' moment Calculate, if this period water level is the most out-of-limit, then suppose that the operating gate moment is t_r-2 Δ T ', the like until current Till period water level is the most out-of-limit.If barrier gate operates, then the barrier gate moment should be postponed.Assuming that the operating gate moment For t_r+ Δ T ', and from t_rProceed by the beginning of period residing for the+Δ T ' moment and calculate by period water level, if this period Water level is the most out-of-limit, then suppose that the operating gate moment is t_r+ 2 Δ T ', the like until present period water level is not got over Only it is limited to.

(c6) if in step (c5) the operating gate moment moved forward to an operating gate moment earlier or after Move to a slower operating gate moment, then should by earlier or slower gatage combination sequence number delete, If the operating gate moment moved forward to schedule periods just or after move to present period end, then terminate the water level of present period Out-of-limit adjustment, enters subsequent period and goes to step (c2)；

(c7) if not traveled through all periods, then enter subsequent period and go to step (c2)；If having traveled through All periods, then exit the out-of-limit adjustment of scheduling process water level.

Table 1: the method for adjustment that scheduling process middle water level is out-of-limit

Sequence number	Previous operating gate situation	Gate adjustable strategies
			1	Barrier gate	Postpone the barrier gate time
2	Open a sluice gate	Switching-off time in advance

See table 2, dispatch the tune that the method for adjustment general principle of end of term reservoir level is out-of-limit with scheduling process reservoir level Adjusting method is similar to, and specifically comprises the following steps that

(d1) according to the deviation requirement of scheduling end of term reservoir level with target water level, to gatage combination sequence Optimize and revise period length and carry out suitable encryption, the time segment length Δ T ' ' adjusted as scheduling end of term reservoir level；

(d2) the comprehensive descision scheduling out-of-limit situation of end of term reservoir level and last operating gate situation, just determining True gate adjustable strategies.If scheduling end of term water level is barrier gate behaviour higher than target water level and last operating gate Make, or the water level in the scheduling end of term less than target water level and last operating gate for opening a sluice gate to operate, then should Postpone the last operating gate moment.If scheduling end of term water level is higher than target water level and last operating gate For opening a sluice gate operation, or the water level in the scheduling end of term is barrier gate operation less than target water level and last operating gate, Then should shift to an earlier date the last operating gate moment；

(d3) suppose that the last operating gate moment is t_r.If the last operating gate moment should be postponed, First find the period at last operating gate moment place, then from the beginning of at the beginning of this period, suppose respectively A rear operating gate moment is t_r+ΔT′,t_r+ 2 Δ T ' ..., until the scheduling end of term, and utilize water balance Principle carries out calculating by period water level, obtains the scheduling end of term water level corresponding to above-mentioned each operating gate moment, Find wherein with the scheduling end of term water level and corresponding operating gate moment of target water level deviation, be optimization and adjust The operating gate moment after whole；

(d4) if should be shifted to an earlier date, when first finding second from the bottom operating gate the last operating gate moment Carve place period, then from the beginning of at the beginning of this period (if the most for once operating gate, then from first time Start at the beginning of Duan), suppose that the last operating gate moment is t respectively_r-ΔT′,t_r-2 Δ T ' ..., until this time Till at the beginning of Duan, and principle of water balance is utilized to carry out, by the calculating of period water level, obtaining corresponding to above-mentioned each gate The scheduling end of term water level in operation moment, find wherein with the scheduling end of term water level and corresponding of target water level deviation In the operating gate moment, it is the operating gate moment after optimizing and revising.

Table 2: scheduling end of term water level method of adjustment

(5) operating gate plan generates

This part is for generating operating gate plan according to the optimization gatage combination sequence obtained.Substantially Thought is the aperture combination number by period contrast adjacent two periods, by two gatage combined informations Contrast resolves, and draws the operating gate plan of this period.Travel through all periods and by the operating gate of each period Plan is all saved in a sequence, is formed and meets the operating gate plan that all operating gate rules require.

If the gatage combination sequence number of adjacent two periods is identical, then this period does not carry out operating gate.If phase The gatage combination sequence number of adjacent two periods is different, then extract corresponding 2 operating gate rule example, and The operation scheme of period gate is obtained by the comparison of content.Owing to each CB example in combination table is according to lock Door unlatching situation arranges from small to large, therefore can be simply by the gatage combination number of adjacent two periods Size, it is judged that be to be by this period opening a sluice gate to operate also being by barrier gate operation.A period aperture combination after if Sequence number combines sequence number more than previous period aperture, then this period carries out opening a sluice gate operation；Otherwise, then this period is carried out Barrier gate operates.

For convenience of description, still as a example by the aperture shown in Fig. 4 combines.Consider respectively to open a sluice gate and two kinds of feelings of barrier gate Condition, the situation of opening a sluice gate suppose former and later two be the optimization aperture combination of period number be to be 1 and 3 respectively, barrier gate situation Assuming that former and later two are the optimization aperture combinations of period number is to be 3 and 1 respectively, extraction gate aperture combination respectively Information is as follows:

The gatage combination number of Fig. 7 left part is become 3 from 1, therefore for opening a sluice gate operation.By to two The comparison analysis of CB example content can draw, gate 2(gate B) need from aperture 1(0.25 aperture) Regulation is to aperture 2(0.5 aperture), gate 1(gate A) and gate 3(gate C) need from full off state It is opened into aperture 1(0.25 aperture).

In like manner, the gatage combination number of Fig. 7 right part is become 1 from 3, therefore operates for barrier gate.Pass through Comparison analysis to two CB example content can obtain, gate 2(gate B) need from aperture 2(0.5 Aperture) regulation to aperture 1(0.25 aperture), gate 1(gate A) and gate 3(gate C) needs from Aperture 1(0.25 aperture) regulate to closed positions.See Fig. 8.

Knowable to the method that gatage combination sequence optimisation generates, the regulating time of above-mentioned gate all should be this Moment value at the beginning of period.The information of gate regulation can be by gate title, regulation moment, target aperture three metamessage Group describes, when certain period has multiple gate to need to be adjusted, when can be constituted this by multiple three metamessage groups The operating gate plan of section.Recycling said method travels through all periods, can obtain the gate of each period Operation planning, collectively forms the operating gate plan in whole schedule periods.

The general principle of the present invention and principal character and advantages of the present invention have more than been shown and described.The industry Skilled person will appreciate that, the present invention is not restricted to the described embodiments, described in above-described embodiment and specification The principle that the present invention is simply described, without departing from the spirit and scope of the present invention, the present invention also has Various changes and modifications, these changes and improvements both fall within scope of the claimed invention.Application claims Protection domain is defined by appending claims and equivalent thereof.

Claims (8)

1. the power station Flood Dispatching Optimization method considering operating gate rule, it is characterised in that successively Including following step:

(1) operating gate rule structuralisation describes, and uses structurized description method to carry out complete expression gate and allows Aperture, opening sequence and aperture combined information, grasp for operating gate rule instantiation and computer analyzing gate Make rule and basis is provided；

(2) operating gate rule instantiation is carried out for specific power station, according to the true gate in specific power station Allow aperture, opening sequence and aperture combined information, dynamically generate and arrange all kinds of operating gate regular texture body Multiple examples, and set up the incidence relation between different instances；

(3) the gatage combination sequence optimized, the theory inquired into according to traditional dynamic programming algorithm are generated Flood discharge discharge curve, calculates optimum aperture combination sequence number by the period, forms the optimization gate in whole schedule periods Aperture combination sequence, for solving the coupled problem of operating gate rule and tradition Flood Dispatching Optimization method；

(4) gatage combination sequence is optimized adjustment, out-of-limit for reservoir level during Dispatching Flood Period, adjust the previous operating gate moment, until eliminate the out-of-limit problem of reservoir level；For scheduling end of term storehouse water Position and the deviation situation of target water level and previous operating gate situation, adjust the last operating gate moment, directly Closest to scheduling end of term water level；

(5) operating gate plan is generated, by optimizing resolving by the period of gatage combination sequence, dynamic State generates and meets the operating gate plan that all operating gate rules require；

In step (1), it is as follows that operating gate rule structuralisation describes method:

The permission aperture of gate, opening sequence and aperture combined information are respectively adopted OM structure, GT structure Body and CB structure characterize wherein；

OM structure is for storing the information of the specific aperture of gate, and content includes aperture sequence number, opening value, opens The rating curve that degree type, aperture describe and be made up of several water levels, flow binary information group；

GT structure is for storing the information of the specific gate in power station, and content includes that gate sequence number, gate describe With the aperture sequence number allowed, GT structure is associated with OM structure by aperture sequence number；

CB structure is used for storing specific gatage combined information, and content includes number of regulation and by some The aperture combined information that individual gate sequence number, aperture sequence number binary information group are constituted, CB structure passes through gate sequence Number associate with GT structure with aperture sequence number.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 1 rule, It is characterized in that,

In step (2), operating gate rule is instantiated, refer to allow for the gate in actual power station Aperture, opening sequence and operation rules information, respectively instantiation OM structure, GT structure and CB structure Body, and form the aperture table for actual power station, gate table and combination table, specifically includes following steps:

First, allowing aperture according to every kind of gate in actual power station, the OM structure building correspondence respectively is real Example, states this gate and allows opening information, and all of OM structure example collectively forms complete aperture table；

Secondly, according to each floodgate in actual power station, build the GT structure example of correspondence respectively, Stating this gate information, all of GT structure example collectively forms complete gate table；

Finally, according to every kind of gatage combination in actual power station, the CB structure building correspondence respectively is real Example, states this aperture combined information, and all of CB structure collectively forms complete combination table.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 2 rule, It is characterized in that,

Rating curve in each described OM structure example should be arranged in strict accordance with water level progressive law Row, in order to calculate, according to reservoir level Fast Interpolation, the flood discharge flow that this aperture is corresponding in solution procedure；

Each OM structure example in described aperture table can arrange from small to large according to opening value；Described gate Each GT structure example in table can arrange according to the type of gate and physical location；In described combination table Each CB structure example then should enter in strict accordance with the order that gate opening sequencing and aperture combine from small to large Row arrangement.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 3 rule, It is characterized in that, between described combination table each CB structure example, order is the most correct to take following methods to check:

(a1) between flood starting-point detection and the highest reservoir level of permission, several water is selected by equidistant principle Position；

(a2) one of them water level selected, calculates different CB corresponding the letting out of structure example under this water level successively Mighty torrent amount；

(a3) verify whether all CB structure examples meet flood discharge flow with combination sequence number monotonically increasing rule Rule, if meeting, by checking, if being unsatisfactory for, illustrates that the information that wherein must have a CB structure example is deposited In mistake, CB structure sample order or the correctness of CB structure example information need to be verified, repair for problem Change CB structure sample order or CB structure example information till passing through checking；

(a4) if not traveled through all water levels, then step (a2) is jumped to, if having traveled through all water levels, then Terminate checking.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 4 rule, It is characterized in that,

In step (3), the gatage combination sequence of optimization specifically uses following step to generate:

(b1) according to the reservoir level at the beginning of the period, various combination sequence number CB structure example is calculated under this reservoir level Actual flood discharge flow, and record actual flood discharge flow and the expection period optimizes the absolute difference between flood discharge flow Or difference square；

(b2) using square minimum CB structure example of absolute difference or difference as period optimizing decision, will A combination thereof sequence number is as the optimum combination sequence number of this period；

(b3) according to reservoir level at the beginning of the period and this optimum combination sequence number, the optimum flood discharge flow of period is calculated；

(b4) water balance equation is utilized to calculate the reservoir level at the beginning of the next period；

(b5) whether having traveled through all periods, if not traveled through all periods, then having gone to step (b1), if Travel through all periods, then terminate ergodic process；

(b6) the multiple optimum gatage combination number obtained is kept to sequence, the lock of optimization can be formed Door aperture combination sequence.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 5 rule, It is characterized in that,

In step (4), gatage combination sequence optimisation adjusts and includes the adjustment that scheduling process reservoir level is out-of-limit Adjustment with scheduling end of term reservoir level；

The out-of-limit method of adjustment of described scheduling process reservoir level is as follows:

(c1) required precision and the permission requirement of gate regulation frequency of Dispatching Flood process control are taken into account, really Segment length when fixed rational gatage combination sequence optimisation adjusts, and using the first period as present period；

(c2) judge the highest limit of reservoir level of period Mo, if the highest limit, be directly entered the next stage；

(c3) if the highest limit of period end reservoir level, recent operating gate is found the most forward；

(c4) judge that this operating gate is opened a sluice gate or barrier gate, if gatage combination sequence number and front period phase Than increasing, then illustrate that this operating gate is for opening a sluice gate operation；If gatage combination sequence number subtracted compared with the front period Little, then illustrate that this operating gate is barrier gate operation；

(c5) if opening a sluice gate operation, then the moment will should be opened a sluice gate in advance, if operating gate planning optimization adjusts the period A length of Δ T ', the former operating gate moment is t_r, then suppose that the operating gate moment is t_r-Δ T ', and from t_rDuring-Δ T ' Proceed by the beginning of period residing for quarter and calculate by period water level, if this period water level is the most out-of-limit, then suppose gate behaviour It is t as the moment_r-2 Δ T ', the like till present period water level is the most out-of-limit, if barrier gate operation, then The barrier gate moment should be postponed, it is assumed that the operating gate moment is t_r+ Δ T ', and from t_rPeriod residing for the+Δ T ' moment Just proceed by and calculate by period water level, if this period water level is the most out-of-limit, then suppose that the operating gate moment is t_r+ 2 Δ T ', the like till present period water level is the most out-of-limit；

(c6) if in step (c5) the operating gate moment moved forward to an operating gate moment earlier or after Move to a slower operating gate moment, then should by earlier or slower gatage combination sequence number delete, If the operating gate moment moved forward to schedule periods just or after move to present period end, then terminate the water level of present period Out-of-limit adjustment, enters subsequent period and goes to step (c2)；

(c7) if not traveled through all periods, then enter subsequent period and go to step (c2), if having traveled through All periods, then exit the out-of-limit adjustment of scheduling process water level；

The method of adjustment of described scheduling end of term reservoir level is as follows:

(d1) according to the deviation requirement of scheduling end of term reservoir level with target water level, to gatage combination sequence Optimize and revise period length to be encrypted, as the time segment length Δ T of scheduling end of term reservoir level adjustment "；

(d2) the comprehensive descision scheduling out-of-limit situation of end of term reservoir level and last operating gate situation, just determining True gate adjustable strategies；If scheduling end of term water level is barrier gate behaviour higher than target water level and last operating gate Make, or the water level in the scheduling end of term less than target water level and last operating gate for opening a sluice gate to operate, then should Postpone the last operating gate moment；If scheduling end of term water level is higher than target water level and last operating gate For opening a sluice gate operation, or the water level in the scheduling end of term is barrier gate operation less than target water level and last operating gate, Then should shift to an earlier date the last operating gate moment；

(d3) suppose that the last operating gate moment is t_rIf the last operating gate moment should be postponed, First find the period at last operating gate moment place, then from the beginning of at the beginning of this period, suppose respectively A rear operating gate moment is t_r+ΔT″,t_r+ 2 Δ T " ..., until the scheduling end of term, and utilize the water yield to put down Weighing apparatus principle carries out calculating by period water level, obtains the scheduling end of term water level corresponding to above-mentioned each operating gate moment, Find wherein minimum with target water level deviation scheduling end of term water level and corresponding operating gate moment, be excellent Change the operating gate moment after adjusting；

(d4) if should be shifted to an earlier date, when first finding second from the bottom operating gate the last operating gate moment Carve the period at place, then from the beginning of at the beginning of this period, if the most for once operating gate, then from first time Start at the beginning of Duan, suppose that the last operating gate moment is t respectively_r-ΔT″,t_r-2 Δ T " ..., until should Till at the beginning of period, and principle of water balance is utilized to carry out, by the calculating of period water level, obtaining corresponding to above-mentioned each lock The scheduling end of term water level in door operation moment, find wherein minimum with target water level deviation scheduling end of term water level and In the corresponding operating gate moment, it is the operating gate moment after optimizing and revising.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 6 rule, It is characterized in that,

In step (5), operating gate plan uses following methods to generate:

The gatage combination sequence number of adjacent two periods in gatage combination sequence is optimized by period contrast, And by the contrast of two gatage combined informations is resolved, draw the operating gate plan of this period, traversal The operating gate plan of each period is also all saved in a sequence by all periods, forms the lock in schedule periods Door operation plan.

The power station Flood Dispatching Optimization method of consideration operating gate the most according to claim 7 rule, It is characterized in that,

If the gatage combination sequence number of adjacent two periods is identical, then this period does not carry out operating gate；After if One period aperture combination sequence number combines sequence number more than previous period aperture, and this period opens a sluice gate operation；After if One period aperture combination sequence number combines sequence number less than previous period aperture, and this period is carried out barrier gate operation.