CN101387679B - Method for weighting dynamic voltage stability margin index - Google Patents
Method for weighting dynamic voltage stability margin index Download PDFInfo
- Publication number
- CN101387679B CN101387679B CN2008102224459A CN200810222445A CN101387679B CN 101387679 B CN101387679 B CN 101387679B CN 2008102224459 A CN2008102224459 A CN 2008102224459A CN 200810222445 A CN200810222445 A CN 200810222445A CN 101387679 B CN101387679 B CN 101387679B
- Authority
- CN
- China
- Prior art keywords
- load
- power
- dynamic
- growth
- burden
- 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
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a method for measuring dynamic voltage stability margin index, while the dynamic voltage stability margin index is the margin to the voltage collapse point of a power system which load increase slowly in a certain speed; the dynamic voltage stability margin index considers that the system is affected by the dynamic components such as the on-load tap changer (ULTC) and the over excitation limit of engines. The method can calculate the dynamic margins of the system which load increases in various types such as single node, regional or full-network type, to provide a visual measurement presenting the distance between the prior operation point of the system and the voltage collapse point, while the distance and the dynamic margin have linear relationship. Therefore, the dynamic margin index has high practical values.
Description
Technical field
The present invention relates to a kind of method of weighting dynamic voltage stability margin index, belong to the simulation and the computing field of electric system.
Background technology
Stablize the progress of research work along with voltage; Emulation tool and modeling work are progressively perfect; China's operation of power networks and planning department are paid special attention to the stable simulation study work of medium-term and long-term voltage, and the important content of in-depth stability Calculation is classified in also clearly that dynamic electric voltage is stable check research as in " power system safety and stability guide rule ".
Weigh system stability and adopt the nargin index more; So-called nargin index is meant the given running status from system; According to certain pattern; Growth through load growth or through-put power progressively approaches the collapse of voltage point, and the current operating point of system can be used as the index of voltage degree of stability to the distance of collapse of voltage point.Traditional static load nargin is based on the continuous tide method, from present operating point, constantly increases with load, finds the solution trend successively, up to through critical point, tries to achieve load margin thereby obtain node PU curve.But the quiescent voltage analytical approach is based on stable state power flow equation or expansion power flow equation; Ignore Electrical Power System Dynamic and regulated effect of factors; With the behavior of algebraic equation descriptive system, the judgement system keeps the stable ability of voltage under microvariations, and the result who therefore obtains is relatively conservative.
The nargin index has the following advantages: can provide one to represent the current operating point of system measuring to collapse of voltage point distance more intuitively to the operations staff; System's operating point is big or small linear to the distance of collapse of voltage point and nargin index; Can take into account the influence that the distribution of work, load growth mode etc. are arranged of various factors in the transient process such as constraint condition, generator more easily.
When calculating the nargin index, the power of each load bus can increase according to any-mode in the network, to approach collapse point, in order to simplify calculating, often supposes that load power increases with following four kinds of modes:
(1) active power of single load bus and/or reactive power increase, and the power of other load buses remains unchanged;
(2) the meritorious and/or reactive power of the load bus of selection area increases, and other some power of load remain unchanged;
(3) the meritorious and/or reactive power of the load bus of a certain selection area increases, and the active power of the load bus of another selection area and/or reactive power reduce, and other load bus power remain unchanged;
(4) all the active power and/or the reactive power of load bus increase simultaneously.
Dynamic voltage stability margin index is to find the solution electric system when load slowly increases with given pace, and system is apart from the nargin of collapse of voltage point; Dynamic voltage stability margin has been considered the influence of dynamic element such as ULTC (ULTC), the restriction of generator overexcitation in the system; Can calculate the dynamic nargin of system under the multiple load growth modes such as single node, subregion and the whole network; For the operations staff provides one to represent the current operating point of system more intuitively to the measuring of collapse of voltage point distance, system's operating point apart from collapse point with dynamically nargin is linear.Therefore, dynamic voltage stability margin index has actual preferably directive significance.
Electric power system full dynamic simulation is meant electric system electromechanical transient, medium and long term process organically united and carries out Digital Simulation.It described electric system be disturbed after whole continuous dynamic process.The differential algebraic equations algorithm use of full dynamic simulation program becomes GEAR method, network equation and the differential equation simultaneous solution of rank variable step.
The main calculation procedure of GEAR method is following:
For the 1 rank differential equation,
(1) definition Nordsieck vector z
m, it is higher derivative y ', y with back " ... wait the y or f (y, t) data that replace former steps.Promptly
(2) prediction
z
M+1 (0)=Pz
mP is a Pascal triangle battle array,
(3) proofread and correct
z
m+1 (i+1)=z
m+1 (i)+LG
m+1 (i)
Wherein: L is a K dimension constant vector,
(4) truncation error
The truncation error in K rank method m step is:
Wherein: C
K+1The constant relevant with rank.
Be the Z that the m step calculates gained
mZ with back
M-1The value that last component of these two vectors subtracts each other.
The normal relative error of using in the emulation, the relative error that promptly requires per step is less than a value ε who stipulates in advance
0:
Y wherein
MaxThe maximal value of the y that had occurred when being integrated to this step, but when the beginning integration, if y=0 then should get y
Max=1.
Under the situation of differential equation group, more than one of y supposes total N variable y (1), y (2) ... Y (N), then relative error is defined as:
ε≤ε
0, think that then this step is effective, change down the step over to, and consider to become rank, variable step.
ε>=ε
0, think that then this step is invalid, become rank, variable step, recomputate.
(5) become rank, variable step control
After per step integration was accomplished, program was estimated current rank respectively, is risen the truncation error of rank and depression of order variable afterwards.Through certain strategy, determine order and step-length when next step calculates.The constant rank only formula during variable step are following:
New step-length is: h
s=R
k* h.Wherein:
Summary of the invention
The present invention based on institute of system of China Electric Power Research Institute power system analysis software-electric power system full dynamic simulation program, has proposed dynamic voltage stability margin index through after researching and analysing.This index has obtained application in the stable research of Chongqing electricity grid voltage, obtained effect preferably.
The method for solving of existing power system voltage stability margin mostly can not taking into account system in the influence of dynamic element; Particularly actuation times such as ULTC (ULTC), generator excitation restriction are than long protection link; Therefore; Based on the dynamic voltage stability margin index that the electric power system full dynamic simulation program is calculated, this index is the ruuning situation of closing to reality system more.
With reference to Fig. 1, dynamic voltage stability margin index is following in the implementation method of electric power system full dynamic simulation program:
(1) builds network topology structure according to computation requirement, set up the load model that system-computed needs, comprise static load model and induction motor model;
Setting up system's static load model is:
Wherein: P
0Be burden with power, P
1, P
2, P
3, P
4, P
5Be respectively constant impedance burden with power ratio, steady current burden with power ratio, firm power burden with power ratio, relevant burden with power ratio with frequency, the relevant burden with power ratio of voltage index, V is system's virtual voltage, V
0Be system reference voltage, Δ f is a frequency variation, L
DPBe the meritorious percent change that change of frequency 1% causes, N
PBe the active voltage index, load or burden without work model formation parameter is named in like manner, in the following formula: P
1+ P
2+ P
3+ P
4+ P
5=1, Q
1+ Q
2+ Q
3+ Q
4+ Q
5=1.
Set up system's motor load model:
According to known motor stator resistance R
S, motor stator reactance X
S, excitatory reactance X
M, rotor resistance R
R, rotor reactance X
R, motor slippage S, calculating machine moment coefficient: A, B, C;
(2) calculate preparation
Confirm fault type, fault-time, fill in corresponding fault card,, fill in cards such as corresponding ULTC (ULTC) card, generator overexcitation restriction card, relay protection action according to system's actual conditions; Accomplish data preliminary work;
(3) import according to calculation requirement: the zone of the speed of load growth, the time of load growth, load growth or bus title,
Wherein: P is the load growth total amount, and Δ P is a unit interval internal loading increment, and Δ t is the unit interval, T
ENDBe the load growth duration,
According to above-mentioned formula construction load sustainable growth model, get final product through the proportionate relationship of controlling Δ P and Δ t;
(4) simulation calculation
A. read in data and network parameter, inspection Rational Parameters property, the model corresponding relation of formation state variable, intermediate variable and parameter constant and each element carries out initialization process to state variable;
B. form admittance matrix and Jacobi matrix, differentiate the rationality of starting condition;
C. revise the Jacobi matrix element of network equation according to the fault card of having filled in, confirm system state;
D. the GEAR method in the application background technology, calculate the prediction term of GEAR method:
z
M+1 (0)=Pz
m(P is a Pascal triangle battle array)
Calculate the correction term in the GEAR method:
z
m+1 (i+1)=z
m+1 (i)+LG
m+1 (i)
Wherein: L is a K dimension constant vector, and m, i are positive integer,
According to the correction term that obtains, with the differential equation and algebraic equation simultaneous solution;
E. carry out the variable step processing based on calculating actual conditions;
F. simulation calculation finishes, and exports required result;
(5) interpretation of result
Obtain the laggard line data analysis of simulation result, confirm the system power limits to growth,,, obtain by the computes dynamic voltage stability margin index according to the node power limit that calculates and system's normal condition lower node initial power from the node power growth chart:
K
p-dynamic voltage stability margin index (%)
P
0-system initial power (MW)
P
1-system limits power (MW)
P
0Can obtain P by system's normal operating mode statistics
1Result by simulation analysis obtains by whole process simulation calculating back.
Three, dynamic voltage stability margin index characteristics
1. dynamic voltage stability margin index has considered that dynamic element is to factors such as the influence of voltage stability, particularly ULTC (ULTC), generator excitation limiting elements in the system, and the result who obtains more presses close to system's actual conditions.
2. but the different modes of load growth in the computing system, like modes such as single busbar, zone and the whole networks, load growth speed also can be set arbitrarily, can the multiple method of operation of simulation system.
3. assumed load initial, the concluding time that increase arbitrarily during the system loading growth pattern after simulated failure.
Advantage of the present invention:
Dynamically nargin index clear physics conception has good calculability and suit property, can calculate the voltage stability margin of system under the sustainable growth of loading in the various electrical networks easily.
Description of drawings
Fig. 1 is the full dynamic simulation program flow diagram;
Fig. 2 is 3 machines, 9 node system wiring diagrams;
Fig. 3 is a bus A load variations curve map.
Embodiment
With 3 machines, 9 node single systems (wiring diagram such as accompanying drawing 2) is simulation example, calculates the performing step that dynamic nargin further specifies this method:
1. set up the data of 3 machines, 9 node systems, confirm that the load model that system-computed needs is 100% constant-impedance model;
2. calculate and prepare
Confirm that system loading is the constant-impedance load, load growth requires to fill in exciter current of generator protection card and ULTC (ULTC) action card based on reality under system's non-failure conditions.
3. selective system median generatrix A is the load growth zone, and growth rate is that per second increases by 1% under the normal mode, and be 800 seconds rise time.
4. simulation calculation
Through the whole process simulation calculation procedure, 3 machines, 9 node systems obtain the whole network node load and change in voltage curve when calculating non-fault under bus A load sustainable growth situation.(as shown in Figure 3)
5. analytical calculation
As can beappreciated from fig. 3, normal condition Down Highway A load is 125MW (P
0), under load sustainable growth situation, power limit is 355MW (P
1), through formula
Calculate dynamic electric voltage nargin index:
K
p=64.8%。
Dynamically the calculating of nargin index is more convenient, and clear concept has obtained using in actual electric network and checking.
The present invention has been described according to the preferred example of implementing.Obviously, reading and understanding above-mentioned detailed description postscript and can make multiple correction and replacement.What this invention is intended to is that the application is built into all these corrections and the replacement that has comprised within the scope that falls into appended claims book or its equivalent.
Claims (1)
1. the appraisal procedure of a weighting dynamic voltage stability margin index, it mainly may further comprise the steps:
(1) builds network topology structure according to computation requirement, set up the load model that system-computed needs, comprise static load model and induction motor model;
Setting up system's static load model is:
Wherein: P
0Be burden with power, P
1, P
2, P
3, P
4, P
5Be respectively constant impedance burden with power ratio, steady current burden with power ratio, firm power burden with power ratio, relevant burden with power ratio with frequency, the relevant burden with power ratio of voltage index, V is system's virtual voltage, V
0Be system reference voltage, Δ f is a frequency variation, L
DPBe the meritorious percent change that change of frequency 1% causes, N
PBe the active voltage index, load or burden without work model formation parameter is named in like manner, in the following formula: P
1+ P
2+ P
3+ P
4+ P
5=1, Q
1+ Q
2+ Q
3+ Q
4+ Q
5=1;
Set up system's motor load model:
According to known motor stator resistance R
S, motor stator reactance X
S, excitatory reactance X
M, rotor resistance R
R, rotor reactance X
R, motor slippage S, calculating machine moment coefficient: A, B, C;
(2) calculate preparation
Confirm fault type, fault-time, fill in corresponding fault card; According to system's actual conditions, fill in corresponding ULTC (ULTC) card, generator overexcitation restriction card and relay protection action card, accomplish data preliminary work;
(3) import according to calculation requirement: the zone of the speed of load growth, the time of load growth, load growth or bus title obtain:
Wherein: P is the load growth total amount, and Δ P is a unit interval internal loading increment, and Δ t is the unit interval, T
ENDBe the load growth duration,
According to above-mentioned formula construction load sustainable growth model, get final product through the proportionate relationship of controlling Δ P and Δ t;
(4) simulation calculation
A. read in data and network parameter, inspection Rational Parameters property, the model corresponding relation of formation state variable, intermediate variable and parameter constant and each element carries out initialization process to state variable;
B. form admittance matrix and Jacobi matrix, differentiate the rationality of starting condition;
C. revise the Jacobi matrix element of network equation according to the fault card of having filled in, confirm system state;
D. the GEAR method in the application background technology, calculate the prediction term of GEAR method:
z
m+1 (0)=Pz
m,
Wherein, P is a Pascal triangle battle array;
Calculate the correction term in the GEAR method:
z
m+1 (i+1)=z
m+1 (i)+LG
m+1 (i)
Wherein: L is a K dimension constant vector, and m, i are positive integer,
According to the correction term that obtains, with the differential equation and algebraic equation simultaneous solution;
E. carry out the variable step processing based on calculating actual conditions;
F. simulation calculation finishes, and exports required result;
(5) interpretation of result
Obtain the laggard line data analysis of simulation result, confirm the system power limits to growth from the node power growth chart, according to the node power limit that calculates and system's normal condition lower node initial power, by the computes dynamic voltage stability margin index:
K
p-dynamic voltage stability margin index (%)
P
0-system initial power (MW)
P
1-system limits power (MW)
P
0Can obtain P by system's normal operating mode statistics
1Result by simulation analysis obtains by whole process simulation calculating back.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102224459A CN101387679B (en) | 2008-09-17 | 2008-09-17 | Method for weighting dynamic voltage stability margin index |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102224459A CN101387679B (en) | 2008-09-17 | 2008-09-17 | Method for weighting dynamic voltage stability margin index |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101387679A CN101387679A (en) | 2009-03-18 |
CN101387679B true CN101387679B (en) | 2012-05-23 |
Family
ID=40477215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102224459A Active CN101387679B (en) | 2008-09-17 | 2008-09-17 | Method for weighting dynamic voltage stability margin index |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101387679B (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336882B (en) * | 2013-05-23 | 2017-12-26 | 国家电网公司 | A kind of overall process dynamic voltage stability margin appraisal procedure based on time-domain-simulation |
CN103606921B (en) * | 2013-12-10 | 2015-08-19 | 上海交通大学 | Based on the Transient Voltage Stability determination methods of real-time broad sense Thevenin's equivalence |
CN105743071B (en) * | 2014-12-12 | 2018-07-20 | 国家电网公司 | A kind of relay protection safety on line nargin assessment system |
CN106099908B (en) * | 2016-05-20 | 2022-02-01 | 中国电力科学研究院 | Method for evaluating stability of medium-and-long-term voltage of receiving-end power grid |
CN107767003B (en) * | 2016-08-18 | 2022-01-18 | 中国电力科学研究院 | Receiving-end power grid stability assessment method |
CN109873454B (en) * | 2017-12-05 | 2023-09-22 | 中国电力科学研究院有限公司 | Method and device for evaluating group voltage supporting capability in receiving end power grid |
CN109980693A (en) * | 2017-12-27 | 2019-07-05 | 广东电网有限责任公司电力调度控制中心 | A kind of power grid security margin index calculation method |
CN109492851B (en) * | 2018-09-06 | 2021-11-30 | 国网浙江省电力有限公司经济技术研究院 | Grid frame margin evaluation method based on load growth uncertainty of different regions |
CN109167374B (en) * | 2018-09-27 | 2021-09-03 | 广东电网有限责任公司 | Distributed power grid voltage stability control method, device, equipment and storage medium |
CN111276970B (en) * | 2020-02-28 | 2023-04-14 | 国网河南省电力公司 | Method and system for determining voltage stability margin based on load distribution control |
US11371573B2 (en) * | 2020-03-09 | 2022-06-28 | GM Global Technology Operations LLC | Methods and systems for EMI assessment for brake pad wear estimation |
CN114966275A (en) * | 2022-05-20 | 2022-08-30 | 国网江苏省电力有限公司苏州供电分公司 | Non-invasive resident load monitoring method and device |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1835336A (en) * | 2006-03-10 | 2006-09-20 | 天津大学 | Power system voltage stable on-line monitoring and prevention control method based on probability |
-
2008
- 2008-09-17 CN CN2008102224459A patent/CN101387679B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1835336A (en) * | 2006-03-10 | 2006-09-20 | 天津大学 | Power system voltage stable on-line monitoring and prevention control method based on probability |
Non-Patent Citations (1)
Title |
---|
马世英等.线性模拟法求解电力***静态电压稳定裕度.《电网技术》.1998,第22卷(第9期),31-34. * |
Also Published As
Publication number | Publication date |
---|---|
CN101387679A (en) | 2009-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101387679B (en) | Method for weighting dynamic voltage stability margin index | |
CN103336882A (en) | Evaluation method of overall process dynamic voltage stabilization margin based on time domain simulation | |
Pekarek et al. | An efficient multirate simulation technique for power-electronic-based systems | |
Sung et al. | Study on optimal location of a resistive SFCL applied to an electric power grid | |
CN103036230A (en) | Dynamic equivalence method of alternating-current-direct-current serial-parallel large power system based on engineering application | |
CN104009475B (en) | Non-topological static equivalence method taking equivalent element comprehensiveness and parameter physical restriction into consideration | |
CN104201674A (en) | Comprehensive load model modeling method considering load low voltage release features | |
CN101232183B (en) | Method for evaluating limit transmitted power applied in electric network on-line scheduling decision | |
CN101192754A (en) | A realization method for integrated load model of distribution network in power system | |
CN105048468A (en) | Power transmission-distribution network integrating voltage stability assessment method based on distributed calculation | |
CN102510060B (en) | Computing method of frequency characteristic coefficient of electric system | |
Adibi et al. | Solution methods for transient and dynamic stability | |
Kopse et al. | Applying a wide-area measurement system to validate the dynamic model of a part of European power system | |
Overbye et al. | Preliminary consideration of voltage stability impacts of geomagnetically induced currents | |
CN101593976B (en) | Method for constructing physical model of double circuit transmission lines on same tower | |
Plumier | Co-simulation of electromagnetic transients and phasor models of electric power systems | |
CN105529736B (en) | Crow bar puts into double-fed fan motor group of planes fault current Equivalent Model computational methods | |
CN102280884B (en) | Power grid equivalence method | |
CN113884748A (en) | Alternating current transmission line overvoltage measuring method and device based on controllable phase shifter | |
CN106849175A (en) | A kind of double-fed fan motor unit crow bar resistance setting method | |
CN108471121A (en) | A kind of optimization method that induction conductivity unstable equilibrium point solves | |
Gordon | Impact of load behavior on transient stability and power transfer limitations | |
Abdalla et al. | Dynamic load modelling and aggregation in power system simulation studies | |
Li et al. | A method for frequency dynamics analysis and load shedding assessment based on the trajectory of power system simulation | |
JP4584157B2 (en) | Method, apparatus and program for evaluating voltage stability in power system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |