CN107832959A - A kind of voltage stability assessment method for considering part throttle characteristics and power constraints - Google Patents

A kind of voltage stability assessment method for considering part throttle characteristics and power constraints Download PDF

Info

Publication number
CN107832959A
CN107832959A CN201711139795.4A CN201711139795A CN107832959A CN 107832959 A CN107832959 A CN 107832959A CN 201711139795 A CN201711139795 A CN 201711139795A CN 107832959 A CN107832959 A CN 107832959A
Authority
CN
China
Prior art keywords
mrow
mtd
msub
mtr
msup
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
Application number
CN201711139795.4A
Other languages
Chinese (zh)
Other versions
CN107832959B (en
Inventor
娄源媛
黎晓
刘崇茹
段秦刚
钱峰
罗钢
伍双喜
刘俊磊
张子泳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Original Assignee
North China Electric Power University
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University, Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd filed Critical North China Electric Power University
Priority to CN201711139795.4A priority Critical patent/CN107832959B/en
Publication of CN107832959A publication Critical patent/CN107832959A/en
Application granted granted Critical
Publication of CN107832959B publication Critical patent/CN107832959B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of voltage stability assessment method for considering part throttle characteristics and power constraints, it is characterised in that comprises the following steps：Step 1：The simplification admittance matrix Y of calculating network simplification Equivalent Model+；Step 2：Calculate the equivalent network admittance matrix Y for considering power constraintsS；Step 3：Form the load characteristic matrix L related to part throttle characteristics；Step 4：The Thevenin's equivalence parameter of computing system load, and carry out Voltage Stability Evaluation；The consideration part throttle characteristics of the present invention and the voltage stability assessment method of power constraints also can correctly assess the voltage stabilization of load in complication system.

Description

A kind of voltage stability assessment method for considering part throttle characteristics and power constraints
Technical field
The present invention relates to voltage stability assessment technology field, more particularly to the electricity for considering part throttle characteristics and power constraints Press stability assessment method.
Background technology
Since the twenties in last century, stability of power system problem is carefully ground as the important prerequisite of safe operation Study carefully, until the end of the seventies in last century to the beginning of the eighties, due to the power network characterized by interconnected network collapse of voltage occurred in succession Splitting accident, Power System Voltage Stability start to be taken seriously, and one as a stabilization of power grids Journal of Sex Research is important Sub- problem.
With PMU appearance, the voltage stability assessment of system is carried out using PMU data and Thevenin's equivalence method to be turned into An important directions in one voltage stabilization research.Domestic and foreign scholars have been done greatly in terms of Thevenin's equivalence parameter identification method Measure work.Initial document is thought, in a data window, monitoring bus load change is the principal element of grid disturbance, its He partly can approximation regard as and keep constant, therefore from the system side Thevenin's equivalence potential and equivalent impedance entered in terms of monitoring bus Keep constant.Therefore Dai Weinan can be estimated using the monitoring busbar voltage in appropriate length data window, current measurement phasor Equivalent potential and equivalent impedance.When monitoring bus load change, voltage source potential amplitude may be constant, but phase angle is due to power Change can occur to change to a certain degree, and then above-mentioned hypothesis is incorrect, have document to propose Dai Weinan of potential phase angle change etc. again It is worth model.More preferable estimation effect can be obtained when only monitoring that bus load increases.If multi-node load increases, now Dai Weinan Equivalent reactance can also change, and above-mentioned hypothesis is also untenable.There is document utilization system node voltage equation in addition analytically Angle analysis Thevenin's equivalence potential and equivalent reactance composition, it is believed that the processing to other node load coupling terms is equivalent The key of model, but do not propose specific processing method.At present, what is largely used in Voltage Stability Evaluation is still that tradition is worn Tie up southern Equivalent Model.The above method has used the vacation of Thevenin's equivalence parameter constant in adjacent region data or data window length If document points out that this hypothesis is difficult to set up in actual parameter estimation, this is that the above method can not obtain preferable parameter Estimation The basic reason of effect.Later it has been proposed that total differential method, the equivalence at enough algebraic equation solvings each moment is constructed Parameter.The deficiency of this method is extremely to rely on initial value, and when grid disturbance is larger, larger mistake also occurs in estimated result Difference.
It is desirable to have a kind of voltage stability assessment method for considering part throttle characteristics and power constraints to solve in the prior art The problem of existing.
The content of the invention
It is an object of the invention to provide the voltage stability assessment method of a kind of consideration part throttle characteristics and power constraints multiple Also the voltage stabilization of load can be correctly assessed in miscellaneous system.
The present invention provides a kind of voltage stability assessment method for considering part throttle characteristics and power constraints, comprises the following steps：
Step 1：The simplification admittance matrix Y of calculating network simplification Equivalent Model+；
Step 2：Calculate the equivalent network admittance matrix Y for considering power constraintsS；
Step 3：Form the load characteristic matrix L related to part throttle characteristics；
Step 4：The Thevenin's equivalence parameter of computing system load, and carry out Voltage Stability Evaluation.
Preferably, comprise the following steps in the step 1：
(1) the topological structure of electric parameter measured according to wide area measurement system establishes the admittance matrix of network, the section of network Point voltage equation is as follows：
(2) node in network is divided into power supply node, load bus and connecting node, V in step (1)L,VT,VGRespectively Load bus voltage, connecting node voltage and the power supply node voltage for having power to inject are represented, wherein power supply node saves including PV Point and balance nodes；
(3) connecting node is eliminated from the admittance matrix of step (1) and is simplified admittance matrix：
Wherein, Y+ LL=YLL-YLT(YTT)-1YTL、Y+ LG=YLG-YLT(YLT)-1YTG、Y+ GL=YGL-YGT(YTT)-1YTL、Y+ GG= YGG-YGT(YTT)-1YTG；
(4) after the simplified admittance matrix of step (3) being removed into balance nodes associated row and related column, just subsequently calculated Simplify admittance matrix Y used in process+：
Preferably, the load bus is with active power or with reactive power injection rate；The connecting node Refer to the node of the injection rate of no any active power and reactive power.
Preferably, the step 2 confirms the idle output of the PV node by wide area measurement system, according to the letter Change admittance matrix Y+, calculate equivalent network admittance matrix YS：
Wherein, diag(QG) represent with power supply The reactive power sent；QGFor the diagonal matrix of diagonal entry；(Y+ LG)*Represent to Y+ LGEach element of matrix seeks conjugation.
Preferably, the equivalent network admittance matrix YSContain the power supply pact constant to voltage magnitude and active power output Beam, the voltage phase angle of power supply can change, the equivalent network admittance matrix YSThe maintenance active power output for considering generator is permanent Fixed and constant voltage magnitude constraint, and consider generator reactive and contribute change to the stable influence of system voltage.
Preferably, the part throttle characteristics in the step 3 includes：Invariable power model, permanent load model, constant-impedance model with And combinations thereof.
Preferably, for each load bus to voltage, the derivative of electric current forms load characteristic matrix L in the step 3, The expression formula of load characteristic matrix L is as follows：
The diagonal matrix that element therein is all made up of the voltage or electric current of load bus.
Preferably, the step 4 is by the load characteristic matrix L and the equivalent network admittance matrix YSIt is combined, The sensitivity of load voltage and load current to disturbance parameter is solved, that is, solves below equation：
Wherein
Id represents unit matrix, has added subscript d to be distinguished with electric current I；Then by voltage sensibility and current sensitive Degree is compared, and obtains Thevenin's equivalence impedance of the load bus based on sensitivity：
Voltage stabilization Judging index VS based on Thevenin's equivalence
Si=| Zth,i|/|Vi/Ii|
Work as SiCloser to 1, system operating point is closer to collapse of voltage point.
Preferably, the equivalent system matrix Y of the step 2sWith the load characteristic matrix L of step 3 in each state It is required for updating；If the topological structure of system changes, need to update simplified admittance matrix Y+ in step 1；If hair Motor reactive power reaches limitation, then needs to repartition load bus in Y+ and power supply node in step 1.
The invention discloses a kind of voltage stability assessment method for considering part throttle characteristics and power constraints, the present invention utilizes section Point voltage equation, angle calculation Thevenin's equivalence potential and equivalent reactance analytically.Consider part throttle characteristics and generator Static constraint, remain the information of complete power flow equation, it is not necessary to " system side Thevenin's equivalence parameter does not change " this One is irrational it is assumed that fundamentally avoiding above-mentioned technological deficiency.The system Thevenin's equivalence parameter arrived of this method is missed Difference is smaller, and the voltage stabilization of load also can be correctly assessed in complication system.
Brief description of the drawings
Fig. 1 is the flow chart for the voltage stability assessment method for considering part throttle characteristics and power constraints.
Fig. 2 is the node standard examples of I EEE 9.
Fig. 3 is the Voltage Stability Evaluation index curve map of load bus 5.
Fig. 4 is the calculating Thevenin's equivalence impedance plot of load bus 5.
Fig. 5 is the Voltage Stability Evaluation index curve map of load bus 7.
Fig. 6 is the calculating Thevenin's equivalence impedance plot of load bus 7.
Fig. 7 is the Voltage Stability Evaluation index curve map of load bus 9.
Fig. 8 is the calculating Thevenin's equivalence impedance plot of load bus 9.
Embodiment
To make the purpose, technical scheme and advantage that the present invention is implemented clearer, below in conjunction with the embodiment of the present invention Accompanying drawing, the technical scheme in the embodiment of the present invention is further described in more detail.In the accompanying drawings, identical from beginning to end or class As label represent same or similar element or the element with same or like function.Described embodiment is the present invention Part of the embodiment, rather than whole embodiments.The embodiments described below with reference to the accompanying drawings are exemplary, it is intended to uses It is of the invention in explaining, and be not considered as limiting the invention.Based on the embodiment in the present invention, ordinary skill people The every other embodiment that member is obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Embodiments of the invention are described in detail below in conjunction with the accompanying drawings.
Consider power constraints and the Voltage Stability Evaluation index calculating method of part throttle characteristics, comprise the following steps：
Step 1：The simplification admittance matrix Y of calculating network simplification Equivalent Model+
The topological structure of electric parameter measured according to wide area measurement system establishes the admittance matrix of the network.By in network Node is divided into power supply node, load bus and connecting node, and connecting node refers to no any active power and reactive power The node of injection rate.The nodal voltage equation of former network is as follows：
Wherein, VL,VT,VGThe load bus voltage for having power to inject, power supply node voltage (power supply node bag are represented respectively Include PV node and balance nodes, but in fact behind calculating do not need balance nodes, retain balance here in order to represent convenient Node) and connecting node voltage.In addition, when there is power supply node reactive power to reach limits value, this power supply node is from PV node Type is converted to PQ nodes, then this power supply node also translates into load bus processing, is no longer called power supply node.
Eliminate connecting node arrives simplification Equivalent Model.
Wherein, Y+ LL=YLL-YLT(YTT)-1YTL、Y+ LG=YLG-YLT(YLT)-1YTG、Y+ GL=YGL-YGT(YTT)-1YTL、Y+ GG= YGG-YGT(YTT)-1YTG。
After above-mentioned simplified admittance matrix is removed into balance nodes associated row and related column, follow-up calculating process institute is just obtained The simplification admittance matrix Y used+。
Step 2：Calculate the equivalent network admittance matrix Y for considering power constraintsS
The idle output of the PV node obtained according to wide area measurement system, it is possible to according to the admittance after above-mentioned network reduction Matrix, calculate an equivalent network admittance matrix.This equivalent network admittance matrix implicitly contains power supply to voltage amplitude Value and the constant constraint of active power output, the voltage phase angle of power supply can change.If it is assumed that the amplitude and phase angle of power supply are not sent out Changing, it can cause Voltage Stability Evaluation result that larger error occurs, and this computational methods are not this it is assumed that therefore Avoid this error.This is that other appraisal procedures are rare.
YSThe expression formula of matrix is as follows：
Wherein,
Diag (Q in formulaG) represent the reactive power Q that is sent with power supplyGFor the diagonal matrix of diagonal entry.(Y+ LG)* Represent to Y+ LGEach element of matrix seeks conjugation.Such representation is also taken below.
Step 3：Form the load characteristic matrix L related to part throttle characteristics
It is believed that load model has a significant impact to the voltage stability of system, it is contemplated that typical load is quiet Step response：Invariable power model, permanent load model, constant-impedance model and combinations thereof.Other loads can certainly be used Static characteristic expression-form.According to each load bus to voltage, the derivative of electric current forms load characteristic matrix, such as table 1.
If a load bus is the superposition of a variety of load types, then by corresponding element overlaid in table 1.
Element and load type corresponding table in the load characteristic matrix of table 1
Load type f(V,V*,I,I*) df/dV df/dV* df/dI df/dI* df/dλ
Constant-impedance f:Vi/Ii-λZ0i=0 1/Ii 0 -Vi/Ii 2 0 -Z0i
Constant current f:Ii-λI0i=0 0 0 1 0 -I0i
Invariable power f:ViIi *-λS0i=0 Ii * 0 0 Vi -S0i
The expression formula of load characteristic matrix L is as follows：
The diagonal matrix that element therein is all made up of the voltage or electric current of load bus.
In the step 3, it is contemplated that the static characteristic of load.In the analysis of voltage static stability, static load Analysis of the characteristic on voltage stabilization also influences very big.Therefore can easily be added in this computational methods various negative Lotus static characteristic, and a kind of advantage of this computational methods.
Step 4：The Thevenin's equivalence parameter of computing system load, calculate Voltage Stability Evaluation index VS.
By above load characteristic matrix L and system admittance matrix YSCombine, it is possible to solve load voltage, and load Sensitivity of the electric current to disturbance parameter, that is, solve following linear equation：
Wherein
IdUnit matrix is represented, has added subscript d to be distinguished with electric current I, then by voltage sensibility and current sensitive Degree is compared, and what is obtained is exactly the Thevenin's equivalence impedance based on sensitivity of load bus：
Voltage stabilization Judging index VS based on Thevenin's equivalence
Si=| Zth,i|/|Vi/Ii|
Work as SiCloser to 1, system operating point is closer to collapse of voltage point.
Based on power flow equation result, or the result based on state estimation, solve linear equation (1), it is possible to obtain Sensitivity of the electrical quantity to the load variations in system.Linear equation is solved, can disposably be tried to achieve under this state, it is all negative The voltage stability sensitivity index V of lotusSi, the voltage sensibility of wherein certain load can be onlyed demand.This depends on which load assigned Df/d λ be not equal to 0.
As shown in Fig. 2 it is proposed by the present invention consider part throttle characteristics and power constraints voltage stability assessment method be with IEEE9 node modular systems.Example includes three generators, wherein three power supply nodes (node 2,3 and a balance nodes 1), three load buses (node 5,7,9), three connecting nodes (node 4,6,8).Because constant power load model is easiest to that electricity occurs Destabilization problems are pressed, and constant power load model is more more complicated than constant-impedance and constant current in analysis and calculating, so in this example In, three load buses are to take constant power load model.The load power of three load buses increases in proportion, proportionality coefficient λ.This That method calculates is a state index S for assessing voltage stabilization, so the corresponding S value of each λ value.Following calculating process It is the detailed process for asking for state index S.Load characteristic matrix L and equivalent system matrix YsThen need calculating S indexs every time Shi Gengxin.
Step 1：The simplification admittance matrix Y of calculating network simplification Equivalent Model+
The topological structure of electric parameter measured according to wide area measurement system establishes the admittance matrix of the network.
Wherein, VL,VG,VTLoad bus voltage is represented respectively, (power supply node includes PV node peace to power supply node voltage Weigh node) and connecting node voltage.Wherein the order of load bus is 5,7, No. 9 nodes；The order of power supply node is 1,2, No. 3 Node；The order of connecting node is 4,6, No. 8 nodes.Eliminate connecting node arrives simplification Equivalent Model.
Wherein, Y+ LL=YLL-YLT(YTT)-1YTL、Y+ LG=YLG-YLT(YLT)-1YTG、Y+ GL=YGL-YGT(YTT)-1YTL、Y+ GG= YGG-YGT(YTT)-1YTG。
Simplified admittance matrix therein removes the Y after the associated row and related column of balance nodes 1+As subsequently calculated The admittance matrix that journey uses.System topological does not change, and this matrix does not just have to update again.Result of calculation is shown in Table 2.
Table 2 simplifies the result of calculation of admittance matrix Y+ (5 × 5) each element
Node number 5 7 8 2 3
5 2.046-12.110i -0.463+1.690i -0.678+3.092i 0 -0.453+3.000i
7 -0.463+1.690i 1.463-15.046i -0.550+2.298i -0.239+6.202i -0.218+5.209i
8 -0.678+3.092i -0.550+2.298i 1.713-12.924i -0.321+2.723i 0
2 0 -0.239+6.202i -0.321+2.723i 0.568-8.823i 0
3 -0.453+3.000i -0.218+5.209i 0 0 0.683-8.0598i
Step 2：Calculate the equivalent system admittance matrix Ys for considering power constraints
Wherein,
This admittance matrix changes with network operation state change, so needing to change and update with systematic parameter λ. In this example, 3 can be shown in Table in the hope of Ys as λ=1, Ys such as tables 4 are tried to achieve as λ=2.255.
, can be with Ys matrix each element values during the λ of table 3=1
1.878-11.567i -0.649+2.649i -0.678+3.092i 0.094-0.560i 0.265-0.940i 0
-0.649+2.6485i 1.152-11.159i -0.7017+3.266i 0.057-0.975i 1.013-3.750i 0.251-0.948i
-0.678+3.092i -0.702+3.267i 1.612-12.502i 0 0.396-0.897i 0.143-0.409i
0.094+0.561i 0.265+0.940i 0 1.878+11.567i -0.649-2.649i -0.678-3.092i
0.057+0.9753i 1.013+3.749i 0.25090+0.94769i -0.649-2.649i 1.152+11.160i -0.7017-3.267i
0 0.396+0.897i 0.143376+0.40856i -0.678-3.092i -0.7017-3.266i 1.612+12.502i
Table 4 is as λ=2.255, Ys matrix each element values
Step 3：Load characteristic matrix L is formed according to part throttle characteristics.
Three load buses are all constant power load model.Sodf/dIL=0,So Arrive：
In this example, as λ=1,5, table 6 can be shown in Table in the hope of L matrix each element values.Work as λ When=2.255, the value for trying to achieve L matrix each elements is shown in Table 7, table 8.
During the λ of table 5=1, L1 matrix each element values
0.899+0.371i 0 0 0 0 0
0 1.018+0.344i 0 0 0 0
0 0 1.262+0.620i 0 0 0
0 0 0 0.899-0.371i 0 0
0 0 0 0 1.018-0.344i
0 0 0 0 1.262-0.620i
During the λ of table 6=1, L2 matrix each element values
0 0 0 0.973-0.068i 0 0
0 0 0 0 0.986+0.011i 0
0 0 0 0 0 0.955-0.073i
0.973+0.068i 0 0 0 0 0
0 0.986-0.011i 0 0 0
0 0 0.955+0.073i 0 0 0
During the λ of table 7=2.255, L1 matrix each element values
1.397+2.668i 0 0 0 0 0
0 0.745+2.740i 0 0 0 0
0 0 1.592+4.437i 0 0 0
0 0 0 1.397-2.668i 0 0
0 0 0 0 0.745-2.740i 0
0 0 0 0 0 1.592-4.437i
During the λ of table 8=2.255, L2 matrix each element values
0 0 0 0.512-0.493i 0 0
0 0 0 0 0.477-0.693i 0
0 0 0 0 0 0.427-0.482i
0.512+0.493i 0 0 0 0 0
0 0.477+0.693i 0 0 0 0
0 0 0.427+0.482i 0 0 0
Step 4：The Thevenin's equivalence parameter of computing system load, and carry out Voltage Stability Evaluation.
By above load characteristic matrix L and system admittance matrix YSCombine, spirit of the solving system electrical quantity to parameter lambda Sensitivity, that is, solve below equation：
Then voltage sensibility is compared with current sensitivity, that obtain is exactly the Dai Wei based on sensitivity of load bus Southern equivalent impedance：
Voltage Stability Evaluation index VS based on Thevenin's equivalence
Si=| Zth,i|/|Vi/Ii|
Work as VSiCloser to 1, system operating point is closer to collapse of voltage point.
State under different load growth coefficients, corresponding different Voltage Stability Evaluation, is shown in Table 9.
Voltage Stability Evaluation VS values corresponding to the different load growth coefficient lambda of table 9
As can be seen here, for system in load growth coefficient lambda=2.255, there is collapse of voltage in system, at this moment, three loads Node voltage Stability Assessment index VS is equal to 1.Also the computational methods for illustrating this assessment voltage stabilization are correct and accurate 's.
As shown in Fig. 3,5 and 7, load bus 5, the voltage of load bus 7 and load bus 9 and VS indexs with λ curve Figure.
As shown in Fig. 4,6 and 8, the Static Equivalent impedance of load bus 5, load bus 7 and load bus 9, that is, VLi/ ILi, the impedance magnitude with system, that is, ZTh, with λ change curve.
The result of calculation that these curves embody VS indexs meets system time of day.
It is last it is to be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations.To the greatest extent The present invention is described in detail with reference to the foregoing embodiments for pipe, it will be understood by those within the art that：It is still Technical scheme described in foregoing embodiments can be modified, or which part technical characteristic is equally replaced Change；And these modifications or replacement, the essence of appropriate technical solution is departed from the essence of various embodiments of the present invention technical scheme God and scope.

Claims (9)

1. a kind of voltage stability assessment method for considering part throttle characteristics and power constraints, it is characterised in that comprise the following steps：
Step 1：The simplification admittance matrix Y of calculating network simplification Equivalent Model+；
Step 2：Calculate the equivalent network admittance matrix Y for considering power constraintsS；
Step 3：Form the load characteristic matrix L related to part throttle characteristics；
Step 4：The Thevenin's equivalence parameter of computing system load, and carry out Voltage Stability Evaluation.
2. the voltage stability assessment method according to claim 1 for considering part throttle characteristics and power constraints, it is characterised in that： Comprise the following steps in the step 1：
(1) the topological structure of electric parameter measured according to wide area measurement system establishes the admittance matrix of network, the node electricity of network Press equation as follows：
<mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>I</mi> <mi>L</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>G</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>Y</mi> <mrow> <mi>L</mi> <mi>L</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>L</mi> <mi>T</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>L</mi> <mi>G</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mrow> <mi>T</mi> <mi>L</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>T</mi> <mi>T</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>T</mi> <mi>G</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mrow> <mi>G</mi> <mi>L</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>G</mi> <mi>T</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Y</mi> <mrow> <mi>G</mi> <mi>G</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>V</mi> <mi>L</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>V</mi> <mi>T</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>V</mi> <mi>G</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
(2) node in network is divided into power supply node, load bus and connecting node, V in step (1)L,VT,VGRepresent respectively Have power inject load bus voltage, connecting node voltage and power supply node voltage, wherein power supply node include PV node and Balance nodes；
(3) connecting node is eliminated from the admittance matrix of step (1) and is simplified admittance matrix：
<mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>I</mi> <mi>L</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>I</mi> <mi>G</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>L</mi> <mi>L</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>L</mi> <mi>G</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>G</mi> <mi>L</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>G</mi> <mi>G</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>V</mi> <mi>L</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>V</mi> <mi>G</mi> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
Wherein, Y+ LL=YLL-YLT(YTT)-1YTL、Y+ LG=YLG-YLT(YLT)-1YTG、Y+ GL=YGL-YGT(YTT)-1YTL、Y+ GG=YGG- YGT(YTT)-1YTG；
(4) after the simplified admittance matrix of step (3) being removed into balance nodes associated row and related column, follow-up calculating process is just obtained Used simplified admittance matrix Y+：
<mrow> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>L</mi> <mi>L</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>L</mi> <mi>G</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>G</mi> <mi>L</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <msup> <mi>Y</mi> <mo>+</mo> </msup> <mrow> <mi>G</mi> <mi>G</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
3. the voltage stability assessment method according to claim 2 for considering part throttle characteristics and power constraints, it is characterised in that： The load bus is with active power or with reactive power injection rate；The connecting node refers to no any active The node of the injection rate of power and reactive power.
4. the voltage stability assessment method according to claim 3 for considering part throttle characteristics and power constraints, it is characterised in that： The step 2 confirms the idle output of the PV node by wide area measurement system, according to the simplified admittance matrix Y+, meter Calculate equivalent network admittance matrix YS：
<mrow> <mi>Y</mi> <mi>s</mi> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mover> <mi>Y</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>L</mi> <mi>L</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>M</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>M</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mover> <mi>Y</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>L</mi> <mi>L</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
Wherein, diag(QG) represent with power supply The reactive power sent；QGFor the diagonal matrix of diagonal entry；(Y+ LG)*Represent to Y+ LGEach element of matrix seeks conjugation.
5. the voltage stability assessment method according to claim 4 for considering part throttle characteristics and power constraints, it is characterised in that： The equivalent network admittance matrix YSContain the power supply constraint constant to voltage magnitude and active power output, the voltage phase angle of power supply It can change, the equivalent network admittance matrix YSThe maintenance active power output for considering generator is constant constant with voltage magnitude Constraint, and consider generator reactive and contribute change to the stable influence of system voltage.
6. the voltage stability assessment method according to claim 1 for considering part throttle characteristics and power constraints, it is characterised in that： Part throttle characteristics in the step 3 includes：Invariable power model, permanent load model, constant-impedance model and combinations thereof.
7. the voltage stability assessment method according to claim 6 for considering part throttle characteristics and power constraints, it is characterised in that： Each load bus forms load characteristic matrix L to the derivative of voltage, electric current in the step 3, load characteristic matrix L Expression formula is as follows：
<mrow> <mi>L</mi> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <mi>d</mi> <mi>V</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <msup> <mi>dV</mi> <mo>*</mo> </msup> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <mi>d</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <msup> <mi>dI</mi> <mo>*</mo> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <msup> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <msup> <mi>dV</mi> <mo>*</mo> </msup> <mo>)</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <msup> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <mi>d</mi> <mi>V</mi> <mo>)</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <msup> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <msup> <mi>dI</mi> <mo>*</mo> </msup> <mo>)</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <msup> <mrow> <mo>(</mo> <mi>d</mi> <mi>f</mi> <mo>/</mo> <mi>d</mi> <mi>I</mi> <mo>)</mo> </mrow> <mo>*</mo> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
The diagonal matrix that element therein is all made up of the voltage or electric current of load bus.
8. the voltage stability assessment method according to claim 7 for considering part throttle characteristics and power constraints, it is characterised in that： The step 4 is by the load characteristic matrix L and the equivalent network admittance matrix YSBe combined, solve load voltage and Sensitivity of the load current to disturbance parameter, that is, solve below equation：
<mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>L</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>L</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>s</mi> </msub> </mtd> <mtd> <msub> <mi>I</mi> <mi>d</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>DV</mi> <mi>L</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>DI</mi> <mi>L</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <mi>D</mi> <mi>f</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>
Wherein
IdUnit matrix is represented, has added subscript d to be distinguished with electric current I；Then by voltage sensibility and current sensitivity phase Than obtaining Thevenin's equivalence impedance of the load bus based on sensitivity：
<mrow> <msub> <mi>Z</mi> <mrow> <mi>T</mi> <mi>h</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>dV</mi> <mi>L</mi> </msub> <mo>/</mo> <msubsup> <mi>dI</mi> <mi>L</mi> <mo>*</mo> </msubsup> </mrow>
Voltage stabilization Judging index VS based on Thevenin's equivalence
Si=| Zth,i|/|Vi/Ii|
Work as SiCloser to 1, system operating point is closer to collapse of voltage point.
9. the voltage stability assessment method according to claim 8 for considering part throttle characteristics and power constraints, it is characterised in that： The equivalent system matrix Y of the step 2sIt is required for updating in each state with the load characteristic matrix L of step 3；If The topological structure of system changes, then needs to update simplified admittance matrix Y in step 1+；If generator reactive power reaches To limitation, then need Y in step 1+Middle load bus and power supply node are repartitioned.
CN201711139795.4A 2017-11-16 2017-11-16 Voltage stability evaluation method considering load characteristics and power supply constraints Active CN107832959B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711139795.4A CN107832959B (en) 2017-11-16 2017-11-16 Voltage stability evaluation method considering load characteristics and power supply constraints

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711139795.4A CN107832959B (en) 2017-11-16 2017-11-16 Voltage stability evaluation method considering load characteristics and power supply constraints

Publications (2)

Publication Number Publication Date
CN107832959A true CN107832959A (en) 2018-03-23
CN107832959B CN107832959B (en) 2021-12-03

Family

ID=61651858

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711139795.4A Active CN107832959B (en) 2017-11-16 2017-11-16 Voltage stability evaluation method considering load characteristics and power supply constraints

Country Status (1)

Country Link
CN (1) CN107832959B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110414816A (en) * 2019-07-19 2019-11-05 福建睿思特科技股份有限公司 One kind being based on least square power system state estimation method
CN112421619A (en) * 2020-11-09 2021-02-26 广东电网有限责任公司江门供电局 Intelligent low-voltage distribution network voltage stability real-time monitoring method
GB2616276A (en) * 2022-03-01 2023-09-06 Univ Warwick Network identification

Citations (14)

* Cited by examiner, † Cited by third party
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
CN102420427A (en) * 2011-12-02 2012-04-18 余杭供电局 Method for calculating voltage stability margin of regional power grid by considering outer network equivalence
CN103258103A (en) * 2013-05-29 2013-08-21 国家电网公司 Thevenin equivalent parameter identification method based on partial least squares regression
CN103606921A (en) * 2013-12-10 2014-02-26 上海交通大学 Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence
CN103760400A (en) * 2014-01-08 2014-04-30 东北电力大学 Voltage stabilization on-line monitoring method based on wide area measuring system
CN104124684A (en) * 2014-07-25 2014-10-29 国家电网公司 Wide-area measurement information based power system voltage stability online evaluation method
CN104734148A (en) * 2015-03-16 2015-06-24 河海大学 Three-phrase power-distributing network continuation power flow analysis of distributed power supply
CN104778367A (en) * 2015-04-20 2015-07-15 山东大学 On-line wide-area Thevenin equivalent parameter calculation method based on single state section
CN104917174A (en) * 2015-06-08 2015-09-16 北京交通大学 Static voltage stability judgment method for active power distribution network
CN105356447A (en) * 2015-10-27 2016-02-24 河海大学 Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system
CN106099914A (en) * 2016-06-29 2016-11-09 燕山大学 A kind of active distribution network safety status classification method based on power supply capacity
CN106295160A (en) * 2016-08-04 2017-01-04 山东大学 AC-DC interconnecting power network Thevenin's equivalence parameter on-line calculation method
CN106682407A (en) * 2016-12-19 2017-05-17 国网上海市电力公司 Voltage stability assessment method based on thevenin equivalence and branch transmission power limits
CN106684857A (en) * 2016-09-05 2017-05-17 国网江苏省电力公司电力科学研究院 Linear optimal power flow model containing unified power flow controller

