CN105914795A - Hydropower station automatic generation active output control method - Google Patents

Abstract

The invention discloses a hydropower station automatic generation active output control system and a hydropower station automatic generation active output control method. The hydropower station automatic generation active output control system comprises an AGC functional block, a PID functional block and a preprocessing functional block, wherein equipment parameters and dynamic parameters are directly input to, or sent to the AGC functional block and the PID functional block for updating after being processed by the preprocessing functional block. The hydropower station automatic generation active output control method can calculate a single-unit active setting value of each generator unit according to a whole station active target value of a hydropower station, performs PID closed-loop adjustment on a single-unit active actually generated value, and meets the operation requirements of a generator unit avoiding a vibration region and hydropower station primary frequency modulation; the hydropower station automatic generation active output control method adopts the method of pre-establishing a combined output model of input AGC units, converts data existing in a unit form into data in a station form, increases solution efficiency, and can meet requirements for automatic active output control of different multi-vibration-region and multi-unit type large and medium-sized hydropower stations.

Description

A kind of power station automatic generation is gained merit force control method

Technical field

The invention belongs to hydroelectric generation and control technical field, relate to a kind of power station automatic generation and gain merit force control method.

Background technology

Power station first full factory gained merit desired value commonly used for the meritorious control exerted oneself controls function by automatic generation amount Distribution was to each turbine-generator units (abbreviation unit), so (Automatic Generation Control is called for short AGC) After by each unit respectively by computer supervisory control system programmable logic controller (PLC) (Programmable Logic Controller, is called for short PLC) or speed regulator carry out the mode of PID (be called for short PID) closed loop regulation.But mesh Not yet have till before one contain and guide Hydropower Plant to go out force modeling, there is the distribution of work at complete station, unit is meritorious regulates and one The systemic algorithm of all real power controls link such as secondary frequency modulation.Automatic generation real power control merit in the actual application in each power station Can, even generally have that algorithm is coarse, integrality lacks, each functional block connecting is poor, logic does not mate conflict between strategy Etc. shortcoming, cause power station in the face of complex operating mode (such as: a frequency modulation frequency modulation conflict, pass through vibrating area etc.) time easy Various abnormal conditions occur.For solving the various problems that occur in producing, or the assistant service appraisal standards of response electrical network and Its particular/special requirement, owing to not having unified algorithm modifying automatic generation real power control function as guidance, power station Time, it not using the entirety of function as optimization object, but ignore the global configuration between strategy and parameter, simple and crudely Partial strategy and parameter are carried out isolated amendment, while solution problem, then brings new problem, add water on the contrary The unordered degree of power station Automatic Generation Control function.

Document 1 " a kind of hydropower station unit AGC control method " (application publication number CN 102664430 A) discloses A kind of AGC control method of hydropower station unit.In literary composition propose be determined by power station gain merit adjustable range, distribution send out Group of motors meritorious exert oneself, allows unit enter vibrating area in turn in the case of inevitable runs and ensures Hydropower Plant Safety and stability.But document 1 has a problem in that: only relate to vibrating area evasion, the vibration by single unit is not proposed District draws specifically building of multiple stage unit associating vibrating area or Operational Zone (in addition to associating vibrating area other gain merit exert oneself region) Mould or calculation；The proposition unit operation time in vibrating area must not exceed 10 minutes, but not to " 10 minutes " Vibrating area run duration carry out any demonstration；Propose to allow unit wheel flow to according to operation cumulative time length and unit priority Enter vibrating area run, but do not propose any relate to the time statistics, priority calculate, power distribution concrete way.

Document 2 " a kind of Automatic Generation Control power station prevents the method that head signal changes the power swing caused " (application Publication No. CN 104638681 A) disclose a kind of Automatic Generation Control power station and prevent from head signal from changing causing power waves Dynamic method, proposes in literary composition, when the change of head signal value causes the meritorious upper limit of exerting oneself of full factory to change, to perform distribution by refusal Value or exit the mode of AGC and prevent that head is abnormal causes full factory power swing.But document 2 has a problem in that: only relate to Relevant issues in terms of head, are not directed to other link of automatic generation real power control；Only propose to prevent from gaining merit during head exception The strategy of fluctuation, does not advises to the reliability how improving power station head measurement function.

Document 3 " Hydropower Stations in Large Scale group's Short-term Optimal Operation method II: high water head many vibrating areas problem " (Journal of Hydraulic Engineering 2011 Year volume 42 the 10th phase page 1168) disclose the associating vibrating area of a kind of high water head, Large Copacity, many vibrating areas unit Modeling pattern, has inquired into the combination of unit, vibrating area combination, vibrating area avoids strategy and solve flow process, to Gao Shui The short-term electricity generation optimization problem of the Huge Power Station group of vibrating area more than has directive significance.But the focal point of document 3 is The cascade operation of GROUP OF HYDROPOWER STATIONS, is not directed to single power station unit in automatic generation real power control and evades vibrating area or pass through and shake The relevant issues in dynamic district.

(power system is automatic for document 4 " project analysis of hydroelectric units primary frequency modulation and AGC typical control strategy and optimization " Change the 3rd phase of volume 39 in 2015 page 146) disclose the matching problem that primary frequency modulation controls with AGC, inquire into emphatically Under the speed regulator power control mode and aperture control model, the Different Strategies that primary frequency modulation coordinates with AGC.But it is civilian Offer 4 and unresolved under the aperture control model of current main flow, primary frequency modulation and the matching problem of AGC.

Above-mentioned document, the most completely provides the systematic approach of power station automatic generation real power control.How to ensure unit Optimal scheme target is exerted oneself under prescribed conditions, makes power station realize the maximization of economic benefit, safety benefit, is still one There is problem to be solved.

Summary of the invention

Present invention solves the problem in that a kind of power station of offer automatic generation is gained merit force control method, can be according to power station The unit that the meritorious desired value in full station calculates each unit is gained merit setting value, and each unit is carried out closed loop feedback regulation respectively, Meet unit simultaneously and evade the service requirement of vibrating area and primary frequency regulation of power network.

The present invention is to be achieved through the following technical solutions:

A kind of power station automatic generation is gained merit Force control system, including AGC functional block, pid function block and pretreatment merit Energy block, device parameter, dynamic parameter directly input or are sent to AGC functional block, PID after the process of preprocessed functional block Functional block updates；

AGC functional block, pid function block receive the pretreatment values that preprocessing function block sends respectively；

AGC functional block is set up the combination of hydropower station unit according to received data and is exerted oneself model, and calculates associating fortune Row region, associating optimized operation region and combine restriction operation area；Gain merit apportioning cost at associating optimum according to full station AGC Operation area or combining restriction operation area, starts corresponding AGC allocation flow and obtains each and put into AGC unit Unit AGC gains merit apportioning cost be issued to unit；

Put into each generating set of AGC further according to unit AGC gain merit apportioning cost and unit PID function throw state of moving back and AGC throws and moves back state, determines that unit is gained merit setting value；Each the meritorious unit that is superimposed to of generating set unit frequency modulation correction is gained merit Setting value, show that unit is meritorious and sets correction value and be sent to pid function block；

Pid function block carries out PID closed loop regulation according to the meritorious correction value that sets of unit to each unit.

A kind of power station automatic generation is gained merit force control method, including following operation:

1) Core Feature of power station automatic generation real power control be divided in two functional blocks of AGC, PID process, It is separated to non-core function in plug-in functional block process；

2) AGC, pid function block receive the pretreatment values that device parameter, dynamic parameter and pretreatment module send respectively；

3) AGC functional block is set up the combination of hydropower station unit according to received data and is exerted oneself model and fortune is combined in calculating Row district, associating optimized operation district, combine restriction Operational Zone；Gain merit apportioning cost in associating optimized operation district according to full station AGC Territory or combining restriction operation area, start corresponding AGC allocation flow obtain each put into AGC unit unit AGC Meritorious apportioning cost is also issued to unit；

Put into each generating set of AGC according to gain merit apportioning cost and unit PID function of the unit AGC assigned to throw and move back State and AGC throw and move back state, determine that unit is gained merit setting value；Again each generating set unit frequency modulation correction is gained merit superposition Gain merit setting value to unit, show that unit is meritorious and set correction value and be sent to pid function block；

4) each generating set PID performs equipment and sets correction value according to unit is meritorious, carries out PID closed loop regulation.

The combination that described AGC functional block is set up model of exerting oneself includes:

The optimum combination putting into AGC generating set is exerted oneself model, utilizes this model according to putting into AGC generating set respectively Plant the combination optimized operation district under Optimal Distribution combination, determine that input AGC generating set is in associating optimized operation district Respectively exert oneself the available Optimal Distribution combination under interval；

Model of exerting oneself is combined in the restriction putting into AGC generating set, utilizes this model according to putting into AGC generating set respectively Combined running district under kind of distributed combination mode, determines that input AGC generating set is respectively exerted oneself district within combining restriction Operational Zone Available distributed combination mode under between.

The unit AGC of described each input AGC generating set gains merit being retrieved as of apportioning cost:

1) calculate full station AGC to gain merit apportioning cost, and judge to be in and put into the associating optimized operation district of AGC generating set also It is to combine restriction Operational Zone；

2) AGC allocation flow is started when meeting following condition: meritorious point of the unit AGC of all input AGC units Join value summation and full station AGC and gain merit the difference of apportioning cost more than AGC meritorious regulation dead band；Or put into AGC unit Combination exert oneself model or cooperation district, associating optimized operation district, combine restriction Operational Zone and change；

3) when carrying out AGC distribution, if the AGC that entirely stands gains merit, apportioning cost is optimum in the associating putting into AGC generating set Operational Zone, then start optimum allocation algorithm；If full AGC apportioning cost of gaining merit in station combines limit put into AGC generating set Operational Zone processed, then start and limit allocation algorithm；

4) calculate and put into the unit AGC of each generating set of AGC and gain merit apportioning cost be issued to generating set.

Described PID closed loop regulation, makes the meritorious real value of sending out of unit trend towards the meritorious correction value that sets of unit, and is maintained at dead In the range of district, when shaping modes is speed regulator aperture pattern, it is computer supervisory control system machine set PLC that PID performs equipment； When shaping modes is speed regulator power mode, it is speed regulator that PID performs equipment.

Compared with prior art, the present invention has a following useful technique effect:

1, the power station automatic generation of the present invention gain merit exert oneself control system and method, have employed pre-build input AGC The combination of unit is exerted oneself the mode of model, the data existed with unit form is converted to the data of power station form, improves and ask Solve efficiency, and calculate associating optimized operation district, the cooperation district in power station, combine the parameters such as restriction Operational Zone, with this The reference of intuitive is provided for operations staff.

2, the power station automatic generation of the present invention is gained merit the system and method controlled of exerting oneself, and gains merit apportioning cost at unit AGC The distribution method of multi-form is used according to the meritorious desired value of unit output model and full station, it is sufficient to meet and shake more during calculating Move the most meritorious demand for control of exerting oneself of the different large medium-size stations of district, multicomputer type.

3, the power station automatic generation of the present invention gain merit exert oneself control system and method, add in pid function block and make Gain merit with unit frequency modulation correction and revise unit and gain merit the way of setting value, drawn by the adjusted open loop vector that primary frequency modulation is caused Enter PID closed loop regulation target, solve rushing of unit primary frequency modulation and AGC, PID function under various real power control pattern Prominent problem, it is achieved that power station primary frequency modulation superposes with the preferable of frequency modulation frequency modulation.

4, the power station automatic generation of the present invention gain merit exert oneself control system and method, considered AGC of Hydropower Stations, PID, the functional characteristics of primary frequency modulation, cover power station gain merit exert oneself control all links, for power station real power control Cooperating of each functional block provides Unified Algorithm, it is therefore prevented that produce between the internal various strategies of function does not mates or conflict Problem, also solves the most ever-present, power station revise AGC or PID strategy time, ignore strategy and parameter it Between global configuration, partial strategy and parameter are carried out isolated amendment simple and crudely, thus add power station automatic generation The problem controlling the unordered degree of function.

5, the present invention solves active power distribution, evades Hydropower Plant vibrating area (or claiming to forbid Operational Zone), compatibility one The problems such as secondary frequency modulation, carry out function division, parameter combing, strategy design, and make merit power station automatic generation real power control Can exchange, by data, the method coordinated and form an organic whole with strategy between block.

Accompanying drawing explanation

Fig. 1 is each functional block structure and the schematic diagram of data stream of the present invention；

Fig. 2 is that the power station automatic generation of the present invention is gained merit the main flow chart controlled of exerting oneself；

Fig. 3 is the meritorious real flow chart sending out value pretreatment of unit of the present invention；

Fig. 4 is the flow chart of the head value pretreatment of the present invention；

Fig. 5 is the flow chart of the full station meritorious desired value pretreatment of the present invention；

Fig. 6 is that the unit of the present invention is gained merit the flow chart of setting value pretreatment；

