CN114357785A - Evaluation parameter determination method for stability analysis applicability of wind power plant equivalent model - Google Patents
Evaluation parameter determination method for stability analysis applicability of wind power plant equivalent model Download PDFInfo
- Publication number
- CN114357785A CN114357785A CN202210020686.5A CN202210020686A CN114357785A CN 114357785 A CN114357785 A CN 114357785A CN 202210020686 A CN202210020686 A CN 202210020686A CN 114357785 A CN114357785 A CN 114357785A
- Authority
- CN
- China
- Prior art keywords
- wind
- power plant
- wind power
- equivalent
- determining
- 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
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012545 processing Methods 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 17
- 238000006116 polymerization reaction Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 abstract description 8
- 230000005540 biological transmission Effects 0.000 abstract description 3
- 238000011158 quantitative evaluation Methods 0.000 abstract description 2
- 230000000875 corresponding effect Effects 0.000 description 30
- 238000011160 research Methods 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010949 in-process test method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Landscapes
- Wind Motors (AREA)
Abstract
The invention relates to an evaluation parameter determination method for stability analysis applicability of an equivalent model of a wind power plant, belonging to the technical field of power transmission and distribution, determining the wind speed of each wind turbine in the original wind power plant and the wind speed of each equivalent wind turbine in the equivalent wind power plant according to the wind speed of the region where the wind power plant is located, further determining the characteristic values of the wind power plant detailed model and the wind power plant equivalent model according to the wind speed, calculating the average Euclidean distance between the wind power plant detailed model and the wind power plant equivalent model and the corresponding characteristic values of the wind power plant detailed model and the wind power plant equivalent model on the key mode to be analyzed, the method comprises the steps of determining mathematical expectation and variance of the average Euclidean distances according to the average Euclidean distances corresponding to the wind speeds of the areas where the wind power plants are located, and enabling the mathematical expectation and variance to be evaluation parameters of stability analysis applicability of the equivalent model of the wind power plants to be evaluated, so that whether the equivalent model of the wind power plants is suitable for quantitative evaluation of oscillation stability analysis or not is achieved.
Description
Technical Field
The invention relates to the technical field of power transmission and distribution, in particular to an evaluation parameter determination method for stability analysis applicability of an equivalent model of a wind power plant.
Background
In recent years, with the rapid development of offshore wind power, a large number of offshore wind power plants adopt alternating current or flexible direct current power transmission grid connection, the problem of oscillation instability generated by grid connection of wind power plants also appears successively, sub-synchronous oscillation of about 30Hz occurs in the wind power plants in the flexible direct current grid connection engineering, and resonance accidents occur in the grid connection engineering of Germany north sea wind power plants in 2014, so that a large number of devices are damaged, and the system is stopped for half a year. These problems have attracted the attention of a large number of experts and scholars, and scholars at home and abroad have conducted a large amount of research on the problem of grid connection stability of wind power plants. The small-signal modeling of the wind power plant is the basis and premise for stability analysis of the wind power grid-connected system, compared with an electromagnetic transient model of the wind power plant, a complete wind power plant state space model is more complex, and detailed modeling is not practical, so that equivalence processing of the wind power plant becomes indispensable, and the wind power plant equivalence model is widely applied to wind power grid-connected stability analysis.
In the research of the grid-connected stability of the wind power plant, a single-machine equivalent model is adopted in most documents, namely the whole wind power plant is equivalent to one equivalent unit for research, and the whole wind power plant is equivalent to two fans in some documents. However, the influence of the equivalent model of the wind power plant on the stability of the grid-connected system of the wind power plant is not researched in the documents.
The existing article for researching the equivalent model of the wind power plant takes the transient response curves of active power, reactive power, voltage and current at the grid-connected position of the wind power plant as evaluation indexes, shows that the complete fitting of the transient response curves is pursued on the premise of continuously increasing modeling complexity, but the change of the equivalent wind power plant on the system stability is not concerned, so that the problems that errors are possible to exist when the existing equivalent model is used for researching the stability of the grid-connected system of the wind power plant, how to quantify the errors which exist and how to evaluate whether the equivalent model of the wind power plant is suitable for oscillation stability analysis are needed to be solved.
Disclosure of Invention
The invention aims to provide an evaluation parameter determination method for stability analysis applicability of a wind power plant equivalent model, so as to quantitatively evaluate whether the wind power plant equivalent model is suitable for oscillation stability analysis.
In order to achieve the purpose, the invention provides the following scheme:
an evaluation parameter determination method for stability analysis applicability of wind power plant equivalent models comprises the following steps:
determining the wind speed of each wind turbine in the original wind power plant based on the wake effect according to the wind speed of the region where the wind power plant is located;
determining the wind speed of each equivalent wind turbine in an equivalent wind power plant according to the wind speed of each wind turbine in the original wind power plant;
carrying out load flow calculation on an original wind power plant system according to the wind speed of each wind turbine in the original wind power plant to obtain a state variable initial value of the original wind power plant system, and carrying out load flow calculation on an equivalent wind power plant system according to the wind speed of each equivalent wind turbine in the equivalent wind power plant to obtain a state variable initial value of the equivalent wind power plant system;
determining a characteristic value of the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is brought in, and determining a characteristic value of the equivalent model of the wind power plant to be evaluated after the initial value of the state variable of the equivalent wind power plant system is brought in;
determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed; the corresponding characteristic value is a strong correlation variable of which one characteristic value of the wind power plant detailed model is similar to one characteristic value of the wind power plant equivalent model;
presetting the wind speeds of areas where a plurality of wind power plants are located, repeating the steps, and obtaining the average Euclidean distance corresponding to the wind speed of the area where each wind power plant is located;
determining mathematical expectation and variance of the average Euclidean distance according to the average Euclidean distance corresponding to the wind speed of the region where the wind power plants are located; the mathematical expectation and the variance are evaluation parameters of the stability analysis applicability of the wind power plant equivalent model to be evaluated.
Optionally, the determining, according to the wind speed of the area where the wind farm is located, the wind speed of each wind turbine in the original wind farm based on the wake effect specifically includes:
considering the wake effect, the wind speed of a wind turbine facing the wind direction in the original wind power plant is equal to the wind speed of the area where the wind power plant is located; the wind speed of the area where the wind power plant is located is between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator;
using formulas
Determining the wind speeds of wind turbines except the wind turbine facing the wind direction in the original wind power plant; wherein, VnIs the wind speed V of the nth wind turbine generator set except the wind turbine generator set facing the wind direction in the original wind power plantn-1The wind speed of the wind turbine generator set at the wind inlet of the nth wind turbine generator set.
Optionally, the determining the wind speed of each equivalent wind turbine in the equivalent wind farm according to the wind speed of each wind turbine in the original wind farm specifically includes:
according to the wind speed of each wind turbine in the original wind power plant, a weighted polymerization method is adopted, and a formula is utilized
Determining the wind speed of each equivalent wind turbine generator in the equivalent wind power plant;
wherein,
for equating m wind turbines of an original wind power plant to the wind speed of the ith equivalent wind turbine, viThe wind speed of the ith wind turbine generator set in the m wind turbine generator sets of the original wind power plant is obtained.
Optionally, the determining a feature value of the wind farm detailed model after the initial value of the state variable of the original wind farm system is introduced, and determining a feature value of the wind farm equivalent model to be evaluated after the initial value of the state variable of the equivalent wind farm system is introduced specifically include:
carrying out linearization processing on the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is introduced, and carrying out linearization processing on the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is introduced;
extracting a first state matrix from the wind power plant detailed model after linearization processing, and extracting a second state matrix from the wind power plant equivalent model after linearization processing;
and determining a characteristic value of the wind power plant detailed model according to the first state matrix, and determining a characteristic value of the wind power plant equivalent model according to the second state matrix.
Optionally, the determining an average euclidean distance between the wind farm detailed model and the wind farm equivalent model and corresponding feature values on the key modality to be analyzed specifically includes:
using formulas
Determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed;
wherein,
is the mean Euclidean distance, λjFor equivalent model of wind power plant in-process analysis gateJ characteristic value on key mode, lambdakFor the kth characteristic value, lambda, of the wind power plant detailed model on the key mode to be analyzedjAnd λkN is the number of feature roots, Re () represents the real part and Im () represents the imaginary part for the corresponding feature root.
An evaluation parameter determination system for stability analysis applicability of wind power plant equivalent models, the system comprising:
the first wind speed determining module is used for determining the wind speed of each wind turbine in the original wind power plant based on the wake effect according to the wind speed of the area where the wind power plant is located;
the second wind speed determination module is used for determining the wind speed of each equivalent wind turbine in the equivalent wind power plant according to the wind speed of each wind turbine in the original wind power plant;
the state variable initial value calculation module is used for carrying out load flow calculation on the original wind power plant system according to the wind speed of each wind turbine in the original wind power plant to obtain a state variable initial value of the original wind power plant system, and carrying out load flow calculation on the equivalent wind power plant system according to the wind speed of each equivalent wind turbine in the equivalent wind power plant to obtain a state variable initial value of the equivalent wind power plant system;
the characteristic value determination module is used for determining the characteristic value of the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is brought in, and determining the characteristic value of the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is brought in;
the average Euclidean distance determining module is used for determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed; the corresponding characteristic value is a strong correlation variable of which one characteristic value of the wind power plant detailed model is similar to one characteristic value of the wind power plant equivalent model;
the circulation module is used for presetting the wind speeds of the areas where the wind power plants are located, repeating the steps and obtaining the average Euclidean distance corresponding to the wind speed of the area where each wind power plant is located;
the evaluation parameter determining module is used for determining mathematical expectation and variance of the average Euclidean distance according to the average Euclidean distance corresponding to the wind speed of the area where the wind power plants are located; the mathematical expectation and the variance are evaluation parameters of the stability analysis applicability of the wind power plant equivalent model to be evaluated.
Optionally, the first wind speed determination module specifically includes:
the third wind speed determination submodule is used for considering the wake effect, and the wind speed of a wind turbine generator facing the wind direction in the original wind power plant is equal to the wind speed of the area where the wind power plant is located; the wind speed of the area where the wind power plant is located is between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator;
a fourth wind speed determination submodule for utilizing a formula
Determining the wind speeds of wind turbines except the wind turbine facing the wind direction in the original wind power plant; wherein, VnIs the wind speed V of the nth wind turbine generator set except the wind turbine generator set facing the wind direction in the original wind power plantn-1The wind speed of the wind turbine generator set at the wind inlet of the nth wind turbine generator set.
Optionally, the second wind speed determination module specifically includes:
a second wind speed determination submodule for adopting a weighted polymerization method and a formula according to the wind speed of each wind turbine in the original wind power plant
Determining the wind speed of each equivalent wind turbine generator in the equivalent wind power plant;
wherein,
for equating m wind turbines of an original wind power plant to the wind speed of the ith equivalent wind turbine, viThe wind speed of the ith wind turbine generator set in the m wind turbine generator sets of the original wind power plant is obtained.
Optionally, the characteristic value determining module specifically includes:
the linearization processing submodule is used for carrying out linearization processing on the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is introduced, and carrying out linearization processing on the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is introduced;
the state matrix extraction submodule is used for extracting a first state matrix from the wind power plant detailed model after linearization processing and extracting a second state matrix from the wind power plant equivalent model after linearization processing;
and the characteristic value determination submodule is used for determining the characteristic value of the wind power plant detailed model according to the first state matrix and determining the characteristic value of the wind power plant equivalent model according to the second state matrix.
Optionally, the average euclidean distance determining module specifically includes:
an average Euclidean distance determining submodule for utilizing the formula
Determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed;
wherein,
is the mean Euclidean distance, λjThe j characteristic value, lambda, of the wind power plant equivalent model on the key mode to be analyzedkFor the kth characteristic value, lambda, of the wind power plant detailed model on the key mode to be analyzedjAnd λkN is the number of feature roots, Re () represents the real part and Im () represents the imaginary part for the corresponding feature root.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses an evaluation parameter determination method for stability analysis applicability of a wind power plant equivalent model, which comprises the steps of determining the wind speed of each wind power unit in an original wind power plant and the wind speed of each equivalent wind power unit in the equivalent wind power plant according to the wind speed of a region where the wind power plant is located, further determining the characteristic values of a wind power plant detailed model and the wind power plant equivalent model according to the wind speeds, calculating the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model on a key mode to be analyzed, determining the mathematical expectation and variance of the average Euclidean distance according to the average Euclidean distances corresponding to the wind speeds of the regions where the wind power plants are located, and taking the mathematical expectation and variance as evaluation parameters of stability analysis applicability of the wind power plant equivalent model to be evaluated to realize quantitative evaluation of whether the wind power plant equivalent model is suitable for oscillation 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 evaluation parameter determination method for stability analysis applicability of an equivalent model of a wind power plant provided by the invention;
FIG. 2 is a schematic diagram of an evaluation parameter determination method for stability analysis applicability of a wind power plant equivalent model provided by the 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 evaluation parameter determination method for stability analysis applicability of a wind power plant equivalent model, so as to quantitatively evaluate whether the wind power plant equivalent model is suitable for oscillation stability analysis.
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.
The modeling process of the equivalent model of the wind power plant to be evaluated in the invention is as follows:
considering factors such as wind power plant topology, wind turbine generator parameters and structure, collection network parameters, alternating current system parameters and the like which a complete wind power plant model should include, thereby determining system state variables and input variables, and respectively establishing mathematical models of an equivalent wind power plant system and an original wind power plant system, as shown in the following formula (1):
wherein X is a state variable of the researched wind power plant system; u is an input variable; y is an output variable;
is a plurality of numbers.
The invention provides an evaluation parameter determination method for stability analysis applicability of a wind power plant equivalent model based on a built wind power plant equivalent model, as shown in a figure 1-2, the method comprises the following steps:
The method specifically comprises the following steps:
firstly, the wind speed of the area where the wind power plant is located is considered to be between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator (the phenomenon of wind abandon caused by overhigh wind speed and overlow wind speed is prevented)
Vcut_in＜Vw_farm＜Vcut_out
Wherein, Vw farmThe wind speed is the wind speed of the area where the wind power plant is located without attenuation; vcut_inThe cutting-in rotating speed; vcut outThe cut-out speed.
Then consider the simple wake effect (due to the effect of the wake, the captured wind speed for a downwind located wind turbine will be lower than an upwind located wind turbine, and operational experience with California wind farms in the United states indicates that the typical value of energy loss due to the wake is 10%).
Considering the wake effect, the wind speed of a wind turbine facing the wind direction in the original wind power plant is equal to the wind speed of the area where the wind power plant is located; the wind speed of the area where the wind power plant is located is between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator;
using formulas
Determining the wind speeds of wind turbines except the wind turbine facing the wind direction in the original wind power plant; wherein, VnIs the wind speed V of the nth wind turbine generator set except the wind turbine generator set facing the wind direction in the original wind power plantn-1The wind speed of the wind turbine generator set at the wind inlet of the nth wind turbine generator set.
And step 102, determining the wind speed of each equivalent wind turbine in the equivalent wind power plant according to the wind speed of each wind turbine in the original wind power plant.
The equivalent wind speed calculation methods of the equivalent models of different wind power plants may be different, and one calculation method is listed as follows: passing wind speed to the power of 3 v3Characteristic curve between output mechanical power of permanent magnet synchronous wind power generator, using v3Obtaining equivalent wind speed v by weighting polymerization methodeq. The triple root-mean-square wind speed weighting polymerization method does not take the influence of the wind energy utilization coefficient on equivalence into account, and the equivalence accuracy can be improved. The equivalent wind speed input by the equivalent fan is as follows:
wherein,
for equating m wind turbines of an original wind power plant to the wind speed of the ith equivalent wind turbine, viThe wind speed of the ith wind turbine generator set in the m wind turbine generator sets of the original wind power plant is obtained.
And 103, carrying out load flow calculation on the original wind power plant system according to the wind speed of each wind turbine in the original wind power plant to obtain a state variable initial value of the original wind power plant system, and carrying out load flow calculation on the equivalent wind power plant system according to the wind speed of each equivalent wind turbine in the equivalent wind power plant to obtain the state variable initial value of the equivalent wind power plant system.
And step 104, determining the characteristic value of the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is brought, and determining the characteristic value of the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is brought.
The method specifically comprises the following steps:
carrying out linearization processing on the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is introduced, and carrying out linearization processing on the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is introduced;
extracting a first state matrix from the wind power plant detailed model after linearization processing, and extracting a second state matrix from the wind power plant equivalent model after linearization processing;
and determining a characteristic value of the wind power plant detailed model according to the first state matrix, and determining a characteristic value of the wind power plant equivalent model according to the second state matrix.
The method comprises the following steps of carrying out linearization processing on an original wind power plant model and a wind power plant equivalent model, extracting respective state matrixes A, and analyzing the state matrixes A to obtain characteristic values of all key modes of an original system and an equivalent system, wherein the characteristic values are shown in the following formulas:
the key modality is the modality which the person using the model wants to analyze, and can be subsynchronous or high-frequency, and the evaluation method is also used for evaluating the equivalent model, so that the error of the equivalent model on the several modalities which the person using the equivalent model wants to research is not large.
And 105, determining the average Euclidean distance of the corresponding characteristic values of the wind power plant detailed model and the wind power plant equivalent model on the key mode to be analyzed.
The corresponding characteristic value is a characteristic value (or a group of characteristic values) of the wind power plant detailed modelThe characteristic root which is close in position and has the same change rule) has a similar strong correlation variable with a characteristic value in the equivalent model of the wind power plant. Namely, the strong correlation variable of the equivalent model feature root should be located in the same link in the equipment model as the strong correlation variable of the corresponding feature root in the original model. For example, the strongly correlated variable of a pair of complex characteristic roots of the equivalent model is direct current voltage udcState variable x corresponding to DC voltage control loop1If u is equal to the set in the detailed model, u is also equal to the set in the detailed modeldcAnd x1The one (or set) of feature roots that are strongly correlated.
The method specifically comprises the following steps:
using formulas
Determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed;
wherein,
is the mean Euclidean distance, λjThe j characteristic value, lambda, of the wind power plant equivalent model on the key mode to be analyzedkFor the kth characteristic value, lambda, of the wind power plant detailed model on the key mode to be analyzedjAnd λkN is the number of feature roots, Re () represents the real part and Im () represents the imaginary part for the corresponding feature root.
And taking the average Euclidean distance of the corresponding characteristic roots of the original model and the equivalent model as an index for measuring the accuracy of the equivalent model.
And 106, presetting the wind speeds of the areas where the wind power plants are located, and repeating the steps to obtain the average Euclidean distance corresponding to the wind speed of the area where each wind power plant is located.
The average Euclidean distance under different output conditions of the wind power plant means that the wind speed V of the wind power plant in the step (2) is changedw_farmThereby obtaining the fresh air at a speed
The wind speed V of the wind power plant is obtainedw_farmEach at a change between cut-in wind speed and cut-out wind speed
Further obtain
Mathematical expectation of
Sum variance
And further providing a reference for the wind power plant equivalent model to be suitable for stability analysis.
Through the steps, the mean value and the variance of errors of different equivalent models and the detailed model on the key mode when the total output of the wind power plant changes can be obtained, the larger the mean value and the variance are, the larger the error of the equivalent model on the key mode is, and the wind power plant equivalent model with the smaller mean value and the smaller variance has better practicability.
According to the method, the error of the equivalent model on stability analysis is quantified by solving the average Euclidean distance between the equivalent model of the wind power plant and the corresponding characteristic value of the detailed model, and a reference is provided for the application of the equivalent model of the wind power plant to oscillation stability analysis.
The invention also provides an evaluation parameter determination system for the stability analysis applicability of the wind power plant equivalent model, which comprises the following steps:
the first wind speed determining module is used for determining the wind speed of each wind turbine in the original wind power plant based on the wake effect according to the wind speed of the area where the wind power plant is located;
the second wind speed determination module is used for determining the wind speed of each equivalent wind turbine in the equivalent wind power plant according to the wind speed of each wind turbine in the original wind power plant;
the state variable initial value calculation module is used for carrying out load flow calculation on the original wind power plant system according to the wind speed of each wind turbine in the original wind power plant to obtain a state variable initial value of the original wind power plant system, and carrying out load flow calculation on the equivalent wind power plant system according to the wind speed of each equivalent wind turbine in the equivalent wind power plant to obtain a state variable initial value of the equivalent wind power plant system;
the characteristic value determination module is used for determining the characteristic value of the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is brought in, and determining the characteristic value of the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is brought in;
the average Euclidean distance determining module is used for determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed; the corresponding characteristic value is a strong correlation variable of which one characteristic value of the wind power plant detailed model is similar to one characteristic value of the wind power plant equivalent model;
the circulation module is used for presetting the wind speeds of the areas where the wind power plants are located, repeating the steps and obtaining the average Euclidean distance corresponding to the wind speed of the area where each wind power plant is located;
the evaluation parameter determining module is used for determining mathematical expectation and variance of the average Euclidean distance according to the average Euclidean distance corresponding to the wind speed of the area where the wind power plants are located; the mathematical expectation and the variance are evaluation parameters of the stability analysis applicability of the wind power plant equivalent model to be evaluated.
The first wind speed determination module specifically includes:
the third wind speed determination submodule is used for considering the wake effect, and the wind speed of a wind turbine generator facing the wind direction in the original wind power plant is equal to the wind speed of the area where the wind power plant is located; the wind speed of the area where the wind power plant is located is between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator;
a fourth wind speed determination submodule for utilizing a formula
Determining the wind speeds of wind turbines except the wind turbine facing the wind direction in the original wind power plant; wherein, VnIs the wind speed V of the nth wind turbine generator set except the wind turbine generator set facing the wind direction in the original wind power plantn-1The wind speed of the wind turbine generator set at the wind inlet of the nth wind turbine generator set.
The second wind speed determination module specifically includes:
a second wind speed determination submodule for adopting a weighted polymerization method and a formula according to the wind speed of each wind turbine in the original wind power plant
Determining the wind speed of each equivalent wind turbine generator in the equivalent wind power plant;
wherein,
for equating m wind turbines of an original wind power plant to the wind speed of the ith equivalent wind turbine, viThe wind speed of the ith wind turbine generator set in the m wind turbine generator sets of the original wind power plant is obtained.
The characteristic value determination module specifically comprises:
the linearization processing submodule is used for carrying out linearization processing on the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is introduced, and carrying out linearization processing on the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is introduced;
the state matrix extraction submodule is used for extracting a first state matrix from the wind power plant detailed model after linearization processing and extracting a second state matrix from the wind power plant equivalent model after linearization processing;
and the characteristic value determination submodule is used for determining the characteristic value of the wind power plant detailed model according to the first state matrix and determining the characteristic value of the wind power plant equivalent model according to the second state matrix.
The average euclidean distance determining module specifically includes:
a mean Euclidean distance determination submodule for utilizing the commonFormula (II)
Determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed;
wherein,
is the mean Euclidean distance, λjThe j characteristic value, lambda, of the wind power plant equivalent model on the key mode to be analyzedkFor the kth characteristic value, lambda, of the wind power plant detailed model on the key mode to be analyzedjAnd λkN is the number of feature roots, Re () represents the real part and Im () represents the imaginary part for the corresponding feature root.
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 evaluation parameter determination method for wind power plant equivalent model stability analysis applicability is characterized by comprising the following steps:
determining the wind speed of each wind turbine in the original wind power plant based on the wake effect according to the wind speed of the region where the wind power plant is located;
determining the wind speed of each equivalent wind turbine in an equivalent wind power plant according to the wind speed of each wind turbine in the original wind power plant;
carrying out load flow calculation on an original wind power plant system according to the wind speed of each wind turbine in the original wind power plant to obtain a state variable initial value of the original wind power plant system, and carrying out load flow calculation on an equivalent wind power plant system according to the wind speed of each equivalent wind turbine in the equivalent wind power plant to obtain a state variable initial value of the equivalent wind power plant system;
determining a characteristic value of the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is brought in, and determining a characteristic value of the equivalent model of the wind power plant to be evaluated after the initial value of the state variable of the equivalent wind power plant system is brought in;
determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed; the corresponding characteristic value is a strong correlation variable of which one characteristic value of the wind power plant detailed model is similar to one characteristic value of the wind power plant equivalent model;
presetting the wind speeds of areas where a plurality of wind power plants are located, repeating the steps, and obtaining the average Euclidean distance corresponding to the wind speed of the area where each wind power plant is located;
determining mathematical expectation and variance of the average Euclidean distance according to the average Euclidean distance corresponding to the wind speed of the region where the wind power plants are located; the mathematical expectation and the variance are evaluation parameters of the stability analysis applicability of the wind power plant equivalent model to be evaluated.
2. The method for determining the evaluation parameter of the wind power plant equivalent model stability analysis applicability according to claim 1, is characterized in that the determining of the wind speed of each wind turbine in the original wind power plant based on the wake effect according to the wind speed of the region where the wind power plant is located specifically comprises:
considering the wake effect, the wind speed of a wind turbine facing the wind direction in the original wind power plant is equal to the wind speed of the area where the wind power plant is located; the wind speed of the area where the wind power plant is located is between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator;
using formulas
Determining the wind speeds of wind turbines except the wind turbine facing the wind direction in the original wind power plant; wherein, VnIs the wind speed V of the nth wind turbine generator set except the wind turbine generator set facing the wind direction in the original wind power plantn-1The wind speed of the wind turbine generator set at the wind inlet of the nth wind turbine generator set.
3. The method for determining the evaluation parameter of the stability analysis applicability of the wind farm equivalent model according to claim 1, is characterized in that the determination of the wind speed of each equivalent wind turbine in the equivalent wind farm according to the wind speed of each wind turbine in the original wind farm specifically comprises:
according to the wind speed of each wind turbine in the original wind power plant, a weighted polymerization method is adopted, and a formula is utilized
Determining the wind speed of each equivalent wind turbine generator in the equivalent wind power plant;
wherein,
for equating m wind turbines of an original wind power plant to the wind speed of the ith equivalent wind turbine, viThe wind speed of the ith wind turbine generator set in the m wind turbine generator sets of the original wind power plant is obtained.
4. The method for determining the evaluation parameter of the wind power plant equivalent model stability analysis applicability according to claim 1, is characterized in that the determining of the feature value of the wind power plant detailed model after the initial state variable value of the original wind power plant system is introduced, and the determining of the feature value of the wind power plant equivalent model to be evaluated after the initial state variable value of the equivalent wind power plant system is introduced specifically comprises:
carrying out linearization processing on the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is introduced, and carrying out linearization processing on the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is introduced;
extracting a first state matrix from the wind power plant detailed model after linearization processing, and extracting a second state matrix from the wind power plant equivalent model after linearization processing;
and determining a characteristic value of the wind power plant detailed model according to the first state matrix, and determining a characteristic value of the wind power plant equivalent model according to the second state matrix.
5. The method for determining the evaluation parameter of the stability analysis applicability of the wind farm equivalent model according to claim 1, characterized in that the determining of the average Euclidean distance between the wind farm detailed model and the wind farm equivalent model in the key mode to be analyzed and corresponding to the feature value specifically comprises:
using formulas
Determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed;
wherein,
is the mean Euclidean distance, λjThe j characteristic value, lambda, of the wind power plant equivalent model on the key mode to be analyzedkFor the kth characteristic value, lambda, of the wind power plant detailed model on the key mode to be analyzedjAnd λkN is the number of feature roots, Re () represents the real part and Im () represents the imaginary part for the corresponding feature root.
6. An evaluation parameter determination system for stability analysis applicability of wind power plant equivalent models is characterized by comprising the following components:
the first wind speed determining module is used for determining the wind speed of each wind turbine in the original wind power plant based on the wake effect according to the wind speed of the area where the wind power plant is located;
the second wind speed determination module is used for determining the wind speed of each equivalent wind turbine in the equivalent wind power plant according to the wind speed of each wind turbine in the original wind power plant;
the state variable initial value calculation module is used for carrying out load flow calculation on the original wind power plant system according to the wind speed of each wind turbine in the original wind power plant to obtain a state variable initial value of the original wind power plant system, and carrying out load flow calculation on the equivalent wind power plant system according to the wind speed of each equivalent wind turbine in the equivalent wind power plant to obtain a state variable initial value of the equivalent wind power plant system;
the characteristic value determination module is used for determining the characteristic value of the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is brought in, and determining the characteristic value of the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is brought in;
the average Euclidean distance determining module is used for determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed; the corresponding characteristic value is a strong correlation variable of which one characteristic value of the wind power plant detailed model is similar to one characteristic value of the wind power plant equivalent model;
the circulation module is used for presetting the wind speeds of the areas where the wind power plants are located, repeating the steps and obtaining the average Euclidean distance corresponding to the wind speed of the area where each wind power plant is located;
the evaluation parameter determining module is used for determining mathematical expectation and variance of the average Euclidean distance according to the average Euclidean distance corresponding to the wind speed of the area where the wind power plants are located; the mathematical expectation and the variance are evaluation parameters of the stability analysis applicability of the wind power plant equivalent model to be evaluated.
7. The evaluation parameter determination system for the wind farm equivalent model stability analysis applicability according to claim 6, wherein the first wind speed determination module specifically comprises:
the third wind speed determination submodule is used for considering the wake effect, and the wind speed of a wind turbine generator facing the wind direction in the original wind power plant is equal to the wind speed of the area where the wind power plant is located; the wind speed of the area where the wind power plant is located is between the cut-in rotating speed and the cut-out rotating speed of the wind turbine generator;
fourth wind speed determinationStator module for utilizing formula
Determining the wind speeds of wind turbines except the wind turbine facing the wind direction in the original wind power plant; wherein, VnIs the wind speed V of the nth wind turbine generator set except the wind turbine generator set facing the wind direction in the original wind power plantn-1The wind speed of the wind turbine generator set at the wind inlet of the nth wind turbine generator set.
8. The evaluation parameter determination system for the wind farm equivalent model stability analysis applicability according to claim 6, wherein the second wind speed determination module specifically comprises:
a second wind speed determination submodule for adopting a weighted polymerization method and a formula according to the wind speed of each wind turbine in the original wind power plant
Determining the wind speed of each equivalent wind turbine generator in the equivalent wind power plant;
wherein,
for equating m wind turbines of an original wind power plant to the wind speed of the ith equivalent wind turbine, viThe wind speed of the ith wind turbine generator set in the m wind turbine generator sets of the original wind power plant is obtained.
9. The evaluation parameter determination system for the wind farm equivalent model stability analysis applicability according to claim 6, characterized in that the characteristic value determination module specifically comprises:
the linearization processing submodule is used for carrying out linearization processing on the wind power plant detailed model after the initial value of the state variable of the original wind power plant system is introduced, and carrying out linearization processing on the wind power plant equivalent model to be evaluated after the initial value of the state variable of the equivalent wind power plant system is introduced;
the state matrix extraction submodule is used for extracting a first state matrix from the wind power plant detailed model after linearization processing and extracting a second state matrix from the wind power plant equivalent model after linearization processing;
and the characteristic value determination submodule is used for determining the characteristic value of the wind power plant detailed model according to the first state matrix and determining the characteristic value of the wind power plant equivalent model according to the second state matrix.
10. The evaluation parameter determination system for the wind power plant equivalent model stability analysis applicability according to claim 6, characterized in that the average Euclidean distance determination module specifically comprises:
an average Euclidean distance determining submodule for utilizing the formula
Determining the average Euclidean distance of the wind power plant detailed model and the wind power plant equivalent model corresponding to the characteristic value on the key mode to be analyzed;
wherein,
is the mean Euclidean distance, λjThe j characteristic value, lambda, of the wind power plant equivalent model on the key mode to be analyzedkFor the kth characteristic value, lambda, of the wind power plant detailed model on the key mode to be analyzedjAnd λkN is the number of feature roots, Re () represents the real part and Im () represents the imaginary part for the corresponding feature root.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210020686.5A CN114357785A (en) | 2022-01-10 | 2022-01-10 | Evaluation parameter determination method for stability analysis applicability of wind power plant equivalent model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210020686.5A CN114357785A (en) | 2022-01-10 | 2022-01-10 | Evaluation parameter determination method for stability analysis applicability of wind power plant equivalent model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114357785A true CN114357785A (en) | 2022-04-15 |
Family
ID=81107037
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210020686.5A Pending CN114357785A (en) | 2022-01-10 | 2022-01-10 | Evaluation parameter determination method for stability analysis applicability of wind power plant equivalent model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114357785A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115750230A (en) * | 2022-11-24 | 2023-03-07 | 盛东如东海上风力发电有限责任公司 | Method and system for dividing digital simulation model of wind turbine generator |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016023527A1 (en) * | 2014-08-14 | 2016-02-18 | 国家电网公司 | Method based on anemometer tower measurement data for determining wind farm discarded wind power |
| CN111651939A (en) * | 2020-06-08 | 2020-09-11 | 华北电力大学 | Dynamic equivalent modeling method for permanent magnet wind power plant considering converter control parameter difference |
-
2022
- 2022-01-10 CN CN202210020686.5A patent/CN114357785A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016023527A1 (en) * | 2014-08-14 | 2016-02-18 | 国家电网公司 | Method based on anemometer tower measurement data for determining wind farm discarded wind power |
| CN111651939A (en) * | 2020-06-08 | 2020-09-11 | 华北电力大学 | Dynamic equivalent modeling method for permanent magnet wind power plant considering converter control parameter difference |
Non-Patent Citations (3)
| Title |
|---|
| ZETIAN ZHENG 等: "Evaluation method for equivalent medels of PMSG-based wind farms considering randomness", IEEE TRANSACTIONS ON SUSTAINABLE ENERGY, vol. 10, no. 2, 22 November 2018 (2018-11-22), pages 1565 - 1574 * |
| 安之 等: "考虑风电随机性的直驱风机风电场等值模型评价方法", 中国电机工程学报, vol. 38, no. 22, 28 September 2018 (2018-09-28), pages 6511 - 6520 * |
| 蔡文超 等: "基于RTDS风电场等值建模若干问题的解决方案", 电气技术, vol. 38, no. 06, 15 June 2017 (2017-06-15), pages 6511 - 6520 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115750230A (en) * | 2022-11-24 | 2023-03-07 | 盛东如东海上风力发电有限责任公司 | Method and system for dividing digital simulation model of wind turbine generator |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102012956B (en) | A kind of wind energy turbine set equivalence method considering the random fluctuation of wind energy turbine set input wind speed and direction | |
| CN103887815B (en) | Based on wind energy turbine set parameter identification and the Dynamic Equivalence of service data | |
| CN109063276B (en) | Wind power plant dynamic equivalent modeling method suitable for long-time domain wind speed fluctuation | |
| CN107947228B (en) | Stochastic stability analysis method for power system containing wind power based on Markov theory | |
| Zhang et al. | Control optimisation for pumped storage unit in micro‐grid with wind power penetration using improved grey wolf optimiser | |
| CN115017854A (en) | Method for calculating maximum allowable capacity of DG (distributed generation) of power distribution network based on multidimensional evaluation index system | |
| CN108899930A (en) | Wind-powered electricity generation station equivalent modeling method based on Principal Component Analysis Method and hierarchical clustering algorithm | |
| CN105762816A (en) | Method for inhibiting subsynchronous interaction | |
| CN114357785A (en) | Evaluation parameter determination method for stability analysis applicability of wind power plant equivalent model | |
| Zheng et al. | Transient stable region of power systems incorporating stochasticity of variable renewable energies and system contingencies | |
| CN106410862A (en) | Wind power plant single machine equivalent method based on active recovery slope correction | |
| CN111651939B (en) | Permanent magnet wind power plant dynamic equivalent modeling method considering control parameter difference of converter | |
| CN112039051A (en) | Real-time modeling method for accessing double-fed wind driven generator into substation bus load | |
| CN109802440B (en) | Offshore wind farm equivalence method, system and device based on wake effect factor | |
| CN117375068A (en) | Nonlinear affine-based wind farm voltage state calculation and optimization method and computer-readable medium | |
| CN115085260B (en) | Oscillation mode fidelity evaluation method of new energy station aggregation model | |
| Li et al. | Wind power forecasting based on time series and neural network | |
| CN109586300B (en) | Method and system for acquiring change interval of tide variable in wind power tide model | |
| Yan et al. | Transient modelling of doubly‐fed induction generator based wind turbine on full operation condition and rapid starting period based on low voltage ride‐through testing | |
| CN104268316A (en) | Off-line calculation and online matching based doubly-fed wind power station probability equivalent modeling method | |
| CN109657380A (en) | A kind of double-fed fan motor field Dynamic Equivalence based on Extended Kalman filter | |
| CN114357787A (en) | Equivalent modeling method and system for offshore wind plant | |
| CN114465280A (en) | Dynamic equivalent modeling method for new energy grid-connected system | |
| Fang et al. | Dynamic equivalence of wind farm considering operational condition of wind turbines | |
| CN104167735B (en) | A kind of non-mechanism equivalent modeling method of wind energy turbine set and device |
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 |