CN111049158A - Method and system for determining broadband oscillation stability of power system based on spectrum radius - Google Patents

Method and system for determining broadband oscillation stability of power system based on spectrum radius Download PDF

Info

Publication number
CN111049158A
CN111049158A CN201911293719.8A CN201911293719A CN111049158A CN 111049158 A CN111049158 A CN 111049158A CN 201911293719 A CN201911293719 A CN 201911293719A CN 111049158 A CN111049158 A CN 111049158A
Authority
CN
China
Prior art keywords
open
doubly
loop
power system
operation state
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
CN201911293719.8A
Other languages
Chinese (zh)
Other versions
CN111049158B (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.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201911293719.8A priority Critical patent/CN111049158B/en
Publication of CN111049158A publication Critical patent/CN111049158A/en
Application granted granted Critical
Publication of CN111049158B publication Critical patent/CN111049158B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Abstract

The invention discloses a method and a system for determining broadband oscillation stability of a power system based on spectrum radius, and belongs to the technical field of power system modeling. The method comprises the following steps: acquiring open-loop frequency characteristics of a power system of the doubly-fed wind turbine under the condition of open-loop operation; determining the open-loop frequency characteristic of the doubly-fed fan under the open-loop operation condition; determining an equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed wind turbine in the open-loop operation state; determining a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to an equivalent open-loop process equation; determining matrix elements of the spectrum radius function matrix, and generating a spectrum radius function curve according to the spectrum radius matrix elements; and judging the stability of the broadband oscillation of the electric power system of the doubly-fed fan in the closed-loop operation state according to the spectrum radius function curve characteristic. The method can be directly applied to the online analysis platform of the power system, and has important engineering application value.

Description

Method and system for determining broadband oscillation stability of power system based on spectrum radius
Technical Field
The invention relates to the technical field of power system modeling, in particular to a method and a system for determining broadband oscillation stability of a power system based on spectrum radius.
Background
New energy power generation represented by wind power and photovoltaic power generation has been widely applied to a large scale in an electric power system. With the increasing proportion of power electronic power supplies such as new power units, the problem of broadband oscillation caused by the complex interaction among power electronic equipment, traditional equipment and networks is increasingly prominent. The method for researching the grid-connected broadband oscillation analysis of the doubly-fed wind turbine has important theoretical significance and practical value for effectively solving the existing oscillation problem and ensuring the safe and stable operation of a power system containing high-proportion new energy.
At present, a great deal of research is carried out on a method for judging the grid-connected broadband oscillation stability of the doubly-fed wind turbine, and a mode analysis method and an impedance analysis method are mainly adopted. However, most of the existing analysis methods are relatively complex and lack of intuition, and the method for judging the grid-connected broadband oscillation stability of the doubly-fed fan is effectively and intuitively researched and has important theoretical practical value. Therefore, in order to visually and effectively determine the new energy broadband oscillation, an intuitive graphical new energy broadband oscillation stability determination method is needed.
Disclosure of Invention
In order to solve the above problems, the present invention provides a method for determining broadband oscillation stability of an electric power system based on a spectrum radius, comprising:
the method comprises the steps that an electric power system connected into the doubly-fed fan is taken as a target object, and the open-loop frequency characteristic of the electric power system under the open-loop operation condition of the doubly-fed fan is obtained;
determining the open-loop frequency characteristic of the doubly-fed wind turbine under the open-loop operation condition according to the open-loop frequency characteristic of the power system;
determining an equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed wind turbine in the open-loop operation state;
determining a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to an equivalent open-loop process equation;
determining matrix elements of the spectrum radius function matrix, and generating a spectrum radius function curve according to the spectrum radius matrix elements;
and judging the stability of the broadband oscillation of the electric power system of the doubly-fed fan in the closed-loop operation state according to the spectrum radius function curve characteristic.
Optionally, the equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state is as follows:
Figure BDA0002319931260000021
satisfies the following conditions:
Figure BDA0002319931260000022
the delta U is a voltage vector of a common connection point of the power system and the fan, and the delta I is a current vector of the common connection point of the power system and the fan;
i ═ 1 or 2, j ═ 1 or 2;
ΔU1and Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000023
to remove GsysMiddle gijThe ith row vector;
Figure BDA0002319931260000024
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000025
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000031
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method is an equivalent open-loop frequency characteristic matrix of the electric power system in the open-loop operation state of the doubly-fed fan.
Optionally, the matrix elements of the spectral radius function matrix are determined according to the following formula:
Figure BDA0002319931260000032
wherein the content of the first and second substances,
parameter pij(s) characterizing the broadband oscillation frequency domain characteristics of the system;
i 1 or 2, j 1 or 2, Δ U1And Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000033
to remove GsysMiddle gijThe ith row vector;
Figure BDA0002319931260000034
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000035
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000036
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method comprises the steps of obtaining an equivalent open-loop frequency characteristic matrix of an alternating current system in the open-loop operation state of the doubly-fed fan;
gwji(s) is equivalent open-loop frequency characteristic G of the doubly-fed fan in the open-loop operation state of the doubly-fed fanwRow j, column i, where s is a complex parameter indicating that the element is a function with respect to s;
ρ is a spectral radius function.
Optionally, the stability of the broadband oscillation of the power system of the doubly-fed wind turbine in the closed-loop operation state is judged according to the spectral radius function curve characteristic, and specifically:
and determining a characteristic value of a spectrum radius function curve according to the spectrum radius function curve, and if the characteristic value of the spectrum radius function curve is greater than 1, determining that the broadband oscillation of the electric power system is in a destabilization state under the closed-loop operation state of the doubly-fed fan.
The invention also provides a system for determining the broadband oscillation stability of the power system based on the spectrum radius, which comprises the following steps:
the first parameter acquisition module is used for acquiring the open-loop frequency characteristic of the electric power system under the open-loop operation condition of the doubly-fed fan by taking the electric power system accessed to the doubly-fed fan as a target object;
the second parameter acquisition module is used for determining the open-loop frequency characteristic of the doubly-fed fan under the open-loop operation condition according to the open-loop frequency characteristic of the power system;
the equivalent module is used for determining an equivalent open-loop process equation of the electric power system of the doubly-fed fan in a closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed fan in the open-loop operation state;
the calculation module is used for determining a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to an equivalent open-loop process equation;
the drawing module is used for determining matrix elements of the spectrum radius function matrix and generating a spectrum radius function curve according to the spectrum radius matrix elements;
and the judging module is used for judging the stability of the broadband oscillation of the electric power system of the double-fed fan in the closed-loop running state according to the spectrum radius function curve characteristic.
Optionally, the equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state is as follows:
Figure BDA0002319931260000041
satisfies the following conditions:
ΔU＝[ΔU1ΔU2]T
ΔI＝[ΔI1ΔI2]T
wherein, Delta U is a PCC point voltage vector connected with the fan, Delta I is a PCC point current vector connected with the fan,
Figure BDA0002319931260000042
For removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe resulting matrix after row i and column j,
Figure BDA0002319931260000043
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe resulting matrix after row i and column j,
Figure BDA0002319931260000044
to remove GsysMiddle gijThe ith row vector sum
Figure BDA0002319931260000045
To remove GsysMiddle gijThe next j-th column vector.
Optionally, the matrix elements of the spectral radius function matrix are determined according to the following formula:
Figure BDA0002319931260000046
wherein the content of the first and second substances,
parameter pij(s) characterizing the broadband oscillation frequency domain characteristics of the system;
i 1 or 2, j 1 or 2, Δ U1And Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000051
to remove GsysMiddle gijThe ith row vector;
Figure BDA0002319931260000052
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000053
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000054
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method comprises the steps of obtaining an equivalent open-loop frequency characteristic matrix of an alternating current system in the open-loop operation state of the doubly-fed fan;
gwji(s) is equivalent open-loop frequency characteristic G of the doubly-fed fan in the open-loop operation state of the doubly-fed fanwRow j, column i, where s is a complex parameter indicating that the element is a function with respect to s;
ρ is a spectral radius function.
Optionally, the stability of the broadband oscillation of the power system of the doubly-fed wind turbine in the closed-loop operation state is judged according to the spectral radius function curve characteristic, and specifically:
and determining a characteristic value of a spectrum radius function curve according to the spectrum radius function curve, and if the characteristic value of the spectrum radius function curve is greater than 1, determining that the broadband oscillation of the electric power system is in a destabilization state under the closed-loop operation state of the doubly-fed fan.
The method is more visual and clear in determining the double-fed wind power grid-connected broadband oscillation stability, meets the monitoring requirement of the power system, can be directly applied to an on-line analysis platform of the power system, and has important engineering application value.
Drawings
FIG. 1 is a flow chart of a method for determining broadband oscillation stability of a power system based on spectral radius according to the present invention;
FIG. 2 is a flowchart of an embodiment of a method for determining broadband oscillation stability of a power system based on spectral radius according to the present invention;
FIG. 3 is a diagram of a two-machine system model according to an embodiment of the present invention;
FIG. 4 is a frequency response curve of a spectrum radius matrix of a method for determining broadband oscillation stability of a power system based on spectrum radius according to the present invention;
FIG. 5 is a time domain simulation waveform of the fan stator voltage current based on the method of determining the broadband oscillation stability of the power system in accordance with the present invention;
FIG. 6 is a diagram of a system for determining broadband oscillation stability of a power system based on spectral radius according to the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the embodiments described herein, which are provided for complete and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, the same units/elements are denoted by the same reference numerals.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Further, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
The invention provides a method for determining broadband oscillation stability of a power system based on spectrum radius, as shown in figure 1, comprising the following steps:
the method comprises the steps that an electric power system connected into the doubly-fed fan is taken as a target object, and the open-loop frequency characteristic of the electric power system under the open-loop operation condition of the doubly-fed fan is obtained;
determining the open-loop frequency characteristic of the doubly-fed wind turbine under the open-loop operation condition according to the open-loop frequency characteristic of the power system;
determining an equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed wind turbine in the open-loop operation state;
the equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state is as follows:
Figure BDA0002319931260000071
satisfies the following conditions:
Figure BDA0002319931260000072
the delta U is a voltage vector of a common connection point of the power system and the fan, and the delta I is a current vector of the common connection point of the power system and the fan;
i ═ 1 or 2, j ═ 1 or 2;
ΔU1and Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000073
to remove GsysMiddle gijThe ith row vector;
Figure BDA0002319931260000074
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000075
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000076
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method is an equivalent open-loop frequency characteristic matrix of the electric power system in the open-loop operation state of the doubly-fed fan.
Determining a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to an equivalent open-loop process equation;
determining matrix elements of the spectrum radius function matrix, and generating a spectrum radius function curve according to the spectrum radius matrix elements;
determining matrix elements of a spectral radius function matrix according to the following formula:
Figure BDA0002319931260000077
wherein the content of the first and second substances,
parameter pij(s) characterizing the broadband oscillation frequency domain characteristics of the system;
i 1 or 2, j 1 or 2, Δ U1And Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000081
to remove GsysMiddle gijThe ith row vector;
Figure BDA0002319931260000082
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000083
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000084
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method comprises the steps of obtaining an equivalent open-loop frequency characteristic matrix of an alternating current system in the open-loop operation state of the doubly-fed fan;
gwji(s) is equivalent open-loop frequency characteristic G of the doubly-fed fan in the open-loop operation state of the doubly-fed fanwRow j, column i, where s is a complex parameter indicating that the element is a function with respect to s;
ρ is a spectral radius function.
The method for judging the stability of the broadband oscillation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state according to the spectrum radius function curve characteristic specifically comprises the following steps:
and determining a characteristic value of a spectrum radius function curve according to the spectrum radius function curve, and if the characteristic value of the spectrum radius function curve is greater than 1, determining that the broadband oscillation of the electric power system is in a destabilization state under the closed-loop operation state of the doubly-fed fan.
The present invention is further described in detail by taking the example shown in fig. 2;
building a simulation system shown in the figure 3 by using an electromagnetic transient software platform, adjusting the working state of the system, and judging the stability of the system model shown in the figure 3 in the current working state by using the method provided by the invention;
based on actual operating characteristics, the open-loop state space expression for the power system model in the DFIG open-loop operating state can be described as:
Figure BDA0002319931260000091
wherein, Δ XsysTo remove all state variables of the remaining subsystems after all fans, AsysIs a state matrix, BsysAs an output matrix, Csys、DsysThe output matrix and the direct transmission matrix are respectively, and the delta U and the delta I are respectively the PCC point voltage and current vector connected with the power system and the fan under the d q coordinate system, and satisfy the following conditions:
ΔU＝[ΔU1ΔU2]T
ΔI＝[ΔI1ΔI2]T
according to the formula (1), under the condition of open-loop operation of the DFIG, an open-loop frequency characteristic model of the power system under the condition of grid-connected open-loop operation of the doubly-fed wind turbine is as follows:
Figure BDA0002319931260000092
step two: according to the actual operation characteristics, the frequency characteristic model of the double-fed fan of the electric power system under the condition of the grid-connected open-loop operation of the double-fed fan is as follows:
Figure BDA0002319931260000093
in the formula, delta U represents a DFIG port voltage vector under an d q coordinate system and is an input vector; delta I is a DFIG port current vector under an d q coordinate system and is an output vector; a. thewIs a state matrix; b iswIs an output matrix; cw、DwRespectively an output matrix and a direct transmission matrix; Δ XwIs the state variable of DFIG.
The open-loop frequency characteristic function of the DFIG shown in fig. 2 is, according to equation (4):
Figure BDA0002319931260000094
the equivalent open-loop frequency characteristic model of the electric power system in the grid-connected closed-loop running state of the DFIG grid-connected double-fed fan is as follows:
Figure BDA0002319931260000095
wherein the content of the first and second substances,
Figure BDA0002319931260000096
indicates removal of GsysThe resulting matrix after row i and column j,
Figure BDA0002319931260000097
indicates removal of GwThe resulting matrix after row i and column j,
Figure BDA0002319931260000098
and
Figure BDA0002319931260000099
respectively represents removal of GsysMiddle gijThe ith row vector and the jth column vector.
According to the formula (6), determining a spectrum radius model of the power system under the grid-connected closed-loop operation state of the doubly-fed wind turbine:
Figure BDA0002319931260000101
according to the spectrum radius matrix of the DFIG grid-connected closed-loop system obtained by calculation of the formula (7), a frequency response curve of the spectrum radius of the closed-loop system in the current working state is drawn and shown in FIG. 3.
It can be seen from the characteristics of the spectrum radius curve that the frequency response curve of the closed-loop system spectrum radius matrix element crosses the straight line representing the system stability critical warning line in fig. 4, and the system is judged to be unstable in the current working state according to the characteristics of the spectrum radius matrix frequency response curve. In order to verify the analysis conclusion of the stability judging method, time domain simulation is carried out on the system in the current working state, and a voltage and current waveform diagram of the DFIG stator of the system is obtained and is shown in FIG. 5.
The present invention further provides a system 200 for determining broadband oscillation stability of an electric power system based on spectrum radius, as shown in fig. 6, including:
the first parameter obtaining module 201 is configured to obtain an open-loop frequency characteristic of an electric power system of a doubly-fed wind turbine under an open-loop operation condition, with the electric power system of the doubly-fed wind turbine as a target object;
the second parameter obtaining module 202 determines the open-loop frequency characteristic of the doubly-fed wind turbine under the open-loop operation condition according to the open-loop frequency characteristic of the power system;
the equivalent module 203 determines an equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed wind turbine in the open-loop operation state;
the calculation module 204 determines a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to the equivalent open-loop process equation;
the drawing module 205 determines matrix elements of the spectrum radius function matrix and generates a spectrum radius function curve according to the spectrum radius matrix elements;
the following formula determines the matrix elements of the spectral radius function matrix:
Figure BDA0002319931260000102
wherein the content of the first and second substances,
parameter pij(s) characterizing the broadband oscillation frequency domain characteristics of the system;
i 1 or 2, j 1 or 2, Δ U1And Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000111
to remove GsysMiddle gijThe ith row vector;
Figure BDA0002319931260000112
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000113
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000114
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method comprises the steps of obtaining an equivalent open-loop frequency characteristic matrix of an alternating current system in the open-loop operation state of the doubly-fed fan;
gwji(s) is equivalent open-loop frequency characteristic G of the doubly-fed fan in the open-loop operation state of the doubly-fed fanwRow j, column i, where s is a complex parameter indicating that the element is a function with respect to s;
ρ is a spectral radius function.
And the judging module 206 is used for judging the stability of the broadband oscillation of the power system of the doubly-fed fan in the closed-loop operation state according to the spectrum radius function curve characteristic.
The equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state is as follows:
Figure BDA0002319931260000115
satisfies the following conditions:
ΔU＝[ΔU1ΔU2]T
ΔI＝[ΔI1ΔI2]T
the delta U is a voltage vector of a common connection point of the power system and the fan, and the delta I is a current vector of the common connection point of the power system and the fan;
i ═ 1 or 2, j ═ 1 or 2;
ΔU1and Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure BDA0002319931260000121
to remove GsysMiddle gijLast ith rowVector quantity;
Figure BDA0002319931260000122
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure BDA0002319931260000123
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure BDA0002319931260000124
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method is an equivalent open-loop frequency characteristic matrix of the electric power system in the open-loop operation state of the doubly-fed fan.
The method for judging the stability of the broadband oscillation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state according to the spectrum radius function curve characteristic specifically comprises the following steps:
and determining a characteristic value of a spectrum radius function curve according to the spectrum radius function curve, and if the characteristic value of the spectrum radius function curve is greater than 1, determining that the broadband oscillation of the electric power system is in a destabilization state under the closed-loop operation state of the doubly-fed fan.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims (8)

1. A method of determining power system wide frequency oscillation stability based on spectral radius, the method comprising:
the method comprises the steps that an electric power system connected into the doubly-fed fan is taken as a target object, and the open-loop frequency characteristic of the electric power system under the open-loop operation condition of the doubly-fed fan is obtained;
determining the open-loop frequency characteristic of the doubly-fed wind turbine under the open-loop operation condition according to the open-loop frequency characteristic of the power system;
determining an equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed wind turbine in the open-loop operation state;
determining a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to an equivalent open-loop process equation;
determining matrix elements of the spectrum radius function matrix, and generating a spectrum radius function curve according to the spectrum radius matrix elements;
and judging the stability of the broadband oscillation of the electric power system of the doubly-fed fan in the closed-loop operation state according to the spectrum radius function curve characteristic.
2. The method of claim 1, wherein an equivalent open-loop process equation of the electric power system of the doubly-fed wind turbine in the closed-loop operation state is as follows:
Figure FDA0002319931250000011
satisfies the following conditions:
Figure FDA0002319931250000012
the delta U is a voltage vector of a common connection point of the power system and the fan, and the delta I is a current vector of the common connection point of the power system and the fan;
i ═ 1 or 2, j ═ 1 or 2;
ΔU1and Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure FDA0002319931250000021
to remove GsysMiddle gijThe ith row vector;
Figure FDA0002319931250000022
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure FDA0002319931250000023
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure FDA0002319931250000024
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method is an equivalent open-loop frequency characteristic matrix of the electric power system in the open-loop operation state of the doubly-fed fan.
3. The method of claim 1, determining matrix elements of the spectral radius function matrix according to the following formula:
Figure FDA0002319931250000025
wherein the content of the first and second substances,
parameter pij(s) characterizing the broadband oscillation frequency domain characteristics of the system;
i 1 or 2, j 1 or 2, Δ U1And Δ I1Voltage vector DeltaU and current vector DeltaI on d axisComponent of, Δ U2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure FDA0002319931250000026
to remove GsysMiddle gijThe ith row vector;
Figure FDA0002319931250000027
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure FDA0002319931250000028
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure FDA0002319931250000029
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method comprises the steps of obtaining an equivalent open-loop frequency characteristic matrix of an alternating current system in the open-loop operation state of the doubly-fed fan;
gwji(s) is equivalent open-loop frequency characteristic G of the doubly-fed fan in the open-loop operation state of the doubly-fed fanwRow j, column i, where s is a complex parameter indicating that the element is a function with respect to s;
ρ is a spectral radius function.
4. The method according to claim 1, wherein the characteristic of the spectrum radius function curve is used for judging the stability of the broadband oscillation of the power system of the doubly-fed wind turbine in the closed-loop operation state, and specifically comprises the following steps:
and determining a characteristic value of a spectrum radius function curve according to the spectrum radius function curve, and if the characteristic value of the spectrum radius function curve is greater than 1, determining that the broadband oscillation of the electric power system is in a destabilization state under the closed-loop operation state of the doubly-fed fan.
5. A system for determining wide frequency oscillation stability of a power system based on spectral radius, the system comprising:
the first parameter acquisition module is used for acquiring the open-loop frequency characteristic of the electric power system under the open-loop operation condition of the doubly-fed fan by taking the electric power system accessed to the doubly-fed fan as a target object;
the second parameter acquisition module is used for determining the open-loop frequency characteristic of the doubly-fed fan under the open-loop operation condition according to the open-loop frequency characteristic of the power system;
the equivalent module is used for determining an equivalent open-loop process equation of the electric power system of the doubly-fed fan in a closed-loop operation state according to the open-loop frequency characteristic of the doubly-fed fan in the open-loop operation state;
the calculation module is used for determining a spectrum radius function matrix of the electric power system of the doubly-fed wind turbine in a closed-loop operation state according to an equivalent open-loop process equation;
the drawing module is used for determining matrix elements of the spectrum radius function matrix and generating a spectrum radius function curve according to the spectrum radius matrix elements;
and the judging module is used for judging the stability of the broadband oscillation of the electric power system of the double-fed fan in the closed-loop running state according to the spectrum radius function curve characteristic.
6. The system of claim 5, wherein the equivalent open-loop process equation of the power system of the doubly-fed wind turbine in the closed-loop operation state is as follows:
Figure FDA0002319931250000031
satisfies the following conditions:
Figure FDA0002319931250000032
the delta U is a voltage vector of a common connection point of the power system and the fan, and the delta I is a current vector of the common connection point of the power system and the fan;
i ═ 1 or 2, j ═ 1 or 2;
ΔU1and Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure FDA0002319931250000041
to remove GsysMiddle gijThe ith row vector;
Figure FDA0002319931250000042
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure FDA0002319931250000043
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure FDA0002319931250000044
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method is an equivalent open-loop frequency characteristic matrix of the electric power system in the open-loop operation state of the doubly-fed fan.
7. The system of claim 5, the matrix elements of the spectral radius function matrix being determined according to the following formula:
Figure FDA0002319931250000045
wherein the content of the first and second substances,
parameter pij(s) characterizing the broadband oscillation frequency domain characteristics of the system;
i 1 or 2, j 1 or 2, Δ U1And Δ I1The components of the voltage vector DeltaU and the current vector DeltaI on the d-axis, DeltaU2And Δ I2The components of the voltage vector delta U and the current vector delta I on the q axis respectively;
gijis GsysRow i and column j in (1);
Figure FDA0002319931250000046
to remove GsysMiddle gijThe ith row vector;
Figure FDA0002319931250000047
representing and removing equivalent open-loop frequency characteristic G of the doubly-fed fan under the open-loop operation state of the doubly-fed fanwThe matrix obtained after the jth row and the ith column;
i is a unit array with a preset dimension;
Figure FDA0002319931250000048
for removing equivalent open-loop frequency characteristic G of electric power system in open-loop operating state of doubly-fed fansysThe matrix obtained after the ith row and the jth column;
Figure FDA0002319931250000049
to remove GsysMiddle gijThe subsequent jth column vector;
Gsysthe method comprises the steps of obtaining an equivalent open-loop frequency characteristic matrix of an alternating current system in the open-loop operation state of the doubly-fed fan;
gwji(s) is under the open-loop operation state of the doubly-fed fanEquivalent open-loop frequency characteristic G of doubly-fed fanwRow j, column i, where s is a complex variable indicating that the element is related tosA function of (a);
ρ is a spectral radius function.
8. The system of claim 5, wherein the spectral radius function curve characteristic is used for judging the stability of the broadband oscillation of the power system of the doubly-fed wind turbine in the closed-loop operation state, and specifically comprises:
and determining a characteristic value of a spectrum radius function curve according to the spectrum radius function curve, and if the characteristic value of the spectrum radius function curve is greater than 1, determining that the broadband oscillation of the electric power system is in a destabilization state under the closed-loop operation state of the doubly-fed fan.
CN201911293719.8A 2019-12-16 2019-12-16 Method and system for determining broadband oscillation stability of power system based on spectrum radius Active CN111049158B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911293719.8A CN111049158B (en) 2019-12-16 2019-12-16 Method and system for determining broadband oscillation stability of power system based on spectrum radius

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911293719.8A CN111049158B (en) 2019-12-16 2019-12-16 Method and system for determining broadband oscillation stability of power system based on spectrum radius

Publications (2)

Publication Number Publication Date
CN111049158A true CN111049158A (en) 2020-04-21
CN111049158B CN111049158B (en) 2022-07-01

Family

ID=70236656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911293719.8A Active CN111049158B (en) 2019-12-16 2019-12-16 Method and system for determining broadband oscillation stability of power system based on spectrum radius

Country Status (1)

Country Link
CN (1) CN111049158B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112636324A (en) * 2020-11-09 2021-04-09 上海电力大学 Stability judgment method and system for feeding multiple power electronic devices into complex weak power grid

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004088937A (en) * 2002-08-28 2004-03-18 Fuji Electric Fa Components & Systems Co Ltd Load machine analysis device for motor
CN1551647A (en) * 2003-04-30 2004-12-01 �ź㴫 Load control device and process with power control
CN109617096A (en) * 2018-12-27 2019-04-12 北京四方继保自动化股份有限公司 A kind of regional power grid broadband disturbance Method of Stability Analysis based on traversal impedance
CN110137941A (en) * 2019-04-10 2019-08-16 中国电力科学研究院有限公司 A kind of network stability control method and device based on double-fed blower stability influence

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004088937A (en) * 2002-08-28 2004-03-18 Fuji Electric Fa Components & Systems Co Ltd Load machine analysis device for motor
CN1551647A (en) * 2003-04-30 2004-12-01 �ź㴫 Load control device and process with power control
CN109617096A (en) * 2018-12-27 2019-04-12 北京四方继保自动化股份有限公司 A kind of regional power grid broadband disturbance Method of Stability Analysis based on traversal impedance
CN110137941A (en) * 2019-04-10 2019-08-16 中国电力科学研究院有限公司 A kind of network stability control method and device based on double-fed blower stability influence

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈骁等: "开环模式谐振条件下直驱风机接入引发电力系统宽频振荡的研究", 《中国电机工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112636324A (en) * 2020-11-09 2021-04-09 上海电力大学 Stability judgment method and system for feeding multiple power electronic devices into complex weak power grid
CN112636324B (en) * 2020-11-09 2023-09-26 上海电力大学 Stability judging method and system for feeding multiple power electronic devices into complex weak power grid

Also Published As

Publication number Publication date
CN111049158B (en) 2022-07-01

Similar Documents

Publication Publication Date Title
JP7194868B1 (en) Methods and apparatus for detecting yaw anomalies with respect to wind, and devices and storage media thereof
CN109217371B (en) Voltage source type converter grid-connected system stability analysis method, device and system considering phase-locked loop influence
CN112564117B (en) Voltage sag management method, voltage sag management device and terminal equipment
CN113708389B (en) Wind farm primary frequency modulation model parameter identification method and system based on actual power response
CN109723609A (en) A kind of fault early warning method and system of paddle change system of wind turbines
CN111049158B (en) Method and system for determining broadband oscillation stability of power system based on spectrum radius
CN110137948B (en) Method and device for evaluating sub-synchronous oscillation risk of doubly-fed wind power plant
CN113032968B (en) Double-fed fan controller parameter identification method and device
CN111987750B (en) Online monitoring method and system for safety phase-advancing capability margin of generator
Zheng et al. Transient stable region of power systems incorporating stochasticity of variable renewable energies and system contingencies
CN113868580A (en) Method for determining minimum peak regulation output of industrial steam supply working condition of extraction condensing heat supply unit
CN113094887A (en) Optimization method and device for frequency shift electromagnetic transient simulation and electronic equipment
KR20130047115A (en) Maximum power tracking emtdc model for wind generators depending on wind velocity
CN113471972B (en) Method and system for monitoring inertia of power system
CN116184828A (en) Online real-time optimization method and system for high-speed valve characteristic curve of self-adaptive steam turbine
CN115549093A (en) Method and system for online modeling and oscillation analysis of new energy power system
CN113794198B (en) Method, device, terminal and storage medium for suppressing broadband oscillation
CN115528684A (en) Ultra-short-term load prediction method and device and electronic equipment
CN111082454B (en) Method and system for obtaining oscillation risk quantitative index of power system
CN113722881A (en) Generator phase-advancing limit calculation method, device and medium based on parameter identification
CN115561668A (en) Method and system for testing limit short-circuit ratio of new energy power generation unit
CN107706948B (en) Multi-dimensional order controlled multi-step Taylor series transient stability analysis method
CN113742906B (en) Method for dynamically equalizing load of induction motor of large-scale industrial enterprise based on dominant dynamic similarity
CN114123336B (en) Control parameter identification method and system for direct-drive fan converter
CN113190982B (en) Method and system for determining key intermediate variables of converter control equipment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant