CN111769575B - Fan parameter optimization oscillation suppression system and method based on modal stability domain - Google Patents

Fan parameter optimization oscillation suppression system and method based on modal stability domain Download PDF

Info

Publication number
CN111769575B
CN111769575B CN202010682091.7A CN202010682091A CN111769575B CN 111769575 B CN111769575 B CN 111769575B CN 202010682091 A CN202010682091 A CN 202010682091A CN 111769575 B CN111769575 B CN 111769575B
Authority
CN
China
Prior art keywords
power grid
wind power
connected system
modal
hyperplane
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.)
Active
Application number
CN202010682091.7A
Other languages
Chinese (zh)
Other versions
CN111769575A (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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202010682091.7A priority Critical patent/CN111769575B/en
Publication of CN111769575A publication Critical patent/CN111769575A/en
Application granted granted Critical
Publication of CN111769575B publication Critical patent/CN111769575B/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
    • 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/28—The renewable source being wind energy
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The invention relates to a fan parameter optimization oscillation suppression system and method based on a modal stability domain, belongs to the technical field of wind power generation, and solves the problems that in the prior art, a wind power grid-connected system has poor oscillation suppression effect and is not beneficial to stable operation of a power grid. The system comprises: the data acquisition module is used for acquiring related data of the wind power grid-connected system; the system aperiodic dynamic energy acquisition module is used for acquiring the total aperiodic dynamic energy of the wind power grid-connected system; the system modal stability domain construction module is used for constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the total non-periodic dynamic energy of the wind power grid-connected system; the modal stability domain boundary fitting module is used for performing hyperplane fitting on the modal stability domain to obtain a fitted hyperplane; and the wind turbine generator control parameter optimization module is used for optimizing the control parameters of the double-fed wind turbine generator. The system can quickly and effectively inhibit oscillation in the wind power grid-connected system and ensure stable operation of the system.

Description

Fan parameter optimization oscillation suppression system and method based on modal stability domain
Technical Field
The invention relates to the technical field of wind power generation, in particular to a system and a method for inhibiting fan parameter optimization oscillation based on a modal stability domain.
Background
The new energy strategy in China sets the key point for the vigorous development of wind power generation. The fluctuation of the wind power generation power, the phase-locked loop, the virtual inertia and the like are additional control links for improving the performance of the wind power generation system, so that the operation mode, the dynamic characteristic and the like of the power system are essentially changed, and particularly after the system is disturbed, the oscillation problem is aggravated by the behavior of a fan different from that of a traditional synchronous generator. Therefore, it is necessary to intensively study a method for suppressing oscillation caused by disturbance of a power system with high wind power permeability.
At present, a low-frequency oscillation mode research method for a wind power grid-connected system is mainly an analysis method based on real-time measurement information, and the method can only process single variables such as a generator power angle, generator active power or tie line flowing power in a certain local part and a short time, and cannot perform overall analysis on the system from a global perspective so as to effectively inhibit oscillation of the wind power grid-connected system. In addition, the existing fan oscillation suppression strategy research in the traditional damping control field such as a Power System Stabilizer (PSS) and a Flexible Alternating Current Transmission System (FACTS) is referred to, and the influence of the fan parameters on the oscillation mode is ignored.
The prior art has at least the following defects that firstly, the oscillation of a wind power grid-connected system can be analyzed and inhibited only by aiming at a single factor, the inhibition efficiency is low, and the effect is poor; secondly, influence of parameters of the fan on the oscillation mode is neglected, and when the parameters are unreasonably set, the wind generating set is possibly in a low damping level or a negative damping state, and stable operation of a power grid is threatened.
Disclosure of Invention
In view of the analysis, the invention aims to provide a system and a method for suppressing fan parameter optimization oscillation based on a modal stability domain, so as to solve the problems of poor effect, low efficiency and easy threat to the stable operation of a power grid in the existing wind power grid-connected system oscillation suppression.
On one hand, the invention provides a fan parameter optimization oscillation suppression system based on a modal stability domain, which is used for suppressing oscillation of a wind power grid-connected system by adjusting control parameters of a doubly-fed wind generating set, and comprises the following components:
the data acquisition module is used for acquiring the operation data, the oscillation mode data and the control parameters of the double-fed wind generating set of the wind power grid-connected system;
the system aperiodic dynamic energy acquisition module is used for acquiring the total aperiodic dynamic energy of the wind power grid-connected system based on the operation data and the oscillation modal data of the wind power grid-connected system and the control parameters of the doubly-fed wind generating set;
the system modal stability domain construction module is used for constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the obtained total aperiodic dynamic energy of the wind power grid-connected system;
the modal stability domain boundary fitting module is used for performing hyperplane fitting on a modal stability domain of the wind power grid-connected system to obtain a fitting hyperplane;
and the wind turbine generator set control parameter optimization module is used for optimizing the control parameters of the doubly-fed wind turbine generator set based on the fitting obtained hyperplane, by taking the maximum modal stability distance of the wind power grid-connected system as a target function and the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system as constraint conditions.
Further, the system aperiodic dynamic energy obtaining module obtains the total aperiodic dynamic energy of the wind power grid-connected system by using a system aperiodic dynamic energy model based on the operation data and the oscillation modal data of the wind power grid-connected system and the control parameters of the doubly-fed wind generating set;
the non-periodic dynamic energy model of the system is as follows:
Figure BDA0002586229260000021
Figure BDA0002586229260000022
Figure BDA0002586229260000023
wherein, UsIs the grid-connected point bus voltage amplitude, U, of the wind power grid-connected systemDFIGFor doubly-fed wind generator terminal voltage amplitude, xDFIGIs the equivalent impedance, theta, of a doubly-fed wind generatorpll0The doubly-fed wind generator phase locking angle theta corresponding to the initial disturbance moments0The initial phase angle of the grid-connected point bus voltage at the initial disturbance moment, a is the fluctuation amplitude of the grid-connected point bus voltage after the phase angle disturbance, omega is the fluctuation angular velocity of the grid-connected point bus voltage after the phase angle disturbance, t is time, IsIs the current amplitude of the grid-connected point bus,
Figure BDA0002586229260000031
is the initial phase angle of the phase-locked angle, KωFor controlling parameters for virtual inertia, KpThe proportional parameter of a phase-locked loop in the wind power grid-connected system is obtained; kIIs an integration parameter of the phase locked loop.
Further, the system modal stability domain building module builds a modal stability domain of the wind power grid-connected system by the following method:
constructing a small disturbance stability domain in a node active power injection space of the wind power grid-connected system as the bottom surface of a modal stability domain;
applying disturbance with different energies to an operation point of any wind power grid-connected system in the bottom surface, carrying out oscillation mode analysis on the operation point to obtain modal frequency and amplitude information of the operation point under each oscillation mode, judging whether the operation point is a critical stable point under a corresponding oscillation mode according to modal stability criteria based on the modal frequency and amplitude information of the operation point under each oscillation mode, and if so, acquiring a unit-time total non-periodic dynamic energy value of the wind power grid-connected system at the operation point as a critical point on a modal stable domain boundary under the corresponding oscillation mode;
traversing all the operating points of the wind power grid-connected system in the bottom surface to obtain a critical point on the boundary of the modal stability domain under each oscillation mode, and fitting the critical point under each oscillation mode to obtain the modal stability domain corresponding to each mode.
Further, the modal stability domain boundary fitting module performs hyperplane fitting on the modal stability domain of the wind power grid-connected system by using the following method:
fitting and obtaining a hyperplane of a modal stability domain of the wind power grid-connected system under each oscillation mode according to the following formula:
Figure BDA0002586229260000032
estimating the fitting accuracy of the obtained hyperplane according to the following formula:
Figure BDA0002586229260000041
wherein, PiActive power, alpha, injected into the ith node of the wind power grid-connected systemiIs a hyperplane parameter, i is 1,2, k is the number of nodes of the wind power grid-connected system, and σ isESRFitting accuracy of the hyperplane;
and when the fitting precision of the hyperplane is greater than a preset threshold value, changing the hyperplane parameters to perform refitting until the obtained fitting precision of the hyperplane is less than the preset threshold value.
Further, the objective function in the wind turbine generator control parameter optimization module is as follows:
Figure BDA0002586229260000042
wherein m represents the dominant oscillation mode of the wind power grid-connected system, PimActive power alpha injected into the ith node of the wind power grid-connected system under m dominant oscillation modesimAnd (3) corresponding hyperplane parameters under the m dominant oscillation mode, wherein i is 1,2, k is the number of nodes of the wind power grid-connected system.
Further, the stability margin constraint condition is:
Sm(P1m,P2m,...,Pkm)≥H0，
the constraint conditions of the tidal current static safety and stability condition are as follows:
Figure BDA0002586229260000043
the frequency range constraint conditions of the wind power grid-connected system are as follows:
Figure BDA0002586229260000044
wherein H0M is a dominant oscillation mode for a preset mode stable distance, J (K) is the virtual inertia of the wind power grid-connected system, JsetTo preset virtual inertia, PimActive power injected into the ith node of the wind power grid-connected system under the m-dominant oscillation mode, wherein i is 1,2P_PLLFor proportional control parameters of the phase-locked loop, KI_PLLFor integral control parameters of phase-locked loops, KωFor virtual inertia control parameters, Ω1Is KP_PLLValue range of (1), omega2Is KI_PLLValue range of (1), omega3Is KωThe value range of (a).
On the other hand, the invention provides a fan parameter optimization oscillation suppression method based on a modal stability domain, which is used for suppressing oscillation of a wind power grid-connected system by adjusting control parameters of a doubly-fed wind generating set, and comprises the following steps:
acquiring operation data, oscillation mode data and control parameters of a double-fed wind generating set of a wind power grid-connected system;
acquiring total aperiodic dynamic energy of the wind power grid-connected system by utilizing a system aperiodic dynamic energy model based on the operation data and the oscillation mode data of the wind power grid-connected system and the control parameters of the double-fed wind generating set;
constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the obtained total aperiodic dynamic energy of the wind power grid-connected system;
performing hyperplane fitting on a modal stability region of the wind power grid-connected system to obtain a fitting hyperplane;
and optimizing control parameters of the double-fed wind generating set by taking the maximum modal stability distance of the wind power grid-connected system as a target function and the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system as constraint conditions on the basis of the hyperplane obtained by fitting.
Further, a modal stability domain of the wind power grid-connected system is constructed in the following way:
constructing a small disturbance stability domain in a node active power injection space of the wind power grid-connected system as the bottom surface of a modal stability domain;
applying disturbance with different energies to an operation point of any wind power grid-connected system in the bottom surface, carrying out oscillation mode analysis on the operation point to obtain modal frequency and amplitude information of the operation point under each oscillation mode, judging whether the operation point is a critical stable point under a corresponding oscillation mode according to modal stability criteria based on the modal frequency and amplitude information of the operation point under each oscillation mode, and if so, acquiring a unit-time total non-periodic dynamic energy value of the wind power grid-connected system at the operation point as a critical point on a modal stable domain boundary under the corresponding oscillation mode;
traversing all the operating points of the wind power grid-connected system in the bottom surface to obtain a critical point on the boundary of the modal stability domain under each oscillation mode, and fitting the critical point under each oscillation mode to obtain the modal stability domain corresponding to each mode.
Further, hyperplane fitting is carried out on the modal stability domain of the wind power grid-connected system in the following mode:
fitting and obtaining a hyperplane of a modal stability domain of the wind power grid-connected system under each oscillation mode according to the following formula:
Figure BDA0002586229260000061
estimating the fitting accuracy of the obtained hyperplane according to the following formula:
Figure BDA0002586229260000062
wherein, PiActive power, alpha, injected into the ith node of the wind power grid-connected systemiIs a hyperplane parameter, i is 1,2, k is the number of nodes of the wind power grid-connected system, and σ isESRFitting accuracy of the hyperplane;
and when the fitting precision of the hyperplane is greater than a preset threshold value, changing the hyperplane parameters to perform refitting until the obtained fitting precision of the hyperplane is less than the preset threshold value.
Further, the objective function is:
Figure BDA0002586229260000063
the stability margin constraint conditions are as follows:
Sm(P1m,P2m,...,Pkm)≥H0，
the constraint conditions of the tidal current static safety and stability condition are as follows:
Figure BDA0002586229260000064
the frequency range constraint conditions of the wind power grid-connected system are as follows:
Figure BDA0002586229260000071
wherein m represents the dominant oscillation mode of the wind power grid-connected system, PimActive power alpha injected into the ith node of the wind power grid-connected system under m dominant oscillation modesimCorresponding hyperplane parameters under m dominant oscillation modes, i is 1,2, k is the number of nodes of the wind power grid-connected system, and H is the number of the nodes of the wind power grid-connected system0A preset modal stability distance, J (K) is a virtual inertia of the wind power grid-connected system, JsetTo preset virtual inertia, KP_PLLFor proportional control parameters of the phase-locked loop, KI_PLLFor integral control parameters of phase-locked loops, KωFor virtual inertia control parameters, Ω1Is KP_PLLValue range of (1), omega2Is KI_PLLValue range of (1), omega3Is KωThe value range of (a).
Compared with the prior art, the invention can realize at least one of the following beneficial effects:
1. the fan parameter optimization oscillation suppression system provided by the invention optimizes the control parameters of the double-fed wind generating set by taking the maximum modal stability distance of the wind power grid-connected system as a target function and the frequency range, stability margin and tidal current static safety and stability conditions of the wind power grid-connected system as constraint conditions, takes the whole factors influencing the oscillation of the wind power grid-connected system into consideration, enables the wind power grid-connected system to quickly recover a stable operation state after disturbance disappears through parameter optimization, and avoids the problems of poor oscillation suppression effect and low suppression efficiency caused by only considering the influence of a single factor on the oscillation of the wind power grid-connected system.
2. The fan parameter optimization oscillation suppression system provided by the invention takes the influence of the fan parameters on the oscillation of the wind power grid-connected system as the constraint condition for optimizing the control parameters of the doubly-fed wind generating set, and avoids the problem that the wind generating set is possibly in a low damping level or a negative damping state to threaten the stable operation of a power grid when the self-set parameters are unreasonable.
In the invention, the technical schemes can be combined with each other to realize more preferable combination schemes. Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and drawings.
Drawings
The drawings are only for purposes of illustrating particular embodiments and are not to be construed as limiting the invention, wherein like reference numerals are used to designate like parts throughout.
FIG. 1 is a schematic diagram of a fan parameter optimization oscillation suppression system based on a modal stability domain according to an embodiment of the present invention;
FIG. 2 is a flowchart of a method for damping oscillation of a wind turbine parameter optimization based on a modal stability domain according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a structure of an IEEE 4 machine 11 node wind power grid-connected system according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a power angle response variation curve of each unit of the wind power grid-connected system before the control parameters of the doubly-fed wind generating set are optimized under a stable working condition according to the embodiment of the invention;
FIG. 5 is a schematic diagram of a power angle response variation curve of each unit of the wind power grid-connected system after control parameters of the doubly-fed wind generating set are optimized under a stable working condition according to the embodiment of the invention;
FIG. 6 is a schematic diagram of a power angle response variation curve of each unit of the wind power grid-connected system before optimization of control parameters of the doubly-fed wind generating set under the instability condition in the embodiment of the present invention;
fig. 7 is a schematic diagram of a power angle response change curve of each unit of the wind power grid-connected system after control parameters of the doubly-fed wind generating set are optimized under the instability condition in the embodiment of the invention.
Detailed Description
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate preferred embodiments of the invention and together with the description, serve to explain the principles of the invention and not to limit the scope of the invention.
System embodiment
The embodiment of the invention discloses a fan parameter optimization oscillation suppression system based on a modal stability domain, and is shown in fig. 1. The system is used for carrying out oscillation suppression on the wind power grid-connected system by adjusting the control parameters of the double-fed wind generating set, and the double-fed wind generating set is connected with a power grid to form the wind power grid-connected system.
This fan parameter optimization oscillation suppression system includes:
and the data acquisition module is used for acquiring the operation data, the oscillation mode data and the control parameters of the double-fed wind generating set of the wind power grid-connected system.
And the system aperiodic dynamic energy acquisition module is used for acquiring the total aperiodic dynamic energy of the wind power grid-connected system based on the operation data and the oscillation mode data of the wind power grid-connected system and the control parameters of the double-fed wind generating set.
And the system modal stability domain construction module is used for constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the obtained total non-periodic dynamic energy of the wind power grid-connected system.
And the modal stability domain boundary fitting module is used for performing hyperplane fitting on a modal stability domain of the wind power grid-connected system to obtain a fitting hyperplane.
And the wind turbine generator set control parameter optimization module is used for optimizing the control parameters of the doubly-fed wind turbine generator set based on the fitting obtained hyperplane, by taking the maximum modal stability distance of the wind power grid-connected system as an objective function and taking the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system as constraint conditions, so that the wind power grid-connected system can quickly recover a stable operation state after disturbance disappears.
Preferably, the system aperiodic dynamic energy obtaining module obtains the total aperiodic dynamic energy of the wind power grid-connected system by using a system aperiodic dynamic energy model based on the operation data and the oscillation mode data of the wind power grid-connected system and the control parameters of the doubly-fed wind generating set.
The non-periodic dynamic energy model of the system is as follows:
Figure BDA0002586229260000091
Figure BDA0002586229260000092
Figure BDA0002586229260000093
wherein, UsIs the grid-connected point bus voltage amplitude, U, of the wind power grid-connected systemDFIGFor doubly-fed wind generator terminal voltage amplitude, xDFIGIs the equivalent impedance, theta, of a doubly-fed wind generatorpll0The doubly-fed wind generator phase locking angle theta corresponding to the initial disturbance moments0The initial phase angle of the grid-connected point bus voltage at the initial disturbance moment, a is the fluctuation amplitude of the grid-connected point bus voltage after the phase angle disturbance, omega is the fluctuation angular velocity of the grid-connected point bus voltage after the phase angle disturbance, t is time, IsIs the current amplitude of the grid-connected point bus,
Figure BDA0002586229260000101
is the initial phase angle of the phase-locked angle, KωFor controlling parameters for virtual inertia, KpThe proportional parameter of a phase-locked loop in the wind power grid-connected system is obtained; kIIs an integration parameter of the phase locked loop.
The total aperiodic dynamic energy of the wind power grid-connected system obtained based on the system aperiodic dynamic energy model is as follows:
Figure 1
Figure BDA0002586229260000103
Figure BDA0002586229260000104
Figure BDA0002586229260000105
wherein n is the number of fans, m is the number of synchronous generators, UDFIGiIs the terminal voltage amplitude, U, of the ith fansFor the grid-connected point bus voltage amplitude, xDFIGiIs the equivalent impedance of the ith fan, thetapll0For disturbing the phase angle, theta, of the fan at the initial moments0The initial phase angle of the bus voltage of the grid-connected point at the initial disturbance moment, a is the fluctuation amplitude after disturbance of the phase angle, omega is the fluctuation angular velocity after disturbance, t is time, IsFor the magnitude of the bus current of the grid-connected point,
Figure BDA0002586229260000106
is the initial phase angle of the phase-locked angle variation value, KωFor controlling parameters for virtual inertia, KpiIs the proportional parameter, K, of the phase-locked loop of the ith fan in the wind power grid-connected systemIiIs an integral parameter, U, of a phase-locked loop of the ith fanGjIs the voltage amplitude of the terminal of the jth synchronous generator, xGjIs the equivalent impedance of the jth synchronous generator, delta0To disturb the power angle of the synchronous generator at the initial moment, DjDamping coefficient, T, of the jth synchronous generatorGjIs the inertia time constant of the jth synchronous generator.
Wherein, the voltage amplitude U of the machine end of the fanDFIGiAnd the voltage amplitude U of the grid-connected point bussAnd the fan phase-locked angle theta at the initial disturbance momentpll0Initial phase angle theta of grid-connected point bus voltage at initial disturbance moments0Grid-connected point bus current amplitude IsInitial phase angle of phase-locked angle variation value
Figure BDA0002586229260000114
Terminal voltage amplitude U of synchronous generatorGjAnd the power angle delta of the synchronous generator at the initial disturbance moment0The operation data is the operation data of the wind power grid-connected system; the phase angle disturbed fluctuation amplitude a and the disturbed fluctuation angular velocity omega are oscillation mode data; proportional parameter K of phase-locked loop of fanpiIntegral parameter K of phase-locked loop of fanIiAnd virtual inertia control parameter KωThe control parameters are control parameters of the double-fed wind generating set.
When the oscillation of the wind power grid-connected system converges, the voltage amplitude U of the grid-connected point busSGradually decreases, at the moment, the non-periodic energy variation of the wind power grid-connected system decreases along with time, namely
Figure BDA0002586229260000111
The accumulation of total non-periodic dynamic energy in unit time of the wind power grid-connected system is represented as a negative value, namely the system continuously consumes disturbance injected energy, and finally the stability is achieved; when the oscillation of the wind power grid-connected system is dispersed, the voltage amplitude U of the grid-connected point busSGradually increases, and the non-periodic energy variation of the wind power grid-connected system increases along with time at the moment, namely
Figure BDA0002586229260000112
Representing total non-periodic dynamic energy of system in single timeThe accumulation of (1) is a positive value, namely the system energy is continuously increased and finally is unstable; when the wind power grid-connected system is in constant-amplitude oscillation, namely
Figure BDA0002586229260000113
The total non-periodic dynamic energy accumulation in the unit time of the representation system is 0, namely the energy consumed by the system is just equal to the energy injected into the system, and at the moment, the wind power grid-connected system is in a critical stable state.
Preferably, the system modal stability domain construction module constructs a modal stability domain of the wind power grid-connected system in the following manner:
a small disturbance stability region is constructed in a node active power injection space of a wind power grid-connected system and serves as the bottom surface of a modal stability region, and within the range of the bottom surface, the wind power grid-connected system has the capability of bearing certain disturbance under the current operation working condition, so that the oscillation instability phenomenon cannot occur.
The node active power injection space of the wind power grid-connected system can be expressed as follows:
ΔWh＝f(P1,P2,P3,…,Pk)，
and k represents a node, namely the active power injected by each node is an independent variable, the total non-periodic dynamic energy is a dependent variable, and the active power injected by the node is the operation data of the wind power grid-connected system.
Applying different energy disturbance, P, to any operating point of the wind power grid-connected system in the bottom surface1,P2,P3,…,PkThe operation point is subjected to oscillation mode analysis to obtain the modal frequency and amplitude information of the operation point under each oscillation mode, whether the operation point is a critical stable point under the corresponding oscillation mode is judged according to the modal stability criterion based on the modal frequency and amplitude information of the operation point under each oscillation mode, and if the operation point is the critical stable point, the total non-periodic dynamic energy value of the unit time of the wind power grid-connected system at the operation point is obtained to be used as the critical point on the boundary of the modal stable domain under the corresponding oscillation mode.
Wherein the mode stability criterion is as follows:
Figure BDA0002586229260000121
according to the modal frequency and amplitude information, finding a scene corresponding to the non-periodic energy variation of each mode from negative to positive in unit time in two sections of unit time before and after the operation point, namely a scene corresponding to the modal stability criterion, and if the modal stability criterion is met in a certain mode, judging that the operation point is a critical stability point of the wind power grid-connected system in the mode.
Traversing all the operating points of the wind power grid-connected system in the bottom surface to obtain a critical point on the boundary of the modal stability domain under each oscillation mode, and fitting the critical point under each oscillation mode to obtain the modal stability domain corresponding to each mode.
Preferably, the modal stability domain boundary fitting module performs hyperplane fitting on the modal stability domain of the wind power grid-connected system by using the following method:
fitting and obtaining a hyperplane of a modal stability domain of the wind power grid-connected system under each oscillation mode according to the following formula:
Figure BDA0002586229260000122
estimating the fitting accuracy of the obtained hyperplane according to the following formula:
Figure BDA0002586229260000123
wherein, PiActive power, alpha, injected into the ith node of the wind power grid-connected systemiFor the hyperplane parameter, i is 1,2ESRThe fitting accuracy of the hyperplane is obtained.
And when the fitting precision of the obtained hyperplane is greater than a preset threshold value, changing the parameters of the hyperplane for re-fitting until the fitting precision of the obtained hyperplane is less than the preset threshold value.
After the wind power grid-connected system is disturbed, the position of the system operation point in each modal stability domain changes. Accordingly, the distance between the system operating point and the boundary of the modal stability domain and the amount of disturbance energy that the system operating point can bear also vary: when the distance between the wind power grid-connected system and the modal stability domain boundary is farther within the modal stability domain range, the stability of the wind power grid-connected system under the mode is stronger, and the capacity of bearing disturbance is larger. Define the modal stability distance as:
Sm(P1m,P2m,...,Pkm)＝Sm0(P1m0,P2m0,...,Pkm0)-Sdst，
in the formula, Sm0(P1m0,P2m0,...,Pkm0) The distance from a disturbed operation point to a stable boundary of a dominant oscillation mode m along a certain specific direction is Sm(P1m,P2m,...,Pkm) And the size of the distance value can represent the stability margin of the disturbed wind power grid-connected system under the dominant oscillation mode m. Definition of SdstThe variation of the distance from the operation point before and after the disturbance to the stable boundary of the dominant oscillation mode m along a specific direction is adopted.
Each oscillation mode corresponds to one mode stability domain, specifically, a dominant oscillation mode is determined by the following method, for example, the operation point is located outside a mode stability domain of the oscillation mode 1, outside a mode stability domain of the oscillation mode 2, and within a mode stability domain of the oscillation mode 3, the operation point is unstable for the oscillation mode 1 and the oscillation mode 2, and is stable for the oscillation mode 3, and a critical distance from the operation point to the oscillation mode 1 is greater than a critical distance to the oscillation mode 2, at this time, for a wind power grid-connected system at the operation point, the oscillation mode 1 is the dominant oscillation mode, the optimization purpose is to optimize control parameters of a doubly-fed wind turbine, so that after disturbance, the operation point can be in the mode stability domain of the oscillation mode 1, and thus the grid-connected system can recover a safe and stable operation state.
When the operating point of the wind power grid-connected system is within the stable boundaryWhen S is presentm(P1m,P2m,...,Pkm) Is greater than 0; when the system operating point is on the stable boundary, Sm(P1m,P2m,...,Pkm) 0; when the system operating point is outside the stable boundary, Sm(P1m,P2m,...,Pkm)＜0。
Preferably, the objective function in the wind turbine generator control parameter optimization module is as follows:
Figure BDA0002586229260000141
wherein m represents the dominant oscillation mode of the wind power grid-connected system under certain disturbance, PimActive power alpha injected into the ith node of the wind power grid-connected system under m dominant oscillation modesimAnd (3) corresponding hyperplane parameters under the m dominant oscillation mode, wherein i is 1,2, k is the number of nodes of the wind power grid-connected system.
In particular, in practice, Sm(P1m,P2m,...,Pkm) The wind power grid-connected system is operated in a safe state by meeting a certain stability margin, namely a certain dynamic stability capacity is achieved, so that the stability margin constraint condition is met:
Sm(P1m,P2m,...,Pkm)≥H0。
preferably, the constraint conditions of the tidal current static safety and stability condition are as follows:
Figure BDA0002586229260000142
the frequency range constraint conditions of the wind power grid-connected system are as follows:
Figure BDA0002586229260000143
wherein H0A preset modal stability distance, m is a dominant oscillation mode, J (K) is a wind power generation modeVirtual inertia of the net system, JsetTo preset virtual inertia, PimActive power injected into the ith node of the wind power grid-connected system under the m-dominant oscillation mode, wherein i is 1,2P_PLLFor proportional control parameters of the phase-locked loop, KI_PLLFor integral control parameters of phase-locked loops, KωFor virtual inertia control parameters, Ω1Is KP_PLLValue range of (1), omega2Is KI_PLLValue range of (1), omega3Is KωThe value range of (a).
Method embodiment
The invention further discloses a fan parameter optimization oscillation suppression method based on a modal stability domain, which is used for suppressing oscillation of a wind power grid-connected system by adjusting control parameters of a double-fed wind generating set. As shown in fig. 2, the method includes:
step 1, collecting operation data, oscillation mode data and control parameters of a double-fed wind generating set of a wind power grid-connected system.
And 2, obtaining the total aperiodic dynamic energy of the wind power grid-connected system by utilizing a system aperiodic dynamic energy model based on the operation data and the oscillation mode data of the wind power grid-connected system and the control parameters of the doubly-fed wind generating set.
And 3, constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the obtained total non-periodic dynamic energy of the wind power grid-connected system.
And 4, performing hyperplane fitting on the modal stability region of the wind power grid-connected system to obtain a fitting hyperplane.
And 5, optimizing control parameters of the doubly-fed wind generating set based on the fit-obtained hyperplane by taking the maximum modal stability distance of the wind power grid-connected system as a target function and the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system as constraint conditions.
Preferably, in step 3, a modal stability region of the wind power grid-connected system is specifically constructed in the following manner:
a small disturbance stability region is constructed in a node active power injection space of a wind power grid-connected system and serves as the bottom surface of a modal stability region, and within the range of the bottom surface, the wind power grid-connected system has the capability of bearing certain disturbance under the current operation working condition, so that the oscillation instability phenomenon cannot occur.
The node active power injection space of the wind power grid-connected system can be expressed as follows:
ΔWh＝f(P1,P2,P3,…,Pk)，
k represents nodes, namely active power injected by each node is an independent variable, and total non-periodic dynamic energy is a dependent variable.
Applying different energy disturbance, P, to any operating point of the wind power grid-connected system in the bottom surface1,P2,P3,…,PkThe operation point is subjected to oscillation mode analysis to obtain the modal frequency and amplitude information of the operation point under each oscillation mode, whether the operation point is a critical stable point under the corresponding oscillation mode is judged according to the modal stability criterion based on the modal frequency and amplitude information of the operation point under each oscillation mode, and if the operation point is the critical stable point, the total non-periodic dynamic energy value of the unit time of the wind power grid-connected system at the operation point is obtained to be used as the critical point on the boundary of the modal stable domain under the corresponding oscillation mode.
Wherein the mode stability criterion is as follows:
Figure BDA0002586229260000161
according to the modal frequency and amplitude information, finding a scene corresponding to the non-periodic energy variation of each mode from negative to positive in unit time in two sections of unit time before and after the operation point, namely a scene corresponding to the modal stability criterion, and if the modal stability criterion is met in a certain mode, judging that the operation point is a critical stability point of the wind power grid-connected system in the mode.
Traversing all the operating points of the wind power grid-connected system in the bottom surface to obtain a critical point on the boundary of the modal stability domain under each oscillation mode, and fitting the critical point under each oscillation mode to obtain the modal stability domain corresponding to each mode.
Preferably, in step 4, the hyperplane fitting is performed on the modal stability domain of the wind power grid-connected system specifically by the following method:
fitting and obtaining a hyperplane of a modal stability domain of the wind power grid-connected system under each oscillation mode according to the following formula:
Figure BDA0002586229260000162
estimating the fitting accuracy of the obtained hyperplane according to the following formula:
Figure BDA0002586229260000163
wherein, PiActive power, alpha, injected into the ith node of the wind power grid-connected systemiIs a hyperplane parameter, i is 1,2, k is the number of nodes of the wind power grid-connected system, and σ isESRAnd the fitting precision of the hyperplane is obtained.
And when the fitting precision of the hyperplane is greater than a preset threshold value, changing the hyperplane parameters to perform refitting until the obtained fitting precision of the hyperplane is less than the preset threshold value.
Preferably, the objective function is:
Figure BDA0002586229260000164
the stability margin constraint conditions are as follows:
Sm(P1m,P2m,...,Pkm)≥H0，
the constraint conditions of the tidal current static safety and stability condition are as follows:
Figure BDA0002586229260000171
the frequency range constraint conditions of the wind power grid-connected system are as follows:
Figure BDA0002586229260000172
wherein m represents the dominant oscillation mode of the wind power grid-connected system under certain disturbance, PimActive power alpha injected into ith node of wind power grid-connected system under m dominant oscillation modesimCorresponding hyperplane parameters under m dominant oscillation modes, i is 1,2, k is the node number of the wind power grid-connected system, and H is0For a predetermined modal stability distance, J (K) is the virtual inertia of the grid-connected system, JsetTo preset virtual inertia, KP_PLLFor proportional control parameters of the phase-locked loop, KI_PLLFor integral control parameters of phase-locked loops, KωFor virtual inertia control parameters, Ω1Is KP_PLLValue range of (1), omega2Is KI_PLLValue range of (1), omega3Is KωThe value range of (a).
The advantageous effects of the present invention are now better demonstrated by the following examples.
The wind power grid-connected system based on the IEEE 4 machine 11 node verifies the energy model and the dynamic characteristics of the doubly-fed wind turbine generator, and the structure of the wind power grid-connected system is shown in figure 3. Wherein, G1, G3 and G4 are synchronous generators, active power is 700MW, and G3 is a reference unit. G2 is replaced by a wind power plant DFIG composed of equal-capacity double-fed wind turbine generators, the wind power plant is formed by connecting 350 double-fed wind turbine generators in parallel, the rated capacity of each double-fed wind turbine generator is 2MW, and the total active power of the wind power plant is 700 MW. And the buses 7 and 9 are respectively connected with loads L1 and L2, and the total load of the wind power grid-connected system is 2734 MW.
1) Verifying stable working condition of wind power grid-connected system
The control parameters of the doubly-fed wind turbine generator set are as follows: kP_PLL＝1p.u.，KI_PLL＝1p.u.，Kω1p.u. Aiming at the wind power grid-connected system, the stability margin H0The set value of (c) is 0.1MW · s.
The bus 8 sends out in t-10 sThree-phase short circuit is generated, and the fault is cut off after the fault duration is 0.1 s. The optimized parameters are calculated as follows: kP_PLL＝0.85p.u.，KI_PLL＝0.17p.u.，Kω0.72p.u. Comparison of power-angle response curves of each unit of the wind power grid-connected system before and after optimization of control parameters of the doubly-fed wind turbine generator (taking G3 as a reference machine) is shown in FIGS. 4 and 5.
As can be seen from comparison between fig. 4 and fig. 5, in this scenario, the oscillation duration of the disturbed wind power grid-connected system before and after parameter optimization is shortened from 15s to about 8s, and the power angle swing amplitude is significantly reduced. Therefore, the power angle response characteristic of the system can be improved through parameter optimization, the wind power grid-connected system can be quickly recovered to be stable after disturbance disappears, and the stability of the system is improved.
2) Verification of instability working condition of wind power grid-connected system
In order to further verify the effectiveness of the parameter optimization strategy under different disturbance energies, a three-phase short-circuit fault occurs when t is 52s at the bus 8, and the fault duration is 0.8 s. The optimized parameters are calculated as follows: kP_PLL＝0.64p.u.，KI_PLL＝0.13p.u.，Kω0.55p.u. Active power curves of all the units of the system before and after optimization of the control parameters of the doubly-fed wind turbine generator are shown in fig. 6 and 7.
As can be seen from fig. 6, in this scenario, after the disturbance disappears, each generator is periodically destabilized and cannot automatically recover to a stable operation state. As can be seen from fig. 7, by parameter optimization, the system can be restored to a stable state after about 10s of oscillation after the disturbance disappears, the generators G1 and DFIG can be restored to the initial operating state, and the synchronous machine G4 is transitioned to a new operating state. The parameter optimization strategy can achieve the effect of suppressing oscillation.
Compared with the prior art, the system and the method for restraining the fan parameter optimization oscillation based on the modal stability domain provided by the invention have the advantages that firstly, the maximum modal stability distance of the wind power grid-connected system is taken as a target function, the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system are taken as constraint conditions, the control parameters of the double-fed wind generating set are optimized, the overall factors influencing the oscillation of the wind power grid-connected system are considered, the wind power grid-connected system can quickly recover the stable running state after the disturbance disappears through parameter optimization, and the problems of poor oscillation restraining effect and low restraining efficiency caused by only considering the influence of a single factor on the oscillation of the wind power grid-connected system are avoided; secondly, the influence of the parameters of the fan on the oscillation of the wind power grid-connected system is used as a constraint condition for optimizing the control parameters of the double-fed wind generating set by the fan parameter optimization oscillation suppression system and the fan parameter optimization oscillation suppression method, so that the problem that the wind generating set is possibly in a low damping level or a negative damping state and the stable operation of a power grid is threatened when the self-set parameters are unreasonable is solved.
Those skilled in the art will appreciate that all or part of the flow of the method implementing the above embodiments may be implemented by a computer program, which is stored in a computer readable storage medium, to instruct related hardware. The computer readable storage medium is a magnetic disk, an optical disk, a read-only memory or a random access memory.
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.

Claims (9)

1. The utility model provides a fan parameter optimization oscillation suppression system based on mode stability territory which characterized in that for carry out the oscillation suppression to wind-powered electricity generation grid-connected system through adjusting double-fed wind generating set's control parameter, fan parameter optimization oscillation suppression system includes:
the data acquisition module is used for acquiring the operation data, the oscillation mode data and the control parameters of the double-fed wind generating set of the wind power grid-connected system;
the system aperiodic dynamic energy acquisition module is used for acquiring the total aperiodic dynamic energy of the wind power grid-connected system by utilizing an aperiodic dynamic energy model based on the operation data and the oscillation modal data of the wind power grid-connected system and the control parameters of the doubly-fed wind generating set;
the system modal stability domain construction module is used for constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the obtained total aperiodic dynamic energy of the wind power grid-connected system;
the modal stability domain boundary fitting module is used for performing hyperplane fitting on a modal stability domain of the wind power grid-connected system to obtain a fitting hyperplane;
the wind turbine generator control parameter optimization module is used for optimizing the control parameters of the doubly-fed wind turbine generator based on the fitting obtained hyperplane, the maximum modal stability distance of the wind power grid-connected system is a target function, and the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system are used as constraint conditions;
the non-periodic dynamic energy model of the system is as follows:
Figure FDA0003196474150000011
Figure FDA0003196474150000012
Figure FDA0003196474150000013
wherein, UsIs the grid-connected point bus voltage amplitude, U, of the wind power grid-connected systemDFIGIs the terminal voltage amplitude, x, of the doubly-fed wind generatorDFIGIs the equivalent impedance, theta, of a doubly-fed wind generatorpll0The doubly-fed wind generator phase locking angle theta corresponding to the initial disturbance moments0The initial phase angle of the grid-connected point bus voltage at the initial disturbance moment, a is the fluctuation amplitude of the grid-connected point bus voltage after the phase angle disturbance, omega is the fluctuation angular velocity of the grid-connected point bus voltage after the phase angle disturbance, t is time, IsElectricity for grid-connected point busThe magnitude of the flow is such that,
Figure FDA0003196474150000021
is the initial phase angle of the phase-locked angle, KωFor controlling parameters for virtual inertia, KpThe proportional parameter of a phase-locked loop in the wind power grid-connected system is obtained; kIIs an integration parameter of the phase locked loop.
2. The fan parameter optimization oscillation suppression system according to claim 1, wherein the system modal stability domain construction module constructs a modal stability domain of the wind power grid-connected system by:
constructing a small disturbance stability domain in a node active power injection space of the wind power grid-connected system as the bottom surface of a modal stability domain;
applying disturbance with different energies to an operation point of any wind power grid-connected system in the bottom surface, carrying out oscillation mode analysis on the operation point to obtain modal frequency and amplitude information of the operation point under each oscillation mode, judging whether the operation point is a critical stable point under a corresponding oscillation mode according to modal stability criteria based on the modal frequency and amplitude information of the operation point under each oscillation mode, and if so, acquiring a unit-time total non-periodic dynamic energy value of the wind power grid-connected system at the operation point as a critical point on a modal stable domain boundary under the corresponding oscillation mode;
traversing all the operating points of the wind power grid-connected system in the bottom surface to obtain a critical point on the boundary of the modal stability domain under each oscillation mode, and fitting the critical point under each oscillation mode to obtain the modal stability domain corresponding to each mode.
3. The fan parameter optimization oscillation suppression system of claim 1, wherein the modal stability domain boundary fitting module performs hyperplane fitting on a modal stability domain of the wind power grid-connected system by:
fitting and obtaining a hyperplane of a modal stability domain of the wind power grid-connected system under each oscillation mode according to the following formula:
Figure FDA0003196474150000022
estimating the fitting accuracy of the obtained hyperplane according to the following formula:
Figure FDA0003196474150000023
wherein, PiActive power, alpha, injected into the ith node of the wind power grid-connected systemiIs a hyperplane parameter, i is 1,2, k is the number of nodes of the wind power grid-connected system, and σ isESRFitting accuracy of the hyperplane;
and when the fitting precision of the hyperplane is greater than a preset threshold value, changing the hyperplane parameters to perform refitting until the obtained fitting precision of the hyperplane is less than the preset threshold value.
4. The wind turbine parameter optimizing oscillation suppression system of claim 1 or 3, wherein the objective function in the wind turbine control parameter optimizing module is:
Figure FDA0003196474150000031
wherein m represents a dominant oscillation mode of the wind power grid-connected system, S represents a distance from a current operation state of the wind power grid-connected system to a stable domain boundary of the dominant oscillation mode under the dominant oscillation mode, and P represents a distance between the current operation state of the wind power grid-connected system and the stable domain boundary of the dominant oscillation modeimActive power alpha injected into the ith node of the wind power grid-connected system under m dominant oscillation modesimAnd (3) corresponding hyperplane parameters under the m dominant oscillation mode, wherein i is 1,2, k is the number of nodes of the wind power grid-connected system.
5. The fan parameter optimization oscillation suppression system of claim 4, wherein the stability margin constraints are:
Sm(P1m,P2m,...,Pkm)≥H0，
the constraint conditions of the tidal current static safety and stability condition are as follows:
Figure FDA0003196474150000032
the frequency range constraint conditions of the wind power grid-connected system are as follows:
Figure FDA0003196474150000033
wherein H0M is a dominant oscillation mode for a preset mode stable distance, J (K) is the virtual inertia of the wind power grid-connected system, JsetTo preset virtual inertia, PimActive power injected into the ith node of the wind power grid-connected system under the m-dominant oscillation mode, i is 1,2, …, K is the number of the nodes of the wind power grid-connected system, and K isP_PLLFor proportional control parameters of the phase-locked loop, KI_PLLFor integral control parameters of phase-locked loops, KωFor virtual inertia control parameters, Ω1Is KP_PLLValue range of (1), omega2Is KI_PLLValue range of (1), omega3Is KωThe value range of (a).
6. A fan parameter optimization oscillation suppression method based on a modal stability domain is characterized by being used for carrying out oscillation suppression on a wind power grid-connected system by adjusting control parameters of a doubly-fed wind generating set, and comprising the following steps:
acquiring operation data, oscillation mode data and control parameters of a double-fed wind generating set of a wind power grid-connected system;
acquiring total aperiodic dynamic energy of the wind power grid-connected system by utilizing a system aperiodic dynamic energy model based on the operation data and the oscillation mode data of the wind power grid-connected system and the control parameters of the double-fed wind generating set;
constructing a modal stability domain of the wind power grid-connected system according to the operation data of the wind power grid-connected system and the obtained total aperiodic dynamic energy of the wind power grid-connected system;
performing hyperplane fitting on a modal stability region of the wind power grid-connected system to obtain a fitting hyperplane;
optimizing control parameters of the double-fed wind generating set by taking the maximum modal stability distance of the wind power grid-connected system as a target function and the frequency range, the stability margin and the static safety and stability condition of the power flow of the wind power grid-connected system as constraint conditions on the basis of the fit obtained hyperplane;
the non-periodic dynamic energy model of the system is as follows:
Figure FDA0003196474150000041
Figure FDA0003196474150000042
Figure FDA0003196474150000043
wherein, UsIs the grid-connected point bus voltage amplitude, U, of the wind power grid-connected systemDFIGIs the terminal voltage amplitude, x, of the doubly-fed wind generatorDFIGIs the equivalent impedance, theta, of a doubly-fed wind generatorpll0The doubly-fed wind generator phase locking angle theta corresponding to the initial disturbance moments0The initial phase angle of the grid-connected point bus voltage at the initial disturbance moment, a is the fluctuation amplitude of the grid-connected point bus voltage after the phase angle disturbance, omega is the fluctuation angular velocity of the grid-connected point bus voltage after the phase angle disturbance, t is time, IsIs the current amplitude of the grid-connected point bus,
Figure FDA0003196474150000051
is the initial phase angle of the phase-locked angle, KωFor controlling parameters for virtual inertia, KpThe proportional parameter of a phase-locked loop in the wind power grid-connected system is obtained; kIIs an integration parameter of the phase locked loop.
7. The fan parameter optimization oscillation suppression method according to claim 6, wherein a modal stability domain of the wind power grid-connected system is constructed by:
constructing a small disturbance stability domain in a node active power injection space of the wind power grid-connected system as the bottom surface of a modal stability domain;
applying disturbance with different energies to an operation point of any wind power grid-connected system in the bottom surface, carrying out oscillation mode analysis on the operation point to obtain modal frequency and amplitude information of the operation point under each oscillation mode, judging whether the operation point is a critical stable point under a corresponding oscillation mode according to modal stability criteria based on the modal frequency and amplitude information of the operation point under each oscillation mode, and if so, acquiring a unit-time total non-periodic dynamic energy value of the wind power grid-connected system at the operation point as a critical point on a modal stable domain boundary under the corresponding oscillation mode;
traversing all the operating points of the wind power grid-connected system in the bottom surface to obtain a critical point on the boundary of the modal stability domain under each oscillation mode, and fitting the critical point under each oscillation mode to obtain the modal stability domain corresponding to each mode.
8. The fan parameter optimization oscillation suppression method according to claim 6, wherein the hyperplane fitting is performed on the modal stability domain of the wind power grid-connected system by:
fitting and obtaining a hyperplane of a modal stability domain of the wind power grid-connected system under each oscillation mode according to the following formula:
Figure FDA0003196474150000052
estimating the fitting accuracy of the obtained hyperplane according to the following formula:
Figure FDA0003196474150000053
wherein, PiActive power, alpha, injected into the ith node of the wind power grid-connected systemiIs a hyperplane parameter, i is 1,2, k is the number of nodes of the wind power grid-connected system, and σ isESRFitting accuracy of the hyperplane;
and when the fitting precision of the hyperplane is greater than a preset threshold value, changing the hyperplane parameters to perform refitting until the obtained fitting precision of the hyperplane is less than the preset threshold value.
9. The method for suppressing oscillation in a wind turbine according to any of claims 6 to 8, wherein the objective function is:
Figure FDA0003196474150000061
the stability margin constraint conditions are as follows:
Sm(P1m,P2m,...,Pkm)≥H0，
the constraint conditions of the tidal current static safety and stability condition are as follows:
Figure FDA0003196474150000062
the frequency range constraint conditions of the wind power grid-connected system are as follows:
Figure FDA0003196474150000063
wherein m represents a dominant oscillation mode of the wind power grid-connected system, S represents a distance from a current operation state of the wind power grid-connected system to a stable domain boundary of the dominant oscillation mode under the dominant oscillation mode, and P represents a distance between the current operation state of the wind power grid-connected system and the stable domain boundary of the dominant oscillation modeimActive power alpha injected into the ith node of the wind power grid-connected system under m dominant oscillation modesimCorresponding hyperplane parameters under m-dominant oscillation modes, i is 1,2, …, k is the number of nodes of the wind power grid-connected system, and H is0A preset modal stability distance, J (K) is a virtual inertia of the wind power grid-connected system, JsetTo preset virtual inertia, KP_PLLFor proportional control parameters of the phase-locked loop, KI_PLLFor integral control parameters of phase-locked loops, KωFor virtual inertia control parameters, Ω1Is KP_PLLValue range of (1), omega2Is KI_PLLValue range of (1), omega3Is KωThe value range of (a).
CN202010682091.7A 2020-07-15 2020-07-15 Fan parameter optimization oscillation suppression system and method based on modal stability domain Active CN111769575B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010682091.7A CN111769575B (en) 2020-07-15 2020-07-15 Fan parameter optimization oscillation suppression system and method based on modal stability domain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010682091.7A CN111769575B (en) 2020-07-15 2020-07-15 Fan parameter optimization oscillation suppression system and method based on modal stability domain

Publications (2)

Publication Number Publication Date
CN111769575A CN111769575A (en) 2020-10-13
CN111769575B true CN111769575B (en) 2021-09-28

Family

ID=72725776

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010682091.7A Active CN111769575B (en) 2020-07-15 2020-07-15 Fan parameter optimization oscillation suppression system and method based on modal stability domain

Country Status (1)

Country Link
CN (1) CN111769575B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113569464A (en) * 2021-06-21 2021-10-29 国网山东省电力公司电力科学研究院 Wind turbine generator oscillation mode prediction method and device based on deep learning network and multi-task learning strategy
CN113783232B (en) * 2021-11-09 2022-02-11 中国电力科学研究院有限公司 Method and device for monitoring synchronization safety and stability of new energy grid-connected system
CN114759618B (en) * 2021-12-31 2024-08-06 华北电力大学 Method and system for determining stability of doubly-fed wind power grid-connected system based on dynamic energy
CN114759578B (en) * 2022-02-25 2024-08-06 华北电力大学 Method for constructing low-frequency oscillation suppression model of wind power grid-connected system

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104680323B (en) * 2015-03-02 2018-05-25 华北电力大学 Electric system energy stabilization domain builds system and method
CN105186502B (en) * 2015-09-18 2017-11-28 天津大学 The power system transient stability analysis method of blower fan containing double-fed based on security domain
CN109193752B (en) * 2018-08-22 2021-01-08 华北电力大学 Method for optimizing low-frequency oscillation control parameters of virtual inertia-containing doubly-fed wind turbine grid-connected system
CN109713685B (en) * 2018-11-07 2021-01-29 华北电力大学 Online positioning method suitable for VSC access-induced subsynchronous oscillation
CN110350588B (en) * 2019-07-30 2020-09-04 华北电力大学 Dynamic energy stability evaluation method and system for doubly-fed fan grid-connected system
CN110350587B (en) * 2019-07-30 2020-09-04 华北电力大学 Doubly-fed fan subsynchronous oscillation suppression method and system based on converter control
CN110311414B (en) * 2019-07-30 2020-11-20 华北电力大学 Doubly-fed wind turbine control parameter optimization method and system based on energy stability domain
CN110417054B (en) * 2019-07-30 2020-10-16 华北电力大学 System and method for adjusting stability of double-fed wind power grid-connected system
CN110309625B (en) * 2019-07-30 2021-09-07 华北电力大学 Method and system for determining energy stability domain of double-fed wind power grid-connected system

