A New MCP Method of Wind Speed Temporal Interpolation and Extrapolation Considering Wind Speed Mixed Uncertainty

In this paper, a missing wind speed data temporal interpolation and extrapolation method in the wind energy industry was investigated. Given that traditional methods have previously ignored part of mixed uncertainty of wind speed, a concrete granular computing method is constructed and a new Measure–Correlate–Predict (MCP) method of wind speed data temporal interpolation and extrapolation considering all mixed uncertainties is proposed, based on granular computing theory by adopting the cloud model method, support vector regression method, artificial neural network, genetic algorithm, and fuzzy c-means clustering algorithm as tools. The importance of considering mixed wind speed uncertainty and the suitability of using granular computing method are illustrated, and wind speed mixed uncertainty analysis is implemented, then, recommended values and estimation tools for wind speed measurement uncertainty and combined uncertainty are provided. An interpolation case of two practical meteorological sites in central Southern China was used to implement and validate the method proposed in this paper. The following conclusions are reached: (a) by using the method proposed in this paper, mixed uncertainty of wind speed can be considered, comparing to other MCP methods used for purposes of comparison, a better estimation of the wind speed is provided, and most evaluation metrics employed in this analysis were superior to other methods, that is to say, the accuracy of the wind resource assessment improved, and the risks of wind farm construction were reduced; (b) granular computing method is suitable for the issue of wind speed data interpolation and extrapolation considering wind speed mixed uncertainty; (c) mixed uncertainty of wind speed can be divided into three levels, and recommended values of granularity are minimum interval of records, 0.3–0.8 m/s, and 1–3 m/s, respectively.


Introduction
With increasing global warming and environmental problems, wind energy has now become one of the most important sources of green energy.Wind resource assessment is a key procedure in wind farm construction, and wind speed observation is the foundation of wind resource assessment [1,2].However, in many practical cases, the use of wind speed observation records is restricted by two factors, one is gaps in the wind speed records caused by instrument failure or the destruction of the meteorological mast, the other is lack of long-term wind speed records.To overcome these Energies 2017, 10, 1231 2 of 21 drawbacks, one or more reference sites have been chosen and the relationship between the target and reference sites have been constructed using statistical methods (Measure-Correlate-Predict (MCP) methods [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]), or physical methods (for example, Wind Atlas Analysis and Application Program (WAsP) [18][19][20]).As per references [13,20], statistical methods tend to provide higher accuracy and are therefore used widely in practical engineering.The MCP methods model correlations between the target and reference sites using simultaneous records; then, those correlations are used with reference site data to predict missing wind data or long-term wind data of the target site.Previous researchers have provided some effective MCP methods that used different types of functions to model correlation: Linear regression [3][4][5][6][7][8][9] has been widely used to model correlation; Derrick [3] proposed MCP method to assess long term wind resources by using linear regression method, and similar works have been done by other researchers [4][5][6][7].MCP methods based on linear regression have been validated in some practical applications [8,9].Support vector machine (SVM) model was also used for correlation, Mohandes et al. [10] introduced the SVM for wind speed prediction and compared it with the multilayer perceptron (MLP) neural networks, the results proved that the SVM model had less root mean square errors than the MLP model.Ji et al. [11] did further research on SVM, a support vector classifier was utilized to estimate the forecasting error and lower mean square error and mean absolute percentage error than traditional SVM method were obtained.Artificial neural networks (ANN) [12][13][14][15] method was another tool to model correlation, Addison et al. [13] investigated the feasibility of using neural networks to make predictions of long term energy yield at a potential wind farm site.Saavedra-Moreno et al. [14] provided very fast training neural-computation techniques for real measure-correlate-predict wind operations in wind farms.Probabilistic method [16,17] also was used to model correlation, Carta et al. [16] proposed the use of a probabilistic Measure-Correlate-Predict (MCP) method to estimate the long-term wind speed characteristics at a potential wind energy conversion site.García-Rojo [17] indicated an automatic method of calculation of a long-term representative wind climate based on the calculation of the joint probability distribution of the wind at a local station and at a long-term reference meteorological mast.
However, there are some shortages and limitations of these existing methods.Wind speed uncertainty is a mixed uncertainty which mainly consist of randomness and fuzziness, specifically, randomness is an aleatory uncertainty due to stochasticity of nature wind speed, and fuzziness is an epistemic uncertainty mainly caused by measurement.The existing methods have been used to treat wind speed data as accurate quantity, and only randomness has been dealt with probability theory or artificial intelligence method, such as regression or ANN, fuzziness of wind speed would always be ignored.Therefore, inaccurate and unscientific predictions have been obtained by ignoring some components of uncertainties, and then errors of wind resource assessment may occur and risks of wind farm construction are raised.
Granular computing theory (GrC) references and imitates human cognition, and extracts the common theory of problem solving; therefore, it could be used to exploit the tolerance for imprecision, uncertainty, and partial truth to achieve tractability, robustness, low cost solutions, and better rapport with reality [21].All components of wind speed mixed uncertainty can be considered and the drawbacks of traditional MCP methods can be overcome by using MCP method based on granular computing theory.
For the purpose of considering wind speed mixed uncertainty, this paper constructed a concrete granular computing method, and proposed a new MCP method of wind speed temporal interpolation and extrapolation.As the foundation of the new method, wind speed uncertainty analysis was implemented, recommended values and estimation tools for wind speed measurement uncertainty and combined uncertainty were provided, and the uncertainty hierarchy was determined.The new method used cloud model, support vector regression (SVR), artificial neural network (ANN), genetic algorithm (GA), and fuzzy c-means clustering algorithm (FCM) as tools, and illustrated a detailed procedure of the new method.In order to implement and validate the new method, a wind data temporal interpolation case was studied where four metrics were used: mean relative error, root mean square error, correlation coefficient, and mean relative error of energy production; and also a comparison of these metrics with those obtained from other methods was implemented.

