CN104333020B - A kind of real-time Low Frequency Oscillation Analysis of power system and optimum Corrective control method - Google Patents

Abstract

The invention relates to a method for real-time low-frequency oscillation analysis and optimal correction control of a power system. The real-time operation data report of the power grid is obtained through the data acquisition module and the unit parameter data report is formed, and the original data report and e format of the CIM model file are analyzed by the data analysis module. Analyze the original data report of the file; use the power flow program and low-frequency oscillation analysis program to perform power flow calculation and static stability calculation on the analyzed data report through the calculation execution module; if the static stability domain is insufficient, start the optimal correction control program for low-frequency oscillation analysis Adjust and correct; output the calculation result report to guide the operation personnel. The invention can make the establishment of the static stable operation mode of the power grid simple and easy, shorten the work flow and formulation time, reduce the work intensity and pressure of operating personnel, and comprehensively improve the voltage quality, safety and stability of the power grid operation.

Description

A Real-time Low-Frequency Oscillation Analysis and Optimal Correction Control Method of Power System

technical field

The invention relates to the field of low-frequency oscillation analysis and optimization of power systems, in particular to a power system real-time low-frequency oscillation analysis and optimal correction control method including optimal power flow constrained by damping ratio constraints and eigenvalue optimization.

Background technique

With the rapid growth of economy and population, people's demand for electric energy continues to increase, which makes the power system load continue to increase. Power system interconnection has many advantages, but is also prone to complex stability issues. At the same time, the power system has entered a new stage of large power grids, large units, and cross-regional networking, and these factors are likely to cause low-frequency oscillations. If the low-frequency oscillation is not effectively suppressed, it will cause more serious faults, large-scale power outages, and even system disconnection, resulting in huge economic losses.

The main harm of low-frequency oscillation to the power grid is that the out-of-control oscillation leads to out-of-synchronization of the power grid. At present, commonly used measures to suppress low-frequency oscillations include: installing series compensators to shorten the electrical distance, installing static var compensators (SVC) to provide dynamic voltage support, using DC transmission, installing unit PSS, DC PSS/PSD (power swing damping ) and other additional control functions, multi-DC coordinated control system and so on. However, the investment of these devices cannot guarantee that low-frequency oscillations will not occur during system operation. How to correctly predict whether low-frequency oscillations will occur during power system operation, and how to correct and control methods based on power generation rescheduling have become research hotspots at home and abroad in recent years.

Corrective control is through the rearrangement of controllable variables in the system. For the problem of improving security through correction control, there are two methods commonly used at present: sensitivity analysis method and optimization method. From the open literature, foreign research in this area began in 2004. The first type of method is based on the sensitivity calculation of the damping ratio in the ground state, and linearly adjusts the output of each generator to meet the requirements of small disturbance stability; The second type is the expected safety cost optimal power flow model with small disturbance stability constraints, which considers the operating conditions before the fault and a series of expected fault sets; the third type transforms the small disturbance stability constraints into matrix minimum singular value constraints , through some approximation techniques, the derivative of the smallest singular value is obtained for each iteration, and then integrated into the interior point algorithm to solve; the fourth type uses the first-order Taylor series expansion, and the first-order sensitivity of the active power output of the unit through the spectral abscissa function is explicit Express the small disturbance stability constraint and add it to the traditional OPF model to realize the solution process of the alternate calculation of the optimal power flow and the eigenvalue sensitivity of the small disturbance stability constraint. This method also considers the voltage stability under heavy load conditions. Its calculated results not only ensure the small disturbance stability, but also ensure the voltage stability. The calculations of several large-scale systems show that this method has fast calculation speed and good robustness. It has a wider application prospect.

Contents of the invention

The technical problem to be solved by the present invention is to provide a preventive adjustment measure for the low-frequency oscillation of the power system by using the method of optimal correction control aiming at whether the existing power system can return to a stable state after a period of time after being subjected to a small disturbance. To improve the stability of the evaluated power system, so that it is always in a safe and stable grid operation mode.

In order to solve the above technical problems, the present invention adopts the following technical solutions:

A real-time low-frequency oscillation analysis and optimal correction control system of a power system, including a data acquisition module, a data analysis module, a calculation execution module, an optimal correction control module, and a system management module; the data acquisition module obtains real-time raw data of CIM model files Reports, raw data reports of e-format files and formation of unit parameter data reports; the data analysis module converts the original data reports of CIM model files and the original data reports of e-format files into the format data reports required by the calculation program by using the custom data interface; calculation The execution module uses the power flow calculation program and the low-frequency oscillation analysis program to calculate the self-defined format data report; the optimal correction control module uses the optimal power flow with damping ratio constraints and the low-frequency oscillation correction control program optimized by the eigenvalue to calculate the original data report Carry out optimization and correction; the system management module completes the setting of relevant parameters.

A control method adapted to the above control system, comprising the following steps:

Step 1, obtaining real-time CIM model file raw data report, e-format file raw data report and unit parameter raw data report;

Step 2, use the data interface to analyze the original data report, and convert the original data report of the CIM model file, the original data report of the e-format file and the original data report of the unit parameters into the power flow program, the low-frequency oscillation analysis program and the optimal correction control program. custom format data report;

Step 3, use the power flow program to perform power flow calculation on the converted real-time format data report;

Step 4, based on the power flow calculation results, use the low-frequency oscillation analysis program to calculate the static stability of the system;

Step 5, check the static stability calculation result file, and judge whether there is an eigenvalue with a real part greater than zero or a value with an insufficient damping ratio in the result; if it does not exist, the system is considered to be statically stable; if it exists, the system is prone to low-frequency oscillation phenomenon, start the low-frequency oscillation analysis optimal correction control program;

Step 6: Start the optimal correction control program for low-frequency oscillation analysis, correct and set the control variables or constraint conditions of the calculated power system, and use the optimal power flow and eigenvalue optimization method with damping ratio constraints to analyze and solve the power system. Obtain correction control data; wherein, the control variables include the active power of the generator and the PV node voltage, and the constraints include the upper and lower limits of the active power and reactive power of the generator and the upper and lower limits of the node voltage;

Step 7, check the corrected calculation result file, and output the calculation result report, the result can ensure that the evaluated power system is in a safe and stable operation mode, and the power grid operation mode is the safe and stable operation limit of the evaluated power system;

Preferably, in step 1, the acquired data is the real-time operation data of the power grid;

Preferably, in step 2, while performing node analysis on the network structure and operating data including the evaluated power system, node equivalent processing is also performed on the network structure and operating data including the evaluated power system.

Preferably, in step 6, the constraints further include system power flow constraints, which are used to ensure that the evaluated power system always meets the system power flow constraints.

Preferably, in step 7, when the optimal calibration is completed, a calculation result report is output, which includes the optimized eigenvalue calculation results and the adjusted active power of the generator for reference and adjustment by the operator.

This program calculates the real-time operation data of the current power system to analyze whether the current system can maintain stability after experiencing small disturbances. If the system loses stability, optimize the control variables or constraints of the power system being evaluated, and evaluate The unit active output, unit reactive output, and PV node voltage are coordinated and controlled, and an optimal operation configuration strategy to ensure the safe and stable operation of the evaluated power system is proposed, which can effectively prevent the low-frequency oscillation of the evaluated power system, and then improve the safety and stability of the evaluated power system Operation level, to ensure the safe, stable and reliable operation of the power grid of the evaluated power system.

Compared with the prior art, the present invention has the following beneficial effects:

1. It can prevent the low-frequency oscillation phenomenon that may occur during the operation of the power system under evaluation, and can optimize the stability limit calculation of the power grid of the power system under evaluation, so that the power grid of the power system under evaluation can pass through The damping ratio constrained optimal power flow correction calculation program obtains the grid operation mode optimization data required for the static, safe and stable operation of the evaluated power system, and can use the grid operation mode optimization data for automatic calculation and verification, and the verification results can be directly used for guidance Compile the grid operation mode for static security and stability of the evaluated power system to ensure the safe, stable and reliable operation of the grid.

2. It can make the formulation of the static, safe and stable operation mode of the power grid simple and easy, shorten the work process and formulation time, reduce the work intensity and pressure of the operators, and comprehensively improve the voltage quality and safety and stability of the power grid operation.