Also Published As

Publication number Publication date
CN111769575A (en) 2020-10-13

Similar Documents

Publication Publication Date Title
CN111769575B (en) Fan parameter optimization oscillation suppression system and method based on modal stability domain
Aguilar et al. Multi-objective PSO applied to PI control of DFIG wind turbine under electrical fault conditions
US12119649B2 (en) Oscillation active damping control method and system for grid-tied type-4 wind turbine generator
CN111864767B (en) Doubly-fed wind turbine generator subsynchronous oscillation active damping control method and system
Li et al. Energy-shaping L2-gain controller for PMSG wind turbine to mitigate subsynchronous interaction
CN110311414B (en) Doubly-fed wind turbine control parameter optimization method and system based on energy stability domain
Liu et al. MMC-STATCOM supplementary wide-band damping control to mitigate subsynchronous control interaction in wind farms
Bhukya et al. Mathematical modelling and stability analysis of PSS for damping LFOs of wind power system
CN112993991A (en) Fan two-channel damping control low-frequency oscillation wide-area damping control method
Senapati et al. Lagrange interpolating polynomial–based deloading control scheme for variable speed wind turbines
Wu et al. Mitigating subsynchronous oscillation using model-free adaptive control of DFIGs
Shabani et al. Investigation of the closed-loop control system on the DFIG dynamic models in transient stability studies
de Oliveira et al. Synchronverter-based frequency control technique applied in wind energy conversion systems based on the doubly-fed induction generator
Ma et al. A low frequency oscillation suppression method for grid-connected DFIG with virtual inertia
Ebrahim et al. Transient stability enhancement of a wind energy distributed generation system by using fuzzy logic stabilizers
Wang et al. Linear parameter varying control of a doubly fed induction generator based wind turbine with primary grid frequency support
Hemeyine et al. Robust takagi sugeno fuzzy models control for a variable speed wind turbine based a DFI-generator
CN112290590B (en) PMSG fan control method based on super-capacitor virtual inertia and short-term primary frequency modulation
CN115800296B (en) Voltage frequency collaborative supporting method for open sea wind power through VSC-MTDC grid-connected system
Shanoob et al. Wind turbine transient response and fault ride‐through improvements with optimal control
CN116418011A (en) Net-structured fan control method based on self-adaptive power correction compensation
Guo et al. Subsynchronous oscillation suppression strategy of wind turbine’s phase-locked loop based on optimal damping ratio
CN114825425A (en) New energy acceptance capacity assessment method and device for voltage drop induced frequency safety
Yuan et al. Active power control strategy of DFIGs in flexible interconnected distribution network
Abdollahi Chirani et al. Investigation of the Impact of SSSC‐Based FLC on the Stability of Power Systems Connected to Wind Farms

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