Granular Computing Theory
GrC references and imitates human cognition, and may be regarded as a label for theories, methodologies, techniques, and tools that make use of granules, such as groups, classes, or clusters of a universe in the process of solving a problem, and is a more philosophical way of thinking than a practical methodology to problem solving [21].Zadeh [22] identified three basic concepts that underlie human cognition: granulation, organization, and causation.Granulation involves the decomposition of a whole into parts; organization involves the integration of parts into a whole; and causation involves the association of cause and effect.There are four steps in the process of using GrC for problem solving, which are partitioning objects into granules, constructing levels, and hierarchy, representing granules and computing with granules [23,24].
The benefits of GrC are evident from its basic guiding principle, which was stated concisely by Zadeh [22] as to "exploit the tolerance for imprecision, uncertainty and partial truth to achieve tractability, robustness, low cost solutions and better rapport with reality".
Granules, levels, and hierarchy are the basic elements of granular computing.Specifically, (a) a granule is a clump of points (objects) drawn together by indistinguishability, similarity, proximity or functionality [22,25]; (b) levels consist of granules, and the properties of granules collectively characterize a level of description and understanding; and (c) hierarchy is the sum of the relationships of levels.

Uncertainty Analysis of Wind Speed
Uncertainty, an inherent characteristic of the nature of wind, turns wind speed into granules, and then granules constitute levels, and then levels constitute hierarchy.Finally, granules, levels, as well as hierarchy all together constitute the granular world.As the foundation of the granular world is based on uncertainty, an uncertainty analysis of wind speed is the foundation of the proposed method.In this section, the hierarchy is determined, and recommended values and estimation tools for wind speed measurement uncertainty as well as combined uncertainty are provided.

Method for Calculating Uncertainty
Uncertainty is defined as a parameter associated with the results of a measurement that characterizes the dispersion of values that can reasonably be attributed to the measurement, and uncertainties are expressed as standard deviations and are denoted standard uncertainties [26].
The Bins method was used to determine the uncertainties as it is a data reduction procedure where wind speed data are divided into 1 m/s continuous intervals (bins) centered on multiples of 1 m/s, with uncertainties calculated in each bin [27].
Uncertainty consists of some components, and in general, the combined standard uncertainties in bin i can be expressed by where c k,i is the sensitivity factor of component k in bin i; u k,i is the standard uncertainty of component k in bin i; M is the number of uncertainty components in each bin; ρ k,l is the correlation coefficient between component k in bin i, and component l in bin i.
To simplify the above expression, the independence of all uncertainty components was assumed, so ρ k,l = 0 when k = l, and ρ k,l = 1 when k = l.As the uncertainty components and combined uncertainty are in the same dimension, c k,i = 1.Combined uncertainty can be simply expressed by Combined standard uncertainties may additionally be expressed by expanded uncertainties.Referring to the International Organization for Standardization (ISO) guide [26] and assuming normal distributions, combined standard uncertainties were multiplied by a coverage factor of 2 in a 95% confidence level.

Hierarchy of Uncertainty
Some sources of wind speed uncertainty exist in the issue of MCP method, they are wind speed measurement, the influence of local obstacles, meso-micro scale weather patterns, stochastic variations in wind speed and direction over time and distance, as well as flight delays and atmospheric stability [9].We should note here that wind speed uncertainty is a mixed uncertainty, which is consisting of randomness and fuzziness.
In practical cases, fuzziness is decided by the campaign of measurement, and they are independent of specific conditions and are most influenced by the instrument of measuring; on the other side, randomness is decided by other sources of wind speed uncertainty, it is strongly dependent on one or multiple specific conditions.
The construction of certain uncertainty levels indicated that the uncertainty components contained in the level are formulated and calculated, so continuity, uniformity, computability, and independence of uncertainty components are demanded.Uncertainty of wind speed measurement, which is represented by fuzziness, meets the requirements and is used to construct a single level; on the other side, randomness is considered in the combined uncertainty rather than a single level.
In this paper, three levels of mixed uncertainty were constructed: bottom level, which ignores all uncertainties; middle level, which includes the measurement uncertainties; and top level, which considers all components of mixed uncertainty.

Uncertainty of Wind Speed Measurement
To determine and characterize every level of uncertainty, some practical cases were collected.Measurement uncertainty represents fuzziness of wind speed.There are four sources of wind speed measurement uncertainty: anemometer operational characteristics, mounting effects, anemometer calibration, and data acquisition.

Operational Characteristics Uncertainty
Anemometers which were used to wind speed measurement are calibrated in wind tunnel, and these wind tunnel calibrations are performed under controlled environmental conditions with a smooth and low turbulent flow.However, during measurements and operation in the field, these instruments are subject to turbulent flow and environmental conditions that can deviate significantly from wind tunnel calibration conditions.The field conditions may significantly influence instrument characteristics and cause operational characteristics uncertainty.
To calculate operational uncertainties, a classification method of an anemometer was provided [27].Known general influence parameters on cup anemometer measurements are: turbulence, air temperature, air density, and average flow inclination angle, these influence parameters must be considered in the classification, hence the classification of an anemometer type is divided into two classes, which is dependent on influence parameters.Details are shown in Table 1.The classification number k was defined and was used in determination of the operational uncertainty, and a classification, referring to a certain class, is expressed by kA, kB, for example 1.7A or 2.5B.The operational standard uncertainty of a cup anemometer may be derived from the classification number assuming a rectangular uncertainty distribution by using Equation (3).
where u m1,i is the operational standard uncertainties in bin i; U i is the center of bin i; and k is the classification number of the anemometer.In this paper, the classification numbers of five popular cup anemometers, which are widely used in practical cases and represent the average uncertainty of cup anemometers, were collected to determine the approximate level of operational uncertainties, and these cup anemometers were measured by ACCUWIND [28].The uncertainties are presented in Table 2 and Figure 1.The classification number k was defined and was used in determination of the operational uncertainty, and a classification, referring to a certain class, is expressed by kA, kB, for example 1.7A or 2.5B.The operational standard uncertainty of a cup anemometer may be derived from the classification number assuming a rectangular uncertainty distribution by using Equation (3).
  where u is the operational standard uncertainties in bin i ; i U is the center of bin i ; and k is the classification number of the anemometer.In this paper, the classification numbers of five popular cup anemometers, which are widely used in practical cases and represent the average uncertainty of cup anemometers, were collected to determine the approximate level of operational uncertainties, and these cup anemometers were measured by ACCUWIND [28].The uncertainties are presented in Table 2 and Figure 1.

Mounting Effects
A lattice mast with a side mounted anemometer is widely used in practical applications, therefore, anemometers are always influenced by flow distortion from both the mast and the boom.Thrust coefficient, C T , which depends upon the porosity of the mast and the drag on the individual members, is assumed to evaluate influence.The average C T of the lattice mast with a square cross-section and shape edge members, square cross-section and round members, triangular cross-section and round members were 0.92, 0.42, and 0.34, respectively.The mounting standard uncertainty can be estimated by Equation ( 4) below [27].
where u m2,i is the mounting standard uncertainty in bin i; C T is the thrust coefficient of mast; L M is the mast leg length; and R is the distance between the observation position and the center of the mast.
The results of the center-line relative wind speed of lattice mast with average C T are shown in Figure 2.
Energies 2017, 10, 1231 6 of 22 Figure 1 shows that the ranges of operational standard uncertainty are about 0.04-0.21m/s in Class A and are 0.13-0.96m/s in Class B, and the average operational standard uncertainties are about 0.1 m/s in Class A and 0.4 m/s in Class B.

Mounting Effects
A lattice mast with a side mounted anemometer is widely used in practical applications, therefore, anemometers are always influenced by flow distortion from both the mast and the boom.
Thrust coefficient, T C , which depends upon the porosity of the mast and the drag on the individual members, is assumed to evaluate influence.The average T C of the lattice mast with a square cross-section and shape edge members, square cross-section and round members, triangular cross-section and round members were 0.92, 0.42, and 0.34, respectively.The mounting standard uncertainty can be estimated by Equation ( 4) below [27].In practical cases, the distance between the observation position and mast center was 3-5 times the mast leg length, therefore, we can conclude from Figure 3 that the range of center-line relative wind speed was 0.969-0.996 of free wind speed.

Uncertainty of the Anemometer Calibration
To obtain an overview of the cup anemometer calibration uncertainty, calibration uncertainties of ten commonly used cup anemometers, which have been quantified by previous researchers [29], were collected.The details are shown in Table 3, and the linear regression of average relative uncertainties of ten anemometers is proposed in Figure 3.In practical cases, the distance between the observation position and mast center was 3-5 times the mast leg length, therefore, we can conclude from Figure 3 that the range of center-line relative wind speed was 0.969-0.996 of free wind speed.

Uncertainty of the Anemometer Calibration
To obtain an overview of the cup anemometer calibration uncertainty, calibration uncertainties of ten commonly used cup anemometers, which have been quantified by previous researchers [29], were collected.The details are shown in Table 3, and the linear regression of average relative uncertainties of ten anemometers is proposed in Figure 3.The relative calibration uncertainties in Table 3 were quantified in the wind tunnels at Otech Engineering, Inc. located in Davis, California [29].
The calibration uncertainty of a commonly used anemometer can be concluded from Figure 3 and take the following expression.
where u m3,i is the calibration standard uncertainties in bin i.
Energies 2017, 10, 1231 8 of 22 The relative calibration uncertainties in Table 3 were quantified in the wind tunnels at Otech Engineering, Inc. located in Davis, California [29].
The calibration uncertainty of a commonly used anemometer can be concluded from Figure 3 and take the following expression.

 
3, 0.000431 0.021578 / 3 where 3, m i u is the calibration standard uncertainties in bin i .

Uncertainty of Data Acquisition
Uncertainty of Data Acquisition may come from transmission, signal conditioning, analog to digital conversion, and data processing in the data acquisition system.
Assuming that the uncertainty of data acquisition system was 0.1% of the wind speed measured range, the standard uncertainty from data acquisition is where , m i u is the wind speed measurement uncertainty in bin i .
In practical conditions, the recommended value of wind speed measurement uncertainty was about 0.3-0.8m/s.

Uncertainty of Data Acquisition
Uncertainty of Data Acquisition may come from transmission, signal conditioning, analog to digital conversion, and data processing in the data acquisition system.
Assuming that the uncertainty of data acquisition system was 0.1% of the wind speed measured range, the standard uncertainty from data acquisition is where u m4,i is the data acquisition standard uncertainties in bin i; and U r is the wind speed measured range.

Combined Uncertainty of Wind Speed Measurement
The combined uncertainty of wind speed measurement can be derived from Equation (3) to Equation ( 6), and take the following expression.
where u m,i is the wind speed measurement uncertainty in bin i.
In practical conditions, the recommended value of wind speed measurement uncertainty was about 0.3-0.8m/s.

Combined Uncertainty
Combined uncertainty of wind speed includes all components of wind speed mixed uncertainty.In this paper, three operational condition classes with two practical meteorological sites were selected to obtain an overview of the combined uncertainty, location, topography, distance, and correlation of the chosen classes are shown in Table 4.The combined uncertainty of correlation can be characterized by standard deviation, which was calculated from a set of simultaneous reference wind speed data at a specific wind speed of target data, for example, when the wind speed of target data was 5.0 m/s, the set of contemporary reference wind speed had a standard deviation of 1.6 m/s, so the combined uncertainty at the wind speed of 5.0 m/s was 1.6 m/s.Results of the combined uncertainty are shown in Figure 4.

Combined Uncertainty
Combined uncertainty of wind speed includes all components of wind speed mixed uncertainty.In this paper, three operational condition classes with two practical meteorological sites were selected to obtain an overview of the combined uncertainty, location, topography, distance, and correlation of the chosen classes are shown in Table 4.The combined uncertainty of correlation can be characterized by standard deviation, which was calculated from a set of simultaneous reference wind speed data at a specific wind speed of target data, for example, when the wind speed of target data was 5.0 m/s, the set of contemporary reference wind speed had a standard deviation of 1.6 m/s, so the combined uncertainty at the wind speed of 5.0 m/s was 1.6 m/s.Results of the combined uncertainty are shown in Figure 4. Hence, it can be concluded from Figure 4 that the range of combined uncertainty at each specific wind speed was 1-3 m/s.

Proposed Method
In this section, a concrete granular computing method is constructed, and a new Measure-Correlate-Predict (MCP) method of wind speed data temporal interpolation and extrapolation considering uncertainties are proposed.The detailed procedure of proposed method is illustrated as follows.

Procedure of Proposed Method
The 5 steps of the proposed method are presented below and details are shown in Figure 5: (a  Hence, it can be concluded from Figure 4 that the range of combined uncertainty at each specific wind speed was 1-3 m/s.

Proposed Method
In this section, a concrete granular computing method is constructed, and a new Measure-Correlate-Predict (MCP) method of wind speed data temporal interpolation and extrapolation considering uncertainties are proposed.The detailed procedure of proposed method is illustrated as follows.

Procedure of Proposed Method
The 5 steps of the proposed method are presented below and details are shown in Figure 5

Determining the Granular Hierarchy
The foundation of wind speed granules is composed of uncertainties.In this paper, the hierarchy of wind speed granules was decided by the structure of uncertainty, of which three levels of granules were set.
Firstly, in the bottom level, all uncertainties were ignored and wind speed data were supposed to be precise, and then this level was made of basic granules, which were the minimum interval of wind speed recording by anemometer.So, actually, the granule computing in the bottom level was equal to traditional methods.In practical conditions, minimum interval of most anemometer records was 0.1 m/s, so, the recommended value of granularity magnitude in this level was 0.1 m/s.

Determining the Granular Hierarchy
The foundation of wind speed granules is composed of uncertainties.In this paper, the hierarchy of wind speed granules was decided by the structure of uncertainty, of which three levels of granules were set.
Firstly, in the bottom level, all uncertainties were ignored and wind speed data were supposed to be precise, and then this level was made of basic granules, which were the minimum interval of wind speed recording by anemometer.So, actually, the granule computing in the bottom level was equal to traditional methods.In practical conditions, minimum interval of most anemometer records was 0.1 m/s, so, the recommended value of granularity magnitude in this level was 0.1 m/s.
Secondly, in the middle level, only the continuous, uniform, computable, and independent uncertainties, which were measurement uncertainties, were considered.Therefore, the magnitude of granularity was calculated by Equation (7).In practical conditions, the recommended value was approximately 0.3-0.8m/s.
Finally, at the top level, all uncertainties in the wind speed interpolation and extrapolation problem were considered, including measurement uncertainty and correlation uncertainty.The magnitude of granularity equal to combined uncertainty, which can be characterized by standard deviation and be calculated from a set of simultaneous reference wind speed data at a specific wind speed of target data (see Section 3.4), was around 1-3 m/s.

Granulation
This step drew together a clump of objects (points) by indistinguishability, similarity, proximity, or functionality to induce granules.The fuzzy c-means clustering algorithm (FCM), which is an efficient and validated clustering algorithm [30] was adopted in this paper.
The steps of the FCM algorithms are presented below [31]: Secondly, in the middle level, only the continuous, uniform, computable, and independent uncertainties, which were measurement uncertainties, were considered.Therefore, the magnitude of granularity was calculated by Equation (7).In practical conditions, the recommended value was approximately 0.3-0.8m/s.
Finally, at the top level, all uncertainties in the wind speed interpolation and extrapolation problem were considered, including measurement uncertainty and correlation uncertainty.The magnitude of granularity equal to combined uncertainty, which can be characterized by standard deviation and be calculated from a set of simultaneous reference wind speed data at a specific wind speed of target data (see Section 3.4), was around 1-3 m/s.

Granulation
This step drew together a clump of objects (points) by indistinguishability, similarity, proximity, or functionality to induce granules.The fuzzy c-means clustering algorithm (FCM), which is an efficient and validated clustering algorithm [30] was adopted in this paper.
The steps of the FCM algorithms are presented below [31]: Normally, the granulations of the target and reference data are independent.At each granule level, the cluster count was identified based on the analysis results of the granular hierarchy, and the initial center of the clusters was generated randomly.Next, the FCM algorithm was used to partition the data into clusters.As the results of the FCM algorithm are sensitive to the initial center of clusters [32,33], the genetic algorithm was adopted to optimize the initial center of clusters and was avoided a partial optimum solution (Figure 6).Normally, the granulations of the target and reference data are independent.At each granule level, the cluster count was identified based on the analysis results of the granular hierarchy, and the initial center of the clusters was generated randomly.Next, the FCM algorithm was used to partition the data into clusters.As the results of the FCM algorithm are sensitive to the initial center of clusters [32,33], the genetic algorithm was adopted to optimize the initial center of clusters and was avoided a partial optimum solution (Figure 6).

Granules Representation
In this study, the cloud model theory was adopted to represent granules.The cloud model transfers the procedure of uncertainty between the qualitative concept and quantitative data representation by using natural language [34], and mainly reflects the fuzziness and randomness of the concept within the affairs and human knowledge in the objective world [35].Suppose that T is the language value of domain u, and mapping

Granules Representation
In this study, the cloud model theory was adopted to represent granules.The cloud model transfers the procedure of uncertainty between the qualitative concept and quantitative data representation by using natural language [34], and mainly reflects the fuzziness and randomness of the concept within the affairs and human knowledge in the objective world [35].
Suppose that T is the language value of domain u, and mapping C T (x) : u →[0, 1], ∀x ∈ u, x → C T (x), then the distribution of C T (x) in u is called the membership cloud of T, or cloud for short.If the distribution of C T (x) is normal, it is named the normal cloud model, which is a random number set that obeys the normal distributive rule and has a stable tendency.It is determined by expectation Ex, entropy En, and super entropy He, and reflect the quantitative characteristics of concept C T (x) [36], Ex determines the center of the cloud, En determines the range of the cloud and, according to "3En" rule [34], about 99.74% of the total cloud drops distribute between [Ex − 3En, Ex + 3En]; He determines the cloud drops' dispersive degree, which means the larger the He is, the more dispersively the cloud drops locate.The sketch map of three digital characteristics of a normal cloud model is shown in Figure 7.For each granule, every wind speed data record was a drop, and the cloud model digital characteristics   , , Ex En He can be obtained from all the drops by using the back cloud generator, thereby the three digital characteristics were representative of each granule.

Granular Computing
Randomness of wind speed was considered in this step.Granular computing is a method of finding the correlation between two worlds of granules, which is complex and non-linear data fitting.There are two alternative methods, artificial neural networks method and support vector regression method, are very suitable for this kind of issue [10][11][12][13][14][15], both of them were adopted as alternative method in this paper and compared in the Section 5. A brief introduction of the two methods was provided below.

Artificial neural networks
The topography of the ANN is shown in Figure 10.Very fast training neural-computation techniques [14] were chosen for ANN training.The ANN was trained by the data of two worlds of granules for each level, and three trained ANNs, corresponding to three levels, were the results of this step.For each granule, every wind speed data record was a drop, and the cloud model digital characteristics (Ex, En, He) can be obtained from all the drops by using the back cloud generator, thereby the three digital characteristics were representative of each granule.

Granular Computing
Randomness of wind speed was considered in this step.Granular computing is a method of finding the correlation between two worlds of granules, which is complex and non-linear data fitting.There are two alternative methods, artificial neural networks method and support vector regression method, are very suitable for this kind of issue [10][11][12][13][14][15], both of them were adopted as alternative method in this paper and compared in the Section 5. A brief introduction of the two methods was provided below.

• Artificial neural networks
The topography of the ANN is shown in Figure 10.Very fast training neural-computation techniques [14] were chosen for ANN training.The ANN was trained by the data of two worlds of granules for each level, and three trained ANNs, corresponding to three levels, were the results of this step.

 Support vector regression
The basic idea of support vector machines (SVM) for regression is to map the data x into a high dimensional feature space via nonlinear mapping and to perform a linear regression in this feature space; more detailed description of SVM for regression can be found in Ref. [37].Radial basis function(RBF) kernels were adopted and five-fold cross validation method was used to determine the c and g .
It is important to note that the correlation constructed was the relationship between two worlds of granules, and the wind speed was in granular format, not in quantitative format.

Interpolation and Extrapolation
For each level, the contemporary reference data were input into trained ANNs to calculate the target data, which meant that three interpolating or extrapolating data were provided by ANNs.
Another ANN was trained using the contemporary data of reference and target granules for

• Support vector regression
The basic idea of support vector machines (SVM) for regression is to map the data x into a high dimensional feature space via nonlinear mapping and to perform a linear regression in this feature space; more detailed description of SVM for regression can be found in Ref. [37].Radial basis function(RBF) kernels were adopted and five-fold cross validation method was used to determine the c and g.
It is important to note that the correlation constructed was the relationship between two worlds of granules, and the wind speed was in granular format, not in quantitative format.

Interpolation and Extrapolation
For each level, the contemporary reference data were input into trained ANNs to calculate the target data, which meant that three interpolating or extrapolating data were provided by ANNs.
Another ANN was trained using the contemporary data of reference and target granules for synthesizing the three results provided by the ANNs, and the topography of the synthesized ANN is shown in Figure 11.The Levenberg-Marquardt (LM) algorithm was chosen for this training.
The final results were provided by the synthesized ANN, and the wind data acquired in this step were in granular format.To transfer the granular format wind data to quantitative data, the Y condition cloud generator were used, and finally, the quantitative results were obtained.
dimensional feature space via nonlinear mapping and to perform a linear regression in this feature space; more detailed description of SVM for regression can be found in Ref. [37].Radial basis function(RBF) kernels were adopted and five-fold cross validation method was used to determine the c and g .
It is important to note that the correlation constructed was the relationship between two worlds of granules, and the wind speed was in granular format, not in quantitative format.

Interpolation and Extrapolation
For each level, the contemporary reference data were input into trained ANNs to calculate the target data, which meant that three interpolating or extrapolating data were provided by ANNs.
Another ANN was trained using the contemporary data of reference and target granules for synthesizing the three results provided by the ANNs, and the topography of the synthesized ANN is shown in Figure 11.The Levenberg-Marquardt (LM) algorithm was chosen for this training.
The final results were provided by the synthesized ANN, and the wind data acquired in this step were in granular format.To transfer the granular format wind data to quantitative data, the Y condition cloud generator were used, and finally, the quantitative results were obtained.

Case Study Description
To implement and validate the granular computing method, a missing wind data interpolation case was utilized.Two meteorological sites located in Hubei province, Central Southern China were chosen; the longitude and latitude, start and end date, equipped instrument, elevation, and local terrain of two sites are shown in Table 5.The two sites are both located in ridgeline with no obstacles around, and the distance between them is 40 km, evaluation map of region of interest is shown in Figure 12.
Data were collected as 1-h values from 1 January 2012 to 31 December 2012, resulting in a total of 8760 records.The correlation coefficient was 0.733, scattered plot of two sites is shown in Figure 13.
The chosen wind farm belongs to a subtropical monsoon climate, located in a complex terrain, and equipped with widely used anemometer, so, gaps in wind speed records are easily emerged due to icing, lightning stroke, equipment failure, and so on.It is a common practical case which is frequently need wind speed interpolation and extrapolation, so, it is a fine representation and the results obtained in this case study can be extended to alternative cases.
In particular, the ten-fold cross validation method was adopted.All records were randomly split into ten mutually exclusive subsets (D 1 , D 2 ...D 10 ) of equal size, the proposed method was trained and tested ten times, each time t ∈ [1,10], it was trained on D/D t and tested on D t , the evaluation metrics of each time was averaged as final evaluation metrics.results obtained in this case study can be extended to alternative cases.
In particular, the ten-fold cross validation method was adopted.All records were randomly split into ten mutually exclusive subsets ( 1 2 10 , ... D D D ) of equal size, the proposed method was trained and tested ten times, each time [1,10] t  , it was trained on / t D D and tested on t D , the evaluation metrics of each time was averaged as final evaluation metrics.In addition, a commercial pitch-regulated 2000 kW rated power wind turbine (GW-2000) was used for the estimations in this study (Figure 14); the power curve was published by the manufacturers in catalogues in table form.Since the power curve data are a discrete pair of values of wind speed and power output, the cubic spline interpolation method was used to calculate the output power from wind speed.In addition, a commercial pitch-regulated 2000 kW rated power wind turbine (GW-2000) was used for the estimations in this study (Figure 14); the power curve was published by the manufacturers in catalogues in table form.Since the power curve data are a discrete pair of values of wind speed and power output, the cubic spline interpolation method was used to calculate the output power from wind speed.In addition, a commercial pitch-regulated 2000 kW rated power wind turbine (GW-2000) was used for the estimations in this study (Figure 14); the power curve was published by the manufacturers in catalogues in table form.Since the power curve data are a discrete pair of values of wind speed and power output, the cubic spline interpolation method was used to calculate the output power from wind speed.It is also important to note that in order to increase accuracy, site analysis was carried out in modal terms with respect to wind direction by ideally subdividing the wind data into 16 wind sectors.

EvaluationMetrics Used
MRE (mean relative error), RMSE (root mean square error), r (correlation coefficient), and MREEP (mean relative error of energy production) were selected as the evaluation metrics using the following expressions.It is also important to note that in order to increase accuracy, site analysis was carried out in modal terms with respect to wind direction by ideally subdividing the wind data into 16 wind sectors.

EvaluationMetrics Used
MRE (mean relative error), RMSE (root mean square error), r (correlation coefficient), and MREEP (mean relative error of energy production) were selected as the evaluation metrics using the following expressions.
where n is the counts of data; V is the interpolated wind speed; V m is the true wind speed; P is the interpolated energy production, kW • h; and P m is the true energy production, kW • h.

Determining the Granular Hierarchy
In this study, the type of anemometer was identified as NRG #40C.As per the NRG #40C white paper [38], the minimum interval was 0.1 m/s; the classification number was 7.7; the range of collected wind speed was 0-17 m/s; and distance to mast center divided by mast leg distance was 0.25.Based on the data, the average measurement uncertainty calculated following the method provided in Section 3.4 was 1.54 m/s, and the measurement and combined uncertainty versus wind speed are shown in Figure 15.

Determining the Granular Hierarchy
In this study, the type of anemometer was identified as NRG #40C.As per the NRG #40C white paper [38], the minimum interval was 0.1 m/s; the classification number was 7.7; the range of collected wind speed was 0-17 m/s; and distance to mast center divided by mast leg distance was 0.25.Based on the data, the average measurement uncertainty calculated following the method provided in Section 3.4 was 1.54 m/s, and the measurement and combined uncertainty versus wind speed are shown in Figure 15.The bottom level was constructed by basic granules, which were the minimum interval of wind speed magnitude recording by anemometer.In this case, the range of collected wind speed was 0-17 m/s, and the magnitude of granularity was 0.1 m/s, therefore the number of clusters was equal to 17 m/s divided by 0.1 m/s, it was 170.
In the middle level, only uncertainties of measurement were considered, and the average measurement standard uncertainty, 0.5 m/s, was adopted for the magnitude of granularity, and the range of collected wind speed was 0-17 m/s, so the number of clusters was equal to 17 m/s divided by 0.5 m/s, it was 34.The bottom level was constructed by basic granules, which were the minimum interval of wind speed magnitude recording by anemometer.In this case, the range of collected wind speed was 0-17 m/s, and the magnitude of granularity was 0.1 m/s, therefore the number of clusters was equal to 17 m/s divided by 0.1 m/s, it was 170.
In the middle level, only uncertainties of measurement were considered, and the average measurement standard uncertainty, 0.5 m/s, was adopted for the magnitude of granularity, and the range of collected wind speed was 0-17 m/s, so the number of clusters was equal to 17 m/s divided by 0.5 m/s, it was 34.
At the top level, all uncertainties were considered, and the average combined uncertainty, 1.54 m/s, was adopted for the magnitude of granularity, and the range of collected wind speed was 0-17 m/s, so the number of clusters was equal to 17 m/s divided by 1.54 m/s, it was 11.

Results Analysis
Missing data in this case study were interpolated based on the granular computing method, and a comparison of the results with traditional methods such as linear regression method, variance ratio method, ANN method, and SVR method was proceeded, results are shown in Figure 16 and Table 6.
It can be observed from the results that GrC with SVR method provided a slight better estimation than GrC with ANN method.To be specific, comparing to GrC with ANN method, the correlation coefficient of the GrC with SVR increase 0.51%, and the MRE, MREEP, and RMSE decrease 6.25%, 0.28% and 2.80%, respectively.
It can be likewise observed from the results that all evaluation indexes of GrC with SVR method provided in this paper are superior to the other methods.The linear regression method, widely used in practical engineering, was chosen as the benchmark, and the correlation coefficient of the GrC with SVR method increase 8.48%, and the MRE, MREEP, and RMSE decrease 61.04%, 60.98% and 15.76%, respectively.It can be concluded from Figure 16 that wind speeds predicted by the granular computing method are very close to the true value in tendency, but prediction error increases in circumstances with small or big wind speed at the target site.
The results of the case study show that the accuracy of the GrC with SVR method meet the needs of practical engineering, and has an advantage over other methods.Therefore, the granular computing method is a fine alternative method in the wind speed temporal interpolation and extrapolation problem.

Conclusions
In this paper, wind speed mixed uncertainty analysis is implemented, and recommended values and estimation tools for wind speed measurement uncertainty and combined uncertainty are provided, as well as a determination of the uncertainty hierarchy.For the purpose of considering mixed uncertainty of randomness and fuzziness, a concrete granular computing method is constructed and a new MCP method of wind speed temporal interpolation and extrapolation based on granular computing theory is proposed; and a wind speed data temporal interpolation case is provided to implement and validate the new method.Hence, the following conclusions are reached: 1.
By using the MCP method proposed in this paper, mixed uncertainty of wind speed had already been considered, thus, a better estimation of the wind speed is provided compared to other methods selected for comparison.In the case study, almost evaluation metrics of interpolation with the proposed method were superior to other methods used in comparison.
In comparison to the linear method, the correlation coefficient of the proposed method increased 8.48%, and the MRE, MREEP, and RMSE decreased 61.04%, 60.98% and 15.76%, respectively.The proposed method improved wind resource assessment accuracy and reduced the risks of wind farm construction.2.
Suitability of using granular computing methods for the issue of wind speed data interpolation and extrapolation is proved.By using GrC method, wind speed mixed uncertainty can be taken into account; accurate results and low cost solutions can be derived.

3.
Mixed uncertainty of wind speed can be divided into three levels, and recommended values of granularity are minimum interval of records, 0.3-0.8m/s, and 1-3 m/s, respectively.Also, estimation tools for wind speed measurement uncertainty and combined uncertainty are provided.

Figure 1 .
Figure 1.Operational standard uncertainty of five anemometers: (a) in class A operational ranges; (b) in class B operational ranges.

Figure 1 . 21 Figure 1
Figure 1.Operational standard uncertainty of five anemometers: (a) in class A operational ranges; (b) in class B operational ranges.

Figure 2 .
Figure 2. Center-line relative wind speed of lattice mast with average T C .
mast centre divided by mast leg distance Center-line relative wind speed Square cross-section and shape edge members Square cross-section and round members Triangular cross-section and round members

Figure 2 .
Figure 2. Center-line relative wind speed of lattice mast with average C T .

Figure 3 .
Figure 3. Average relative calibration uncertainty of ten commonly used cup anemometers with best linear regression fit.

Figure 3 .
Figure 3. Average relative calibration uncertainty of ten commonly used cup anemometers with best linear regression fit.

Figure 4 .
Figure 4. Standard uncertainties of three operational condition classes.

Figure 4 .
Figure 4. Standard uncertainties of three operational condition classes.

Energies 2017, 10 , 1231 10 of 22 Figure 5 .
Figure 5.A schematic view of wind speed data interpolation and extrapolation based on the granular computing method.

Figure 5 .
Figure 5.A schematic view of wind speed data interpolation and extrapolation based on the granular computing method.

Figure 6 .
Figure 6.Flow chart of genetic algorithm procedure.

Figure 6 .
Figure 6.Flow chart of genetic algorithm procedure.
produces a random number with normal distribution whose expectation is En and standard deviation is He .

Figure 9 .
Figure 9. Flow chart of forward cloud generator.

Figure 9 .
Figure 9. Flow chart of forward cloud generator.

Figure 10 .
Figure 10.Topological structure of predicted Artificial Neural Network (ANN) in single granular level.

Figure 11 .
Figure 11.Topological structure of the synthesized ANN.Figure 11.Topological structure of the synthesized ANN.

Figure 11 .
Figure 11.Topological structure of the synthesized ANN.Figure 11.Topological structure of the synthesized ANN.

Figure 12 .
Figure 12.Evaluation map of region of interest.Solid dot show the target and reference site.Figure 12. Evaluation map of region of interest.Solid dot show the target and reference site.Energies 2017, 10, 1231 16 of 22

Figure 13 .
Figure 13.Scattered plot of wind speed of target site versus reference site with best fit line from 1 January 2012 to 31 December 2012. 2000

Figure 13 .
Figure 13.Scattered plot of wind speed of target site versus reference site with best fit line from 1 January 2012 to 31 December 2012.

Figure 13 .
Figure 13.Scattered plot of wind speed of target site versus reference site with best fit line from 1 January 2012 to 31 December 2012.

Figure 14 .
Figure 14.Output power curve of the wind turbine used in the case study.

Figure 14 .
Figure 14.Output power curve of the wind turbine used in the case study.

Figure 15 .
Figure 15.Combined and measurement uncertainty in this case.

Figure 15 .
Figure 15.Combined and measurement uncertainty in this case.

Table 1 .
Class A and Class B operational ranges.

Table 2 .
Classification number of five popular anemometers.

Table 1 .
Class A and Class B operational ranges.

Table 2 .
Classification number of five popular anemometers.

Table 3 .
Relative calibration uncertainty of ten commonly used anemometers.

Table 4 .
Topography, distance, and correlation coefficient of the classes.

Table 4 .
Topography, distance, and correlation coefficient of the classes.
T , or cloud for short.If the distribution of is normal, it is named the normal cloud model, which is a random number set that obeys the normal distributive rule and has a stable tendency.It is determined by expectation Ex, entropy En , and super entropy He , [36]x T C x[36], Ex determines the center of

Table 5 .
Description of the two meteorological sites.

Table 5 .
Description of the two meteorological sites.Evaluation map of region of interest.Solid dot show the target and reference site.