Evaluation Method of Wind Speed Time-Shifting Characteristics at Multiple Scales and Its Application in Wind Power System

(e wind speed sequences at different spatial positions have a certain spatiotemporal coupling relationship. It is of great significance to analyze the clustering effect of the wind farm(s) and reduce the adverse impact of large-scale wind power integration if we can grasp this relationship at multiple scales. At present, the physical method cannot optimize the timeshifting characteristics in real time, and the research scope is concentrated on the wind farm. (e statistical method cannot quantitatively describe the temporal relationship and the speed variation among wind speed sequences at different spatial positions. To solve the above problems, a quantification method of wind speed time-shifting characteristics based on wind process is proposed in this paper. Two evaluation indexes, the delay time and the decay speed, are presented to quantify the time-shifting characteristics. (e effectiveness of the proposed method is verified from the perspective of the correlation between wind speed sequences. (e time-shifting characteristics of wind speed sequences under the wind farms scale and the wind turbines scale are studied, respectively. (e results show that the proposed evaluation method can effectively achieve the quantitative analysis of time-shifting and could improve the results continuously according to the actual wind conditions. Besides, it is suitable for any spatial scale.(e calculation results can be directly applied to the wind power system to help obtain the more accurate output of the wind farm.


Introduction
A high proportion of renewable energy is one of the basic characteristics of the future electric power system, and wind power is an important component of renewable energy [1]. With the rapid expansion of wind power integration, the volatility of wind power has a more adverse impact on the grid dispatching, safe and stable operation of the electric power system, and its power quality [2]. In fact, there is a certain spatiotemporal coupling relationship among the wind speed sequences at different spatial points [3]. If we can grasp this coupling relationship at multiple temporal-spatial scales, it is of great significance to analyze the clustering effect [4], deal with the ramp events and extreme events, obtain the real output of wind farm(s) [5], and so forth, which can effectively reduce the negative effect of the wind power integration [6].
At present, the research algorithms on the spatiotemporal coupling relationship among wind speed time series can be divided into two categories, the model-based physical method and the data-based statistical method, according to the driving mode. e physical method takes the wind speed data at a certain point as input and considers factors such as terrain, roughness, altitude, and flow to establish the aerodynamic model to describe the coupling relationship among the wind speeds at different spatial locations. Zheng [7] assumed that the wind spreads at the speed of the starting point and established the motion equation of wind speed to obtain the average delay time between the wind speed time series at two spatial positions. She et al. [8] used the rotor diameter and geographic coordinate matrix of each wind turbine in the wind farm as the inherent characteristics and took the wind speed and direction of the wind tower, rotor speed, and pitch angle of each wind turbine as inputs to establish the wake effect model and the time-lag effect model; then the dynamic spatial relationship among wind speed time series at each wind turbine location was obtained. Zeng et al. [9] first developed the continuous partial differential equation of each wind turbine by calculating the wake effects, according to the specific spatial locations and distributions of correlated wind turbines. en, the finite volume method is used to differentiate the differential equation at each grid point and derive the spatial correlation matrix of wind speed. Finally, the wind speed at each wind turbine location is obtained by solving the above differential equation through the given boundary conditions. Ye et al. [10] proposed a novel spatial correlation matrix by using computational fluid dynamics to mathematically describe the correlations among wind speeds at different wind turbine locations. e physical method studies the temporal and spatial distribution and the coupling relationship of wind speeds by establishing the aerodynamic model among the wind speeds at different spatial positions. However, this aerodynamic model cannot be optimized in real time according to the actual wind conditions. Besides, when the spatial scale increases, the complexity of the constructed model increases sharply, which leads to the increasing difficulty of modeling and solving the model simultaneously. But it is difficult to ensure the accuracy of the aerodynamic model if the factors affecting wind speed are simplified. erefore, the application of this method is mainly concentrated in the interior of the wind farm, which is difficult to expand to a larger spatial scale. e statistical method uses the actual wind data at different spatial locations to analyze the spatiotemporal coupling relationship among the wind speeds. e research methods can be mainly divided into the following three types: (1) e spatiotemporal coupling relationship is studied by calculating the distribution of wind speeds at different spatial locations, from the perspectives of average wind speed, wind speed frequency distribution, wind shear exponent, and so forth [11]. Liu et al. [12] analyzed the spatiotemporal distributions of the wind resource at the hub height in China from the aspects of average wind speed and average wind power density, based on the large network of wind speed observations, including 2430 meteorological stations from 2006 to 2015. e results show that northeast China has the highest wind speed, while South China has the lowest wind speed. Li and Yu [13] conducted the statistical analysis of the onshore, nearshore, and offshore wind energy conditions from the aspects of Weibull shape and scale factors, average turbulence intensity, and average wind power density. e results indicate that the average wind speed at the offshore location is 1-2 times higher than that nearby onshore area at the height of 50 m above ground level. Bosch et al. [14] estimated the distributions of global offshore wind energy potential from the aspects of average wind speed, mean capacity factor which could be built across the viable offshore area, and generation potential per year. e results show that the wind speed increases with the distance to the coast.
(2) e spatiotemporal coupling relationship is studied by establishing the mapping relationship among wind speed time series at different spatial positions through interpolation, neural network, principal component analysis, and so forth [15]. Philippopoulos and Deligiorgi [16] used two feedforward neural network architectures to establish the relationship among hourly wind speeds at six meteorological sites. An insight into the underlying inputoutput function approximation of the neural networks is obtained by examining their ability to incorporate the mean wind variability characteristics. Miranda and Dunn [17] proposed a spatially correlated wind speed model to characterize the wind resources between various sites by using the multivariate time-series approach with actual wind speed data in different zones. Park and Jang [18] presented a modified Korean parameter-elevation regression on independent slopes model for spatial interpolation of monthly wind speeds. e results show that the proposed model can perform more appropriately to represent the phenomena where similar wind speeds appear continuously along ridges and coastlines, compared with the cokriging model.
(3) e correlation coefficient is used to study the correlation characteristics of wind speed time series at different spatial positions, which is an important part of the spatiotemporal coupling relationship. e calculation method can be divided into the following categories. One is to directly calculate the correlation coefficient based on the wind speed time series [19]. Another is to use the Copula function to fit the joint probability distribution of the wind speed sequences and then calculate the correlation coefficient based on the fitted Copula function [20]. e commonly used correlation coefficients include the Pearson linear correlation coefficient [21], Kendall rank correlation coefficient [22], and cross-correlation coefficient [23]. e statistical method studies the spatiotemporal coupling relationship of wind speed sequences at different spatial positions based on the actual wind data. e obtained results of this method can be continuously improved according to actual conditions, and the research method can be applied to any spatial scale. However, only the correlation relationship between wind speed sequences is quantified in this method. Compared with the physical method, the temporal relationship and the law of speed variation are both implicit in the spatial distribution, and the mapping relationship is not quantified yet.
In conclusion, when the physical method is used to study the spatiotemporal coupling relationship among wind speeds at different spatial positions, the established aerodynamic model cannot be optimized in real time according to the actual wind conditions, and the application of this method mainly concentrated in the interior of the wind farm; it is difficult to expand to a larger spatial scale. When the statistical method is applied, only the correlation relationship of wind speed sequences is quantified, and the temporal relationship and the law of speed variation among wind speed sequences at different spatial positions have not been quantified yet. In order to solve the above problems, a quantification method of time-shifting characteristics based on wind process is proposed in this paper. e results of the proposed method can be continuously improved according to the actual wind conditions and could be applied to any spatial scale. e main contributions of this paper are as follows: (1) e concept of wind process is introduced to study the time-shifting characteristics. In order to cope with the complexity and variability of wind conditions, the wind sequences at different spatial positions are divided according to the certain time window at first, and then the time-shifting characteristics are studied based on the wind processes. (2) Two evaluation indexes, the delay time and the decay speed, are proposed to quantify the time-shifting characteristics of the temporal relationship and the speed variation in each wind process. (3) e time-shifting characteristics are studied at two spatial scales, the wind farms scale and the wind turbines scale. (4) e improved clustering method of k-means, xmeans, is employed to extract the typical wind processes to explore the time-shifting characteristics among wind speeds under various typical wind processes. e obtained results can be directly used in the wind power system. (5) e validity and practicability of the proposed method are verified by applying the calculation results of time-shifting among wind speed sequences at different spatial positions to the wind power system. e remainder of this paper is organized as follows. Section 2 describes the research process of wind speed timeshifting. Section 3 studies the wind speed time-shifting characteristics under two spatial scales, the wind farms scale and the wind turbines scale. Section 4 verifies the effectiveness and practicability of the proposed time-shifting evaluation method. Section 5 concludes this paper.

Research Process of Wind Speed Time-Shifting
is section proposes an evaluation method of the timeshifting characteristics among wind speed sequences at different spatial positions. e research process is as follows.

Outline of the Research Process on Wind Speed Time-
Shifting. First, the wind sequence at each point is preprocessed. e wind sequence is divided according to a certain time window to obtain the wind processes at each point. en some wind processes are screened out to compose a new wind sequence that can be used to study the time-shifting characteristics by judging whether the main wind direction in the wind process is consistent with the relative position of points. For example, the wind sequence at point p can be expressed as a composed sequence of wind process α p , wind process β p , wind process m p , and so forth. en, two evaluation indexes, the delay time (Δt) and the decay speed (Δv), are proposed to study the time-shifting characteristics between wind processes at two points. e firefly optimization algorithm is used to calculate the proposed indexes, in which the constraint conditions are the distance between two points and the actual wind conditions; the optimal objective is the minimum Euclidean distance between the corresponding wind processes. Taking the wind process α p and the wind process α q as an example, the wind process α p at point p and the wind process α q at point q can be expressed as v α respectively. e distance between two points and the wind conditions of the wind process α p at point p are set as the constraint conditions, and the wind process α q is offset through the time axis and the speed axis to achieve the objective of the minimum Euclidean distance between the wind process α p at point p and the new wind process α q ′ at point q. en the corresponding delay time Δt α , decay speed Δv α , and the new wind process Finally, the time-shifting characteristics of wind speed sequences at two points are studied based on the proposed indexes. e delay time and the decay speed of each wind process at point p and point q are calculated, and then the statistical analysis method is used to analyze the timeshifting. Furthermore, the improved clustering method of kmeans, x-means, is employed to extract the typical wind processes, and the time-shifting characteristics of different typical wind processes are studied so that the research results can be applied to the actual planning and operation of the electric power system.
In order to verify the effectiveness of the proposed timeshifting analysis method, the wind speed sequence at point q is reconstructed after considering the time-shifting characteristics. t-Copula function is used to fit the joint probability distribution between wind speed sequences under two conditions, that is, whether the time-shifting is considered. en the Kendall rank correlation coefficient τ, the Spearman rank correlation coefficient ρ s , and the tail correlation coefficient λ are calculated based on the t-Copula function.
e correlation coefficients increase with the increasing similarity between two wind speed sequences. e outline of the research process on wind speed timeshifting is shown in Figure 1.
It should be noted that, due to the fact that the wind speed can gradually recover within a certain distance after passing through the wind turbine(s) or the wind farm(s), the evaluation indexes for time-shifting analysis should be different according to the distance and the installed capacity of wind power, that is, whether the decay speed should be Mathematical Problems in Engineering 3 contained when studying the time-shifting characteristics of wind speed sequences at different spatial positions. According to the research results of Risø Laboratory at the Technical University of Denmark, even for a large-scale wind farm of 30 km × 30 km at the upwind direction, the weakening effect on the wind speed at the downwind direction could vanish after a distance of 30 km to 60 km [24]. erefore, it can be considered that there is no influence of the decay speed among wind speed sequences at different locations when the distance between the wind farms exceeds 60 km.

Calculation Method of the Proposed Indexes Based on Firefly Optimization
Algorithm. e firefly optimization algorithm [25] is applied to calculate the proposed indexes, the delay time and decay speed. Each firefly is regarded as a point in the algorithm, and the fireflies are distributed randomly in the search space at first.
ere are two parameters in the firefly optimization algorithm, the brightness and the attractiveness. e brightness of the firefly is decided by the objective function, which indicates the position and decides the next moving direction. e attractiveness between two fireflies is inversely proportional to the distance, which decides the next moving distance. e firefly is attracted and moves toward the relatively brighter one according to the brightness and attractiveness; then the fireflies gather through the movement. In the optimization process, the brightness and attractiveness are updated iteratively until the optimization objective is achieved. e specific steps of the firefly algorithm are as follows.
Step (a): to set the initial parameters of the algorithm, including the maximum attraction β0, the light attraction coefficient c, the step factor α, the maximum number of iterations, and the number of fireflies.
Step (b): to initialize the position of fireflies.
Step (c): to calculate the objective function of each firefly and regard it as the maximum fluorescent brightness I 0 .
Step (d): to define the relative brightness between fireflies, as shown in the following equation: where r ij is the relative distance between firefly i and firefly j and the light attraction coefficient c reflects the weakening characteristics with the increase of the distance and the absorption of the propagation medium. to define the attractiveness between fireflies, as shown in the following equation: to calculate the relative brightness and attractiveness between fireflies according to equations (1) and (2) and sort them by the relative brightness. To propose evaluation indexes Taking wind process α p and wind process α q as an example To offset αp through the time axis and the speed axis

Mathematical Problems in Engineering
Step (e): to update the position of fireflies, as shown in the following equation: where x i and x j are the positions of firefly i and firefly j, respectively. e step factor α is a disturbance term, which is added to enlarge the search area and avoid the optimization algorithm falling into partial convergence prematurely. In general, α ∈ [0, 1] and ε i ∈ [− 0.5, 0.5].
Step (f ): to repeat Step (c) to Step (e) until the maximum iterations number is reached, and then the optimal results can be obtained.

Correlation Analysis Based on t-Copula
Function. e correlation analysis between wind speed time series at different spatial points is an important method to test the validity of the proposed time-shifting indexes and evaluate the effect of the time-shifting characteristics. At present, most researchers use the linear correlation coefficient to measure the correlation relationship between wind speed sequences. However, when the distribution of the variable follows the nonnormal distribution or nonelliptic distribution, it may get wrong conclusions [26]. e distribution of wind speed is more consistent with Weibull distribution, so it is more reasonable to use the rank correlation coefficient to measure the correlation relationship between wind speed sequences. Furthermore, compared with the normal weather conditions, the wind power output in the extreme wind conditions has a greater impact on the electric power system, and the output correlation among wind turbines or wind farms should be paid more attention in these weather conditions. t-Copula [27] function is one of the commonly used Copula functions. It has a thicker tail and is more sensitive to the changes of the tail correlation between random variables, which can capture the symmetric tail correlation. erefore, the t-Copula function is used to fit the joint probability density of wind speed sequences under two conditions in this paper, including the original sequences and the sequences considering time-shifting characteristics. en, the Kendall rank correlation coefficient, Spearman rank correlation coefficient, and tail correction coefficient can be calculated to verify the validity of the proposed time-shifting indexes and evaluate the effect of the time-shifting characteristics. e specific steps of the correlation analysis based on t-Copula function are as follows.
e binary t-Copula function is shown in equation (4), where the linear correlation coefficient between two variables is ρ and the degree of freedom is k. e density function and the distribution function of t-Copula are depicted in Figure 2: where t − 1 k is the inverse function of the distribution function of unary t distribution with freedom degree of k. e calculation equations of Kendall rank correlation coefficient τ, Spearman rank correlation coefficient ρ s , and tail correlation coefficient λ based on t-Copula function are shown in the following equations: where t k+1 ((

Extraction of Typical Wind Processes Based on the x-Means
Clustering Method. e improved algorithm of k-means algorithm, x-means [28], is adopted to extract the typical wind processes when analyzing the time-shifting characteristics of each typical wind process. Compared with kmeans algorithm, x-means algorithm is not required to give the specific value of the clustering number n in advance but only needs to give a range within which the true n reasonably lies. e clustering result is not only the set of centroids but also the value for n which scores best in this range. e specific steps of x-means are as follows.
Step (a): to input m wind processes to be clustered and set the range of clustering number as [n min , n max ], where n max ≤ m.
Step (b): to select n min initial clustering centers from original wind processes randomly.
Step (c): to calculate the Euclidean distance between other wind processes and each clustering center and assign the corresponding wind process to the nearest category.
Step (d): to calculate the average value of wind processes in each clustering set and use the results as the new clustering centers.
Step (e): to repeat Step (c) and Step (d) until the clustering centers no longer change, and the optimal clustering results can be obtained with the clustering number of n min .

Mathematical Problems in Engineering
Step (b) to Step (e) are the main steps of k-means algorithm when the clustering number is n min .
Step (f ): to apply the 2-means algorithm in each clustering set, and the clustering results are evaluated according to the Bayesian Information Criterion (BIC). e calculation method of BIC can be seen in [28]. If the BIC score is larger when 2-means algorithm is used in the clustering set, the set will be divided into two sets; otherwise, the clustering set will be retained.
Step (g): to repeat Step (f ) in each new clustering set until the clustering number n is no longer increased or n ≥ n max .

Case Study and Analysis
is section analyzes the time-shifting characteristics of wind speed time series in different points under the wind farms scale and the wind turbines scale, respectively. e data includes wind speed and wind direction. In the wind farms scale, the data is the mesoscale data located in South Gobi of Mongolian; the time range is from January 1, 2004, to October 31, 2018, and the temporal resolution is 1 h. All points are located in the flat terrain in this case. In the wind turbines scale, the data is the real operation data of one wind farm in Northern China; the time range is from January 1, 2014, to December 30, 2014, and the temporal resolution is 1 min. All points are located in the hilly terrain in this case.

Time-Shifting under Wind Farms Scale.
e wind sequence at each point is preprocessed at first. e wind sequence is divided into different wind processes with a time window of 1 day in this case, and then the wind process sets are obtained at each point. e relative distribution of each point is shown in Figure 3, so the wind processes with the main wind direction between 0265°and 275°at the corresponding point should be screened out for the time-shifting study. For example, when studying the time-shifting characteristics of wind speed sequences from point 1 to other points (2)(3)(4)(5), the wind processes with the main direction between 265°and 275°at point 1 and the contemporaneous wind processes at other points should be screened out simultaneously. en, the new wind speed sequence that can be used to study the time-shifting characteristics at each point is reconstructed. e average speeds of new wind speed sequences at points 1 to 4 are 8.32 m/s, 8.44 m/s, 6.68 m/s, and 6.06 m/s, respectively. en, the time-shifting characteristics of new wind speed sequences between point 1 and points 2 to 5, between point 2 and points 3 to 5, between point 3 and points 4 to 5, and between point 4 and point 5 are studied, respectively. Since the distances between two points are all more than 100 km, only the delay time is used as the time-shifting evaluation index in this case. e delay time distributions of wind speed time series between different points are shown in Figure 4. It can be seen that, with the increase of distance, the median delay time between wind speed sequences at two points increases, and the frequency distribution becomes more dispersed. e delay time of the wind speed sequence from point 1 to other points (2-5) varies in 0 h-26 h, 0 h-31 h, 0 h-49 h, and 0 h-78 h, respectively. e average delay time between wind speed sequences at two points has a correlation relationship with the average wind speed at the starting point and the distance between two points, but it is not completely equal to the quotient between these two variables. Take the wind speed sequence at point 1 as an example; if the distances between point 1 and points 2 to 5 are divided by the average wind speed at point 1, the corresponding results are 3 h, 7 h, 10 h, and 13 h, respectively. When the proposed time-shifting analysis method is applied, the average delay time of wind speed sequence from point 1 to other points (2-5) is 2 h, 5 h, 8 h, and 14 h, respectively.
is is because even if the average speed of two wind processes is the same, the wind conditions can also be quite  different, and the time-shifting characteristics are various under different wind processes. erefore, it is more reasonable to calculate the delay time between wind speed sequences at two points based on the similarity of wind processes.
In order to preliminarily analyze the mapping relationship among the time-shifting characteristics of wind speed sequences at different points, the distances between points, and the actual wind conditions, the wind processes at point 1 are divided according to the average wind speed of each process, with the interval of 1 m/s. e average delay time of point 1 to points 2-5 is calculated in each wind speed bin, respectively, and the results are depicted in Figure 5. As can be seen, the delay time increases with the distance when the wind speed is almost the same and decreases with the wind speed when the distance is basically the same. However, since the wind speed is inevitably affected by various factors such as terrain and obstacles when it propagates from one point to another point in space, the propagation time is not equal to the quotient of distance and its value. Furthermore, the wind process is the time series composed of different wind speeds, so the propagation law of the wind process is more complicated.
is is why some dots in Figure 5 look like the "abnormal dots," such as the average delay time from point 1 to point 2 when the average speed of the wind process is between 8 m/s and 9 m/s and the average delay time from point 1 to point 3 when the average speed of the wind process is between 11 m/s and 12 m/s. erefore, it is necessary to further analyze the time-shifting characteristics of different wind processes at different spatial scales.
x-Means algorithm is applied to cluster the wind processes at point 1 to extract typical wind processes, and then the typical wind process sets are constructed. e range of clustering number is set as in [3,15] when using x-means algorithm. e optimal clustering number is 24 in this case. e average delay time of each wind process set from point 1 to other points (2-5) is calculated, respectively, as shown in Figure 6. Six typical wind process sets are selected randomly, as depicted in Figure 7. e red solid line is the clustering center and the black-dashed lines are the actual wind processes of each set. e average delay time from point 1 to other points and the average wind speed of these six typical wind process sets are listed in Table 1.
e new wind speed sequence at each point considering the time-shifting characteristics is reconstructed, respectively, according to the calculation results of the delay time. e wind speed time series of whether considering the timeshifting is shown in Figure 8. As can be seen, the wind speed sequences at different points could be more similar if the time-shifting characteristics are considered.
In order to further verify the validity of the proposed method, the correlation analysis is used to test the effectiveness of the time-shifting characteristics. t-Copula function is used to fit the original wind speed sequences of two points and the new wind speed sequences considering the time-shifting characteristics of two points, respectively. en the Kendall correlation coefficient τ, Spearman correlation coefficient ρ s , and tail correlation coefficient λ are

1#-2#
1#   respectively, when the Kendall correlation coefficient τ is used as the evaluation index; the correlations are increased by 0.046, 0.107, 0.149, and 0.288, respectively, when the Spearman correlation coefficient ρ s is used as the evaluation index. Furthermore, the tail correlation of wind speed sequences is also improved when the time-shifting characteristics are taken into account, which means the appearance times under high wind conditions and small wind conditions are more consistent. e tail correlation coefficients of wind speed sequences between point 1 and points 2 to 5 are increased by 0.090, 0.050, 0.025, and 0.045, respectively.

Time-Shifting under Wind Turbines Scale Containing
Wake Effect. e wind sequence at each point is preprocessed at first. e wind sequence is divided into different wind processes with a time window of 1 hour in this case, and then the wind process sets are obtained at each point. e relative distribution of each point is shown in Figure 9, so the wind processes with the main wind direction between 130°and 140°at the corresponding point should be screened out for the time-shifting study. e average speeds of new wind speed sequences at points 1 to 3 are 11.45 m/s, 10.17 m/ s, and 10.49 m/s, respectively. en, the time-shifting characteristics of new wind speed sequences between point 1 and points 2 to 4, between point 2 and points 3 to 4, and between point 3 and point 4 are studied, respectively. Since the data used in this case is the actual operation data of the wind farm, the factors that affect wind speed sequences at different locations should include both the delay time and the decay speed. e distributions of delay time and decay speed of wind speed time series between different points are shown in Figure 10. It can be seen that, with the increase of distance, the median delay time between wind speed sequences at two points increases, and the frequency distribution becomes more dispersed. e wind speed at the downstream location is not always lower than that at the upstream location. e average decay speeds of the wind speed sequence from point 1 to other points is is because the points are all located in the hilly terrain in this case, and the propagation of wind speed between different points is greatly affected by the terrain. Furthermore, the wake effect of wind turbines is weakened with the increase of the distance and the relief of the terrain.
In order to preliminarily analyze the mapping relationship among the time-shifting characteristics of wind speed sequences at different points, the distances between points, and the actual wind conditions, the wind processes at point 1 are divided according to the average wind speed of each process, with the interval of 1 m/s. e average delay time and average decay speed of wind processes from point 1 to points 2-4 are calculated in each wind speed bin, respectively, and the results are depicted in Figure 11. As can be seen, the delay time increases with the distance when the wind speed is almost the same and decreases with the wind speed when the distance is basically the same. However, due to the influence of the small distance between two points and the low time resolution of the data, the delay time of wind speed sequence between two points does not decrease with the increase of wind speed in some situations. is situation happens more frequently with the shortening of the distance between two points. e decay speed of wind speed sequences between different points varies with the value of wind speed. e improved algorithm of k-means, x-means, is employed to cluster the wind processes at point 1 to extract typical wind processes, and then the typical wind process sets are constructed. e range of clustering number is set as in [4,20], and the optimal clustering number is 32 in this case. e average delay time and the average decay speed of each wind process set from point 1 to other points (2)(3)(4) are calculated, respectively, and the results are shown in Figure 12. Six typical wind process sets are selected randomly, as depicted in Figure 13. e average delay time and the average decay speed from point 1 to other points and the average wind speeds of these six typical wind process sets are listed in Table 3.
e new wind speed sequence at each point considering the time-shifting characteristics is reconstructed, respectively, according to the calculation results of the delay time and decay speed. e wind speed time series of whether considering the time-shifting is shown in Figure 14. As can be seen, the wind speed sequences at different points could be more similar if the time-shifting characteristics are considered.
In order to further verify the validity of the proposed method, the correlation analysis is used to test the effectiveness of the time-shifting characteristics. t-Copula function is used to fit the original wind speed sequences of two points and the new wind speed sequences considering the timeshifting characteristics of two points, respectively. en the Kendall correlation coefficient τ, Spearman correlation coefficient ρ s , and tail correlation coefficient λ are calculated based on the fitted function, and the results are listed in Table 4. It can be seen that the correlation of wind speed sequences at different points is enhanced significantly when the time-shifting characteristics are considered. For example, taking point 1 as the starting point, if the time-shifting characteristics are taken into account, the correlations of wind speed time series between point 1 and points 2 to 4 are increased by 0.142, 0.103, and 0.202, respectively, when the Kendall correlation coefficient τ is used as the evaluation index; the correlations are increased by 0.065, 0.046, and 0.124, respectively, when the Spearman correlation coefficient ρ s is used as the evaluation index. Furthermore, the tail correlations of wind speed sequences between point 1 and points 2 to 4 are increased by 0.374, 0.147, and 0.802, respectively. It should be noted that, different from the case in wind farms scale, the points are all located in the hilly terrain in this case, the mapping relationship among the wind speed sequences at different points is more complicated, and the time-shifting effect does not show a more obvious trend with the increase of distance.

Verification of the Proposed Time-Shifting Evaluation Method
is section takes the data on December 31, 2014, under the wind turbines scale as an example to verify the calculation results of time-shifting and illustrate the effectiveness of the  point 1 is between 130°and 140°. If so, it should be judged at first which wind process set the wind process belongs to, and then the wind process at points 2 to 4 should be constructed according to the above calculation results of time-shifting. If not, the corresponding wind speed at points 2 to 4 can be considered the same as that at point 1. e actual wind speed sequence and wind speed sequences of whether considering time-shifting at points 2 to 4 are shown in Figure 15. It can be seen that the wind speed sequence is more similar to the actual sequence at each point if the time-shifting characteristics are considered.
To further verify the calculation results of timeshifting, the average absolute deviations between the calculated wind speed sequences of whether considering the time-shifting and the actual wind speed sequence are used as the evaluation index to quantify the effectiveness of time-shifting. e average absolute deviations between the calculated wind speed sequence and the actual wind speed sequence at points 2 to 4 are 1.27 m/s, 1.18 m/s, and 1.26 m/s, respectively, if the time-shifting is not taken into account. If the time-shifting is considered, the average absolute deviations are 1.06 m/s, 1.02 m/s, and 1.04 m/s, respectively. e average absolute deviations are reduced by 16.5%, 13.6%, and 17.5%, respectively, when the time-shifting is considered. e average absolute deviations between wind speed sequences of whether considering the time-shifting and the actual wind speed sequence in each time window at points 2 to 4 are depicted in Figure 16. It can be seen that the wind speed     sequences are more similar to the actual wind speed sequences in all time windows if the time-shifting is considered.
In order to study the influence of time-shifting on the wind power generation system, the measured wind speed sequence and the calculated wind speed sequences under two conditions (considering time-shifting or not) at each point are converted to the power sequences based on the theoretical power curve of the wind turbine, respectively; then the total output of wind turbines at points 1 to 4 is obtained by adding the output of wind turbine at each point. e total power sequences of four wind turbines calculated by the measured wind speed sequence and the wind speed sequences of whether considering time-shifting are shown in Figure 17. It can be seen that the calculated total power is more similar to the actual power when the time-shifting is considered. To further illustrate the necessity of timeshifting, the influence of time-shifting is quantified from two perspectives: one is the volatility of the total power sequence and the other is the generation of four wind turbines. Fluctuation ratio (FR) [29] is used to evaluate the volatility of the power sequence in this paper. e calculation method of FR is shown in the following equation: where P (t) is the power of wind turbines at time t and P cap is the rated power of wind turbines, 8.02 MW in this study. e total power generation and the volatility of the power sequence of four wind turbines are calculated through the wind speed sequences under two conditions, and the results are compared with the results when using the measured wind speed sequences, as shown in Table 5. As can be seen from the table, if the time-shifting is taken into account, the calculated daily power generation is closer to the actual generation, and the volatility of the calculated power sequence is more similar to that of the actual power sequence. e aggregation characteristics of wind turbines output are reflected. Compared with the actual daily power generation, the deviation of the calculated generation is 8.6% if the time-shifting is considered and is 10.6% if the time-shifting is not considered. Compared with the volatility of the actual power sequence, the deviation of the volatility of the calculated power sequence is 70.2% if the time-shifting is considered and is 77.2% if the time-shifting is not considered.

Conclusions
A quantification method of wind speed time-shifting at different spatial positions based on wind process is proposed in this paper. As a result of this study, the following conclusions can be drawn: (1) e proposed evaluation method can effectively achieve the quantitative analysis of the time-shifting characteristics of wind speed sequences at different spatial positions. Furthermore, the results of the proposed method can be continuously improved according to the actual wind conditions and could be applied to any spatial scale. (2) e time-shifting effect between wind speed sequences increases with the distance when points are located in the flat terrain. But when points are located in the hilly terrain, the time-shifting effect does not show a more obvious trend with the increase of distance. (3) e median delay time of wind speed sequences increases and the frequency distribution of the delay time becomes more dispersed with the increase of distance, no matter what terrain the points are located in. (4) e wind speed at the downstream location could be larger than that at the upstream location in the wind farm when points are located in the hilly terrain. (5) x-Means clustering method can achieve the effective division of wind processes, and the time-shifting characteristics of different wind processes are quite various. (6) e calculated output of wind turbines is more similar to the actual output in the wind power system if the time-shifting characteristics are taken into account.
e time-shifting characteristics of wind speeds are studied when the wind direction is consistent with the relative position of points in this paper. In future research, the proposed evaluation method will be applied to analyze the time-shifting characteristics of wind speeds in full wind directions to better grasp the spatiotemporal coupling relationship of wind speeds at multiple scales.

Data Availability
e data used to support the findings of this paper are available upon request from the corresponding author.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.