CN106655201B - A kind of safety-optimized control method of electric system thermostabilization based on security domain - Google Patents
A kind of safety-optimized control method of electric system thermostabilization based on security domain Download PDFInfo
- Publication number
- CN106655201B CN106655201B CN201610838264.3A CN201610838264A CN106655201B CN 106655201 B CN106655201 B CN 106655201B CN 201610838264 A CN201610838264 A CN 201610838264A CN 106655201 B CN106655201 B CN 106655201B
- Authority
- CN
- China
- Prior art keywords
- route
- load
- synchronous generator
- operating point
- security domain
- 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.)
- Active
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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The safety-optimized control method of electric system thermostabilization based on thermal stability security domain that the invention discloses a kind of, determines the N-1 failure that power grid may occur according to electrical network parameter and structural analysis, and obtains the corresponding hyperplane coefficient of each route;Then, the transimission power of each route in electric system is calculated by electric system initial launch point and hyperplane coefficient, whether the transimission power for judging each route is more than its limiting value, judge whether initial launch point is unstability operating point accordingly, and filters out the route of over-limit condition and the Overflow RateHT of computation overload route;Again, filter out the maximum route of Overflow RateHT and its corresponding hyperplane coefficient, the minimum generation adjustment amount of unstability operating point is calculated, if each route adjusted no longer overloads up to stable operating point, cutting load adjustment otherwise is carried out to the Operation of Electric Systems point adjusted.The method of the present invention can refer to conductive Force system operations staff and carry out more efficiently control to the thermal stability unstability operating status of electric system.
Description
Technical field
The invention belongs to electrical network analysis fields, and are related to field of power system control.
Background technique
Economic fast development is continuously increased power demand, and electric power networks scale constantly expands, and gradually forms transregional
Interconnection even transnational interconnected network.China has formed the pattern of all regions Power System Interconnection in addition to Taiwan, European network at present
Plan to further strengthen main grid structure construction within 10 years futures, realizes the interconnection of power grid adjacent with southern east.Simultaneously as one
The secondary energy and growth area distribution are uneven, load center apart from power center farther out so that a large amount of electric power need to pass through electricity
The remote transmission of net.The transimission power of interregional transmission cross-section increasingly close to its transmission limit, threatens the safety and stability of power grid
Operation.Further, the extensive transregional interconnection of power grid increases the factor for influencing power network safety operation, an electric power member
Part breaks down and will likely involve the whole network and in turn result in major accident, brings baptism to the safe and stable operation of power grid.
2003, " 8.14 " beautiful ash moss had occurred in the world in succession, " 8.28 " London has a power failure on a large scale, " 9.1 " Sydney and horse
Come that West Asia has a power failure on a large scale and " 9.28 " Italy has a power failure on a large scale, wherein " 8.14 " beautiful ash moss is North America most serious since the dawn of human civilization
Power outage, being affected by power failure population is more than 50,000,000.The nearly half area electric power paralysis in the whole nation occurs in July, 2012, India
Large-scale blackout is affected by power failure number up to 600,000,000.These serious power outages are all mainly by certain local power at first
Accident causes, and then causes other regional electric power accidents, eventually causes the whole network accident.Massive blackout accident not only results in huge
Huge economic loss, and will seriously affect people's production and living.For complicated interconnected electric power system, force device is many
More, how effective monitoring system, which runs and prevent electric power accident, becomes particularly important.
Power system security domain can provide a large amount of security information, and the appearance of the theory specifies for methods of risk assessment
New direction.It characterizes static system voltage stability domain boundary and complex power injection on heat-staple thermal stability security domain boundary is critical
The set of point can be indicated using hyperplane.In addition, security domain is capable of providing the nargin information of system operating point, to table
It is horizontal to levy system risk, system risk is quantified, computation burden is small, as a result accurately, can be used in quickly calculating safety
Transition probability.And the nargin information of security domain can be used to as management and running personnel in failure process middle filling machine cutting load, with
And it avoids risk during scheduling operation and reliable foundation is provided.
Although having carried out many researchs around security domain and having achieved a large amount of practical achievement, security domain is in warm
It is less using carrying out in the coordinated control field of stability, and the thermal stability for carrying out electric system using security domain there is no to assist
The theory and algorithm of regulation system aspect.
Summary of the invention
The case where losing thermal stability for electric system proposes one kind based on thermostabilization and successively carries out generator safely
The method of power output adjustment and cutting load adjustment, this method can not only be effectively reduced the risk for losing load, while can help
The case where route in system coordinator's guarantee system is not in overlond running.First according to electrical network parameter and structure to being
System carries out N-1 accident analysis, calculates the thermal stability security domain of the corresponding electric system of various possible N-1 malfunctions, obtains
The corresponding hyperplane coefficient of each route of system under each N-1 failure;Then it by system initial launch point and is calculated
The transimission power of each route of hyperplane coefficient computing system, and judge whether each route overloads.Based on thermal stability security domain, to out
The operating point for now overloading route successively calculates synchronous machine power output adjustment amount and load adjustment amount, obtains the adjustable strategies of system, from
And guarantee the stable operation of system.
In order to solve the above-mentioned technical problem, it is safety-optimized to invent a kind of electric system thermostabilization based on thermal stability security domain
Control method, comprising the following steps:
Step 1: determining the N-1 failure that power grid may occur according to electrical network parameter and structural analysis, calculate various possible
The thermal stability security domain of the corresponding electric system of malfunction obtains under each N-1 failure each route pair in electric system
The hyperplane coefficient answered;
Step 2: calculating each route in electric system by the hyperplane coefficient in electric system initial launch point and step 1
Transimission power, and judge whether the transimission power of each route is more than its limiting value, if skipping to step not above limiting value
Five, otherwise, which is unstability operating point, filters out the route of over-limit condition, the route of over-limit condition is
Overload route, and the Overflow RateHT of computation overload route;
Step 3: the maximum route of Overflow RateHT in step 2 and its corresponding hyperplane coefficient are filtered out, then, based on heat
The calculation method for stablizing the minimum generation adjustment amount of security domain calculates the minimum generation adjustment amount of unstability operating point, and judges to adjust
Whether each route no longer overloads in electric system afterwards, if no longer overloading, skips to step 5, otherwise executes step 4;
Step 4: the minimum tangential load amount calculation method based on thermal stability security domain is to power train adjusted in step 3
The operating point of system carries out cutting load adjustment;
Step 5: to obtain stable operating point.
Further, in the present invention, according to the adjustment of synchronous generator power output and electric power in the step 3 and step 4
The sequencing of system cutting load adjustment, and based on the minimum tangential load amount of thermal stability security domain and minimum generation adjustment amount
The particular content of calculation method, the minimum tangential load amount and minimum generation adjustment amount that calculate unstability operating point is as follows:
If step 1) power grid there is a phenomenon where there is circuit overload after N-1 failure, by synchronous generator power output adjustment and
Cutting load adjustment stablizes power system recovery, based on the thermal stability security domain hyperplane in active power injecting power space,
Its mathematic(al) representation are as follows:
α1P1+α2P2+α3P3+…αnPn=1 (1)
In formula (1), PiIt is the active power injection of electric system node i, i=1,2,3 ..., n;αiIt is electric system node
The hyperplane coefficient of i;N is the number of electric system node;
If operating point P (P1, P2..., Pn) be located at except thermal stability security domain boundary, then operating point P (P1, P2..., Pn)
Distance to thermal stability security domain bounding hyperplane HP is minimum tangential load amount or minimum generation adjustment amount Δ P, is equal to and loses
The difference of each active power injection, remembers Δ P=[Δ P between steady operating point and stable operating point adjusted1,ΔP2,…,ΔPn],
Minimum tangential load amount or minimum generation adjustment amount Δ P are orthogonal with thermal stability security domain bounding hyperplane HP, therefore:
ΔP1:ΔP2:...:ΔPn=α1:α2:...:αn (2)
Cutting load operating point P ' (P adjusted1', P2' ..., Pn'), straight line where PP ' indicates at this time are as follows:
Simultaneous solution formula (1) and formula (3) obtain cutting load operating point P ' (P adjusted1', P2' ..., Pn');
It is 0 that step 2), which enables the flag bit Flag of each synchronous generator node,;Synchronous generator power output is adjusted, is to be based on
It is carried out in the thermal stability security domain space of synchronous generator node active power injection, is based on synchronous generator node wattful power
The boundary representation of the thermal stability security domain of rate injection are as follows:
∑αgPg=1- ∑ αlPl (4)
In formula (4), PgIt is injected for the active power of synchronous generator node, PlFor the burden with power numerical value of load bus, αg
For the thermal stability security domain bounding hyperplane coefficient of synchronous generator node, αlIt is super flat for the stabilization security domain boundaries of load bus
Face coefficient;
G is synchronous generator node, and l is load bus;
And only choose the item injected in formula (3) containing synchronous generator active power, i.e.,
NgSynchronous generator number in representation formula (3), Pg(Ng) indicate NgThe wattful power of a synchronous generator node
Rate injection, P ' (Ng) indicate cutting load N adjustedgThe active power of a synchronous generator node is injected, αg(Ng) indicate the
NgThe thermal stability security domain bounding hyperplane coefficient of a synchronous generator node.
Step 3) successively judges whether the flag bit Flag of each synchronous generator node is 0, while judging synchronous generator
Whether the power output of node is greater than 0, chooses synchronous generator node that first meets the two conditions as reference unit ref,
And joint type (4) and formula (5):
Pg(ref) the active power injection of synchronous generator in reference unit ref, α are indicatedg(ref) reference unit is indicated
The thermal stability security domain bounding hyperplane coefficient of synchronous generator in ref.
By formula (6) obtain it is adjusted after operating point in synchronous generator node active power inject vector Pg';
Whether step 4) successively examines the power output of synchronous generator node out-of-limit, if power output is out-of-limit, flag bit setting
It is 1, and the power output of the synchronous generator node is fixed as its limiting value, continues to examine next;If it is no out-of-limit to contribute,
Examine it is next, until having inspected;
Step 5) is by the synchronous generator node in the operating point of result after the adjustment obtained in step 3) and step 4)
Active power injects vector Pg' in be updated;Judge the sum of flag bit of each synchronous generator node sum (Flag) whether etc.
In 0, while judging whether the flag bit Flag (ref) of reference unit is equal to 1, if there is a condition to meet, jump to step 6),
If two conditions are not satisfied, step 3) is skipped to;
Step 6) judges whether the sum of flag bit of each synchronous generator node sum (Flag) is equal to 0, if meeting this
Part has then reached operating point by adjusting synchronous generator power output and has restored to stablize, so far obtained the controlling party of corresponding operating point P
Method;Otherwise, judge whether the flag bit Flag (ref) of reference unit is equal to 1, cutting load adjustment is carried out if meeting the condition;
Step 7) filters out the N that hyperplane coefficient is negative because the load in electric system cannot controllably increaselA load
Node and corresponding flag bit flag, NlThe active power of a load bus injects vector Pl;
If the load bus of step 8) operating point hangs with load, the active power of the load bus is injected to negative value,
It is 0 if without load, successively judges that the power output of each load bus whether less than 0, while judging the flag bit of each load bus
Whether it is 0, chooses first and meet the load bus of the two above-mentioned conditions as reference load ref simultaneously;It is calculated by formula (7)
The active power injection of each node in operating point after cutting load;
Pl(Nl) indicate NlThe active power of a load bus is injected, P 'l(Nl) indicate N adjustedlA load section
The active power injection of point, αl(Nl) NlThe thermal stability security domain bounding hyperplane coefficient of a load bus.Pl(ref) table
Show the active power injection of load bus in reference load ref, αl(ref) reference load node in reference load ref is indicated
Thermal stability security domain bounding hyperplane coefficient.
It is 0 that step 9), which enables the flag bit flag of each load bus,;Successively examining the active power injection of load bus is
No more lower limit, if power output gets over lower limit, the flag bit of load bus is set as 1, and the active power of the load bus is injected
It is fixed as its lower limit value, continues to examine next;If power output is no out-of-limit, examine it is next, until by all load sections
Until point has inspected;
The active power of node each in operating point obtained in step 9) is infused in operating point adjusted by step 10)
Load bus active power inject vector Pl' be updated;Judging the sum of the flag bit of each load bus sum (flag) is
It is no to be equal to 0, while judging whether without circuit overload, if two conditions all meet, reaches operating point and has restored to stablize,
So far the control method of corresponding operating point P is obtained;Otherwise return step 8).
Compared with prior art, the beneficial effects of the present invention are:
The present invention is based on thermal stability security domain theories, carry out to the measure of minimum tangential load and minimum generated output adjustment
It is further perfect, for the thermostabilization instability status of electric system, synchronous generator is controlled orderly, in phase and is carried out
Cutting load ensure that the normal operation of load in system.At the same time, utilization of the invention reduces the risk level of system.
Detailed description of the invention
Fig. 1 is the flow chart of the synchronous machine power output adjustment provided by the invention based on thermal stability security domain;
Fig. 2 (a) and Fig. 2 (b) are minimum generation adjustment amount provided by the invention, minimum tangential load amount schematic diagram respectively;
Fig. 3 is the flow chart of the cutting load adjustment provided by the invention based on thermal stability security domain;
Fig. 4 (a), Fig. 4 (b) reflect 300 node modular system of New England in example of the present invention with Fig. 4's (c) jointly
Wiring diagram;
Fig. 5 (a), Fig. 5 (b) and Fig. 5 (c) reflect route 255-33 in example of the present invention jointly needs the hair carried out after disconnecting
Motor power output regulating measures histogram.
Specific embodiment
Technical solution of the present invention is described in further detail with specific implementation example with reference to the accompanying drawing.
A kind of safety-optimized control method of electric system thermostabilization based on thermal stability security domain of the present invention, is described in detail such as
Under, comprising:
Step 1: the N-1 failure that power grid may occur is determined according to electrical network parameter and structural analysis, is calculated various possible
The thermal stability security domain of the corresponding electric system of N-1 malfunction obtains under each N-1 failure each route in electric system
Corresponding hyperplane coefficient.
By taking 300 node modular system of New England as an example, as shown in Fig. 4 (a), Fig. 4 (b) and Fig. 4 (c), including 69 hairs
Motor node and 62 load bus.411 routes are shared in system, therefore have 411 kinds of N-1 malfunctions, each N-1 event
Remaining 410 routes respectively correspond its thermal stability security domain border coefficient when barrier state.Interface journey is write by MATLAB
Sequence calls the data of 300 node systems in matpower to carry out the calculating on electric system thermal stability security domain boundary.It completes
After the search of critical point, MATLAB, which needs to fit electric system thermostabilization under given N-1 malfunction using critical point, to be pacified
The hyperplane coefficient of gamut boundary is prepared to calculate unstability operating point control method for coordinating.
Since 411 N-1 malfunctions are excessive, 10 N-1 malfunctions shown in table 1 are only enumerated in the patent
For be illustrated, hyperplane coefficient calculating is carried out respectively to this 10 N-1 failures.
1 N-1 failure of table
N-1 failure | N-1 failure |
Route 274-276 is disconnected | Route 260-53 is disconnected |
Route 264-118 is disconnected | Route 98-100 is disconnected |
Route 255-33 is disconnected | Route 96-97 is disconnected |
Route 263-109 is disconnected | Route 116-119 is disconnected |
Route 265-145 is disconnected | Route 39-52 is disconnected |
Step 2: each route in electric system is calculated by the hyperplane coefficient in electric system initial launch point and step 1
Transimission power, and judge whether the transimission power of each route is more than its limiting value, if skipping to step not above limiting value
Five, otherwise, which is unstability operating point, filters out the route of over-limit condition, the route of over-limit condition is
Overload route, and the Overflow RateHT of computation overload route.
To the N-1 failure enumerated in step 1, route 274-276 is disconnected, route 260-53 disconnection not will cause route
Power is out-of-limit, then is adjusted to step 5;Route 264-118 disconnection will cause route 116-119 (Overflow RateHT 52.43%), route
(Overflow RateHT is by 119-161 (Overflow RateHT 47.81%), route 161-118 (Overflow RateHT 4.05%), route 257-43
35.71%) it overloads;Route 255-33 disconnection will cause route 16-36 (Overflow RateHT 3.85%) overload;Route 263-109 is disconnected
Meeting causes route 3-129 (Overflow RateHT 19.30%), route 31-43 (Overflow RateHT 13.66%), route 112-116 (mistake
Load rate be 10.30%), route 109-129 (Overflow RateHT 19.30%), route 257-43 (Overflow RateHT 143.35%) overload;
Route 265-145 disconnection will cause route 257-43 (Overflow RateHT 8.99%) overload;Route 98-100 disconnection will cause route
(Overflow RateHT is by 96-97 (Overflow RateHT 40.12%), route 97-98 (Overflow RateHT 45.72%), route 136-138
27.15%), route 101-136 (Overflow RateHT 3.90%), route 138-96 (Overflow RateHT 40.12%) overload;Route 96-
97 disconnections will cause route 98-100 (Overflow RateHT 19.31%) overload;Route 116-119 disconnection will cause route 4-16 (mistake
Load rate be 20.06%), route 18-36 (Overflow RateHT 34.49%), route 112-116 (Overflow RateHT 27.28%), route 3-
4 (Overflow RateHTs 20.06%) overload;Route 39-52 disconnection will cause route 4-16 (Overflow RateHT 6.09%), route 16-36
(Overflow RateHT 13.16%), route 112-116 (Overflow RateHT 7.92%), route 3-4 (Overflow RateHT 6.09%) overload.
Step 3: filtering out the maximum route of Overflow RateHT in step 2 and its corresponding hyperplane coefficient, then, based on heat
The calculation method for stablizing the minimum generation adjustment amount of security domain calculates the minimum generation adjustment amount of unstability operating point, and judges to adjust
Whether each route no longer overloads in electric system afterwards, if no longer overloading, skips to step 5, otherwise executes step 4.
Step 3-1) if power grid there is a phenomenon where there is circuit overload after N-1 failure, is contributed by synchronous generator and is adjusted
Stablize power system recovery with cutting load adjustment, the thermal stability security domain hyperplane based on active power Injection Space,
Mathematic(al) representation are as follows:
α1P1+α2P2+α3P3+…αnPn=1 (1)
In formula (1), PiIt is the active power injection of electric system node i, i=1,2,3 ..., n;αiIt is electric system node
The hyperplane coefficient of i;N is the number of electric system node;
If operating point P (P1, P2..., Pn) be located at except thermal stability security domain boundary, as shown in Fig. 2 (a) and Fig. 2 (b), then
Operating point P (P1, P2..., Pn) to the distance of thermal stability security domain bounding hyperplane HP be minimum tangential load amount or minimum hair
Electric adjustment amount Δ P is equal to the difference of each active power injection between unstability operating point and stable operating point adjusted, remembers Δ P
=[Δ P1,ΔP2,…,ΔPn], minimum tangential load amount or minimum generation adjustment amount Δ P and thermal stability security domain bounding hyperplane
HP is orthogonal, therefore:
ΔP1:ΔP2:...:ΔPn=α1:α2:...:αn (2)
Cutting load operating point P ' (P adjusted1', P2' ..., Pn'), straight line where PP ' indicates at this time are as follows:
Simultaneous solution formula (1) and formula (3) obtain cutting load operating point P ' (P adjusted1', P2' ..., Pn');
Step 3-2) as shown in Figure 1, enabling the flag bit Flag of each synchronous generator node is 0.Adjust synchronous generator
Power output, is carried out in the thermal stability security domain space injected based on synchronous generator node active power, based on synchronous hair
The boundary representation of the thermal stability security domain of motor node active power injection are as follows:
∑αgPg=1- ∑ αlPl (4)
In formula (4), PgIt is injected for the active power of synchronous generator node, PlFor the burden with power numerical value of load bus, αg
For the thermal stability security domain bounding hyperplane coefficient of synchronous generator node, αlIt is super flat for the stabilization security domain boundaries of load bus
Face coefficient;
In formula (4), g is synchronous generator node, and l is load bus;
And only choose the item injected in formula (3) containing synchronous generator active power, i.e.,
NgSynchronous generator number in representation formula (3), Pg(Ng) indicate NgThe wattful power of a synchronous generator node
Rate injection, P ' (Ng) indicate cutting load N adjustedgThe active power of a synchronous generator node is injected, αg(Ng) indicate the
NgThe thermal stability security domain bounding hyperplane coefficient of a synchronous generator node.
Step 3-3) successively judge whether the flag bit Flag of each synchronous generator node is 0, while judging synchronous generator
Whether the power output of machine node is greater than 0, chooses first and meets the synchronous generator node of the two conditions as reference unit
Ref, and joint type (4) and formula (5):
Pg(ref) the active power injection of synchronous generator in reference unit ref, α are indicatedg(ref) reference unit is indicated
The thermal stability security domain bounding hyperplane coefficient of synchronous generator in ref.By formula (6) obtain it is adjusted after operating point in it is same
It walks generator node active power and injects vector Pg';
Step 3-4) successively examine the power output of synchronous generator node whether out-of-limit, if power output is out-of-limit, flag bit is set
It is set to 1, and the power output of the synchronous generator node is fixed as its limiting value, continues to examine next;If it is no out-of-limit to contribute,
Then examine it is next, until having inspected;
Step 3-5) by the synchronous generator in the operating point of result after the adjustment obtained in step 3-3) and step 3-4)
Machine node active power injects vector Pg' in be updated;Judge the sum of flag bit of each synchronous generator node sum (Flag)
Whether it is equal to 0, while judges whether the flag bit Flag (ref) of reference unit is equal to 1, if there is a condition to meet, jumps to
Step 6) skips to step 3-3 if two conditions are not satisfied);
Step 3-6) judge whether the sum of flag bit of each synchronous generator node sum (Flag) is equal to 0, if meeting this
Part has then reached operating point by adjusting synchronous generator power output and has restored to stablize, so far obtained the controlling party of corresponding operating point P
Method;Otherwise, judge whether the flag bit Flag (ref) of reference unit is equal to 1, cutting load adjustment is carried out if meeting the condition;
After route 255-33 is disconnected, the load factor of route 16-36 is maximum, adjust synchronous generator power output can so that
Operating point restores to stablize, and route 16-36 is just fully loaded at this time, so far obtains the control method of the operating point;
After route 265-145 is disconnected, the load factor of route 257-43 is maximum, and the power output for adjusting synchronous generator can just make
It obtains operating point to restore to stablize, route 257-43 is just fully loaded at this time, so far obtains the control method of the operating point;
After route 96-97 is disconnected, the load factor of route 98-100 is maximum, and the power output for adjusting synchronous machine can be so that operation
Point restores to stablize, and route 98-100 is just fully loaded at this time, so far obtains the control method of the operating point;
After route 116-119 is disconnected, the load factor of route 18-36 is maximum, and the power output for adjusting synchronous machine can be so that fortune
Row point restores to stablize, and route 4-16, route 18-36, route 112-116, route 3-4 can be operated normally at this time, so far be obtained
The control method of the operating point.
The synchronized machine adjustment of four above-mentioned route N-1 failures can restore stable, not need to carry out cutting load adjustment again.
For after route 255-33 disconnection, the control method of acquisition such as Fig. 5 (a), Fig. 5 (b) and Fig. 5 (c) are shown.
After route 264-118 is disconnected, the load factor of route 116-119 is maximum, after synchronous generator power output adjustment, line
Road 116-119 full-load run, route 263-109 heavy service, however route 119-161 (Overflow RateHT 44.27%), route
161-118 (Overflow RateHT 0.51%) still overloads, therefore after route 264-118 disconnection, is only contributed by synchronous generator
Adjustment cannot make system restore to stablize, and need to carry out cutting load adjustment;
After route 263-109 is disconnected, the load factor of route 257-43 is maximum, after synchronous generator power output adjustment, line
Road 3-129, route 109-129, route 257-43, route 31-43 are operated normally, however (Overflow RateHT is route 112-116
23.10%) it still overloads, therefore after route 263-109 disconnection, only system cannot be made extensive by synchronous generator power output adjustment
It is multiple to stablize, it needs to carry out cutting load adjustment;
After route 98-100 is disconnected, the load factor of route 97-98 is maximum, after synchronous generator power output adjustment, route
96-97, route 97-98, route 136-138, route 138-96 are operated normally, however (Overflow RateHT is route 101-136
3.68%) it still overloads, therefore after route 98-100 disconnection, only system cannot be made to restore by synchronous generator power output adjustment
Stablize, needs to carry out cutting load adjustment;
After route 116-119 is disconnected, the load factor of route 18-36 is maximum, after synchronous generator power output adjustment, route
4-16, route 16-36, route 3-4 are operated normally, however route 112-116 (Overflow RateHT 2.90%) still overloads, therefore
After route 116-119 is disconnected, only system cannot be made to restore to stablize by synchronous generator power output adjustment, need to carry out cutting load tune
It is whole.
Step 4: the minimum tangential load amount calculation method based on thermal stability security domain is to power train adjusted in step 3
The operating point of system carries out cutting load adjustment.As shown in Figure 3.
Step 4-1) because the load in electric system cannot controllably increase, filter out the N that hyperplane coefficient is negativelIt is a negative
Lotus node and corresponding flag bit flag, NlThe active power of a load bus injects vector Pl;
Step 4-2) if the load bus of operating point hangs with load, the active power injection of the load bus is negative
Value, is 0 if without load, successively judges that the power output of each load bus whether less than 0, while judging the mark of each load bus
Whether position is 0, chooses first and meets the load bus of the two above-mentioned conditions as reference load ref simultaneously;It is counted by formula (7)
The active power injection of each node in operating point after calculating cutting load;
Pl(Nl) indicate NlThe active power of a load bus is injected, P 'l(Nl) indicate N adjustedlA load section
The active power injection of point, αl(Nl) NlThe thermal stability security domain bounding hyperplane coefficient of a load bus.Pl(ref) table
Show the active power injection of load bus in reference load ref, αl(ref) reference load node in reference load ref is indicated
Thermal stability security domain bounding hyperplane coefficient.
Step 4-3) to enable the flag bit flag of each load bus be 0;Successively examine the active power injection of load bus
Whether lower limit is got over, if power output gets over lower limit, the flag bit of load bus is set as 1, and the active power of the load bus is infused
Enter to be fixed as its lower limit value, continues to examine next;If power output is no out-of-limit, examine it is next, until by all loads
Until node has inspected;
Step 4-4) active power of node each in operating point obtained in step 4-3) is infused in operation adjusted
The active power of load bus in point injects vector Pl' be updated;Judge the sum of the flag bit of each load bus sum
(flag) whether it is equal to 0, while judges whether to have reached operating point if two conditions all meet without circuit overload
Restore to stablize, so far obtains the control method of corresponding operating point P;Otherwise return step 4-2).
Synchronous generator power output is only adjusted in step 3 and is unable to four stable route N-1 failures of recovery system, passes through step
Stable operation may be implemented after rapid four cutting load adjustment.By taking route 264-118 disconnection as an example, operating point needs the hair carried out
Motor is contributed, and regulating measures are as shown in table 2, and the cutting load measure needed is as shown in table 3, and the load bus that do not list does not need
Cutting load adjustment is carried out, i.e. cutting load amount is 0.
2 route 264-118 of table needs the generator output regulating measures carried out after disconnecting
3 route 264-118 of table needs the cutting load regulating measures carried out after disconnecting
Step 5: to finally obtain stable operating point.
To sum up, the present invention is based on thermal stability security domain theories, the measure to minimum tangential load and minimum generated output adjustment
Carried out it is further perfect, for the thermostabilization instability status of electric system, orderly, in phase control synchronous generator, with
And cutting load is carried out, it ensure that the normal operation of load in system.At the same time, utilization of the invention reduces load in system
Loss, to reduce the risk level of system.
Although above in conjunction with figure, invention has been described, and the invention is not limited to above-mentioned specific embodiment parties
Formula, the above mentioned embodiment is only schematical, rather than restrictive, and those skilled in the art are in this hair
Under bright enlightenment, without deviating from the spirit of the invention, many variations can also be made, these belong to guarantor of the invention
Within shield.
Claims (2)
1. a kind of safety-optimized control method of electric system thermostabilization based on thermal stability security domain, which is characterized in that including with
Lower step:
Step 1: determining the N-1 failure that power grid may occur according to electrical network parameter and structural analysis, various possible N-1 are calculated
The thermal stability security domain of the corresponding electric system of malfunction obtains under each N-1 failure each route pair in electric system
The hyperplane coefficient answered;
Step 2: calculating the biography of each route in electric system by the hyperplane coefficient in electric system initial launch point and step 1
Defeated power, and judge whether the transimission power of each route is more than its limiting value, if skipping to step 5 not above limiting value,
Otherwise, which is unstability operating point, filters out the route of over-limit condition, and the route of over-limit condition is to overload
Route, and the Overflow RateHT of computation overload route;
Step 3: filtering out the maximum route of Overflow RateHT in step 2 and its corresponding hyperplane coefficient, then, it is based on thermostabilization
The calculation method of the minimum generation adjustment amount of security domain calculates the minimum generation adjustment amount of unstability operating point, and judges electricity after adjustment
Whether each route no longer overloads in Force system, if no longer overloading, skips to step 5, otherwise executes step 4;
Step 4: the minimum tangential load amount calculation method based on thermal stability security domain is to electric system adjusted in step 3
Operating point carries out cutting load adjustment;
Step 5: to obtain stable operating point.
2. the safety-optimized control method of electric system thermostabilization based on thermal stability security domain according to claim 1, special
Sign is, is contributed the successive of adjustment and the adjustment of electric system cutting load in the step 3 and step 4 according to synchronous generator
Sequentially, the calculation method of minimum tangential load amount and minimum generation adjustment amount and based on thermal stability security domain, calculates unstability operation
The particular content of the minimum tangential load amount of point and minimum generation adjustment amount is as follows:
If step 1) power grid there is a phenomenon where there is circuit overload after N-1 failure, is contributed by synchronous generator and adjusts and cut negative
Lotus adjusts to stablize power system recovery, the thermal stability security domain hyperplane based on active power Injection Space, mathematical table
Up to formula are as follows:
α1P1+α2P2+α3P3+…αnPn=1 (1)
In formula (1), PiIt is the active power injection of electric system node i, i=1,2,3 ..., n;αiIt is electric system node i
Hyperplane coefficient;N is the number of electric system node;
If operating point P (P1, P2..., Pn) be located at except thermal stability security domain boundary, then operating point P (P1, P2..., Pn) arrive heat
The distance for stablizing security domain boundaries hyperplane HP is minimum tangential load amount or minimum generation adjustment amount Δ P, is equal to unstability and transports
The difference of each active power injection, remembers Δ P=[Δ P between row point and stable operating point adjusted1,ΔP2,…,ΔPn], it is minimum
Cutting load amount or minimum generation adjustment amount Δ P are orthogonal with thermal stability security domain bounding hyperplane HP, therefore:
ΔP1:ΔP2:...:ΔPn=α1:α2...:...:αn (2)
Cutting load operating point P ' (P adjusted1', P2' ..., Pn'), straight line where PP ' indicates at this time are as follows:
Simultaneous solution formula (1) and formula (3) obtain cutting load operating point P ' (P adjusted1', P2' ..., Pn');
It is 0 that step 2), which enables the flag bit Flag of each synchronous generator node,;Synchronous generator power output is adjusted, is based on synchronization
It carries out in the thermal stability security domain space of generator node active power injection, is infused based on synchronous generator node active power
The boundary representation of the thermal stability security domain entered are as follows:
∑αgPg=1- ∑ αlPl (4)
In formula (4), PgIt is injected for the active power of synchronous generator node, PlIt is injected for the active power of load bus, αgIt is same
Walk the thermal stability security domain bounding hyperplane coefficient of generator node, αlFor the stabilization security domain boundaries hyperplane system of load bus
Number;
And only choose the item injected in formula (3) containing synchronous generator active power, i.e.,
NgSynchronous generator number in representation formula (3), Pg(Ng) indicate NgThe active power of a synchronous generator node is infused
Enter, P ' (Ng) indicate cutting load N adjustedgThe active power of a synchronous generator node is injected, αg(Ng) indicate NgIt is a
The thermal stability security domain bounding hyperplane coefficient of synchronous generator node;
Step 3) successively judges whether the flag bit Flag of each synchronous generator node is 0, while judging synchronous generator node
Power output whether be greater than 0, choose first and meet the synchronous generator node of the two conditions as reference unit ref, parallel connection
Vertical (4) and formula (5):
Pg(ref) the active power injection of synchronous generator in reference unit ref, α are indicatedg(ref) it indicates in reference unit ref
The thermal stability security domain bounding hyperplane coefficient of synchronous generator;
By formula (6) obtain it is adjusted after operating point in synchronous generator node active power inject vector Pg';
Whether step 4) successively examines the power output of synchronous generator node out-of-limit, if power output is out-of-limit, flag bit is set as 1,
And the power output of the synchronous generator node is fixed as its limiting value, continue to examine next;If it is no out-of-limit to contribute, examine
It is next, until having inspected;
Step 5) is active by the synchronous generator node in the operating point of result after the adjustment obtained in step 3) and step 4)
Power injects vector Pg' in be updated;Judge whether the sum of flag bit of each synchronous generator node sum (Flag) is equal to 0,
Judge whether the flag bit Flag (ref) of reference unit is equal to 1 simultaneously, if there is a condition to meet, jumps to step 6), if two
A condition is not satisfied, then skips to step 3);
Step 6) judges whether the sum of flag bit of each synchronous generator node sum (Flag) is equal to 0, if meeting the condition,
Reach operating point by adjusting synchronous generator power output to restore to stablize, so far obtains the control method of corresponding operating point P;It is no
Then, judge whether the flag bit Flag (ref) of reference unit is equal to 1, cutting load adjustment is carried out if meeting the condition;
Step 7) filters out the N that hyperplane coefficient is negative because the load in electric system cannot controllably increaselA load bus
And corresponding flag bit flag, NlThe active power of a load bus injects vector Pl;
If the load bus of step 8) operating point hangs with load, the active power of the load bus is injected to negative value, if not having
Having load is then 0, successively judge the power output of each load bus whether less than 0, while judge each load bus flag bit whether
It is 0, chooses first and meet the load bus of the two above-mentioned conditions as reference load ref simultaneously;It is cut by formula (7) calculating negative
The active power injection of each synchronous generator node in operating point after lotus;
Pl(Nl) indicate NlThe active power of a load bus is injected, P 'l(Nl) indicate N adjustedlA load bus
Active power injection, αl(Nl) indicate NlThe thermal stability security domain bounding hyperplane coefficient of a load bus;Pl(ref) table
Show the active power injection of load bus in reference load ref, αl(ref) reference load node in reference load ref is indicated
Thermal stability security domain bounding hyperplane coefficient;
It is 0 that step 9), which enables the flag bit flag of each load bus,;Successively examine whether the active power injection of load bus gets over
Lower limit, if power output gets over lower limit, the flag bit of load bus is set as 1, and the active power of the load bus is injected and is fixed
For its lower limit value, continue to examine next;If power output is no out-of-limit, examine next, is examined until by all load bus
Until having tested;
The active power of synchronous generator node each in operating point obtained in step 9) is infused in fortune adjusted by step 10)
The active power of load bus in row point injects vector Pl' be updated;Judge the sum of the flag bit of each load bus sum
(flag) whether it is equal to 0, while judges whether to have reached operating point if two conditions all meet without circuit overload
Restore to stablize, so far obtains the control method of corresponding operating point P;Otherwise return step 8).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610838264.3A CN106655201B (en) | 2016-09-20 | 2016-09-20 | A kind of safety-optimized control method of electric system thermostabilization based on security domain |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610838264.3A CN106655201B (en) | 2016-09-20 | 2016-09-20 | A kind of safety-optimized control method of electric system thermostabilization based on security domain |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106655201A CN106655201A (en) | 2017-05-10 |
CN106655201B true CN106655201B (en) | 2019-11-22 |
Family
ID=58852516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610838264.3A Active CN106655201B (en) | 2016-09-20 | 2016-09-20 | A kind of safety-optimized control method of electric system thermostabilization based on security domain |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106655201B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107508289B (en) * | 2017-09-04 | 2020-12-08 | 天津大学 | Practical and fast dynamic security domain algorithm considering tripping time difference at two ends of line |
CN108493921B (en) * | 2018-01-19 | 2021-04-23 | 中国南方电网有限责任公司 | Method for quickly constructing thermal stability security domain based on key node injection power space |
CN109542972B (en) * | 2018-09-29 | 2023-01-03 | 天津大学 | Multi-dimensional space thermal stability security domain representation method based on correlation model |
CN109193813B (en) * | 2018-09-30 | 2021-11-02 | 中国南方电网有限责任公司 | Radar map drawing method for thermal stability security domain of electric power system |
CN110365048B (en) * | 2019-07-19 | 2022-11-29 | 昆明理工大学 | Black-start net rack dynamic stability assessment method based on optimized trend |
CN113097994A (en) * | 2021-03-15 | 2021-07-09 | 国网浙江省电力有限公司 | Power grid operation mode adjusting method and device based on multiple reinforcement learning agents |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102231174A (en) * | 2011-06-29 | 2011-11-02 | 天津大学 | Method for calculating minimum tangential load |
CN103106338A (en) * | 2013-01-17 | 2013-05-15 | 天津大学 | Method of fast generating boundary of electric system thermal stability security domain in decision space |
CN103607375A (en) * | 2013-10-28 | 2014-02-26 | 天津大学 | Network N-1 security-region-boundary calculation and security evaluation method |
CN104899798A (en) * | 2015-06-30 | 2015-09-09 | 天津大学 | Transient risk control method based on consideration on rotary standby wind power integration system |
-
2016
- 2016-09-20 CN CN201610838264.3A patent/CN106655201B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102231174A (en) * | 2011-06-29 | 2011-11-02 | 天津大学 | Method for calculating minimum tangential load |
CN103106338A (en) * | 2013-01-17 | 2013-05-15 | 天津大学 | Method of fast generating boundary of electric system thermal stability security domain in decision space |
CN103607375A (en) * | 2013-10-28 | 2014-02-26 | 天津大学 | Network N-1 security-region-boundary calculation and security evaluation method |
CN104899798A (en) * | 2015-06-30 | 2015-09-09 | 天津大学 | Transient risk control method based on consideration on rotary standby wind power integration system |
Non-Patent Citations (2)
Title |
---|
An improved risk assessment considering spinning reserve of power system with wind farm integration;Yuan Zeng等;《TENCON 2015 - 2015 IEEE Region 10 Conference》;20151130;第28-33页 * |
基于安全域的电网最小切负荷计算方法;王菲等;《中国电机工程学报》;20100505;第30卷(第13期);第1-4页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106655201A (en) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106655201B (en) | A kind of safety-optimized control method of electric system thermostabilization based on security domain | |
CN109936133B (en) | Power system vulnerability analysis method considering information and physics combined attack | |
CN106469908B (en) | A kind of comprehensive aid decision-making method of the static security that multiple faults is coordinated | |
Ustun et al. | Differential protection of microgrids with central protection unit support | |
CN106549425B (en) | A kind of failure isolated island fast control method containing flexible direct current | |
CN104810853B (en) | Digraph-based islanding method of distribution network including distributed generations | |
CN107134785B (en) | Power transmission network voltage coordination optimization control method considering network structure optimization | |
CN110896231B (en) | Distributed photovoltaic capacity receiving calculation method and system for power distribution network in poverty alleviation area | |
CN105281344A (en) | Smart distribution network self-restoration optimization method considering power quality and uncertainty constraint thereof | |
Yakout et al. | Enhancement of frequency stability of power systems integrated with wind energy using marine predator algorithm based PIDA controlled STATCOM | |
CN114336674B (en) | Distributed toughness frequency control method for alternating-current micro-grid | |
CN107248752B (en) | high-permeability photovoltaic distributed voltage control method based on network topology identification | |
CN110932281A (en) | Multi-section cooperative correction method and system based on quasi-steady-state sensitivity of power grid | |
CN105186496B (en) | A kind of microgrid group running status appraisal procedure and system | |
CN107332277B (en) | Active power distribution network island operation method considering source load storage operation characteristics | |
CN110165667A (en) | Meter and the transmission & distribution collaboration idle work optimization method and system of Static Voltage Security constraint | |
CN108400597B (en) | Direct-current fault static safety analysis method considering security control and primary frequency modulation characteristics | |
CN113890039B (en) | Multi-terminal flexible direct-current power distribution network power flow scheduling optimization method | |
CN102638036A (en) | Assistant decision-making method for fast branch overloads in plan security check | |
CN107451743B (en) | Active power distribution network island operation method considering source load operation characteristics | |
CN110912199A (en) | Point distribution and constant volume method and device for multi-direct-current feed-in receiving-end power grid synchronous phase modifier | |
CN109377020B (en) | Power transmission network planning method considering load transfer capacity of power distribution network | |
Yu et al. | Inverse time overcurrent optimization of distribution network with DG based on improved gray wolf algorithm | |
CN107069757B (en) | A kind of mixing multi-infeed HVDC system secondary voltage control method | |
CN110768301B (en) | Micro-grid frequency synchronization anti-attack cooperative control method |
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 |