CN108808701A - A kind of system inertia real-time computing technique based on response - Google Patents
A kind of system inertia real-time computing technique based on response Download PDFInfo
- Publication number
- CN108808701A CN108808701A CN201810750396.XA CN201810750396A CN108808701A CN 108808701 A CN108808701 A CN 108808701A CN 201810750396 A CN201810750396 A CN 201810750396A CN 108808701 A CN108808701 A CN 108808701A
- Authority
- CN
- China
- Prior art keywords
- inertia
- moment
- active
- loss
- response
- 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.)
- Pending
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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The system inertia real-time computing technique based on response that the present invention relates to a kind of, according to the real-time inertia center frequency information of electric system and the busbar information of load bus, when electric system disturbs, frequency can change, and the frequency change rate at inertia center is calculated using real-time inertia centre frequency；At two moment for taking corresponding disturbance initial voltage to change greatly respectively, take the two moment corresponding inertia centre frequency change rate and burden with power；The equation that two moment generator amatures move is subtracted each other, phase shift is divided by obtain system equivalent inertia.The present invention can calculate system equivalent inertia without carrying out preset cutting load action directly according to the response message after the system failure.
Description
Technical field
The present invention relates to technical field of power systems, especially a kind of system inertia real-time computing technique based on response.
Background technology
When generator failure occurs for current power system, generator bus is disconnected with system, and often occurrence frequency is stablized
Sex chromosome mosaicism, traditional low-frequency load shedding method promote the frequency retrieval of system attached to rated value by fixed low frequency load shedding equipment
Closely.And adaptive off-load scheme is made that based on frequency response models and starts sex work, is deduced using frequency response models
The quantitative ratio relation of moment system inertia centre frequency rate of descent and the active vacancy of system after disturbance, proportionality coefficient are system etc.
It is worth inertia, the core algorithm of all adaptive low-frequency load shedding schemes after this proportionate relationship becomes.
Invention content
The system inertia real-time computing technique based on response that the purpose of the present invention is to provide a kind of, to overcome the prior art
Present in defect.
To achieve the above object, the technical scheme is that：A kind of system inertia real-time computing technique based on response,
It realizes in accordance with the following steps：
Step S1：When electric system disturbs, system inertia centre frequency f is obtained by wide-area monitoring systemscoi,
Obtain inertia centre frequency change rateAnd two values are subjected to standardization；
Step S2：When remembering that t1, t2 are respectively that two of disturbance initial voltage variation in preset range occur for electric system
It carves, and obtains the two moment corresponding inertia center change rate and burden with power, burden with power is subjected to standardization；It is corresponding
Two moment of t1, t2, the equation of motion of generator amature are as follows：
Wherein：It is the inertia centre frequency change rate at t1 and t2 moment respectively；ΔPT=t1, Δ PT=t2
It is the active vacancy of system at t1 and t2 moment, MeqThe equivalent inertia based on response that is system；
The active vacancy of system is shown such as following formula：
ΔPT=t1=PL, t=t1+PLoss, t=t1-PM, t=t1
ΔPT=t2=PL, t=t2+PLoss, t=t2-PM, t=t2
Wherein：PL, t=t1, PM, t=t1, PLoss, t=t1It is that t1 moment system loadings are active, machinery is active and loss respectively；
PL, t=t2, PM, t=t2, PLoss, t=t2It is that t2 moment system loadings are active, machinery is active and loss respectively；
Step S3：The generator amature equation at two moment is subtracted each other, goes to indicate with active, the mechanical active and loss of load
Active vacancy, obtains following formula：
Due to the limitation of governor inertial element, lead to (PM, t=t1-PM, t=t2) ignore；Due to the spy of loss itself
Property, lead to (PLoss, t=t1-PLoss, t=t2) ignore；Then ignore after the active variable quantity of machinery and Dissipation change amount system is used
Amount is：
In an embodiment of the present invention, the t1 and intervals t2 500ms.
Compared to the prior art, the invention has the advantages that：A kind of system based on response proposed by the present invention
Inertia real-time computing technique can be directly according to the response message after the system failure, meter without carrying out preset cutting load action
Calculation system equivalence inertia.With existing adaptive low-frequency load shedding scheme comparison, after the present invention can be calculated quickly and is out of order
System equivalent inertia, and then the active vacancy of system is quickly estimated, to preventing the further unstability of frequency from playing an important role.
Description of the drawings
Fig. 1 is IEEE39 node diagrams in the embodiment of the present invention.
Specific implementation mode
Below in conjunction with the accompanying drawings, technical scheme of the present invention is specifically described.
The present invention provides a kind of system inertia real-time computing technique based on response, is surveyed according to wide-area monitoring systems (WAMS)
Inertia center frequency information and the active information of load are measured, generator amature equation of motion computing system inertia is utilized.Specific implementation
Method is：
Step S1：When electric system disturbs, system inertia center frequency is obtained using wide-area monitoring systems (WAMS)
Rate fcoi, and then calculate inertia centre frequency change rateWherein, two values are subjected to standardization.
Step S2：It takes t1, t2 to correspond to electric system respectively and two moment that disturbance initial voltage changes greatly occurs, generally
Take t1, t2 differs 500ms, and takes the two moment corresponding inertia center change rate and burden with power, by burden with power into
Row standardization；At the two corresponding moment, there is the equation of motion of generator amature as follows：
In formula：It is the inertia centre frequency change rate at t1 and t2 moment respectively；ΔPT=t1, Δ PT=t2
It is the active vacancy of system at t1 and t2 moment.MeqFor the equivalent inertia of system.
The active vacancy of system can be indicated such as following formula：
ΔPT=t1=PL, t=t1+PLoss, t=t1-PM, t=t1
ΔPT=t2=PL, t=t2+PLoss, t=t2-PM, t=t2
In formula：PL, t=t1, PM, t=t1, PLoss, t=t1It is that t1 moment system loadings are active, machinery is active and loss respectively；
PL, t=t2, PM, t=t2, PLoss, t=t2It is that t2 moment system loadings are active, machinery is active and loss respectively.
Step S3：The generator amature equation at two moment is subtracted each other, goes to indicate with active, the mechanical active and loss of load
Active vacancy can obtain following formula：
Wherein：The limitation of governor inertial element leads to (PM, t=t1-PM, t=t2) very little, the characteristic of loss itself causes
(PLoss, t=t1-PLoss, t=t2) also very little.
Therefore, ignoring the active variable quantity of machinery and Dissipation change amount can obtain：
The system inertia real-time computing technique based on response thereby is achieved.
Further, in order to allow skilled in the art realises that technical solution proposed by the present invention, with reference to specific reality
Example is applied to illustrate.
IEEE39 node systems are selected, as shown in Figure 1,0.2s moment G3, G9 units fall off simultaneously, while 25-26 lines
Road, 18-3 circuits, 14-4 circuits, 11-6 circuit breakings, the i.e. right half part in IEEE39 node systems establish isolated island, load
Including node 12,15,16,18,20,21,23,24,26,27,28,29, generator includes fired power generating unit G4, G5, G6, G7.Respectively
Generated power and inertia are as shown in table 1, and each load bus is active as shown in table 2.The practical inertia of system is 231.6s at this time.
1 each generated power of table and inertia
| Generator | Active (p.u.) | Inertia (s) |
| G4 | 6.32 | 57.2 |
| G5 | 5.08 | 52.0 |
| G6 | 6.50 | 69.6 |
| G7 | 5.60 | 52.8 |
2 each load bus of table is active
The first step：By system inertia centre frequency, inertia centre frequency change rate is calculated, in t=2.78s, system
Inertia centre frequency has reached the threshold value 49Hz of cutting load action.
Second step：Before system inertia centre frequency drops to 49Hz, it is 500ms to take t1 and t2 time intervals, and
To the burden with power information at two moment.
Third portion：Subtracted each other using two moment generator amature equations, system inertia center frequency value can be calculated and fallen
Before dropping down onto 49Hz, corresponding system equivalent inertia of each moment is obtained.Concrete outcome is shown in Table 3.
The corresponding estimation inertia of 3 different time node of table and its error
As t=2.78s, the inertia of estimation is 236.2s, and the practical inertia of system is 231.6s at this time, error 2%, ten
Split-phase is close.
In the present embodiment, it is only necessary to obtain system inertia centre frequency and burden with power information, occur in electric system
It, can information calculates system equivalent inertia according to response immediately after failure.And existing adaptive low-frequency load shedding scheme, it is being
System breaks down after disturbance, when frequency, which drops to low-frequency load shedding, acts threshold value 49Hz, needs to carry out a preset round
Cutting load acts, and system equivalent inertia can be just calculated according to the secondary cutting load amount.
The above are preferred embodiments of the present invention, all any changes made according to the technical solution of the present invention, and generated function is made
When with range without departing from technical solution of the present invention, all belong to the scope of protection of the present invention.
Claims (2)
1. a kind of system inertia real-time computing technique based on response, which is characterized in that realize in accordance with the following steps：
Step S1：When electric system disturbs, system inertia centre frequency f is obtained by wide-area monitoring systemscoi, obtain
Inertia centre frequency change rateAnd two values are subjected to standardization；
Step S2：Remember that t1, t2 are respectively electric system generation disturbance initial voltage two moment of the variation in preset range, and
The two moment corresponding inertia center change rate and burden with power are obtained, burden with power is subjected to standardization；Corresponding t1, t2
Two moment, the equation of motion of generator amature are as follows：
Wherein：It is the inertia centre frequency change rate at t1 and t2 moment respectively；ΔPT=t1, Δ PT=t2It is t1
With the active vacancy of system at t2 moment, MeqThe equivalent inertia based on response that is system；
The active vacancy of system is shown such as following formula：
ΔPT=t1=PL, t=t1+PLoss, t=t1-PM, t=t1
ΔPT=t2=PL, t=t2+PLoss, t=t2-PM, t=t2
Wherein：PL, t=t1, PM, t=t1, PLoss, t=t1It is that t1 moment system loadings are active, machinery is active and loss respectively；PL, t=t2,
PM, t=t2, PLoss, t=t2It is that t2 moment system loadings are active, machinery is active and loss respectively；
Step S3：The generator amature equation at two moment is subtracted each other, goes to indicate active with active, the mechanical active and loss of load
Vacancy obtains following formula：
Due to the limitation of governor inertial element, lead to (PM, t=t1-PM, t=t2) ignore；Due to the characteristic of loss itself, lead
Cause (PLoss, t=t1-PLoss, t=t2) ignore；Then ignore after the active variable quantity of machinery and Dissipation change amount system inertia is：
2. a kind of system inertia real-time computing technique based on response according to claim 1, which is characterized in that t1 and t2
It is spaced 500ms.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810750396.XA CN108808701A (en) | 2018-07-10 | 2018-07-10 | A kind of system inertia real-time computing technique based on response |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810750396.XA CN108808701A (en) | 2018-07-10 | 2018-07-10 | A kind of system inertia real-time computing technique based on response |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108808701A true CN108808701A (en) | 2018-11-13 |
Family
ID=64075897
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810750396.XA Pending CN108808701A (en) | 2018-07-10 | 2018-07-10 | A kind of system inertia real-time computing technique based on response |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108808701A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110943451A (en) * | 2019-12-12 | 2020-03-31 | 华北电力大学 | System partition area inertia evaluation method based on disturbance data |
| CN111293685A (en) * | 2020-02-28 | 2020-06-16 | 华北电力大学 | System partition inertia evaluation method based on coherent recognition |
| CN111756039A (en) * | 2020-07-09 | 2020-10-09 | 西安交通大学 | New energy power system inertia estimation method based on probability statistics |
| CN111756038A (en) * | 2020-07-09 | 2020-10-09 | 西安交通大学 | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics |
| CN113346482A (en) * | 2021-05-19 | 2021-09-03 | 电子科技大学 | Method for predicting frequency space-time distribution of wide area power system based on SFR model |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102508061A (en) * | 2011-10-22 | 2012-06-20 | 东北电力大学 | Power shortage estimation method based on wide area phasor measurement system |
| CN103956746A (en) * | 2014-03-28 | 2014-07-30 | 西安交通大学 | Automatic-adaptation under-frequency load shedding method based on frequency changing rate response |
-
2018
- 2018-07-10 CN CN201810750396.XA patent/CN108808701A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102508061A (en) * | 2011-10-22 | 2012-06-20 | 东北电力大学 | Power shortage estimation method based on wide area phasor measurement system |
| CN103956746A (en) * | 2014-03-28 | 2014-07-30 | 西安交通大学 | Automatic-adaptation under-frequency load shedding method based on frequency changing rate response |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110943451A (en) * | 2019-12-12 | 2020-03-31 | 华北电力大学 | System partition area inertia evaluation method based on disturbance data |
| CN110943451B (en) * | 2019-12-12 | 2021-05-25 | 华北电力大学 | System partition area inertia evaluation method based on disturbance data |
| CN111293685A (en) * | 2020-02-28 | 2020-06-16 | 华北电力大学 | System partition inertia evaluation method based on coherent recognition |
| CN111756039A (en) * | 2020-07-09 | 2020-10-09 | 西安交通大学 | New energy power system inertia estimation method based on probability statistics |
| CN111756038A (en) * | 2020-07-09 | 2020-10-09 | 西安交通大学 | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics |
| CN111756038B (en) * | 2020-07-09 | 2021-08-13 | 西安交通大学 | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics |
| CN111756039B (en) * | 2020-07-09 | 2021-11-16 | 西安交通大学 | New energy power system inertia estimation method based on probability statistics |
| CN113346482A (en) * | 2021-05-19 | 2021-09-03 | 电子科技大学 | Method for predicting frequency space-time distribution of wide area power system based on SFR model |
| CN113346482B (en) * | 2021-05-19 | 2022-10-18 | 电子科技大学 | Method for predicting wide area power system frequency space-time distribution based on SFR model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108808701A (en) | A kind of system inertia real-time computing technique based on response | |
| EP3446133B1 (en) | Method of determining mutual voltage sensitivity coefficients between a plurality of measuring nodes of an electric power network | |
| US7096175B2 (en) | Stability prediction for an electric power network | |
| CN112257212B (en) | Modeling method and system of comprehensive load model considering distributed power generation system | |
| EP1381132B1 (en) | Power network parameter estimation | |
| CN103248055B (en) | A kind of solve regional power grid off-the-line after correct choosing cut the method for generator and load | |
| JP2016103972A (en) | System and method for modelling load in electrical power network | |
| US11038352B1 (en) | Control system with load selection strategy that minimizes load overshed | |
| CN110912151B (en) | On-pole reactive compensation control method and device for alternating current transmission line | |
| JP7077250B2 (en) | Power system stabilization system | |
| JP2017099187A (en) | Power generation amount estimation device and method for distributed power source | |
| CN107767060B (en) | Theoretical line loss calculation system and method for distribution network line | |
| Adewole et al. | Synchrophasor-based online coherency identification in voltage stability assessment | |
| CN108879668A (en) | A kind of active vacancy real-time computing technique of system based on response | |
| CN109274098B (en) | Simple power system emergency control method considering damping | |
| CN108808700A (en) | A kind of system inertia computational methods based on load local information | |
| US20150120228A1 (en) | Power-system short-circuit capacity monitoring method and system thereof | |
| Adewole et al. | Extended synchrophasor‐based online voltage stability assessment using synchronous generator‐derived indices | |
| CN103944158A (en) | Online power control method based on small interference stable risk detection and system thereof | |
| CN109494718B (en) | Damping-considered emergency control method for complex power system | |
| CN108923434A (en) | A kind of centralized low-frequency load shedding control method based on response | |
| Purwar et al. | Novel adaptive algorithm for optimal relay setting with improved coordination | |
| US11146102B2 (en) | Method and control center arrangement for the control of an electrical energy transmission grid, and computer program product | |
| CN116404923B (en) | Control parameter adjusting device and method of excitation system and storage medium | |
| CN115000924B (en) | Line admittance protection criterion construction method and device for high-proportion new energy system |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181113 |
|
| RJ01 | Rejection of invention patent application after publication |