CN104868478A - Method for starting dynamic partitioning scheme under condition of power grid emergency - Google Patents

Method for starting dynamic partitioning scheme under condition of power grid emergency Download PDF

Info

Publication number
CN104868478A
CN104868478A CN201510305457.8A CN201510305457A CN104868478A CN 104868478 A CN104868478 A CN 104868478A CN 201510305457 A CN201510305457 A CN 201510305457A CN 104868478 A CN104868478 A CN 104868478A
Authority
CN
China
Prior art keywords
power grid
msub
mover
impedance
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
Application number
CN201510305457.8A
Other languages
Chinese (zh)
Other versions
CN104868478B (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201510305457.8A priority Critical patent/CN104868478B/en
Publication of CN104868478A publication Critical patent/CN104868478A/en
Application granted granted Critical
Publication of CN104868478B publication Critical patent/CN104868478B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method for starting a dynamic partitioning scheme under the condition of power grid emergency, and the method comprises the following steps: 10) measuring and calculating an impedance modulus index value of a weak node in each partitioned power grid; 20) building a relation between a power margin and an impedance modulus index; 30) measuring and calculating a threshold value of the impedance modulus index according to a set power margin limit and a formula (1); 40) comparing the impedance modulus index values, obtained at step 10), of all weak node in the partitioned power grids with the threshold value, measured at step 30), of the impedance modulus index: starting the dynamic partitioning scheme if the impedance modulus index value of any one weak node in the partitioned power grid is less than the threshold value of the impedance modulus index; or else, indicating that the voltage of the partitioned power grids are safe at this moment, and returning to step 10). The method can achieve the starting of dynamic partitioning under the condition of power grid emergency, and guarantees the reliability of power supply.

Description

Method for starting dynamic partitioning scheme in power grid emergency state
Technical Field
The invention belongs to the technical field of emergency control of a large power grid, and particularly relates to a method for starting a dynamic partitioning scheme in an emergency state of the power grid.
Background
The electrical connection among large power grids is greatly enhanced by the interconnection of the regional power grids, the rapid development of the extra-high voltage alternating current and direct current system and the basic formation of a long-distance large-capacity power transmission pattern. The development of the global and integrated power pattern leads to more complex power grid stability characteristics and aggravates dynamic interaction influence among regional power grids.
Under the background, whether the load can be quickly controlled under the emergency condition of the power grid and even under the fault condition is the key for determining that the power grid keeps safe and stable operation by reasonably recombining the subarea power supply network by using the original line, further controlling the power flow of related elements and improving the subarea reactive voltage support function. For the power grid which realizes the layered partition operation, the dynamic partition technology is undoubtedly an effective means for preventing and controlling the power grid operation and enabling the power grid to exit the emergency state. In addition, the national institute of Electrical safety Accident No. 599 Emergency handling and survey handling regulations (hereinafter, 599 Command) make clear regulations on the load shedding proportion of the power grid in different areas under the condition of different degrees of accidents of the power system, and the reasonable application of the power grid dynamic partitioning technology can actively respond to 599 Command, reduce the load shedding proportion to the maximum extent and ensure the reliability of power supply.
The dynamic partitioning of the power grid is a comprehensive problem, and needs to have an overall grasp on the overall situation of the power grid (including aspects of safety and stability, short-circuit current, power supply reliability and the like). The core of the dynamic partitioning is to judge whether a partitioning scheme needs to be adopted or not according to the serious condition of the system, and to formulate a reasonable starting criterion for starting the dynamic partitioning scheme. At present, related researches are relatively lacked in the aspect of partition scheme starting criteria.
Disclosure of Invention
The technical problem is as follows: the technical problem to be solved by the invention is as follows: the method for starting the dynamic partitioning scheme in the power grid emergency state is provided, so that the dynamic partitioning is started in the power grid emergency state, and the power supply reliability is ensured.
The technical scheme is as follows: in order to solve the technical problems, the invention adopts the technical scheme that:
a method for starting a dynamic partitioning scheme in an emergency state of a power grid comprises the following steps:
step 10), measuring and calculating impedance mode index values of weak nodes in each subarea power grid in the power grid;
step 20) establishing a relational expression between the power margin and the impedance mode index as shown in the formula (1);
formula (1)
In the above formula, KPRepresenting a power margin; vindexRepresenting an impedance mode index; theta represents the equivalent impedance angle of the network;representing the load impedance angle;
step 30), according to the set power margin limit value and the formula (1), measuring and calculating an impedance mode index threshold value;
step 40) comparing the impedance mode index values of all weak nodes in the partitioned power grid measured in the step 10) with the impedance mode index threshold measured in the step 30), and if the impedance mode index value of any weak node in the partitioned power grid is smaller than the impedance mode index threshold, starting a dynamic partitioning scheme; otherwise, the voltage of the subarea power grid at the moment is safe, and the step 10) is returned.
Further, the impedance mode index threshold is 0.358.
Has the advantages that: compared with the prior art, the invention has the following advantages:
the method comprises the steps of firstly determining that node voltage is subjected to stability analysis by using impedance mode indexes based on Thevenin equivalence, then determining the relation between the impedance mode indexes and power margins, determining an impedance mode index threshold value in an emergency state of a power grid according to a limit value of the power margins determined in advance, finally performing index calculation on weak nodes in partitions, and starting a dynamic partition scheme once the impedance mode index value of the node is lower than the impedance mode index threshold value. The method skillfully associates the voltage index value with the given power margin and calculates the impedance mode index threshold value, so that the starting threshold value of the partition scheme is more real and credible. In addition, the method of the invention considers the situation that the power supply side has insufficient active/reactive response supporting capability and the load sudden change occurs under the heavy load condition. These factors will mainly cause frequency, voltage, short circuit current problems in the partitioned grid. For the existing interconnected power grid pattern of the regional power grid, the probability of the occurrence of general frequency problems is not very high unless island operation occurs, and the general frequency problems are mainly the voltage problems. Therefore, the voltage index is considered in the starting criterion of the dynamic partitioning scheme in the emergency state of the power grid. This enables more efficient and accurate start-up. In the method, the static voltage indexes are adopted for analysis, the voltage stability level of each node in the subarea can be accurately reflected, whether the subarea power grid is in an emergency state or not is judged, and the method has the advantages of simplicity in calculation, strong operability and the like.
Drawings
FIG. 1 is a block flow diagram of the present invention;
FIG. 2 is a system diagram of node voltage Thevenin equivalence.
Detailed Description
The invention is described in further detail below with reference to embodiments and with reference to the drawings. The invention is not limited to the examples given.
The emergency state of the power grid mainly refers to the state that the power grid enters an unstable state after various operation constraint conditions are broken, and the dynamic partitioning technology mainly aims at the state that the power grid enters an abnormal operation state due to the vulnerability of the power grid after the power grid is disturbed from the inside and the outside. Such a state may not necessarily cause the constraint condition to be out of limit. When the margin left by the operation constraint is smaller than a certain threshold value, the power grid can be considered to enter an emergency state. Factors which can cause the power grid to enter an emergency state when the power supply side and the load side change excessively are considered, such as insufficient active/reactive response supporting capacity of the power supply, severe fluctuation of intermittent new energy power of wind power and the like, sudden load change under heavy load and the like.
As shown in fig. 1, a method for starting a dynamic partitioning scheme in an emergency state of a power grid according to an embodiment of the present invention includes the following steps:
and step 10) measuring and calculating impedance mode index values of weak nodes in each subarea power grid in the power grid.
Step 20) establishing a relational expression between the power margin and the impedance mode index as shown in the formula (1);
formula (1)
In the above formula, KPRepresenting a power margin; vindexRepresenting an impedance mode index; theta represents the equivalent impedance angle of the network;representing the load impedance angle.
Step 30), according to the set power margin limit value and the formula (1), measuring and calculating an impedance mode index threshold value;
step 40) comparing the impedance mode index values of all weak nodes in the partitioned power grid measured in the step 10) with the impedance mode index threshold measured in the step 30), and if the impedance mode index value of any weak node in the partitioned power grid is smaller than the impedance mode index threshold, starting a dynamic partitioning scheme; otherwise, the voltage of the subarea power grid at the moment is safe, and the step 10) is returned.
And step 10), measuring and calculating an impedance mode index value of a weak node, and firstly, carrying out Thevenin equivalent parameter calculation on the weak node in the partitioned power grid to obtain the equivalent impedance of the Thevenin system, thereby calculating the impedance mode index value of the weak node.
A real-time thevenin equivalent circuit diagram as shown in fig. 2 is established. In FIG. 2, all the electromotive forces in the system are equalized to one equivalent electromotive forceRepresents; similarly, the line impedance except the node is used as an equivalent impedanceTo indicate that the user is not in a normal position,representing the node load impedance. Aiming at any weak node of an electric field of a certain subarea, a complex system is converted into a simple alternating current system through thevenin equivalence. The method comprises the following specific steps:
step 101) when the real-time Thevenin equivalent parameter of the system is calculated, firstly, the initial value of the parameter is calculated by using the data obtained by the previous two times of sampling and processing. According to the formula (201), obtaining the initial value of thevenin equivalent impedance based on the traditional two-point method
Formula (201)
Obtaining the initial value of the system Thevenin equivalent potential according to the formula (202)
<math> <mrow> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mi>eq</mi> <mn>0</mn> </mrow> </msub> <mo>=</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mi>eq</mi> <mn>0</mn> </mrow> </msub> <msub> <mi>I</mi> <mrow> <mi>L</mi> <mn>0</mn> </mrow> </msub> <mo>+</mo> <msub> <mover> <mi>U</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mi>L</mi> <mn>0</mn> </mrow> </msub> </mrow> </math> Formula (202)
In the formula,representing the first sampled voltage phasor, UL0Representing the magnitude of the voltage sampled for the first time,representing the voltage phase angle obtained by the first sampling;representing the second sampled voltage phasor, UL1Representing the magnitude of the voltage obtained by the second sampling,representing the voltage phase angle obtained by the second sampling; i isL0Representing the current amplitude obtained by the first sampling; i isL1Representing the current magnitude obtained by the second sampling.
Calculating the corresponding equivalent impedance prediction value according to the formula (203)
Formula (203)
Wherein,showing that the ith iteration obtains the predicted value of the equivalent impedance,representing the predicted value of the equivalent resistance obtained by the ith iteration,representing the predicted value of equivalent reactance obtained in the ith iteration, Eeqi-1Representing the calculated equivalent potential amplitude for the (i-1) th iteration,eqi-1representing the i-1 th iteratively calculated equivalent potential phase angle, ULiRepresenting the magnitude of the voltage obtained from the ith sample,representing the phase angle of the voltage obtained by the ith sampling, ILiRepresenting the current magnitude obtained from the ith sample.
For the ith sampling, the voltage phasor obtained by the ith sampling is utilizedCurrent data ILiAnd the amplitude E of the equivalent potential obtained by the i-1 th calculationeqi-1Phase angle of sum potentialeqi-1Calculating to obtain the equivalent potential correction quantityThe respective intermediate judgment variable Δ Z and the angle correction amount Δ are calculated according to equation (204):
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>&Delta;Z</mi> <mo>=</mo> <mo>|</mo> <msub> <mover> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mo>~</mo> </mover> <mi>eqi</mi> </msub> <mo>|</mo> <mo>-</mo> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mi>eqi</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>|</mo> </mtd> </mtr> <mtr> <mtd> <mi>&Delta;&delta;</mi> <mo>=</mo> <mo>|</mo> <mi>arctan</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mover> <mi>X</mi> <mo>~</mo> </mover> <mi>eqi</mi> </msub> <mo>-</mo> <msub> <mi>X</mi> <mrow> <mi>eqi</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mover> <mi>R</mi> <mo>~</mo> </mover> <mi>eqi</mi> </msub> <mo>-</mo> <msub> <mi>R</mi> <mrow> <mi>eqi</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>|</mo> </mtd> </mtr> </mtable> </mfenced> </math> formula (204)
Thereby obtaining the correction quantity of the equivalent potential
If it is U Li I Li > U Li - 1 I Li - 1 , And Δ Z>0, then <math> <mrow> <mi>&Delta;</mi> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>=</mo> <mo>-</mo> <mn>0.01</mn> <mo>&angle;</mo> <mi>&Delta;&delta;</mi> <mo>;</mo> </mrow> </math>
If it is U Li I Li > U Li - 1 I Li - 1 , And Δ Z<0, then <math> <mrow> <mi>&Delta;</mi> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>=</mo> <mn>0.01</mn> <mo>&angle;</mo> <mi>&Delta;&delta;</mi> <mo>;</mo> </mrow> </math>
If it is U Li I Li < U Li - 1 I Li - 1 , And Δ Z>0, then <math> <mrow> <mi>&Delta;</mi> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>=</mo> <mn>0.01</mn> <mo>&angle;</mo> <mi>&Delta;&delta;</mi> <mo>;</mo> </mrow> </math>
If it is U Li I Li < U Li - 1 I Li - 1 , And Δ Z<0, then <math> <mrow> <mi>&Delta;</mi> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>=</mo> <mo>-</mo> <mn>0.01</mn> <mo>&angle;</mo> <mi>&Delta;&delta;</mi> <mo>;</mo> </mrow> </math>
If it is U Li I Li = U Li - 1 I Li - 1 , Then <math> <mrow> <mi>&Delta;</mi> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>=</mo> <mn>0</mn> <mo>.</mo> </mrow> </math>
In the formula,represents the i-1 th iteration to obtain the equivalent impedance Reqi-1Representing the equivalent resistance value, X, obtained in the i-1 st iterationeqi-1Represents the equivalent reactance value, U, obtained in the i-1 th iterationLi-1Represents the voltage amplitude, I, obtained by the I-1 th samplingLi-1The current amplitude obtained by the i-1 th sampling is shown.
Therefore, the Thevenin equivalent potential phasor obtained by the ith sampling is calculated according to the formula (205)
<math> <mrow> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eqi</mi> </msub> <mo>=</mo> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mrow> <mi>eqi</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mi>&Delta;</mi> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> </mrow> </math> Formula (205)
Wherein,representing thevenin equivalent potential phasor obtained by sampling at the (i) -1 st time,an equivalent potential correction amount is indicated,and expressing thevenin equivalent potential phasor obtained by the ith sampling.
Obtaining thevenin equivalent potential phasor corresponding to the ith sampling in calculationThen, substituting the equation into a formula (211) to calculate the equivalent impedance phasor after the ith sampling
<math> <mrow> <msub> <mover> <mi>E</mi> <mo>&CenterDot;</mo> </mover> <mi>eqi</mi> </msub> <mo>=</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>eqi</mi> </msub> <msub> <mi>I</mi> <mi>Li</mi> </msub> <mo>+</mo> <msub> <mover> <mi>U</mi> <mo>&CenterDot;</mo> </mover> <mi>Li</mi> </msub> </mrow> </math> Formula (211)
Wherein,representing the thevenin equivalent potential phasor obtained by the ith sampling,represents the equivalent impedance phasor, I, after the ith samplingLiRepresenting the current amplitude, U, obtained from the ith sampleLiRepresenting the voltage amplitude obtained by the ith sampling.
And thus, the calculation of the real-time Thevenin equivalent parameters of the ith sampling data is completed, the ith +1 sampling data is continuously obtained, and the calculation of the real-time Thevenin equivalent parameters is repeated.
Step 102) measuring and calculating equivalent impedance of the Thevenin system:
as can be seen from the Thevenin equivalent circuit diagram shown in FIG. 2, the apparent power obtained by the load node in the polar coordinate system is as shown in equation (206):
<math> <mrow> <msub> <mover> <mi>S</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <mo>=</mo> <msub> <mover> <mi>U</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <msubsup> <mover> <mi>I</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> <mo>*</mo> </msubsup> <mo>=</mo> <msub> <mover> <mi>I</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <msubsup> <mover> <mi>I</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> <mo>*</mo> </msubsup> <mo>=</mo> <mfrac> <msubsup> <mi>E</mi> <mi>eq</mi> <mn>2</mn> </msubsup> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>+</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> </mrow> </math> formula (206)
As shown in equation (207), the load node absorbs the active power PLReactive power QLAre respectively as
Formula (207)
In the formula, theta represents the equivalent impedance angle of the network,representing the load impedance angle, EeqRepresenting the amplitude of thevenin equivalent electromotive force; zeqRepresenting Thevenin equivalent impedance modulus, ZLRepresenting the magnitude of the node load impedance modulus,representing a node load impedance;representing thevenin equivalent impedance;
PLand QLTo ZLDerivative to obtain ZL=ZeqWhen the temperature of the water is higher than the set temperature,at the moment, the active power transmitted by the line reaches the maximum value PLmaxAnd reactive power Q obtained at the load sideLmaxAlso, it reaches a maximum, as shown in equation (208):
formula (208)
Step 103) measuring and calculating an impedance model index test value:
from step 102), when the transmission power of the line is maximum, the voltage is in a critical stable state, and the load impedance and the thevenin equivalent impedance value are equal, that is, the equation (209):
<math> <mrow> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <mo>|</mo> <mo>=</mo> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>|</mo> </mrow> </math> formula (209)
And judging the voltage stability by using an impedance mode index, wherein the index test value is defined as:
<math> <mrow> <msub> <mi>V</mi> <mi>index</mi> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <mo>|</mo> </mrow> </mfrac> </mrow> </math> formula (210)
The impedance model index test value accurately reflects the difference between the current operation state point and the maximum transmission power point of the system. When the index test value is closer to 0, the operation point is closer to the voltage stability critical point, and the stability of the node voltage is poorer.
In step 20), the process of establishing the relation between the power margin and the impedance mode index is as follows:
power margin KPIs defined as
K P = P L max P L - 1 Formula (301)
Combining formula (207) with formula (301) to obtain KPAs shown in formula (302):
formula (302)
From formula (210)
<math> <mrow> <mfrac> <mrow> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>eq</mi> </msub> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <msub> <mover> <mi>Z</mi> <mo>&CenterDot;</mo> </mover> <mi>L</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>=</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>V</mi> <mi>index</mi> </msub> </mrow> </math> Formula (303)
Then, combining the formula (302) and the formula (303) to obtain the power margin KPAnd an impedance mode index VindexThe relationship between
In step 30), the power margin limit value is set in advance, and may be defined according to the operation rule in the field. For example, according to the technical Specification for calculating safety and stability of Power System (Standard number DL/T1234-2013) on Power margin KPAnd deducing a corresponding impedance mode index threshold value as a criterion whether the node voltage reaches an emergency state. The technical specification of safety and stability calculation of the electric power system indicates that K is used when the regional load is maximum or the section current reaches the maximumPNot less than 8%; under N-1 failure, KPNot less than 5%. This gives the formula (401):
formula (401)
This gives:
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>V</mi> <mi>index</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0.413</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>V</mi> <mi>index</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0.358</mn> </mtd> </mtr> </mtable> </mfenced> </math>
thus, V can be adjustedindex<0.358 as a dynamic partition launch criterion.
An example is illustrated below.
Taking a weak node of 230kV Shanghai Kyoho 1# as an example, carrying out load flow calculation and stable file operation through BPA software, wherein the simulation time is 10s, obtaining real-time information (including active power, reactive power, voltage amplitude and phase angle) of the weak node, importing the data into MATLAB, and carrying out programming and voltage index calculation to obtain results shown in Table 1.
TABLE 1230 kV Shanghai Kyoho No. 1 Voltage index results
Time/s Voltage index value Time/s Voltage index value
0.5 0.763099818763312 5.5 0.662943865068431
1.0 0.679590478635123 6.0 0.667358883729164
1.5 0.498000794003520 6.5 0.675677440313633
2.0 0.457948311073828 7.0 0.685469870716686
2.5 0.462524250660278 7.5 0.695403741868622
3.0 0.465589621643638 8.0 0.703210457591864
3.5 0.453258840128252 8.5 0.707505750033555
4.0 0.404039885804003 9.0 0.708775548352132
4.5 0.617469536821706 9.5 0.709931902035739
5.0 0.572478419419662 10.0 0.711083592870503
As can be seen from the data in Table 1, the voltage index value V of the node is within 10s of the simulationindex>0.358, which shows the stability of the bus voltage of No. 1 Hu Kyoff at 230 kV. If the partition scheme needs to be started, the same calculation needs to be performed on other weak nodes of the whole partition, and the conclusion can be reached only if the observation index value exceeds 0.358. The index calculation of other nodes is the same as the calculation method of the weak nodes, and is not repeated.
It should be noted that: the described embodiments are only some embodiments of the present application and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.

Claims (2)

1. A method for starting a dynamic partitioning scheme in an emergency state of a power grid is characterized by comprising the following steps:
step 10), measuring and calculating impedance mode index values of weak nodes in each subarea power grid in the power grid;
step 20) establishing a relational expression between the power margin and the impedance mode index as shown in the formula (1);
formula (1)
In the above formula, KPRepresenting a power margin; vindexRepresenting an impedance mode index; theta represents the equivalent impedance angle of the network;representing the load impedance angle;
step 30), according to the set power margin limit value and the formula (1), measuring and calculating an impedance mode index threshold value;
step 40) comparing the impedance mode index values of all weak nodes in the partitioned power grid measured in the step 10) with the impedance mode index threshold measured in the step 30), and if the impedance mode index value of any weak node in the partitioned power grid is smaller than the impedance mode index threshold, starting a dynamic partitioning scheme; otherwise, the voltage of the subarea power grid at the moment is safe, and the step 10) is returned.
2. A method for starting a dynamic partitioning scheme for grid emergency situations as claimed in claim 1, wherein said impedance mode indicator threshold is 0.358.
CN201510305457.8A 2015-06-05 2015-06-05 A kind of method of the startup dynamic partition scheme under the power network state of emergency Expired - Fee Related CN104868478B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510305457.8A CN104868478B (en) 2015-06-05 2015-06-05 A kind of method of the startup dynamic partition scheme under the power network state of emergency

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510305457.8A CN104868478B (en) 2015-06-05 2015-06-05 A kind of method of the startup dynamic partition scheme under the power network state of emergency

