CN113935247A - Partition virtual inertia estimation method considering wind speed randomness and correlation - Google Patents

Partition virtual inertia estimation method considering wind speed randomness and correlation Download PDF

Info

Publication number
CN113935247A
CN113935247A CN202111302789.2A CN202111302789A CN113935247A CN 113935247 A CN113935247 A CN 113935247A CN 202111302789 A CN202111302789 A CN 202111302789A CN 113935247 A CN113935247 A CN 113935247A
Authority
CN
China
Prior art keywords
wind speed
wind
inertia
speed
correlation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111302789.2A
Other languages
Chinese (zh)
Other versions
CN113935247B (en
Inventor
马喜平
甄文喜
姚寅
沈渭程
王臣赟
陈柏旭
董开松
朱宏毅
赵霖
魏博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
State Grid Gansu Electric Power Co Ltd
Original Assignee
STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
State Grid Gansu Electric Power 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 STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE, State Grid Gansu Electric Power Co Ltd filed Critical STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
Priority to CN202111302789.2A priority Critical patent/CN113935247B/en
Publication of CN113935247A publication Critical patent/CN113935247A/en
Application granted granted Critical
Publication of CN113935247B publication Critical patent/CN113935247B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/06Wind turbines or wind farms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • 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
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Water Supply & Treatment (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Human Resources & Organizations (AREA)
  • Artificial Intelligence (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Wind Motors (AREA)

Abstract

The invention discloses a partition virtual inertia estimation method considering wind speed randomness and correlation, which comprises the following steps of: step 1: considering wake flow and time delay effects of wind speed, processing historical wind speed data by adopting a segmented aggregation approximation method, and establishing daily variation curves of the wind speed of each unit in the field; step 2: clustering the daily variation curve of the wind speed of each unit by adopting a dual-scale spectral clustering algorithm to realize the partition of the units in the field and form a plurality of clusters; and step 3: fitting the wind speed probability distribution of each cluster center unit by adopting a nonparametric kernel density estimation method, and constructing an optimal Copula function to analyze the correlation of the wind speed between each cluster center and a wind measuring tower; and 4, step 4: and (4) considering uncertainty between wind speed and inertia, and estimating virtual inertia storage of each cluster in the wind power plant based on the internal structure and wind speed data of the actual wind power plant.

Description

Partition virtual inertia estimation method considering wind speed randomness and correlation
Technical Field
The invention belongs to the technical field of wind power generation, and particularly relates to a partition virtual inertia estimation method considering wind speed randomness and correlation.
Background
With the large-scale grid connection of wind power generation, the problems caused by the intermittency and the fluctuation of wind energy are increasingly prominent. Because the wind turbine generator is connected with the power grid through the converter, the wind turbine generator is decoupled from the frequency of the power grid, the wind turbine generator cannot directly participate in system frequency modulation, and effective inertia support cannot be provided. The available inertia of the entire power system will continue to decrease and the frequency response capability will continue to decrease. Therefore, the wind power plant with the power of 30MW or more must have a primary frequency modulation function in the detailed rules for implementation of auxiliary service management of grid-connected power plants in central China, published by the national energy agency in 2020. The frequency response of the fan is mainly realized by the technologies of virtual inertia control, droop control, primary frequency modulation and the like. With the popularization and application of the technologies, the wind turbine generator set has equivalent inertia supporting capacity. However, an accurate estimation method of the available inertia of each unit in the wind power plant is still lacking. The accurate estimation of the available inertia can provide reference for the distribution and the scheduling of the auxiliary frequency modulation output of each unit after disturbance, and also can provide theoretical support for the frequency modulation control strategy setting of the units in the region according to the reserve of the available inertia in each region.
The main influencing factors of the available inertia of the wind farm are the wind speed and the operating state in the wind farm. The wind power generation sets in the large wind power plant are large in number, the distance between the wind power generation sets is close, and certain correlation exists between the wind speeds of the wind power generation sets and the running state of the wind power generation sets. Documents, namely a wind power plant wind speed dependent model constructed by rattan Copula and a scene generation method considering the relevant relation of the Copula output of the multi-wind power plant, adopt Copula functions to analyze the relevance of wind speeds in the wind power plant. In the literature, the multi-wind-field output scene generation method considering the time-space correlation considers various wind speed influence factors, establishes a conditional Copula joint distribution function for wind power prediction, and realizes interval prediction of wind power field output. The literature, "wind power correlation analysis based on a mixed Copula function" considers the space-time correlation of wind speed, and a wind speed space-time distribution correlation model is established based on a multivariate normal distribution function and a Copula function. According to the literature, tail characteristics among fan output are considered in a wind power plant virtual inertia multi-machine cooperative control strategy considering frequency modulation capability, and a wind power output correlation analysis method based on a hybrid Copula function is provided. Research shows that the mixed Copula function has obvious advantages in representing wind speed correlation, but when constructing the mixed multi-element Copula function, related parameters are difficult to estimate, the quality of the correlation representation is influenced by weight value selection, the intermediate calculation process is complex, and for the reason, when the Copula function is used for describing the correlation between variables, the Copula function is used for calculating the correlation between the variables. On the premise of ensuring no distortion, the intermediate calculation process can be greatly simplified by selecting a single Copula function.
The current research on wind power plant inertia response focuses on virtual inertia control optimization, relatively few researches on wind power plant available inertia evaluation are carried out, and the wind power plant available inertia evaluation is evaluation of the whole wind power plant inertia. A novel method for evaluating potential help of wind capacity to system frequency stability provides an evaluation method for frequency modulation capability of a plurality of wind power plants based on the principle that rotating kinetic energy of wind turbine generators is equal. According to the literature, the probability method for evaluating comprehensive droop and inertial response of the wind power plant evaluates the joint inertial response capability of the wind power plant at a specific average wind speed by processing a Gaussian distribution model of atmospheric turbulence in blocks. On the basis, the literature, "analysis modeling of wind power plant inertia and droop response for short-term frequency regulation of a power system" comprehensively evaluates the joint frequency modulation effect of virtual inertia control and droop control. However, the wind power plants are taken as a whole by the above evaluation methods, and the distribution difference of the available inertia in the wind power plants is ignored, so that a large error may exist between the evaluation result and the actual value, and the problem of the distribution of the available inertia in the wind power plants cannot be solved.
The estimation of the available inertia of the wind farm typically uses the wind speed of a particular set of locations or the average wind speed of the wind farm to estimate the level of rotational inertia of the wind farm as a whole. An anti-intersection to copulas document proposes a load shedding model based on VSWT speed control, and deduces An inertia time constant from a single variable speed constant frequency wind generating set. However, due to the problem of computational efficiency, the model cannot accurately estimate the available inertia of the large wind farm.
Disclosure of Invention
The invention provides a partition virtual inertia estimation method considering wind speed randomness and correlation, and takes the wind speed correlation of each unit in a wind power plant, the position distribution of each unit, the running state and other factors into consideration.
The technical scheme adopted by the invention is as follows:
the partition virtual inertia estimation method considering the randomness and the correlation of the wind speed comprises the following steps of:
step 1: considering wake flow and time delay effects of wind speed, processing historical wind speed data by adopting a segmented aggregation approximation method, and establishing daily variation curves of the wind speed of each unit in the field;
step 2: clustering the daily variation curve of the wind speed of each unit by adopting a dual-scale spectral clustering algorithm to realize the partition of the units in the field and form a plurality of clusters;
and step 3: fitting the wind speed probability distribution of each cluster center unit by adopting a nonparametric kernel density estimation method, and constructing an optimal Copula function to analyze the correlation of the wind speed between each cluster center and a wind measuring tower;
and 4, step 4: and (4) considering uncertainty between wind speed and inertia, and estimating virtual inertia storage of each cluster in the wind power plant based on the internal structure and wind speed data of the actual wind power plant.
Preferably, in step 2, the partitioning process of the wind farm is implemented as follows:
step 201: distance-based similarity metric
Set X (X) with M wind speed data1,x2,…,xM) Wherein each wind speed curve xiHas a data length of L; the Euclidean distance is the real distance between two data points in the L-dimensional space, and the wind speed curve xiAnd xjThe euclidean distance between them is defined as:
Figure BDA0003338986730000031
step 202: similarity measure based on morphological features
The morphological characteristics represent the fluctuation degree of the curve, and the information entropy is used as a measurement index of the fluctuation degree of the wind speed curve; setting a wind speed curve x in a certain time periodiThere are n possible values of x1,x2,…,xn(ii) a Probability of each possible value being p1,p2,…,pn(ii) a The wind speed curve x for that time periodiEntropy value H ofnIs defined as:
Figure BDA0003338986730000041
the size of the entropy reflects the fluctuation degree of the curve, and the larger the entropy is, the larger the fluctuation degree of the curve is; similarly, the closer the information entropy values of the two wind speed curves are, the higher the fluctuation similarity of the two curves is;
the wind speed similarity matrix of each fan is constructed by extracting similarity measurement features based on distance and form, and daily variation curves of wind speeds close to each other in distance and similar in fluctuation degree are gathered into a class by utilizing a spectral clustering algorithm, so that the partition of the wind power plant is realized.
Preferably, in step 3, after the in-field unit clustering and dimensionality reduction is performed in step 2, the whole wind power plant is divided into N unit areas, a wind speed sequence joint distribution probability histogram between each unit area and a wind measuring tower is drawn, tail characteristics and symmetry characteristics of the wind speed sequence joint distribution probability histogram are analyzed, and an optimal Copula function is selected according to a method with the minimum Euclidean distance from an empirical Copula function.
Preferably, the specific process of constructing the optimal Copula function to analyze the correlation between the wind speed of each cluster center and the anemometer tower is as follows:
step 301: first, an empirical Copula function, a random variable X, Y, and empirical distribution functions of sum F, respectively, are introducedn(x) And Hn(x) The empirical Copula function is defined as:
Figure BDA0003338986730000042
in which u, v is E [0,1 ]];Fn(xi) When the content is less than or equal to u,
Figure BDA0003338986730000043
otherwise, the value is 0;
obtaining relevant unknown parameters of each theoretical Copula function by using a maximum likelihood estimation method; according to the formula (1), the euclidean distance squared of the theoretical Copula function and the empirical Copula function is obtained by combining the defined formula of the euclidean distance and the empirical Copula function as follows:
Figure BDA0003338986730000044
wherein P is a theoretical Copula function type,
Figure BDA0003338986730000051
reflects the squared Euclidean distance between a theoretical Copula function and an empirical Copula function,
Figure BDA0003338986730000052
the Copula function corresponding to the minimum value is an optimal Copula function;
after the optimal Copula function is selected, the binary joint probability density function between the cluster i and the anemometer tower m is as follows:
fi(vi,vm)=Ci(uvi,vvm)f(vi)f(vm) (5)
when inertia estimation is carried out, detailed wind speed data measured by the anemometer tower are known, an optimal Copula function is established to describe the wind speed correlation between the anemometer tower and each cluster, and therefore the instantaneous wind speed probability density function of each cluster is obtained as follows:
Figure BDA0003338986730000053
preferably, in step 4, the available inertia estimation process between the partitions based on the wind speed is as follows:
according to the load shedding control mode of the wind power plant, the wind power plant is divided into a low wind speed area, a medium wind speed area and a high wind speed area;
in the low wind speed region, according to the related theory of aerodynamics, when the wind turbine runs in the MPPT operation mode, the mechanical power captured by the wind turbine can be represented by the following formula:
Figure BDA0003338986730000054
where ρ is the air density, CP maxThe wind energy capture coefficient is shown as lambda is the blade tip speed ratio, beta is the pitch angle, R is the fan blade radius, and V is the wind speed; wherein, the wind energy capture coefficient CP maxAssociated with the tip speed ratio λ and pitch angle β:
Figure BDA0003338986730000055
Figure BDA0003338986730000056
in the above formula, omega is the rotating speed of the generator; p is the number of pole pairs of the generator; g is the number of pole pairs of the gear box;
when the load shedding ratio is d%, the power value is shown in the following formula:
Figure BDA0003338986730000061
combining the two formulas to obtain:
CP-de=(1-d%)CP max (11)
the tip speed ratio when the load shedding rate is d% can be obtained according to a bisection method, and the relation between the fan rotating speed and the wind speed when the load shedding rate is d% can be obtained through the definition of the tip speed ratio:
Figure BDA0003338986730000062
when the fan operates in a medium wind speed and constant rotating speed area, the wind speed basically reaches the upper limit, overspeed control cannot continue to act, the rotating speed of a rotor of the fan cannot obviously change along with the wind speed, but the power captured by the fan is continuously improved through variable pitch control, the rotating speed equation of the part can be approximately linearized, and the wind speed and rotating speed value (v) of two end points of a curve is obtainedn,wmax)、(vws,ws) Constructing a straight line, and obtaining a relation of the rotating speed and the wind speed as follows:
Figure BDA0003338986730000063
when the wind speed of the fan reaches the high wind speed constant power region, the rotating speed of the fan reaches the constant maximum rotating speed omegamaxIf the fan needs to be in d% load reduction operation, only the reference power P is neededrefSwitching to (1-d%) PnThe rotating speed cannot be changed;
the wind speed-rotating speed relational expression under different wind speeds is solved by utilizing the wind speed partition method, and then the rotating speed-inertia relational expression is combined, namely
Figure BDA0003338986730000064
Wherein J represents inertia, ωiAnd representing the actual rotating speed of the fan i to obtain a relation between the available inertia and the wind speed as follows:
Figure BDA0003338986730000071
preferably, the wind speed interval is divided according to the following steps: the preset load shedding target can be realized only by converting the running state of the fan from MPPT (maximum power point tracking) to a d% load shedding curve through overspeed control, and the region is called as a low wind speed region; d% load shedding rate is required to be realized through rotating speed and pitch angle control, and the d% load shedding rate is called as a medium wind speed area; the load reduction target can only be achieved through variable pitch control, and the area is called a high wind speed area.
Preferably, according to the available inertia evaluation process, the wind farm available inertia confidence interval is calculated as follows:
(1) inputting a Copula function and a wind speed-inertia relational expression;
(2) constructing a probability density function f (E) of the available inertia of each cluster, namely obtaining the probability density function according to a formula (5) and a formula (6) and combining a relational expression (15) of wind speed and inertia:
Figure BDA0003338986730000072
(2) integrating f (E) to obtain a probability distribution function phi (E) of the available inertia, i.e.
Figure BDA0003338986730000073
(3) Obtaining inertia values of a probability distribution function phi (E) corresponding to the wind turbine group at alpha/2 and 1-alpha/2 probability values by adopting a binary search-numerical integration method, and respectively recording the upper limit and the lower limit of an available inertia interval of the wind turbine group;
(4) superposing the available inertia of each cluster of the wind power plant to obtain a confidence interval of the available inertia of the full wind power plant, and analyzing the wind speed of each cluster in the cluster by using a wake effect to further obtain the distribution of the available inertia of each cluster;
(5) output available inertia E, and distribution of inertia.
The invention has the beneficial effects that:
the invention provides a wind power plant partition virtual inertia estimation method based on Copula function and clustering algorithm by considering wind power plant wind speed randomness and inter-unit wind speed correlation, and the wind power plant partition virtual inertia estimation method analyzes the internal structure and wind speed data of an actual wind power plant. And obtaining the total available inertia confidence interval of the wind power plant and the available inertia distribution of each unit. A simulation example is constructed according to the actual wind speed and output data of a wind field in Gansu, and the simulation result shows that the algorithm provided by the invention can effectively perform partition clustering on the wind power plant and the accuracy of available inertia estimation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of a partitioned virtual inertia estimation method considering wind speed randomness and correlation according to the present invention.
Fig. 2 is a schematic diagram showing a simplified model of the wake effect.
Fig. 3 shows a processing flow of the wind speed data of the unit.
FIG. 4 shows a fan speed-power curve based on load shedding.
Fig. 5 shows the fan operating state.
FIG. 6 is a diagram illustrating a geographical location profile of a wind farm machine.
FIG. 7 is a graph showing cluster clustering results and daily variation curves of wind speed at the clustering center; (a) the clustered unit geographical position distribution map is obtained, and the unit with the same color is a cluster; (b) a first machine group; (c) a second machine group; (d) and (5) a cluster III.
FIG. 8 shows the histograms of wind speed distributions and the fitted probability density function for W1 and W2; (a) a central fan W1; (b) anemometer tower W2.
FIG. 9 shows a two-bin binary joint wind velocity distribution histogram.
FIG. 10 is a graph of probability density function and distribution function of Gumbel-copula; (a) a probability density function graph; (b) a distribution function graph.
FIG. 11 is a graph of a probability density function and a distribution function of Clayton-copula; (a) a probability density function graph; (b) a distribution function graph.
FIG. 12 is a graph of the probability density function and distribution function of Frank-copula; (a) a probability density function graph; (b) a distribution function graph.
FIG. 13 is a diagram showing a comparison of inertia available from an equivalent machine method and a real measurement method.
Fig. 14 shows the actual values of the available inertia and the theoretical confidence intervals of the available inertia (confidence is 90%).
FIG. 15 is a graph showing a usable inertia profile; (a) actual value of available inertia; (b) the available inertia desired value profile.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, 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 some embodiments of the present invention, but not all embodiments. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. 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 embodiment specifically provides a partition virtual inertia estimation method considering wind speed randomness and correlation, as shown in fig. 1, including the following steps:
step 1: considering wake flow and time delay effects of wind speed, processing historical wind speed data by adopting a segmented aggregation approximation method, and establishing daily variation curves of the wind speed of each unit in the field;
step 2: clustering the daily variation curve of the wind speed of each unit by adopting a dual-scale spectral clustering algorithm to realize the partition of the units in the field and form a plurality of clusters;
and step 3: fitting wind speed probability distribution of each cluster center unit by adopting a nonparametric kernel density estimation method (solving wind speed edge distribution after fitting the probability distribution of the center units by adopting the nonparametric kernel density estimation method, providing necessary conditions for solving parameters in a Copula function later), and constructing an optimal Copula function to analyze the correlation of the wind speed between each cluster center and a anemometer tower;
and 4, step 4: and (4) considering uncertainty between wind speed and inertia, and estimating virtual inertia storage of each cluster in the wind power plant based on the internal structure and wind speed data of the actual wind power plant.
In step 1, the difference in wind speed between different units in the wind farm is mainly affected by the wake effect between the units.
The simplified wake flow model is shown in fig. 2, and the wind speed relationship between the upstream and downstream units is shown in the following formula.
Figure BDA0003338986730000101
In the formula VoThe wind speed is the wind speed blowing to the wind turbine generator set; vxThe wind speed d leaving the wind turbine set is the wind speed reduction coefficient, the size and the thrust coefficient C of the wind turbine setTThe radius R of the blades of the wind turbine generator, the distance X between adjacent wind turbine generators and the wake descent coefficient K are related. The land wind farm generally takes K as 0.075.
In a large wind power plant, there is an obvious time delay in the real-time wind condition transmitted from an upstream fan to a downstream fan, and the time delay phenomenon of the wind speed is particularly obvious under the condition of considering the wake effect. Assuming that the delay time of the wind speed from the wind farm end to the wind turbine j is tau, the wind speed of the upstream wind turbine considering the delay is:
vj(t)=vo(t-τ)
as shown in fig. 3, when the wind speed data information amount is large, firstly, a piecewise aggregation approximation method is adopted to perform dimensionality reduction characterization on the wind speed curve. The basic form of high-dimensional data is approximately expressed by low-dimensional data, the dimension reduction processing of wind speed data can be realized, and a wind speed sequence Y with the length of n is { Y ═ Y1,y2,…,ynUsing the wind speed number series with the length of m
Figure BDA0003338986730000111
To approximate it, where m<n and m is divisible by n. Then
Figure BDA0003338986730000112
The ith element in (b) can be calculated by:
Figure BDA0003338986730000113
secondly, clustering is carried out on the wind speed change curve based on a dual-scale spectral clustering algorithm, and the specific process is as follows:
step 201: distance-based similarity metric
Set X (X) with M wind speed data1,x2,…,xM) Wherein each wind speed curve xiHas a data length of L; the Euclidean distance is the real distance between two data points in the L-dimensional space, and the wind speed curve xiAnd xjThe euclidean distance between them is defined as:
Figure BDA0003338986730000114
step 202: similarity measure based on morphological features
The morphological characteristics represent the fluctuation degree of the curve, and the information entropy is used as a measurement index of the fluctuation degree of the wind speed curve; setting a wind speed curve x in a certain time periodiThere are n possible values of x1,x2,…,xn(ii) a Probability of each possible value being p1,p2,…,pn(ii) a The wind speed curve x for that time periodiEntropy value H ofnIs defined as:
Figure BDA0003338986730000115
the size of the entropy reflects the fluctuation degree of the curve, and the larger the entropy is, the larger the fluctuation degree of the curve is; similarly, the closer the information entropy values of the two wind speed curves are, the higher the fluctuation similarity of the two curves is;
the wind speed similarity matrix of each fan is constructed by extracting similarity measurement features based on distance and form, and daily variation curves of wind speeds close to each other in distance and similar in fluctuation degree are gathered into a class by utilizing a spectral clustering algorithm, so that the partition of the wind power plant is realized.
The wind turbine clustering algorithm divides all the turbines in the farm into a plurality of turbine groups based on wind speed distribution in the farm, and establishes an available inertia evaluation model for each divided turbine group, so that the comprehensive reflection capability of the characteristics of the whole farm of the wind power plant is improved, and the accuracy of the available inertia evaluation of the whole farm is improved.
And (4) integrating the similarity measurement based on the distance and morphological characteristics, and clustering the wind speed curves of the units by using a dual-scale spectral clustering method to realize cluster grouping of the units in the field. The calculation process is as follows:
inputting: inputting wind speed data sets of each unit
Clustering:
1. determining a width parameter of a Gaussian kernel function so as to construct a similarity matrix H;
2. calculating a Laplace matrix of the similarity matrix and determining an optimal classification number k;
3. obtaining the eigenvectors corresponding to the first k maximum eigenvalues;
4. forming a feature matrix for the selected feature vectors;
5. and carrying out K-means clustering on the characteristic matrix to obtain a curve clustering result.
In step 3, the non-parametric kernel density estimation method is a non-parametric estimation method. Let { z i1,2, …, n is a sample space of a random variable z, the probability density of z is u (z), and the kernel density of u (z) is estimated as:
Figure BDA0003338986730000121
Figure BDA0003338986730000122
wherein z is a kernel density function argument; z is a radical ofiData for sample point i; l is the window width; k (delta) is a Gaussian kernel function, and delta is a kernel function argument.
After the in-field unit clustering and dimensionality reduction is carried out in the step 2, the whole wind power plant is divided into N unit areas after the in-field unit clustering and dimensionality reduction, a wind speed sequence joint distribution probability histogram between each unit area and a wind measuring tower is drawn, tail characteristics and symmetry characteristics of the histogram are analyzed, and an optimal Copula function is selected according to a method of minimizing the Euclidean distance from an empirical Copula function.
The specific process is as follows:
step 301: first, an empirical Copula function, a random variable X, Y, and empirical distribution functions of sum F, respectively, are introducedn(x) And Hn(x) The empirical Copula function is defined as:
Figure BDA0003338986730000131
in which u, v is E [0,1 ]];Fn(xi) When the content is less than or equal to u,
Figure BDA0003338986730000132
otherwise, the value is 0;
obtaining relevant unknown parameters of each theoretical Copula function by using a maximum likelihood estimation method; according to the formula (1), the euclidean distance squared of the theoretical Copula function and the empirical Copula function is obtained by combining the defined formula of the euclidean distance and the empirical Copula function as follows:
Figure BDA0003338986730000133
wherein P is a theoretical Copula function type,
Figure BDA0003338986730000134
reflects the squared Euclidean distance between a theoretical Copula function and an empirical Copula function,
Figure BDA0003338986730000135
the Copula function corresponding to the minimum value is an optimal Copula function;
after the optimal Copula function is selected, the binary joint probability density function between the cluster i and the anemometer tower m is as follows:
fi(vi,vm)=Ci(uvi,vvm)f(vi)f(vm) (5)
when inertia estimation is carried out, detailed wind speed data measured by the anemometer tower are known, an optimal Copula function is established to describe the wind speed correlation between the anemometer tower and each cluster, and therefore the instantaneous wind speed probability density function of each cluster is obtained as follows:
Figure BDA0003338986730000136
in step 4, the wind speed-inertia uncertainty analysis process of the wind power plant is as follows:
the wind driven generator makes the running state deviate from the maximum power point by controlling the rotation speed or the pitch angle of the wind driven generator, and then the reserved standby power control mode is the load shedding control of the fan. The load shedding control mode mainly comprises overspeed control and pitch angle control.
The available inertia estimation process between the partitions based on the wind speed is as follows:
wind speed interval division basis: as shown in fig. 4, in the area surrounded by ABB 'a', the predetermined load shedding target can be realized only by converting the operating state of the fan from MPPT to a d% load shedding curve through overspeed control, which is called as a low wind speed area; the BCB' area is limited by the maximum rotating speed of the fan, and d% load shedding rate is required to be realized through rotating speed and pitch angle control, and the area is called a medium wind speed area; in the figure, the C' D interval reaches the maximum rotating speed of 1.2p.u, and the load reduction target can be achieved only through variable pitch control, and the area is called as a high wind speed area.
TABLE 1 interval of wind speeds
Figure BDA0003338986730000141
In addition to the first three wind speed zones, as shown in fig. 5, in the startup zone, when the wind speed is below the cut-in wind speed or when the wind speed is above the cut-out wind speed: the wind speed is limited, and the wind turbine cannot normally run in a grid-connected mode or has no rotating speed reduction capacity, so that the wind turbine generator does not participate in frequency response of the power system in the wind speed area, and has no available inertia.
In the low wind speed region, according to the related theory of aerodynamics, when the wind turbine runs in the MPPT operation mode, the mechanical power captured by the wind turbine can be represented by the following formula:
Figure BDA0003338986730000142
where ρ is the air density, CP maxThe wind energy capture coefficient is shown as lambda is the blade tip speed ratio, beta is the pitch angle, R is the fan blade radius, and V is the wind speed; wherein, the wind energy capture coefficient CP maxAssociated with the tip speed ratio λ and pitch angle β:
Figure BDA0003338986730000151
Figure BDA0003338986730000152
in the above formula, omega is the rotating speed of the generator; p is the number of pole pairs of the generator; g is the number of pole pairs of the gear box;
when the load shedding ratio is d%, the power value is shown in the following formula:
Figure BDA0003338986730000153
combining the two formulas to obtain:
CP-de=(1-d%)CP max (11)
the tip speed ratio when the load shedding rate is d% can be obtained according to a bisection method, and the relation between the fan rotating speed and the wind speed when the load shedding rate is d% can be obtained through the definition of the tip speed ratio:
Figure BDA0003338986730000154
when the wind turbine operates in a medium wind speed constant-speed area, the wind speed basically reaches the upper limit, overspeed control cannot continue to function, the rotating speed of a rotor of the wind turbine does not obviously change along with the wind speed, but the power captured by the wind turbine is continuously improved through variable pitch control, and the wind power plant available inertia evaluation method based on a mixed Copula function is based on documents]The partial rotating speed equation can be approximately linearized, and the wind speed and rotating speed value (v) passing through two end points of the curven,wmax)、 (vws,ws) Constructing a straight line, and obtaining a relation of the rotating speed and the wind speed as follows:
Figure BDA0003338986730000155
when the wind speed of the fan reaches the high wind speed constant power region, the rotating speed of the fan reaches the constant maximum rotating speed omegamaxIf necessaryD% load reduction operation is carried out on the fan, and only the reference power P is requiredrefSwitching to (1-d%) PnAnd the rotation speed does not change.
The available inertia of the fan refers to the capacity of the fan, which can participate in inertia response or primary frequency modulation response. When the frequency changes, the fan changes the output of the unit by releasing the kinetic energy of the rotor, reduces the power difference with the load and inhibits the rapid change of the frequency. Let the rotor be an object with inertia J and angular velocity ω rotation, and its rotational kinetic energy be:
Figure BDA0003338986730000161
for a single fan in a wind power plant, when inertia response is carried out, the kinetic energy of a rotor cannot be completely released, and the rotor has the limit of the lowest rotating speed, omega, in the inertia response processmin0.75 (p.u.). Therefore, if the actual rotational speed of the fan i is ωiThen, the actual available inertia of this fan is:
Figure BDA0003338986730000162
according to the equations (13) and (14), the relation between the available inertia and the wind speed is obtained as follows:
Figure BDA0003338986730000163
according to the available inertia probability evaluation process, the calculation steps of the wind power plant available inertia confidence interval are as follows:
(1) inputting a Copula function and a wind speed-inertia relational expression;
(2) constructing a probability density function f (E) of the available inertia of each cluster, namely obtaining the probability density function according to a formula (5) and a formula (6) and combining a relational expression (15) of wind speed and inertia:
Figure BDA0003338986730000164
(2) integrating f (E) to obtain a probability distribution function phi (E) of the available inertia, i.e.
Figure BDA0003338986730000171
(3) Obtaining inertia values of a probability distribution function phi (E) corresponding to the wind turbine group at alpha/2 and 1-alpha/2 probability values by adopting a binary search-numerical integration method, and respectively recording the upper limit and the lower limit of an available inertia interval of the wind turbine group;
(4) superposing the available inertia of each cluster of the wind power plant to obtain a confidence interval of the available inertia of the full wind power plant, and analyzing the wind speed of each cluster in the cluster by using a wake effect to further obtain the distribution of the available inertia of each cluster;
(5) output available inertia E, and distribution of inertia.
Example simulation
Taking a certain wind farm in Gansu as an example, a schematic layout diagram of fans in the wind farm is shown in FIG. 6: the horizontal and vertical coordinates in the figure are longitude and latitude coordinates of the position of the fan, and the lower right corner is the position of the wind measuring tower. And (3) evaluating the available inertia interval of the whole wind power plant by adopting the accurate wind speed data at the anemometer tower, and selecting the wind speed data with the height of 70m of the anemometer tower when modeling the wind speed of the fan, wherein the fan is of a uniform type, and the detailed data are shown in a table 2.
TABLE 2 wind turbine parameters
Figure BDA0003338986730000172
1.1 Cluster clustering
Clustering is carried out on the daily change curves of the wind speeds of the 33 fans by using a dual-scale spectral clustering algorithm, the clustering number is 3, and the result is shown in a figure 7 below. And obtaining the position distribution maps of the three clustered clusters by combining the positions of the clusters according to the curve clustering result.
1.2 optimal Copula function
According to the method for analyzing the cluster wind speed correlation by the Copula function, the wind speed data of a certain cluster center fan W1 and a wind measuring tower W2 are adopted for example analysis. The wind speed data is instantaneous wind speed data measured at 10-minute intervals on a certain day by the unit.
In the embodiment, nonparametric kernel density estimation method is adopted to determine wind speed edge distribution U and V of two units, and frequency histograms corresponding to wind speeds W1 and W2 are plotted as shown in FIG. 8.
From the binary distribution histogram (fig. 9), the tail of the probability density function is asymmetric, and an archimedes-Copula function can be selected to describe the relationship between the two. The parameters of three copula functions obtained by the maximum likelihood estimation method are as follows: gumbel-copula: 5.7320, respectively; clayton-copula: 4.4871, respectively; frank-copula: 20.1128.
as can be seen from the probability density function graphs (FIG. 10, FIG. 11 and FIG. 12) of the three copula functions, Gumbel-copula and Frank-copula functions have higher similarity in shape, which indicates that the two functions can be used for characterizing the correlation between the wind speeds of the two machine sets. For more accurate selection of the appropriate copula function, the above euclidean distance method is used for comparison. See table 3:
TABLE 3 correlation measurement coefficients of three copula functions and their Euclidean distances
Figure BDA0003338986730000181
Among the three, the Gumbel-Copula function has a small squared Euclidean distance, and a Clayton-Copula function can be selected in the combined distribution of the Gumbel-Copula function and the Clayton-Copula function, and can better show the relationship between the Gumbel-Copula function and the Clayton-Copula function.
1.3 available inertia confidence intervals and distributions
And describing the wind speed relativity of the three clusters by using the optimal copula function between the three clusters and the anemometer tower obtained in the step 1.2, and deducing the wind speed of each cluster by using the known anemometer tower wind speed data. And then solving the inertia of each cluster by using a wind speed-inertia relational expression. The wind speed distribution inside the cluster can be obtained by combining the wake effect and the time delay effect, and further the available inertia distribution can be obtained.
(1) Wind power plant available inertia evaluation curve based on time sequence
A wind power plant in Gansu is selected as an example, the method provided by the invention is used for evaluating the upper and lower limits of the total available inertia of the wind power plant, and the upper and lower limits are compared with historical data and actual available inertia calculated by average wind speed. The effectiveness of the method in practical application is illustrated.
And in actual operation, performing available inertia evaluation by using the historical wind speed data of 0:00-24:00 of a certain day in summer of the wind power plant M. And compares the available inertia interval estimated by the method of the invention with the available inertia curve based on the measured value. The lower graph is an available inertia change curve obtained by actually measuring the rotating speed value of the fan.
And comparing the available inertia confidence interval obtained by the method with the actual value.
Fig. 13 is a time-dependent variation curve of total available inertia of a wind farm in one day, and in combination with a 1.1 daily variation curve of wind speed, it can be seen that the use of the equivalent machine method has higher accuracy when the wind speed is higher, because when the wind speed is too high, the wind turbine enters a constant rotation speed or constant power region, and the rotation speed of the wind turbine does not change with the change of the wind speed any more, so that the rotation speeds of the units are basically consistent and are the highest rotation speeds, so the equivalent machine method works, however, as the wind speed continuously decreases, within a time period of 0.5 to 0.8 days, the accuracy of the equivalent machine method is lower, even the available inertia is 0, because many units cannot provide inertia support when the wind speed is too low. This is also the reason why the equivalent machine method cannot accurately estimate the available inertia when the wind speed is low.
FIG. 14 shows the upper and lower estimated useful inertia limits when the confidence level is 90%, and it can be seen that the probability method used herein to estimate the confidence interval of the useful inertia can better contain the actual useful inertia value, especially in the case of medium and high wind speed, but sometimes inaccurate. In 144 groups of wind speed time sequence data in one day, 10 groups with actual inertia exceeding the upper limit of a confidence interval or being lower than the lower limit of the confidence interval exist, and the error rate is 6.944%. Meanwhile, the method for solving the available inertia has good accuracy in a low wind speed interval. When the wind speed is too low, the available inertia cannot be estimated by the equivalent machine method, but the method can accurately identify the available inertia value, but also has the problem of overhigh available inertia estimation at a certain time point, and in practical application, the principle that the estimation value needs to be kept as much as possible is considered, so that the lower limit curve of the available inertia estimation has better reference value.
(2) Available spatial distribution of inertia
The distribution of the available inertia of the wind power plant is calculated by using the wind speed data of each unit of the wind power plant at the time of 10:00 of the M wind power plant and is shown in FIG. 15: fig. 15(a) is an actual value of the available inertia, and fig. 15(b) is a distribution diagram of an expected value of the available inertia obtained by the method of the present invention, which is obtained by comparing the two diagrams, and the distribution of the available inertia obtained by the method of the present invention is substantially consistent with the actual value, which cannot be realized by the equivalent machine method. The average error of the inertia distribution calculated by the method is 4.0552% by calculating the error between the theoretical value and the actual value of the available inertia of each unit. The effectiveness of the method provided by the invention in available inertia distribution estimation is fully proved. Meanwhile, dispatching personnel can master the potential of output provided by each unit when the inertia is supported according to the graph, and a basis is provided for adjusting the frequency response and the operation state of the unit.
The above description is only for the purpose of illustrating the technical solutions of the present invention and not for the purpose of limiting the same, and other modifications or equivalent substitutions made by those skilled in the art to the technical solutions of the present invention should be covered within the scope of the claims of the present invention without departing from the spirit and scope of the technical solutions of the present invention.

Claims (7)

1. The partition virtual inertia estimation method considering the randomness and the correlation of the wind speed is characterized by comprising the following steps of:
step 1: considering wake flow and time delay effects of wind speed, processing historical wind speed data by adopting a segmented aggregation approximation method, and establishing daily variation curves of the wind speed of each unit in the field;
step 2: clustering the daily variation curve of the wind speed of each unit by adopting a dual-scale spectral clustering algorithm to realize the partition of the units in the field and form a plurality of clusters;
and step 3: fitting the wind speed probability distribution of each cluster center unit by adopting a nonparametric kernel density estimation method, and constructing an optimal Copula function to analyze the correlation of the wind speed between each cluster center and a wind measuring tower;
and 4, step 4: and (4) considering uncertainty between wind speed and inertia, and estimating virtual inertia storage of each cluster in the wind power plant based on the internal structure and wind speed data of the actual wind power plant.
2. The method for estimating virtual inertia of a partition in consideration of randomness and correlation of wind speed according to claim 1, wherein in step 2, the partition process of the wind farm is implemented as follows:
step 201: distance-based similarity metric
Set X (X) with M wind speed data1,x2,…,xM) Wherein each wind speed curve xiHas a data length of L; the Euclidean distance is the real distance between two data points in the L-dimensional space, and the wind speed curve xiAnd xjThe euclidean distance between them is defined as:
Figure FDA0003338986720000011
step 202: similarity measure based on morphological features
The morphological characteristics represent the fluctuation degree of the curve, and the information entropy is used as a measurement index of the fluctuation degree of the wind speed curve; setting a wind speed curve x in a certain time periodiThere are n possible values of x1,x2,…,xn(ii) a Probability of each possible value being p1,p2,…,pn(ii) a The wind speed curve x for that time periodiEntropy value H ofnIs defined as:
Figure FDA0003338986720000012
the size of the entropy reflects the fluctuation degree of the curve, and the larger the entropy is, the larger the fluctuation degree of the curve is; similarly, the closer the information entropy values of the two wind speed curves are, the higher the fluctuation similarity of the two curves is;
similarity measurement characteristic extraction based on distance and form is adopted to construct a similarity matrix of wind speed of each fan, and wind speed daily variation curves close to the distance and similar in fluctuation degree are gathered into one type by utilizing a spectral clustering algorithm, so that the partition of the wind power plant is realized.
3. The method for estimating virtual inertia of partitions considering randomness of wind speed and correlation according to claim 2, wherein in step 3, after clustering and dimensionality reduction of the units in the wind farm are performed in step 2, the whole wind farm is divided into N unit areas, a wind speed sequence joint distribution probability histogram between each unit area and a wind measuring tower is drawn, tail characteristics and symmetry characteristics of the histogram are analyzed, and an optimal Copula function is selected according to a method of minimizing an Euclidean distance from an empirical Copula function.
4. The method for estimating virtual inertia of a partition in consideration of randomness and correlation of wind speed according to claim 3, wherein the specific process of constructing the optimal Copula function to analyze the correlation of wind speed between each cluster center and the anemometer tower is as follows:
step 301: first, an empirical Copula function, a random variable X, Y, and empirical distribution functions of sum F, respectively, are introducedn(x) And Hn(x) The empirical Copula function is defined as:
Figure FDA0003338986720000021
in which u, v is E [0,1 ]];Fn(xi) When the content is less than or equal to u,
Figure FDA0003338986720000022
otherwise, the value is 0;
obtaining relevant unknown parameters of each theoretical Copula function by using a maximum likelihood estimation method; according to the formula (1), the euclidean distance squared of the theoretical Copula function and the empirical Copula function is obtained by combining the defined formula of the euclidean distance and the empirical Copula function as follows:
Figure FDA0003338986720000023
wherein P is a theoretical Copula function type,
Figure FDA0003338986720000024
reflects the squared Euclidean distance between a theoretical Copula function and an empirical Copula function,
Figure FDA0003338986720000031
the Copula function corresponding to the minimum value is an optimal Copula function;
after the optimal Copula function is selected, the binary joint probability density function between the cluster i and the anemometer tower m is as follows:
fi(vi,vm)=Ci(uvi,vvm)f(vi)f(vm) (5)
when inertia estimation is carried out, detailed wind speed data measured by the anemometer tower are known, and an optimal Copula function is established to describe the wind speed correlation between the anemometer tower and each cluster, so that the probability density function of the instantaneous wind speed of each cluster is obtained as follows:
Figure FDA0003338986720000032
5. the partition virtual inertia estimation method considering wind speed randomness and correlation according to claim 4, wherein in step 4, the available inertia estimation process between partitions based on wind speed is as follows:
according to the load shedding control mode of the wind power plant, the wind power plant is divided into a low wind speed area, a medium wind speed area and a high wind speed area;
in the low wind speed region, according to the related theory of aerodynamics, when the wind turbine runs in the MPPT operation mode, the mechanical power captured by the wind turbine can be represented by the following formula:
Figure FDA0003338986720000033
where ρ is the air density, CPmaxThe wind energy capture coefficient is shown as lambda is the blade tip speed ratio, beta is the pitch angle, R is the fan blade radius, and V is the wind speed; wherein, the wind energy capture coefficient CPmaxAssociated with the tip speed ratio λ and pitch angle β:
Figure FDA0003338986720000034
Figure FDA0003338986720000035
in the above formula, omega is the rotating speed of the generator; p is the number of pole pairs of the generator; g is the number of pole pairs of the gear box;
when the load shedding ratio is d%, the power value is shown in the following formula:
Figure FDA0003338986720000041
combining the two formulas to obtain:
CP-de=(1-d%)CPmax (11)
the tip speed ratio when the load shedding rate is d% can be obtained according to a bisection method, and the relation between the fan rotating speed and the wind speed when the load shedding rate is d% can be obtained through the definition of the tip speed ratio:
Figure FDA0003338986720000042
when the fan operates in a medium-wind-speed constant-speed area, the wind speed basically reaches the upper limit, overspeed control cannot continue to act, and the rotating speed of a rotor of the fan cannot follow the wind at the momentThe speed is obviously changed, but the power captured by the fan is continuously improved through variable pitch control, the partial rotating speed equation can be approximately linearized, and the wind speed and rotating speed value (v) of two end points of the curve is usedn,wmax)、(vws,ws) Constructing a straight line, and obtaining a relation of the rotating speed and the wind speed as follows:
Figure FDA0003338986720000043
when the wind speed of the fan reaches the high wind speed constant power region, the rotating speed of the fan reaches the constant maximum rotating speed omegamaxIf the fan needs to be in d% load reduction operation, only the reference power P is neededrefSwitching to (1-d%) PnThe rotating speed cannot be changed;
the wind speed-rotating speed relational expression under different wind speeds is solved by utilizing the wind speed partition method, and then the rotating speed-inertia relational expression is combined, namely
Figure FDA0003338986720000044
Wherein J represents inertia, ωiAnd representing the actual rotating speed of the fan i to obtain a relation between the available inertia and the wind speed as follows:
Figure FDA0003338986720000051
6. the partition virtual inertia estimation method considering wind speed randomness and correlation as claimed in claim 5, wherein the wind speed interval division is based on: the preset load shedding target can be realized only by converting the running state of the fan from MPPT (maximum power point tracking) to a d% load shedding curve through overspeed control, and the region is called as a low wind speed region; the d% load shedding rate is required to be realized through the control of the rotating speed and the pitch angle, and the d% load shedding rate is called as a medium wind speed area; the load reduction target can only be achieved through variable pitch control, and the area is called a high wind speed area.
7. The method for estimating virtual inertia of a partition taking into account randomness and correlation of wind speed according to claim 5, wherein the confidence interval of available inertia of the wind farm according to the available inertia estimation procedure is calculated as follows:
(1) inputting a Copula function and a wind speed-inertia relational expression;
(2) constructing a probability density function f (E) of the available inertia of each cluster, namely obtaining the probability density function according to a formula (5) and a formula (6) and combining a relational expression (15) of wind speed and inertia:
Figure FDA0003338986720000052
(2) integrating f (E) to obtain a probability distribution function phi (E) of the available inertia, i.e.
Figure FDA0003338986720000053
(3) Obtaining inertia values of a probability distribution function phi (E) corresponding to the wind turbine group at alpha/2 and 1-alpha/2 probability values by adopting a binary search-numerical integration method, and respectively recording the upper limit and the lower limit of an available inertia interval of the wind turbine group;
(4) superposing the available inertia of each cluster of the wind power plant to obtain a confidence interval of the available inertia of the full wind power plant, and analyzing the wind speed of each cluster in the cluster by using a wake effect to further obtain the distribution of the available inertia of each cluster;
(5) output available inertia E, and distribution of inertia.
CN202111302789.2A 2021-11-05 2021-11-05 Partition virtual inertia estimation method considering wind speed randomness and correlation Active CN113935247B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111302789.2A CN113935247B (en) 2021-11-05 2021-11-05 Partition virtual inertia estimation method considering wind speed randomness and correlation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111302789.2A CN113935247B (en) 2021-11-05 2021-11-05 Partition virtual inertia estimation method considering wind speed randomness and correlation