Fig. 7 is the flow chart that the AGC of the present invention is sequentially allocated pattern；

Fig. 8 is that the unit that calculates of the present invention is gained merit the flow chart of setting value.

Detailed description of the invention

Below in conjunction with specific embodiment, the present invention is described in further detail, described in be explanation of the invention rather than Limit.

Seeing Fig. 1, Fig. 2, a kind of power station automatic generation is gained merit Force control system, including AGC functional block, PID merit Energy block and preprocessing function block, device parameter, dynamic parameter directly input or are sent to AGC after the process of preprocessed functional block Functional block, pid function block update；

AGC functional block, pid function block receive the pretreatment values that preprocessing function block sends respectively；

AGC functional block is set up the combination of hydropower station unit according to received data and is exerted oneself model, and according to mould of exerting oneself Type calculates cooperation region, associating optimized operation region and combines restriction operation area；The distribution of work is had according to full station AGC Value is in associating optimized operation region or combines restriction operation area, starts corresponding AGC allocation flow and obtains each input The unit AGC of AGC unit gains merit apportioning cost be issued to unit；

Put into each generating set of AGC further according to unit AGC gain merit apportioning cost and unit PID function throw state of moving back and AGC throws and moves back state, determines that unit is gained merit setting value；Each the meritorious unit that is superimposed to of generating set unit frequency modulation correction is gained merit Setting value, show that unit is meritorious and sets correction value and be sent to pid function block；

Pid function block carries out PID closed loop regulation according to the meritorious correction value that sets of unit to each unit.

A kind of power station automatic generation is gained merit force control method, including following operation:

1) Core Feature of power station automatic generation real power control be divided in two functional blocks of AGC, PID process, It is separated to non-core function in plug-in functional block process；

2) AGC, pid function block receive the pretreatment values that device parameter, dynamic parameter and pretreatment module send respectively； AGC, pid function block calculate according to dynamic parameter and device parameter and obtain intermediate parameters；

3) AGC functional block is set up the combination of hydropower station unit according to received data and is exerted oneself model and fortune is combined in calculating Row district, associating optimized operation district, combine restriction Operational Zone；Gain merit apportioning cost in associating optimized operation district according to full station AGC Territory or combine restriction operation area, start corresponding AGC allocation flow obtain each put into AGC unit unit AGC Meritorious apportioning cost is also issued to unit；

Put into each generating set of AGC according to gain merit apportioning cost and unit PID function of the unit AGC assigned to throw and move back State and AGC throw and move back state, determine that unit is gained merit setting value；Again each generating set unit frequency modulation correction is gained merit superposition Gain merit setting value to unit, show that unit is meritorious and set correction value and be sent to pid function block；

4) each generating set PID performs equipment and sets correction value according to unit is meritorious, carries out PID closed loop regulation.

The pretreatment values that described AGC is received include the pretreatment values of head value, the pretreatment values of the meritorious desired value in full station, Unit is gained merit the pretreatment values of setting value and the pretreatment values of unit bad condition running priority level；The pretreatment that PID is received Value includes that unit is gained merit the pretreatment values of setting value and the meritorious real pretreatment values sending out value of unit；

Described device parameter is the parameter characteristic of each generating set in power station, and the device parameter that AGC functional block receives includes: Under different heads, under the unit output upper limit, different head, under unit each unit vibrating area scope, different head, each unit of unit limits Operational Zone processed scope；The device parameter that pid function block receives includes: unit meritorious regulation dead band, unit be meritorious set difference limit, Unit meritorious sudden change boundary, coefficient of frequency modulation；

Described dynamic parameter is equipment real-time status or operation conditions real-time parameter；

The calculating of described intermediate parameters includes: AGC meritorious regulation dead band, AGC meritorious setting difference limit, AGC gain merit Sudden change boundary, unit AGC have distribution of work step-length, unit frequency modulation correction to gain meritWith each unit operation of unit under different heads District.

The combination that described AGC functional block is set up model of exerting oneself includes:

The optimum combination putting into AGC generating set is exerted oneself model, utilizes this model according to putting into AGC generating set respectively Plant the combination optimized operation district under Optimal Distribution combination, determine that input AGC generating set is in associating optimized operation district Respectively exert oneself the available Optimal Distribution combination under interval；

Model of exerting oneself is combined in the restriction putting into AGC generating set, utilizes this model according to putting into AGC generating set respectively Combined running district under kind of distributed combination mode, determines that input AGC generating set is respectively exerted oneself district within combining restriction Operational Zone Available distributed combination mode under between.

Described input AGC generating set is respectively exerted oneself the available Optimal Distribution combination side under interval in associating optimized operation district Being retrieved as of formula:

1) Operational Zone scope is limited, to throwing according to the upper limit of exerting oneself under each head of unit, each unit vibrating area scope, each unit The generating set entering AGC is grouped；

2) under current head, Operational Zone scope, the upper limit of exerting oneself, meter are limited according to unit unit vibrating area scope, unit Calculate the unit optimized operation district of each group of generating set；

3) for each group of generating set, exert oneself in the distribution situation in each unit optimized operation district according to generating set, calculate each Group generating set combination optimized operation district under various Optimal Distribution modes；

4) for the generating set of all input AGC, the difference in unit optimized operation district it is according to each group of generating set Distribution mode and each combination optimized operation district organizing generating set of correspondence, calculate owning under various Optimal Distribution combination Put into the combination optimized operation district of AGC generating set；

5) according to putting into the AGC unit union combining optimized operation district under all Optimal Distribution combinations, calculate Put into the associating optimized operation district of AGC generating set；

6) according to putting into AGC generating set combination optimized operation district under various Optimal Distribution combinations, throwing is determined Enter the available Optimal Distribution combination that AGC generating set is respectively exerted oneself under interval in associating optimized operation district.

Obtaining of the available distributed combination mode that described input AGC generating set is respectively exerted oneself under interval within combining restriction Operational Zone It is taken as:

1) Operational Zone scope is limited, to throwing according to the upper limit of exerting oneself under each head of unit, each unit vibrating area scope, each unit The generating set entering AGC is grouped；

2) under current head, according to unit unit vibrating area, the upper limit of exerting oneself, the unit operation of each group of generating set is calculated District；

3) for each group of generating set, exert oneself in the distribution situation in each unit operation district according to generating set, calculate each group and send out Group of motors combined running district under various distribution modes；

4) for the generating set of all input AGC, the different distributions in unit operation district it is according to each group of generating set Mode and each combined running district organizing generating set of correspondence, all input AGC calculated under various distributed combination mode send out The combined running district of group of motors；

5) calculate the cooperation district putting into AGC generating set and combine restriction Operational Zone；

6) according to putting into AGC generating set combined running district under various distributed combination modes, determine that input AGC sends out The available distributed combination mode that group of motors is respectively exerted oneself under interval within combining restriction Operational Zone.

The unit AGC of described each input AGC generating set gains merit being retrieved as of apportioning cost:

1) calculate full station AGC to gain merit apportioning cost, and judge to be in and put into the associating optimized operation district of AGC generating set also It is to combine restriction Operational Zone；

2) AGC allocation flow is started when meeting following condition: meritorious point of the unit AGC of all input AGC units Join value summation and full station AGC and gain merit the difference of apportioning cost more than AGC meritorious regulation dead band；Or put into AGC unit Combination exert oneself model or cooperation district, associating optimized operation district, combine restriction Operational Zone and change；

3) when carrying out AGC distribution, if the AGC that entirely stands gains merit, apportioning cost is optimum in the associating putting into AGC generating set Operational Zone, then start optimum allocation algorithm；If full AGC apportioning cost of gaining merit in station combines limit put into AGC generating set Operational Zone processed, then start and limit allocation algorithm；

4) calculate and put into the unit AGC of each generating set of AGC and gain merit apportioning cost be issued to generating set.

Described AGC distribution, when starting optimum allocation algorithm, determines the target distribution combination putting into AGC generating set Mode and target are exerted oneself combination；Under target exerts oneself combination, if not having the unit operation district of generating set to become Change then to use and be sequentially allocated pattern；If the unit operation district having generating set changes, then use proportional allocations pattern or Use and transregional be sequentially allocated pattern；

Described AGC distribution, when starting restriction allocation algorithm, determines the target distribution combination putting into AGC generating set Mode and target are exerted oneself combination；Under target exerts oneself combination, if not having the unit operation district of generating set to become Change, be sequentially allocated pattern have the unit operation district of generating set to change if then using, can use proportional allocations pattern or Use and transregional be sequentially allocated pattern；

Described unit AGC gains merit being retrieved as of apportioning cost: calculate unit AGC active power adjustment value；If unit AGC has Merit adjusted value has distribution of work step-length less than or equal to unit AGC, and desired value of directly being gained merit by unit AGC is as unit AGC Meritorious apportioning cost, and it is issued to generating set；Otherwise, desired value of being gained merit by unit AGC substep is converted into unit AGC Meritorious apportioning cost, and it is issued to generating set.

Described unit is gained merit being retrieved as of setting value:

1) when unit PID function exits, take the meritorious real value of sending out of unit unit and gain merit setting value as unit；

2) when unit PID function puts into and unit does not puts into AGC, taking unit setting value of gaining merit has as unit AGC Distribution of work value, and is less than or equal to if manually setting the value real difference sending out value meritorious with unit in the range of unit meritorious setting difference limit The unit output upper limit, then take operations staff and manually set value and gain merit setting value as unit, otherwise then keeps unit to gain merit setting value Constant；

3) as unit PID function input and unit commitment AGC, if unit AGC gains merit, apportioning cost is meritorious with unit real Send out the difference being worth and and be less than or equal to the unit output upper limit in the range of unit meritorious setting difference limit, then take unit AGC and gain merit Apportioning cost is gained merit setting value as unit, otherwise then keeps unit setting value of gaining merit constant.

Concrete include following operating procedure:

S1000) reject the higher unnecessary function of complexity, and automatic generation real power control is divided into AGC, PID and Plug-in functional block；

S2000) parameter and the data stream of each functional block are determined；

S3000) the input parameter of each functional block is pre-processed；

S4000) the abnormal exit strategy of each functional block is determined,

S5000) set up the combination of Hydropower Plant to exert oneself model calculate cooperation district, associating optimized operation district, associating Limit Operational Zone；

S6000) calculate and put into the unit AGC of AGC unit and gain merit apportioning cost；

S7000) according to put into each unit of AGC unit AGC gain merit apportioning cost and and unit PID function throw move back state, AGC throws and moves back state etc., determines that unit is gained merit setting value；

S8000) meritorious for each the unit unit frequency modulation correction unit that is superimposed to is gained merit setting value, show that unit setting of gaining merit is repaiied On the occasion of；

S9000) each unit PID performs equipment and sets correction value according to unit is meritorious, carries out PID closed loop regulation, makes list The meritorious real value of sending out of machine trends towards unit meritorious setting correction value, and is maintained in dead zone range.

1, AGC functional block, the setting up of pid function block

The power station automatic generation that the present invention provides is gained merit Force system and control method, be reject complexity higher need not Want function, and automatic generation real power control is divided into AGC, PID and plug-in functional block, be separately operable in plant stand level control In device processed and unit level controller.

In order to reduce program complexity, the present invention does not consider that rate of water consumed in electric generator Curve Matching, the automatic startup-shutdown of unit etc. are the most multiple Miscellaneous function, and be separated to non-core function in plug-in functional block process, to promote Procedure Haleness.

In order to make the complexity of each functional block be controlled, need each functional block is undertaken in automatic generation real power control Function carry out classifying rationally.Do not produce horizontal function between the pid function block of each unit to rely on or data exchange；To the greatest extent Function may be distributed in the pid function block of each unit, reduce the computational load of AGC functional block；Reduce AGC and The computational load of pid function block, AGC functional block and pid function block avoid participating in the pretreatment of input parameter as far as possible.

1.1, plant stand level functional block:

1) AGC functional block: intermediate parameters calculates, Unit Combination goes out force modeling, cooperation district calculates, meritorious point of unit Join, abnormal exit strategy etc.；

2) head value preprocessing function block；

3) the meritorious desired value preprocessing function block in full station；

4) unit bad condition running priority level preprocessing function block.

1.2, unit level functional block:

1) pid function block: intermediate parameters calculating, meritorious closed loop regulation, abnormal exit strategy etc.；

2) meritorious real the sending out of unit is worth preprocessing function block；

3) unit is gained merit setting value preprocessing function block.

2, the structure of each functional block and data stream

See Fig. 1, the structure of each functional block that present invention determine that and data stream, including:

S2100) dynamic parameter is determined, i.e. consersion unit real-time status or the parameter that inputted in real time by operations staff, including: The unit of each unit is meritorious real sends out value, setting value gained merit by head value, the meritorious setting value in full station, the unit of each unit, respectively The state of platform unit, each unit have zero defects, Quan Chang to have zero defects, frequency departure, each unit bad condition running priority Level；

