CN104240151A - Transient stability optimal correcting and control system and method for power system - Google Patents

Transient stability optimal correcting and control system and method for power system Download PDF

Info

Publication number
CN104240151A
CN104240151A CN201410460340.2A CN201410460340A CN104240151A CN 104240151 A CN104240151 A CN 104240151A CN 201410460340 A CN201410460340 A CN 201410460340A CN 104240151 A CN104240151 A CN 104240151A
Authority
CN
China
Prior art keywords
bpa
original data
msub
data report
transient
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
CN201410460340.2A
Other languages
Chinese (zh)
Other versions
CN104240151B (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.)
Guangxi University
Original Assignee
Guangxi 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 Guangxi University filed Critical Guangxi University
Priority to CN201410460340.2A priority Critical patent/CN104240151B/en
Publication of CN104240151A publication Critical patent/CN104240151A/en
Application granted granted Critical
Publication of CN104240151B publication Critical patent/CN104240151B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a transient stability optimal correcting and control system and method for a power system. The system comprises a data collecting module, a computation execution module, a data analysis module and a correction execution module, wherein the data collecting module obtains a BPA load flow format original data report and a BPA transient stability format original data report in a current operation mode; a BPA load flow computation program and a BPA transient stability computation program are used for computing the BPA load flow format original data report and the BPA transient stability format original data report; the BPA load flow format original data report and the BPA transient stability format original data report are analyzed by means of a self-defining data interface module; the original data reports are optimized and corrected by means of a modern interior point algorithm. By means of the system and method, programming of a transient security and stability operation mode of a power grid is easy, work procedures are reduced, programming time is shortened, labor intensity and pressure of workers are reduced, and voltage quality, safety and stability of operation of the power grid are comprehensively improved.

Description

