CN103646294A - Hydropower-station flood optimization dispatching method taking gate operation rules into consideration - Google Patents

Abstract

The invention discloses a hydropower-station flood optimization dispatching method which takes gate operation rules into consideration. The method includes five steps: (1) structural description of the gate operation rules, which solves an expression problem of the complex gate operation rules; (2) instantiation of the gate operation rules, which solves a problem of how to use the structural-body examples of the gate operation rules to describe specific hydropower-station gate operation rules; (3) optimization generation of a gate opening combination number sequence, which solves a problem of coupling of the gate operation rules and traditional hydropower-station flood optimization dispatching methods; (4) optimization adjustment of the gate opening combination number sequence, which solves problems that in a flood dispatching process, during specific periods, water level exceeds the limit and at the end of a dispatching period, deviation exists between the water level and a target water level; (5) generation of a gate operation plan, which solves a problem of generating the gate operation plan dynamically according to the gate opening combination number sequence. The hydropower-station flood optimization dispatching method which takes the gate operation rules into consideration improves prediction accuracy of a later-period flood-discharge flow process and improves the rapid response capability of a flood prevention decision of a hydropower station.

Description

A kind of power station Flood Dispatching Optimization method of 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 power station Flood Dispatching Optimization method of considering operating gate rule.

Background technology

Traditional power station Flood Dispatching Optimization method is considered as a scheduler object by all gates, only consider the whole discharge capacity of all gates, do not consider the operating gate rules such as permission aperture, opening sequence, aperture combination of each gate, only can calculate the total flood discharge discharge curve of optimization of gate of hydropower station group integral body, and cannot automatically make the Optimum Operation plan of each gate.In actual motion link, need according to this, optimize total flood discharge discharge curve, the operation planning of manually formulating each gate by dispatcher.This mode is difficult to meet power station become more meticulous scheduling and the business demand of rolling scheduling in real time, has also had a strong impact on the efficiency of flood decision simultaneously.In addition, owing to must following the working rule of each gate in actual motion, uncontrollable gate is according to the total flood discharge discharge curve operation of the optimization of expection, cause the deviation between the actual total letdown flow of gate and the total letdown flow of optimization of expection, and then cause later stage actual library water level with expection reservoir level deviation, further aggravate the extent of deviation between actual flood discharge flow and expection flood discharge flow, had a strong impact on the practical value of traditional power station Flood Dispatching Optimization method.

Power station Flood Dispatching Optimization mainly solves the optimization of flood into reservoir and holds the problem of letting out, and is guaranteeing, under the inherently safe prerequisite of power station, to alleviate as far as possible the damaging effect of flood to downstream, and the recycling of flood with due regard to.Tradition power station Flood Dispatching Optimization method is considered as an object by all gates, calculates the optimization flood discharge discharge curve in whole power station by Optimized Operation, does not consider the working rule of each gate in this process.In practical implementation, need to rely on dispatcher's experience and constantly tentative calculation, with reference to the expection having obtained, optimize flood discharge discharge curve and actual reservoir level, determine and meet the operation planning that each gate allows the working rules such as aperture, opening sequence, aperture combination, and the artificial gatage that regulates in real time release floodwatering facility, make actual flood discharge flow approaching as far as possible with expection flood discharge flow.Owing to having ignored the working rule of each gate, the true overflow performance of optimization flood discharge discharge curve and power station that tradition power station Flood Dispatching Optimization method obtains is not inconsistent, thereby cause the substantial deviation of carrying out link and plan link, significantly reduced classic method result of calculation value and can referential.

The final purpose of Flood Dispatching Optimization is by the rational management to flood, makes national economy minimization of loss.Yet the national economy loss that flood causes itself is difficult to quantification calculates, so power station Flood Dispatching Optimization adopts the mode of simplification to consider flood losses conventionally.Generally, national economy loss and maximum being proportionate property of flood discharge flow that flood causes, can be converted into national economy minimization of loss target maximum flood discharge flow and minimize target.Utilize mathematical derivation to prove, while using definite storage capacity in the Dispatching Flood phase, in fact the minimized scheme of maximum flood discharge flow is exactly the most uniform scheme of flood discharge flow.Because initial water level and the last water level of Dispatching Flood are all determined in advance, and the water yield of flood warehouse-in is also determined in advance, therefore according to the principle of water balance, can determine in advance the total flood discharge water yield in schedule periods, this means that the flood discharge flow sum of day part in schedule periods is a fixed value.According to mathematical principle, when a plurality of independent variable sums one timings, can explain its degree of uniformity by the quadratic sum of each independent variable.Therefore, can set up following objective function:

$\mathrm{Min}{\&Integral;}_{t_0}^{t_d}{q}^2\left(t\right)\mathrm{dt},t\&Element;[{\mathrm{<figure-callout\; id="0"\; label="t"\; filenames="HDA0000412922920000021.png,HDA0000412922920000031.png"\; state="\{\{state\}\}">t</figure-callout>}}_0,t_d]$

For the ease of using computer solving, can carry out discretize to above-mentioned objective function:

$\mathrm{Min}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{q}_i^2\&CenterDot;\mathrm{\&Delta;t},i=\mathrm{1,2},L,N;N=\frac{t_d-{\mathrm{<figure-callout\; id="0"\; label="t"\; filenames="HDA0000412922920000021.png,HDA0000412922920000031.png"\; state="\{\{state\}\}">t</figure-callout>}}_0}{\mathrm{\&Delta;t}};$

Constraint condition:

Storage capacity constraint $\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}(Q_i-q_i)\&CenterDot;\mathrm{\&Delta;t}={\mathrm{\&Delta;V}}_f,i=\mathrm{1,2},L,N$

Traffic constraints q when nothing is let 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 _ibe the flood discharge flow of i period, Q _maxfor the maximum reservoir inflow (crest discharge) of flood, segment length when Δ t is, hop count when N is total, Δ V _ffor the storage capacity of intending employing, z _ibe the reservoir level of i period, q _max(z _i) be that all release floodwatering facilities are at z _imaximum discharge capacity under reservoir level, can obtain according to maximum discharge capacity curve interpolation V _ibe the storage capacity at the beginning of i period.

The storage capacity size delta V that plan is employed _fshould be less than or equal to total storage capacity of reservoir.Δ V _fconcrete size should and determine in conjunction with artificial experience according to later stage inflow forecast situation.When the later stage, water was more, should set less Δ V _f, make schedule periods interior with flood discharge water under larger peak clipping flow, so that being reduced to as early as possible flood in flood water-break process, reservoir restricts water supply below position.When the later stage, water was less, can set larger Δ V _f, make schedule periods interior with flood discharge water under less peak clipping flow, reduce downstream flood loss.