S2200) determining device parameter, be different from dynamic parameter, the characteristic that device parameter embodies power station difference unit is poor Different, typically arranged according to dispatch from the factory parameter or result of the test of unit by attendant, and the loading when program starts, including: no Limit with each unit of unit under unit each unit vibrating area scope, different head under the unit output upper limit, different head under head Operational Zone scope, unit meritorious regulation dead band, unit is meritorious sets difference limit, unit meritorious sudden change boundary, coefficient of frequency modulation etc.； The setting principle of device parameter is:

2210) the unit output upper limit under different heads, determines or by testing according to the meritorious control characteristic of unit under different heads Draw；

2220) unit each unit vibrating area scope under different heads, determines according to the stable characteristic of unit under different heads Or drawn by test；

2230) under different heads, each unit of unit limits Operational Zone scope, according to the stable characteristic of unit under different heads Determine or drawn by test；

2240) the gain merit ratio of regulation dead band and the unit output upper limit of unit unit is necessarily less than or is equal to the AGC of electrical network and examines Core precision, the meritorious control characteristic further according to water turbine set is determined

I.e. unit unit meritorious regulation dead bandMeet:

Wherein ε is electrical network AGC assessment precision, and i is machine group number,For the unit output upper limit；

2250) the meritorious difference limit that sets of unit unit must be considerably larger than unit meritorious regulation dead band: unit unit is meritorious sets difference limitMeet:

2260) unit unit meritorious sudden change boundary have to be larger than unit meritorious regulation difference limit: unit unit meritorious sudden change boundaryMeet:

2270) coefficient of frequency modulation should be greater than or examine the coefficient of correspondence of formula equal to primary frequency regulation of power network device electricity, to ensure The regulation of PID closed loop can the meritorious adjusted open loop deviation that produces of compatible primary frequency modulation.

S2300) determine intermediate parameters, i.e. calculate the parameter of gained according to dynamic parameter and device parameter, including: AGC has Merit regulation dead band, AGC meritorious setting difference limit, AGC meritorious sudden change boundary, unit AGC have distribution of work step-length, unit Frequency modulation correction is meritorious, unit each unit operation district model under unit each unit optimized operation district scope, different head under different heads Enclose；The computing formula of intermediate parameters and follow principle, including:

S2310) AGC gain merit regulation dead band more than or equal to all inputs AGC unit units gain merit regulate dead band maximum Value, the summation in regulation dead band of gaining merit less than or equal to all input AGC unit units；I.e. AGC meritorious regulation dead band p_d Meet:

$max\[p_d^1\×{\mathrm{\θ}}^1,p_d^2\×{\mathrm{\θ}}^2,\mathrm{...}{p}_f^n\×{\mathrm{\θ}}^n\]\≤p_d=f\[\underset{i=1}{\overset{n}{\Σ}}(p_d^i\×{\mathrm{\θ}}^i)\]\≤\underset{i=1}{\overset{n}{\Σ}}(p_d^i\×{\mathrm{\θ}}^i)$

Wherein θⁱFor unit commitment AGC variable, as unit commitment AGC, θ value is 1, when unit exits AGC, θ value is 0.

S2320) AGC gain merit set difference limit more than or equal to all inputs AGC unit units gain merit set difference limit maximum Value, gains merit less than or equal to all input AGC unit units and sets the summation of difference limit, dead much larger than the meritorious regulation of AGC District, i.e. AGC are meritorious sets difference limit p_spMeet:

$\left\{\begin{array}{c}max\[p_{sp}^1\×{\mathrm{\θ}}^1,p_{sp}^2\×{\mathrm{\θ}}^2,\mathrm{...}{p}_{sp}^n\×{\mathrm{\θ}}^n\]\≤p_{sp}=f\[\underset{i=1}{\overset{n}{\Σ}}(p_{sp}^i\×{\mathrm{\θ}}^i)\]\≤\underset{i=1}{\overset{n}{\Σ}}(p_{sp}^i\×{\mathrm{\θ}}^i)\\ p_{sp}>>p_d\end{array}\right.$

S2330) AGC gain merit sudden change boundary more than or equal to all generating set units gain merit sudden change boundary maximum, little In or the summation of sudden change boundary of gaining merit equal to all generating set units, meritorious more than AGC set difference limit；I.e. AGC has Merit sudden change boundary p_errMeet:

$\left\{\begin{array}{c}max\[p_{err}^1,p_{err}^2,\mathrm{...}{p}_{err}^n\]\≤p_{err}=f\left(\underset{i=1}{\overset{n}{\Σ}}{p}_{err}^i\right)\≤\underset{i=1}{\overset{n}{\Σ}}{p}_{err}^i\\ p_{err}>p_{sp}\end{array}\right.$

S2340) unit AGC has distribution of work step-length set difference limit meritorious less than unit to deduct unit meritorious regulation dead band, i.e. Unit AGC has distribution of work step-lengthMeet:

S2350) when unit primary frequency modulation action and unit PID function put into, unit frequency modulation correction is meritorious inclined equal to frequency Difference and the product of coefficient of frequency modulation and this unit output upper limit under current head subtract unit and gain merit the minimum in setting value gained difference Value, otherwise the correction of unit frequency modulation is gained merit is 0, i.e. unit frequency modulation correction is gained meritMeet:

$p_f^i=\min [\mathrm{\Δf}\×k^i,{\stackrel{-}{p}}_h^i-p_{\mathrm{set}}^i]\×{\mathrm{\σ}}^i\×{\mathrm{\μ}}^i$

Wherein Δ f is frequency departure, kⁱFor unit coefficient of frequency modulation,For the unit output upper limit under current head,For unit Unit is gained merit setting value, μⁱFor unit primary frequency modulation action variable, when unit primary frequency modulation action, μ value is 1, when one When secondary frequency modulation is not operating, μ value is 0, σⁱVariable, when unit unit PID puts into, σ value is put into for unit unit PID Being 1, when unit unit PID exits, σ value is 0；

S2360) unit each unit optimized operation district and not under unit each unit operation district scope be difference head under different heads Union with each unit of unit under head restriction Operational Zone:

Operational Zone=optimized operation district U limits Operational Zone

S2400) data of exchange between data stream, i.e. each functional block are determined, including: (PID is defeated for meritorious real the value of unit Go out to AGC functional block), unit AGC gains merit, and (AGC is defeated for the unit output upper limit under apportioning cost (two-way), current head Go out to pid function block), unit frequency modulation correction meritorious (PID exports to AGC functional block) etc..

3, the parameter of input is pre-processed by each preprocessing function block, including:

S3100) the real pretreatment of sending out value meritorious to unit, flow process is as it is shown on figure 3, comprise the following steps:

S3110) priority in active power measurement sources multiple to unit is ranked up, if the highest active power measurement source of priority Measured value quality preferably, takes the measured value in the highest active power measurement source of priority as meritorious real the value of unit；

S3120) if the measured value quality in the highest active power measurement source of priority is bad, then high to priority second meritorious survey The measured value quality in amount source judges；

S3130) if the measured value quality in the high active power measurement source of priority second preferably, takes the meritorious survey that priority second is high The measured value in amount source, as the meritorious real value of sending out of unit, otherwise continues the measured value quality in the active power measurement source judging that priority the 3rd is high, And according to active power measurement source priority from high to low order successively recursion；

S3140) when the active power measurement source that recursion to priority is minimum, the measured value quality in the source of measurement is not judged, directly Access the measured value in this active power measurement source as meritorious real the value of unit.

S3200) pretreatment to head value, flow process as shown in Figure 4, comprises the following steps:

S3210) determine automatic water header value, have following steps:

1) the measured value h to all level measurement devices_iFilter, if h_i>h_maxOr h_i<h_minThen by this measured value Get rid of mark γ_iSet to 0, on the contrary γ_iValue is 1, wherein h_maxFor highest possible water level, h_minFor minimum possible water level.

2) mean value of all head measured values not being excluded is sought

$\stackrel{\&OverBar;}{h}=\underset{i=1}{\overset{n}{\&Sigma;}}(h_i\&times;{\mathrm{\&gamma;}}_i)\&divide;\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{\&gamma;}}_i$

3) ask all head measured values of not being excluded withThe maximum Δ h of absolute difference

$\mathrm{\&Delta;}h=max\&lsqb;\left|h_1-\stackrel{\&OverBar;}{h}\right|\&times;{\mathrm{\&gamma;}}_1,\left|h_2-\stackrel{\&OverBar;}{h}\right|\&times;{\mathrm{\&gamma;}}_2,...\left|h_n-\stackrel{\&OverBar;}{h}\right|\&times;{\mathrm{\&gamma;}}_n\&rsqb;$

4) if Δ h is more than acceptable head measure error scope h_err, then by the row of head measured value the highest for difference degree Except mark γ_iSet to 0；

5) 2 to 4 steps are repeated, until Δ h≤h_err；

6) intermediate value of all head measured values not being excluded is taken as automatic water header value.

S3220) determine manual water header value, comprise the following steps:

S3221) when head setting means is manual, the artificial setting value of operations staff is taken as manual water header value；

S3222) when head setting means is automatic, final head value is taken as manual water header value.

S3230) determine final head value, comprise the following steps:

S3231) when head setting means is manual, manual water header value is taken as final head value；

S3232) when head setting means is automatic, if all head measured values are not the most in effective range, the most automatically will Head setting means is cut to manually；

S3233) when head setting means is automatic, automatic water header value is compared with final head value, if difference In normal fluctuation range, then take from kinetic head value as final head value, otherwise, then head setting means is switched to hand Dynamic；

S3234) after finally head value changes, if the meritorious upper limit of exerting oneself in the full station under this head is less than meritorious real the value in station entirely, Or have the unit output upper limit real send out value meritorious less than this unit unit, then return to the final head value before changing.

S3300) pretreatment of desired value meritorious to full station, flow process is as it is shown in figure 5, comprise the following steps:

S3310) when AGC function exits, take all unit unit AGC and gain merit apportioning cost summation as the meritorious mesh in full station Scale value, electric network active setting value, collection control meritorious setting value, power plant gains merit setting value；

S3320) when AGC function input and regulation source are at electrical network, control is meritorious to be set as collecting to take the meritorious desired value in full station Value and power plant gain merit setting value, if the electric network active setting value real difference of sending out value summation meritorious with all unit units exists AGC is meritorious to be set in the range of difference limits and is in Quan Zhan cooperation district, then take electric network active setting value as the meritorious mesh in full station Scale value, on the contrary then keep the meritorious desired value in full station constant；

S3330) put into when AGC function and regulation source is when collection control, take the meritorious desired value in full station and set as electric network active Value and power plant gain merit setting value, if the meritorious setting value of the collection control real difference of sending out value summation meritorious with all unit units exists AGC is meritorious sets in the range of difference limits and is in Quan Zhan cooperation district, then take the meritorious setting value of collection control as complete meritorious mesh of standing Scale value, on the contrary then keep the meritorious desired value in full station constant；

S3340) when AGC function input and regulation source are in power plant, the meritorious desired value in full station is taken as electric network active setting Value and the meritorious setting value of collection control, if power plant gains merit, the setting value real difference of sending out value summation meritorious with all unit units exists AGC is meritorious sets in the range of difference limits and is in Quan Zhan cooperation district, then power taking factory gains merit setting value as the meritorious mesh in full station Scale value, on the contrary then keep the meritorious desired value in full station constant；

S3350) put into when AGC function and regulation source at electrical network time, if power plant and the communication system of power grids are interrupted, or electrical network without Setting value issues and exceedes certain time interval T₁, then AGC regulation source automatically switches to power plant, T₁Not less than meritorious setting value curve Issue the time interval of adjacent 2 on mode lower curve；

S3360) put into when AGC function and regulation source is when collection control, if power plant and collection control communication disruption, or collection control nothing Setting value issues and exceedes certain time interval T₂(T₂=T₁), then AGC regulation source automatically switches to power plant.

The pretreatment of setting value of S3400) unit being gained merit, flow process as shown in Figure 6, comprises the following steps:

S3410) when unit PID function exits, take the meritorious real value of sending out of unit unit and gain merit setting value as unit, take list Machine setting value of gaining merit is gained merit apportioning cost as unit AGC；

S3420) when unit PID function puts into and unit does not puts into AGC, take unit and gain merit setting value as unit AGC Meritorious apportioning cost, if operations staff manually sets value and the meritorious real difference of value of sending out of unit in the range of unit meritorious setting difference limit And less than or equal to the unit output upper limit, then take operations staff and manually set value and gain merit setting value as unit, otherwise then keep singly Machine setting value of gaining merit is constant；

S3430) as unit PID function input and unit commitment AGC, if unit AGC gains merit, apportioning cost has with unit The real difference of value of sending out of merit in the range of unit meritorious setting difference limit and less than or equal to the unit output upper limit, then takes unit AGC Meritorious apportioning cost is gained merit setting value as unit, otherwise then keeps unit setting value of gaining merit constant.

S3500) pretreatment to unit bad condition running priority level, comprises the following steps:

The operation time that S3510) each unit be in after turn(a)round last time restriction Operational Zone and vibrating area is weighted system Meter, wherein vibrating area runs the statistical weight of time far above the statistical weight limiting the Operational Zone operation time；

S3520) time after each unit weighted statistical is ranked up, the most from short to long, sets the most from high to low Fixed automatic priority；

S3530) operations staff is to each group setup manual priority level；

S3540) when priority set-up mode is automatic, dynamic priority is taken from as unit bad condition running priority level, When priority set-up mode is manual, take manual priority level as unit bad condition running priority level.

4, the exception policy that each functional block exits is determined

Specified conditions are set the most in a program, when specified conditions trigger, it is believed that power station or Hydropower Plant are in exception State, beyond program can the condition range of normal process, thus exit whole merits of power station automatic generation real power control Can or partial function transfer Artificial Control to, abnormal exit strategy is rationally set and contributes to by suitably reducing program adaptability Make program complexity be controlled, specifically include:

S4100) abnormality causing unit PID function to exit is determined

S4110) the gain merit difference of meritorious superposition value real the value meritorious with unit of setting value and unit frequency modulation correction of unit exceeds Unit meritorious regulation dead band, and continue for some time t₁:

$\left|p_{set}^i+p_f^i-p^i\right|>p_d^i$

Wherein pⁱIt is worth for meritorious real the sending out of unit.

S4120) unit gains merit setting value real the send out difference of value meritorious with unit beyond unit meritorious sudden change boundary:

$\left|p_{set}^i-p^i\right|>p_{err}^i$

S4130) there is saltus step in the meritorious real value of sending out of unit, and hop value is beyond unit meritorious sudden change boundary:

S4140) unit is in non-power generating state；

S4150) unit all active power measurements source measured value quality is bad；

S4160) unit has accident.

S4200) abnormality causing power station AGC function to exit is determined, including:

S4210) AGC performs equipment and performs equipment communication abnormality or interruption with certain unit PID；

S4220) the gain merit superposition value of summation of the meritorious desired value in full station and all unit frequency modulation corrections is gained merit with all unit units The real difference of value summation of sending out is beyond AGC meritorious regulation dead band, and continues for some time t₂(t₂>=t₁):

$\left|p_{set}+\underset{i=1}{\overset{n}{\&Sigma;}}{p}_f^i-\underset{i=1}{\overset{n}{\&Sigma;}}{p}^i\right|>p_d$

Wherein p_setFor the meritorious desired value in full station.

S4230) the meritorious desired value in full station real the send out difference of value summation meritorious with all unit units is beyond the meritorious sudden change of AGC Boundary:

$\left|p_{set}-\underset{i=1}{\overset{n}{\&Sigma;}}{p}^i\right|>p_{err}$

S4240) meritorious real the sending out of all unit units is worth summation generation saltus step, and hop value is beyond AGC meritorious sudden change boundary；

S4250) the meritorious desired value in full station is not in Quan Zhan cooperation district；

S4260) certain unit all active power measurements source measured value quality is bad；

S4270) there is accident in full factory.

S4300) abnormality causing unit to exit AGC is determined, including:

S4310) AGC function in power station is in the state of exiting；

S4320) unit PID function is in the state of exiting；

S4330) unit setting value of gaining merit is not equal to unit AGC and gains merit apportioning cost, and continues for some time t₃:

WhereinGain merit apportioning cost for the unit AGC of unit.

S4340) the meritorious real value of sending out of unit is in unit vibrating area and exceedes certain time t₄；

S4350) the meritorious real value of sending out of unit is in unit and limits Operational Zone and exceed certain time t₅(t₅>>t₄)。

5, set up the combination of Hydropower Plant to exert oneself model calculate cooperation district, associating optimized operation district, combine restriction Operational Zone, comprises the following steps:

S5100) when S2200 arranges device parameter, to each unit each unit vibrating area scope under different heads and list Machine limits Operational Zone scope and is modified；

S5200) set up and put into the optimum combination of AGC unit and exert oneself model, and calculate that to put into the associating of AGC unit optimum Operational Zone；

S5300) model of exerting oneself is combined in the restriction setting up input AGC unit, and calculates the cooperation putting into AGC unit District limits Operational Zone with combining；

S5400) the associating optimized operation district of AGC unit will be put into plus all unit AGC not putting into AGC unit Meritorious apportioning cost, associating optimized operation district of entirely being stood, the meritorious setting value in full station is set for operations staff reference is provided；

S5500) the cooperation district putting into AGC unit is not put into AGC unit with combining restriction Operational Zone plus all Unit AGC gain merit apportioning cost, obtain Quan Zhan cooperation district and complete station and combine restriction Operational Zone, for operations staff's setting The meritorious setting value in full station provides reference.

Described S5100) when S2000 arranges device parameter, to each unit each unit vibrating area scope under different heads Limit Operational Zone scope with unit to be modified, including:

S5110) dimensionality reduction correction, by suitably reducing unit regulation performance, reaches to simplify the purpose of real power control strategy；

S5111) unit unit vibrating area and unit are limited Operational Zone different heads relatively and carry out merger process, its The merger priority of middle unit vibrating area limits Operational Zone higher than unit.Such as power station unit when 204 meters of head, unit Vibrating area is 0 to 120MW, 210MW to 420MW, and it is 120MW to 210MW that unit limits Operational Zone, During 206 meters of head, unit vibrating area is 0 to 120MW, 240MW to 480MW, and it is 120MW that unit limits Operational Zone To 240MW, then the unit vibrating area of 204 meters to 206 meters of head and unit can be limited in the case of allowing and run District carries out merger process, and after merger, unit vibrating area is 0 to 120MW, 210MW to 480MW, and unit limits and runs District is 120MW to 210MW；

S5112) unit vibrating area and unit are limited Operational Zone different units relatively and carry out merger process, method with Step is with reference to S5111.

S5120) correction is expanded, the meritorious deviation that regulation dead band or primary frequency modulation produce because unit is gained merit in order to prevent unit Slip into vibrating area or limit Operational Zone, the scope of vibrating area and restriction Operational Zone suitably being expanded, the wherein expansion of unit vibrating area Big priority limits Operational Zone higher than unit.Such as power station unit when 204 meters of head, unit vibrating area be 0 to 120MW, it is 120MW to 210MW that unit limits Operational Zone, then can expand correction unit vibrating area is 0 to 130MW, Expand correction unit limiting Operational Zone is 130MW to 220MW.

Described S5200) set up and put into the optimum combination of unit of AGC and exert oneself model, and calculate and put into AGC unit Associating optimized operation district, including:

S5210) Operational Zone scope is limited according to the upper limit of exerting oneself under each head of unit, each unit vibrating area scope, each unit, The unit putting into AGC is grouped, can be by the method for dimensionality reduction merger approximation parameters by different types of when packet Unit merger is same type, it is also possible to from the angle simplifying algorithm (but may reduce executing efficiency), will be every Platform unit is individually grouped；

S5220) under current head, Operational Zone scope, the upper limit of exerting oneself are limited according to unit unit vibrating area scope, unit, Calculate the unit optimized operation district of each group of unit:

Unit optimized operation district=[0, the upper limit of exerting oneself]-unit limits Operational Zone-unit vibrating area

Such as certain unit is exerted oneself upper limit 700MW under current head, unit vibrating area be 0 to 120MW, 210MW extremely 420MW, it is 120MW to 210MW that unit limits Operational Zone, then unit optimized operation district is 420MW to 700MW.

S5230) for each group of unit, according to unit output in the distribution situation in each unit optimized operation district, each group of machine is calculated Group combination optimized operation district under various Optimal Distribution modes, as a example by a certain group of unit:

S5231) list this group unit to gain merit all Optimal Distribution modes in each unit optimized operation district of exerting oneself, it is assumed that this point Group has n platform unit, every unit to have m unit optimized operation district, and the most all of Optimal Distribution mode is following equations group Disaggregation:

$\left\{\begin{array}{c}{x}_1+x_2+...+x_m=n\\ x_1,x_2,...,x_m\&le;n\\ x_1,x_2,...,x_m\&Element;N\end{array}\right.$

Wherein x₁,x₂,…,x_mRepresent successively and be in the 1st, 2 ..., the unit number of units in m unit optimized operation district, by equation group Solution represent with the form of m dimensional vector, it is assumed that total n_AKind Optimal Distribution mode:

$a_1=\left[\begin{array}{c}{a}_{11}\\ a_{21}\\ .\\ .\\ .\\ a_{m1}\end{array}\right]=\left[\begin{array}{c}{x}_1^1\\ x_2^1\\ .\\ .\\ .\\ x_m^1\end{array}\right],a_2=\left[\begin{array}{c}{a}_{12}\\ a_{22}\\ .\\ .\\ .\\ a_{m2}\end{array}\right]=\left[\begin{array}{c}{x}_1^2\\ x_2^2\\ .\\ .\\ .\\ x_m^2\end{array}\right],\mathrm{...}{a}_{n_A}=\left[\begin{array}{c}{a}_{1n_A}\\ a_{2n_A}\\ .\\ .\\ .\\ a_{{\mathrm{mn}}_A}\end{array}\right]=\left[\begin{array}{c}{x}_1^{n_A}\\ x_2^{n_A}\\ .\\ .\\ .\\ x_m^{n_A}\end{array}\right]$

Then solution of equations collection can be expressed as matrix a:

$a=\&lsqb;a_1,a_2,...a_{n_A}\&rsqb;=\left[\begin{array}{cccc}{a}_{11}& a_{12}& \mathrm{...}& a_{{\ln }_A}\\ a_{21}& a_{22}& \mathrm{...}& a_{2n_A}\\ .& .& .& .\\ .& .& .& .\\ .& .& .& .\\ a_{m1}& a_{m2}& \mathrm{...}& a_{{\mathrm{mn}}_A}\end{array}\right]$

Calculating this group unit combination optimized operation district under every kind of Optimal Distribution mode, all of combination optimized operation district is permissible It is expressed as row vector A, the vector element of the set A in optimized operation districtSet a with Optimal Distribution mode Vector elementIn correspondence with each other:

$A=\&lsqb;(\underset{\&OverBar;}{p_1},\stackrel{\&OverBar;}{p_1}),(\underset{\&OverBar;}{p_2},\stackrel{\&OverBar;}{p_2}),\mathrm{...},(\underset{\&OverBar;}{p_m},\stackrel{\&OverBar;}{p_m})\&rsqb;\&times;a=\&lsqb;A_1,A_2,...A_{n_A}\&rsqb;$

Whereinp ₁ 、p ₂ 、…p _m It is followed successively by this group unit the 1st, the lower limit in m the unit optimized operation district in 2 ..., It is followed successively by this group unit the 1st, the upper limit in m the unit optimized operation district in 2 ....

S5240) for the unit of all input AGC, it is in the different optimum of unit optimized operation district according to each group of unit Distribution mode and each combination optimized operation district organizing unit of correspondence, calculate all inputs under various Optimal Distribution combination The combination optimized operation district of AGC unit, including:

S5241) all Optimal Distribution combinations of input AGC unit are enumerated, it is assumed that the unit putting into AGC is divided into A, b, c ... group, Optimal Distribution mode number is respectively n_A,n_B,n_C..., comprise the following steps:

1) by a, b, c ... set a, b, c of the Optimal Distribution mode of group unit ... be converted to the collection of Optimal Distribution mode sequence number Close Δ a, Δ b, Δ c ...:

$\begin{array}{c}a=\&lsqb;a_1,a_2,\mathrm{...},a_{n_A}\&rsqb;\&DoubleRightArrow;\mathrm{\&Delta;}a\&lsqb;1,2,\mathrm{...},n_A\&rsqb;\\ b=\&lsqb;b_1,b_2,\mathrm{...},b_{n_B}\&rsqb;\&DoubleRightArrow;\mathrm{\&Delta;}b\&lsqb;1,2,\mathrm{...},n_B\&rsqb;\\ c=\&lsqb;c_1,c_2,\mathrm{...},c_{n_C}\&rsqb;\&DoubleRightArrow;\mathrm{\&Delta;}c\&lsqb;1,2,\mathrm{...},n_C\&rsqb;\\ .\\ .\\ .\end{array}$

2) calculate matrix delta ab of all Optimal Distribution combination elements composition of a group, b group, and expand into row vector Δ ab:

$\mathrm{\&delta;}ab={\mathrm{\&Delta;a}}^T\&lsqb;1,1,\mathrm{...},1\&rsqb;\&times;10+\left[\begin{array}{c}1\\ 1\\ .\\ .\\ .\\ 1\end{array}\right]\mathrm{\&Delta;}b$

Δ ab={ [1,0 ..., 0] × δ ab, [0,1 ..., 0] × δ ab ..., [0,0 ..., 1] × δ ab}

3) calculate matrix delta abc of all Optimal Distribution combination elements composition of a group, b group, c group, and expand into Row vector Δ abc:

$\mathrm{\&delta;}abc={\mathrm{\&Delta;ab}}^T\&lsqb;1,1,...,1\&rsqb;\&times;10+\left[\begin{array}{c}1\\ 1\\ .\\ .\\ .\\ 1\end{array}\right]\mathrm{\&Delta;}c$

Δ abc={ [1,0 ..., 0] × δ abc, [0,1 ..., 0] × δ abc ..., [0,0 ..., 1] × δ abc}

4) above procedure is continued, until drawing the row vector of all Optimal Distribution combinations composition putting into AGC unit Δ abc ..., in row vector, each element digital is from a high position to low level, is followed successively by a, b, the Optimal Distribution of c ... group unit Mode sequence number.

S5242) input AGC unit combination optimized operation district under every kind of Optimal Distribution combination is calculated, including:

1) matrix delta of a group, the b group combination optimized operation p-block element p composition under all Optimal Distribution combinations is calculated AB, and expand into row vector AB:

$\mathrm{\&delta;}AB=A^T\&lsqb;1,1,...,1\&rsqb;+\left[\begin{array}{c}1\\ 1\\ .\\ .\\ .\\ 1\end{array}\right]B$

AB={ [1,0 ..., 0] × δ AB, [0,1 ..., 0] × δ AB ..., [0,0 ..., 1] × δ AB}

2) square of a group, b group, the c group combination optimized operation p-block element p composition under all Optimal Distribution combinations is calculated Battle array δ ABC, and expand into row vector ABC:

$\mathrm{\&delta;}ABC={\mathrm{AB}}^T\&lsqb;1,1,...,1\&rsqb;+\left[\begin{array}{c}1\\ 1\\ .\\ .\\ .\\ 1\end{array}\right]C$

ABC={ [1,0 ..., 0] × δ ABC, [0,1 ..., 0] × δ ABC ..., [0,0 ..., 1] × δ ABC}

3) above procedure is continued, until drawing the combination optimum putting into AGC unit under all Optimal Distribution combinations Row vector ω of Operational Zone composition, and the row vector Δ abc formed with Optimal Distribution combination ... in correspondence with each other.

S5250) according to putting into the AGC unit union combining optimized operation district under all Optimal Distribution combinations, Draw put into AGC unit associating optimized operation district Ω:

$\mathrm{\&omega;}=\&lsqb;{\mathrm{\&omega;}}_1,{\mathrm{\&omega;}}_2,{\mathrm{\&omega;}}_3...\&rsqb;=\&lsqb;(\underset{\&OverBar;}{{\mathrm{\&omega;}}_1},\stackrel{\&OverBar;}{{\mathrm{\&omega;}}_1}),(\underset{\&OverBar;}{{\mathrm{\&omega;}}_2},\stackrel{\&OverBar;}{{\mathrm{\&omega;}}_2}),(\underset{\&OverBar;}{{\mathrm{\&omega;}}_3},\stackrel{\&OverBar;}{{\mathrm{\&omega;}}_3}),...\&rsqb;$

$\mathrm{\&Omega;}=(\underset{\&OverBar;}{{\mathrm{\&omega;}}_1},\stackrel{\&OverBar;}{{\mathrm{\&omega;}}_1})\&cup;(\underset{\&OverBar;}{{\mathrm{\&omega;}}_2},\stackrel{\&OverBar;}{{\mathrm{\&omega;}}_2})\&cup;(\underset{\&OverBar;}{{\mathrm{\&omega;}}_3},\stackrel{\&OverBar;}{{\mathrm{\&omega;}}_3})\&cup;...$

S5260) according to putting into AGC unit combination optimized operation district under various Optimal Distribution combinations, throwing is determined Enter the available Optimal Distribution combination that AGC unit is respectively exerted oneself under interval in associating optimized operation district, comprise the following steps:

S5261) bound of optimum combination Operational Zone corresponding for every kind of Optimal Distribution mode described in S5242 is entered by size Row sequence, identical then the merging of size；

S5262) according to the bound after sequence, the associating optimized operation district putting into AGC unit is split, draw many Individual interval of exerting oneself, it is assumed that interval number of exerting oneself is n:

$\mathrm{\&Omega;}=\&lsqb;(\underset{\&OverBar;}{{\mathrm{\&Omega;}}_1},\stackrel{\&OverBar;}{{\mathrm{\&Omega;}}_1}),(\underset{\&OverBar;}{{\mathrm{\&Omega;}}_2},\stackrel{\&OverBar;}{{\mathrm{\&Omega;}}_2}),...,(\underset{\&OverBar;}{{\mathrm{\&Omega;}}_n},\stackrel{\&OverBar;}{{\mathrm{\&Omega;}}_n})\&rsqb;$

S5263) by combination optimized operation district that is interval and that put into corresponding to AGC various Optimal Distribution combination of respectively exerting oneself Contrast, draw the available Optimal Distribution combination under interval of respectively exerting oneself, it is assumed that Optimal Distribution combination has m Kind, corresponding combination optimized operation district is respectively as follows:

ω=[ω₁,ω₂,…,ω_m]

Then:

$\left\{\begin{array}{cc}{f}_i^j=1,& \begin{array}{cc}if& (\underset{\&OverBar;}{{\mathrm{\&Omega;}}_i},\stackrel{\&OverBar;}{{\mathrm{\&Omega;}}_i})\&Subset;{\mathrm{\&omega;}}_j\end{array}\\ f_i^j=0,& \begin{array}{cc}if& (\underset{\&OverBar;}{{\mathrm{\&Omega;}}_i},\stackrel{\&OverBar;}{{\mathrm{\&Omega;}}_i})\&NotSubset;{\mathrm{\&omega;}}_j\end{array}\\ i\&Element;(1,2,\mathrm{...},n),& j\&Element;(1,2,\mathrm{...},m)\\ f_i=\&lsqb;f_i^1,f_i^2,\mathrm{...},f_i^m\&rsqb;& \end{array}\right.$

Wherein f_iIt is the row vector of the available Optimal Distribution combination representing that i-th exerts oneself under interval, described in S5241 Row vector Δ abc ... corresponding, when row vector f_iElement value when being 1, corresponding row vector Δ abc ... representated by middle element Optimal Distribution combination can use.

Concrete, 1) assume that the unit of all input AGC is divided into 2 groups, often 2 units of group, 2 lists of every unit Machine Operational Zone, the unit operation district of the 1st group of unit is 100MW to 200MW, 300MW to 400MW, Qi Zhongqian Person be unit limit Operational Zone, the latter is unit optimized operation district, the unit operation district of the 2nd group of unit be 50MW extremely 150MW, 300MW to 350MW, both unit optimized operation district, then Optimal Distribution combination includes:

{[(0,2),(2,0)],[(0,2),(1,1)],[(0,2),(0,2)]}

2) the combined running district that every kind of distributed combination mode is corresponding is:

{(300,400)×2+(50,150)×2,(300,400)×2+(50,150)+(300,350),

(300,400) × 2+ (300,350) × 2}={ (700,1100), (950,1300), (1200,1500) }

3) union that cooperation district is all combined running districts of AGC unit, 700MW to 1500MW are put into.

4) bound putting into all combination optimized operation districts of AGC unit is ranked up {700,950,1100,1200,1300,1500}；

5) according to sequence after bound to put into AGC unit associating optimized operation district split, draw multiple go out Power is interval:

{(700,950),(950,1100),(1100,1200),(1200,1300),(1300,1500)}

6) the combination optimized operation district respectively exerted oneself under interval and various Optimal Distribution combination is compared, draw input AGC unit is respectively exerted oneself the available Optimal Distribution combination under interval in associating optimized operation district, is respectively

{[(0,2),(2,0)]},{[(0,2),(2,0)],[(0,2),(1,1)]},{[(0,2),(1,1)]},

{[(0,2),(1,1)],[(0,2),(0,2)]},{(0,2),(0,2)]}

Described S5300) set up and put into the restriction combination of AGC unit and exert oneself model, and calculate the connection putting into AGC unit Close Operational Zone and combine restriction Operational Zone, comprising the following steps:

S5310) according to the packet mode of S5210, the unit putting into AGC is grouped.

S5320) under current head, according to unit unit vibrating area, the upper limit of exerting oneself, the unit operation of each group of unit is calculated District: unit operation district=[0, the upper limit of exerting oneself]-unit vibrating area

Such as certain unit is exerted oneself upper limit 700MW under current head, unit vibrating area be 0 to 120MW, 210MW extremely 420MW, then unit operation district is 120MW to 210MW, 420MW to 700MW.

S5330) for each group of unit, according to the unit output distribution situation in each unit operation district, calculate each group of unit and exist Combined running district under various distribution modes, method and step are with reference to S5230.

S5340) for the unit of all input AGC, the different distributions mode in unit operation district it is according to each group of unit And each combined running district organizing unit of correspondence, calculate the combination of all input AGC units under various distributed combination mode Operational Zone, method and step are with reference to S5240；

S5350) calculate the cooperation district putting into AGC unit and combine restriction Operational Zone, comprise the following steps:

S5351) according to the union in input AGC unit combined running district under all distributed combination modes, input is drawn The cooperation district of AGC unit, method and step are with reference to S5250；

S5352) from the cooperation district putting into AGC unit, deduct the associating optimized operation district that S5250 draws, draw Put into AGC unit combines restriction Operational Zone.

S5360) according to putting into AGC unit combined running district under various distributed combination modes, input AGC machine is determined The available distributed combination mode that group is respectively exerted oneself under interval within combining restriction Operational Zone, method and step are with reference to S5260:

1) assume that the unit of all input AGC is divided into 2 groups, often 2 units of group, 2 unit operation districts of every unit, The unit operation district of the 1st group of unit is 100MW to 200MW, 300MW to 400MW, and wherein the former is unit limit Operational Zone processed, the latter is unit optimized operation district, and the unit operation district of the 2nd group of unit is 50MW to 150MW, 300MW To 350MW, both unit optimized operation districts, then distributed combination mode includes:

{[(2,0),(2,0)],[(2,0),(1,1)],[(2,0),(0,2)],[(1,1),(2,0)],[(1,1),(1,1)],

[(1,1),(0,2)],[(0,2),(2,0)],[(0,2),(1,1)],[(0,2),(0,2)]}

2) the combined running district that every kind of distributed combination mode is corresponding is:

{(100,200)×2+(50,150)×2,(100,200)×2+(50,150)+(300,350),

(100,200)×2+(300,350)×2,(100,200)+(300,400)+(50,150)×2,

(100,200)+(300,400)+(50,150)+(300,350),(100,200)+(300,400)+(300,350)×2,

(300,400)×2+(50,150)×2,(300,400)×2+(50,150)+(300,350),

(300,400) × 2+ (300,350) × 2}={ (300,700), (550,900), (800,1100),

(500,900),(750,1100),(1000,1300),(700,1100),(950,1300),(1200,1500)}

3) union that cooperation district is all combined running districts of AGC unit, 300MW to 1500MW are put into.

4) according to S5260, the associating optimized operation district putting into AGC unit is 700MW to 1500MW.

5) from the cooperation district putting into AGC unit, deduct associating optimized operation district, obtain putting into the connection of AGC unit Close and limit Operational Zone, 300MW to 700MW.

6) bound putting into all combined running districts of AGC unit is ranked up:

{300,500,550,700,750,800,900,950,1000,1100,1200,1300,1500}

7) according to sequence after bound to put into AGC unit combine restriction Operational Zone split, draw multiple go out Power is interval:

{(300,500),(500,550),(550,700)}

8) the combined running district respectively exerted oneself under interval and various distributed combination mode is contrasted, draw input AGC machine Group is combining the available distributed combination mode that restriction Operational Zone is respectively exerted oneself under interval, is respectively as follows:

{[(2,0),(2,0)]},{[(2,0),(2,0)],[(1,1),(2,0)]},

{[(2,0),(2,0)],[(2,0),(1,1)],[(1,1),(2,0)]}

Described S5400) the associating optimized operation district of AGC unit will be put into plus all units not putting into AGC unit AGC gains merit apportioning cost, associating optimized operation district of entirely being stood:

$\mathrm{\&Omega;}+\underset{i=1}{\overset{n}{\&Sigma;}}(p_{AGC}^i\&times;{\mathrm{\&theta;}}^{\left|i-1\right|})$

Described S5500) the cooperation district putting into AGC unit is not put into AGC with combining restriction Operational Zone plus all The unit AGC of unit gains merit apportioning cost, obtains Quan Zhan cooperation district and complete station and combines restriction Operational Zone, computing formula ginseng Examine S5400.

6, calculate and put into the unit AGC of AGC unit and gain merit apportioning cost, including:

S6100) calculate full station AGC to gain merit apportioning cost, and judge to be in and put into the associating optimized operation district of AGC unit also It is to combine restriction Operational Zone；

S6200) when meeting specified conditions, AGC allocation flow is started；

S6300) if the AGC that entirely stands gains merit, apportioning cost in the associating optimized operation district putting into AGC unit, then starts optimal sorting Join algorithm；

S6400) if the AGC that entirely stands gains merit, apportioning cost limits Operational Zone in combining of input AGC unit, then start to limit and divide Join algorithm；

S6500) calculate unit AGC gain merit apportioning cost according to the unit AGC putting into each unit of AGC desired value of gaining merit, And it is issued to unit.

Further, described S6100), comprise the following steps:

1) calculate full station AGC to gain merit apportioning cost p_AGC:

$p_{AGC}=p_{set}-\underset{i=1}{\overset{n}{\&Sigma;}}(p_{AGC}^i\&times;{\mathrm{\&theta;}}^{\left|i-1\right|})$

2) with put into AGC unit combine optimized operation district, combine restriction Operational Zone compare, it is assumed that put into AGC The restriction Operational Zone of combining of unit is 300MW to 700MW, and associating optimized operation district is 700MW to 1500MW, Then when full station AGC gain merit apportioning cost at 300MW to 700MW between time, be in and combine restriction Operational Zone, stand AGC entirely When meritorious apportioning cost is between 700MW to 1500MW, it is in associating optimized operation district.

Described S6200) meeting following condition for the moment, startup AGC allocation flow:

S6210) gain merit apportioning cost summation and full station AGC of the unit AGC of all input AGC units gains merit apportioning cost Difference is more than AGC meritorious regulation dead band:

$|\underset{i=1}{\overset{n}{\&Sigma;}}(p_{AGC}^i\&times;{\mathrm{\&theta;}}^i)-p_{AGC}|>p_d$

S6220) put into AGC unit combination exert oneself model or cooperation district, combine restriction Operational Zone, associating optimum Operational Zone changes, with one below as trigger condition:

1) there is unit commitment AGC or have unit to exit AGC:

Wherein Δ θⁱFor unit former input AGC variable, when unit commitment AGC, Δ θ value is 1, otherwise Δ θ value is 0.

2) change of water level causes the input upper limit of exerting oneself of AGC unit, unit vibrating area, unit restriction Operational Zone to become Change.

Described S6300) if the AGC that entirely stands gains merit, apportioning cost in the associating optimized operation district putting into AGC unit, then starts Optimum allocation algorithm, comprises the following steps:

S6310) determine the target distribution combination of input AGC unit, comprise the following steps:

S6311) determine that can meet full station AGC gains merit all Optimal Distribution combination sides of input AGC unit of apportioning cost Formula:

$\left\{\begin{array}{c}\begin{array}{cc}f=f_i,& p_{AGC}\&Element;(\underset{\&OverBar;}{{\mathrm{\&Omega;}}_i},\stackrel{\&OverBar;}{{\mathrm{\&Omega;}}_i})\end{array}\\ i\&Element;(1,2,\mathrm{...},n)\end{array}\right.$

Wherein f is to represent that meeting full station AGC gains merit the row vector of all Optimal Distribution combinations of apportioning cost, with S5241 Described row vector Δ abc ... corresponding, when row vector f_iElement value when being 1, corresponding row vector Δ abc ... middle element Representative Optimal Distribution combination can be used.

S6312) full station AGC can be met gain merit the Optimal Distribution combination of apportioning cost with current by putting into AGC unit Distributed combination mode compare, select unit to pass through vibrating area platform time minimum Optimal Distribution combination and divide as target Cloth combination, is minimum and identical, the most entirely if there being multiple Optimal Distribution combination unit to pass through vibrating area platform Portion, as target distribution combination, comprises the following steps:

1) unit optimized operation district adjacent for input AGC unit and unit being limited Operational Zone merger is unit non-vibration district, By unit in unit optimized operation district, unit limits Operational Zone, the Optimal Distribution combination in unit operation district is converted to unit Optimal Distribution combination in unit non-vibration district, it is assumed that under current head unit unit limit Operational Zone be 0 to 120MW, 210MW to 420MW, unit vibrating area is 120MW to 210MW, and unit optimized operation district is 420MW To 700MW, then the unit non-vibration district of unit is 0 to 120MW, 210MW to 700MW, when unit gains merit When power is between 420MW to 700MW, unit is in the 3rd unit operation district and is in the 2nd unit non-vibration district.

2) according to the unit Optimal Distribution combination in unit non-vibration district, packet calculates unit and passes through the platform of vibrating area, Assume that certain group unit distribution mode in unit non-vibration district current is (x₁, x₂, x₃... x_m), wherein x_mFor being in The unit number of units in m non-vibration district, can meet full station AGC in jth kind and gain merit the Optimal Distribution combination of apportioning cost Under, this group unit distribution mode in unit non-vibration district isThen this group unit passes through vibrating area Platform time be:

$\left|\underset{i=1}{\overset{1}{\&Sigma;}}{x}_i^j-\underset{i=1}{\overset{1}{\&Sigma;}}{x}_i\right|+\left|\underset{i=1}{\overset{2}{\&Sigma;}}{x}_i^j-\underset{i=1}{\overset{2}{\&Sigma;}}{x}_i\right|+...+\left|\underset{i=1}{\overset{m}{\&Sigma;}}{x}_i^j-\underset{i=1}{\overset{m}{\&Sigma;}}{x}_i\right|$

3) by under jth kind Optimal Distribution combination, all packet units pass through the platform time of vibrating area and are added, and are input AGC unit passes through vibrating area platform time under jth kind Optimal Distribution combination；

4), the Optimal Distribution combination of apportioning cost of gaining merit from all satisfied full station AGC of input AGC unit, machine is selected Group passes through time minimum Optimal Distribution combination of vibrating area platform as target distribution combination, if there being multiple Optimal Distribution Combination unit passes through vibrating area platform and is minimum and identical, the most all as target distribution combination.

S6320) determine that the target putting into AGC unit is exerted oneself combination, comprise the following steps:

S6321) enumerate input AGC unit and can meet all combinations of exerting oneself of target distribution combination, it is assumed that The unit putting into AGC has 2 groups, and the 1st group has 3 units, and the 2nd group has 2 units, and target distribution combination is: { [(2,1), (1,1)] }, i.e. the 1st group unit has 2 to be in the 1st unit optimized operation district, and 1 is in the 2nd unit Excellent Operational Zone, the 2nd group of unit has 1 to be in the 1st unit optimized operation district, and 1 is in the 2nd unit optimized operation District, then combination of exerting oneself includes:

{[(1,1,2),(1,2)],[(1,2,1),(1,2)],[(2,1,1),(1,2)],[(1,1,2),(2,1)],[(1,2,1),(2,1)],[(2,1,1),(2,1)]}

Wherein first element in array [(1,1,2), (1,2)] represent, the 1st group the 1st of unit the, 2 units be in the 1st Individual unit optimized operation district, the 3rd unit is in the 2nd unit optimized operation district, at the 1st unit of the 2nd group of unit In the 1st unit optimized operation district, the 2nd unit is in the 2nd unit optimized operation district, it is also possible to combination side of exerting oneself Formula simplification is shown as:

{(1,1,2,1,2),(1,2,1,1,2),(2,1,1,1,2),(1,1,2,2,1),(1,2,1,2,1),(2,1,1,2,1)}

S6322) combination of exerting oneself that can meet target distribution combination is presently in putting into each unit of AGC Operational Zone compare, the combination of exerting oneself selecting unit to pass through vibrating area platform time minimum is exerted oneself combination side as target Formula, including:

1) unit optimized operation district adjacent for input AGC unit and unit being limited Operational Zone merger is unit non-vibration district, By input AGC unit in unit optimized operation district, unit limits Operational Zone, the combination of exerting oneself in unit operation district is changed For the combination of exerting oneself in unit non-vibration district, method is with reference to S6312；

2) according to putting into the AGC unit combination of exerting oneself in unit non-vibration district, the platform passing through vibrating area is calculated, Assuming there is m platform unit commitment AGC, the combination of exerting oneself in former unit non-vibration district is for (x₁, x₂, x₃... x_m), What jth kind can meet target distribution combination exerts oneself under combination, puts into AGC unit in unit non-vibration district Combination of exerting oneself isThen put into AGC unit and pass through the platform of vibrating area and time be:

$\underset{i=1}{\overset{m}{\&Sigma;}}\left|x_i^j-x_i\right|$

3) from all meet target distribution combination exert oneself combination, select to put into AGC unit and pass through vibration Combination that what district's platform was time minimum exert oneself is exerted oneself combination as target.

S6323) if having multiple combination of exerting oneself to pass through vibrating area platform and being minimum and identical, then input is weighted Compare again after the bad condition running priority level of AGC unit, it is assumed that total n combination of exerting oneself passes through vibrating area Platform is minimum and identical, then computing formula is:

$min\{\underset{i=1}{\overset{m}{\&Sigma;}}\&lsqb;|x_i^1-x_i|\&times;f\left({\mathrm{\&eta;}}_i\right)\&rsqb;,\&lsqb;\underset{i=1}{\overset{m}{\&Sigma;}}|x_i^2-x_i|\&times;f\left({\mathrm{\&eta;}}_i\right)\&rsqb;,...,\underset{i=1}{\overset{m}{\&Sigma;}}\&lsqb;|x_i^n-x_i|\&times;f\left({\mathrm{\&eta;}}_i\right)\&rsqb;\}$

Wherein η is unit bad condition running priority level.

S6330) under target exerts oneself combination, the unit operation district without unit changes, then use successively Allocation model, allocation flow schematic diagram is as it is shown in fig. 7, comprises following steps:

S6331) calculate gain merit apportioning cost and all input AGC unit unit AGC of full station AGC to gain merit apportioning cost summation Difference, as value to be allocated, it is assumed that have the n platform unit commitment AGC, value Δ p the most to be allocated to be:

$\mathrm{\&Delta;}p=p_{AGC}-\underset{i=1}{\overset{n}{\&Sigma;}}{p}_{AGC}^i$

S6332) if value to be allocated is more than 0, then calculate input AGC each unit unit AGC and gain merit apportioning cost with current The absolute difference of the unit optimized operation district upper limit, as unit can apportioning cost, and be arranged according to size, if value to be allocated is little In 0, then calculate that to put into AGC each unit unit AGC apportioning cost of gaining merit exhausted with the difference of current unit optimized operation district lower limit To value, as unit can apportioning cost, and be arranged according to size, it is assumed that have n platform unit commitment AGC, current optimized operation district Scope is respectivelyThen unit can apportioning cost pⁱ' it is:

$\left\{\begin{array}{cc}{p}^{i^{\&prime;}}=\left|\stackrel{\&OverBar;}{p_i}-p_{AGC}^i\right|,& \mathrm{\&Delta;}p>0\\ p^{i^{\&prime;}}=\left|\underset{\&OverBar;}{p_i}-p_{AGC}^i\right|,& \mathrm{\&Delta;}p<0\end{array}\right.$

S6333) select S6332 unit can the unit of apportioning cost sequence the 1st, take unit can apportioning cost, value to be allocated absolute Value, unit AGC have after the minimum of a value in distribution of work step-length, and travel direction correction as unit apportioning cost, value to be allocated Obtaining new value to be allocated after deduction unit apportioning cost, unit can obtain new list after apportioning cost deduction unit apportioning cost absolute value Machine can apportioning cost；

S6334) if value to be allocated is not equal to 0, then select S6332 unit can apportioning cost sequence the 2nd unit, take unit Can apportioning cost, value absolute value to be allocated, unit AGC make after having the minimum of a value in distribution of work step-length, and travel direction correction For unit apportioning cost, obtaining new value to be allocated after value to be allocated deduction unit apportioning cost, unit can divide by apportioning cost deduction unit Obtaining new unit after joining value absolute value can apportioning cost；

S6335) if value to be allocated is not equal to 0, then above procedure is continued, until value to be allocated is equal to 0；

S6336) if after all unit commitments, value to be allocated is not equal to 0, then repeat S6333, S6334, S6335, Until value to be allocated is equal to 0.

S6337) the unit apportioning cost putting into each unit of AGC is superimposed to unit AGC gain merit apportioning cost, obtains each The unit AGC of unit gains merit desired value

$p_{agc}^i={\mathrm{\&Delta;p}}^i+p_{AGC}^i$

Wherein Δ pⁱFor being assigned to the unit apportioning cost of unit i.

S6340) under target exerts oneself combination, if there being the unit operation district of unit to change, then equal proportion is used Allocation model, comprises the following steps:

S6341) calculate target to exert oneself the intermediate value in the combination optimized operation district putting into AGC unit under combination, it is assumed that have n Platform unit commitment AGC, target unit optimized operation district scope is respectivelyThen combine optimum The intermediate value of Operational Zone is:

$\underset{i=1}{\overset{n}{\&Sigma;}}(\underset{\&OverBar;}{p_n}\&divide;2+\stackrel{\&OverBar;}{p_n}\&divide;2)$

S6342) calculate full station AGC and gain merit apportioning cost and the difference combining optimized operation district intermediate value, as value Δ p to be allocated:

$\mathrm{\&Delta;}p=p_{AGC}-\underset{i=1}{\overset{n}{\&Sigma;}}(\underset{\&OverBar;}{p_n}\&divide;2+\stackrel{\&OverBar;}{p_n}\&divide;2)$

S6343) exert oneself according to target and put into the size of AGC each unit unit optimized operation district scope under combination, will Value proportional allocations to be allocated to each unit, obtains the unit apportioning cost Δ p of each unitⁱ:

${\mathrm{\&Delta;p}}^i=\mathrm{\&Delta;}p\&times;(\stackrel{\&OverBar;}{p_i}-\underset{\&OverBar;}{p_i})\&divide;\underset{i=1}{\overset{n}{\&Sigma;}}(\stackrel{\&OverBar;}{p_i}-\underset{\&OverBar;}{p_i})$

S6344) the unit apportioning cost putting into each unit of AGC is superimposed in each unit target unit optimized operation district Value, obtains the unit AGC of each unit and gains merit desired value

$p_{agc}^i={\mathrm{\&Delta;p}}^i+(\stackrel{\&OverBar;}{p_i}-\underset{\&OverBar;}{p_i})\&divide;2$

S6350) under target exerts oneself combination, if there being the unit operation district of unit to change, it is possible to use transregional It is sequentially allocated pattern, comprises the following steps:

1) Operational Zone that combination of target being exerted oneself is presently in the input each unit of AGC compares, for unit The unit that Operational Zone changes, using in target unit optimized operation district upper lower limit value near the limit value in former unit operation district as repairing Positive unit AGC gains merit apportioning cost, unit unchanged for unit operation district, apportioning cost that unit AGC is gained merit as Revise unit AGC to gain merit apportioning cost；

2) gaining merit apportioning cost and target unit optimized operation district based on putting into AGC each unit correction unit AGC, employing depends on Sub-distribution pattern, method and step are with reference to S6330, but need to use the former to replace input AGC each unit unit AGC Meritorious apportioning cost, uses the latter to replace current unit optimized operation district.

Described S6400) if the AGC that entirely stands gains merit, apportioning cost limits Operational Zone in combining of input AGC unit, then start Limit allocation algorithm, comprise the following steps:

S6410) determine the target distribution combination of input AGC unit, comprise the following steps:

S6411) the input AGC unit drawn according to S5360 combine limit Operational Zone in respectively exert oneself under interval available Distributed combination mode, determines that can meet full station AGC gains merit all distributed combination sides of input AGC unit of apportioning cost Formula, method and step are with reference to S6311；

All distributed combination modes of apportioning cost of S6412) can gaining merit with satisfied full station AGC from input AGC unit select Minimum unit is in unit and limits the combination of Operational Zone, as a example by S5360, when full station AGC apportioning cost of gaining merit is more than 500MW, less than 550MW time, distributed combination mode has { [(2,0), (2,0)], [(1,1), (2,0)] }, the wherein distribution of [(2,0), (2,0)] There are 2 units to be in unit under combination and limit Operational Zone, have at 1 unit under the distributed combination mode of [(1,1), (2,0)] Limit Operational Zone in unit, therefore select the latter as target distribution combination；

S6413) if the unit that input AGC unit has multiple distributed combination mode to be in unit restriction Operational Zone is Less and identical, then compare with current distributed combination mode, select unit to pass through the distribution group that vibrating area platform time is minimum Conjunction mode as target distribution combination, if having multiple distributed combination mode unit pass through vibrating area platform time be minimum also And identical, the most all as target distribution combination, method and step are with reference to S6312.

S6420) determine that the target putting into AGC unit is exerted oneself combination:

S6421) enumerate input AGC unit and can meet all combinations of exerting oneself of target distribution combination, method With step with reference to S6321；

S6422) can meet target distribution combination exert oneself combination to all, weighting puts into AGC unit After bad condition running priority level, the combination of exerting oneself selecting minimum unit to be in unit restriction Operational Zone is exerted oneself as target Combination:

1) by the combination conversion display of exerting oneself of unit, character γ in this unit operation district of serial number of unit operation district is replaced, When this unit operation district is unit restriction Operational Zone, γ value is 1, the γ value when this unit operation district is unit optimized operation district It is 0, it is assumed that have m platform unit, then the power combination after conversion is:

(γ₁,γ₂,γ₃,...,γ_m)

2), after weighting unit bad condition running priority level, select minimum unit to be in unit and limit the combination of exerting oneself of Operational Zone Mode is exerted oneself combination as target, it is assumed that total n kind combination of exerting oneself meets target distribution combination, then calculate Formula is:

$min\{\underset{i=1}{\overset{m}{\&Sigma;}}\&lsqb;{\mathrm{\&gamma;}}_i^1\&times;f\left({\mathrm{\&eta;}}_i\right)\&rsqb;,\underset{i=1}{\overset{m}{\&Sigma;}}\&lsqb;{\mathrm{\&gamma;}}_i^2\&times;f\left({\mathrm{\&eta;}}_i\right)\&rsqb;,...,\underset{i=1}{\overset{m}{\&Sigma;}}\&lsqb;{\mathrm{\&gamma;}}_i^n\&times;f\left({\mathrm{\&eta;}}_i\right)\&rsqb;\}$

Wherein η is unit bad condition running priority level.

S6423) if the target that S6422 draws is exerted oneself, combination, more than a kind, the most therefrom selects unit to pass through vibrating area The secondary minimum combination of exerting oneself of platform exerts oneself combination, method and step with reference to S6322 as target；

S6424) if the target that S6423 draws is exerted oneself, combination is more than a kind, then weighting unit bad condition runs excellent Comparing after first level, method and step are with reference to S6323 again.

S6430) under target exerts oneself combination, the unit operation district without unit changes, then use successively Allocation model, method and step are with reference to S6330.

S6440) under target exerts oneself combination, if there being the unit operation district of unit to change, equal proportion can be used Allocation model, method and step are with reference to S6340.

S6450) under target exerts oneself combination, if there being the unit operation district of unit to change, it is possible to use transregional Being sequentially allocated pattern, method and step are with reference to S6350.

Described S6500) according to put into each unit of AGC unit AGC gain merit desired value calculate unit AGC have the distribution of work Value, and it is issued to unit, comprise the following steps:

S6510) unit AGC active power adjustment value is calculated

${\mathrm{\&Delta;p}}_{agc}^i=\left|p_{agc}^i-p^i\right|;$

S6520) if unit AGC active power adjustment value has distribution of work step-length less than or equal to unit AGC, the most directly will be single Machine AGC desired value of gaining merit is gained merit apportioning cost as unit AGC, and is issued to unit:

$\begin{array}{cc}if& {\mathrm{\&Delta;p}}_{agc}^i\&le;{\mathrm{\&Delta;p}}_{sp}^i\end{array}$

$p_{AGC}^i=p_{agc}^i$

S6530) if unit AGC active power adjustment value has distribution of work step-length more than unit AGC, then substep is by unit AGC Meritorious target value conversion is that unit AGC gains merit apportioning cost, and is issued to unit, flow process as shown in Figure 8, including following step Rapid:

S6531) desired value real the send out value meritorious more than unit if unit AGC gains merit, then add meritorious for unit real value of sending out Unit AGC has distribution of work step-length, is issued to unit as unit AGC apportioning cost of gaining merit, otherwise then by meritorious for unit real Send out value to deduct unit AGC and have distribution of work step-length, be issued to unit and gain merit apportioning cost as unit AGC:

$\left\{\begin{array}{cc}{p}_{AGC}^i=p^i+{\mathrm{\&Delta;p}}_{sp}^i,& p_{agc}^i>p^i\\ p_{AGC}^i=p^i-{\mathrm{\&Delta;p}}_{sp}^i,& p_{agc}^i<p^i\end{array}\right.$

S6532) set time t is waited₆:

t₃=t₆×n n∈N₊

S6533) unit AGC active power adjustment value is recalculated:

${\mathrm{\&Delta;p}}_{agc}^i=\left|p_{agc}^i-p^i\right|$

S6534) if unit AGC active power adjustment value has distribution of work step-length less than or equal to unit AGC, the most directly will be single Machine AGC desired value of gaining merit is gained merit apportioning cost as unit AGC, and is issued to unit, otherwise repeat S6531, S6532, S6533, S6534.

Described S7000) according to put into each unit of AGC unit AGC gain merit apportioning cost and and unit PID function throw Move back state, AGC throws and moves back state etc., determines that unit gains merit setting value, method and step with reference to Fig. 6.

Described S8000) meritorious for each the unit unit frequency modulation correction unit that is superimposed to is gained merit setting value, show that unit is meritorious and set Periodical repair on the occasion of

${\mathrm{\&Delta;p}}_{set}^i=p_{set}^i+p_f^i$

Described S9000) each unit PID perform equipment and set correction value according to unit is meritorious, carries out PID closed loop regulation, Make the meritorious real value of sending out of unit trend towards unit meritorious setting correction value, and be maintained in dead zone range, when shaping modes is speed governing During device aperture pattern, it is computer supervisory control system machine set PLC that PID performs equipment；When shaping modes is speed regulator power mode Time, it is speed regulator that PID performs equipment.

The general principle of the present invention and principal character and advantages of the present invention have more than been shown and described.The technology people of the industry Member is it should be appreciated that the present invention is not restricted to the described embodiments, and described in above-described embodiment and specification is the present invention's Principle, without departing from the spirit and scope of the present invention, the present invention also has various changes and modifications, these change and Improve the most within the scope of the claimed invention.

Claims (17)

1. a power station automatic generation is gained merit Force control system, it is characterized in that, including AGC functional block, pid function block and preprocessing function block, device parameter, dynamic parameter directly input or are sent to AGC functional block after the process of preprocessed functional block, pid function block updates；

AGC functional block, pid function block receive the pretreatment values that preprocessing function block sends respectively；

AGC functional block is set up the combination of hydropower station unit according to received data and is exerted oneself model, and calculates cooperation region, associating optimized operation region and combine restriction operation area；In associating optimized operation region or combining restriction operation area according to full station AGC apportioning cost of gaining merit, starting corresponding AGC allocation flow and obtain each the unit AGC putting into AGC unit and gain merit apportioning cost be issued to unit；

Put into each generating set of AGC and throw state of moving back and AGC throws and moves back state further according to gain merit apportioning cost and unit PID function of unit AGC, determine that unit is gained merit setting value；Each generating set unit frequency modulation correction be meritorious to be superimposed to unit and gains merit setting value, show that unit is meritorious and sets correction value and be sent to pid function block；

Pid function block carries out PID closed loop regulation according to the meritorious correction value that sets of unit to each unit.

2. power station as claimed in claim 1 automatic generation is gained merit Force control system, it is characterized in that, described AGC functional block receives the pretreatment values including that setting value preparation block gained merit by head value preparation block, the meritorious desired value preparation block in full station, unit and unit bad condition running priority level preparation block sends；Pid function block receives and includes that unit is gained merit the meritorious real pretreatment values sending out value preparation block of setting value preparation block and unit；Pid function block sends to AGC functional block and includes that unit meritorious real sends out value, unit frequency modulation correction has work value, unit to gain merit regulation dead band and unit is gained merit and set the parameter value that difference limits.

3. power station as claimed in claim 1 automatic generation is gained merit Force control system, it is characterised in that the combination that described AGC functional block is set up model of exerting oneself includes:

The optimum combination putting into AGC generating set is exerted oneself model, utilize this model according to putting into AGC generating set combination optimized operation district under various Optimal Distribution combinations, determine and put into the available Optimal Distribution combination that AGC generating set is respectively exerted oneself under interval in associating optimized operation district；

Model of exerting oneself is combined in the restriction putting into AGC generating set, utilize this model according to putting into AGC generating set combined running district under various distributed combination modes, determine and put into the available distributed combination mode that AGC generating set is respectively exerted oneself under interval within combining restriction Operational Zone.

4. power station as claimed in claim 1 automatic generation is gained merit Force control system, it is characterized in that, the AGC allocation flow of described AGC functional block is: if the AGC that entirely stands gains merit, apportioning cost is in the associating optimized operation district putting into AGC generating set, then start optimum allocation algorithm；If full station AGC gains merit, apportioning cost limits Operational Zone in combining of input AGC generating set, then start restriction allocation algorithm；Target value conversion of being gained merit by the unit AGC obtained under AGC allocation flow is that unit AGC gains merit apportioning cost, and is issued to generating set.