Optimal correction control system and method for transient stability of power system
Technical Field
The invention relates to the field of transient stability analysis and optimization of a power system, in particular to a transient stability optimal correction control method of the power system based on BPA and transient stability constraint optimal power flow.
Background
In current competitive power market environments, the operation of a power system is often close to its steady limit state. Once a fault or disturbance occurs, a large-scale power failure accident and even system breakdown are often caused. Therefore, the power system cannot operate in the past conservative mode when the operation stability analysis (particularly transient stability analysis) of the power system and the optimal correction are found in the power market environment, and the transient safety of the system cannot be guaranteed under the condition of failure. This requires that the operation planning must optimize the operating conditions under the premise of satisfying the stability, i.e. adding stability constraints, especially transient stability constraints, to the optimal power flow.
Disclosure of Invention
The invention aims to solve the technical problem that aiming at the defect that the existing power system can not lead the evaluated power system to be in the transient stable power grid operation mode under the sudden failure, the invention provides the optimal correction control system and the optimal correction control method for the transient safety and stability of the power system, so as to improve the power grid transient stability of the evaluated power system, lead the power grid to be always in the transient safe and stable power grid operation mode, carry out the power flow optimization calculation on the power grid of the evaluated power system, realize that the power grid of the evaluated power system can obtain the power grid operation mode optimization data required by the transient safe and stable operation of the evaluated power system through a BPA system under the specified operation mode, automatically calculate and verify the power grid operation mode optimization data, and the verification result can be directly used for guiding the establishment of the transient safe and stable operation mode of the evaluated power system, so as to ensure the safe, stable and reliable operation of the power grid.
In order to solve the technical problems, the invention adopts the following technical scheme:
an optimal correction control system for transient stability of a power system comprises a data acquisition module, a calculation execution module, a data analysis module and a correction execution module. The data acquisition module is used for acquiring a BPA trend format original data report and a BPA transient stable format original data report in the current operation mode; the calculation execution module calculates the BPA power flow file original data report and the BPA transient stability original data report by adopting a BPA power flow calculation program and a BPA transient stability calculation program; the data analysis module is used for analyzing the BPA trend format original data report and the BPA transient stable format original data report by using the user-defined data interface module; and the correction execution module optimizes and corrects the original data report by adopting a modern interior point algorithm.
The control method adaptive to the control system comprises the following steps:
step one, providing a BPA trend format original data report and a BPA transient stability format original data report in an operation mode;
carrying out BPA power flow calculation on data in the BPA power flow format original data report by using BPA power system calculation analysis software;
setting a fault line, a fault occurrence time and a fault removal time, and performing BPA transient stability calculation on data in the BPA transient stability original data report by using BPA power system calculation and analysis software;
judging whether transient stability calculation of the BPA is unstable or not through the limit of the relative swing angle of the generator rotor; if the stable region is insufficient and the system is unstable, starting a transient stability constraint optimal power flow correction calculation program;
analyzing a BPA trend format original data report and a BPA transient stability format original data report containing the network structure and the operation data of the evaluated power system by using a custom data interface module;
starting a transient stability constraint optimal power flow correction calculation program, correcting and setting control variables or constraint conditions of the evaluated power system, and analyzing and solving the corrected and set evaluated power system by adopting an optimal power flow model and an optimal power flow algorithm, wherein the control variables comprise generator active power, generator reactive power and PV node voltage, and the constraint conditions comprise interconnection line constraint and node voltage constraint;
step seven, writing the corrected and set data back to the BPA trend format original data report;
step eight, outputting a BPA flow format optimized data report and a BPA transient stable format optimized data report which enable calculation convergence;
and step nine, covering the data in the BPA flow format original data report with the data in the BPA flow format optimized data report, namely ensuring that the evaluated power system is in a transient safe and stable operation mode, wherein the power grid operation mode is the transient safe and stable operation limit of the evaluated power system.
Preferably, in the fifth step, BPA network node analysis is performed on the network configuration and operation data including the power system under evaluation, and BPA network node equivalent processing is also performed on the network configuration and operation data including the power system under evaluation.
Preferably, in step six, the constraint condition further includes a system power flow constraint for ensuring that the evaluated power system always satisfies the system power flow constraint.
Preferably, in step six, the constraint condition further comprises a generator rotor relative swing angle limit constraint for ensuring that the generator is not unstable in synchronous operation.
Preferably, in the seventh step, if the BPA transient stability calculation is unstable, the original data report is corrected by using the optimal power flow model containing the transient stability constraint, and the data after the correction setting is written back to the BPA power flow format original data report.
According to the scheme, whether a current stabilizing system can be kept stable after experiencing the fault is judged by presetting the system fault, if the system is unstable, the control variable or constraint condition of the evaluated power system is optimized, the unit active power output, the unit reactive power output and the PV node voltage of the evaluated power system are coordinately controlled, an optimized operation configuration strategy for transient state safe and stable operation of the evaluated power system is provided, the transient state safe and stable operation level of the evaluated power system is further improved, and the safe, stable and reliable operation of a power grid of the evaluated power system is ensured. Compared with the prior art, the invention has the following beneficial effects:
1. the method can improve that the power grid of the evaluated power system is always in the transient state safe and stable power grid operation mode, and can perform power flow limit optimization calculation on the power grid of the evaluated power system, so that the power grid of the evaluated power system can obtain power grid operation mode optimization data required by the transient state safe and stable operation of the evaluated power system through a transient state stability constraint optimal power flow correction calculation program in any operation mode, and can perform automatic calculation and verification by using the power grid operation mode optimization data, and the verification result can be directly used for guiding and compiling the power grid operation mode of the transient state safe and stable of the evaluated power system so as to ensure the safe, stable and reliable operation of the power grid.
2. The method can make the establishment of the transient safe and stable operation mode of the power grid simple and easy, shorten the working process and the formulation time, reduce the working intensity and pressure of operators, and comprehensively improve the voltage quality and the safety and stability of the operation of the power grid.
Drawings
Fig. 1 is a block diagram of a transient stability optimal correction control system for an electric power system according to the present invention.
Fig. 2 is a schematic block diagram of a flow of the transient stability optimal correction control method of the power system according to the present invention.
FIG. 3 is a plot of rotor relative yaw before correction of the fault lines FUS 1M-FUS 2M in an exemplary embodiment.
FIG. 4 is a corrected rotor yaw relative curve for the exemplary embodiment of the fault lines FUS 1M-FUS 2M.
FIG. 5 is a rotor yaw relative curve of the embodiment before correction of fault line WENL 2M-LONGB 1M.
FIG. 6 is a corrected rotor yaw relative curve for the embodiment with fault line WENL 2M-LONGB 1M.
FIG. 7 is a rotor yaw curve of the fault line YUZ 1M-YUZ 2M before correction in an exemplary embodiment.
FIG. 8 is a corrected rotor yaw relative curve for fault line YUZ 1M-YUZ 2M in an exemplary embodiment.
FIG. 9 is a rotor yaw curve of the fault line DONGF 2M-LONGB 1M before calibration in an exemplary embodiment.
FIG. 10 is a corrected rotor yaw relative curve for the fault line DONGF 2M-LONGB 1M in an exemplary embodiment.
Fig. 11 is a rotor relative rocking curve of the fault line LUOD2M — length 2M before correction in the exemplary embodiment.
Fig. 12 is a corrected rotor relative rocking curve for the fault line LUOD 2M-long b2M in an example embodiment.
Detailed Description
The technical solution and the advantages of the present invention will be further explained with reference to the following embodiments and the accompanying drawings.
The model is roughly divided into two types, one is based on time domain simulation, a rotor motion equation is converted into an equivalent algebraic equation, and the algebraic equation and an angle criterion are used as constraints and are directly added into the optimal power flow model; or based on the Lyapunov direct method, the equation and the stability criterion are directly added. The former method has accurate model and accurate and reliable result, but has larger calculation amount and longer time consumption; while the latter results are somewhat conservative.
On the basis of the optimal power flow, the differential of the motion equation of the rotor of the generator under the fault is used as equality constraint by using an implicit trapezoidal integral rule, and the limit of the relative swing angle of the rotor is used as inequality constraint and added into the optimal power flow, so that an optimal power flow model of transient stability constraint is established. The objective is to solve an initial system state that exhibits transient stability in the event of a fault, while also meeting economic requirements. The optimal power flow of transient stability constraint belongs to a typical dynamic optimization problem, is also a research hotspot of an electric power system, and is solved by adopting a modern interior point algorithm.
The terms to which the invention relates are to be interpreted as follows:
BPA: the method refers to a power system calculation analysis software package introduced and developed by the Chinese institute of electrical science.
And (3) load flow calculation: the method is characterized in that under a certain operation mode and a certain wiring mode, the voltage, the current and the direction of the power system from a power supply to various places of a load and the distribution situation of the power are calculated.
Transient stability calculation: the method is used for calculating whether the power system can reach a new steady-state operation state or recover to an original state through a transient process after the power system is suddenly subjected to large interference under a certain operation condition. The large disturbance generally means a short-circuit fault, a sudden disconnection of a line, a generator, or the like.
Optimal power flow with transient stability constraints: when the structural parameters and the load condition of the system are given, the available control variables are adjusted, and transient stability constraint conditions are added to find the load flow distribution which can meet all the operation constraint conditions and enable a certain performance index of the system to reach the optimal value.
The mathematical model used in the invention is as follows:
1. objective function
Taking the minimum active power output correction quantity of the generator as an objective function:
<math><mrow> <mi>min</mi> <mi> f</mi> <mrow> <mo>(</mo> <mo>&CenterDot;</mo> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mrow> </munder> <msup> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>G</mi> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow></math>
wherein, PG1iFor the current generator i active power, PGiTo adjust the active power of the rear generator i.
2. Equality constraint
1) Node injection power flow balance equation
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>Di</mi> </msub> <mo>-</mo> <msub> <mi>e</mi> <mi>i</mi> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>ij</mi> </msub> <msub> <mi>e</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>B</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>B</mi> <mi>ij</mi> </msub> <msub> <mi>e</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Q</mi> <mi>Ri</mi> </msub> <mo>-</mo> <msub> <mi>Q</mi> <mi>Di</mi> </msub> <mo>-</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>ij</mi> </msub> <msub> <mi>e</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mi>B</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>e</mi> <mi>i</mi> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>B</mi> <mi>ij</mi> </msub> <msub> <mi>e</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>n</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow></math>
Wherein,
n: the number of nodes;
PDi: the active power of load i;
QRi: reactive power of a reactive power source i;
QDi: the reactive power of the load is;
ei,fi: the real part and the imaginary part of the voltage of the node i;
ej,fj: the voltage amplitude and phase angle of node j;
Gij,Bij: the node admittance real part and the imaginary part of the branch ij;
Sn: a system node set;
2) equation of motion of generator rotor
Differentiating the rotor motion equation by using an implicit trapezoidal integration method to obtain a rotor equation under an expected fault:
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mfrac> <mi>&Delta;t</mi> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msubsup> <mi>&omega;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&omega;</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> <msub> <mi>&omega;</mi> <mi>N</mi> </msub> <mo>=</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&omega;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&omega;</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>-</mo> <mfrac> <mi>&Delta;t</mi> <mrow> <mn>2</mn> <msub> <mi>M</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> <msubsup> <mi>&omega;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>+</mo> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mi>ei</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> <msubsup> <mi>&omega;</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <mo>-</mo> <msubsup> <mi>P</mi> <mi>ei</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> <mo>,</mo> <mi>t</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>t</mi> </msub> <mo>,</mo> <mi>k</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow></math>
wherein, <math><mrow> <msubsup> <mi>P</mi> <mi>ei</mi> <mi>t</mi> </msubsup> <mo>=</mo> <msub> <mi>E</mi> <mi>i</mi> </msub> <munder> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mrow> </munder> <msub> <mi>E</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>ij</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&delta;</mi> <mi>j</mi> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>B</mi> <mi>ij</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&delta;</mi> <mi>j</mi> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>;</mo> </mrow></math>
k: predicting the number of faults;
t: a time period;
ωN: synchronizing the rotating speed;
Δ t: numerical integration step length;
Di: generator i damping coefficient;
: the rotation angle of the generator i at each time interval after the fault is removed;
: the speed of the generator i at each time interval after the fault is removed;
: the electromagnetic power of the generator i at each time interval after the fault is removed;
Ei: the potential amplitude of the generator i;
Ej: the magnitude of the potential of generator j;
Gij: only the real part of the existing admittance array of the generator node is contained after the fault;
Bij: only containing the virtual part of the existing admittance array of the generator node after the fault;
St: a set of time periods.
3) Initial value equation of generator
In steady state operation, the generator potential satisfies the formula:
<math><mrow> <msub> <mi>E</mi> <mi>i</mi> </msub> <mo>&angle;</mo> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mn>0</mn> </msubsup> <mo>=</mo> <msub> <mover> <mi>V</mi> <mo>&CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>+</mo> <mi>j</mi> <mfrac> <mrow> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>Q</mi> <mi>Gi</mi> </msub> </mrow> <msub> <mover> <mover> <mi>V</mi> <mo>&CenterDot;</mo> </mover> <mo>^</mo> </mover> <mi>i</mi> </msub> </mfrac> <msubsup> <mi>X</mi> <mi>di</mi> <mo>&prime;</mo> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow></math>
namely:
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>E</mi> <mi>i</mi> </msub> <msub> <mi>e</mi> <mi>i</mi> </msub> <mi>sin</mi> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msub> <mi>E</mi> <mi>i</mi> </msub> <msub> <mi>f</mi> <mi>i</mi> </msub> <mi>cos</mi> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <msubsup> <mi>X</mi> <mi>di</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>e</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>f</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>E</mi> <mi>i</mi> </msub> <msub> <mi>e</mi> <mi>i</mi> </msub> <mi>cos</mi> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msub> <mi>E</mi> <mi>i</mi> </msub> <msub> <mi>f</mi> <mi>i</mi> </msub> <mi>sin</mi> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mn>0</mn> </msubsup> <mo>+</mo> <msub> <mi>Q</mi> <mi>Gi</mi> </msub> <msubsup> <mi>X</mi> <mi>di</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow></math>
wherein,
Ei: the potential amplitude of the generator i;
QGi: reactive power of generator i;
X′d: transient reactance of the generator;
ei,fi: the real part and the imaginary part of the end voltage of the generator i;
: and (5) issuing the angle of the generator i in an initial state.
3. Inequality constraint condition
1) Operating constraints
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <munder> <mi>P</mi> <mo>&OverBar;</mo> </munder> <mi>Gi</mi> </msub> <mo>&le;</mo> <msub> <mi>P</mi> <mi>Gi</mi> </msub> <mo>&le;</mo> <msub> <mover> <mi>P</mi> <mo>&OverBar;</mo> </mover> <mi>Gi</mi> </msub> </mtd> <mtd> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <munder> <mi>Q</mi> <mo>&OverBar;</mo> </munder> <mi>Ri</mi> </msub> <mo>&le;</mo> <msub> <mi>Q</mi> <mi>Ri</mi> </msub> <mo>&le;</mo> <msub> <mover> <mi>Q</mi> <mo>&OverBar;</mo> </mover> <mi>Ri</mi> </msub> </mtd> <mtd> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>R</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <munder> <mi>V</mi> <mo>&OverBar;</mo> </munder> <mi>i</mi> </msub> <mo>&le;</mo> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>&le;</mo> <msub> <mover> <mi>V</mi> <mo>&OverBar;</mo> </mover> <mi>i</mi> </msub> </mtd> <mtd> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow></math>
Wherein,
the upper limit of the active power of the generator i;
P Gi: the active power lower limit of the generator i;
the upper limit of the reactive power source i;
Q Ri: the lower reactive power limit of the reactive power source i;
the upper voltage limit of node i;
V i: the lower voltage limit of node i;
SR: a reactive power supply set;
Sn: a system node set;
2) transient stability constraint
Taking the inertia center COI of the system as a reference, the relative swing angle limit of each generator rotor is as follows:
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <munder> <mi>&delta;</mi> <mo>&OverBar;</mo> </munder> <mo>&le;</mo> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&delta;</mi> <mi>COI</mi> <mn>0</mn> </msubsup> <mo>&le;</mo> <mover> <mi>&delta;</mi> <mo>&OverBar;</mo> </mover> </mtd> </mtr> <mtr> <mtd> <munder> <mi>&delta;</mi> <mo>&OverBar;</mo> </munder> <mo>&le;</mo> <msubsup> <mi>&delta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&delta;</mi> <mi>COI</mi> <mi>t</mi> </msubsup> <mo>&le;</mo> <mover> <mi>&delta;</mi> <mo>&OverBar;</mo> </mover> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>i</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> <mo>,</mo> <mi>t</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow></math>
wherein, <math><mrow> <msubsup> <mi>&delta;</mi> <mi>COI</mi> <mn>0</mn> </msubsup> <mo>=</mo> <mfrac> <mrow> <munder> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mrow> </munder> <msubsup> <mi>&delta;</mi> <mi>j</mi> <mn>0</mn> </msubsup> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> <mrow> <munder> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mrow> </munder> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>,</mo> <msubsup> <mi>&delta;</mi> <mi>COI</mi> <mi>t</mi> </msubsup> <mo>=</mo> <mfrac> <mrow> <munder> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mrow> </munder> <msubsup> <mi>&delta;</mi> <mi>j</mi> <mi>t</mi> </msubsup> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> <mrow> <munder> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>&Element;</mo> <msub> <mi>S</mi> <mi>G</mi> </msub> </mrow> </munder> <msub> <mi>M</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>;</mo> </mrow></math>
issuing a generator inertia center angle in an initial state;
predicting the inertia center angle of each time interval after the fault is removed;
Mj: the inertia time constant of the jth generator;
δ: a lower limit of the generator rotor angle;
and (4) upper limit of the angle of the generator rotor.
Referring to fig. 1 and 2, the present embodiment is implemented according to the following specific method: setting an expected system fault, a fault occurrence time and a fault removal time in the current stable state of the system through a power system network structure and a load flow format original data report and a transient stability original data report of operation data, and performing BPA load flow calculation and transient stability calculation on data in the original data report; if oscillation occurs in BPA transient stability calculation, correcting and setting control variables or constraint conditions by using a transient stability constraint optimal power flow correction calculation program, and writing corrected data back to a BPA power flow format original data report; outputting a BPA flow format optimized data report and a BPA transient stable format optimized data report when the calculation is converged; and covering the data in the BPA flow format original data report with the data in the BPA flow format optimized data report, so that the evaluated power system can be ensured to be in a transient safe and stable operation mode, wherein the power grid operation mode is the transient safe and stable operation limit of the evaluated power system.
Now, by comprehensively testing a 195-node system of a power grid in a certain area, the minimum active power output correction quantity of the generator after a three-phase short-circuit fault occurs in the middle section of a line of 141 lines is obtained, and a test result of 5 representative lines is given. The analysis steps are as follows:
1. the evaluated power system comprises 195 nodes, 141 lines and 29 generators, and BPA load flow calculation is carried out;
2. setting an expected system fault, and performing BPA transient stability calculation;
3. and if the transient stability calculation cannot be kept stable, performing transient stability constraint optimal power flow optimal correction, and simultaneously setting equivalent control variables and constraint conditions for joint optimization, wherein the algorithm of the joint optimization adopts a modern interior point algorithm.
4. When the power flow calculation is converged and the transient stability calculation is stable, outputting a BPA power flow format optimization data report and a BPA transient stability format optimization data report when the calculation is converged;
5. and covering the data in the BPA flow format original data report with the data in the BPA flow format optimized data report, so that the evaluated power system can be ensured to be in a transient safe and stable operation mode, the power grid operation mode is an optimal transient safe and stable operation mode of the evaluated power system, and a calculation analysis report is output.
In the transient safe and stable operation optimization of the actual power system, the method has very obvious effect on the transient safe and stable operation of the evaluated power system. The following is illustrated by way of example:
2. test result with minimum active output correction of generator as objective function
To show the characteristics of the algorithm, the integral step length is selected to be 0.02s, the minimum generator active output correction is calculated according to the objective function of the minimum generator active output correction, and the method is described in detail through the test results of 5 representative lines. Table 1 is a table of test data with integration step size of 0.02 s.
Table 1195 node system 5 test data table with representative line target of minimum active correction of generator
Reference capacity: 100MVA
The test result curves are shown in fig. 3-12, wherein fig. 3 and 4 are the relative swing curves of the rotor before and after the calibration of the fault lines FUS 1M-FUS 2M, respectively; FIGS. 5 and 6 are the rotor yaw curves before and after correction of fault line WENL 2M-LONGB 1M, respectively; FIGS. 7 and 8 are graphs of rotor yaw before and after calibration for fault lines YUZ 1M-YUZ 2M, respectively; FIGS. 9 and 10 are the rotor relative rocking curves before and after calibration of the fault line DONGF 2M-LONGB 1M, respectively; fig. 11 and 12 are rotor relative rocking curves before and after correction of the fault line LUOD 2M-LONGB 2M, respectively.
And (3) calculating a curve before correction by transient stability, keeping the current operation state of the system, generating an expected fault at 0s, and cutting off the fault at 0.15s to obtain a simulation curve. And the corrected curve is obtained by calculating an objective function taking the active power output correction quantity of the generator as the minimum. The system changes the current operation state, an expected fault occurs at 0s, and the fault is removed in 0.15 s.
As can be seen from the results of table 1 and fig. 3 to fig. 12, the optimal power flow calculation module with transient stability constraint can effectively perform preventive control on the system. The active output of the generator set is reasonably adjusted; the control of each line of the evaluated power system can meet the safety and stability guide rule of the power system and the voltage quality and reactive power management regulation of the power system. Transient stability of the evaluated power system is effectively controlled.
The embodiments of the present invention have been described above with reference to the accompanying drawings, but the implementation is not limited to the above-described embodiments, and those skilled in the art can make various changes or modifications within the scope of the appended claims.