Publications (2)

Publication Number Publication Date
CN104868478A true CN104868478A (en) 2015-08-26
CN104868478B CN104868478B (en) 2017-07-07

Family

ID=53914111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510305457.8A Expired - Fee Related CN104868478B (en) 2015-06-05 2015-06-05 A kind of method of the startup dynamic partition scheme under the power network state of emergency

Country Status (1)

Country Link
CN (1) CN104868478B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105811404A (en) * 2016-03-22 2016-07-27 山东大学 Stable situation monitoring method for quiescent voltage of distribution network with synergic transmission and distribution
CN106684859A (en) * 2016-11-24 2017-05-17 东南大学 Power grid short circuit current suppression method based on dynamic partition technique
CN108493979A (en) * 2018-03-28 2018-09-04 国电南瑞科技股份有限公司 A kind of extra-high voltage direct-current access provincial power network partition zone optimizing progress control method
CN108631290A (en) * 2017-03-21 2018-10-09 国家电网公司 A kind of multi-infeed DC receiving end power grid direct current drop point screening technique
CN116937576A (en) * 2023-09-15 2023-10-24 国网山西省电力公司运城供电公司 Large-scale active power distribution network stability evaluation method, system and equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1324458A2 (en) * 2001-12-24 2003-07-02 ABB Technology AG A method of the identification of weak and/or strong branches of an electric power transmission system
CN101118265A (en) * 2007-09-17 2008-02-06 重庆大学 Process for real time recognizing voltage stability of electrified wire netting trough recognizing weak links of electric network
CN102055191A (en) * 2010-12-31 2011-05-11 重庆大学 Wide area load shedding control method for quiescent voltage stabilization considering external power grid equivalence
CN103944170A (en) * 2014-04-16 2014-07-23 西安交通大学 Self-adaption urgent load shedding method based on electric system response information

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1324458A2 (en) * 2001-12-24 2003-07-02 ABB Technology AG A method of the identification of weak and/or strong branches of an electric power transmission system
CN101118265A (en) * 2007-09-17 2008-02-06 重庆大学 Process for real time recognizing voltage stability of electrified wire netting trough recognizing weak links of electric network
CN102055191A (en) * 2010-12-31 2011-05-11 重庆大学 Wide area load shedding control method for quiescent voltage stabilization considering external power grid equivalence
CN103944170A (en) * 2014-04-16 2014-07-23 西安交通大学 Self-adaption urgent load shedding method based on electric system response information

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯光明: "电力系统电压控制区域的研究", 《CNKI优秀硕士论文全文数据库》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105811404A (en) * 2016-03-22 2016-07-27 山东大学 Stable situation monitoring method for quiescent voltage of distribution network with synergic transmission and distribution
CN105811404B (en) * 2016-03-22 2018-01-19 山东大学 A kind of distribution static voltage stability situation monitoring method of transmission & distribution collaboration
CN106684859A (en) * 2016-11-24 2017-05-17 东南大学 Power grid short circuit current suppression method based on dynamic partition technique
CN106684859B (en) * 2016-11-24 2019-07-12 东南大学 A kind of grid short circuit electric current suppressing method based on Dynamic partitioning techniques
CN108631290A (en) * 2017-03-21 2018-10-09 国家电网公司 A kind of multi-infeed DC receiving end power grid direct current drop point screening technique
CN108493979A (en) * 2018-03-28 2018-09-04 国电南瑞科技股份有限公司 A kind of extra-high voltage direct-current access provincial power network partition zone optimizing progress control method
CN116937576A (en) * 2023-09-15 2023-10-24 国网山西省电力公司运城供电公司 Large-scale active power distribution network stability evaluation method, system and equipment
CN116937576B (en) * 2023-09-15 2023-11-24 国网山西省电力公司运城供电公司 Large-scale active power distribution network stability evaluation method, system and equipment