According to above-mentioned mathematical model, use traditional dynamic programming algorithm can calculate the optimization flood discharge discharge curve in whole power station.

From above-mentioned mathematical model, can find out, prior art scheme has only been considered the maximum discharge capacity constraint of all release floodwatering facilities in power station, do not consider the impact that the operating gate rules such as permission aperture, opening sequence, aperture combination of each gate produce flood discharge process, only can obtain 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 manually to monitor in real time that reservoir level changes, determine whether open or closed shutter, and calculate corresponding gatage according to the optimization flood discharge flow of Current Library water level and expection, workload is very large.In addition,, because the aperture of release floodwatering facility gate is all confined to several discrete fixedly opening values, actual gate flood discharge flow is optimized flood discharge flow with expection can there is larger deviation.This deviation all exists in whole schedule periods, and the continuous accumulation of deviation makes the optimization flood discharge graph of schedule periods posterior part lose prediction and directive significance.Traditional rolling scheduling method also cannot thoroughly address this problem, make dispatcher accurately anticipation go out the development trend of reservoir level and letdown flow in schedule periods, even can cause the dispatcher who lacks experience to judge flood-proofing terrain by accident, to flood control, bring larger potential safety hazard.

Summary of the invention

The deficiency existing for prior art, the object of the invention is to provide a kind of power station Flood Dispatching Optimization method of considering operating gate rule, solved the problem that classic method actual execution link and plan link exist serious deviation, improved the practice guiding action of power station optimization flood discharge discharge curve, dispatcher's working strength and error probability have been reduced, meanwhile, auxiliary dispatching personnel study and judge later stage flood-proofing terrain variation more accurately.

To achieve these goals, the present invention realizes by the following technical solutions:

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

(1) operating gate rule structuralisation is described, and adopts structurized describing method to come complete expression gate to allow the complex informations such as aperture, opening sequence and aperture combination, for operating gate rule instantiation and computer analyzing operating gate rule provide basis;

(2) for specific power station, carry out the instantiation of operating gate rule, according to the true gate in actual power station, allow the information such as aperture, opening sequence and aperture combination, according to certain rule, dynamically generate and arrange the Multi-instance of all kinds of operating gate regular texture bodies, and set up the incidence relation between different instances;

(3) generate the gatage combination sequence of optimizing, the theoretical flood discharge discharge curve of inquiring into according to traditional power station Flood Dispatching Optimization method, by the period, calculate optimum aperture and combine sequence number, form the optimization gatage combination sequence in whole schedule periods, for solving the coupled problem of operating gate rule and traditional Flood Dispatching Optimization method;

(4) gatage combination sequence is optimized to adjustment, for indivedual out-of-limit periods of reservoir level in Dispatching Flood process, by designation method, adjusts the last time operating gate moment, until eliminate the out-of-limit problem of reservoir level.For the deviation situation of scheduling end of term reservoir level and target water level and operating gate situation last time, by the last operating gate moment of designation method adjustment, until it is the most approaching to dispatch end of term water level;

(5) generate operating gate plan, by optimizing resolving by the period of gatage combination sequence, dynamically generation meets the operating gate plan of all operating gate rule requests.

In step (1), described operating gate rule comprises permission aperture, opening sequence and the aperture combined information of gate, and its structural description method is as follows:

Adopt respectively aperture structure (hereinafter to be referred as OM structure), gate structure body (hereinafter to be referred as GT structure) and aperture unitized construction body (hereinafter to be referred as CB structure) to characterize permission aperture, opening sequence and the aperture combined information of gate.

Described OM structure is for storing the information of the specific aperture of gate, and its content comprises the rating curve that aperture sequence number, opening value, aperture type, aperture are described and consisted of several water levels, flow binary information group;

Described GT structure is for the information of the specific gate in storage of water power station, and its content comprises the aperture sequence number that gate sequence number, gate are described and allowed, and GT structure is associated with OM structure by aperture sequence number;

Described CB structure is used for storing specific gatage combined information, its content comprises number of regulation and the aperture combined information consisting of several gate sequence numbers, aperture sequence number binary information group, and CB structure is associated with GT structure with aperture sequence number by gate sequence number.

In step (2), operating gate rule is carried out to instantiation, refer to that the gate for actual power station allows aperture, opening sequence and working rule information, difference instantiation OM structure, GT structure and CB structure, and form aperture table, gate table and the combination table for this power station.Specifically comprise the following steps:

First, according to every kind of gate in actual power station, allow aperture, build respectively corresponding OM structure example, explain the relevant information that this gate allows aperture.The aperture table of all common completes of OM structure example.

Secondly, according to each floodgate in actual power station, build respectively corresponding GT structure example, explain the relevant information of this gate.The gate table of all common completes of GT structure example.

Finally, every kind of gatage combination according to actual power station, builds respectively corresponding CB structure example, explains the relevant information of this aperture combination.The combination table of all common completes of CB structure.

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

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

Take following methods to check between each CB structure example of described combination table, whether order is correct:

(a1) at flood, rise between water transfer position and the highest reservoir level of permission by selected several water levels of equidistant principle;

(a2) selected one of them water level, calculates flood discharge flow corresponding to different CB structure examples under this water level successively;

(a3) verify whether whether all CB structure examples meet flood discharge flow with its combination sequence number monotonically increasing rule, if satisfied by checking, if do not meet explanation, wherein must there is not the information of a CB structure example to have mistake, need the correctness of verification CB structure example order or CB structure example information, for problem, revise CB structure example order or CB structure example information until by checking;

(a4) if do not traveled through all water levels, jump to step (a2), if traveled through all water levels, finish checking.

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

(b1) according to the reservoir level at the beginning of the period, calculate the actual flood discharge flow of various combination sequence number CB structure example under this reservoir level, and record absolute value difference between actual flood discharge flow and expection period optimization flood discharge flow or difference square;

(b2) using square minimum CB structure example of absolute value difference or difference as period optimizing decision, combined the optimum combination sequence number of sequence number as this period;

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

(b4) utilize the reservoir level at the beginning of water balance equation calculates the next period;

(b5) whether travel through all periods, if do not traveled through all periods, gone to step (b1), if traveled through all periods, stopped ergodic process;