5. a power station automatic generation is gained merit force control method, it is characterised in that include following operation:

1) Core Feature of power station automatic generation real power control is divided in two functional blocks of AGC, PID processes, be separated to non-core function in plug-in functional block process；

2) AGC, pid function block receive the pretreatment values that device parameter, dynamic parameter and pretreatment module send respectively；AGC, pid function block calculate according to dynamic parameter and device parameter and obtain intermediate parameters；

3) AGC functional block is set up the combination of hydropower station unit according to received data and is exerted oneself model calculate cooperation district, associating optimized operation district, combine restriction Operational Zone；In associating optimized operation region or combining restriction operation area according to full station AGC apportioning cost of gaining merit, starting corresponding AGC allocation flow and obtain each the unit AGC putting into AGC unit and gain merit apportioning cost be issued to unit；

Put into each generating set of AGC according to gain merit apportioning cost and unit PID function of the unit AGC assigned to throw state of moving back and AGC throws and moves back state, determine that unit is gained merit setting value；Again meritorious for each the generating set unit frequency modulation correction unit that is superimposed to is gained merit setting value, show that unit is meritorious and set correction value and be sent to pid function block；

4) each generating set PID performs equipment and sets correction value according to unit is meritorious, carries out PID closed loop regulation.

6. power station as claimed in claim 5 automatic generation is gained merit force control method, it is characterized in that, the pretreatment values that described AGC is received includes that the pretreatment values of head value, the meritorious pretreatment values of desired value in full station, unit are gained merit the pretreatment values of setting value and the pretreatment values of unit bad condition running priority level；The pretreatment values that PID is received includes that unit is gained merit the pretreatment values of setting value and the meritorious real pretreatment values sending out value of unit.

7. power station as claimed in claim 5 automatic generation is gained merit force control method, it is characterized in that, described device parameter is the parameter characteristic of each generating set in power station, and the device parameter that AGC functional block receives includes: under different heads, under the unit output upper limit, different head, under unit each unit vibrating area scope, different head, each unit of unit limits Operational Zone scope；The device parameter that pid function block receives includes: unit meritorious regulation dead band, unit meritorious setting difference limit, unit meritorious sudden change boundary, coefficient of frequency modulation；

Described dynamic parameter is equipment real-time status or operation conditions real-time parameter；

The calculating of described intermediate parameters includes: AGC meritorious regulation dead band, AGC meritorious setting difference limit, AGC meritorious sudden change boundary, unit AGC have distribution of work step-length, unit frequency modulation correction meritoriousWith unit each unit operation district under different heads.

8. power station as claimed in claim 6 automatic generation is gained merit force control method, it is characterised in that the acquisition of described pretreatment values includes:

Being retrieved as of the meritorious real pretreatment values sending out value of a, unit: the priority in active power measurement sources multiple to unit is ranked up, takes the highest measured value of the measured value quality in active power measurement source priority preferably as meritorious real the value of unit；

B, being retrieved as the pretreatment values of head value: obtain final head value according to head setting means and head value situation of change；

C, being retrieved as of pretreatment values of desired value meritorious to full station: put into lower different situations according to AGC function and combine the meritorious desired value in full station, electric network active setting value, the meritorious setting value of collection control and power plant setting value of gaining merit and obtain；

D, unit are gained merit being retrieved as of pretreatment values of setting value: gain merit setting value and unit PID function throws state of moving back and AGC throws and moves back state according to the unit of input, determine that unit is gained merit setting value；

F, being retrieved as of pretreatment values of unit bad condition running priority level: 1) each unit be in after turn(a)round last time limit operation time of Operational Zone and vibrating area and be weighted statistics, wherein vibrating area runs the statistical weight of time far above the statistical weight limiting Operational Zone and running the time；2) time after each unit weighted statistical is ranked up, the most from short to long, sets automatic priority the most from high to low；3) to each group setup manual priority level；4) when priority set-up mode is automatic, take from dynamic priority as unit bad condition running priority level, when priority set-up mode is manual, take manual priority level as unit bad condition running priority level.

9. power station as claimed in claim 6 automatic generation is gained merit force control method, it is characterised in that the combination that described AGC functional block is set up model of exerting oneself includes:

The optimum combination putting into AGC generating set is exerted oneself model, utilize this model according to putting into AGC generating set combination optimized operation district under various Optimal Distribution combinations, determine and put into the available Optimal Distribution combination that AGC generating set is respectively exerted oneself under interval in associating optimized operation district；

Model of exerting oneself is combined in the restriction putting into AGC generating set, utilize this model according to putting into AGC generating set combined running district under various distributed combination modes, determine and put into the available distributed combination mode that AGC generating set is respectively exerted oneself under interval within combining restriction Operational Zone.

10. power station as claimed in claim 9 automatic generation is gained merit force control method, it is characterised in that being retrieved as of the available Optimal Distribution combination that described inputs AGC generating set is respectively exerted oneself under interval within combining optimized operation district:

1) limit Operational Zone scope according to the upper limit of exerting oneself under each head of unit, each unit vibrating area scope, each unit, the generating set putting into AGC is grouped；

2) under current head, limit Operational Zone scope, the upper limit of exerting oneself according to unit unit vibrating area scope, unit, calculate the unit optimized operation district of each group of generating set；

3) for each group of generating set, exert oneself in the distribution situation in each unit optimized operation district according to generating set, calculate each group of generating set combination optimized operation district under various Optimal Distribution modes；

4) for the generating set of all input AGC, it is in each combination optimized operation district organizing generating set of the different distributions mode in unit optimized operation district and correspondence according to each group of generating set, calculates the combination optimized operation district of all input AGC generating sets under various Optimal Distribution combination；

5) according to putting into the AGC unit union combining optimized operation district under all Optimal Distribution combinations, the associating optimized operation district putting into AGC generating set is calculated；

6) according to putting into AGC generating set combination optimized operation district under various Optimal Distribution combinations, determine and put into the available Optimal Distribution combination that AGC generating set is respectively exerted oneself under interval in associating optimized operation district.

11. power station as claimed in claim 9 automatic generations are gained merit force control method, it is characterised in that described input AGC generating set combine limit Operational Zone in being retrieved as of available distributed combination mode respectively exerted oneself under interval:

1) limit Operational Zone scope according to the upper limit of exerting oneself under each head of unit, each unit vibrating area scope, each unit, the generating set putting into AGC is grouped；

2) under current head, according to unit unit vibrating area, the upper limit of exerting oneself, the unit operation district of each group of generating set is calculated；

3) for each group of generating set, exert oneself in the distribution situation in each unit operation district according to generating set, calculate each group of generating set combined running district under various distribution modes；

4) for the generating set of all input AGC, it is in each combined running district organizing generating set of the different distributions mode in unit operation district and correspondence according to each group of generating set, calculates the combined running district of all input AGC generating sets under various distributed combination mode；

5) calculate the cooperation district putting into AGC generating set and combine restriction Operational Zone；

6) according to putting into AGC generating set combined running district under various distributed combination modes, determine and put into the available distributed combination mode that AGC generating set is respectively exerted oneself under interval within combining restriction Operational Zone.

12. power station as claimed in claim 6 automatic generations are gained merit force control method, it is characterised in that the unit AGC of described each input AGC generating set gains merit being retrieved as of apportioning cost:

1) calculate full station AGC to gain merit apportioning cost, and judge to be in and put into the associating optimized operation district of AGC generating set and still combine restriction Operational Zone；

2) AGC allocation flow is started when meeting following condition: gain merit apportioning cost summation and full station AGC of the unit AGC of all input AGC units gains merit the difference of apportioning cost more than AGC meritorious regulation dead band；Or the combination putting into AGC unit exert oneself model or cooperation district, associating optimized operation district, combine restriction Operational Zone and change；

3) when carrying out AGC distribution, if the AGC that entirely stands gains merit, apportioning cost is in the associating optimized operation district putting into AGC generating set, then start optimum allocation algorithm；If full station AGC gains merit, apportioning cost limits Operational Zone in combining of input AGC generating set, then start restriction allocation algorithm；

4) calculate and put into the unit AGC of each generating set of AGC and gain merit apportioning cost be issued to generating set.

13. power station as claimed in claim 6 automatic generations are gained merit force control method, it is characterised in that described AGC distribution, when starting optimum allocation algorithm, determines and puts into the target distribution combination of AGC generating set and target is exerted oneself combination；Under target exerts oneself combination, if not having the unit operation district of generating set to change, using and being sequentially allocated pattern；If the unit operation district having generating set changes, then use proportional allocations pattern or use and transregional be sequentially allocated pattern；

Described AGC distribution start limit allocation algorithm time, determine and put into the target distribution combination of AGC generating set and target is exerted oneself combination；Under target exerts oneself combination, if not having the unit operation district of generating set to change, being sequentially allocated pattern have the unit operation district of generating set to change if then using, proportional allocations pattern can be used or use and transregional be sequentially allocated pattern；

Described unit AGC gains merit being retrieved as of apportioning cost: calculate unit AGC active power adjustment value；If unit AGC active power adjustment value has distribution of work step-length less than or equal to unit AGC, desired value of directly being gained merit by unit AGC is gained merit apportioning cost as unit AGC, and is issued to generating set；Otherwise, desired value of being gained merit by unit AGC substep is converted into unit AGC and gains merit apportioning cost, and is issued to generating set.

14. power station according to claim 6 automatic generations are gained merit force control method, it is characterised in that described unit is gained merit being retrieved as of setting value:

1) when unit PID function exits, take the meritorious real value of sending out of unit unit and gain merit setting value as unit；

2) when unit PID function puts into and unit does not puts into AGC, take unit setting value of gaining merit to gain merit apportioning cost as unit AGC, if manually setting value real the send out difference of value meritorious with unit in the range of unit is gained merit and set difference limit and less than or equal to the unit output upper limit, then take operations staff manually to set value and gain merit setting value as unit, otherwise then keep unit setting value of gaining merit constant；

3) as unit PID function input and unit commitment AGC, if unit AGC gains merit, apportioning cost real the send out difference of value meritorious with unit is in the range of unit is gained merit and set difference limit and less than or equal to the unit output upper limit, then take unit AGC apportioning cost of gaining merit to gain merit setting value as unit, otherwise then keep unit setting value of gaining merit constant.

15. power station according to claim 6 automatic generations are gained merit force control method, it is characterized in that, described PID closed loop regulation, the meritorious real value of sending out of unit is made to trend towards unit meritorious setting correction value, and be maintained in dead zone range, when shaping modes is speed regulator aperture pattern, it is computer supervisory control system machine set PLC that PID performs equipment；When shaping modes is speed regulator power mode, it is speed regulator that PID performs equipment.

16. power station according to claim 6 automatic generations are gained merit force control method, it is characterised in that be additionally provided with the abnormal exit strategy of each functional block in described AGC functional block, pid function block, including:

100) abnormality causing unit PID function to exit is determined；

200) abnormality causing power station AGC function to exit is determined；

300) abnormality causing unit to exit AGC is determined.

17. power station according to claim 16 automatic generations are gained merit force control method, it is characterised in that the abnormality that described determination causes unit PID function to exit includes:

110) the unit meritorious superposition value real difference of sending out value meritorious with unit of setting value and unit frequency modulation correction of gaining merit regulates dead band beyond unit is meritorious, and continues for some time t₁；

120) unit gains merit setting value real the send out difference of value meritorious with unit beyond unit meritorious sudden change boundary；

130) there is saltus step in the meritorious real value of sending out of unit, and hop value is beyond unit meritorious sudden change boundary；

140) unit is in non-power generating state；

150) unit all active power measurements source measured value quality is bad；

160) unit has accident；

Described determine that the abnormality causing power station AGC function to exit includes:

210) AGC performs equipment and performs equipment communication abnormality or interruption with certain unit PID；

220) the gain merit superposition value real difference of sending out value summation meritorious with all unit units of summation of the meritorious desired value in full station and all unit frequency modulation corrections regulates dead band beyond AGC is meritorious, and continues for some time t₂, t₂>=t₁；

230) the meritorious desired value in full station real the send out difference of value summation meritorious with all unit units is beyond AGC meritorious sudden change boundary；

240) meritorious real the sending out of all unit units is worth summation generation saltus step, and hop value is beyond AGC meritorious sudden change boundary；

250) the meritorious desired value in full station is not in Quan Zhan cooperation district；

260) certain unit all active power measurements source measured value quality is bad；

270) there is accident in full factory；

Described determine that the abnormality causing unit to exit AGC includes:

310) AGC function in power station is in the state of exiting；

320) unit PID function is in the state of exiting；

330) unit setting value of gaining merit is not equal to unit AGC and gains merit apportioning cost, and continues for some time t₃；

340) the meritorious real value of sending out of unit is in unit vibrating area and exceedes certain time t₄；

350) the meritorious real value of sending out of unit is in unit and limits Operational Zone and exceed certain time t₅, t₅>>t₄。