This application claims all benefits accruing under 35 U.S.C. §119 from China Patent Application 201410247853.5, filed on Jun. 5, 2014 in the China Intellectual Property Office, the disclosure of which is incorporated herein by reference.
1. Technical Field
The present disclosure relates to an approach to assess available wind resource, especially to an approach to assess available wind resource distribution based on interpolation method.
2. Description of the Related Art
In large-scale wind power grid background, wind power has taken huge challenge to traditional electrical systems in safe and stable operation because of its volatility and randomness. As a basis for various research fields of wind power, performance analysis and modeling of wind resources distribution has become particularly important.
Interpolation model is usually used for meteorological data downscaling, data assimilation and initial field construction. The inverse distance weighting method is most widely used in interpolation model. However, this method relies on the physical similarity and anti-correlation of the distances. Because the impact of external terrain and weather factors on models is ignored, and the introduction of the local topography parameter description is merely considered, the accuracy is limited. Furthermore, the complex sequence features will greatly increase the complexity of the model and has more requirements to the sequence. Thus, statistical models are often difficult to describe the terrain, weather and other factors of the law on the resource characteristics, and fail to character the intrinsic physical principles.
What is needed, therefore, is an approach to assess available wind resource distribution that can overcome the above-described shortcomings.
Many aspects of the embodiments can be better understood with reference to the following drawings. The components in the drawings are not necessarily drawn to scale, the emphasis instead being placed upon clearly illustrating the principles of the embodiments. Moreover, in the drawings, like reference numerals designate corresponding parts throughout the several views.
The disclosure is illustrated by way of example and not by way of limitation in the figures of the accompanying drawings in which like references indicate similar elements. It should be noted that references to “an” or “one” embodiment in this disclosure are not necessarily to the same embodiment, and such references mean at least one.
Referring to
Step S1, calculating a correlation coefficient between each two of a plurality of anemometer towers, and forming a plurality of first groups by grouping the plurality of anemometer towers for the first time with the shortest distance clustering method, wherein values of wind velocity and wind direction of each two of the plurality of anemometer towers are taken as two different random variables;
Step S2, forming a plurality of second groups by dividing the plurality of anemometer towers in each of the plurality of first groups for the second time, wherein the plurality of anemometer towers are grouped with absolute value of differences of mean wind velocity between each two of the plurality of anemometer towers in each of the plurality of first groups; and
Step S3, obtaining a plurality of values of wind velocity and wind direction at an target point through at least two anemometer towers by the inverse distance weighting method, wherein the at least two anemometer towers are nearest to the target point.
In step S1, the plurality of first groups can be obtained by clustering the plurality of anemometer towers through the plurality of wind velocity data of the plurality of anemometer towers for the past 6 hours.
The shortest distance clustering method belongs to Hierarchy Cluster. The basic principle is to calculate levels of similarity (or distance) based on a certain criterion among different categories. Then two categories that have the highest similarity and recalculate levels of similarity among different categories are combined. The all steps will be redid until the number of categories becomes the same with a set value. The method of clustering with beeline criterion can realize grouping of data, considering it is well-defined and has a very simple model. The specific method of this kind of clustering is shown in followings.
According to the shortest distance clustering method, the distance between two different categories, for example, the distance Dpq between and Gp and Gq, is the minimum value of distance between unit from Gp and unit from Gq. That is,
Then a distance matrix (m×m) is constructed, in which the number of elements to be clustered is m. After mergering Gp and Gq between which distance is minimum among the matrix elements except the diagonal matrix elements, a first category Gr can be obtained. Namely, Gr={Gp, Gq}. In Gr,
d
rk=min{dpk,dqk} (k≠p,q) (2).
Then the distance between the first categories Gr and other categories such as Gp and Gq is calculated. As a result, a second matrix can be obtained, of which the number of elements is (m−1). Again, the minimum distance dij among the second matrix, merger Gi and Gj is found out, and then the distance between the categories in the second matrix and other categories such as Gi and Gj is calculate. The circulation can be stopped when the number of categories reaches to the set value.
The values of wind velocity and wind direction of two of the plurality of anemometer towers can be taken as two different random variable, such as X and Y. Then the correlation coefficient between X and Y is:
wherein μx is the expected value of X, σx is the standard deviation of X, xi is a plurality of individual data in X; μy is the expected value of Y, σy is the standard deviation of the Y, yi is a plurality of individual data in Y.
In addition, while grouping the plurality of anemometer towers with correlation coefficients, the minimum distance should correspondent to maximum correlation coefficient. As correlation coefficients reflect the similarity of fluctuation between different values of wind velocity and wind direction, and the anemometer towers with similar fluctuation will be grouped together. If the differences of wind resource are obvious, the number of groups should be added, otherwise the number of groups should be reduced.
In step S2, because the correlation coefficient does not reflect the value relationships between wind velocity, the anemometer towers classified as one group may exist differences in wind velocity values. Thus the plurality of anemometer towers in the first groups need to be grouped for the second time. In addition, because there may be great wind velocity magnitude volatility, the number of groups is not pre-specified here and is judged in the grouping process until the maximum differences between mean wind velocity in the same second group less than a threshold.
In Step S3, referring to
Step S3.1, calculating a geographical distance r between the target point and each of the plurality of anemometer towers;
Step S3.2, finding out a first anemometer tower i and a second anemometer tower j, wherein the geographical distances r between the target point and the first anemometer tower i and the second anemometer tower j are minimum;
Step S3.3, confirming a plurality of third anemometer towers from the plurality anemometer towers, wherein the plurality of third anemometer towers, the first anemometer tower i, and the second anemometer tower j belong to the same first group;
Step S3.4, obtaining the plurality of values of velocity and wind direction at target points through the first anemometer tower, the second anemometer tower, and the plurality of third anemometer towers based on the method of inverse distance weighting interpolation.
In step S3.1, after finding out the first anemometer tower i, the second anemometer tower j, the plurality of third anemometer towers within the same first group can be obtained based on step S1.
In step S3.2, inverse distance weighted interpolation method need to select an anemometer tower closest to the target point to be estimated as a grouping standard. Furthermore, select only the most recent one anemometer tower could cause instability in results. Thus by selecting the closest two anemometer towers, the calculation accuracy will be improved through the wind velocity and wind direction interpolation consolidated results.
In step S3.4, a weighting function W(r) is set to show the distance from the target point to the known points such as the first anemometer tower, the second anemometer tower, and the plurality of third anemometer towers, and then the values of wind velocity and wind direction through weighting data of the first anemometer tower, the second anemometer tower, and the plurality of third anemometer towers can be obtained.
The W(r) can be set as:
wherein r is the distance between the target point and the known points, R is influence radius which shows the known point has no reference value to the target point if the distance r is larger than R, m is an integer which is larger than 1. Apparently, the closer the distance between the known point to the target point is, the larger the value of weight will be.
In one embodiment, the model is constructed based on the function (5), which shows the closer the distance is, the higher the level of similarity of resources will be. Then a vector interpolation of first wind velocity U in a first direction and a second wind velocity V in a second direction are:
wherein the first direction U are substantially perpendicular to the second direction V. In one embodiment, the first direction U represents the direction from west to the east, and the second direction V represents the direction from south to the north.
The approach to assess available wind resource distribution based on interpolation method has following advantages. First, the shortest distance method is introduced into the approach to assess available wind resource distribution, and the anemometer towers are grouped based on historical actual sequence of real-time packet masts, thus a higher similarity is guaranteed between the fluctuations on the law and values of the packet mast. Second, during performing interpolation calculation, the nearest at least two masts are selected, thus the calculation accuracy is improved through the wind velocity and direction interpolation consolidated results. Third, the accuracy is improved, and the overall level of the model error is limited by the inherently accuracy of such statistical interpolation method.
The approach to assess available wind resource is applied in the district Jiuquan, GanSu province, China. 22 anemometer towers are selected.
First, referring to
Second, referring to
Third, the wind velocity and wind direction at the target point is obtained through at least two anemometer towers in the same group after grouped for the second time by the inverse distance weighting method.
In comparison, one anemometer tower selected to calculate the wind velocity and wind direction is labeled as No. 1. Two anemometer towers in the same group selected to calculate the wind velocity and wind direction is labeled as No. 2. All anemometer towers in the same group selected to calculate the wind velocity and wind direction is labeled as No. 3. Referring to
Depending on the embodiment, certain of the steps of methods described may be removed, others may be added, and that order of steps may be altered. It is also to be understood that the description and the claims drawn to a method may include some indication in reference to certain steps. However, the indication used is only to be viewed for identification purposes and not as a suggestion as to an order for the steps.
It is to be understood that the above-described embodiments are intended to illustrate rather than limit the disclosure. Variations may be made to the embodiments without departing from the spirit of the disclosure as claimed. It is understood that any element of any one embodiment is considered to be disclosed to be incorporated with any other embodiment. The above-described embodiments illustrate the scope of the disclosure but do not restrict the scope of the disclosure.
Number | Date | Country | Kind |
---|---|---|---|
201410247853.5 | Jun 2014 | CN | national |