(b6) a plurality of optimum gatage combination number obtaining is retained in sequence, can forms the gatage combination sequence of optimization.

In step (4), the adjustment of gatage combination sequence optimisation comprises the adjustment that scheduling process reservoir level is out-of-limit and the adjustment of dispatching end of term reservoir level;

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

(c1) take into account the accuracy requirement of Dispatching Flood process control and the permission requirement that gate regulates frequency, determine segment length when rational gatage combination sequence optimisation is adjusted, and using the first period as the current period;

(c2) judge the whether higher limit of reservoir level of period end, if not higher limit directly enters the next stage;

(c3) if the higher limit of period Mo reservoir level is found recent operating gate forward;

(c4) judge that this gate adjustment open a sluice gate or barrier gate, if gatage combination sequence number is compared increase with the front period, illustrate that this operating gate operates for opening a sluice gate; If gatage combination sequence number was compared and reduced with the front period, illustrate that this operating gate is barrier gate operation;

(c5), if open a sluice gate operation, should will open a sluice gate constantly in advance.If segment length is Δ T ' when gatage combination sequence optimisation is adjusted, former operating gate is t constantly _r, suppose that operating gate is t constantly _r-Δ T ', and from t _rat the beginning of the residing period in-Δ T ' moment, start to carry out to calculate by period water level, if this period water level is still out-of-limit, suppose that operating gate is t constantly _r-2 Δ T ', the like until current period water level is not out-of-limit.If barrier gate operation, should postpone barrier gate constantly.Suppose that operating gate is t constantly _r+ Δ T ', and from t _rat the beginning of the residing period in+Δ T ' moment, start to carry out to calculate by period water level, if this period water level is still out-of-limit, suppose that operating gate is t constantly _r+ 2 Δ T ', the like until current period water level is not out-of-limit;

(c6) if in step (c5) operating gate constantly moved forward to an operating gate more early constantly or after move to a slower operating gate constantly, gatage combination sequence number more early or slower should be deleted, if operating gate has constantly moved forward to schedule periods just or after move to current period Mo, stop the out-of-limit adjustment of water level of current period, enter next period and go to step (c2);

(c7) if do not traveled through all periods, enter next period and go to step (c2), if traveled through all periods, exiting the out-of-limit adjustment of scheduling process water level.

The out-of-limit method of adjustment of the method for adjustment ultimate principle of described scheduling end of term reservoir level and scheduling process reservoir level is similar, and concrete steps are as follows:

(d1), according to the deviation requirement of scheduling end of term reservoir level and target water level, gatage combination sequence optimisation is adjusted to the suitable encryption of period progress row, as the time segment length Δ T ' ' that dispatches end of term reservoir level adjustment;

(d2) comprehensively the judgement out-of-limit situation of scheduling end of term reservoir level and for the last time operating gate situation, determine correct gate adjustment strategy.If scheduling end of term water level is barrier gate operation higher than target water level and last operating gate, or the water level in the scheduling end of term lower than target water level and last operating gate for opening a sluice gate operation, 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 lower than target water level and last operating gate, should shift to an earlier date the last operating gate moment;

(d3) suppose that last operating gate is t constantly _r.If should postpone last operating gate constantly, first find the period at last operating gate moment place, then, from the beginning of this period, suppose that respectively last operating gate is t constantly _r+ Δ T ' ', t _r+ 2 Δ T ' ', until the scheduling end of term, and utilize principle of water balance to carry out calculating by period water level, obtain corresponding to above-mentioned each operating gate scheduling end of term water level constantly, find wherein with the scheduling end of term water level of target water level deviation minimum and corresponding operating gate constantly, be operating gate after optimizing and revising constantly;

(d4) if should shift to an earlier date last operating gate constantly, first find the period at second from the bottom operating gate moment place, then from the beginning of this period (if altogether operating gate for once, from the beginning of first period), supposes that respectively last operating gate is t constantly _r-Δ T ' ', t _r-2 Δ T ' ', until at the beginning of this period, and utilize principle of water balance to carry out calculating by period water level, obtain corresponding to above-mentioned each operating gate scheduling end of term water level constantly, find wherein with the scheduling end of term water level of target water level deviation minimum and corresponding operating gate constantly, be operating gate after optimizing and revising constantly.

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

By period contrast, optimize the gatage combination sequence number of adjacent two periods in gatage combination sequence, and by the contrast of two gatage combined informations is resolved, draw the operating gate plan of this period.Travel through all periods and the operating gate plan of each period is all saved in a sequence, forming the operating gate plan in schedule periods.

If the gatage of adjacent two periods combination sequence number is identical, this period is not carried out operating gate; If rear period aperture combination sequence number is greater than last period aperture combination sequence number, this period is opened a sluice gate to operate; If after period aperture combination sequence number be less than last period aperture combination sequence number, this period is carried out to barrier gate operation.

The present invention has considered the working rules such as the permission aperture, opening sequence, aperture combination of gate of hydropower station in Optimum Flood Dispatching Model, can directly export the Optimum Operation plan of each gate between overflow stage, significantly reduce dispatcher's workload and work difficulty; Compare with conventional solution, the actual discharge capacity characteristic that flood discharge discharge curve meets power station is more optimized in the power station that the present invention calculates, can predict more exactly the Evolution of reservoir level and downstream flood discharge flow during Dispatching Flood, the risk of auxiliary dispatching personnel during to Dispatching Flood better studied and judged.

Accompanying drawing explanation

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

Fig. 2 is the gate flood discharge flow of stepped variation;

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

Fig. 4 is that OM structure, GT structure and CB structure all adopt text structure piece to represent;

Fig. 5 is operating gate rule example table;

Fig. 6 is the classification of correctness verification;

Fig. 7 is that aperture combination number changes example;

Fig. 8 is that operating gate plan generates.

Embodiment

For technological means, creation characteristic that the present invention is realized, reach object and effect is easy to understand, below in conjunction with embodiment, further set forth the present invention.

Referring to Fig. 1, Flood Dispatching Optimization method of the present invention comprises that the description of operating gate rule structuralisation, the instantiation of operating gate rule, the generation of gatage combination sequence optimisation, the adjustment of gatage combination sequence optimisation, operating gate plan generate five steps.

Wherein, operating gate rule structuralisation is described and have been adopted structurized describing method to come complete expression gate to allow the complex informations such as aperture, opening sequence and aperture combination, for operating gate rule instantiation and computer analyzing operating gate rule provide basis.

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

