CN107147114B - A kind of method for visualizing for realizing power system steady state voltage stability domain - Google Patents
A kind of method for visualizing for realizing power system steady state voltage stability domain Download PDFInfo
- Publication number
- CN107147114B CN107147114B CN201710419843.9A CN201710419843A CN107147114B CN 107147114 B CN107147114 B CN 107147114B CN 201710419843 A CN201710419843 A CN 201710419843A CN 107147114 B CN107147114 B CN 107147114B
- Authority
- CN
- China
- Prior art keywords
- load
- voltage stability
- node
- power
- power system
- 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
-
- 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]
-
- 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/002—Flicker reduction, e.g. compensation of flicker introduced by non-linear load
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of method for visualizing for realizing power system steady state voltage stability domain, are related to power system voltage stabilization technical field.The method for visualizing is by being calculated online the corresponding static voltage stability load margin of whole forecast failure situations, and it is arranged by ascending order, using the corresponding power system load nargin of forecast failure situation as a point in load level reference axis, another point on using electric system actual negative charge values as load level reference axis, then centered on load margin, its left side forms Static Voltage Stability Region, its right side forms quiescent voltage unstable domain, finally whole fault conditions are sorted by the size of voltage stabilization load margin, it is listed in the visualization Static Voltage Stability Region that forecast failure can be obtained on a figure；The visualization in power system steady state voltage stability domain may be implemented in the present invention, according to the difference of forecast failure collection, can also obtain the corresponding Static Voltage Stability Region of different forecast failure collections.
Description
Technical field
The invention belongs to power system voltage stabilization technical field more particularly to a kind of realization power system static voltage are steady
The method for visualizing of localization.
Background technique
The theoretical basis of static electric voltage stability analysis is the performable theory of power flow equation, usually with the function of electric power networks
Rate delivery limits are as the index for judging power system steady state voltage stability.Wu Zhengqiu, Li Bo, Zhong Hao, Zeng Xingjia, Li Lian
The big power system steady state voltage stability limit and nargin, which calculate, summarizes [J] Power System and its Automation journal, and 2010,22
(1): 126-132, to common static voltage stabilities such as load margin, Jacobian matrix singular value, sensitivity index, modulus of impedance
Index is reviewed, and it is normal to compared Continuation Method, direct method, Nonlinear Programming Method, singular value decomposition method, Sensitivity Method etc.
The advantage and disadvantage of voltage stability index calculation method.
Static Voltage Stability Region is a kind of effective tool for describing power system steady state voltage stability boundary, existing static state
Voltage stability domain construction method is all based on the coordinate system of electric system carrying out practically parameter building.NiuBen, Yu Yixin, Jia Hong
Outstanding person, Yang Yanbin, He Nanqiang, Tang Zhiyu, Zhang Yiming, three-dimensional visualization technique [J] power grid of Fu Hongjun Static Voltage Stability Region
Technology, 2005,29 (7): the second approximation of 56-59, Wang Gang, Zhang Xuemin, Mei Shengwei Static Voltage Stability Region Boundary analyzes [J]
Proceedings of the CSEE, 2008,28 (19): 30-35, waiting researchers to define the boundary of Static Voltage Stability Region is power flow equation
Singular point obtain the secondary of static voltage stability boundary in route space using the power flow equation at electric system singular point
Approximate equation realizes the visualization of stable region to obtain the geometry description of a kind of Static Voltage Stability Region Boundary.
Existing quiescent voltage domain construction method is all based on the coordinate system of electric system carrying out practically parameter building, usually
Complex, power system dispatcher is difficult to intuitively grasp the static voltage stability state of electric system.Therefore, such as
It is problem to be solved that, which designs a Static Voltage Stability Region that is more succinct, being easy to understand and realize,.
Summary of the invention
In view of the deficiencies of the prior art, the present invention provides a kind of visualization side for realizing power system steady state voltage stability domain
Method, can it is more succinct, be easy to understand depict power system steady state voltage stability domain.
The present invention is to solve above-mentioned technical problem by the following technical solutions: a kind of realization power system static electricity
The method for visualizing of voltage stabilization field, including the following steps:
Step 1, the operating parameter for acquiring electric system construct electric system according to the demand that Operation of Electric Systems is analyzed
N-n forecast failure set SC；
Step 2, based on electric system current operating status and parameter, consider forecast failure set S one by oneCIn whole
Forecast failure situation establishes the optimization mould of power system steady state voltage stability load margin calculating in the case of m-th of forecast failure
Type PmAre as follows:
Pm:min-λm
s.t.f(x,λm)=0
g(x)≤0
Wherein, λmIt is the corresponding power system load nargin of m-th of forecast failure situation, m ∈ SC, f (x, λm)=0 is to contain
The augmentation power flow equation of load margin, g (x)≤0 are the physical limits of power system steady stability operation constraint and power equipment
Constraint；X is the variable of electric system augmentation power flow equation；
Step 3, using optimization algorithm calculation optimization model P one by onem, obtain all pre- in electric system forecast failure set
Think the corresponding load margin λ of fault conditionm, m ∈ SC；Judge whether to complete the corresponding load margin λ of whole forecast failure situationsm
Calculating, if so, enter in next step；If it is not, then continuing calculated load nargin λm；
Step 4 establishes a two-dimensional coordinates, and abscissa is load level, and ordinate is forecast failure situation；
Step 5 arranges load margin λ by ascending orderm, m ∈ SC, corresponding with forecast failure situation in the two-dimensional coordinates
Power system load nargin λmAs a point in load level reference axis；
Step 6, centered on load margin, left side formed Static Voltage Stability Region, right side formed quiescent voltage not
Stable region；
The visualized graphs of step 7, output power static system voltage stability domain and unstable domain.
Further, in the step 1 electric system N-n forecast failure set SCThe most commonly used is N-1 or N-2 anticipation events
Hinder set SC。
Further, the optimization algorithm in the step 3 uses original dual interior point, has good fast convergence
Property and robustness, be very suitable for solve large-scale nonlinear constrained minimization problem.
Further, the solution procedure of the original dual interior point includes: 1. to introduce slack variable, will be in step 2
Inequality constraints be converted to equality constraint；2. the Optimized model without inequality constraints is configured to Lagrangian；③
According to Lagrangian, disturbance KKT conditional equation group is derived；4. disturbing KKT conditional equation using Newton Algorithm, can obtain
Obtain the solution of Optimized model.
Further, the method for visualizing for realizing power system steady state voltage stability domain uses MATLAB Programming with Pascal Language
To realize.
Compared with prior art, the method for visualizing provided by the present invention for realizing power system steady state voltage stability domain,
By sorting to the corresponding static voltage stability load margin of forecast failure situations whole in forecast failure set, to obtain quiet
State voltage stability boundary, the method for visualizing are easily achieved, can be used for electric system in line computation；The present invention obtains quiet
State voltage stability domain is more succinct, intuitive, understands convenient for power system dispatcher and uses.
Detailed description of the invention
It, below will be to attached drawing needed in embodiment description in order to illustrate more clearly of technical solution of the present invention
It is briefly described, it should be apparent that, the accompanying drawings in the following description is only one embodiment of the present of invention, general for this field
For logical technical staff, without creative efforts, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of flow chart for the method for visualizing for realizing power system steady state voltage stability domain of the present invention；
Fig. 2 is that IEEE-30 node power static system voltage stability domain figure of the present invention is shown.
Specific embodiment
With reference to the attached drawing in the embodiment of the present invention, the technical solution in the present invention is clearly and completely described,
Obviously, described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based in the present invention
Embodiment, those of ordinary skill in the art's every other embodiment obtained without creative labor,
It shall fall within the protection scope of the present invention.
As shown in Figure 1, a kind of method for visualizing for realizing power system steady state voltage stability domain provided by the present invention is adopted
It is realized with MATLAB Programming with Pascal Language, by taking IEEE-30 node power static system voltage stabilization domain analysis as an example, including it is following
Step:
Step 1, the operating parameter for acquiring electric system, according to the demand that Operation of Electric Systems is analyzed, building includes 10
The electric system forecast failure set S of N-1 forecast failure situationC, specific as shown in table 1；
1 IEEE-30 system N-1 forecast failure collection of table
| Failure number | Headend node | Endpoint node |
| Branch 1 | 1 | 2 |
| Branch 2 | 1 | 3 |
| Branch 3 | 2 | 4 |
| Branch 4 | 3 | 4 |
| Branch 5 | 2 | 5 |
| Branch 6 | 2 | 6 |
| Branch 7 | 4 | 6 |
| Branch 8 | 5 | 7 |
| Branch 9 | 6 | 7 |
| Branch 10 | 6 | 8 |
Step 2, based on electric system current operating status and parameter, consider forecast failure set S one by oneCIn whole
Forecast failure situation establishes the optimization mould of power system steady state voltage stability load margin calculating in the case of m-th of forecast failure
Type PmAre as follows:
Pm:min-λm
Wherein, λmIt is the corresponding power system load nargin of m-th of forecast failure situation, m ∈ SC；SBFor node set；SG
For generator node set；SRFor reactive source set；PGiThe active power issued for node i generator；QRiFor all kinds of nothings of node i
The reactive power that function source issues；PLiAnd QLiRespectively node i load is active and reactive power；bPiAnd bQiRespectively node i is negative
Lotus increases direction；ViAnd δiThe respectively voltage magnitude and phase angle of node i；YijFor node admittance matrix element；δij=δi-δj-
αij, αijFor node admittance matrix respective element phase angle； P GiRespectively PGiCorresponding bound； Q RiRespectively QRiIt is right
Answer bound； V iIt is node i voltage magnitude bound respectively；The variable x of electric system augmentation power flow equation includes generator
The active-power P of sendingGi, reactive source issue reactive power QRi, node load active-power PLi, node load reactive power
QLi, node voltage amplitude ViWith node voltage phase angle δi；
Step 3, using original dual interior point calculation optimization model P one by onem, obtain IEEE-30 node power system
The corresponding load margin λ of whole forecast failure situation in forecast failure setm, as shown in table 2；Judge whether to complete whole anticipations
The corresponding load margin λ of fault conditionmCalculating, if so, enter in next step；If it is not, then continuing calculated load nargin
λm；
2 IEEE-30 node power system forecast failure collection corresponding load nargin of table
| Failure number | Load margin |
| Branch 1 | 1.6926 |
| Branch 2 | 1.6915 |
| Branch 3 | 1.7324 |
| Branch 4 | 1.6999 |
| Branch 5 | 1.7326 |
| Branch 6 | 1.6974 |
| Branch 7 | 1.7101 |
| Branch 8 | 1.7273 |
| Branch 9 | 1.7922 |
| Branch 10 | 1.7529 |
Step 4 establishes a two-dimensional coordinates, and abscissa is load level, and ordinate is forecast failure situation；
Step 5 arranges load margin λ by ascending orderm, in the two-dimensional coordinates, with the corresponding electricity of forecast failure situation
Force system load margin λmAs a point in load level reference axis；
Step 6, centered on load margin, left side formed Static Voltage Stability Region, right side formed quiescent voltage not
Stable region, as shown in Figure 2；
The visualized graphs of step 7, output power static system voltage stability domain and unstable domain, as shown in Fig. 2, from figure
In as can be seen that the present invention can be succinct depict power system steady state voltage stability domain.
Above disclosed is only a specific embodiment of the invention, but scope of protection of the present invention is not limited thereto,
Anyone skilled in the art in the technical scope disclosed by the present invention, can readily occur in variation or modification,
It is covered by the protection scope of the present invention.
Claims (7)
1. a kind of method for visualizing for realizing power system steady state voltage stability domain, which is characterized in that including the following steps:
Step 1, the operating parameter for acquiring electric system construct the N-n of electric system according to the demand that Operation of Electric Systems is analyzed
Forecast failure set SC；
Step 2, based on electric system current operating status and parameter, consider forecast failure set S one by oneCIn whole anticipations
Fault condition establishes the Optimized model P of power system steady state voltage stability load margin calculating in the case of m-th of forecast failurem
Are as follows:
Pm:min-λm
s.t.f(x,λm)=0
g(x)≤0
Wherein, λmIt is the corresponding power system load nargin of m-th of forecast failure situation, m ∈ SC, f (x, λm)=0 is containing load
The augmentation power flow equation of nargin, g (x)≤0 are power system steady stability operation constraint and the physical limit constraint of power equipment；
X is the variable of electric system augmentation power flow equation；
Step 3, using optimization algorithm calculation optimization model P one by onem, obtain all anticipation events in electric system forecast failure set
Hinder the corresponding load margin λ of situationm, m ∈ SC；Judge whether to complete the corresponding load margin λ of whole forecast failure situationsmMeter
It calculates, if so, entering in next step；If it is not, then continuing calculated load nargin λm；
Step 4 establishes a two-dimensional coordinates, and abscissa is load level, and ordinate is forecast failure situation；
Step 5 arranges load margin λ by ascending orderm, m ∈ SC, in the two-dimensional coordinates, with the corresponding electricity of forecast failure situation
Force system load margin λmAs a point in load level reference axis；
Step 6, centered on load margin, left side formed Static Voltage Stability Region, right side formed quiescent voltage it is unstable
Domain；
The visualized graphs of step 7, output power static system voltage stability domain and unstable domain.
2. realizing the method for visualizing in power system steady state voltage stability domain as described in claim 1, which is characterized in that described
The N-n forecast failure set S of electric system in step 1CThe most commonly used is N-1 or N-2 forecast failure set SC。
3. realizing the method for visualizing in power system steady state voltage stability domain as described in claim 1, which is characterized in that described
The variable x of electric system augmentation power flow equation includes the active-power P that generator issues in step 2Gi, reactive source issue it is idle
Power QRi, node load active-power PLi, node load reactive power QLi, node voltage amplitude ViWith node voltage phase angle δi。
4. realizing the method for visualizing in power system steady state voltage stability domain as claimed in claim 3, which is characterized in that described
The Optimized model P of load margin is calculated in step 2mAre as follows:
Pm:min-λm
Wherein, λmIt is the corresponding power system load nargin of m-th of forecast failure situation, m ∈ SC；SBFor node set；SGFor hair
Motor node set；SRFor reactive source set；PGiThe active power issued for node i generator；QRiFor all kinds of reactive sources of node i
The reactive power of sending；PLiAnd QLiRespectively node i load is active and reactive power；bPiAnd bQiRespectively the load of node i increases
Add direction；ViAnd δiThe respectively voltage magnitude and phase angle of node i；YijFor node admittance matrix element；δij=δi-δj-αij, αij
For node admittance matrix respective element phase angle； P GiRespectively PGiCorresponding bound； Q RiRespectively QRiAbove and below corresponding
Limit； V iIt is node i voltage magnitude bound respectively；VjFor the voltage magnitude of node j.
5. realizing the method for visualizing in power system steady state voltage stability domain as described in claim 1, which is characterized in that described
Optimization algorithm in step 3 uses original dual interior point.
6. realizing the method for visualizing in power system steady state voltage stability domain as claimed in claim 5, which is characterized in that described
The solution procedure of original dual interior point includes: 1. to introduce slack variable, and the inequality constraints in the step 2 is converted to
Equality constraint；2. the Optimized model without inequality constraints is configured to Lagrangian；3. being pushed away according to Lagrangian
Lead disturbance KKT conditional equation group；4. disturbing KKT conditional equation using Newton Algorithm, the solution of Optimized model can be obtained.
7. realizing the method for visualizing in power system steady state voltage stability domain as described in claim 1, which is characterized in that use
MATLAB Programming with Pascal Language is realized.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710419843.9A CN107147114B (en) | 2017-06-06 | 2017-06-06 | A kind of method for visualizing for realizing power system steady state voltage stability domain |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710419843.9A CN107147114B (en) | 2017-06-06 | 2017-06-06 | A kind of method for visualizing for realizing power system steady state voltage stability domain |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107147114A CN107147114A (en) | 2017-09-08 |
| CN107147114B true CN107147114B (en) | 2019-08-13 |
Family
ID=59779622
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710419843.9A Active CN107147114B (en) | 2017-06-06 | 2017-06-06 | A kind of method for visualizing for realizing power system steady state voltage stability domain |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107147114B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103580022A (en) * | 2013-11-05 | 2014-02-12 | 国家电网公司 | Electrical power system dynamic reactive storage computing method |
| CN103795058A (en) * | 2014-02-12 | 2014-05-14 | 国家电网公司 | Static voltage stability margin analyzing and system fault ordering method of power system |
| CN104734151A (en) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | Electric system static security assessment method |
| CN105406480A (en) * | 2015-12-31 | 2016-03-16 | 国网冀北电力有限公司电力科学研究院 | Voltage stability preventive control method and device |
| CN106558879A (en) * | 2015-09-30 | 2017-04-05 | 中国电力科学研究院 | A kind of electricity net safety stable analysis method based on safety and stability index |
-
2017
- 2017-06-06 CN CN201710419843.9A patent/CN107147114B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103580022A (en) * | 2013-11-05 | 2014-02-12 | 国家电网公司 | Electrical power system dynamic reactive storage computing method |
| CN103795058A (en) * | 2014-02-12 | 2014-05-14 | 国家电网公司 | Static voltage stability margin analyzing and system fault ordering method of power system |
| CN104734151A (en) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | Electric system static security assessment method |
| CN106558879A (en) * | 2015-09-30 | 2017-04-05 | 中国电力科学研究院 | A kind of electricity net safety stable analysis method based on safety and stability index |
| CN105406480A (en) * | 2015-12-31 | 2016-03-16 | 国网冀北电力有限公司电力科学研究院 | Voltage stability preventive control method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107147114A (en) | 2017-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hans et al. | Hierarchical distributed model predictive control of interconnected microgrids | |
| Zhou et al. | Data-adaptive robust unit commitment in the hybrid AC/DC power system | |
| US8606420B2 (en) | Solving large-scale security-constrained economic dispatch problem in real-time | |
| Kiani et al. | Wholesale energy market in a smart grid: Dynamic modeling and stability | |
| Razmi et al. | A literature review of the control challenges of distributed energy resources based on microgrids (MGs): past, present and future | |
| Abido | Multiobjective optimal VAR dispatch using strength pareto evolutionary algorithm | |
| Wang et al. | Distributed optimal power flow of DC microgrids: A penalty based ADMM approach | |
| CN104573906A (en) | System and method for analyzing oscillatory stability in electrical power transmission systems | |
| CA2838724A1 (en) | Direct current transmission and distribution system and method of operating the same | |
| Ilyushin et al. | State-of-the-Art Literature Review of Power Flow Control Methods for Low-Voltage AC and AC-DC Microgrids | |
| JP2017028915A (en) | Arrangement calculation apparatus for power storage device, arrangement calculation method for power storage device, arrangement calculation system for power storage device, and program | |
| Zhang et al. | Frequency-constrained unit commitment for power systems with high renewable energy penetration | |
| CN113610432A (en) | Data acquisition analysis formula photovoltaic station accuse system based on BS version | |
| Cheng et al. | Distributed economic power dispatch and bus voltage control for droop-controlled DC microgrids | |
| CN101534010B (en) | Method for solving the local boundary of cut-set voltage stability regions on the basis of perturbation | |
| WO2020041903A1 (en) | Nano/micro-grid power distribution system | |
| CN107147114B (en) | A kind of method for visualizing for realizing power system steady state voltage stability domain | |
| CN106483393A (en) | Collocation method and system are unified in harmonic wave and voltage dip monitoring point | |
| CN115136438A (en) | Distributed resource management device and distributed resource management method | |
| CN117977704A (en) | Method and system for constructing source network charge storage integrated collaborative system control architecture | |
| CN109063253B (en) | Reliability modeling method for aviation power supply system based on Bayesian network | |
| CN110649712A (en) | Urban power grid operation visual monitoring system based on electricity consumption information acquisition system | |
| Yang et al. | Emergency Power Supply Restoration Strategy for Distribution Network Considering Support of Microgrids with High-Dimensional Dynamic Correlations | |
| Vilaisarn et al. | Linear three-phase power flow model for autonomous microgrid systems | |
| Li et al. | Distributed Stochastic Scheduling of Massive Backup Batteries in Cellular Networks for Operational Reserve and Frequency Support Ancillary Services |
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 |