Publications (2)

Publication Number Publication Date
CN113935247A true CN113935247A (en) 2022-01-14
CN113935247B CN113935247B (en) 2024-05-31

Family

ID=79285788

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111302789.2A Active CN113935247B (en) 2021-11-05 2021-11-05 Partition virtual inertia estimation method considering wind speed randomness and correlation

Country Status (1)

Country Link
CN (1) CN113935247B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115238933A (en) * 2022-09-23 2022-10-25 西安德纳检验检测有限公司 Wind turbine generator inertia response detection method, device and system based on multipoint measurement
CN115310388A (en) * 2022-10-13 2022-11-08 南京理工大学 Method for calculating three-dimensional asymmetric double-Gaussian wake flow wind speed of wind turbine with space change
CN117788165A (en) * 2024-02-27 2024-03-29 薪海数字科技(烟台)有限公司 Enterprise supply chain transaction risk prediction method and system based on artificial intelligence

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106528912A (en) * 2016-09-19 2017-03-22 国网浙江省电力公司经济技术研究院 Method for estimating frequency regulation capacity of wind power plant
CN110826927A (en) * 2019-11-11 2020-02-21 国网山东省电力公司日照供电公司 Wind power plant available inertia evaluation method
CN111709454A (en) * 2020-05-22 2020-09-25 东南大学 Multi-wind-field output clustering evaluation method based on optimal copula model
CN112117768A (en) * 2020-03-26 2020-12-22 广西大学 Wind generating set subsection frequency modulation control method based on power tracking curve switching
WO2021164112A1 (en) * 2020-02-18 2021-08-26 山东大学 Frequency control method and system during using wind farm as black-start power source by means of optimal configuration of energy storage
CN113346517A (en) * 2021-05-11 2021-09-03 国网甘肃省电力公司电力科学研究院 Damping support strategy of virtual synchronous machine

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106528912A (en) * 2016-09-19 2017-03-22 国网浙江省电力公司经济技术研究院 Method for estimating frequency regulation capacity of wind power plant
CN110826927A (en) * 2019-11-11 2020-02-21 国网山东省电力公司日照供电公司 Wind power plant available inertia evaluation method
WO2021164112A1 (en) * 2020-02-18 2021-08-26 山东大学 Frequency control method and system during using wind farm as black-start power source by means of optimal configuration of energy storage
CN112117768A (en) * 2020-03-26 2020-12-22 广西大学 Wind generating set subsection frequency modulation control method based on power tracking curve switching
CN111709454A (en) * 2020-05-22 2020-09-25 东南大学 Multi-wind-field output clustering evaluation method based on optimal copula model
CN113346517A (en) * 2021-05-11 2021-09-03 国网甘肃省电力公司电力科学研究院 Damping support strategy of virtual synchronous machine

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
XIPING MA, ET AL.: "Secondary control strategy of islanded micro-grid based on multi-agent consistency", IEEE, 4 January 2018 (2018-01-04), pages 1 - 6 *
柯贤波;张文朝;李朋旺;牛拴保;盛四清;杨俊炜;: "高风电渗透率系统的模糊自适应虚拟惯量控制", 电网技术, no. 06, 5 June 2020 (2020-06-05) *
马喜平等: "考虑相关性的风光抽蓄互补发电系统优化运行", 电力系统及其自动化学报, vol. 31, no. 11, 19 March 2019 (2019-03-19), pages 92 - 102 *
魏博;邵冲;张柏林;汤文;沈渭程;: "基于下垂特性的风电场参与电网快速频率调整实测分析", 电气技术, no. 06, 15 June 2020 (2020-06-15) *
魏博等: "基于下垂特性的风电场参与电网快速频率调整实测分析", 电气技术, no. 6, 30 June 2020 (2020-06-30), pages 39 - 44 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115238933A (en) * 2022-09-23 2022-10-25 西安德纳检验检测有限公司 Wind turbine generator inertia response detection method, device and system based on multipoint measurement
CN115238933B (en) * 2022-09-23 2022-12-09 西安德纳检验检测有限公司 Wind turbine generator inertia response detection method, device and system based on multipoint measurement
CN115310388A (en) * 2022-10-13 2022-11-08 南京理工大学 Method for calculating three-dimensional asymmetric double-Gaussian wake flow wind speed of wind turbine with space change
CN115310388B (en) * 2022-10-13 2022-12-23 南京理工大学 Method for calculating three-dimensional asymmetric double-Gaussian wake flow wind speed of wind turbine with space variation
CN117788165A (en) * 2024-02-27 2024-03-29 薪海数字科技(烟台)有限公司 Enterprise supply chain transaction risk prediction method and system based on artificial intelligence
CN117788165B (en) * 2024-02-27 2024-05-14 薪海数字科技(烟台)有限公司 Enterprise supply chain transaction risk prediction method and system based on artificial intelligence