Gatage combination sequence optimisation generates the theoretical flood discharge discharge curve of inquiring into according to traditional Flood Dispatching Optimization method, by the period, calculate optimum aperture and combine sequence number, form the optimization gatage combination sequence in whole schedule periods, for solving the coupled problem of operating gate rule and traditional Flood Dispatching Optimization method.

The adjustment of gatage combination sequence optimisation, for the out-of-limit period of reservoir level in Dispatching Flood process, is adjusted the last time operating gate moment by designation method, until eliminate the out-of-limit problem of reservoir level; For the deviation situation of scheduling end of term reservoir level and target water level and operating gate situation last time, by the last operating gate moment of designation method adjustment, until it is the most approaching to dispatch end of term water level.

Operating gate plan generates by optimizing resolving by the period of gatage combination sequence, and dynamically generation meets the operating gate plan of all operations rule request.

As can be seen here, only gatage combination sequence optimisation generating portion belongs to the improvement to traditional Flood Dispatching Optimization method, and other four parts propose first by the present invention, for further supplementing and improve traditional Flood Dispatching Optimization method system.

Operating gate rule structuralisation is described as all the other four parts provides basic standard to support, and all the other four parts form the complete step of power station Flood Dispatching Optimization method.

(1) operating gate rule structuralisation is described

Gate of hydropower station working rule mainly comprises permission aperture, opening sequence and the aperture combination of gate.Allow aperture to refer to the permission aperture position of each gate, can adopt height (definitely aperture) or the ratio (relative opening degree) of actual Lift and maximum Lift of the actual unlatching of gate to represent.In Dispatching Flood process, the actual gatage of choosing is necessary for one of several discrete apertures of permission.Therefore, gatage can only strengthen step by step, and corresponding flood discharge flow is step-like variation, cannot accurately follow the tracks of expection flood discharge discharge process.It is the principal element (see figure 2) that traditional Flood Dispatching Optimization method is carried out link flood discharge flow deviation that gate allows aperture.Owing to allowing aperture all for specific gate, so its structuring expression is fairly simple, only need to be by the set of aperture type (definitely aperture, relative opening degree) and permission aperture, as { 0,0.25,0.5,0.75,1}, is saved in corresponding gate information.

Part gate of hydropower station scheduling rule has simply been arranged the priority opening sequence between a plurality of gates, and could open the gate that priority is taken second place after the highest that gate standard-sized sheet of regulation priority.When the flood discharge of needs starting gate, that the highest gate of first operator precedence level, along with flood into reservoir flow increases gradually and strengthens the aperture of this gate.When this gate reaches fully open position and flood into reservoir continuation increase, choose one or more gates that priority is taken second place, and in time increase or reduce gatage.In the flood water-break stage, adopt contrary adjusting order, the gate of opening afterwards first closes, and after the gate of first opening, closes.For this situation, only need to adopt " priority sequence " can meet structured storage requirement.

Yet, also have part gate of hydropower station scheduling rule to require to regulate in turn between a plurality of gates, to optimize the HYDRODYNAMIC CONDITION RELATING TO in downstream, power station.For example, some hydropower station has A(left side), in the middle of B(), C(right side) three gates, requirement is along with the increasing gradually of flood into reservoir flow, first open B gate to 0.25 aperture and 0.5 aperture, then keep B gate in 0.5 aperture, open A, two gate to 0.25 apertures of C and 0.5 aperture, simultaneously when flood discharge flow further increases, again strengthen B gate to 0.75 aperture and standard-sized sheet, and then strengthen successively A, C gate to 0.75 aperture and standard-sized sheet.In this case, because low priority gate can just be opened in advance before high priority gate standard-sized sheet, above-mentioned " priority sequence " method cannot meet the gate opening sequential organization memory requirement in this situation.

For addressing the above problem, the present invention stores different gatage combined situation, and by the ordered arrangement of different gatage combined informations, implies the complicated gate opening order information of storage, realizes the indirect statement of complicated gate opening order.The method not only can solve the complicated opening sequence statement problem of above-mentioned regulating gate in turn well, and can effectively compatible and alternative simple " priority sequence " expression method.In addition, the method has also realized the integrated storage of different gate opening orders and aperture combined information, has reduced the redundance of information storage.

According to above-mentioned analysis, referring to Fig. 3, adopt respectively OM, GT, tri-kinds of structures of CB to come permission aperture, opening sequence and the aperture combined information of characterized gate.OM structure is for storing the information of the specific aperture of gate, and content comprises aperture sequence number, opening value, aperture type, aperture description, and the overflow performance curve consisting of several water levels, flow binary information group.GT structure is for the information of the specific gate in storage of water power station, and content comprises the sequence number that gate sequence number, gate are described and allowed aperture, associated with OM structure by aperture sequence number.CB structure is used for storing specific gatage combined information, and content comprises number of regulation, and the aperture combined information consisting of several gate sequence numbers, aperture sequence number binary information group, associated with GT structure with aperture sequence number by gate sequence number.This by the mode of opening information and gate unpack storage, can avoid pattern and Installation Elevation all identical a plurality of gates re-define identical opening information, redundant information is minimized.

(2) operating gate rule instantiation

Operating gate rule structuralisation is described the statement problem that only solves different types of information.For specific power station, need to different structure body, carry out instantiation according to permission aperture and the aperture combined situation of its concrete gate quantity, each gate, to utilize structuring expression method to describe all information in specific power station.

Because gate quantity, permission aperture, the aperture combination in different power stations are all different, therefore the quantity of all kinds of structure instantiations is also different.The quantity of OM structure example depends on the quantity of gate of hydropower station and allows the quantity of aperture.In principle, should all set up an OM structure example to all permission apertures of each gate.But when there are a plurality of gates in power station, and the overflow performance curve homogeneous phase under its permission aperture and corresponding aperture simultaneously, can not re-define the OM structure example of identical aperture.Take above-mentioned A, B, tri-gates of C is example, if three gate patterns, Installation Elevation are all identical, allow aperture to be { 0.25,0.5,0.75,1 }, and the overflow performance curve under arbitrary aperture is also identical, and actual need to define 4 OM structure examples, opening value is respectively 0.25,0.5,0.75,1.The quantity of GT structure example equates with the quantity of gate of hydropower station, carries out instantiation respectively for each gate.The quantity of CB structure example equates with the quantity of each gatage combination of power station, carries out instantiation respectively for each aperture combination.

