CN109494765A - Alternating current-direct current combined hybrid system transient stability control method based on EEAC - Google Patents
Alternating current-direct current combined hybrid system transient stability control method based on EEAC Download PDFInfo
- Publication number
- CN109494765A CN109494765A CN201811346098.0A CN201811346098A CN109494765A CN 109494765 A CN109494765 A CN 109494765A CN 201811346098 A CN201811346098 A CN 201811346098A CN 109494765 A CN109494765 A CN 109494765A
- Authority
- CN
- China
- Prior art keywords
- machine
- generator
- power
- rotor angle
- generator rotor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
The alternating current-direct current combined hybrid system transient stability control method based on EEAC that the present invention provides a kind of, comprising: the generator rotor angle and power data of each generator after breaking down are received by control centre；Machine group identification is carried out to each generator；According to generator rotor angle data and power data, machine group after identification is equivalent to a double machine set system using corresponding local inertia center, then double machine set systems are equivalent to an one machine infinity bus system, and are the generator rotor angle relationship matched curve after failure by one machine infinity bus system curve matching；Transient stability differentiation is carried out using extended equal area criterion strategy.This method can after being latched failure the effectively amount of unbalance between balancing machine power and electromagnetic power, cut machine capacity after reducing failure, have it is biggish can practicality, reliable guarantee can be provided for safe and stable operation of power system.
Description
Technical field
The present invention relates to power grid power supply technique field more particularly to a kind of alternating current-direct current combined hybrid system transient state based on EEAC are steady
Qualitative Control.
Background technique
With the development of power grid power supply technique, China's extra-high voltage operation scale will just have reached " six friendships five in the end of the year 2016
Directly ", major regional power grid has passed through the realization interconnection of extra-high voltage alternating current-direct current transmission line of electricity, and alternating current-direct current mixed connection has become China's power grid
Characteristic feature.With the reduction of electrical distance between the increase and direct current drop point of direct current scale, when receiving end fault in ac transmission system,
It may result in multiple-circuit line simultaneously or sequentially commutation failure and be easy to cause dc power to transmit when multiple-circuit line commutation failure
It interrupts, the serious safe and stable operation for threatening direct current system.This interconnected power system operation variation is big, and alternating current-direct current hybrid connected structure is multiple
Miscellaneous and existing some weak links bring a series of choose to the stability analysis, control and operation of AC and DC power system
War.
In order to realize that the emergent control to the significantly more efficient protection of bulk power system, for alternating current-direct current mixing bulk power grid includes
When DC line is latched, the prior art of use includes: to be supported by the DC line of other branches blocked line power；
Or the method using HVDC Modulation, such as frequency modulation(PFM), the modulation of direct current active power in case of emergency carry out power support, mention
The transient stability of high direct current both ends AC system；Or under large disturbances by cut machine control change cut node inertia,
Power distribution in kinetic energy and reconfiguration system maintains electrical power system transient power-angle stability, such as pre- based on section contact heat input
The transient stability of survey promptly cuts machine strategy and cuts machine strategy etc. based on control moment angular speed deviation.The peace of the above prior art
Full Stable Control Strategy can guarantee the transient stability of electric system to a certain extent, but the above-mentioned prior art is all from list
Aspect is analyzed and researched, less to guarantee power system transient stability progress quantitative analysis, seldom considers to occur sternly when system
, there is relatively high power vacancy in weight failure, and scheduling normal DC route has been unable to meet power to the support of faulty line power and lacks
When volume, one-side control strategy cannot maintain the transient stability of electric system.Need to then end cutting machine, HVDC Modulation sent at this time
Ac and dc systems transient stability is improved etc. the emergency control policy of more resource coordinatings.
On the other hand, controlled at present for direct-current emergency and the coordination control strategy of power generation machine-cut machine is based on time domain mostly
Emulation carries out qualitative analysis, the deficiencies of it is low that there are computational efficiencies, and time-consuming, can not theoretically obtain direct-current emergency control amount and
Generator cuts machine amount, do not solved the problems, such as from mechanism between more resources how coordinated control.
Summary of the invention
Of the invention provides a kind of alternating current-direct current combined hybrid system transient stability control method based on EEAC, is handed over improving
The stability of direct current hybrid system.
To achieve the goals above, this invention takes following technical solutions.
The alternating current-direct current combined hybrid system transient stability control method based on EEAC that the present invention provides a kind of, feature exist
In, comprising:
The generator rotor angle and power data of each generator after breaking down are received by control centre；
According to the generator rotor angle data, machine group identification is carried out to each generator；
According to the generator rotor angle data and power data by the machine group after identification using corresponding local inertia center etc.
Effect is a double machine set system, double machine set systems is then equivalent to an one machine infinity bus system, and will be described
One machine infinity bus system curve matching is the generator rotor angle relationship matched curve after failure；
According to the power angular dependence matched curve after the failure, carried out using extended equal area criterion EEAC strategy temporary
State stability distinguishing.
Further, carrying out transient stability differentiation using extended equal area criterion EEAC strategy includes: when acceleration area
When domain Δ ac > do not take the retardation area Δ RED under any control strategy, system unstability takes direct-current emergency control strategy, meter
It calculates by the direct-current emergency control increased retardation area Δ DC of institute otherwise, as Δ ac < Δ RED, system keeps transient stability；
As Δ ac > Δ RED+ Δ DC, system continues unstability, calculates by cutting the machine control increased retardation area Δ gt of institute, and take and cut
The control strategy of machine.
Further, according to the generator rotor angle data, machine group identification is carried out to each generator, comprising: use
EEAC is arranged from small to large according to each generator's power and angle, is differentiated by maximal clearance, and each generator is divided into critical machine
Group and remaining machine group.
Further, the criticality benchmark group and remaining machine group are adopted according to the generator rotor angle data and power data
A double machine set system is equivalent to corresponding local inertia center, comprising:
It is equivalent to the criticality benchmark group to lose the machine set system of transient stability according to following formula (1):
The remaining machine group is equivalent to according to following formula (2) to keep the machine set system of transient stability:
Wherein, MsIt is the equivalent inertia time constant of critical machine, MiIt is the inertial time of every generator in criticality benchmark group
Between constant, δsIt is the equivalent generator rotor angle of criticality benchmark group, δiIt is the generator rotor angle of every generator in critical machine, MAIt is remaining machine group
Equivalent inertia time constant, MjIt is the inertia time constant of every generator in a remaining group of planes, δAIt is the equivalent of remaining machine group
Generator rotor angle, δjIt is the generator rotor angle that remaining machine group concentrates every generator.
Further, two machine set systems are equivalent to an one machine infinity bus system, comprising: according to following public affairs
Two machine set systems are equivalent to an one machine infinity bus system by formula (3):
Wherein, M is the inertia time constant of one machine infinity bus system, and δ is the generator rotor angle of one machine infinity bus system, PmIt is single machine
The mechanical output of Infinite bus system, PmiIt is the mechanical output of every generator in criticality benchmark group, PmjIt is every in remaining machine group
The mechanical output of platform generator.
It further, is the power angular dependence matched curve after failure by the one machine infinity bus system curve matching,
It specifically includes: according to the electromagnetic power of the one machine infinity bus system and generator, obtaining in image one machine infinity bus system
Shown in power angular dependence matched curve such as following formula (4):
Wherein, p, b, c are unknown parameter, PmIt is the mechanical output of one machine infinity bus system, Pe is the electromagnetic work of generator
Rate,For in equation of rotor motion, the second dervative of generator rotor angle δ.
Further, as acceleration region area Δ ac > do not take the retardation area Δ RED under any control strategy,
System unstability takes direct-current emergency control strategy, calculates by the direct-current emergency control increased retardation area Δ DC of institute otherwise, works as Δ
When ac < Δ RED, system keeps transient stability；As Δ ac > Δ RED+ Δ DC, system continues unstability, calculates by cutting machine control
The increased retardation area Δ gt of system institute, and take the control strategy for the machine of cutting, comprising:
The acceleration region area Δ ac is calculated according to following formula (5):
The retardation area that do not take under any control strategy is calculated according to following formula (6):
It is calculated according to following formula (7) by the direct-current emergency control increased retardation area Δ DC of institute:
It is calculated according to following formula (8) by cutting the machine control increased retardation area Δ gt of institute:
Wherein,
δ0Generator rotor angle when occurring for failure, δcGenerator rotor angle when for fault clearance, δqTo execute function when emergency control policy
Angle, δ 'hTo implement the critical clearing angle after direct-current emergency power control, δ "hIt is that the limit after further implementing to cut machine is taken to cut
Except angle, Δ PG is to cut machine amount, PIFor electromagnetic power under normal circumstances, PIIIFor the electromagnetic power after fault clearance.
The technical side provided by the alternating current-direct current combined hybrid system transient stability control method based on EEAC of aforementioned present invention
Case can be seen that the present invention technical characterstic big using the distinctive quick response of direct current transportation itself and pondage, in EEAC
On the basis of derived direct current active power modulation principle, and urgent dc power control is proposed based on EEAC principle and cuts machine phase
In conjunction with alternating current-direct current Transient Stability Control strategy, existing cut is aided with emergency DC power raising in machine measure in systems, can
With effectively the amount of unbalance between balancing machine power and electromagnetic power, the machine of cutting after reducing failure hold after being latched failure
Amount.The method of the present invention in alternating current-direct current mixed connection catenet have it is biggish can practicality, can be applied to on-line operation, adapt to
The development trend of the following power grid can provide reliable guarantee for safe and stable operation of power system.
The additional aspect of the present invention and advantage will be set forth in part in the description, these will become from the following description
Obviously, or practice through the invention is recognized.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill of field, without creative efforts, it can also be obtained according to these attached drawings others
Attached drawing.
Fig. 1 is a kind of alternating current-direct current combined hybrid system transient stability control method based on EEAC of the embodiment of the present invention 1
Flow chart；
Fig. 2 is a kind of alternating current-direct current combined hybrid system transient stability control method based on EEAC of the embodiment of the present invention 1
Flow diagram；
Fig. 3 is that direct-current emergency power control and machine of cutting control the schematic diagram matched in the embodiment of the present invention 1；
Fig. 4 is that the three machine systems of use in the embodiment of the present invention 2 connect infinite busbar topological structure schematic diagram；
Fig. 5 is the dc power emergency lifting of the embodiment of the present invention 2 and the simulation result diagram of emergency descent；
Fig. 6 is criticality benchmark group and generator's power and angle curve graph after the failure of remaining machine group in the embodiment of the present invention 2；
Fig. 7 is generator electromagnetic power curve graph after the failure of the embodiment of the present invention 2；
Fig. 8 be embodiment 2 it is equivalent after one machine infinity bus system power angle characteristics curve graph；
Fig. 9 is the simulation result diagram of the embodiment of the present invention 2.
Specific embodiment
Embodiments of the present invention are described below in detail, the example of the embodiment is shown in the accompanying drawings, wherein from beginning
Same or similar element or element with the same or similar functions are indicated to same or similar label eventually.Below by ginseng
The embodiment for examining attached drawing description is exemplary, and for explaining only the invention, and is not construed as limiting the claims.
Those skilled in the art of the present technique are appreciated that unless expressly stated, singular " one " used herein, " one
It is a ", " described " and "the" may also comprise plural form.It is to be further understood that being arranged used in specification of the invention
Diction " comprising " refers to that there are the feature, integer, step, operations, but it is not excluded that in the presence of or addition it is one or more other
Feature, integer, step, operation.
Those skilled in the art of the present technique are appreciated that unless otherwise defined, all terms used herein (including technology art
Language and scientific term) there is meaning identical with the general understanding of those of ordinary skill in fields of the present invention.Should also
Understand, those terms such as defined in the general dictionary, which should be understood that, to be had and the meaning in the context of the prior art
The consistent meaning of justice, and unless defined as here, it will not be explained in an idealized or overly formal meaning.
In order to facilitate understanding of embodiments of the present invention, it is further explained below in conjunction with attached drawing with specific embodiment
It is bright, and do not constitute the restriction to the embodiment of the present invention.
Embodiment 1
Fig. 1 is the flow chart of the alternating current-direct current combined hybrid system transient stability control method of the invention based on EEAC, and Fig. 2 is
The flow diagram of alternating current-direct current combined hybrid system transient stability control method based on EEAC of the invention, referring to Figures 1 and 2,
This method comprises:
S1 receives the generator rotor angle and power data of each generator after breaking down by control centre.
S2 carries out machine group identification according to generator rotor angle data, to each generator.
Preferably, it is arranged from small to large using extended equal area criterion EEAC according to each generator's power and angle, by between maximum
Gap differentiates, each generator is divided into criticality benchmark group S and remaining machine group A.
EEAC is to be differentiated using the gap of generator's power and angle δ to criticality benchmark group.Specific method is, when being
After system receives impact, the generator rotor angle of each generator after breaking down by arranging from small to large, between front and back in this arrangement
Generator is divided into criticality benchmark S groups and remaining unit A groups by the maximum place of gap.
Schematically, it is assumed that system collects one group of generator rotor angle data after by failure:
δk=[δk1, δk2... δkn] k=1,2,3 ..., n
K-th of generating set presses the generator rotor angle sampled value of Δ t time interval after a failure in system, and after a failure certain is for the moment
It carves, the generator rotor angle of generating set each in system is ranked up, as a result as follows:
δ2<δ3<δ1<δ6<δ4<δ5
The interval sorted above is ranked up again, result is as follows:
(δ3-δ2)<(δ6-δ1)<(δ5-δ4)<(δ4-δ6)<(δ1-δ3)
Then it is the largest interval of this group of generator's power and angle between generator 3 and 1, is sorted according to generator rotor angle size, pass through the 3rd
Generating set is divided into the criticality benchmark being made of generator 1,6,4,5 by a and the 1st generator, and is made of generator 2,3
Remaining unit.
S3 is according to the generator rotor angle data and power data by the machine group after identification using corresponding local inertia center
It is equivalent to a double machine set system, double machine set systems are then equivalent to an one machine infinity bus system, and will be described
One machine infinity bus system curve matching be failure after power angular dependence matched curve.
Further, the criticality benchmark group and remaining machine group are adopted according to the generator rotor angle data and power data
A double machine set system is equivalent to corresponding local inertia center, comprising:
It is equivalent to the criticality benchmark group to lose the machine set system of transient stability according to following formula (1):
The remaining machine group is equivalent to according to following formula (2) to keep the machine set system of transient stability:
Wherein, MsIt is the equivalent inertia time constant of critical machine, Mi is the inertia of every generator in criticality benchmark group
Time constant, δsIt is the equivalent generator rotor angle of criticality benchmark group, δ i is the generator rotor angle of every generator in critical machine, MAIt is remaining unit
The equivalent inertia time constant of group, MjIt is the inertia time constant of every generator in a remaining group of planes, δABe remaining machine group etc.
Imitate generator rotor angle, δjIt is the generator rotor angle that remaining machine group concentrates every generator.
Further, two machine set systems are equivalent to an one machine infinity bus system, comprising: according to following public affairs
Two machine set systems are equivalent to an one machine infinity bus system by formula (3):
Wherein, M is the inertia time constant of one machine infinity bus system, and δ is the generator rotor angle of one machine infinity bus system, PmIt is single machine
The mechanical output of Infinite bus system, PmiIt is the mechanical output of every generator in criticality benchmark group, PmjIt is every in remaining machine group
The mechanical output of platform generator.
It further, is the power angular dependence matched curve after failure by the one machine infinity bus system curve matching,
Then it is solved using least square method, specifically includes: according to the electromagnetic power of the one machine infinity bus system and generator, obtaining
Into image one machine infinity bus system shown in power angular dependence matched curve such as following formula (4):
Wherein, p, b, c are unknown parameter, PmIt is the mechanical output of one machine infinity bus system, PeIt is the electromagnetic work of generator
Rate,For the second dervative of generator rotor angle δ in equation of rotor motion.
It should be noted that the power that control centre receives each generator after breaking down calculates according to the following formula
It arrives:
Wherein:
PeiFor the electromagnetic power of i-th generator in m platform generator.Wherein,For the voltage vector of node i,For section
The conjugate complex number of point i Injection Current,The conjugate complex number of the j element of admittance matrix i-th,To be total to for the voltage vector of node j
Yoke plural number, EiAnd EjFor the interior electromotive force of generator, GiiFor the self-conductance of the i-th node, GijFor the i-th j in admittance matrix members
The conductance of element, DijIt is during failure and after failure for the constant being not zero, Y with CijijTo be retracted to generator interior nodes
The modulus value of the j element of admittance matrix i-th, BijFor the susceptance of the i-th j in admittance matrix elements, δijFor node i and node j it
Between angle, δiAnd δjFor generator's power and angle.
S4 carries out transient stability using EEAC strategy and sentences according to the power angular dependence matched curve after the failure
Not.
Preferably, carrying out transient stability differentiation using EEAC strategy includes: when acceleration region Δ ac > do not take is any
When retardation area Δ RED under control strategy, system unstability takes direct-current emergency control strategy, and calculating is controlled by direct-current emergency
The increased retardation area of institute, Δ DC, otherwise, as Δ ac < Δ RED, system keeps transient stability；When Δ ac > Δ RED+ Δ
When DC, system continues unstability, calculates by cutting the machine control increased retardation area Δ gt of institute, and take the control strategy for the machine of cutting.
Specifically, Fig. 3 is that direct-current emergency power control and machine of cutting control the schematic diagram matched in the present embodiment, referring to figure
3, δ0Generator rotor angle when occurring for failure, δcGenerator rotor angle when for fault clearance, Δ ac are acceleration regions, Δ RED be do not take it is any
Retardation area under control strategy, δqTo execute generator rotor angle when emergency control policy, from generator rotor angle δcTo δqBetween be and meanwhile consider
Data transmission period calculates time and operating delay time, δ 'hTo implement the excision of the limit after direct-current emergency power control
Angle, δ "hIt is to take further to implement to cut the critical clearing angle after machine.Wherein, it is mechanical to be equivalent to reduction for dc power Emergency Assistance
Power increases retardation area under this effect.Add if urgent dc power support cannot be such that the retardation area of unit is greater than
Fast area then needs to be implemented and cuts machine strategy.Δ gt is needed by cutting machine strategy come increased retardation area.It is assumed that dc power
Emergency Assistance amount is Δ PDC, cutting machine amount is Δ PG。
Accelerating area is the area in the region of Δ ac shown in Fig. 3, and the acceleration region is calculated according to following formula (5)
Area Δ ac:
Retardation area is the area in the region Δ RED shown in Fig. 3, calculates described not taking according to following formula (6) and appoints
Retardation area under what control strategy:
As Δ ac < Δ RED, system can keep transient stability, and the system unstability as Δ ac > Δ RED needs at this time
Direct-current emergency control strategy is taken, is calculated by the direct-current emergency control increased retardation area Δ DC of institute.According to following formula (7)
It calculates by the direct-current emergency control increased retardation area Δ DC of institute:
The system unstability as Δ ac > Δ RED+ Δ DC, needs further to take at this time the control strategy for the machine of cutting, according to following
Formula (8) is calculated by cutting the machine control increased retardation area Δ gt of institute:
Wherein,
According to formula (7)-(13), it can acquire and cut machine amount Δ PG, PIFor electromagnetic power under normal circumstances, PIIIFor fault clearance
Electromagnetic power afterwards.
Embodiment 2
Fig. 8 be embodiment 2 it is equivalent after one machine infinity bus system power angle characteristics curve graph；
Fig. 9 is the simulation result diagram of the embodiment of the present invention 2
Fig. 4 is that the three machine systems of use in the embodiment of the present invention 2 connect infinite busbar topological structure schematic diagram, referring to Fig. 4,
Including three fired power generating units.Wherein, generator excitation model uses IEEE-AC1A, and governor model uses IEEE-GOV4, bears
Load is modeled with constant impedance.Direct current system voltage class is 500kV, and under normal circumstances, the control mode of direct current system is: whole
Stream side is used using constant current control and inverter side determines gamma kick.When 10s, occur on the alternating current circuit between b2 and b3
Three-phase shortcircuit.When 10.55s, failure is removed, reference power 100MVA.Using the direct-current emergency control based on EEAC theory
Measure processed, to improve the transient stability of ac and dc systems.Specifically control process includes:
Direct-current emergency control ability is tested.
As in the system of fig. 4, test is made to the emergent control of direct current first.Due to direct-current emergency power control
Fast with speed, the big advantage of regulation capacity, is the important control measure of guarantee system safe and stable operation, so when adjacent
When alternating current circuit is broken down, power support is carried out to faulty line first with direct-current emergency control.Fig. 5 is that the present invention is implemented
The dc power emergency lifting of example 2 and the simulation result diagram of emergency descent, according to Fig. 5, for the mould in the 500Kv direct current system 5s
The emergency lifting of quasi- dc power and decline, Fig. 9 is promoted from 0.5pu. to 0.7pu., 0.5 respectively drops to 0.3pu.,
The emulation controlled direct-current emergency is realized in PSCAD, power ascension or decline result are more accurate.
Carry out machine group identification.It is short that three-phase occurs for an AC line in system as shown in Figure 4,10s between b2 and b3
Road, the maximal clearance after EEAC is arranged using each generator's power and angle from small to large differentiate criticality benchmark.Table 1 is the system of acquisition
Each generator's power and angle data after breaking down.
Table 1: generator's power and angle data
Generator rotor angle data are the opposite data that presteady state generator rotor angle occurs with failure in table 1.According to table 1, after taking failure to occur compared with
The short time, unit generator rotor angle has had the tendency that separated when 0.2s, so 0.2s is selected to carry out machine group identification, can therefore
After barrier occurs when 0.2s, generator is with respect to the ascending arrangement of generator rotor angle are as follows: δ1< δ3< δ2, maximal clearance is in δ3And δ1Between, most
2.7 ° are divided between big.Because of δ3And δ1Between interval it is maximum, so all units are separated between 1 and unit 3.δ3's
Generator rotor angle ratio δ1Greatly, more likely unstability, so being greater than δ for all1Unit be divided into critical machine, be less than or equal to δ1Machine component
For a remaining group of planes.Therefore it can differentiate that G2, G3 are criticality benchmark group, G1 is remaining machine group.
Further it is equivalent to one machine infinity bus system.Fig. 6 is criticality benchmark group and remaining machine group in the embodiment of the present invention 2
Failure after generator's power and angle curve graph, Fig. 7 be the embodiment of the present invention 2 failure after generator electromagnetic power curve graph, reference
Fig. 6 and Fig. 7, failure is at the 10s moment, and 10.55s terminates and open-circuit line two sides switch.
Criticality benchmark group and remaining machine group are equivalent to a dual systems with respective center inertia machine, then by two-shipper
System is equivalent to infinitely great one-of-a-kind system.The power angle characteristics curve of one machine infinity bus system after equivalent is as shown in Figure 8.
Matched curve in Fig. 8 is to be fitted by the generator rotor angle measured of 10.55~10.7s in Fig. 6 and Fig. 7 and electromagnetic power
The fired power generating unit generator rotor angle relation curve arrived.Experiment curv is the generator rotor angle relation curve of practical fired power generating unit.It can be seen that there are micro- for the two
Small difference, but the implementation of control strategy is not influenced.
It is calculated according to EEAC and accelerates area Δ ac=86.15, retardation area Δ RED=62.68, due to Δ ac > Δ
RED, system are not able to maintain transient stability, in order to stablize system, need to take corresponding control measure.
The characteristics of carrying out direct-current emergency control first, playing direct current quick response promotes 100MW, increased retardation area Δ
DC=9.388 accelerates area to be still greater than retardation area due to Δ ac > Δ RED+ Δ DC, and system cannot be made to stablize.
Since system is not able to maintain transient stability still after direct-current emergency control, using machine strategy is cut, Δ gt is calculated
=14.08, being calculated and cutting machine amount is 68MW.Therefore, the thermal motor of 1 100MW is cut off in 10.55s.Fig. 9 is embodiment 2
Simulation result diagram, as shown in Figure 9, it can be seen that only relying on HVDC Modulation or cutting machine strategy has generating set generator rotor angle flyweight
Certain inhibiting effect, but system is insufficient to allow to restore to stablize.The control strategy combined with the machine of cutting using HVDC Modulation is effectively
The generator rotor angle flyweight for inhibiting fired power generating unit, makes system restore stabilization.
In conclusion the embodiment of the present invention passes through extraction system stability direct from electric system actually disturbed track
Qualitative information and quantitative information ensure that the integrality and accuracy of information；This method is not to the complexity of mathematical model simultaneously
The multi- cluster characteristic of degree and dynamic process proposes too many limitation, and calculating speed is fast, high-efficient, can theoretically be controlled
Amount has certain engineering application value.
Those of ordinary skill in the art will appreciate that: attached drawing is the schematic diagram of one embodiment, and the process in attached drawing is simultaneously
It is not necessarily necessary to the implementation present invention.
As seen through the above description of the embodiments, those skilled in the art can be understood that the present invention can
It realizes by means of software and necessary general hardware platform.Based on this understanding, technical solution of the present invention essence
On in other words the part that contributes to existing technology can be embodied in the form of software products, the computer software product
It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer equipment
(can be personal computer, server or the network equipment etc.) executes the certain of each embodiment or embodiment of the invention
Method described in part.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
In the technical scope disclosed by the present invention, any changes or substitutions that can be easily thought of by anyone skilled in the art,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with scope of protection of the claims
Subject to.
Claims (7)
1. a kind of alternating current-direct current combined hybrid system transient stability control method based on EEAC characterized by comprising
The generator rotor angle and power data of each generator after breaking down are received by control centre；
According to the generator rotor angle data, machine group identification is carried out to each generator；
The machine group after identification is equivalent to using corresponding local inertia center according to the generator rotor angle data and power data
Then double machine set systems are equivalent to an one machine infinity bus system by one double machine set system, and by the single machine
Infinite bus system curve matching is the generator rotor angle relationship matched curve after failure；
According to the power angular dependence matched curve after the failure, it is steady that transient state is carried out using extended equal area criterion EEAC strategy
Qualitative discrimination.
2. the method according to claim 1, wherein described is carried out using extended equal area criterion EEAC strategy
Transient stability differentiation includes:
System unstability, takes direct-current emergency control strategy, calculates by the direct-current emergency control increased retardation area Δ DC of institute otherwise, as Δ ac
When < Δ RED, system keeps transient stability；As Δ ac > Δ RED+ Δ DC, system continues unstability, calculates by cutting machine control
The increased retardation area Δ gt of institute, and take the control strategy for the machine of cutting.
3. the method according to claim 1, wherein described according to the generator rotor angle data, to described each
Generator carries out machine group identification, comprising: is arranged using EEAC according to each generator's power and angle, is sentenced from small to large by maximal clearance
Not, each generator is divided into criticality benchmark group and remaining machine group.
4. according to the method described in claim 2, it is characterized in that, described will according to the generator rotor angle data and power data
The criticality benchmark group and remaining machine group use corresponding local inertia center to be equivalent to a double machine set system, comprising:
It is equivalent to the criticality benchmark group to lose the machine set system of transient stability according to following formula (1):
The remaining machine group is equivalent to according to following formula (2) to keep the machine set system of transient stability:
Wherein, MsIt is the equivalent inertia time constant of critical machine, MiBe every generator in criticality benchmark group inertia time it is normal
Number, δsIt is the equivalent generator rotor angle of criticality benchmark group, δiIt is the generator rotor angle of every generator in critical machine, MABe remaining machine group etc.
Imitate inertia time constant, MjIt is the inertia time constant of every generator in a remaining group of planes, δAIt is the equivalent function of remaining machine group
Angle, δjIt is the generator rotor angle that remaining machine group concentrates every generator.
5. according to the method described in claim 4, it is characterized in that, described be equivalent to one for two machine set systems
One machine infinity bus system, comprising: two machine set systems are equivalent to an Infinite bus power system system according to following formula (3)
System:
Wherein, M is the inertia time constant of one machine infinity bus system, and δ is the generator rotor angle of one machine infinity bus system, PmIt is that single machine is infinite
The mechanical output of big system, PmiIt is the mechanical output of every generator in criticality benchmark group, PmjIt is every hair in remaining machine group
The mechanical output of motor.
6. according to the method described in claim 5, it is characterized in that, described by the one machine infinity bus system curve matching
It for the power angular dependence matched curve after failure, specifically includes: according to the electromagnetism of the one machine infinity bus system and generator
Power obtains in image one machine infinity bus system shown in power angular dependence matched curve such as following formula (4):
Wherein, p, b, c are unknown parameter, PmIt is the mechanical output of one machine infinity bus system, Pe is the electromagnetic power of generator,
For in equation of rotor motion, the second dervative of generator rotor angle δ.
7. according to the method described in claim 6, it is characterized in that, described as acceleration region area Δ ac > do not take and appoint
When retardation area Δ RED under what control strategy, system unstability takes direct-current emergency control strategy, calculates by direct-current emergency control
Otherwise, as Δ ac < Δ RED, system keeps transient stability to the increased retardation area Δ DC of system institute；When Δ ac > Δ RED+ Δ
When DC, system continues unstability, calculates by cutting the machine control increased retardation area Δ gt of institute, and take the control strategy for the machine of cutting, wraps
It includes:
The acceleration region area Δ ac is calculated according to following formula (5):
The retardation area that do not take under any control strategy is calculated according to following formula (6):
It is calculated according to following formula (7) by the direct-current emergency control increased retardation area Δ DC of institute:
It is calculated according to following formula (8) by cutting the machine control increased retardation area Δ gt of institute:
Wherein,
δ0Generator rotor angle when occurring for failure, δcGenerator rotor angle when for fault clearance, δqTo execute generator rotor angle when emergency control policy, δ 'h
To implement the critical clearing angle after direct-current emergency power control, δ "hIt is to take further to implement to cut the critical clearing angle after machine,
Δ PG is to cut machine amount, PIFor electromagnetic power under normal circumstances, PIIIFor the electromagnetic power after fault clearance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811346098.0A CN109494765B (en) | 2018-11-13 | 2018-11-13 | EEAC-based AC-DC hybrid system transient stability control method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811346098.0A CN109494765B (en) | 2018-11-13 | 2018-11-13 | EEAC-based AC-DC hybrid system transient stability control method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109494765A true CN109494765A (en) | 2019-03-19 |
| CN109494765B CN109494765B (en) | 2022-03-18 |
Family
ID=65695779
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811346098.0A Active CN109494765B (en) | 2018-11-13 | 2018-11-13 | EEAC-based AC-DC hybrid system transient stability control method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109494765B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111030112A (en) * | 2020-01-02 | 2020-04-17 | 北京交通大学 | Method for judging transient stability of alternating current-direct current hybrid power system containing flexible direct current |
| CN112366689A (en) * | 2020-10-27 | 2021-02-12 | 中国电力科学研究院有限公司 | Method and system for determining absolute inertia of alternating current power grid |
| CN114188956A (en) * | 2021-12-03 | 2022-03-15 | 南瑞集团有限公司 | Transient stability emergency control method and device considering direct current successive blocking time difference and storage medium |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007181387A (en) * | 2005-12-01 | 2007-07-12 | Chubu Electric Power Co Inc | Method for determining power system transient stability and its equipment |
| CN102035204A (en) * | 2010-12-15 | 2011-04-27 | 北京交通大学 | Traction substation calculation method for transient-state analysis of power system |
| CN103311960A (en) * | 2013-06-17 | 2013-09-18 | 国家电网公司 | Method for dividing coherent stability regions of strong-connection power grids |
| CN103606922A (en) * | 2013-12-10 | 2014-02-26 | 国家电网公司 | Approximate judgment method for power angle stability of electric power system based on typical fault set |
| CN103683267A (en) * | 2013-10-29 | 2014-03-26 | 南京南瑞集团公司 | Power support capacity distribution method of different DC lines in multi-DC system |
| CN104331837A (en) * | 2014-08-13 | 2015-02-04 | 国网电力科学研究院 | Simplification method for optimal cutting machine control strategy search of electric system transient stability |
| CN104809512A (en) * | 2014-12-18 | 2015-07-29 | 国电南瑞科技股份有限公司 | Rapid assessment method for performance indexes of different emergency control measures in electric power system |
| CN107181254A (en) * | 2016-11-23 | 2017-09-19 | 国电南瑞科技股份有限公司 | The matrix form quantizating index of plurality of stable problem is calculated and alternating current-direct current control method for coordinating |
| CN107994588A (en) * | 2017-12-08 | 2018-05-04 | 甘肃省电力公司风电技术中心 | Lift the idle control strategy of wind power plant of wind power integration light current net transient stability |
| CN108054768A (en) * | 2017-12-15 | 2018-05-18 | 北京交通大学 | Transient stability evaluation in power system method based on principal component analysis |
-
2018
- 2018-11-13 CN CN201811346098.0A patent/CN109494765B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007181387A (en) * | 2005-12-01 | 2007-07-12 | Chubu Electric Power Co Inc | Method for determining power system transient stability and its equipment |
| CN102035204A (en) * | 2010-12-15 | 2011-04-27 | 北京交通大学 | Traction substation calculation method for transient-state analysis of power system |
| CN103311960A (en) * | 2013-06-17 | 2013-09-18 | 国家电网公司 | Method for dividing coherent stability regions of strong-connection power grids |
| CN103683267A (en) * | 2013-10-29 | 2014-03-26 | 南京南瑞集团公司 | Power support capacity distribution method of different DC lines in multi-DC system |
| CN103606922A (en) * | 2013-12-10 | 2014-02-26 | 国家电网公司 | Approximate judgment method for power angle stability of electric power system based on typical fault set |
| CN104331837A (en) * | 2014-08-13 | 2015-02-04 | 国网电力科学研究院 | Simplification method for optimal cutting machine control strategy search of electric system transient stability |
| CN104809512A (en) * | 2014-12-18 | 2015-07-29 | 国电南瑞科技股份有限公司 | Rapid assessment method for performance indexes of different emergency control measures in electric power system |
| CN107181254A (en) * | 2016-11-23 | 2017-09-19 | 国电南瑞科技股份有限公司 | The matrix form quantizating index of plurality of stable problem is calculated and alternating current-direct current control method for coordinating |
| CN107994588A (en) * | 2017-12-08 | 2018-05-04 | 甘肃省电力公司风电技术中心 | Lift the idle control strategy of wind power plant of wind power integration light current net transient stability |
| CN108054768A (en) * | 2017-12-15 | 2018-05-18 | 北京交通大学 | Transient stability evaluation in power system method based on principal component analysis |
Non-Patent Citations (3)
| Title |
|---|
| XIAO-PING ZHANG: "Multiterminal Voltage-Sourced Converter-Based HVDC Models for Power Flow Analysis", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
| 李文勇: "基于EEAC理论的多机电力系统暂态稳定控制策略研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
| 韦肖燕 等: "采用储能电源辅助的暂态稳定紧急控制方法", 《电工技术学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111030112A (en) * | 2020-01-02 | 2020-04-17 | 北京交通大学 | Method for judging transient stability of alternating current-direct current hybrid power system containing flexible direct current |
| CN111030112B (en) * | 2020-01-02 | 2021-10-26 | 北京交通大学 | Method for judging transient stability of alternating current-direct current hybrid power system containing flexible direct current |
| CN112366689A (en) * | 2020-10-27 | 2021-02-12 | 中国电力科学研究院有限公司 | Method and system for determining absolute inertia of alternating current power grid |
| CN112366689B (en) * | 2020-10-27 | 2022-09-20 | 中国电力科学研究院有限公司 | Method and system for determining absolute inertia of alternating current power grid |
| CN114188956A (en) * | 2021-12-03 | 2022-03-15 | 南瑞集团有限公司 | Transient stability emergency control method and device considering direct current successive blocking time difference and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109494765B (en) | 2022-03-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109494765A (en) | Alternating current-direct current combined hybrid system transient stability control method based on EEAC | |
| Lin et al. | Division algorithm and interconnection strategy of restoration subsystems based on complex network theory | |
| CN111030112A (en) | Method for judging transient stability of alternating current-direct current hybrid power system containing flexible direct current | |
| CN103559347A (en) | Method for establishing electromagnetic transient simulation model of large-scale AC-DC (Alternating Current - Direct Current) power system | |
| Ge et al. | A model and data hybrid-driven short-term voltage stability real-time monitoring method | |
| CN102445660B (en) | Gray Verhulst model-based prediction method of power angle of generator | |
| CN111639841B (en) | Lightning trip-out risk comprehensive evaluation method for high-voltage transmission line | |
| Rezvani et al. | A new approach for steady-state analysis of a hybrid AC-DC microgrid | |
| CN109787266A (en) | Extra-high voltage direct-current complex fault on-line analysis decision-making technique and system | |
| Guo et al. | Real time simulation for the study on smart grid | |
| CN107069747B (en) | A kind of minimum start-up mode based on regional voltage stabilization determines method | |
| CN106021943B (en) | A kind of DC Line Fault screening technique of combination GPU software and hardware architecture features design | |
| Liu et al. | An improved phasor method for modelling and analysis of modular multilevel matrix converter | |
| Chatterjee et al. | Application of trajectory sensitivity for the evaluation of the effect of TCSC placement on transient stability | |
| CN105470980A (en) | Alternating current/direct current coordination control method based on failure of direct current commutation of power grid all the time | |
| CN108199418A (en) | A kind of interconnected network transient stability emergent control cuts machine solution formulation method | |
| TW201415757A (en) | High performance micro grid island operation grounding failure analysis method | |
| Omer et al. | Analysis of Commutation Failure in LCC-HVDC Systems During Unbalanced Conditions Using Sequence Components | |
| Prince et al. | A fault detection technique based on line parameters in ring-configured DC microgrid | |
| Yanxia et al. | Research on protective relaying systems based on multi-agent system | |
| Li et al. | The optimization of undervoltage load shedding strategy | |
| Guo et al. | Lyapunov exponents over variable window sizes for prediction of rotor angle stability | |
| Padmawansa et al. | Improving Transient Stability of an Islanded Microgrid Using PV Based Virtual Synchronous Machines | |
| CN103997049B (en) | A kind of electrical network real time emergency control method | |
| CN103995958B (en) | Method for constructing harmonic wave impedance characteristic function |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |