CN114336685A - Oscillation mode identification method for wind power transmission system through MMC-HVDC - Google Patents

Oscillation mode identification method for wind power transmission system through MMC-HVDC Download PDF

Info

Publication number
CN114336685A
CN114336685A CN202210020688.4A CN202210020688A CN114336685A CN 114336685 A CN114336685 A CN 114336685A CN 202210020688 A CN202210020688 A CN 202210020688A CN 114336685 A CN114336685 A CN 114336685A
Authority
CN
China
Prior art keywords
mmc
wind power
hvdc
model
mathematical model
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.)
Pending
Application number
CN202210020688.4A
Other languages
Chinese (zh)
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.)
Huaneng Rudong Baxianjiao Offshore Wind Power Co ltd
Huaneng Clean Energy Research Institute
North China Electric Power University
Huaneng Group Technology Innovation Center Co Ltd
Original Assignee
Huaneng Clean Energy Research Institute
North China Electric Power University
Huaneng Group Technology Innovation Center Co Ltd
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 Huaneng Clean Energy Research Institute, North China Electric Power University, Huaneng Group Technology Innovation Center Co Ltd filed Critical Huaneng Clean Energy Research Institute
Priority to CN202210020688.4A priority Critical patent/CN114336685A/en
Publication of CN114336685A publication Critical patent/CN114336685A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Control Of Eletrric Generators (AREA)
  • Wind Motors (AREA)

Abstract

The invention relates to an oscillation mode identification method and system for a wind power transmission system through MMC-HVDC. The method comprises the following steps: establishing mathematical models of all parts of a system for sending wind power out through MMC-HVDC and connecting to obtain an integrated mathematical model of the system; carrying out linearization processing on the integrated mathematical model to obtain a system small interference model; extracting a characteristic matrix in the small interference model for analysis, obtaining all the modes and corresponding state parameters of the system, carrying out normalization processing, and locking main participation links under the corresponding modes according to the state parameters with the highest correlation degree, thereby positioning an oscillation source and determining an oscillation mode; the oscillation mode includes wind farm internal, wind farm control system, MMC internal, MMC control system and interaction. The method covers all possible modes of the system by establishing a detailed mathematical model of the wind turbine generator and the MMC converter station, and effectively identifies and classifies the oscillation mode of the wind power transmitted out of the system through the MMC-HVDC.

Description

Oscillation mode identification method for wind power transmission system through MMC-HVDC
Technical Field
The invention relates to the technical field of power transmission and distribution, in particular to an oscillation mode identification method for a system for sending wind power out through MMC-HVDC.
Background
In order to solve the problems of global climate change, environmental pollution, fossil energy depletion and the like, energy structure transformation is actively promoted in various countries, and renewable energy power generation represented by wind power is an important component of a future power system. In recent years, Chinese wind power is rapidly developed, and the installed proportion of the wind power is continuously increased. Among various modes of wind power transmission, flexible direct current transmission is widely applied to the aspect of large-scale wind power transmission due to the advantages of economy, capability of providing alternating voltage and frequency support for a weak power grid or a passive system and the like. The modular multilevel converter based high voltage direct current (MMC-HVDC) has many advantages, such as low harmonic content, no commutation failure, active and reactive decoupling control, and is the mainstream topological structure of the current flexible direct current transmission system.
Multiple oscillation phenomena occur in existing wind power sending projects which are built at home and abroad and put into operation, including a Hami wind power plant, a Hebei Staphylou wind power plant, a south Virginian flexible and straight project and a south Australia multi-end flexible and straight project in China, and subsynchronous oscillation phenomena of different degrees occur, so that corresponding economic losses, such as tripping of a thermal power generating unit, damage of a crowbar circuit of a fan, system shutdown and the like, are caused. Therefore, the method has great significance for identifying and classifying the oscillation modes of the wind power plant through the MMC-HVDC sending-out system.
Recent research shows that the oscillation phenomenon easily occurs when a wind power plant is connected with a high series compensation line, a weak alternating current power grid and other power electronic equipment. MMC-HVDC has the problem of multi-band oscillation, wherein the middle-low frequency band oscillation has a direct relation with a circulating current suppression controller, inner and outer ring control parameters of a main controller, particularly the bandwidth of a phase-locked loop, and the high frequency band oscillation is related with the link delay of the controller. In summary, the wind power has oscillation problem through the MMC-HVDC transmission system, so it is necessary to provide a suitable oscillation mode identification method to identify the oscillation mode of the wind power through the MMC-HVDC transmission system.
Disclosure of Invention
The invention aims to provide an oscillation mode identification method of a system for sending wind power out through MMC-HVDC, so as to effectively identify and classify the oscillation mode of the system.
In order to achieve the purpose, the invention provides the following scheme:
an oscillation mode identification method for wind power through an MMC-HVDC sending system comprises the following steps:
establishing a mathematical model of each part of a wind power sending system through MMC-HVDC;
connecting the wind power through mathematical models of all parts of the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind power through the MMC-HVDC sending-out system;
carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out of the MMC-HVDC system;
extracting a characteristic matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through MMC-HVDC;
normalizing each state parameter, and locking a main participation link under a corresponding mode according to the state parameter with the highest correlation degree;
positioning an oscillation source and determining an oscillation mode according to the main participation link; the oscillation modes comprise a wind power plant internal oscillation mode, a wind power plant control system oscillation mode, an MMC internal oscillation mode, an MMC control system oscillation mode and an interaction oscillation mode.
Optionally, the establishing a mathematical model of each part of the wind power output system through MMC-HVDC specifically includes:
dividing the wind power into a wind power generator set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system, respectively modeling each part independently, and establishing a mathematical model of each part of the wind power sending-out system through the MMC-HVDC
Figure BDA0003462425080000021
Wherein X is the state variable of the DC power transmission system, U is the input variable,
Figure BDA0003462425080000022
is the first derivative of the state variable, f () is the first derivative of the state variable
Figure BDA0003462425080000023
A functional relationship with the state variable X and the input variable U; y is the output variable and Y () is the functional relationship between the output variable Y and the state variable X and the input variable U.
Optionally, the connecting the mathematical models of the parts of the MMC-HVDC transmission system of the wind power to obtain the integrated mathematical model of the parts of the MMC-HVDC transmission system of the wind power specifically includes:
and converting the electrical quantities of the wind turbine generator and the MMC converter station to the same coordinate system by adopting an interface conversion formula, and connecting the mathematical models of the wind power transmission system through MMC-HVDC transmission systems to obtain an integrated mathematical model of the wind power transmission system through MMC-HVDC transmission systems.
Optionally, the performing linearization on the integrated mathematical model to obtain a small interference model of the wind power transmitted through the MMC-HVDC transmission system specifically includes:
for the integrated mathematical modelCarrying out linearization processing to obtain a small interference model of the wind power output system through MMC-HVDC
Figure BDA0003462425080000031
Where A is the feature matrix, B is the input matrix, C is the output matrix, and D is the feedforward matrix.
Optionally, the extracting the feature matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through the MMC-HVDC system specifically includes:
extracting a characteristic matrix A in the small interference model for analysis, solving a characteristic value of the characteristic matrix A to obtain characteristic roots, wherein each group of characteristic roots represents a mode, and each mode corresponds to a state parameter; the state quantities comprise the state variables.
An oscillation mode identification system for wind power transmitted by an MMC-HVDC system comprises:
the mathematical model establishing module is used for establishing a mathematical model of each part of the wind power output system through the MMC-HVDC;
the mathematical model connecting module is used for connecting the mathematical models of all parts of the wind power output system through the MMC-HVDC to obtain an integrated mathematical model of the wind power output system through the MMC-HVDC;
the model linearization processing module is used for carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out of the MMC-HVDC system;
the model characteristic matrix analysis module is used for extracting a characteristic matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through the MMC-HVDC;
the state parameter normalization processing module is used for performing normalization processing on each state parameter and locking main participation links under corresponding modes according to the state parameter with the highest correlation degree;
the oscillation positioning and oscillation mode identification module is used for positioning an oscillation source and determining an oscillation mode according to the main participation link; the oscillation modes comprise a wind power plant internal oscillation mode, a wind power plant control system oscillation mode, an MMC internal oscillation mode, an MMC control system oscillation mode and an interaction oscillation mode.
Optionally, the mathematical model building module specifically includes:
a mathematical model establishing unit for dividing the wind power into a wind power set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system, respectively modeling each part separately, and establishing a mathematical model of each part of the wind power sending-out system through the MMC-HVDC
Figure BDA0003462425080000041
Wherein X is the state variable of the DC power transmission system, U is the input variable,
Figure BDA0003462425080000042
is the first derivative of the state variable, f () is the first derivative of the state variable
Figure BDA0003462425080000043
A functional relationship with the state variable X and the input variable U; y is the output variable and Y () is the functional relationship between the output variable Y and the state variable X and the input variable U.
Optionally, the mathematical model connection module specifically includes:
and the mathematical model connecting unit is used for converting the electrical quantities of the wind turbine generator and the MMC converter station to the same coordinate system by adopting an interface conversion formula, and connecting the mathematical models of the wind turbine generator and the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind turbine generator and the MMC-HVDC sending-out system.
Optionally, the model linearization processing module specifically includes:
the model linearization processing unit is used for carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out from the MMC-HVDC system
Figure BDA0003462425080000044
Wherein A is a feature matrix and B isThe input matrix, C the output matrix, and D the feedforward matrix.
Optionally, the model feature matrix analysis module specifically includes:
the model characteristic matrix analysis unit is used for extracting a characteristic matrix A in the small interference model for analysis, solving a characteristic value of the characteristic matrix A to obtain characteristic roots, wherein each group of characteristic roots represents a mode, and each mode corresponds to a state parameter; the state quantities comprise the state variables.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides an oscillation mode identification method and system for wind power transmitted out of a system through MMC-HVDC, wherein the method comprises the following steps: establishing a mathematical model of each part of a wind power sending system through MMC-HVDC; connecting the wind power through mathematical models of all parts of the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind power through the MMC-HVDC sending-out system; carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out of the MMC-HVDC system; extracting a characteristic matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through MMC-HVDC; normalizing each state parameter, and locking a main participation link under a corresponding mode according to the state parameter with the highest correlation degree; positioning an oscillation source and determining an oscillation mode according to the main participation link; the oscillation modes comprise a wind power plant internal oscillation mode, a wind power plant control system oscillation mode, an MMC internal oscillation mode, an MMC control system oscillation mode and an interaction oscillation mode. The method covers all possible modes of the system by establishing a detailed mathematical model of the wind turbine generator and the MMC converter station, effectively identifies and classifies the oscillation mode of the wind power transmitted out of the system through the MMC-HVDC, and the identification result can provide reference for system stability analysis.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of an oscillation mode identification method of a wind power transmission system through MMC-HVDC according to the present invention;
FIG. 2 is a schematic diagram of the method for identifying the oscillation mode of wind power transmitted through MMC-HVDC transmission system according to the present invention;
FIG. 3 is a schematic diagram illustrating a transformation of a wind farm and an MMC rotating coordinate system according to an embodiment of the present invention;
FIG. 4 is a structural diagram of an oscillation mode identification system of a wind power transmission system through MMC-HVDC according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide an oscillation mode identification method of a system for sending wind power out through MMC-HVDC, so as to effectively identify and classify the oscillation mode of the system.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a flow chart of an oscillation mode identification method of a system for sending wind power out through MMC-HVDC, and fig. 2 is a schematic principle diagram of the oscillation mode identification method of the system for sending wind power out through MMC-HVDC. Referring to fig. 1 and 2, the method for identifying the oscillation mode of the wind power transmission system through the MMC-HVDC system according to the present invention includes:
step 101: and establishing a mathematical model of each part of the wind power output system through MMC-HVDC.
The invention divides wind power into a wind power generator set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system (hereinafter referred to as the system for short), and models each part independently. Firstly, considering the detailed structure of the wind turbine generator and the internal harmonic characteristics of the MMC, determining a state variable and an input variable, respectively establishing a state space equation of the wind turbine generator, the MMC converter station and a connecting circuit, and obtaining a mathematical model of each part, wherein the mathematical model is shown as the following formula (1):
Figure BDA0003462425080000061
the steps of establishing the state space model in the three parts of the system are the same, and the formula (1) is a general expression of the state space equation. Wherein X is a state variable of the direct current transmission system; u is an input variable; y is an output variable;
Figure BDA0003462425080000062
for the first derivative of the state variable, f () represents the first derivative of the state variable as a function of the state variable and the input variable, and y () represents the output variable as a function of the state variable and the input variable.
Therefore, the step 101 of establishing a mathematical model of each part of the wind power transmission system through the MMC-HVDC transmission specifically includes:
dividing the wind power into a wind power generator set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system, respectively modeling each part independently, and establishing a mathematical model of each part of the wind power sending-out system through the MMC-HVDC
Figure BDA0003462425080000063
Wherein X is the state variable of the DC power transmission system, U is the input variable,
Figure BDA0003462425080000064
is the first derivative of the state variable, f () is the first derivative of the state variable
Figure BDA0003462425080000065
A functional relationship with the state variable X and the input variable U; y is the output variable and Y () is the functional relationship between the output variable Y and the state variable X and the input variable U.
Step 102: and connecting the wind power through mathematical models of all parts of the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind power through the MMC-HVDC sending-out system.
And establishing model interfaces of all parts, connecting all parts to obtain an integrated mathematical model, and performing linearization processing to obtain a system small interference model.
The process of establishing the model interface of each part is as follows:
and establishing different rotating coordinate systems determined by the wind power plant and the MMC phase-locked loop, wherein the coordinate conversion relation is shown in figure 3. In FIG. 3,. omega.pRepresenting the angular frequency, omega, of the phase-locked loop of a wind farm side converter1Representing the MMC phase-locked loop rotation angular frequency. dq is a coordinate system rotated at a constant angular frequency, with the index m (i.e. d)mqmCoordinate system) at an angular frequency ω1A rotating coordinate system with a phase angle theta11=ω1t), t represents time; subscript w (i.e. d)wqwCoordinate system) at an angular frequency ωpA rotating coordinate system with a phase angle θ p (θ p ═ ω [)pt)。θp-θ1Is the phase angle difference of the two rotating coordinate systems. i.e. isOutputting three-phase current i for wind farmsj(j ═ a, b, c) a spatially synthesized vector in the complex plane; i.e. isThe dq axis component under the coordinate system of the wind farm is isd,isq,isThe dq axis component under the coordinate system of the MMC converter station is Isd,Isq
The specific interface conversion formula is shown in the following formula (2):
Figure BDA0003462425080000071
wherein X represents an MMC side variable and X represents a wind farm side variable. x is the number ofsd、xsqRepresenting dq-axis components of a wind farm side variable x, e.g. isd、isq。Xsd、XsqRepresenting dq-axis components of MMC-side variable X, e.g. Isd、Isq
And (3) converting the electrical quantities of the wind power plant and the MMC converter station to the same coordinate system according to an interface conversion formula shown in the formula (2), so as to connect the mathematical models and obtain an integrated mathematical model of the whole system. The connection of the mathematical models means that the electrical quantities of the wind power plant and the MMC converter station are converted to the same coordinate system, and the process is essentially to convert state variables on a connection line so as to connect all parts.
Therefore, the step 102 of connecting the mathematical models of the parts of the MMC-HVDC transmission system of the wind power to obtain the integrated mathematical model of the parts of the MMC-HVDC transmission system of the wind power specifically includes:
and converting the electrical quantities of the wind turbine generator and the MMC converter station to the same coordinate system by adopting an interface conversion formula, and connecting the mathematical models of the wind power transmission system through MMC-HVDC transmission systems to obtain an integrated mathematical model of the wind power transmission system through MMC-HVDC transmission systems.
Step 103: and carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out of the MMC-HVDC system.
The integrated mathematical model is subjected to linearization processing to obtain a small interference model of the system, which is specifically shown in the following formula (3):
Figure BDA0003462425080000081
where A is the feature matrix, B is the input matrix, C is the output matrix, and D is the feedforward matrix. The expression form of the integrated mathematical model after linearization is shown as formula (3).
Therefore, the step 103 of performing linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted through the MMC-HVDC transmission system specifically includes:
for the integrated mathematical modelLine linearization processing is carried out to obtain a small interference model of the wind power output system through MMC-HVDC
Figure BDA0003462425080000082
Where A is the feature matrix, B is the input matrix, C is the output matrix, and D is the feedforward matrix.
Step 104: and extracting a characteristic matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through MMC-HVDC.
Based on the small interference model (3) in step 103, a stability characterization matrix a is extracted and analyzed to obtain all states of the system and oscillation characteristic information corresponding to the states. The stability characterization matrix A can represent the stability of the system, the state of the system can be obtained by analyzing the matrix A, the state of the system refers to characteristic roots obtained by solving characteristic values of the characteristic matrix A, and each group of characteristic roots represents a mode and contains oscillation characteristic information such as damping ratio, oscillation frequency and the like. Each state parameter corresponding to each modality refers to all parameters of the system, including system parameters and state variables. The system parameters comprise fan, MMC main circuit parameters, connecting circuit parameters and the like.
Therefore, the step 104 of extracting the feature matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through the MMC-HVDC system specifically includes:
extracting a characteristic matrix A in the small interference model for analysis, solving a characteristic value of the characteristic matrix A to obtain characteristic roots, wherein each group of characteristic roots represents a mode, and each mode corresponds to a state parameter; the state quantities comprise the state variables.
Step 105: and carrying out normalization processing on each state parameter, and locking main participation links under corresponding modes according to the state parameter with the highest correlation degree.
And (3) carrying out normalization processing on each state parameter, and locking main participation links (namely an electric quantity, an electric element or a certain molecular system and the like) in the mode according to the parameter with the highest degree of correlation. Each state parameter refers to all parameters of the system, including system parameters and state variables. There are many methods for determining the degree of participation (correlation), such as sensitivity analysis, which are not unique, and this belongs to basic knowledge in the field of power system stability research, and is not described herein again. The state parameter with the highest degree of correlation can be obtained, and the related control link or system structure can be positioned, namely the main participation link in the mode is locked. The main participation links in each mode can provide basis for positioning the oscillation source.
Step 106: and positioning an oscillation source and determining an oscillation mode according to the main participation link.
The oscillation mode is divided into 5 types of the interior of a wind power plant, the control system of the wind power plant, the interior of an MMC, the control system of the MMC and interaction through a main participation link to position an oscillation source. The main participation links are different, the corresponding main participation variables are also different, the participation variables are from a wind power plant and an MMC, and the generation source (namely the oscillation source) of the oscillation can be known to be a subsystem by judging which part the participation variables belong to, so that the corresponding oscillation mode is further determined. For example, the main participating variable obtained through analysis is a wind power plant side phase-locked loop control parameter, and then the oscillating wind power plant side is related to a wind power plant control link and belongs to an oscillating mode of a wind power plant control system.
The method can effectively identify the oscillation mode of the wind power transmitted out of the system through the MMC-HVDC through the step 101 and the step 106, and the identification result can provide reliable reference for the stability analysis of the subsequent system.
The invention provides an oscillation mode identification method of a system for sending wind power out through MMC-HVDC, which comprises the steps of dividing the system into a wind power generator set, an MMC converter station and a connecting circuit, and respectively carrying out independent modeling on each part; establishing an interface, connecting all parts to obtain an integrated mathematical model, and performing linearization processing to obtain a system small interference model; extracting a stability characterization matrix A, and analyzing the stability characterization matrix A to obtain all states of the system and corresponding oscillation characteristic information thereof; normalizing each state parameter, and locking main participation links (namely, an electric quantity, an electric element or a certain molecular system and the like) in the mode according to the parameter with the highest degree of correlation; the oscillation mode is divided into 5 types of the interior of a wind power plant, the control system of the wind power plant, the interior of an MMC, the control system of the MMC and interaction through a main participation link to position an oscillation source. The method covers all possible modes of the system by establishing a detailed mathematical model of the wind turbine generator and the MMC converter station, effectively identifies the oscillation mode of the wind power transmitted out of the system through the MMC-HVDC, and the identification result can provide reference for system stability analysis.
Based on the method for identifying the oscillation mode of the wind power output system through the MMC-HVDC provided by the invention, the invention also provides an identification system for identifying the oscillation mode of the wind power output system through the MMC-HVDC, and referring to fig. 4, the system comprises:
the mathematical model establishing module 401 is used for establishing a mathematical model of each part of the wind power output system through the MMC-HVDC;
the mathematical model connecting module 402 is used for connecting mathematical models of all parts of the MMC-HVDC sending-out system of the wind power to obtain an integrated mathematical model of the MMC-HVDC sending-out system of the wind power;
the model linearization processing module 403 is configured to perform linearization processing on the integrated mathematical model to obtain a small interference model of the wind power passing through the MMC-HVDC transmission system;
a model feature matrix analysis module 404, configured to extract a feature matrix in the small interference model for analysis, so as to obtain all modes and corresponding state parameters of the wind power transmitted from the MMC-HVDC transmission system;
a state parameter normalization processing module 405, configured to perform normalization processing on each state parameter, and lock a main participation link in a corresponding mode according to a state parameter with the highest correlation;
an oscillation positioning and oscillation mode identification module 406, configured to position an oscillation source and determine an oscillation mode according to the main participation link; the oscillation modes comprise a wind power plant internal oscillation mode, a wind power plant control system oscillation mode, an MMC internal oscillation mode, an MMC control system oscillation mode and an interaction oscillation mode.
The system identifies the oscillation mode of the wind power through the MMC-HVDC sending-out system, positions the oscillation source, and can divide the oscillation mode into 5 types of wind power plant interior, wind power plant control system, MMC interior, MMC control system and interaction.
The mathematical model building module 401 specifically includes:
a mathematical model establishing unit for dividing the wind power into a wind power set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system, respectively modeling each part separately, and establishing a mathematical model of each part of the wind power sending-out system through the MMC-HVDC
Figure BDA0003462425080000111
Wherein X is the state variable of the DC power transmission system, U is the input variable,
Figure BDA0003462425080000112
is the first derivative of the state variable, f () is the first derivative of the state variable
Figure BDA0003462425080000113
A functional relationship with the state variable X and the input variable U; y is the output variable and Y () is the functional relationship between the output variable Y and the state variable X and the input variable U.
The mathematical model connection module 402 specifically includes:
and the mathematical model connecting unit is used for converting the electrical quantities of the wind turbine generator and the MMC converter station to the same coordinate system by adopting an interface conversion formula, and connecting the mathematical models of the wind turbine generator and the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind turbine generator and the MMC-HVDC sending-out system.
The model linearization processing module 403 specifically includes:
the model linearization processing unit is used for carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out from the MMC-HVDC system
Figure BDA0003462425080000114
Wherein A is a feature matrix and B isThe input matrix, C the output matrix, and D the feedforward matrix.
The model feature matrix analysis module 404 specifically includes:
the model characteristic matrix analysis unit is used for extracting a characteristic matrix A in the small interference model for analysis, solving a characteristic value of the characteristic matrix A to obtain characteristic roots, wherein each group of characteristic roots represents a mode, and each mode corresponds to a state parameter; the state quantities comprise the state variables.
The invention provides an oscillation mode identification method and system for sending wind power out of a system through MMC-HVDC, which provide a basis for a parameter adjustment method for system stability analysis and oscillation inhibition by classifying system oscillation modes, including an internal oscillation mode of a wind power plant, an oscillation mode of a wind power plant control system, an internal oscillation mode of an MMC, an oscillation mode of an MMC control system and an interaction oscillation mode.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (10)

1. An oscillation mode identification method for wind power through an MMC-HVDC sending system is characterized by comprising the following steps:
establishing a mathematical model of each part of a wind power sending system through MMC-HVDC;
connecting the wind power through mathematical models of all parts of the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind power through the MMC-HVDC sending-out system;
carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out of the MMC-HVDC system;
extracting a characteristic matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through MMC-HVDC;
normalizing each state parameter, and locking a main participation link under a corresponding mode according to the state parameter with the highest correlation degree;
positioning an oscillation source and determining an oscillation mode according to the main participation link; the oscillation modes comprise a wind power plant internal oscillation mode, a wind power plant control system oscillation mode, an MMC internal oscillation mode, an MMC control system oscillation mode and an interaction oscillation mode.
2. The method according to claim 1, wherein the establishing of the mathematical model of each part of the wind power output system through MMC-HVDC specifically comprises:
dividing the wind power into a wind power generator set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system, respectively modeling each part independently, and establishing a mathematical model of each part of the wind power sending-out system through the MMC-HVDC
Figure FDA0003462425070000011
Wherein X is the state variable of the DC power transmission system, U is the input variable,
Figure FDA0003462425070000012
is the first derivative of the state variable, f () is the first derivative of the state variable
Figure FDA0003462425070000013
A functional relationship with the state variable X and the input variable U; y is the output variable and Y () is the functional relationship between the output variable Y and the state variable X and the input variable U.
3. The method according to claim 2, wherein the connecting the mathematical models of the parts of the MMC-HVDC delivery system of the wind power to obtain the integrated mathematical model of the parts of the MMC-HVDC delivery system of the wind power comprises:
and converting the electrical quantities of the wind turbine generator and the MMC converter station to the same coordinate system by adopting an interface conversion formula, and connecting the mathematical models of the wind power transmission system through MMC-HVDC transmission systems to obtain an integrated mathematical model of the wind power transmission system through MMC-HVDC transmission systems.
4. The method according to claim 3, wherein the linearizing the integrated mathematical model to obtain a small interference model of the wind power transmitted through the MMC-HVDC transmission system specifically comprises:
carrying out linearization treatment on the integrated mathematical model to obtain a small interference model of the wind power through an MMC-HVDC sending-out system
Figure FDA0003462425070000021
Where A is the feature matrix, B is the input matrix, C is the output matrix, and D is the feedforward matrix.
5. The method according to claim 4, wherein the extracting of the feature matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted from the MMC-HVDC transmission system specifically comprises:
extracting a characteristic matrix A in the small interference model for analysis, solving a characteristic value of the characteristic matrix A to obtain characteristic roots, wherein each group of characteristic roots represents a mode, and each mode corresponds to a state parameter; the state quantities comprise the state variables.
6. An oscillation mode identification system for wind power through an MMC-HVDC transmission system is characterized by comprising:
the mathematical model establishing module is used for establishing a mathematical model of each part of the wind power output system through the MMC-HVDC;
the mathematical model connecting module is used for connecting the mathematical models of all parts of the wind power output system through the MMC-HVDC to obtain an integrated mathematical model of the wind power output system through the MMC-HVDC;
the model linearization processing module is used for carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out of the MMC-HVDC system;
the model characteristic matrix analysis module is used for extracting a characteristic matrix in the small interference model for analysis to obtain all the modes and corresponding state parameters of the wind power transmitted out of the system through the MMC-HVDC;
the state parameter normalization processing module is used for performing normalization processing on each state parameter and locking main participation links under corresponding modes according to the state parameter with the highest correlation degree;
the oscillation positioning and oscillation mode identification module is used for positioning an oscillation source and determining an oscillation mode according to the main participation link; the oscillation modes comprise a wind power plant internal oscillation mode, a wind power plant control system oscillation mode, an MMC internal oscillation mode, an MMC control system oscillation mode and an interaction oscillation mode.
7. The system of claim 6, wherein the mathematical model building module specifically comprises:
a mathematical model establishing unit for dividing the wind power into a wind power set, an MMC converter station and a connecting circuit through an MMC-HVDC sending-out system, respectively modeling each part separately, and establishing a mathematical model of each part of the wind power sending-out system through the MMC-HVDC
Figure FDA0003462425070000031
Wherein X is the state variable of the DC power transmission system, U is the input variable,
Figure FDA0003462425070000032
is the first derivative of the state variable, f () is the first derivative of the state variable
Figure FDA0003462425070000033
A functional relationship with the state variable X and the input variable U; y is the output variable and Y () is the functional relationship between the output variable Y and the state variable X and the input variable U.
8. The system according to claim 7, wherein the mathematical model connection module specifically comprises:
and the mathematical model connecting unit is used for converting the electrical quantities of the wind turbine generator and the MMC converter station to the same coordinate system by adopting an interface conversion formula, and connecting the mathematical models of the wind turbine generator and the MMC-HVDC sending-out system to obtain an integrated mathematical model of the wind turbine generator and the MMC-HVDC sending-out system.
9. The system of claim 8, wherein the model linearization processing module specifically comprises:
the model linearization processing unit is used for carrying out linearization processing on the integrated mathematical model to obtain a small interference model of the wind power transmitted out from the MMC-HVDC system
Figure FDA0003462425070000034
Where A is the feature matrix, B is the input matrix, C is the output matrix, and D is the feedforward matrix.
10. The system of claim 9, wherein the model feature matrix analysis module specifically comprises:
the model characteristic matrix analysis unit is used for extracting a characteristic matrix A in the small interference model for analysis, solving a characteristic value of the characteristic matrix A to obtain characteristic roots, wherein each group of characteristic roots represents a mode, and each mode corresponds to a state parameter; the state quantities comprise the state variables.
CN202210020688.4A 2022-01-10 2022-01-10 Oscillation mode identification method for wind power transmission system through MMC-HVDC Pending CN114336685A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210020688.4A CN114336685A (en) 2022-01-10 2022-01-10 Oscillation mode identification method for wind power transmission system through MMC-HVDC

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210020688.4A CN114336685A (en) 2022-01-10 2022-01-10 Oscillation mode identification method for wind power transmission system through MMC-HVDC

Publications (1)

Publication Number Publication Date
CN114336685A true CN114336685A (en) 2022-04-12

Family

ID=81026089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210020688.4A Pending CN114336685A (en) 2022-01-10 2022-01-10 Oscillation mode identification method for wind power transmission system through MMC-HVDC

Country Status (1)

Country Link
CN (1) CN114336685A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110707728A (en) * 2019-09-09 2020-01-17 上海电力大学 Subsynchronous oscillation suppression method based on subsynchronous oscillation short-circuit ratio index
CN111049178A (en) * 2019-11-19 2020-04-21 国网甘肃省电力公司 Method for analyzing stability control of direct-drive permanent magnet wind turbine generator through VSC-HVDC grid connection
CN113270874A (en) * 2021-05-19 2021-08-17 国网内蒙古东部电力有限公司 Stability judging device for direct-drive wind power plant through LCC-HVDC sending-out system
CN113612237A (en) * 2021-07-16 2021-11-05 国网江苏省电力有限公司电力科学研究院 Method for positioning resonance-induced subsynchronous oscillation source in offshore wind farm

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110707728A (en) * 2019-09-09 2020-01-17 上海电力大学 Subsynchronous oscillation suppression method based on subsynchronous oscillation short-circuit ratio index
CN111049178A (en) * 2019-11-19 2020-04-21 国网甘肃省电力公司 Method for analyzing stability control of direct-drive permanent magnet wind turbine generator through VSC-HVDC grid connection
CN113270874A (en) * 2021-05-19 2021-08-17 国网内蒙古东部电力有限公司 Stability judging device for direct-drive wind power plant through LCC-HVDC sending-out system
CN113612237A (en) * 2021-07-16 2021-11-05 国网江苏省电力有限公司电力科学研究院 Method for positioning resonance-induced subsynchronous oscillation source in offshore wind farm

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘凯,等: "基于零序环流控制器的海上风电MMC-HVDC并网系统小干扰稳定分析与优化控制", 中国电机工程学报, vol. 41, no. 12, pages 4069 - 4080 *
曾禹村: "《信号与系统》", 30 September 2018, 北京理工大学出版社, pages: 349 - 358 *
王一凡,等: "混合型风电场经MMC-HVDC送出系统的振荡模式分析", 电力系统保护与控制, vol. 48, no. 09, pages 18 - 26 *
陈剑,等: "基于对抗式迁移学习的含柔性高压直流输电的风电系统次同步振荡源定位", 电工技术学报, vol. 36, no. 22, pages 4704 - 4715 *

Similar Documents

Publication Publication Date Title
CN109449958B (en) A kind of double-fed blower grid-connected system method for analyzing stability
CN110601268B (en) Doubly-fed fan grid-connected port output impedance modeling and stability analysis method
CN110994668A (en) Stability analysis method based on output impedance model of doubly-fed wind power plant grid-connected system
CN109586337B (en) Frequency domain modeling-based VSC grid-connected system subsynchronous oscillation risk assessment method
CN111697618B (en) Impedance modeling and stability analysis method for wide frequency domain full model of large-scale wind power plant
CN111327075B (en) Two-stage type photovoltaic grid-connected system low voltage ride through dynamic evaluation method and system
CN114597978A (en) PMSG grid-connected system stability determination method and system
CN112787326A (en) Doubly-fed wind turbine dynamic inertia evaluation method and system considering controller parameter influence
CN103117562A (en) Control method of high-voltage cascade energy feedback converter power module
CN113378347A (en) Wind turbine generator frequency domain impedance modeling method based on modularized multiport
Yu et al. Small-signal modeling of wind farm with direct-drive PMSG using the component connection method
CN114336685A (en) Oscillation mode identification method for wind power transmission system through MMC-HVDC
CN111769597A (en) Dimension reduction modeling analysis method of doubly-fed wind generator
CN107240922A (en) Marine Fractional Frequency Power Transmission System control method and device
Abdelwahed et al. Comparative study between fast terminal and second order sliding mode controls applied to a wind energy conversion system
Wu et al. Eliminating frequency coupling of DFIG system using a complex vector PLL
CN112952901B (en) Distributed stability analysis method for multi-fan grid-connected system
CN109657380A (en) A kind of double-fed fan motor field Dynamic Equivalence based on Extended Kalman filter
CN115632395A (en) Oscillation suppression method for direct current receiving end feed-in station and near-field wind power plant system
CN114792984A (en) Quick positioning method for subsynchronous/supersynchronous oscillation source
CN115021246A (en) Doubly-fed fan grid-connected system stability analysis method based on Gehr circle theorem
CN114465280A (en) Dynamic equivalent modeling method for new energy grid-connected system
CN112994113A (en) Doubly-fed wind generator sequence impedance modeling method and system
CN108667008B (en) Universal small signal modeling method for subsynchronous oscillation analysis of power system
Gao et al. Improved extended kalman filter based dynamic equivalent method of DFIG wind farm cluster

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
TA01 Transfer of patent application right

Effective date of registration: 20240221

Address after: 102206 No. 2 Huilongguan Beinong Road, Changping District, Beijing

Applicant after: NORTH CHINA ELECTRIC POWER University

Country or region after: China

Applicant after: HUANENG RUDONG BAXIANJIAO OFFSHORE WIND POWER CO.,LTD.

Applicant after: HUANENG GROUP TECHNOLOGY INNOVATION CENTER Co.,Ltd.

Applicant after: HUANENG CLEAN ENERGY Research Institute

Address before: 102206 No. 2 Nong Road, Huilongguan, Beijing, Changping District

Applicant before: NORTH CHINA ELECTRIC POWER University

Country or region before: China

Applicant before: HUANENG GROUP TECHNOLOGY INNOVATION CENTER Co.,Ltd.

Applicant before: HUANENG CLEAN ENERGY Research Institute

TA01 Transfer of patent application right