In order to eliminate the dependence to commercial data base management system (DBMS), realize light-weighted information configuration, the present invention adopts the block structure of textual form to represent operating gate Rule Information, and computer program can be resolved these block structures, and the corresponding structure of instantiation.It should be noted that the rating curve in OM example should be in strict accordance with the arrangement of water level progressive law, so that calculate flood discharge flow corresponding to this aperture according to reservoir level Fast Interpolation in solution procedure.Still take above-mentioned A, B, tri-gates of C is example, and its operating gate rule statement as shown in Figure 4.

Adopt the mode of object array that the structure example of three types is all written in calculator memory, set up respectively aperture table, gate table and combination table as shown in Figure 5, each member in each table is a structure example.Aperture table and gate table, to the order no requirement (NR) between each example, still in order more clearly to reflect aperture Changing Pattern and gate situation, should be deposited according to certain rule in order.For example, the structure example in aperture table can be arranged from small to large according to opening value.Structure example in gate table can be arranged according to the type of gate and physical location.Each example in combination table should carry out ordered arrangement in strict accordance with gate opening sequencing and combined situation, to follow the rational operating gate plan of operating gate Rulemaking.

When given reservoir level z, the flood discharge flow that CB example is corresponding can uniquely be determined.Method is all aperture sequence numbers in traversal CB example, find respectively OM example corresponding to this aperture sequence number, then travel through water level, flow information group in this OM example, find respectively two adjacent water level values that are less than or equal to and are more than or equal to z, be designated as z ' and z ' ', corresponding flood discharge flow is respectively q ' and q ' ', and when reservoir level is z, the flow of this OM example can be expressed as q _om=q '+(q ' '-q ') (z-z ')/(z ' '-z ').If do not find z ' and z ' ', when reservoir level is z, the flood discharge flow of this OM example is zero.The flood discharge flow of each OM example in cumulative CB example, obtains total flood discharge flow successively, is the flood discharge flow of this CB example.Therefore, flood discharge flow is the function of reservoir level and aperture combination sequence number, is designated as q _cb=f _rq(z, c), wherein z is reservoir level, c is aperture combination sequence number, q _cbfor the total flood discharge flow in power station under this gatage combination.

All structure examples are written into after internal memory, and reply example information is carried out correctness verification, is divided into three kinds of sequence number verification, associated verification and data verifications, sees Fig. 6.Wherein, sequence number verification is carried out verification to the combination sequence number of the gate sequence number of the aperture sequence number of each OM example, GT example and CB example.Require each sequence number to increase progressively with sequence of natural numbers, must not duplicate or jump.Associated verification is for verifying that whether the relevance between each example is correct, comprise whether the aperture sequence number in gate table exists corresponding OM example in aperture table, whether the gate sequence number in combination table there is corresponding GT example in gate table, and whether corresponding aperture sequence number is present in this GT example.Whether data verification there is mistake for verifying data data self, and this mistake may come from source book, also may come from the artificial error of Data Input Process.Verification content mainly comprises that whether the rating curve of OM example is reasonable, and each CB example of combination table and sequentially whether reasonable.The principle verification OM example larger according to the higher flood discharge flow of reservoir level under same aperture, contrasts each water level, the flow binary information group of OM example one by one, checks whether water level and flow are monotone increasing, do not meet the demands and judge that this OM example exists mistake.Combination table has been stored along with flood into reservoir increases and strengthens step by step the aperture combined information of starting gate, so under same bank water level, and the letdown flow that closes sequence number CB example compared with group should be less than the letdown flow of larger combination sequence number CB example.Based on this principle, can take following steps to check each CB example of combination table and sequentially whether reasonable:

(a1) at flood, rise between water transfer position and the highest reservoir level of permission by selected several water levels of equidistant principle;

(a2) selected one of them water level, calculates flood discharge flow corresponding to different CB examples under this water level successively;

(a3) verify whether all CB structure examples meet flood discharge flow with combination sequence number monotonically increasing rule, if satisfied by checking, if do not meet explanation, wherein must there is not the information of a CB example to have mistake, need the correctness of verification CB structure example order or CB structure example information, for problem, revise CB structure example order or CB structure example information until by checking;

(a4) if do not traveled through all water levels, jump to step (a2), if traveled through all water levels, finish checking.

(3) gatage combination sequence optimisation generates

Gate of hydropower station need to be followed certain opening sequence and aperture combination requirement, and this two category information is all stored in combination table.Single CB example has represented a kind of specific gatage combination, and ordered arrangement between CB example can guarantee that the operating gate scheme that Optimization Solution obtains meets opening sequence requirement all the time.The in the situation that of correct structure aperture table, gate table and combination table, only need in the combination table of ordered arrangement, search for the optimum CB sequence number of different periods, can consider permission aperture, opening sequence and the aperture combination constraint of gate comprehensively, efficiently solve the coupled problem that operating gate rule and Flood Dispatching Optimization calculate.

In addition, power station allows to cross over and regulate between different gatage combinations, for example, when flood increases very soon, gate in complete shut-down position can be crossed over 0.25 aperture and directly open to 0.5 aperture, so the redirect between any two combination sequence numbers, all can meet gate opening order and aperture combination requirement naturally.Therefore, between the combination sequence number of adjacent time interval, do not have interconnection constraint, the optimizing space of gatage combination sequence number is whole combination table, makes the optimizing operation of aperture combination sequence number obtain great simplification.

As previously mentioned, in traditional Flood Dispatching Optimization method, power station discharge capacity is constrained to q _i≤ q _max(z _i), in this process, do not consider the working rule of gate, make model fail truly to reflect the operating gate rule constrain in power station.Structuring by operating gate rule is explained, and flood discharge flow in power station can be expressed as the function of reservoir level and aperture combination sequence number, i.e. q _cb=f _rq(z, c).Because each CB example in combination table is according to gatage situation ordered arrangement, therefore can utilize maximum combination sequence number to try to achieve the power station discharge capacity under particular bin water level.So the power station discharge capacity constraint in traditional Flood Dispatching Optimization method can be revised as q _i≤ f _rq(z _i, c _max), c wherein _maxmaximum opening combination sequence number for power station.Utilize this amended mathematical model and traditional method for solving, try to achieve the optimization flood discharge discharge curve of expection.

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

Objective function:

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

Constraint condition:

Actual letdown flow constraint $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 when nothing is let out in advance $q_i^j\&le;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)\&CenterDot;\mathrm{\&Delta;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)$

The constraint of reservoir peak level $z_{i+1}^j\&le;f_{\mathrm{vz}}({\mathrm{<figure-callout\; id="0"\; label="V"\; filenames="HDA0000412922920000021.png,HDA0000412922920000031.png"\; state="\{\{state\}\}">V</figure-callout>}}_0+V_f),j=\mathrm{1,2},L,R;i=\mathrm{1,2},L,N$

In formula, q _ibe the optimization flood discharge flow of i period,

be the flood discharge flow of j aperture of i choosing period of time combination number, be j gatage combination number, R is gatage combination quantity, f _rqthe funtcional relationship between flood discharge flow and reservoir level, gatage combination number, f _vzfor the funtcional relationship of reservoir level and storage capacity,

be period Mo (at the beginning of the i+1 period) storage capacity in j aperture combination situation of i choosing period of time,

it is period Mo (at the beginning of the i+1 period) reservoir level in j aperture combination situation of i choosing period of time.If objective function is minimum during j=0, this period should be by all gate Close Alls.Because actual letdown flow constraint has reflected the restriction of gatage to letdown flow, so this model does not need to consider maximum discharge capacity constraint.

In period optimizing process, there will be the situation that the constraint of traffic constraints and peak level cannot be simultaneously satisfied.If keep outbound flow to be less than reservoir inflow, the reservoir level of period end can surpass high limit.If maintenance period Mo water level, below high limit, must guarantee that outbound flow is greater than reservoir inflow.Now, should preferentially guarantee that outbound flow is less than reservoir inflow, avoid causing artificial flood, and solve the out-of-limit problem of water level by follow-up operating gate planning optimization method of adjustment.

Adopt the concrete solution procedure of above-mentioned mathematical model solving-optimizing gatage combination sequence as follows:

(b1) according to the reservoir level at the beginning of the period, calculate the actual flood discharge flow of various combination sequence number CB structure example under this reservoir level, and record absolute value difference between actual flood discharge flow and expection period optimization flood discharge flow or difference square;

(b2) using square minimum CB structure example of absolute value difference or difference as period optimizing decision, combined the optimum combination sequence number of sequence number as this period;

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

(b4) utilize the reservoir level at the beginning of water balance equation calculates the next period;

(b5) whether travel through all periods, if do not traveled through all periods, gone to step (b1), if traveled through all periods, stopped ergodic process.

(b6) N obtaining an optimum gatage combination number is retained in sequence, can forms the gatage combination sequence of optimization.

(4) gatage combination sequence optimisation adjustment

As previously mentioned, optimize gate flood discharge flow that gatage combination sequence only pursues each stage close to the optimization flood discharge flow of expecting, cannot avoid occurring in Dispatching Flood process the problem that reservoir level transfinites.In addition, the impact of later stage water situation on reservoir risk considered in part power station, and the demand of flood resource optimization utilization, can require Dispatching Flood end of term reservoir level to return back to designated water level.The adjustment of gatage combination sequence optimisation is for solve the out-of-limit problem of above-mentioned scheduling process reservoir level and scheduling end of term water level adjustment problem simultaneously.

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

(c1) take into account the accuracy requirement of Dispatching Flood process control and the permission requirement that gate regulates frequency, determine segment length's Δ T ' when rational gatage combination sequence optimisation is adjusted, Δ T ' should get 15 minutes or 10 minutes, using the first period as the current period;

(c2) judge the whether higher limit of reservoir level of period end, if not higher limit directly enters the next stage;

(c3) if the higher limit of period Mo reservoir level is found recent operating gate (period that gatage combination sequence number changes) forward;

(c4) judge that this gate adjustment open a sluice gate or barrier gate.If gatage combination sequence number is compared increase with the front period, illustrate that this operating gate is for opening a sluice gate operation; If gatage combination sequence number was compared and reduced with the front period, illustrate that this gate is barrier gate operation;

(c5), if open a sluice gate operation, should will open a sluice gate constantly in advance.If former operating gate is t constantly _r, suppose that operating gate is t constantly _r-Δ T ', and from t _rat the beginning of the residing period in-Δ T ' moment, start to carry out to calculate by period water level, if this period water level is still out-of-limit, suppose that operating gate is t constantly _r-2 Δ T ', the like until current period water level is not out-of-limit.If barrier gate operation, should postpone barrier gate constantly.Suppose that operating gate is t constantly _r+ Δ T ', and from t _rat the beginning of the residing period in+Δ T ' moment, start to carry out to calculate by period water level, if this period water level is still out-of-limit, suppose that operating gate is t constantly _r+ 2 Δ T ', the like until current period water level is not out-of-limit.

(c6) if in step (c5) operating gate constantly moved forward to an operating gate more early constantly or after move to a slower operating gate constantly, gatage combination sequence number more early or slower should be deleted, if operating gate has constantly moved forward to schedule periods just or after move to current period Mo, stop the out-of-limit adjustment of water level of current period, enter next period and go to step (c2);

(c7), if do not traveled through all periods, enter next period and go to step (c2); If traveled through all periods, 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	Operating gate situation last time	Gate is adjusted strategy
			1	Barrier gate	Postpone the barrier gate time
2	Open a sluice gate	Open a sluice gate the time in advance

Referring to table 2, the out-of-limit method of adjustment of the scheduling method of adjustment ultimate principle of end of term reservoir level and scheduling process reservoir level is similar, and concrete steps are as follows:

(d1), according to the deviation requirement of scheduling end of term reservoir level and target water level, gatage combination sequence optimisation is adjusted to the suitable encryption of period progress row, as the time segment length Δ T ' ' that dispatches end of term reservoir level adjustment;

(d2) comprehensively the judgement out-of-limit situation of scheduling end of term reservoir level and for the last time operating gate situation, determine correct gate adjustment strategy.If scheduling end of term water level is barrier gate operation higher than target water level and last operating gate, or the water level in the scheduling end of term lower than target water level and last operating gate for opening a sluice gate operation, 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 lower than target water level and last operating gate, should shift to an earlier date the last operating gate moment;

(d3) suppose that last operating gate is t constantly _r.If should postpone last operating gate constantly, first find the period at last operating gate moment place, then, from the beginning of this period, suppose that respectively last operating gate is t constantly _r+ Δ T ', t _r+ 2 Δ T ', until the scheduling end of term, and utilize principle of water balance to carry out calculating by period water level, obtain corresponding to above-mentioned each operating gate scheduling end of term water level constantly, find wherein with the scheduling end of term water level of target water level deviation and corresponding operating gate constantly, be operating gate after optimizing and revising constantly;

(d4) if should shift to an earlier date last operating gate constantly, first find the period at second from the bottom operating gate moment place, then from the beginning of this period (if altogether operating gate for once, from the beginning of first period), supposes that respectively last operating gate is t constantly _r-Δ T ', t _r-2 Δ T ', until at the beginning of this period, and utilize principle of water balance to carry out calculating by period water level, obtain corresponding to above-mentioned each operating gate scheduling end of term water level constantly, find wherein with the scheduling end of term water level of target water level deviation and corresponding operating gate constantly, be operating gate after optimizing and revising constantly.

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 having obtained.Basic thought is by the period, to contrast the aperture combination number of adjacent two periods, by the contrast of two gatage combined informations is resolved, draws the operating gate plan of this period.Travel through all periods and the operating gate plan of each period is all saved in a sequence, forming the operating gate plan that meets all operating gate rule requests.

If the gatage of adjacent two periods combination sequence number is identical, this period is not carried out operating gate.If the gatage of adjacent two periods combination sequence number is different, extracts corresponding 2 operating gates rule example, and by the comparison of content, obtain the operation scheme of period gate.Because each CB example in combination table is arranged from small to large according to gate opening situation, therefore can be simply the size of gatage combination by adjacent two periods number, judgement is to open a sluice gate this period operate or carry out barrier gate operation.If rear period aperture combination sequence number is greater than last period aperture combination sequence number, this period is opened a sluice gate operation; Otherwise this period is carried out barrier gate operation.

For convenience of explanation, still with the aperture shown in Fig. 4, be combined as example.Consider respectively to open a sluice gate and two kinds of situations of barrier gate, the situation of opening a sluice gate supposes that former and later two are that the optimization aperture combination number of period is to be respectively 1 and 3, barrier gate situation supposes that former and later two are that the optimization aperture combination number of period is to be respectively 3 and 1, and extraction gate aperture combined information is as follows respectively:

The gatage combination of Fig. 7 left part number becomes 3 from 1, therefore for opening a sluice gate, operates.By drawing the compare of analysis of two CB example content, gate 2(gate B) need to be from aperture 1(0.25 aperture) be adjusted to aperture 2(0.5 aperture), gate 1(gate A) and gate 3(gate C) need to be opened into aperture 1(0.25 aperture from complete shut-down state).

In like manner, the gatage of Fig. 7 right side part combination number becomes 1 from 3, therefore for barrier gate operates.By obtaining the compare of analysis of two CB example content, gate 2(gate B) need to be from aperture 2(0.5 aperture) be adjusted to aperture 1(0.25 aperture), gate 1(gate A) and gate 3(gate C) need to be from aperture 1(0.25 aperture) be adjusted to complete shut-down position.See Fig. 8.

The method generating from gatage combination sequence optimisation, the adjusting time of above-mentioned gate all should be the moment value at the beginning of this period.The information that gate regulates can be described by gate title, the adjusting moment, target aperture three metamessage groups, when certain period has a plurality of gates to regulate, can consist of a plurality of three metamessage groups the operating gate plan of this period.Recycling said method travels through all periods, can obtain the operating gate plan of each period, jointly forms the operating gate plan in whole schedule periods.

More than show and described ultimate principle of the present invention and principal character and advantage of the present invention.The technician of the industry should understand; the present invention is not restricted to the described embodiments; that in above-described embodiment and instructions, describes just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the claimed scope of the invention.The claimed scope of the present invention is defined by appending claims and equivalent thereof.

Claims (9)

1. a power station Flood Dispatching Optimization method of considering operating gate rule, is characterized in that, comprises successively following step:

(1) operating gate rule structuralisation is described, and adopts structurized describing method to come complete expression gate to allow aperture, opening sequence and aperture combined information, for operating gate rule instantiation and computer analyzing operating gate rule provide basis;

(2) for specific power station, carry out the instantiation of operating gate rule, according to the true gate in specific power station, allow aperture, opening sequence and aperture combined information, dynamically generate and arrange the Multi-instance of all kinds of operating gate regular texture bodies, and set up the incidence relation between different instances;

(3) generate the gatage combination sequence of optimizing, the theoretical flood discharge discharge curve of inquiring into according to traditional power station Flood Dispatching Optimization method, by the period, calculate optimum aperture and combine sequence number, form the optimization gatage combination sequence in whole schedule periods, for solving the coupled problem of operating gate rule and traditional Flood Dispatching Optimization method;

(4) gatage combination sequence is optimized to adjustment, for the out-of-limit period of reservoir level in Dispatching Flood process, adjusts the last time operating gate moment, until eliminate the out-of-limit problem of reservoir level; For the deviation situation of scheduling end of term reservoir level and target water level and operating gate situation last time, adjust last operating gate constantly, until it is the most approaching to dispatch end of term water level;

(5) generate operating gate plan, by optimizing resolving by the period of gatage combination sequence, dynamically generation meets the operating gate plan of all operating gate rule requests.

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

In step (1), operating gate rule structuralisation describing method is as follows:

Adopt respectively OM structure, GT structure and CB structure to characterize wherein the permission aperture of gate, opening sequence and aperture combined information;

OM structure is for storing the information of the specific aperture of gate, and content comprises the rating curve that aperture sequence number, opening value, aperture type, aperture are described and consisted of several water levels, flow binary information group;

GT structure is for the information of the specific gate in storage of water power station, and content comprises the aperture sequence number that gate sequence number, gate are described and allowed, and GT structure is associated with OM structure by aperture sequence number;

CB structure is used for storing specific gatage combined information, and content comprises number of regulation and the aperture combined information consisting of several gate sequence numbers, aperture sequence number binary information group, and CB structure is associated with GT structure with aperture sequence number by gate sequence number.

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

In step (2), operating gate rule is carried out to instantiation, refer to that the gate for actual power station allows aperture, opening sequence and working rule information, difference instantiation OM structure, GT structure and CB structure, and form aperture table, gate table and the combination table for actual power station, specifically comprise the following steps:

First, according to every kind of gate in actual power station, allow aperture, build respectively corresponding OM structure example, explain this gate and allow opening information, the aperture table of all common completes of OM structure example;

Secondly, according to each floodgate in actual power station, build respectively corresponding GT structure example, explain this gate information, the gate table of all common completes of GT structure example;

