CN108899897A - The power system voltage weak node on-line identification method of data-driven - Google Patents
The power system voltage weak node on-line identification method of data-driven Download PDFInfo
- Publication number
- CN108899897A CN108899897A CN201810764027.6A CN201810764027A CN108899897A CN 108899897 A CN108899897 A CN 108899897A CN 201810764027 A CN201810764027 A CN 201810764027A CN 108899897 A CN108899897 A CN 108899897A
- Authority
- CN
- China
- Prior art keywords
- node
- voltage
- generator
- generator node
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
The power system voltage weak node on-line identification method of data-driven, is related to electric system online steady stability Situation Assessment field.The present invention is the problem in order to solve the existing voltage stabilization weak node discrimination method based on Thevenin's equivalence parameter and impedance matching principle, applied to nonlinear power system accuracy difference.The voltage and current information of each generator node is obtained according to Real-Time Power System Operation States, it is a generator node by each generator node aggregation in electric system, the net power output of generator node after being polymerize, the generator node potential after being polymerize according to the voltage and current information and the net power output；Using the voltage and current information of some load bus in the potential, electric system, equivalent impedance is obtained；Using the ratio of equivalent impedance and the load impedance of corresponding load bus as the load impedance of the load bus than index, the identification of voltage weak node is realized.For recognizing voltage weak node.
Description
Technical field
The present invention relates to the online recognitions of electric system weak node, belong to the online steady stability Situation Assessment of electric system
Field.
Background technique
The power supply of electric system becomes more diverse with load, and power generation characteristics and part throttle characteristics are to power network safety operation
Influence receive significant attention.By impedance matching principle by Thevenin's equivalence in conjunction with Network Voltage Stability, electricity can be improved
The identification efficiency of net Voltage Instability.So far, Thevenin's equivalence parameter identification method is all based on various assumed conditions, has one
Fixed limitation.Moreover, proper restraint of some discrimination method due to not accounting for parameter to be identified itself, causes to ask
There are severe deviations in the result obtained.Therefore, seek that a kind of accuracy is high, adaptable weak node discrimination method seems especially
It is important.
Although there is a large amount of literature research Thevenin's equivalence parameter identification method at present, the time-varying of bulk power grid is all received
Property and nonlinear restriction.So far, existing Thevenin's equivalence parameter identification method mainly has following three classes：
(1) based on the Thevenin's equivalence parameter identification method that more moment locally measure, this kind of methods need equal value nodes
The measured data of discontinuity surface when locating more.By assuming that in sampling time interval, to the Dai Weinan parameter of system in terms of equal value nodes
It does not change, calculates Thevenin's equivalence parameter using the method for curve matching.But such method is the problem is that system is non-
The premise for measuring time interval Nei Daiweinan parameter constant is unable to satisfy when linear stronger.
(2) the networked examination method measured based on single moment the whole network, all nodes of discontinuity surface when this kind of methods need single
Measured data.Although overcoming the parameter drift problem in the method for measurement of more moment parts, such method there are the problem of
It is to change the characteristic of non-linear element or without specific coupling terms processing mode.
(3) the Thevenin's equivalence parameter identification method for locally being measured based on single moment and virtually being disturbed, this kind of methods are needed
The non-of system part is obtained after the practical metric data of discontinuity surface system part when wanting list, with the virtual disturbance of addition practical measures number
According to matching, computing system Thevenin's equivalence parameter.This kind of methods are the problem is that having initial value dependence or having implicit
Premise.
Summary of the invention
The present invention is existing based on the voltage stabilization of Thevenin's equivalence parameter and impedance matching principle weakness in order to solve
Node authentication method, the problem applied to nonlinear power system accuracy difference.The power system voltage of data-driven is now provided
Weak node on-line identification method.
The power system voltage weak node on-line identification method of data-driven, the described method comprises the following steps：
Step 1: the voltage and current information of each generator node is obtained according to Real-Time Power System Operation States, it will
Each generator node aggregation in electric system is a generator node, the net output of the generator node after being polymerize
Power, the generator node potential after being polymerize according to the voltage and current information and the net power output；
Step 2: using after polymerization generator node potential, in electric system some load bus voltage and current
Information obtains equivalent impedance；
Step 3: being hindered the ratio of equivalent impedance and the load impedance of corresponding load bus as the load of the load bus
It is anti-that the identification of voltage weak node in electric system is realized according to impedance ratio index L than index L.
Beneficial effects of the present invention are：
The application is changed system generator node and load bus classification polymerization by abbreviation power system network respectively
Letter obtains the classification polymerization equivalent network that duty value branch is linked in source.By calculate enter in terms of some node etc.
It is worth the modulus of impedance ratio of impedance and the load impedance of the node, identifies power system voltage weak node.The application is worn with existing
The impedance modular character for tieing up southern equivalent parameters discrimination method is compared, the metric data of discontinuity surface when this method is based on list, it is only necessary to be
System generator node metric data, does not need network parameter.Also, the application is demonstrated by being applied to ieee standard system
The validity of the power system voltage weak node on-line identification method of data-driven.
Detailed description of the invention
Fig. 1 is original power network；
Fig. 2 is classification polymerization star equivalent network；
Fig. 3 is classification polymerization triangle equivalent network；
Fig. 4 is the abbreviation of classification polymerization triangle equivalent network；
Fig. 5 is equivalence PQ node connection in the classification polymerization equivalent network of source；
Fig. 6 is the L index of 30 node system of IEEE；
Fig. 7 is the L index of 57 node system of IEEE.
Specific embodiment
Specific embodiment one：Present embodiment is illustrated referring to figs. 1 to Fig. 5, data described in present embodiment are driven
Dynamic power system voltage weak node on-line identification method, the described method comprises the following steps：
Step 1: the voltage and current information of each generator node is obtained according to Real-Time Power System Operation States, it will
Each generator node aggregation in electric system is a generator node, the net output of the generator node after being polymerize
Power, the generator node potential after being polymerize according to the voltage and current information and the net power output；
Step 2: using after polymerization generator node potential, in electric system some load bus voltage and current
Information obtains equivalent impedance；
Step 3: being hindered the ratio of equivalent impedance and the load impedance of corresponding load bus as the load of the load bus
It is anti-that the identification of voltage weak node in electric system is realized according to impedance ratio index L than index L.
In present embodiment, abbreviation is carried out for original power network shown in FIG. 1, by generator node and load bus point
Type of Collective, abbreviation, obtains classification polymerization star equivalent network shown in Fig. 2 respectively.In Fig. 1, G1 to Gm indicates each generator
Node, SLnIndicate the net consumption power of n-th of load bus, PLnIndicate the net consumption active power of n-th of load bus, QLn
Indicate that n-th of load bus consumes reactive power only, j indicates imaginary part；SLiIndicate that each load bus consumes power only,
PLiIndicate that each load bus consumes active power only, i indicates studied load bus.Utilize the voltage and electricity of generator node
Flow measurement information calculates the electrical potential parameters of classification polymerization equivalent network equivalence PV (generator) node.
In step 2, star triangular transformation, decoupling equivalence PV (power generation are carried out to the classification polymerization triangle equivalent network of Fig. 2
Machine) and equivalence PQ node, the classification polymerization star equivalent network of Fig. 3 is obtained, by Fig. 3 medium value PV node to VXImpedance tear open
Point, as shown in figure 4, classification polymerization equivalent network of equivalent PQ (load) node connection in source of Fig. 5 can be obtained, in Fig. 4
VpqeqIndicate the equivalent potential of other load bus.
Specific embodiment two：Present embodiment is the electricity of the electric system to data-driven described in specific embodiment one
Ironed weak bus on-line identification method is described further, in present embodiment, in step 1, according to the voltage and current
The process of generator node potential after information and the net power output are polymerize is：
According to formula：
Obtain all generator node net power outputs and SG,
In formula, m is the number of generator node, SGFor the sum of all generator node net power outputs, SGkFor each hair
Motor node generated output power, SGLkPower, S are consumed for each generator node institute's on-loadCkFor each generator node
Susceptance output power over the ground；
The sum that all generator nodes export electric current only is：
In formula, IGExport the sum of electric current only for all generator nodes,For the conjugation of each generator node voltage,For the conjugation of each generator node susceptance output power over the ground,Function is consumed for each generator node institute's on-load
The conjugation of rate,For the conjugation of each generator node susceptance output power over the ground；
Generator node potential after polymerization is：
In formula,Export the conjugation of the sum of electric current, V only for all generator nodespveqFor the equivalence electricity of equivalent PV node
Gesture.
Specific embodiment three：Present embodiment is the electricity of the electric system to data-driven described in specific embodiment two
Ironed weak bus on-line identification method is described further, and in present embodiment, in step 2, utilizes the generator section after polymerization
The voltage and current information of some load bus in point potential, electric system, obtaining equivalent impedance is：
In formula, ZpveqFor the equivalent impedance of equivalent PV node to some load bus, VLFor the voltage of some load bus,
ILFor the electric current of some load bus.
In present embodiment, for the impedance parameter for seeking classification polymerization equivalent network, classification polymerization equivalent network is carried out
Abbreviation, the equivalent PQ node connection for obtaining Fig. 5 polymerize equivalent network in the classification of source.Joined using the potential of equivalent PV node
Number calculates the equivalent impedance parameter entered in terms of studied load bus.
Although seeking for the equivalent impedance and Thevenin's equivalence impedance, is all in terms of load bus into system, this method
The classification of seeking being built upon of equivalent impedance polymerize on the basis of equivalent network, different from asking for traditional Thevenin's equivalence impedance
It takes, so referred to as class Thevenin's equivalence impedance.
Specific embodiment four：Present embodiment is the electricity of the electric system to data-driven described in specific embodiment three
Ironed weak bus on-line identification method is described further, negative with some according to equivalent impedance in step 3 in present embodiment
The ratio of lotus impedance obtains some load impedance than index L in real time, is according to impedance ratio index L：
In formula, ZLFor the impedance of some load bus.
It is proposed by the present invention to utilize class Thevenin's equivalence impedance ratio index shown in formula 5, Ke Yiyong in present embodiment
To describe the vulnerability of each load bus.Largely simulation results show, which can be with accurate recognition load bus
Vulnerability, calculating speed is fast, and only needs the metric data of system PV node, does not need network parameter, and it is quiet to be suitable for powernet
State stablizes Situation Assessment and decision.
Embodiment 1：
Simulation analysis is carried out by taking 30 node system of IEEE as an example, specific steps and result are as follows：
1) the real-time measurement data of each generator node are replaced using the data of Load flow calculation.It is all in simulation calculation
Load type is taken as constant output characteristic, and the power of load consumption is gradually increased with identical multiple, generated output power also with
Identical multiple does corresponding growth, until trend does not restrain.Each node voltage vector sum current vector is by MATPOWER trend
It is calculated, convergence precision 10-6。
2) electrical potential parameters of classification polymerization equivalent network equivalence PV node are calculated using formula 1 to formula 3.
3) class Thevenin's equivalence impedance parameter is calculated using formula 4.
4) the impedance ratio index of studied node is calculated using formula 5, as a result as shown in Figure 6.
It can be seen from the figure that load bus 3 is most weak node, and it can be seen that as system loading increases, respectively
The L index of load bus all increases, and when close to system crash, the L index of weak node increases very fast, the L of non-weak node
Index increases slower.
5) result above actually matches with electric system with discovery.
It is one embodiment of the present of invention below, carries out simulation analysis, specific steps by taking 57 node system of IEEE as an example
It is as follows with result：
1) the real-time measurement data of each generator node are replaced using the data of Load flow calculation.It is all in simulation calculation
Load type is taken as constant output characteristic, and the power of load consumption is gradually increased with identical multiple, generated output power also with
Identical multiple does corresponding growth, until trend does not restrain.Each node voltage vector sum current vector is by MATPOWER trend
It is calculated, convergence precision 10-6。
2) electrical potential parameters of classification polymerization equivalent network equivalence PV node are calculated using formula 1 to formula 3.
3) class Thevenin's equivalence impedance parameter is calculated using formula 4.
4) the impedance ratio index that studied node is calculated using formula 5, as a result as Fig. 7 shows.
It can be seen from the figure that load bus 17 is most weak node, and it can be seen that as system loading increases,
The L index of each load bus all increases, and when close to system crash, the L index of weak node increases comparatively fast, non-weak node
L index increases slower.
5) result above actually matches with electric system with discovery.
Claims (3)
1. the power system voltage weak node on-line identification method of data-driven, which is characterized in that the method includes following
Step：
Step 1: the voltage and current information of each generator node is obtained, by electric power according to Real-Time Power System Operation States
Each generator node aggregation in system is a generator node, the net output work of the generator node after being polymerize
Rate, the generator node potential after being polymerize according to the voltage and current information and the net power output；
Step 2: using after polymerization generator node potential, in electric system some load bus voltage and current information,
Obtain equivalent impedance；
Step 3: using the ratio of equivalent impedance and the load impedance of corresponding load bus as the load impedance ratio of the load bus
Index L realizes the identification of voltage weak node in electric system according to impedance ratio index L.
2. the power system voltage weak node on-line identification method of data-driven according to claim 1, which is characterized in that
Generator node potential in step 1, after being polymerize according to the voltage and current information and the net power output
Process be：
According to formula：
Obtain all generator node net power outputs and SG,
In formula, m is the number of generator node, SGFor the sum of all generator node net power outputs, SGkFor each generator
Node generated output power, SGLkPower, S are consumed for each generator node institute's on-loadCkOver the ground for each generator node
Susceptance output power；
The sum that all generator nodes export electric current only is：
In formula, IGExport the sum of electric current only for all generator nodes,For the conjugation of each generator node voltage,For
The conjugation of each generator node susceptance output power over the ground,Being total to for power is consumed for each generator node institute's on-load
Yoke,For the conjugation of each generator node susceptance output power over the ground；
Generator node potential after polymerization is：
In formula,Export the conjugation of the sum of electric current, V only for all generator nodespveqFor the equivalent potential of equivalent PV node.
3. the power system voltage weak node on-line identification method of data-driven according to claim 2, which is characterized in that
In step 2, using after polymerization generator node potential, in electric system some load bus voltage and current information, obtain
It is to equivalent impedance：
In formula, ZpveqFor the equivalent impedance of equivalent PV node to some load bus, VLFor the voltage of some load bus, ILFor
The electric current of some load bus.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810764027.6A CN108899897B (en) | 2018-07-12 | 2018-07-12 | Data-driven online identification method for weak voltage nodes of power system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810764027.6A CN108899897B (en) | 2018-07-12 | 2018-07-12 | Data-driven online identification method for weak voltage nodes of power system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108899897A true CN108899897A (en) | 2018-11-27 |
| CN108899897B CN108899897B (en) | 2021-05-18 |
Family
ID=64349600
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810764027.6A Expired - Fee Related CN108899897B (en) | 2018-07-12 | 2018-07-12 | Data-driven online identification method for weak voltage nodes of power system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108899897B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102157938A (en) * | 2011-04-12 | 2011-08-17 | 哈尔滨工业大学 | Method for online identifying lower-voltage stability nodes of power system |
| CN104778367A (en) * | 2015-04-20 | 2015-07-15 | 山东大学 | On-line wide-area Thevenin equivalent parameter calculation method based on single state section |
| CN107658881A (en) * | 2017-09-27 | 2018-02-02 | 国网上海市电力公司 | Voltage stability critical point determination methods based on Thevenin's equivalence method |
-
2018
- 2018-07-12 CN CN201810764027.6A patent/CN108899897B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102157938A (en) * | 2011-04-12 | 2011-08-17 | 哈尔滨工业大学 | Method for online identifying lower-voltage stability nodes of power system |
| CN104778367A (en) * | 2015-04-20 | 2015-07-15 | 山东大学 | On-line wide-area Thevenin equivalent parameter calculation method based on single state section |
| CN107658881A (en) * | 2017-09-27 | 2018-02-02 | 国网上海市电力公司 | Voltage stability critical point determination methods based on Thevenin's equivalence method |
Non-Patent Citations (2)
| Title |
|---|
| YUNFEI WANG等: "Voltage Stability Monitoring Based on the", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
| 邢晨等: "电力系统戴维南等值参数辨识方法综述", 《中国科技论文在线》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108899897B (en) | 2021-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109061506A (en) | Lithium-ion-power cell SOC estimation method based on Neural Network Optimization EKF | |
| CN102074955B (en) | Method based on knowledge discovery technology for stability assessment and control of electric system | |
| CN110492479B (en) | Method for identifying rotational inertia and damping of distributed grid-connected equipment | |
| CN104615882B (en) | Extensive time-lag power system feature value calculating method based on EIGD | |
| CN104778367B (en) | Wide area Thevenin's equivalence parameter on-line calculation method based on a single state section | |
| CN102510263A (en) | Method for identifying practical parameters of synchronous generator on basis of load rejection test and numerical difference | |
| CN108448585B (en) | A kind of electric network swim equation solution method of linearization based on data-driven | |
| CN107069741A (en) | A kind of novel linear tidal current computing method | |
| Zhao et al. | Criterion to evaluate power system online transient stability based on adjoint system energy function | |
| Sharma et al. | Data-driven identification of nonlinear power system dynamics using output-only measurements | |
| Liu et al. | Unscented particle filter for SOC estimation algorithm based on a dynamic parameter identification | |
| Murugan et al. | Determining true harmonic contributions of sources using neural network | |
| CN107516115B (en) | A kind of load model canonical parameter extracting method based on k- central point algorithm | |
| CN109802433A (en) | A kind of gird-connected inverter oscillation of power inhibition system and method | |
| CN105576651B (en) | A kind of medium and small water power group of planes hybrid parallel dynamic equivalent method | |
| CN103986400B (en) | Based on the model reference adaptive system parameters automatic setting method that two dimension fuzzy controls | |
| Huang et al. | Damping inter‐area oscillation using reinforcement learning controlled TCSC | |
| CN108899897A (en) | The power system voltage weak node on-line identification method of data-driven | |
| CN109698505A (en) | The regulation of the online prevention and control of bulk power grid static voltage stability quantifies mapping calculation method | |
| CN110703038B (en) | Harmonic impedance estimation method suitable for fan access power distribution network | |
| CN116131261A (en) | Micro-grid dynamic equivalent modeling method considering model robustness | |
| Shi et al. | An improved adaptive square root cubature kalman filter method for estimating state-of-charge of lithium-ion batteries | |
| CN113794198A (en) | Method, device, terminal and storage medium for suppressing broadband oscillation | |
| Liu et al. | Dynamic aggregation modeling for droop control inverter based on slow coherency algorithm | |
| Sang et al. | Double Inverted Pendulum control based on three-loop PID and improved BP Neural network |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190108 Address after: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin Applicant after: HARBIN INSTITUTE OF TECHNOLOGY Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: NORTHEAST BRANCH OF STATE GRID CORPORATION OF CHINA Address before: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin Applicant before: Harbin Institute of Technology |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210518 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |