CN115085224A - New energy power system frequency stability rapid analysis method - Google Patents

New energy power system frequency stability rapid analysis method Download PDF

Info

Publication number
CN115085224A
CN115085224A CN202210995878.8A CN202210995878A CN115085224A CN 115085224 A CN115085224 A CN 115085224A CN 202210995878 A CN202210995878 A CN 202210995878A CN 115085224 A CN115085224 A CN 115085224A
Authority
CN
China
Prior art keywords
new energy
power system
energy power
matrix
frequency response
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
CN202210995878.8A
Other languages
Chinese (zh)
Other versions
CN115085224B (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.)
Liyang Research Institute of Southeast University
Original Assignee
Liyang Research Institute of Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Liyang Research Institute of Southeast University filed Critical Liyang Research Institute of Southeast University
Priority to CN202210995878.8A priority Critical patent/CN115085224B/en
Publication of CN115085224A publication Critical patent/CN115085224A/en
Application granted granted Critical
Publication of CN115085224B publication Critical patent/CN115085224B/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
    • H02J3/241—The oscillation concerning frequency
    • 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/22—The renewable source being solar energy
    • H02J2300/24—The renewable source being solar energy of photovoltaic origin
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously

Landscapes

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

Abstract

The invention discloses a method for rapidly analyzing the frequency stability of a new energy power system, which comprises the steps of firstly, establishing a frequency response model of the new energy power system, wherein the frequency response model comprises a traditional generator set frequency response link, a wind power plant frequency response link and a photovoltaic power station frequency response link; secondly, performing model dimensionality reduction on the new energy power system frequency response model to obtain a low-order new energy power system frequency response model; and finally, analyzing the frequency stability by using a low-order frequency response model of the new energy power system, and realizing the rapid evaluation of the frequency stability of the new energy power system. The method fully considers the flexible control potential of the new energy multi-frequency modulation resource, can realize the rapid analysis of the stability and the dynamic characteristic of the frequency response high-dimensional model, and has profound significance for improving the operation stability of the new energy power system.

Description

New energy power system frequency stability rapid analysis method
Technical Field
The invention relates to the field of frequency stability analysis of power systems, in particular to a method for rapidly analyzing frequency stability of a new energy power system.
Background
The construction and development of a new energy power system is an important guarantee for realizing the double-carbon target, is an important measure for coping with global climate change and fulfilling international commitments, and has important significance for promoting the adjustment of energy structures in China, increasing energy supply and guaranteeing energy safety. The new energy power system actively participates in the adjustment of the frequency of the power system along with new energy such as wind power, photovoltaic power and the like, and the new energy power system comprises various types of typical elements such as wind power, photovoltaic power, thermal power generating units, hydroelectric generating units and the like. A prime system formed by the speed regulator of the hydroelectric generating set and the prime motor generates negative damping, and frequency oscillation is caused in a region with high water-electricity ratio, so that the system frequency stability is influenced. The new energy transient process with the converter as an interface has multi-time scale interweaving and discrete events, and the continuous process has characteristics of mixing, leading control strategy and the like, so that the high-order analysis of the model is difficult, and leading poles are difficult to distinguish.
In the existing system frequency response model considering that new energy participates in frequency modulation, the frequency modulation link of the new energy is generally simplified and equalized, so that the dimension of the model is reduced, and the calculation efficiency is improved. The new energy centralized and distributed construction trend leads to the existence of various frequency modulation resources with different frequency response dynamic characteristics in a new energy power system, and model errors brought by the traditional equivalent mode are not negligible.
Active frequency control of new energy is an important way for improving system frequency stability under high-power shortage, the existing research does not relate to a new energy power system frequency stability analysis method containing multiple frequency modulation resources, and the problem to be solved is that how to realize rapid analysis of new energy power system frequency stability is to be solved by considering the difference of dynamic characteristics of different frequency modulation resources.
Disclosure of Invention
In order to solve the above problems, the present invention provides a method for rapidly analyzing the frequency stability of a new energy power system, which can effectively determine the frequency stability of the new energy power system and analyze the dynamic change of the frequency.
In order to achieve the purpose, the invention is realized by the following technical scheme:
a method for rapidly analyzing frequency stability of a new energy power system comprises the following steps:
step 1: establishing a new energy power system frequency response model of multi-type frequency modulation resources comprising a traditional generator set frequency response link, a wind power plant frequency response link and a photovoltaic power station frequency response link;
step 2: performing model dimensionality reduction on the frequency response model of the new energy power system established in the step 1 to obtain a low-order frequency response model of the new energy power system;
and step 3: and analyzing the frequency stability according to the low-order new energy power system frequency response model.
The invention is further improved in that: the step 1 of establishing a frequency response link of the traditional generator set comprises the step of obtaining a traditional generator set type and corresponding frequency response parameters in a new energy power system, wherein the traditional generator set type at least comprises one of a thermal power generating unit and a hydroelectric generating unit.
The invention is further improved in that: in the step 1, the wind power plant frequency response link adopts comprehensive inertia control, and a transfer function meets the following formula:
Figure 663077DEST_PATH_IMAGE001
in the formula:
Figure 535087DEST_PATH_IMAGE002
the active variable quantity of the wind power plant;
Figure 824117DEST_PATH_IMAGE003
the angular frequency variation of the wind power plant;
Figure 707760DEST_PATH_IMAGE004
calculating a virtual inertia coefficient of the wind power plant;
Figure 320748DEST_PATH_IMAGE005
the droop coefficient of the wind power plant;
Figure 289841DEST_PATH_IMAGE006
the required time for transmitting the frequency modulation control command of the wind power plant power control system to the inverter is obtained;
Figure 382562DEST_PATH_IMAGE007
the execution time of the wind power plant inverter; s is a complex frequency domain;
the frequency response link of the photovoltaic power station adopts droop control, and a transfer function meets the following formula:
Figure 386290DEST_PATH_IMAGE008
in the formula:
Figure 155532DEST_PATH_IMAGE009
the active variable quantity of the photovoltaic power station is obtained;
Figure 611921DEST_PATH_IMAGE010
the angular frequency variation of the photovoltaic power station;
Figure 242754DEST_PATH_IMAGE011
adjusting difference coefficients of the photovoltaic power station;
Figure 835409DEST_PATH_IMAGE012
the required time for transmitting the frequency modulation control command of the photovoltaic power station power control system to the inverter is obtained;
Figure 41132DEST_PATH_IMAGE013
the required time for the photovoltaic power plant inverter to execute.
The invention is further improved in that: the specific steps of establishing the frequency response model of the new energy power system containing the multi-type frequency modulation resources in the step 1 are as follows: the frequency response links of all frequency modulation resources are connected in parallel to serve as negative feedback of a new energy power system model to form a new energy power system frequency response model, and the new energy power system frequency response model meets the following formula:
Figure 860183DEST_PATH_IMAGE014
in the formula:
Figure 419340DEST_PATH_IMAGE015
the angular frequency variation of the new energy power system is obtained;
Figure 256715DEST_PATH_IMAGE016
the active variable quantity of the new energy power system is obtained;His the inertia time constant of the generator;Dis the system damping coefficient.
The invention is further improved in that: the specific steps of the step 2 are as follows:
step 2.1, converting the frequency response model of the new energy power system established in the step 1 into a state space equation, and satisfying the following formula:
Figure 977547DEST_PATH_IMAGE017
Figure 283894DEST_PATH_IMAGE018
in the formula:
Figure 381163DEST_PATH_IMAGE019
is a roomnThe state vector of the dimension(s) is,
Figure 341553DEST_PATH_IMAGE020
is composed ofnThe differentiation of the dimensional state vector is performed,
Figure 233286DEST_PATH_IMAGE021
is composed ofqThe dimension controls the input vector of the input vector,
Figure 292509DEST_PATH_IMAGE022
is composed ofpThe vector is output in dimension and the vector is output,Ain order to be a matrix of the system,Bin order to input the matrix, the input matrix is,Cto be the output matrix, the output matrix is,D 1 is a direct transmission matrix;
step 2.2, calculating a controllable gram matrix of the new energy power systemPAnd considerable gram matrixQThe expression is:
Figure 193469DEST_PATH_IMAGE023
Figure 739857DEST_PATH_IMAGE024
in the formula:
Figure 68070DEST_PATH_IMAGE025
is a transpose of the system matrix,
Figure 83430DEST_PATH_IMAGE026
in order to be a transpose of the input matrix,
Figure 788081DEST_PATH_IMAGE027
is the transpose of the output matrix;
step 2.3, to the controllable gram matrixPAnd considerable gram matrixQCholesky decomposition is performed, and the expression is as follows:
Figure 454555DEST_PATH_IMAGE028
Figure 94615DEST_PATH_IMAGE029
in the formula:
Figure 721905DEST_PATH_IMAGE030
representing a controllable gram matrixPCholesky decomposition of (1);
Figure 89301DEST_PATH_IMAGE031
representing a controllable gram matrixPTranspose of Cholesky decomposition of (1);
Figure 485648DEST_PATH_IMAGE032
a Cholesky decomposition representing a observable gram matrix Q;
Figure 31030DEST_PATH_IMAGE033
a transpose representing the Cholesky decomposition of the observable gram matrix Q;
step 2.4, for
Figure 145616DEST_PATH_IMAGE034
Singular value decomposition is carried out:
Figure 319633DEST_PATH_IMAGE035
in the formula:Uis a unitary matrix;Σis a semi-positive definite diagonal matrix;
Figure 836065DEST_PATH_IMAGE036
are unitary matrices of different orders;
step 2.5, calculating a balance change matrix:
Figure 552348DEST_PATH_IMAGE037
Figure 154231DEST_PATH_IMAGE038
in the formula:
Figure 129009DEST_PATH_IMAGE039
is a balance change matrix after Cholesky decomposition of the controllable gram matrix P;
Figure 375314DEST_PATH_IMAGE040
is a considerable gram matrixQThe Cholesky decomposed equilibrium change matrix of (1);
step 2.6, carrying out balance transformation on the state space matrix:
Figure 387132DEST_PATH_IMAGE041
in the formula:
Figure 335365DEST_PATH_IMAGE042
a matrix of the system after the balance transformation is represented,
Figure 723621DEST_PATH_IMAGE043
representing the input matrix after the balance transformation,
Figure 90012DEST_PATH_IMAGE044
representing the output matrix after the balance transformation,
Figure 272731DEST_PATH_IMAGE045
the direct transfer matrix after the balance transformation is represented,
Figure 973840DEST_PATH_IMAGE046
is that
Figure 900208DEST_PATH_IMAGE047
Transposing;
and 2.7, truncating the state space matrix after the balance transformation according to the expected model order, calculating a state space equation of the new energy power system after the model dimensionality reduction, and converting the state space equation of the new energy power system after the model dimensionality reduction into a transfer function:
Figure 855525DEST_PATH_IMAGE048
in the formula:
Figure 474726DEST_PATH_IMAGE049
transmission of state space equations representing new energy power systemA function;
Figure 666060DEST_PATH_IMAGE050
，
Figure 271485DEST_PATH_IMAGE051
，
Figure 205943DEST_PATH_IMAGE052
respectively, the transfer function numerator 0 order,m-1order summ-2Coefficients of order;
Figure 120678DEST_PATH_IMAGE053
，
Figure 406166DEST_PATH_IMAGE054
respectively of the order 0 of the denominator of the transfer function,n-1the coefficients of the order of the first order,
Figure 549702DEST_PATH_IMAGE055
、
Figure 604246DEST_PATH_IMAGE056
respectively m-th order of the numerator and n-th order of the denominator.
The invention is further improved in that: the step 3 comprises the following specific steps: calculating a system pole of the low-order new energy power system frequency response model in the step 2, wherein the expression is as follows:
Figure 955462DEST_PATH_IMAGE057
judging the stability of the frequency of the new energy power system according to the positions of the system poles, wherein the new energy power system is stable if and only if all the system poles are on the left half plane of the complex plane; and selecting a conjugate complex root closest to the imaginary axis as a dominant pole, and analyzing the frequency dynamic characteristic through the position of the dominant pole on a complex plane and the damping ratio of the dominant pole, wherein the larger the imaginary part is, the smaller the damping ratio is, and the weaker the frequency stability of the new energy power system is.
The invention has the beneficial effects that: according to the method, a frequency response model of the new energy power system comprising elements such as a traditional generator set, a wind power plant, a photovoltaic power station and the like is established, the model is subjected to dimensionality reduction, and then the pole of a computing system is solved, so that frequency stability analysis is carried out, and the rapid evaluation of the frequency stability of the new energy power system is realized. The method fully considers the flexible control potential of the new energy multi-frequency modulation resource, can realize the rapid analysis of the stability and the dynamic characteristic of the frequency response high-dimensional model, and has profound significance for improving the operation stability of the new energy power system.
Drawings
FIG. 1 is a flow chart of the implementation of the method of the present invention;
FIG. 2 is a frequency response model of a new energy power system according to an embodiment of the present invention;
fig. 3 is a time domain simulation curve of the frequency response model of the new energy power system.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
As shown in fig. 1, the method of the present invention comprises the following steps:
step 1: establishing a new energy power system frequency response model, which comprises establishing a traditional generator set frequency response link, a wind power plant frequency response link and a photovoltaic power plant frequency response link, and establishing a new energy power system frequency response model containing multiple types of frequency modulation resources;
step 2: performing model dimensionality reduction on the frequency response model of the new energy power system established in the step 1 to obtain a low-order frequency response model of the new energy power system;
and step 3: and analyzing the frequency stability according to the low-order new energy power system frequency response model.
In the step 1, establishing a frequency response link of a traditional generator set, specifically comprising the following steps: the method comprises the steps of obtaining the type of a traditional power generation set in a new energy power system and corresponding frequency response parameters, wherein the type of the traditional power generation set comprises at least one of a thermal power generating unit and a hydroelectric generating unit.
When the traditional generating set is a thermal power generating set, the dynamic characteristics of a speed regulator and a prime motor are considered in the frequency response link of the thermal power generating set, and a transfer function meets the following formula:
Figure 462667DEST_PATH_IMAGE058
in the formula:
Figure 409894DEST_PATH_IMAGE059
the active variable quantity of the thermal power generating unit is obtained;
Figure 584523DEST_PATH_IMAGE060
the angular frequency variation of the thermal power generating unit is obtained;
Figure 575482DEST_PATH_IMAGE061
the percentage of the steady-state output power of the high-pressure cylinder to the total output power of the steam turbine is shown;
Figure 835562DEST_PATH_IMAGE062
is the intermediate reheat steam volume effect time constant;
Figure 586480DEST_PATH_IMAGE063
is the high pressure steam volume time constant;Rthe difference adjustment coefficient of the thermal power generating unit is obtained;
Figure 731461DEST_PATH_IMAGE064
is the time constant of the speed regulator of the thermal power generating unit; and S is a complex frequency domain.
When the traditional generator set is a hydroelectric generating set, the dynamic characteristics of a speed regulator and a prime motor are considered in the frequency response link of the hydroelectric generating set, and a transfer function meets the following formula:
Figure 768687DEST_PATH_IMAGE065
in the formula:
Figure 391429DEST_PATH_IMAGE066
the active variable quantity of the hydroelectric generating set;
Figure 805093DEST_PATH_IMAGE067
the angular frequency variation of the hydroelectric generating set;
Figure 813369DEST_PATH_IMAGE068
is the water flow inertia time constant;
Figure 287076DEST_PATH_IMAGE069
the difference adjustment coefficient of the hydroelectric generating set;
Figure 397115DEST_PATH_IMAGE070
，
Figure 348890DEST_PATH_IMAGE071
and
Figure 477252DEST_PATH_IMAGE072
the proportional coefficient, the integral coefficient and the differential coefficient of the speed regulator are respectively;
Figure 856281DEST_PATH_IMAGE073
is the time constant of the speed regulator of the hydroelectric generating set.
In the step 1, a wind power plant frequency response link and a photovoltaic power station frequency response link respectively adopt comprehensive inertia control and droop control, and a transfer function of the wind power plant frequency response link meets the following formula:
Figure 188036DEST_PATH_IMAGE074
in the formula:
Figure 943502DEST_PATH_IMAGE075
the active variable quantity of the wind turbine generator is obtained;
Figure 926371DEST_PATH_IMAGE076
the angular frequency variation of the wind turbine generator is obtained;
Figure 882825DEST_PATH_IMAGE077
the virtual inertia coefficient of the wind power plant is obtained;
Figure 826511DEST_PATH_IMAGE078
the droop coefficient of the wind power plant;
Figure 247652DEST_PATH_IMAGE079
the required time for transmitting the frequency modulation control command of the wind power plant power control system to the inverter is obtained;
Figure 960393DEST_PATH_IMAGE080
the time required for the wind farm inverter to execute.
The frequency response link transfer function of the photovoltaic power station meets the following formula:
Figure 822170DEST_PATH_IMAGE081
in the formula:
Figure 253151DEST_PATH_IMAGE082
the active variable quantity of the photovoltaic power station is obtained;
Figure 740633DEST_PATH_IMAGE083
the angular frequency variation of the photovoltaic power station;
Figure 42302DEST_PATH_IMAGE084
adjusting difference coefficients of the photovoltaic power station;
Figure 340559DEST_PATH_IMAGE085
the required time for transmitting the frequency modulation control command of the photovoltaic power station power control system to the inverter is obtained;
Figure 258837DEST_PATH_IMAGE086
the required time for the photovoltaic power plant inverter to execute.
And multiplying a traditional generator set frequency response link, a wind power plant frequency response link and a photovoltaic power station frequency response link by a proportional coefficient of each power generation unit accounting for all power output, and connecting the frequency response links of all frequency modulation resources in parallel to form a new energy power system frequency response model as negative feedback of the new energy power system model. The power system model satisfies the following formula:
Figure 18851DEST_PATH_IMAGE087
in the formula:
Figure 706184DEST_PATH_IMAGE088
the angular frequency variation of the new energy power system is obtained;
Figure 175343DEST_PATH_IMAGE089
the active variable quantity of the new energy power system is obtained;His the inertia time constant of the generator;Dis the system damping coefficient.
In the embodiment, a new energy power system comprising a photovoltaic power station, a wind power plant, a thermal power generating unit and a hydroelectric generating unit is selected as an object. And constructing a frequency response model of the new energy power system, as shown in fig. 2. The new energy power system is provided with 2 photovoltaic power stations, 2 wind power stations, 1 thermal power unit and 1 hydroelectric power unit. The parameters are set as follows:
photovoltaic power plant: 1 st photovoltaic power station difference adjustment coefficient
Figure 439971DEST_PATH_IMAGE090
=25, execution time of the 1 st photovoltaic power plant inverter
Figure 144622DEST_PATH_IMAGE091
=2, time for sending frequency modulation control command of power control system of 1 st photovoltaic power station to inverter
Figure 561828DEST_PATH_IMAGE092
=0.05, 2 nd photovoltaic power plant adjustment coefficient
Figure 326522DEST_PATH_IMAGE093
Execution time of 2 nd photovoltaic power station inverter =10
Figure 81375DEST_PATH_IMAGE094
Time required for sending frequency modulation control instruction of power control system of No. 2 photovoltaic power station to inverter
Figure 58559DEST_PATH_IMAGE095
=0.02；
Wind power plant: 1 st wind farm virtual inertia
Figure 595850DEST_PATH_IMAGE096
=30, droop coefficient of the 1 st wind farm
Figure 265866DEST_PATH_IMAGE097
=30, time for transmitting frequency modulation control command of 1 st wind power plant power control system to inverter
Figure 505086DEST_PATH_IMAGE098
=2, execution time of the 1 st wind farm inverter
Figure 285961DEST_PATH_IMAGE099
=0.03, virtual inertia of the 2 nd wind farm
Figure 943338DEST_PATH_IMAGE100
=15, droop coefficient of 2 nd wind farm
Figure 908889DEST_PATH_IMAGE101
=15, time for sending frequency modulation control command of 2 nd wind power plant power control system to inverter
Figure 510772DEST_PATH_IMAGE102
=5, execution time of 2 nd wind farm inverter
Figure 236282DEST_PATH_IMAGE103
=0.08；
Thermal power generating unit: percentage of steady state output power of high pressure cylinder to total output power of steam turbine
Figure 607221DEST_PATH_IMAGE104
=0.35, intermediate reheat steam volume effect time constant
Figure 743673DEST_PATH_IMAGE105
=12, high pressure steam volume time constant
Figure 708218DEST_PATH_IMAGE106
=0.3 difference adjustment coefficient of thermal power generating unitR=0.03, thermal power unit speed regulator time constant
Figure 96474DEST_PATH_IMAGE107
=0.2；
A hydroelectric generating set: water flow inertia time constant
Figure 449482DEST_PATH_IMAGE108
=2.5, difference adjustment coefficient of hydroelectric generating set
Figure 632202DEST_PATH_IMAGE109
=0.04, governor ratio
Figure 349622DEST_PATH_IMAGE110
=1, integral coefficient
Figure 275990DEST_PATH_IMAGE111
=2, differential coefficient
Figure 480575DEST_PATH_IMAGE112
=0.7, time constant of speed regulator of hydroelectric generating set
Figure 99775DEST_PATH_IMAGE113
=0.2；
Parameters of the power system: the equivalent inertia after considering the new energy is as follows: inertia time constant of generatorH=4.5, system damping coefficientD=0.1；
Power generation ratio coefficient:
Figure 38912DEST_PATH_IMAGE114
=0.025，
Figure 628026DEST_PATH_IMAGE115
=0.025，
Figure 828063DEST_PATH_IMAGE116
=0.025，
Figure 759110DEST_PATH_IMAGE117
=0.025，
Figure 779018DEST_PATH_IMAGE117
=0.025，
Figure 437402DEST_PATH_IMAGE117
=0.025，
Figure 491945DEST_PATH_IMAGE118
=0.180，
Figure 62735DEST_PATH_IMAGE119
=0.720。
as shown in fig. 2, the original new energy power system is a linear steady control system, the order of the model is 15, and the dimension of the model is reduced, and the specific steps are as follows:
1) converting the frequency response model of the new energy power system established in the step 1 into a state space equation, and satisfying the following formula:
Figure 835519DEST_PATH_IMAGE120
Figure 46662DEST_PATH_IMAGE018
in the formula:
Figure 955712DEST_PATH_IMAGE121
is composed ofnThe state vector of the dimension(s) is,
Figure 228562DEST_PATH_IMAGE122
to representnDifferentiation of the dimensional state vector;
Figure 347697DEST_PATH_IMAGE123
is composed ofqDimension control input vector;
Figure 223249DEST_PATH_IMAGE124
is composed ofpAnd outputting the vector dimension.AIn order to be a matrix of the system,Bin order to input the matrix, the input matrix is,Cin order to output the matrix, the input matrix,D 1 is a direct transfer matrix.
2) Controllable gram matrix for calculating new energy power systemPAnd considerable gram matrixQ：
Figure 127751DEST_PATH_IMAGE125
Figure 164977DEST_PATH_IMAGE126
In the formula:
Figure 36987DEST_PATH_IMAGE127
is a transpose of the system matrix,
Figure 450651DEST_PATH_IMAGE128
is a transpose of the input matrix,
Figure 475239DEST_PATH_IMAGE129
is the transpose of the output matrix;
3) for controllable gram matrixPAnd considerable gram matrixQPerforming Cholesky decomposition:
Figure 683366DEST_PATH_IMAGE130
Figure 777093DEST_PATH_IMAGE131
in the formula:
Figure 994447DEST_PATH_IMAGE132
representing a controllable gram matrixPCholesky decomposition of (1);
Figure 607963DEST_PATH_IMAGE133
representing a controllable gram matrixPTranspose of Cholesky decomposition of (1);
Figure 252571DEST_PATH_IMAGE134
a Cholesky decomposition representing a observable gram matrix Q;
Figure 836523DEST_PATH_IMAGE135
a transpose representing the Cholesky decomposition of the observable gram matrix Q;
4) to pair
Figure 732935DEST_PATH_IMAGE136
Singular value decomposition is carried out:
Figure 591170DEST_PATH_IMAGE137
in the formula:Uis a unitary matrix;Σis a semi-positive definite diagonal matrix; v is a unitary matrix of a different order,
Figure 531313DEST_PATH_IMAGE138
representing transposes of unitary matrices of different orders;
5) calculating an equilibrium change matrix:
Figure 209419DEST_PATH_IMAGE139
Figure 643942DEST_PATH_IMAGE140
in the formula:
Figure 356683DEST_PATH_IMAGE141
is the equilibrium change matrix after Cholesky decomposition of the observable gram matrix Q;
Figure 733307DEST_PATH_IMAGE142
is a balance change matrix after Cholesky decomposition of the controllable gram matrix P;
6) and (3) carrying out balance transformation on the state space matrix:
Figure 164288DEST_PATH_IMAGE041
in the formula:
Figure 136923DEST_PATH_IMAGE143
a matrix of the system after the balance transformation is represented,
Figure 704171DEST_PATH_IMAGE144
representing the input matrix after the balance transformation,
Figure 986117DEST_PATH_IMAGE145
representing the output matrix after the balance transformation,
Figure 638815DEST_PATH_IMAGE045
the direct transfer matrix after the balance transformation is represented,
Figure 680720DEST_PATH_IMAGE146
is that
Figure 495617DEST_PATH_IMAGE147
Transposing;
7) according to the expected model order, truncating the state space matrix after the balance transformation, calculating a state space equation of the system after the dimensionality reduction, and converting the state space equation of the system after the dimensionality reduction into a transfer function:
Figure 823830DEST_PATH_IMAGE148
in the formula:
Figure 839191DEST_PATH_IMAGE149
representing a transfer function of a state space equation of the new energy power system;
Figure 278262DEST_PATH_IMAGE150
，
Figure 210315DEST_PATH_IMAGE151
，
Figure 443850DEST_PATH_IMAGE152
respectively, the transfer function numerator 0 order,m-1order summ-2The coefficient of the order;
Figure 212086DEST_PATH_IMAGE153
，
Figure 454849DEST_PATH_IMAGE154
respectively of the order 0 of the denominator of the transfer function,n-1the coefficients of the order of the first order,
Figure 241408DEST_PATH_IMAGE155
、
Figure 52369DEST_PATH_IMAGE156
respectively m-th order of the numerator and n-th order of the denominator.
Step signals with the model order of 5 and 10 percent after dimensionality reduction are input into the system to obtain a system frequency dynamic response curve, as shown in fig. 3, the process that the system frequency changes and tends to be stable is completely consistent with the original new energy power system, and the model after dimensionality reduction can be well approximated to the dynamic process of the original new energy power system.
And 3, carrying out frequency stability analysis, solving the model after dimensionality reduction, and calculating a pole, wherein the expression is as follows:
Figure 901376DEST_PATH_IMAGE157
and judging the stability of the frequency of the new energy power system according to the positions of the poles of the system, wherein the new energy power system is stable if and only if all the poles of the system are on the left half plane of the complex plane. And selecting a conjugate complex root closest to the imaginary axis as a dominant pole, analyzing the frequency dynamic characteristic through the position of the dominant pole on a complex plane and the damping ratio of the dominant pole, wherein the larger the imaginary part is, the smaller the damping ratio is, and the weaker the frequency stability of the system is.
The dominant characteristic roots of the frequency response model of the new energy power system and the frequency response model of the original new energy power system after dimensionality reduction are calculated to be-0.1959 +/-0.5610iAnd-0.1961. + -. 0.5614iThe damping ratio of the dimension reduction model and the damping ratio of the full-order model are respectively 0.3296 and 0.3298, the leading characteristic roots before and after dimension reduction are basically consistent, but the original new energy power system contains 15 characteristic roots, leading poles are difficult to distinguish, and the frequency response model of the new energy power system after dimension reduction only has 5 poles, so that the leading poles are easy to distinguish. The order of the reduced dimension model is further adjusted, and the absolute error of the obtained reduced order system is shown in table 1:
Figure 72464DEST_PATH_IMAGE158
along with the increase of the order of the model, the absolute error of the reduced-order model and the frequency response model of the original new energy power system is obviously reduced.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (6)

1. A method for rapidly analyzing frequency stability of a new energy power system is characterized by comprising the following steps: the method comprises the following steps:
step 1: establishing a new energy power system frequency response model comprising a traditional generator set frequency response link, a wind power plant frequency response link and a photovoltaic power station frequency response link;
step 2: performing model dimensionality reduction on the new energy power system frequency response model established in the step 1 to obtain a low-order new energy power system frequency response model;
and step 3: and analyzing the frequency stability according to the low-order new energy power system frequency response model.
2. The method for rapidly analyzing the frequency stability of the new energy power system according to claim 1, wherein the method comprises the following steps: in the step 1, establishing a frequency response link of a traditional generator set comprises acquiring a traditional generator set type and corresponding frequency response parameters in a new energy power system, wherein the traditional generator set type at least comprises one of a thermal power generating unit and a hydroelectric generating unit.
3. The method for rapidly analyzing the frequency stability of the new energy power system according to claim 1, wherein the method comprises the following steps: in the step 1, the frequency response link of the wind power plant adopts comprehensive inertia control, and a transfer function meets the following formula:
Figure 472087DEST_PATH_IMAGE001
in the formula:
Figure 506908DEST_PATH_IMAGE002
the active variable quantity of the wind power plant;
Figure 963297DEST_PATH_IMAGE003
the angular frequency variation of the wind power plant;
Figure 594130DEST_PATH_IMAGE004
the virtual inertia coefficient of the wind power plant is obtained;
Figure 452364DEST_PATH_IMAGE005
the droop coefficient of the wind power plant;
Figure 383718DEST_PATH_IMAGE006
the required time for transmitting the frequency modulation control command of the wind power plant power control system to the inverter is obtained;
Figure 327403DEST_PATH_IMAGE007
the required time for the wind farm inverter to execute; s is a complex frequency domain;
the frequency response link of the photovoltaic power station adopts droop control, and a transfer function meets the following formula:
Figure 761927DEST_PATH_IMAGE008
in the formula:
Figure 474668DEST_PATH_IMAGE009
the active variable quantity of the photovoltaic power station is obtained;
Figure 585712DEST_PATH_IMAGE010
the angular frequency variation of the photovoltaic power station;
Figure 751114DEST_PATH_IMAGE011
adjusting difference coefficients of the photovoltaic power station;
Figure 989329DEST_PATH_IMAGE012
the required time for transmitting the frequency modulation control command of the photovoltaic power station power control system to the inverter is obtained;
Figure 946789DEST_PATH_IMAGE013
the required time for the photovoltaic power plant inverter to execute.
4. The method for rapidly analyzing the frequency stability of the new energy power system according to claim 1, wherein the method comprises the following steps: the specific steps of establishing a new energy power system frequency response model comprising a traditional generator set frequency response link, a wind power plant frequency response link and a photovoltaic power plant frequency response link in the step 1 are as follows: the frequency response links of all frequency modulation resources are connected in parallel to serve as negative feedback of a new energy power system model to form a new energy power system frequency response model, and the new energy power system frequency response model meets the following formula:
Figure 369681DEST_PATH_IMAGE014
in the formula:
Figure 897745DEST_PATH_IMAGE015
the angular frequency variation of the new energy power system is obtained;
Figure 798705DEST_PATH_IMAGE016
the active variable quantity of the new energy power system is obtained;His the inertia time constant of the generator;Dthe damping coefficient of the system is; and S is a complex frequency domain.
5. The method for rapidly analyzing the frequency stability of the new energy power system according to claim 1, wherein the method comprises the following steps: the specific steps of the step 2 are as follows:
step 2.1, converting the frequency response model of the new energy power system established in the step 1 into a state space equation, and satisfying the following formula:
Figure 610672DEST_PATH_IMAGE017
Figure 938885DEST_PATH_IMAGE018
in the formula:
Figure 219825DEST_PATH_IMAGE019
is composed ofnThe state vector of the dimension(s) is,
Figure 786460DEST_PATH_IMAGE020
is composed ofnThe differential of the dimensional state vector is taken,
Figure 328300DEST_PATH_IMAGE021
is composed ofqThe dimensions of the input vector are controlled in such a way that,
Figure 702780DEST_PATH_IMAGE022
is composed ofpThe vector is output in dimension and the vector is output,Ain order to be a matrix of the system,Bin order to input the matrix, the input matrix is,Cin order to output the matrix, the input matrix,D 1 is a direct transmission matrix;
step 2.2, calculating a controllable gram matrix of the new energy power systemPAnd considerable gram matrixQThe expression is:
Figure 330071DEST_PATH_IMAGE023
Figure 697467DEST_PATH_IMAGE024
in the formula:
Figure 93813DEST_PATH_IMAGE025
is a transpose of the system matrix,
Figure 904775DEST_PATH_IMAGE026
in order to be a transpose of the input matrix,
Figure 19361DEST_PATH_IMAGE027
is the transpose of the output matrix;
step 2.3, to the controllable gram matrixPAnd considerable gram matrixQCholesky decomposition is performed, and the expression is as follows:
Figure 190448DEST_PATH_IMAGE028
Figure 441301DEST_PATH_IMAGE029
in the formula:
Figure 423164DEST_PATH_IMAGE030
representing a controllable gram matrixPCholesky decomposition of (1);
Figure 149680DEST_PATH_IMAGE031
representing a controllable gram matrixPTranspose of Cholesky decomposition of (1);
Figure 734245DEST_PATH_IMAGE032
a Cholesky decomposition representing a observable gram matrix Q;
Figure 980550DEST_PATH_IMAGE033
a transpose representing the Cholesky decomposition of the observable gram matrix Q;
step 2.4, for
Figure 992368DEST_PATH_IMAGE034
Singular value decomposition is carried out:
Figure 474690DEST_PATH_IMAGE035
in the formula:Uis a unitary matrix;Σis a semi-positive definite diagonal matrix;
Figure 597366DEST_PATH_IMAGE036
are unitary matrices of different orders;
step 2.5, calculating a balance change matrix:
Figure 698178DEST_PATH_IMAGE037
Figure 146476DEST_PATH_IMAGE038
in the formula:
Figure 582006DEST_PATH_IMAGE039
is the equilibrium change matrix after Cholesky decomposition of the observable gram matrix Q;
Figure 508373DEST_PATH_IMAGE040
is the equilibrium change matrix after Cholesky decomposition of the controllable gram matrix P;
step 2.6, carrying out balance transformation on the state space matrix:
Figure 994850DEST_PATH_IMAGE041
in the formula:
Figure 473104DEST_PATH_IMAGE042
a matrix of the system after the balance transformation is represented,
Figure 271296DEST_PATH_IMAGE043
representing the input matrix after the balance transformation,
Figure 876721DEST_PATH_IMAGE044
representing the output matrix after the balance transformation,
Figure 76758DEST_PATH_IMAGE045
the direct transfer matrix after the balance transformation is represented,
Figure 725914DEST_PATH_IMAGE046
is that
Figure 276981DEST_PATH_IMAGE047
Transposing;
and 2.7, truncating the state space matrix after the balance transformation according to the expected model order, calculating a state space equation of the new energy power system after the model dimensionality reduction, and converting the state space equation of the new energy power system after the model dimensionality reduction into a transfer function:
Figure 420518DEST_PATH_IMAGE048
in the formula:
Figure 475061DEST_PATH_IMAGE049
representing a transfer function of a state space equation of the new energy power system;
Figure 94786DEST_PATH_IMAGE050
，
Figure 742936DEST_PATH_IMAGE051
，
Figure 549218DEST_PATH_IMAGE052
respectively, the transfer function numerator 0 order,m-1order summ-2Coefficients of order;
Figure 114061DEST_PATH_IMAGE053
，
Figure 855752DEST_PATH_IMAGE054
respectively of the order 0 of the denominator of the transfer function,n-1the coefficient of the order of the first order,
Figure 115832DEST_PATH_IMAGE055
、
Figure 850438DEST_PATH_IMAGE056
respectively m-th order of the numerator and n-th order of the denominator.
6. The method for rapidly analyzing the frequency stability of the new energy power system according to claim 5, wherein the method comprises the following steps: the specific steps of the step 3 are as follows: calculating a system pole of the low-order new energy power system frequency response model in the step 2, wherein the expression is as follows:
Figure 879574DEST_PATH_IMAGE057
judging the stability of the frequency of the new energy power system according to the positions of the system poles, wherein the new energy power system is stable if and only if all the system poles are on the left half plane of the complex plane; and selecting the conjugate complex root closest to the imaginary axis as a dominant pole, analyzing the frequency dynamic characteristic according to the position of the dominant pole on the complex plane and the damping ratio of the dominant pole, wherein the larger the imaginary part is, the smaller the damping ratio is, and the weaker the frequency stability of the new energy power system is.
CN202210995878.8A 2022-08-19 2022-08-19 New energy power system frequency stability rapid analysis method Active CN115085224B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210995878.8A CN115085224B (en) 2022-08-19 2022-08-19 New energy power system frequency stability rapid analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210995878.8A CN115085224B (en) 2022-08-19 2022-08-19 New energy power system frequency stability rapid analysis method

Publications (2)

Publication Number Publication Date
CN115085224A true CN115085224A (en) 2022-09-20
CN115085224B CN115085224B (en) 2022-11-01

Family

ID=83244242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210995878.8A Active CN115085224B (en) 2022-08-19 2022-08-19 New energy power system frequency stability rapid analysis method

Country Status (1)

Country Link
CN (1) CN115085224B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115764979A (en) * 2022-10-12 2023-03-07 东南大学溧阳研究院 Distributed photovoltaic system stability quantitative evaluation method considering communication delay
CN117311138A (en) * 2023-11-30 2023-12-29 华中科技大学 Method and system for calculating stability margin domain of control parameter of water turbine adjusting system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120083935A1 (en) * 2010-10-04 2012-04-05 Wells Charles H Decoupling controller for power systems
CN112636368A (en) * 2020-12-10 2021-04-09 南京工程学院 Automatic power generation control method for multi-source multi-region interconnected power system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120083935A1 (en) * 2010-10-04 2012-04-05 Wells Charles H Decoupling controller for power systems
CN112636368A (en) * 2020-12-10 2021-04-09 南京工程学院 Automatic power generation control method for multi-source multi-region interconnected power system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115764979A (en) * 2022-10-12 2023-03-07 东南大学溧阳研究院 Distributed photovoltaic system stability quantitative evaluation method considering communication delay
CN115764979B (en) * 2022-10-12 2023-12-08 东南大学溧阳研究院 Distributed photovoltaic system stability quantitative evaluation method considering communication delay
CN117311138A (en) * 2023-11-30 2023-12-29 华中科技大学 Method and system for calculating stability margin domain of control parameter of water turbine adjusting system
CN117311138B (en) * 2023-11-30 2024-02-23 华中科技大学 Method and system for calculating stability margin domain of control parameter of water turbine adjusting system

Also Published As

Publication number Publication date
CN115085224B (en) 2022-11-01

Similar Documents

Publication Publication Date Title
CN115085224B (en) New energy power system frequency stability rapid analysis method
Guo et al. Hydropower system operation stability considering the coupling effect of water potential energy in surge tank and power grid
CN110429648B (en) Small interference stability margin probability evaluation method considering wind speed random fluctuation
CN110750882A (en) Wind power ratio limit value analytical calculation method considering frequency constraint
Liu et al. Damping characteristics analysis of hydropower units under full operating conditions and control parameters: Accurate quantitative evaluation based on refined models
CN107947228B (en) Stochastic stability analysis method for power system containing wind power based on Markov theory
CN111884259A (en) Site-level wind power generating set self-adaptive equivalence method considering system small interference stability characteristics
Morovati et al. Robust output feedback control design for inertia emulation by wind turbine generators
Zou et al. Design of intelligent nonlinear robust controller for hydro-turbine governing system based on state-dynamic-measurement hybrid feedback linearization method
Zou et al. Optimized Robust Controller Design Based on CPSOGSA Optimization Algorithm and H 2/H∞ Weights Distribution Method for Load Frequency Control of Micro-Grid
Gao et al. A fast high-precision model of the doubly-fed pumped storage unit
CN116436042B (en) Wind-water-fire system stability analysis method considering wind turbine frequency modulation dead zone
CN103956767B (en) A kind of wind farm grid-connected method for analyzing stability considering wake effect
CN110048428B (en) Power system probability load flow calculation method based on probability conservation principle
CN114880863B (en) Self-adaptive frequency division order reduction method for distributed renewable energy cluster impedance aggregation model
CN106952180B (en) Method for establishing double-fed distributed wind power system low-order frequency response model
Fang et al. Modeling and simulation of hydraulic transients for hydropower plants
CN113162063B (en) Design method of multi-direct-current coordination controller for inhibiting ultralow frequency oscillation
Li et al. Small Signal stability analysis and optimize control of large-scale wind power collection system
CN114465280A (en) Dynamic equivalent modeling method for new energy grid-connected system
CN114298478A (en) Small disturbance stability identification method and system for wind power grid-connected system
Xia et al. Frequency regulation strategy for AC–DC system during black start
CN113595065A (en) Method for inhibiting ultralow frequency oscillation of water-light complementary system
Xu et al. Stability of hydropower units under full operating conditions considering nonlinear coupling of turbine characteristics
CN111969624A (en) Damping control method and system of wind power grid-connected system containing virtual synchronous generator

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