Description of drawings

Fig. 1 is a structural block diagram of the power system real-time low-frequency oscillation analysis and optimal correction control system according to the present invention.

Fig. 2 is a schematic block diagram of the flow of the power system real-time low-frequency oscillation analysis and optimal correction control method according to the present invention.

detailed description

The terms involved in the present invention are interpreted as:

Low-frequency oscillation: refers to the vibration of the generator's rotor angle, speed, and related electrical quantities, such as line power, bus voltage, etc., with approximately equal amplitude or increased frequency between 0.1 and 2.5 Hz under the action of small disturbances.

Power flow calculation: refers to the calculation of the voltage, current magnitude and direction, and power distribution of the power system from the power source to the load under a certain operating mode and wiring mode.

Optimal Power Flow (OPF): When the structural parameters and load conditions of the system are given, adjust the available control variables (such as generator output power, adjustable transformer taps, etc.) The distribution of power flow under the operating constraints and making a certain performance index of the system (such as power generation cost or network loss) reach the optimal value.

Static stability: refers to the ability of the power system to automatically return to the initial operating state without spontaneous oscillation or non-periodic out-of-step after a small disturbance.

Static stability constrained optimal power flow: refers to the optimal power flow with static stability as the constraint under the constraints of the optimal power flow.

Static stability constrained optimal power flow aims at possible small disturbance unstable points, through optimization calculation, analyzes and gives the optimal operation scheme that fully meets the static stability requirements, so as to improve the safety and stability of power grid operation.

The static stability optimal correction control model is as follows:

1) The objective function to be optimized is: the sum of the adjusted generator active power is the minimum.

$\mathrm{<span\; class="notranslate">\; min</span>}\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; \&Element;</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; G</span>}}}{\mathrm{<span\; class="notranslate">\; \&Sigma;</span>}}\mathrm{<span\; class="notranslate">\; \&Delta;</span>}{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; Gi</span>}}$

Where: S _G is the set of generator nodes

2) Constraints include equality constraints and inequality constraints.

① Equality constraints:

a) The power flow equation of each node in the network.

$\left\{\begin{array}{c}{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; Gi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; Di</span>}}<span\; class="notranslate">\; -</span>\underset{{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; \&Element;</span>\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}}}{\mathrm{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; Y</span>}}_{\mathrm{<span\; class="notranslate">\; ij</span>}}{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}\mathrm{<span\; class="notranslate">\; cos</span>}{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; ij</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\\ {\mathrm{<span\; class="notranslate">\; Q</span>}}_{\mathrm{<span\; class="notranslate">\; Gi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; Q</span>}}_{\mathrm{<span\; class="notranslate">\; Di</span>}}<span\; class="notranslate">\; -</span>\underset{{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; \&Element;</span>\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}}}{\mathrm{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; Y</span>}}_{\mathrm{<span\; class="notranslate">\; ij</span>}}{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}\mathrm{<span\; class="notranslate">\; sin</span>}{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; ij</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\end{array}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; \&Element;</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}}\right.$

in:

V _i : the magnitude of the i-th node voltage;

Y _ij : element amplitude of node admittance matrix;

δ _ij = δ _i - δ _j - α _ij ;

δ _i : the phase angle of the i-th node voltage;

δ _j : the phase angle of the jth node voltage;

α _ij : element phase angle of the node admittance matrix;

S _N : system node set.

b) The initial value equation of the generator.

Taking the fourth-order model of the generator as an example, each generator is equipped with an IEEE-I type DC excitation system, and the following formula can be obtained:

$\left\{\begin{array}{c}{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; Gi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}\mathrm{<span\; class="notranslate">\; sin</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}\mathrm{<span\; class="notranslate">\; cos</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\\ {\mathrm{<span\; class="notranslate">\; Q</span>}}_{\mathrm{<span\; class="notranslate">\; Gi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}\mathrm{<span\; class="notranslate">\; cos</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}\mathrm{<span\; class="notranslate">\; sin</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\\ {\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}\mathrm{<span\; class="notranslate">\; sin</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; the\; si</span>}}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\\ {\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}\mathrm{<span\; class="notranslate">\; cos</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; the\; si</span>}}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\\ {\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; )</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\\ {\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; )</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 0</span>}\end{array}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; \&Element;</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; G</span>}}\right.$

in:

E' _di : Direct-axis transient electromotive force of generator i;

E' _qi : generator i quadrature axis transient electromotive force;

δ _i : rotation angle of generator i;

I _di : the direct-axis component of the stator current of generator i;

I _qi : the quadrature axis component of the stator current of generator i;

X' _di : direct-axis transient reactance of generator i;

X' _qi : quadrature axis transient reactance of generator i;

R _si : armature resistance of generator i;

E _fdi : output voltage of generator i excitation system.

c) Differential equations.

$\left\{\begin{array}{c}\frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; \&delta;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; \&omega;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&omega;</span>}}_{\mathrm{<span\; class="notranslate">\; the\; s</span>}}\\ \frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; \&omega;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span>\frac{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Mi</span>}}}{{\mathrm{<span\; class="notranslate">\; m</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}<span\; class="notranslate">\; -</span>\frac{<span\; class="notranslate">\; [</span>{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}<span\; class="notranslate">\; ]</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}}{{\mathrm{<span\; class="notranslate">\; m</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}<span\; class="notranslate">\; -</span>\frac{<span\; class="notranslate">\; [</span>{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}<span\; class="notranslate">\; ]</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}}{{\mathrm{<span\; class="notranslate">\; m</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}<span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; D.</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&omega;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; \&omega;</span>}}_{\mathrm{<span\; class="notranslate">\; the\; s</span>}}<span\; class="notranslate">\; )</span>}{{\mathrm{<span\; class="notranslate">\; m</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}\\ \frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}<span\; class="notranslate">\; -</span>\frac{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; )</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}<span\; class="notranslate">\; +</span>\frac{{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}\\ \frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; di</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}<span\; class="notranslate">\; +</span>\frac{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; )</span>{\mathrm{<span\; class="notranslate">\; I</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; qi</span>}}^{<span\; class="notranslate">\; \&prime;</span>}}\\ \frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; K</span>}}_{\mathrm{<span\; class="notranslate">\; Ei</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; E.</span>}}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}<span\; class="notranslate">\; )</span>}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Ei</span>}}}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}<span\; class="notranslate">\; +</span>\frac{{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; Ri</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Ei</span>}}}\\ \frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; Ri</span>}}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; Ri</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}}<span\; class="notranslate">\; +</span>\frac{{\mathrm{<span\; class="notranslate">\; K</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}}{\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}<span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; K</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}{\mathrm{<span\; class="notranslate">\; K</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}<span\; class="notranslate">\; +</span>\frac{{\mathrm{<span\; class="notranslate">\; K</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Ai</span>}}}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; refi</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span>\\ \frac{\mathrm{<span\; class="notranslate">\; d</span>}{\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}}{\mathrm{<span\; class="notranslate">\; dt</span>}}<span\; class="notranslate">\; =</span><span\; class="notranslate">\; -</span>\frac{{\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}}{{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}}<span\; class="notranslate">\; +</span>\frac{{\mathrm{<span\; class="notranslate">\; K</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}}{{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; Fi</span>}}<span\; class="notranslate">\; )</span>}^{\mathrm{<span\; class="notranslate">\; 2</span>}}}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; fdi</span>}}\end{array}\right.$

② Inequality constraints are network physical limitations and operational limitations (can be used in combination), including:

a) Node voltage limit, the value of the limit is related to voltage level, node type, area or node normal or emergency conditions.

${\underset{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; V</span>}}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&le;</span>{\stackrel{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; V</span>}}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; \&Element;</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}}$

in:

the upper limit of the amplitude of the i-th node voltage;

i: the lower limit of the amplitude of the i-th node voltage.

b) Unit output limitation, including active output limitation and reactive output limitation.

${\underset{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; P</span>}}}_{\mathrm{<span\; class="notranslate">\; Si</span>}}<span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; Si</span>}}<span\; class="notranslate">\; \&le;</span>{\stackrel{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; P</span>}}}_{\mathrm{<span\; class="notranslate">\; Si</span>}}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; \&Element;</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; SG</span>}}$

${\underset{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; Q</span>}}}_{\mathrm{<span\; class="notranslate">\; Ri</span>}}<span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; Q</span>}}_{\mathrm{<span\; class="notranslate">\; Ri</span>}}<span\; class="notranslate">\; \&le;</span>{\stackrel{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; Q</span>}}}_{\mathrm{<span\; class="notranslate">\; Ri</span>}}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; \&Element;</span>{\mathrm{<span\; class="notranslate">\; S</span>}}_{\mathrm{<span\; class="notranslate">\; R</span>}}$

in:

Generator active power upper and lower limits.

Generator reactive power upper and lower limits.

c) The real part of all eigenvalues of the equilibrium point state matrix is less than a given value:

{Re(λ)|λ∈λ(A)}≤ε (ε≤0)

d) Damping ratio constraint: the damping ratio of the eigenvalues of the state matrix is greater than a certain value (generally a value above 0.03)

η≥ξ (ξ≥0.03), where η is the damping ratio

The technical solutions of the present invention will be further specifically described below through the embodiments and in conjunction with the accompanying drawings.

Embodiment: A power system low-frequency oscillation analysis and optimal correction control method based on real-time data, as shown in FIG. 2 , the steps are as follows:

Step 1, obtaining real-time CIM model file raw data report, e-format file raw data report and unit parameter raw data report;

Step 2, use the data interface to analyze the original data report, and convert the original data report of the CIM model file, the original data report of the e-format file and the original data report of the unit parameters into the power flow program, the low-frequency oscillation analysis program and the optimal correction control program. custom format data report;

Step 3, use the power flow program to perform power flow calculation on the converted real-time format data report;

Step 4, based on the power flow calculation results, use the low-frequency oscillation analysis program to calculate the static stability of the system;

Step 5, check the static stability calculation result file, and judge whether there is an eigenvalue with a real part greater than zero or a value with an insufficient damping ratio in the result; if it does not exist, the system is considered to be statically stable; if it exists, the system is prone to low-frequency oscillation phenomenon, start the low-frequency oscillation analysis optimal correction control program;

Step 6, start the optimal correction control program for low-frequency oscillation analysis, and correct the control variables or constraint conditions of the calculated power system: use the optimal power flow and eigenvalue optimization method with damping ratio constraints to analyze and solve the power system, Obtain the correction control data, wherein the control variables include the active power of the generator, the PV node voltage, etc., and the constraints include the upper and lower limits of the active power and reactive power of the generator and the upper and lower limits of the node voltage;

Step 7, check the corrected calculation result file, and output the calculation result report, the result can ensure that the evaluated power system is in a safe and stable operation mode, and the power grid operation mode is the safe and stable operation limit of the evaluated power system;

Taking the 568 node system of a power grid in a certain area as an example, the real-time low-frequency oscillation operation analysis and optimal correction control of the system are carried out. The analysis steps are as follows:

1. The system has 568 nodes, 489 lines, 63 generators, and 313 transformers;

2. Obtain real-time CIM model file raw data report, e-format file raw data report and unit parameter raw data report;

3. Use the data interface to analyze the original data report, and convert the original data report of the CIM model file, the original data report of the e-format file and the original data report of the unit parameters into the power flow program, low-frequency oscillation analysis program and optimal correction control program. Custom format data report;

4. Use the power flow program to perform power flow calculations on the converted real-time format data reports;

5. Based on the calculation results of the power flow, the static stability calculation of the system is carried out by using the low-frequency oscillation analysis program;

6. Check the static stability calculation result file, and judge whether there are eigenvalues whose real part is greater than zero or values with insufficient damping ratio; if they do not exist, the system is considered to be statically stable; if they exist, the system is prone to low-frequency oscillations , start the low-frequency oscillation analysis optimal correction control program;

7. Start the optimal correction control program for low-frequency oscillation analysis, and correct and set the control variables or constraints of the calculated power system: use the optimal power flow and eigenvalue optimization method with damping ratio constraints to analyze and solve the power system, and obtain Correct the control data, wherein the control variables include generator active power, PV node voltage, etc., and constraints include the upper and lower limits of the active power and reactive power of the generator and the upper and lower limits of the node voltage;

8. Check the corrected calculation result file, coordinate and control the unit output of the system, and propose a configuration strategy for the static and stable operation of the calculated power system.

In the operation of the actual power system, the effect of this method on the stable operation analysis and correction control of the calculated power system is very obvious. Here is an example to illustrate:

Taking the static stability calculation results of the system as the base state, the optimal power flow is constrained by different damping ratios for correction control. The results are shown in the following table: where η represents the damping ratio.

serial number 1 2 3 4 damping ratio constraints Ground state η≥0.042 η≥0.044 η≥0.046 minimum damping ratio 0.04165170 0.04201111 0.04422099 0.04609023 P _G1 601.719 602.1612 576.8442 516.2662 P _G2 663.064 662.5029 634.6277 551.9978 P _G3 53.9325 53.9325 53.9325 53.9325 P _G4 488.689 488.3209 474.5019 452.3960 P _G5 491.162 491.2232 477.6311 464.5942 P _G6 59.0119 59.9563 76.0690 76.4700 P _G7 59.8597 61.5207 54.8828 59.6928 P _G8 276.2517 267.5292 198.3985 146.4452 P _G9 143.6 143.3400 115.2347 109.2728 P _G10 137.4036 135.9359 134.1355 130.8074 P _G11 206.0002 209.5390 225.2033 243.2707 P _G12 137.4036 139.3678 171.7857 210.0001 P _G13 74.0148 74.9380 130.3000 130.3000 P _G14 71.4248 70.7207 85.7741 130.3000 P _G15 636.774 657.2837 667.0000 667.0000 P _G16 633.343 614.6777 425.4178 321.8792 P _G17 605.852 600.7695 585.0646 591.6176 P _G18 598.646 595.2671 577.2408 583.2024 P _G19 252.175 256.4784 297.0540 333.4123 P _G20 256.526 256.4067 289.7756 336.8446 P _G21 81.1407 81.2795 87.4786 88.5387 P _G22 319.825 320.8706 328.9543 340.2272 P _G23 355.169 357.2485 398.7204 422.5860

P _G24 277.831 277.9717 273.9139 245.2894 P _G25 274.88 275.0514 277.9240 277.9083 P _G26 203.791 206.6434 235.0000 235.0000 P _G27 159.262 161.3318 187.8086 198.2184 P _G28 120.408 118.4469 94.4343 91.4639 P _G29 85.3499 84.8570 70.5627 65.9353 P _G30 252.104 248.3396 262.7101 284.4454 P _G31 263.721 268.1907 309.5787 353.141 P _G32 637.395 641.3828 709.5728 744.00 P _G33 100.021 99.0672 120.0154 139.1296 P _G34 29.2494 28.9245 44.7000 44.7000 P _G35 9.8574 9.9362 5.4071 0.2939 P _G36 119.041 119.1388 118.4691 135.0548 P _G37 129.625 131.0796 154.1303 158.80 P _G38 129.659 128.1881 109.5159 100.5506 P _G39 254.991 256.9994 260.6055 257.3609 P _G40 256.477 259.3306 288.7616 323.9847 P _G41 630.681 621.7782 532.8348 480.3806 P _G42 579.095 568.0007 486.8934 455.2222 P _G43 148.308 146.1131 151.4585 155.8615 P _G44 233.096 235.6788 275.7366 312.8406 P _G45 31.2908 35.4583 59.3691 83.9922 P _G46 983.917 986.5852 1028.9622 1016.7915 P _G47 81.3681 81.3678 77.0366 79.3562 P _G48 40.1094 40.1573 40.1572 40.1571 P _G49 82.5724 82.5663 82.5724 82.5724 P _G50 27.8106 28.1422 43.2837 45.6 P _G51 19.0347 22.5372 40.0997 45.6 P _G52 24.3893 24.8949 11.3887 0.0001

P _G53 30.1646 32.6000 32.6000 32.6000 P _G54 29.7373 26.6814 0 0 P _G55 16.734 18.8451 24.5000 24.50 P _G56 18.0163 14.1664 2.8436 6.5597 P _G57 29.5122 29.0492 24.2081 24.2236 P _G58 46.7769 46.7769 46.7769 46.7769 P _G59 24.9042 24.3922 17.1532 1.0022 P _G60 19.3787 19.7717 21.7400 21.7400 P _G61 9.4195 8.4041 12.7200 20.2718 P _G62 0 0 0 0 P _G63 14.2714 16.5109 15.5611 22.8659

Through the above example, it is illustrated that rationally adjusting the active output of the generating set can control the damping ratio of the eigenvalues calculated for the static stability of the calculated power system to be above a certain value, so that the static stability of the system will not cause low-frequency oscillation. Through the overall optimization of the system, the calculated power system can always maintain a static and stable operation mode, making the power grid operate reliably and stably.

The embodiments of the present invention are described above with reference to the accompanying drawings, but the implementation is not limited by the above embodiments, and those skilled in the art can make various changes or modifications within the scope of the appended claims.

Claims (5)

1. A real-time low-frequency oscillation analysis of a power system and a control method for an optimal correction control system, characterized in that, First set up a corresponding control system, the control system includes a data acquisition module, a data analysis module, a calculation execution module, an optimal correction control module and a system management module; the data acquisition module obtains real-time CIM model file raw data reports, The original data report of the e-format file and the original data report of the unit parameters; the data analysis module uses the custom data interface to convert the original data report of the CIM model file and the original data report of the e-format file into the format data report required by the calculation program; calculation execution The module uses the power flow calculation program and the low-frequency oscillation analysis program to calculate the self-defined format data report; the optimal correction control module uses the optimal power flow with damping ratio constraints and the low-frequency oscillation analysis optimal correction control program for the original Optimize and correct the data report; the system management module completes the setting of relevant parameters; Described control method comprises the following steps: Step 1, obtaining real-time CIM model file raw data report, e-format file raw data report and unit parameter raw data report; Step 2: Use the data interface to analyze the original data report, and convert the original data report of the CIM model file, the original data report of the e-format file, and the original data report of the unit parameters into a power flow calculation program, a low-frequency oscillation analysis program, and an optimal correction for low-frequency oscillation analysis Custom format data reports required by the control program; Step 3, use the power flow calculation program to perform power flow calculation on the converted real-time format data report; Step 4, based on the power flow calculation results, use the low-frequency oscillation analysis program to calculate the static stability of the system; Step 5, check the static stability calculation result file, and judge whether there is an eigenvalue with a real part greater than zero or a value with an insufficient damping ratio in the result; if it does not exist, the system is considered to be statically stable; if it exists, the system is prone to low-frequency oscillation phenomenon, start the low-frequency oscillation analysis optimal correction control program; Step 6: Start the optimal correction control program for low-frequency oscillation analysis, correct and set the control variables or constraint conditions of the calculated power system, and use the optimal power flow and eigenvalue optimization method with damping ratio constraints to analyze and solve the power system. Obtain correction control data; wherein, the control variables include the active power of the generator and the PV node voltage, and the constraints include the upper and lower limits of the active power and reactive power of the generator and the upper and lower limits of the node voltage; Step seven, check the corrected calculation result file, and output the calculation result report, the result can ensure that the evaluated power system is in a safe and stable operation mode, which is the safe and stable operation limit of the evaluated power system. 2. The control method according to claim 1, wherein in step 1, the acquired data is real-time operation data of the power grid. 3. The control method according to claim 1, characterized in that, in step 2, while performing node analysis on the network structure and operating data containing the evaluated power system, the network structure and operating data containing the evaluated power system are also analyzed The network structure and operation data are treated as node equivalents. 4. The control method according to claim 1, wherein in step 6, the constraints further include system power flow constraints, which are used to ensure that the evaluated power system always meets the system power flow constraints. 5. The control method according to claim 1, characterized in that, in step 7, when the optimal correction is completed, the output calculation result report contains the optimized characteristic value calculation result and the active power adjustment of the generator The power is adjusted for the operator's reference.