Claims (6)

1. An optimal correction control system for transient stability of a power system comprises a data acquisition module, a calculation execution module, a data analysis module and a correction execution module; the data acquisition module is used for acquiring a BPA trend format original data report and a BPA transient stable format original data report in the current operation mode; the calculation execution module is used for calculating the BPA power flow file original data report and the BPA transient stability original data report by adopting a BPA power flow calculation program and a BPA transient stability calculation program; the data analysis module is used for analyzing the BPA trend format original data report and the BPA transient stable format original data report by using the user-defined data interface module; and the correction execution module optimizes and corrects the original data report by adopting a modern interior point algorithm.
2. A control method adapted to the power system transient stability optimal correction control system of claim 1, comprising the steps of:
giving a BPA trend format original data report and a BPA transient stable format original data report in an operation mode;
carrying out BPA power flow calculation on data in the BPA power flow format original data report by using BPA power system calculation analysis software;
setting a fault line, a fault occurrence time and a fault removal time, and performing BPA transient stability calculation on data in the BPA transient stability original data report by using BPA power system calculation and analysis software;
judging whether transient stability calculation of the BPA is unstable or not through the limit of the relative swing angle of the generator rotor; when the stable region is insufficient and the system is unstable, starting a transient stability constraint optimal power flow correction calculation program;
analyzing a BPA trend format original data report and a BPA transient stability format original data report containing the network structure and the operation data of the evaluated power system by using a custom data interface module;
starting a transient stability constraint optimal power flow correction calculation program, correcting and setting control variables or constraint conditions of the evaluated power system, and analyzing and solving the corrected and set evaluated power system by adopting an optimal power flow model and an optimal power flow algorithm, wherein the control variables comprise generator active power, generator reactive power and PV node voltage, and the constraint conditions comprise interconnection line constraint and node voltage constraint;
step seven, writing the corrected and set data back to the BPA trend format original data report;
step eight, outputting a BPA flow format optimized data report and a BPA transient stable format optimized data report when calculation convergence occurs;
and step nine, covering the data in the BPA flow format original data report with the data in the BPA flow format optimized data report, namely ensuring that the evaluated power system is in a transient safe and stable operation mode, wherein the power grid operation mode is the transient safe and stable operation limit of the evaluated power system.
3. The control method according to claim 2, wherein in the fifth step, BPA network node analysis is performed on the network structure and the operation data including the evaluated power system, and BPA network node equivalent processing is also performed on the network structure and the operation data including the evaluated power system.
4. The control method according to claim 2, wherein in step six, the constraint conditions further include a system power flow constraint for ensuring that the evaluated power system always satisfies the system power flow constraint.
5. The control method of claim 2, wherein in step six, the constraints further include generator rotor relative yaw angle limit constraints for ensuring that the generator does not destabilize during synchronous operation.
6. The control method according to claim 2, wherein in step seven, when the BPA transient stability calculation is unstable, the original data report is corrected by using an optimal power flow model containing transient stability constraints, and the corrected data is written back to the BPA power flow format original data report.
CN201410460340.2A 2014-09-11 2014-09-11 A kind of optimal corrective control of electric power system transient stability and method Expired - Fee Related CN104240151B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410460340.2A CN104240151B (en) 2014-09-11 2014-09-11 A kind of optimal corrective control of electric power system transient stability and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410460340.2A CN104240151B (en) 2014-09-11 2014-09-11 A kind of optimal corrective control of electric power system transient stability and method