Patent Citations (14)

* Cited by examiner, † Cited by third party
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
CN102420427A (en) * 2011-12-02 2012-04-18 余杭供电局 Method for calculating voltage stability margin of regional power grid by considering outer network equivalence
CN103258103A (en) * 2013-05-29 2013-08-21 国家电网公司 Thevenin equivalent parameter identification method based on partial least squares regression
CN103606921A (en) * 2013-12-10 2014-02-26 上海交通大学 Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence
CN103760400A (en) * 2014-01-08 2014-04-30 东北电力大学 Voltage stabilization on-line monitoring method based on wide area measuring system
CN104124684A (en) * 2014-07-25 2014-10-29 国家电网公司 Wide-area measurement information based power system voltage stability online evaluation method
CN104734148A (en) * 2015-03-16 2015-06-24 河海大学 Three-phrase power-distributing network continuation power flow analysis of distributed power supply
CN104778367A (en) * 2015-04-20 2015-07-15 山东大学 On-line wide-area Thevenin equivalent parameter calculation method based on single state section
CN104917174A (en) * 2015-06-08 2015-09-16 北京交通大学 Static voltage stability judgment method for active power distribution network
CN105356447A (en) * 2015-10-27 2016-02-24 河海大学 Analysis method for steady-state security region of electric-pneumatic interconnected integrated energy system
CN106099914A (en) * 2016-06-29 2016-11-09 燕山大学 A kind of active distribution network safety status classification method based on power supply capacity
CN106295160A (en) * 2016-08-04 2017-01-04 山东大学 AC-DC interconnecting power network Thevenin's equivalence parameter on-line calculation method
CN106684857A (en) * 2016-09-05 2017-05-17 国网江苏省电力公司电力科学研究院 Linear optimal power flow model containing unified power flow controller
CN106682407A (en) * 2016-12-19 2017-05-17 国网上海市电力公司 Voltage stability assessment method based on thevenin equivalence and branch transmission power limits

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
付正鑫: "基于广域测量技术戴维南等值参数的计算及在电压稳定性分析中的应用", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 *
周旭: "计及电网拓扑优化的电压稳定性研究", 《中国优秀硕士学位论文全文库》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110414816A (en) * 2019-07-19 2019-11-05 福建睿思特科技股份有限公司 One kind being based on least square power system state estimation method
CN110414816B (en) * 2019-07-19 2021-11-23 福建睿思特科技股份有限公司 Least square-based power system state estimation method
CN112421619A (en) * 2020-11-09 2021-02-26 广东电网有限责任公司江门供电局 Intelligent low-voltage distribution network voltage stability real-time monitoring method
GB2616276A (en) * 2022-03-01 2023-09-06 Univ Warwick Network identification