Finally, every kind of gatage combination according to actual power station, builds respectively corresponding CB structure example, explains this aperture combined information, the combination table of all common completes of CB structure.

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

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

Each OM structure example in described aperture table can be arranged from small to large according to opening value; Each GT structure example in described gate table can be arranged according to the type of gate and physical location; Each CB structure example in described combination table should be arranged in strict accordance with gate opening sequencing and aperture combination order from small to large.

5. the power station Flood Dispatching Optimization method of consideration operating gate rule according to claim 4, is characterized in that, takes following methods to check between each CB structure example of described combination table, whether order is correct:

(a1) at flood, rise between water transfer position and the highest reservoir level of permission by selected several water levels of equidistant principle;

(a2) selected one of them water level, calculates flood discharge flow corresponding to different CB structure examples under this water level successively;

(a3) verify whether all CB structure examples meet flood discharge flow with combination sequence number monotonically increasing rule, if satisfied by checking, if do not meet explanation, wherein must there is not the information of a CB structure example to have mistake, need the correctness of verification CB structure example order or CB structure example information, for problem, revise CB structure example order or CB structure example information until by checking;

(a4) if do not traveled through all water levels, jump to step (a2), if traveled through all water levels, finish checking.

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

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

(b1) according to the reservoir level at the beginning of the period, calculate the actual flood discharge flow of various combination sequence number CB structure example under this reservoir level, and record absolute value difference between actual flood discharge flow and expection period optimization flood discharge flow or difference square;

(b2) using square minimum CB structure example of absolute value difference or difference as period optimizing decision, combined the optimum combination sequence number of sequence number as this period;

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

(b4) utilize the reservoir level at the beginning of water balance equation calculates the next period;

(b5) whether travel through all periods, if do not traveled through all periods, gone to step (b1), if traveled through all periods, stopped ergodic process;

(b6) a plurality of optimum gatage combination number obtaining is retained in sequence, can forms the gatage combination sequence of optimization.

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

In step (4), the adjustment of gatage combination sequence optimisation comprises the adjustment that scheduling process reservoir level is out-of-limit and the adjustment of dispatching end of term reservoir level;

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

(c1) take into account the accuracy requirement of Dispatching Flood process control and the permission requirement that gate regulates frequency, determine segment length when rational gatage combination sequence optimisation is adjusted, and using the first period as the current period;

(c2) judge the whether higher limit of reservoir level of period end, if not higher limit directly enters the next stage;

(c3) if the higher limit of period Mo reservoir level is found recent operating gate forward;

(c4) judge that this operating gate open a sluice gate or barrier gate, if gatage combination sequence number is compared increase with the front period, illustrate that this operating gate operates for opening a sluice gate; If gatage combination sequence number was compared and reduced with the front period, illustrate that this operating gate is barrier gate operation;

(c5) if open a sluice gate operation, should will open a sluice gate constantly in advance, if segment length is Δ T ' when operating gate planning optimization is adjusted, former operating gate is t constantly _r, suppose that operating gate is t constantly _r-Δ T ', and from t _rat the beginning of the residing period in-Δ T ' moment, start to carry out to calculate by period water level, if this period water level is still out-of-limit, suppose that operating gate is t constantly _r-2 Δ T ', the like until current period water level is not out-of-limit, if barrier gate operation should postpone barrier gate constantly, suppose that operating gate is t constantly _r+ Δ T ', and from t _rat the beginning of the residing period in+Δ T ' moment, start to carry out to calculate by period water level, if this period water level is still out-of-limit, suppose that operating gate is t constantly _r+ 2 Δ T ', the like until current period water level is not out-of-limit;

(c6) if in step (c5) operating gate constantly moved forward to an operating gate more early constantly or after move to a slower operating gate constantly, gatage combination sequence number more early or slower should be deleted, if operating gate has constantly moved forward to schedule periods just or after move to current period Mo, stop the out-of-limit adjustment of water level of current period, enter next period and go to step (c2);

(c7) if do not traveled through all periods, enter next period and go to step (c2), if traveled through all periods, exiting 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 and target water level, gatage combination sequence optimisation is adjusted to long being encrypted of period, as the time segment length Δ T ' ' that dispatches end of term reservoir level adjustment;

(d2) comprehensively the judgement out-of-limit situation of scheduling end of term reservoir level and for the last time operating gate situation, determine correct gate adjustment strategy; If scheduling end of term water level is barrier gate operation higher than target water level and last operating gate, or the water level in the scheduling end of term lower than target water level and last operating gate for opening a sluice gate operation, 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 lower than target water level and last operating gate, should shift to an earlier date the last operating gate moment;

(d3) suppose that last operating gate is t constantly _rif, should postpone last operating gate constantly, first find the period at last operating gate moment place, then, from the beginning of this period, suppose that respectively last operating gate is t constantly _r+ Δ T ' ', t _r+ 2 Δ T ' ', until the scheduling end of term, and utilize principle of water balance to carry out calculating by period water level, obtain corresponding to above-mentioned each operating gate scheduling end of term water level constantly, find wherein with the scheduling end of term water level of target water level deviation minimum and corresponding operating gate constantly, be operating gate after optimizing and revising constantly;

(d4) if should shift to an earlier date last operating gate constantly, first find second from the bottom operating gate period at place constantly, then from the beginning of this period, if operating gate for once altogether,, from the beginning of first period, suppose that respectively last operating gate is t constantly _r-Δ T ' ', t _r-2 Δ T ' ', until at the beginning of this period, and utilize principle of water balance to carry out calculating by period water level, obtain corresponding to above-mentioned each operating gate scheduling end of term water level constantly, find wherein with the scheduling end of term water level of target water level deviation minimum and corresponding operating gate constantly, be operating gate after optimizing and revising constantly.

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

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

By period contrast, optimize the gatage combination sequence number of adjacent two periods in gatage combination sequence, and by the contrast of two gatage combined informations is resolved, draw the operating gate plan of this period, travel through all periods and the operating gate plan of each period is all saved in a sequence, forming the operating gate plan in schedule periods.

9. the power station Flood Dispatching Optimization method of consideration operating gate rule according to claim 8, is characterized in that,

If the gatage of adjacent two periods combination sequence number is identical, this period is not carried out operating gate; If rear period aperture combination sequence number is greater than last period aperture combination sequence number, this period is opened a sluice gate to operate; If after period aperture combination sequence number be less than last period aperture combination sequence number, this period is carried out to barrier gate operation.

Images