Also Published As

Publication number Publication date
CN113935247B (en) 2024-05-31

Similar Documents

Publication Publication Date Title
CN113935247B (en) Partition virtual inertia estimation method considering wind speed randomness and correlation
CN107909211B (en) Wind field equivalent modeling and optimization control method based on fuzzy c-means clustering algorithm
CN109751206B (en) Fan blade icing fault prediction method and device and storage medium
CN106875033B (en) Wind power cluster power prediction method based on dynamic self-adaption
Hu et al. Adaptive confidence boundary modeling of wind turbine power curve using SCADA data and its application
Zhu et al. A prediction model for wind farm power generation based on fuzzy modeling
Xu et al. Correlation based neuro-fuzzy Wiener type wind power forecasting model by using special separate signals
CN112186761B (en) Wind power scene generation method and system based on probability distribution
CN113205210A (en) Wind speed and power prediction method, system, equipment and storage medium for wind power plant with complex terrain
CN114386718A (en) Wind power plant output power short-time prediction algorithm combined with particle swarm neural network
CN110991701A (en) Wind power plant fan wind speed prediction method and system based on data fusion
CN115510043A (en) Wind power curve abnormal data eliminating method
CN117967497A (en) Dual-mode control method and system based on wind speed timely monitoring wind generating set
CN112211794B (en) Cabin temperature abnormity early warning method, device, equipment and medium of wind turbine generator
CN109802634A (en) A kind of intelligent O&amp;M method and operational system of the photovoltaic plant based on big data
CN113610285A (en) Power prediction method for distributed wind power
CN106451527B (en) Double-fed fan motor field group of planes polymerization and system based on rotor current dynamic characteristic
CN110273818B (en) Fan blade icing fault monitoring method based on shaft transformation thickness classification
CN117039926A (en) Wind power plant multi-cluster available inertial measurement unit and method
CN117060430A (en) Improved LSTM network-based grid-structured wind power frequency modulation prediction method
CN115879715A (en) Response characteristic considered equivalent aggregation method for polymorphic flexible resources of virtual power plant
CN114320770A (en) Effective wind speed measuring method and device for wind turbine generator
Cascianelli et al. Experimental prediction intervals for monitoring wind turbines: An ensemble approach
CN113269430A (en) Wind power project generating capacity post-evaluation method based on actual measurement data of operation anemometer tower
Yang et al. Wind-storage combined system based on just-in-time-learning prediction model with dynamic error compensation

Legal Events

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