Also Published As

Publication number Publication date
CN107832959B (en) 2021-12-03

Similar Documents

Publication Publication Date Title
Du et al. On-line estimation of transmission line parameters, temperature and sag using PMU measurements
EP1841037B1 (en) System, method, and article of manufacture for determining parameter values associated with an electrical grid
CN104778367B (en) Wide area Thevenin&#39;s equivalence parameter on-line calculation method based on a single state section
EP1841036A2 (en) System, method and article of manufacture for controlling operation of an electrical power generation system
CN107832959A (en) A kind of voltage stability assessment method for considering part throttle characteristics and power constraints
CN106682407B (en) Voltage stability assessment method based on Thevenin&#39;s equivalence and branch transmission power limit
CN103760400B (en) Voltage stabilization on-line monitoring method based on wide area measurement system
CN103595040A (en) Power distribution network comprising direct-current distributed power supply and analysis method of comprehensive load characteristics of power distribution network
CN106443246A (en) Online identification method of small-disturbance stability parameters based on PMU measurement data
EP3041107A1 (en) Method and system for identifying element parameter and power correction factor of electric power system
CN104009475B (en) Non-topological static equivalence method taking equivalent element comprehensiveness and parameter physical restriction into consideration
CN103023418A (en) Online parameter identification method of synchronous generator based on wide-area measurement information
CN107609802A (en) A kind of alternating current-direct current series-parallel connection Power Systems/voltage stability margin appraisal procedure containing multi-infeed HVDC
CN104124684B (en) Based on the Power System Voltage Stability online evaluation method of wide area measurement information
CN104182644A (en) Power distribution network state estimation method fusing distributed generation characteristics
CN104142419B (en) Power grid short-circuit current acquiring method with influences of loads considered
Prasad et al. Assessment on various node voltage stability indices—A review
CN104240036B (en) A kind of Transient Voltage Stability quantitative estimation method based on critical system equiva lent impedance
CN108196150A (en) A kind of asymmetric transmission line parameter measuring method of common-tower double-return
CN111211581B (en) Method and device for determining new energy acceptance capacity in power system
Danish et al. Comparative analysis of load flow calculation methods with considering the voltage stability constraints
CN104934969A (en) Method for calculating parameter of electric power line
CN109342817A (en) A kind of non-fully commutation three-phase power linear electrical parameter estimation method based on PMU measurement
Raju et al. Voltage instability identification with modified L-index using synchrophasor data
CN113536674B (en) Line parameter identification method based on BP neural network and improved SCADA data

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