Also Published As

Publication number Publication date
CN104868478B (en) 2017-07-07

Similar Documents

Publication Publication Date Title
CN104868478B (en) A kind of method of the startup dynamic partition scheme under the power network state of emergency
CN103323742B (en) Based on the electric network fault recognition methods of random distributed points PMU
CN103901319B (en) Method for detecting transient voltage stability of power grid
CN106548265B (en) Power transmission network reliability assessment method based on cascading failure accident chain search
CN105207186B (en) A kind of distance protecting method for the transmission line of electricity containing THE UPFC
CN105917539A (en) Method and apparatus of reactor turn-to-turn protection
CN102403716A (en) Dynamic equalizing method for multi-infeed alternating/direct-current power grid
CN107104420B (en) I segment protection method of route distance suitable for THE UPFC access
CN103050933B (en) Based on the large-sized battery energy-accumulating power station interface protection method of one-terminal current Sudden Changing Rate
CN109066610B (en) Island power grid line fault positioning method
CN103870703A (en) Dynamic short-circuit ratio calculation method based on Thevenin equivalent parameter tracking
CN113203950B (en) Method and system for diagnosing ground fault degree of generator stator of marine nuclear power platform
CN103577698B (en) A kind of on-Line Voltage stabilizing determination method of alternating current-direct current combined hybrid system
Dobson New angles for monitoring areas
CN104993466A (en) Cascading failure fast dynamic simulation method suitable for alternating current-direct current power grid
CN105303454A (en) Vulnerability based power grid early-warning method
CN110472373A (en) A kind of dynamic frequency after Power System Disturbances estimates measurement method
CN107067126A (en) It is a kind of based on power flow transfer than thermally-stabilised key transmission channel recognition method
CN108092272A (en) A kind of voltage stabilization on-line monitoring method based on the Kalman filtering that fades
CN111884178A (en) Active power distribution network single-phase line break fault protection method and system
Zamani et al. A PMU-augmented stability power limit assessment for reliable arming of special protection systems
CN110783923B (en) Transformer substation main transformer neutral point connection processing method and system
CN104578008B (en) A kind of zero-sequence direction component faulty action preventing method
CN113315100A (en) Micro-grid protection method and system based on convolutional neural network
CN102522745A (en) Voltage-gradient-based steepest voltage collapse margin calculation method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170707