Publications (2)

Publication Number Publication Date
CN104240151A true CN104240151A (en) 2014-12-24
CN104240151B CN104240151B (en) 2017-12-08

Family

ID=52228164

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410460340.2A Expired - Fee Related CN104240151B (en) 2014-09-11 2014-09-11 A kind of optimal corrective control of electric power system transient stability and method

Country Status (1)

Country Link
CN (1) CN104240151B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105005831A (en) * 2015-08-06 2015-10-28 河海大学 Quasi direct current dynamic optimal power flow based on power system power flow coupling relation
CN105046588A (en) * 2015-08-13 2015-11-11 河海大学 Improved DC (Direct Current) dynamic optimal power flow calculating method based on network loss iteration
CN106026083A (en) * 2016-06-27 2016-10-12 吉林大学 Power grid transient stabilization analysis method based on MATLAB
CN106972494A (en) * 2017-05-18 2017-07-21 广西大学 A kind of AC-DC hybrid power grid sends section conveying power optimization method outside

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100114400A1 (en) * 2008-11-05 2010-05-06 Abb Research Ltd. Voltage Regulation Optimization
CN103066595A (en) * 2012-12-26 2013-04-24 中国电力科学研究院 Optimization method of extra-high voltage transient stability control
CN103746368A (en) * 2013-04-10 2014-04-23 广西大学 Method of optimizing static safe and stable operation limit of electric power system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100114400A1 (en) * 2008-11-05 2010-05-06 Abb Research Ltd. Voltage Regulation Optimization
CN103066595A (en) * 2012-12-26 2013-04-24 中国电力科学研究院 Optimization method of extra-high voltage transient stability control
CN103746368A (en) * 2013-04-10 2014-04-23 广西大学 Method of optimizing static safe and stable operation limit of electric power system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
夏小琴等: "暂态稳定约束最优潮流的改进降阶内点算法", 《电工技术学报》 *
苏维波等: "一种基于BPA数据的静态安全最优校正分析方法", 《电网与清洁能源》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105005831A (en) * 2015-08-06 2015-10-28 河海大学 Quasi direct current dynamic optimal power flow based on power system power flow coupling relation
CN105005831B (en) * 2015-08-06 2019-04-30 河海大学 A kind of calculation method of the quasi- direct current Dynamic Optimal Power Flow Problem based on electric power system tide coupled relation
CN105046588A (en) * 2015-08-13 2015-11-11 河海大学 Improved DC (Direct Current) dynamic optimal power flow calculating method based on network loss iteration
CN105046588B (en) * 2015-08-13 2019-04-30 河海大学 A kind of calculation method of the improvement direct current Dynamic Optimal Power Flow Problem based on network loss iteration
CN106026083A (en) * 2016-06-27 2016-10-12 吉林大学 Power grid transient stabilization analysis method based on MATLAB
CN106026083B (en) * 2016-06-27 2019-02-19 吉林大学 A kind of Power Network Transient Stability analysis method based on MATLAB
CN106972494A (en) * 2017-05-18 2017-07-21 广西大学 A kind of AC-DC hybrid power grid sends section conveying power optimization method outside

Also Published As

Publication number Publication date
CN104240151B (en) 2017-12-08

Similar Documents

Publication Publication Date Title
Feng et al. A comprehensive approach for preventive and corrective control to mitigate voltage collapse
Sarasúa et al. Dynamic response and governor tuning of a long penstock pumped-storage hydropower plant equipped with a pump-turbine and a doubly fed induction generator
CN102682358B (en) A kind of assessment wind-electricity integration scale and Net Frame of Electric Network adaptive planning simulation method
Kalcon et al. Small-signal stability analysis of multi-terminal VSC-based DC transmission systems
Zhang et al. Coordinated supplementary damping control of DFIG and PSS to suppress inter-area oscillations with optimally controlled plant dynamics
Xi et al. Torsional oscillation damping control for DFIG-based wind farm participating in power system frequency regulation
CN104240151B (en) A kind of optimal corrective control of electric power system transient stability and method
Du et al. Frequency regulation analysis of modern power systems using start-stop peak shaving and deep peak shaving under different wind power penetrations
CN110266039A (en) Consider wind power fluctuation and cuts the governor parameter optimization method of machine-cut load disturbance frequency modulation next time
CN104485670B (en) The control method of voltage sensitivity industrial load time-varying damping characteristic in island network
Wang et al. Small signal stability analysis considering grid-connected wind farms of DFIG type
Martínez–Lucas et al. Combined hydro-wind frequency control scheme: Modal analysis and isolated power system case example
Zhang et al. Online frequency security assessment based on analytical model considering limiting modules
CN104113086B (en) A kind of wind-powered electricity generation thermoelectricity cuts machine policy co-ordination optimization method
Husham et al. Decentralized stability enhancement of DFIG-based wind farms in large power systems: Koopman theoretic approach
CN104767207A (en) Voltage stabilization prevention and control method based on reactive storage sensitivity of generators
Moreira et al. Identification of dynamic simulation models for variable speed pumped storage power plants
CN115940157A (en) Method, device and equipment for automatically generating load flow scene of stability control strategy checking task
Riquelme et al. RoCoF-Minimizing H₂ Norm Control Strategy for Multi-Wind Turbine Synthetic Inertia
CN105186440A (en) Relaying protection fixed value setting method based on electromechanical transient method
KR102348138B1 (en) Apparatus for determining stability margin of online generator special protection system and method thereof
Deshmukh et al. Review on stability analysis of grid connected wind power generating system
CN104638621B (en) Reclosing and emergent control comprehensive optimization method under a kind of power network based on DSR
Li et al. Adaptive robust backstepping design for turbine valve controller
CN106374509B (en) A kind of method of large-scale photovoltaic electricity generation grid-connecting system shortsightedness voltage oscillation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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

Granted publication date: 20171208