Next Article in Journal
Weighted-PSO Applied to Tune Sliding Mode Plus PI Controller Applied to a Boost Converter in a PV System
Next Article in Special Issue
A Flower Pollination Method Based Global Maximum Power Point Tracking Strategy for Point-Absorbing Type Wave Energy Converters
Previous Article in Journal
Experimental Identification of the Optimal Current Vectors for a Permanent-Magnet Synchronous Machine in Wave Energy Converters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of Generalized Interaction Wake Models Systems with ELM Variation for Off-Shore Wind Farms

1
School of Logistics Engineering, Wuhan University of Technology, Wuhan 430063, China
2
National Engineering Research Center for Water Transport Safety (WTS Center), Wuhan 430063, China
3
State Key Laboratory of Ocean Engineering, Shanghai Jiaotong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Energies 2019, 12(5), 863; https://doi.org/10.3390/en12050863
Submission received: 20 January 2019 / Revised: 23 February 2019 / Accepted: 28 February 2019 / Published: 5 March 2019

Abstract

:
This paper reports a novel frandsen generalized wake model and its variation model-frandsen generalized normal distribution wake model for off-shore wind farms. Two different new wake models in off-shore wind farms have been studied comparatively. Their characteristics have been analyzed through mathematical modeling and derivation. Meanwhile, simulation experiments show that the proposed two new wake models have different properties. Furthermore, the distributions of wind speed and wind direction are modeled by the statistical methods and Extreme Learning Machine through the off-shore wind farms of Yangshan Deepwater Harbor in the Port of Shanghai, China. In addition, the data of wind energy are provided to verify and test the correctness and effectiveness of the proposed two models. Wind power has been demonstrated by wind rose and wind resources with real-time data. These techniques contribute to enhance planning, utilization and exploitation for wind power of off-shore wind farms.

1. Introduction

The global warming and the climate changes lead to a gradual shift from conventional to renewable energy sources which are more reliable clean resources. The main reason is the growing trend of global energy demands accompanied with the detrimental effects of overuse conventional fossil fuel sources namely [1,2,3].
Wind is one of the fastest growing energy sources, and it is also pollution-free, renewable and abundant. At present, some researchers have done a lot of research work on wind energy [4,5]. For example, based on short term wind speed forecasting of variable weight, Li et al. provided the research and application of a combined model [6]. Romanic et al. studied wind and tornado climatologies and wind resource modelling for a modern development situated in ‘Tornado Alley’ [7]. Ahmed Shata Ahmed investigated wind energy characteristics and wind park installation in Shark El-Ouinat, Egypt [8]. Using standard exergy and Extended Exergy Accounting (EEA) approaches, Aghbashlo et al. studied performance assessment of a wind power plant [9].
Wind energy has a low carbon footprint, which is a type of renewable energy. Some researchers and experts believe that local power generation sources like micro-grids on wind farm or micro power sources are better, in view of transmission efficiency and grid laying cost [10]. Some of the methods by which sustainable energy can be harvested are solar, wind, vibration, tidal, etc. Ideally, a power source should be sustainable, reliable, economical and eco-friendly [11,12].
Among renewable energy sources, Wind Power (WP) is considered to be one of the most prospective renewable energy technologies, and its usage has been increased immensely over recent decades. Usually, Horizontal Axis Wind Turbines (HAWTs) as types of wind energy converters are a kind of medium to large scale rotating machinery and their technologies have been developed substantially during the last decades, resulting in significant performance improvement. Despite the considerable achievements in the aerodynamic improvement of the HAWTs, wind energy still has certain drawbacks that make the popularity of the wind turbine technology difficult.
Due to the location optimization caused by operation in the unsteady flow condition and dynamic loading exerted on the Wind Turbine (WT), one of its main problems is how to increase performance of the control system as decision-making, which makes the most effective use of the WT controller. The WT location optimization is an interesting topic, which provides a predominant approach to increase the total WP of the Wind Farm (WF) and decreases the wake effect of WTs [13,14,15]. In recent years, the size of WTs has been more extensive from a few kilowatts to several megawatts. Lots of experiences have shown that the larger the WTs, the lower the cost per kilowatt installed [16]. Furthermore, their total costs of production, installation and maintenance are less than the total of smaller WTs achieving the same WP [17]. In 2013, Chen et al. investigated tower height matching optimization for WT positioning in the wind farm [18]. One scheme for properly handling these aerodynamic interactions is to use and promote wake models in the optimization and distributed algorithms control. An alternative approach is to present an online control method where each WT adjusts its own sense model coefficients, which is in line with the communication of local Wind Farms (WFs) [19,20,21]. For example, V. Seshadri Sravan Kumar and D.Thukaram presented accurate modeling of doubly fed induction generator [22] based WFs in load flow analysis. Farajzadeh et al. proposed statistical modeling of the power grid from a wind farm standpoint [23]. Tian et al. developed and verified a new WT wake model with two dimensions in 2015 [24,25,26,27].
It is worth mentioning that Wind Energy (WE) has developed rapidly in China in recent years. For example, Figure 1 and Figure 2 show off-shore wind farms (OSWFs) of Yangshan Deepwater Harbor in the Port of Shanghai, China, which is one of the largest wide range or scope WFs in China [28,29].
From the control point of view, the research of WFs has attracted a great deal of interest from researchers. Recently, Ebrahimi et al. proposed a new optimizing power control scheme based on a centralizing WF control system [30]. Song et al. presented the decision model of WF layout design with three dimensions [31]. Sadegh Ghani Varzaneh et al. have studied a novel simplified model for assessment of power alteration of Doubly-fed Induction Generator (DFIG) [32] based wind farm participating in frequency control system. Hossain has presented a nonlinear controller for transient stability enhancement of DFIG based new bridge type fault current limiter [33] for WFs. Yao et al. studied coordinated control of hybrid WFs [34] in a Permanent Magnet Synchronous Generator (PMSG) based and fixed-speed induction generator (FSIG) for WTs during an asymmetrical grid fault. Li et al. also proposed adaptive fault lenient control of WTs with safeguard transient performance considering active power control of WFs [35]. From met-mast and remote sensing techniques, Chaurasiya et al. proposed comparative analysis of Weibull parameters with wind data measured [36]. Hamid Atighechi et al. presented an effective load shedding remedial action strategy [37] for WF generation. Suganthi et al. proposed an Improved Differential Evolution algorithm [38] based on congestion management in the presence of wind turbine generators.
Recently, there are some machine learning techniques used in such prediction of wind direction and wind speed for OSWFs. For example, Wan et al. proposed an Extreme Learning Machine (ELM) based method of probabilistic forecasting for wind power generation [39]. Vlastimir Nikolić et al. presented a novel wake model based upon ELM for sensor-less computation of wind speed based on WT parameters in WFs [40]. Lazarevska presented an alternative approach to forecasting the wind speed based on Extreme Learning Machine [41]. Wu et al. presented a real-time precise wind speed estimation approach and sensor-less control for variable pitch and variable speed Wind Turbine Power Generation System (WTPGS) [42].
Based on the above discussion, in general, the main contributions of this paper are as follows. Firstly, a Frandsen Generalized Wake Model (FGWM) and its variation model-Frandsen Generalized Normal Distribution Wake Model (FGNDWM) for WFs have been analyzed and presented with mathematical derivation forms. Then, comparisons of these two different wake models of OSWFs have been presented. Furthermore, comparative experiments of both Wake Models have been studied. Finally, focused on the OSWFs of Yangshan Deepwater Harbor in the Port of Shanghai, China, wind rose, wind Weibull probability density distribution and ELM prediction are elaborated and discussed through the OSWFs of Yangshan Harbor. Simulation figures are also provided to show the effectiveness of the proposed approach.
The structure of the paper is as follows: FGWM and FGNDWM are derived and studied in Section 2 and Section 3. Meanwhile, in Section 4 and Section 5, comparative analysis and experiments for the proposed two novel different wake models are studied. Furthermore, in Section 6, the wind rose, wind Weibull probability density distribution and ELM prediction are proposed for OSWFs of Yangshan Deepwater Harbor in the Port of Shanghai, China. Finally, in Section 7, the conclusions are summarized.

2. A Frandsen Generalized Wake Model (FGWM) for OSWFs

Frandsen Generalized Wake Model (FGWM) in the ideal state is shown in Figure 3, where the far wake area is described by the sideways trapezoidal region. The near field is denoted W T 1 i j (radius is r 1 i j ) and can be seen as a turbulent wake. In the down-wind distance x i j , the wind speed denoted P 1 i j and P 2 i j are assumed to be equal v 0 i j , and the wind speed on S x i j is given by v x i j . The circular cross-section radius is r x i j , where i and j are row vector and column vector for Wind Turbine in large-scale Wind Farms, respectively.
In an ideal state, the model assumes that the far wake region spreads with a linear approach, and the distribution of wind speed is homogeneous on every cross-section. In FGWM, the tube includes the near wake area, which is described in the green rectangular region in Figure 3. The fluid inlet mass flow rate in the tube is equal to i = 1 n j = 1 m ρ π r x i j 2 v 0 i j and it goes through the WT. The fluid outlet mass flow rate is equal to i = 1 n j = 1 m ρ π r x i j 2 v x i j . The FGWM assumes
i = 1 n j = 1 m M ˙ i j = i = 1 n j = 1 m M i j t i j = i = 1 n j = 1 m ρ π r x i j 2 v x i j ,
for n = 1 , 2 , , N ; m = 1 , 2 , , N .
By conservation of momentum, we obtain
i = 1 n j = 1 m T i j = i = 1 n j = 1 m M ˙ i j v 0 i M ˙ i j v x i j = i = 1 n j = 1 m M i j t i j v 0 i j M i j t i j v x i j ,
for n = 1 , 2 , , N ; m = 1 , 2 , , N .
We assume that the radius of the row vector ith and column vector jth actuator disk is r r i j , then the area of ith × jth actuator disk is given as A r i j = π r r i j 2 . According to the definition of thrust coefficient C T i j , one can have
i = 1 n j = 1 m C T i j = i = 1 n j = 1 m T T h r u s t F o r c e i j T D y n a m i c F o r c e i j = i = 1 n j = 1 m T i j T i j m a x = i = 1 n j = 1 m 1 2 ρ ( v 0 i j 2 v x i j 2 ) A r i j 1 2 ρ v 0 i j 2 A r i j = i = 1 n j = 1 m 2 ρ v 0 i j 2 A r i j a i j ( 1 a i j ) 1 2 ρ v 0 i j 2 A r i j ,
which is equivalent to
i = 1 n j = 1 m T i j = i = 1 n j = 1 m C T i j T i j m a x .
Substituting Equations (1), (3) and (4) into Equation (2), we obtain the following equation:
i = 1 n j = 1 m v x i j 2 v 0 i v x i j + 1 2 C T i j ( r r i j r x i j ) 2 v 0 i j = 0 .
Solving Equation (5), one can have the following equation:
i = 1 n j = 1 m v x i j = i = 1 n j = 1 m v 0 i j ± v 0 i j 2 2 C T i j ( r r i j r x i j ) 2 v 0 i j 2 2 = i = 1 n j = 1 m 1 2 ± 1 2 1 2 C T i j ( r x i j / r r i ) 2 v 0 i j = i = 1 n j = 1 m 1 1 2 1 ± 1 2 C T i j ( r x i j / r r i j ) 2 v 0 i j .
Using the physical solution of Equation (6), then we obtain
i = 1 n j = 1 m v x i j = i = 1 n j = 1 m 1 1 2 1 1 2 C T i j ( r x i j / r r i j ) 2 v 0 i j .
The above equation gives the main variables and results of the FGWM.

3. A Frandsen Generalized Normal Distribution Wake Model (FGNDWM) for OSWFs

The Frandsen Generalized Normal Distribution Wake Model (FGNDWM) called Frandsen Generalized Gaussian Distribution Wake Model (FGGDWM) for OSWFs is illustrated in Figure 4 where the two dotted lines A i j and C i j are selected to be the boundaries of the FGNDWM tube. The far wake region is confined to the dotted line tube region while the farthest boundary is extended to infinity. The wind speed on P 1 i j and P 2 i j are recognized as v 0 i j and the wind speed on S x i j is v x i j . Here, i and j are row vector and column vector for Wind Turbine in large-scale offshore Wind Farms, respectively.
The conservation of mass does not hold comparing to the FGWM. The following equation is considered as the outlet mass flow rate
i = 1 n j = 1 m M ˙ i j = M i j t = i = 1 n j = 1 m 0 + 2 π r x i j ρ v x i j , r x i j d r x i j , n = 1 , 2 , , N ; m = 1 , 2 , , N .
The FGNDWM satisfies the following equation:
i = 1 n j = 1 m v x i j , r x i j = i = 1 n j = 1 m 1 Ψ ( x i j ) exp r i j 2 2 σ i j 2 v 0 i j ,
where σ i j is the standard deviation, which is also called the characteristic width of FGNDWM and Ψ ( x i j ) is a coefficient related to x i j .
According to the momentum conservation law, one can obtain the following result
i = 1 n j = 1 m 0 + 2 π r x i j ρ v x i j , r x i j v 0 i v x i j , r x i j d r x i j = i = 1 n j = 1 m T i j ,
which is equivalent to
i = 1 n j = 1 m lim b + 0 b 2 π r x i j ρ v x i j , r x i j v 0 i j v x i j , r x i j d r x i j = i = 1 n j = 1 m T i j .
From the definition of the thrust coefficient, one can have the following equation:
i = 1 n j = 1 m T i j = C T i j T i j m a x = 1 2 i = 1 n j = 1 m ρ A r i j v 0 i 2 C T i j = 1 2 i = 1 n j = 1 m ρ π r r i j 2 v 0 i j 2 C T i j .
Substituting Equations (9) and (12) into Equation (11), we obtain
i = 1 n j = 1 m lim b + 0 b 2 π r x i j ρ v x i j , r x i j ( v 0 i v x i j , r x i j ) d r x i j = i = 1 n j = 1 m C T i j T i j m a x .
From this equation, we obtain:
i = 1 n j = 1 m lim b + 0 b 2 π r x i j ρ v 0 i j 2 1 Ψ ( x i j ) exp r i j 2 2 σ i j 2 d r i j i = 1 n j = 1 m lim b + 0 b 2 π r x i j ρ v 0 i j 2 1 Ψ ( x i j ) exp r i j 2 2 σ i j 2 2 d r i j = i = 1 n j = 1 m 1 2 ρ π r r i j 2 v 0 i j 2 C T i j , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which is also equivalent to
i = 1 n j = 1 m lim b + 0 b 2 r x i j 1 Ψ ( x i j ) exp r i j 2 2 σ i j 2 d r i j i = 1 n j = 1 m lim b + 0 b 2 r x i j 1 Ψ ( x i j ) exp r i j 2 2 σ i j 2 2 d r i j = i = 1 n j = 1 m 1 2 r r i j 2 C T i j , n = 1 , 2 , , N ; m = 1 , 2 , , N .
Computing this equation, then we have
i = 1 n j = 1 m lim b + r x i j 2 | 0 b i = 1 n j = 1 m lim b + 0 b Ψ ( x i j ) exp r i j 2 2 σ i j 2 r i j d r i j i = 1 n j = 1 m lim b + r x i j 2 | 0 b + 2 i = 1 n j = 1 m lim b + 0 b Ψ ( x i j ) exp r i j 2 2 σ i j 2 r i j d r i j i = 1 n j = 1 m lim b + 0 b Ψ 2 ( x i j ) exp 2 r i j 2 2 σ i j 2 r i j d r i j = i = 1 n j = 1 m 1 4 r r i j 2 C T i j , n = 1 , 2 , , N ; m = 1 , 2 , , N
or also
i = 1 n j = 1 m lim b + 0 b Ψ ( x i j ) exp r i j 2 2 σ i 2 r i j d r i j i = 1 n j = 1 m lim b + 0 b Ψ 2 ( x i j ) exp r i j 2 σ i 2 r i j d r i j = i = 1 n j = 1 m 1 4 r r i j 2 C T i j , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which can be also written as
i = 1 n j = 1 m lim b + Ψ ( x i j ) σ i j 2 exp r x i j 2 2 σ i j 2 | 0 b + i = 1 n j = 1 m lim b + Ψ 2 ( x i j ) σ i j 2 2 exp r i j 2 σ i j 2 | 0 b = i = 1 n j = 1 m 1 4 r r i j 2 C T i j , n = 1 , 2 , , N ; m = 1 , 2 , , N .
Then, we obtain
i = 1 n j = 1 m Ψ ( x i j ) σ i j 2 Ψ 2 ( x i j ) σ i j 2 2 = i = 1 n j = 1 m 1 4 C T i j r r i j 2 , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which is
i = 1 n j = 1 m Ψ 2 ( x i j ) 2 Ψ ( x i j ) + 1 2 C T i j r r i j σ i j 2 = 0 , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
or equivalently
i = 1 n j = 1 m Ψ 2 ( x i j ) 2 Ψ ( x i j ) + C T i j 2 σ i j r r i j 2 = 0 , n = 1 , 2 , , N ; m = 1 , 2 , , N .
By solving Equation (21), we obtain
i = 1 n j = 1 m Ψ ( x i j ) = i = 1 n j = 1 m 2 ± 4 2 C T i j ( σ i j / r r i j ) 2 2 = i = 1 n j = 1 m 1 ± 1 C T i j 2 ( σ i j / r r i j ) 2 = i = 1 n j = 1 m 1 ± 1 C T i j 2 r r i j σ i j 2 , n = 1 , 2 , , N ; m = 1 , 2 , , N .
From Equation (22), the physical solution is
i = 1 n j = 1 m Ψ ( x i j ) = i = 1 n j = 1 m 1 1 C T i j 2 r r i j σ i j 2 , n = 1 , 2 , , N ; m = 1 , 2 , , N .
By substituting this solution into Equation (9), we obtain
i = 1 n j = 1 m v x i j , r x i j = i = 1 n j = 1 m 1 1 1 C T i j 2 r r i j σ i j 2 exp r i j 2 2 σ i j 2 v 0 i j , n = 1 , 2 , , N ; m = 1 , 2 , , N .
Here, σ i j is recognized as a linear function of x i j in FGNDWM. In FGWM and FGNDWM, owing to every plane is perpendicular to the axis, the rates of mass flow are equal to each other in both FGWM and FGNDWM. According to the law of mass conservation, we can calculate and get:
i = 1 n j = 1 m 0 + 2 π r i j ρ i j v x i j , r x i j d r i j = i = 1 n j = 1 m 0 r x i j 2 π r i j ρ i j v x i j d r i j + i = 1 n j = 1 m r x i j + 2 π r i j ρ i v 0 i j d r i j , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which is
i = 1 n j = 1 m lim b + 0 b 2 π r i j ρ i j v x i j , r x i j d r i j = i = 1 n j = 1 m 0 r x i j 2 π r i j ρ i j v x i j d r i j + i = 1 n j = 1 m lim b + r x i j b 2 π r i j ρ i j v 0 i j d r i j , n = 1 , 2 , , N ; m = 1 , 2 , , N .
Substituting Equations (7) and (24) into Equation (26), we obtain
i = 1 n j = 1 m lim b + 0 b 2 π r i j ρ i j 1 1 1 C T i j 2 r r i j σ i j 2 exp r i j 2 2 σ i j 2 v 0 i j d r i j = i = 1 n j = 1 m 0 r x i j 2 π r i j ρ i j 1 1 2 1 1 2 C T i j ( r x i j / r r i j ) 2 v 0 i j d r i j + i = 1 n j = 1 m lim b + r x i j b 2 π r i j ρ i j v 0 i j d r i j , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
or equivalently
i = 1 n j = 1 m lim b + 0 b r i j 1 1 1 C T i j 2 r r i j σ i j 2 exp r i j 2 2 σ i j 2 d r i j = i = 1 n j = 1 m 0 r x i j r i j 1 1 2 1 1 2 C T i j ( r x i j / r r i j ) 2 d r i j + i = 1 n j = 1 m lim b + r x i j b r i j d r i j , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which can be also written as
i = 1 n j = 1 m lim b + 0 b r i j d r i j i = 1 n j = 1 m 1 1 C T i j 2 r r i j σ i 2 lim b + 0 b r i j exp r i j 2 2 σ i j 2 d r i j = i = 1 n j = 1 m 0 r x i j r i j 1 2 + 1 2 1 2 C T i j ( r x i j / r r i j ) 2 d r i j + i = 1 n j = 1 m lim b + r x i j b r i j d r i j n = 1 , 2 , , N ; m = 1 , 2 , , N ,
or equivalently
i = 1 n j = 1 m lim b + r i j 2 2 | 0 b i = 1 n 1 1 C T i j 2 r r i j σ i 2 lim b + 0 b r i j exp r i j 2 2 σ i j 2 d r i j = i = 1 n j = 1 m 1 2 r i 2 2 | 0 r x i j + i = 1 n j = 1 m 1 2 1 2 C T i j ( r x i j / r r i j ) 2 r i j 2 2 | 0 r x i j + i = 1 n j = 1 m lim b + r i j 2 2 | r r i j b n = 1 , 2 , , N ; m = 1 , 2 , , N .
From Equation (30), we have
i = 1 n j = 1 m 1 1 C T i j 2 r r i j σ i j 2 lim b + 0 b r i j exp r i j 2 2 σ i j 2 d r i j = i = 1 n j = 1 m 1 2 r x i j 2 2 + i = 1 n j = 1 m 1 2 1 2 C T i j ( r x i j / r r i j ) 2 r x i j 2 2 i = 1 n j = 1 m r x i j 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which can be expressed by the following simple equation:
i = 1 n j = 1 m 1 1 C T i j 2 r r i j σ i j 2 σ i 2 = i = 1 n j = 1 m r x i j 2 4 + i = 1 n j = 1 m r x i j 2 4 1 2 C T i j ( r x i j / r r i j ) 2 n = 1 , 2 , , N ; m = 1 , 2 , , N
or also
i = 1 n j = 1 m 1 1 C T i j 2 r r i j σ i j 2 σ i j 2 = i = 1 n j = 1 m r x i j 2 4 i = 1 n j = 1 m r x i j 2 4 1 2 C T i j ( r x i j / r r i j ) 2 n = 1 , 2 , , N ; m = 1 , 2 , , N .
Now, we can deduce the following results from Equation (33):
i = 1 n j = 1 m σ i j 2 σ i j σ i j 2 C T i j r r i j 2 2 = i = 1 n j = 1 m r x i j 2 4 i = 1 n j = 1 m r x i j 2 4 1 2 C T i j ( r x i j / r r i j ) 2 n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which leads to
i = 1 n j = 1 m σ i j 2 r x i j 2 4 1 1 2 C T i j ( r x i j / r r i j ) 2 = i = 1 n j = 1 m σ i j σ i j 2 C T i j r r i j 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N .
Now, by squaring the two sides of Equation (35), we obtain:
i = 1 n j = 1 m σ i j 4 r x i j 2 2 σ i j 2 1 1 2 C T i j ( r x i j / r r i j ) 2 + r x i j 4 16 1 1 2 C T i j ( r x i j / r r i j ) 2 2 = i = 1 n j = 1 m σ i j 2 σ i j 2 C T i j r r i j 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N .
By arranging the different terms, we obtain:
i = 1 n j = 1 m r x i j 2 2 1 1 2 C T i j ( r x i j / r r i j ) 2 σ i j 2 + C T i j r r i j 2 2 σ i j 2 = i = 1 n j = 1 m r x i j 4 16 1 1 2 C T i j ( r x i j / r r i j ) 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which gives the following equality:
i = 1 n j = 1 m σ i j 2 = i = 1 n j = 1 m r x i j 4 16 1 1 2 C T i j ( r x i j / r r i j ) 2 2 r x i j 2 2 1 1 2 C T i j ( r x i j / r r i j ) 2 + C T i j r r i j 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which can be rewritten as:
i = 1 n j = 1 m σ i j 2 = i = 1 n j = 1 m r x i j 2 4 × r x i j 2 4 1 2 1 2 C T i j ( r x i j / r r i j ) 2 + 1 2 C T i j ( r x i j / r r i j ) 2 r x i j 2 2 + r x i j 2 2 1 2 C T i j ( r x i j / r r i j ) 2 + C T i j r r i j 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which is equivalent to:
i = 1 n j = 1 m σ i j 2 = i = 1 n j = 1 m r x i j 2 4 r x i j 2 2 + r x i j 2 2 1 2 C T i j ( r x i j / r r i j ) 2 + C T i j r r i j 2 2 r x i j 2 2 + r x i j 2 2 1 2 C T i j ( r x i j / r r i j ) 2 + C T i j r r i j 2 2 n = 1 , 2 , , N ; m = 1 , 2 , , N .
Then, we deduce
i = 1 n j = 1 m σ i j 2 = i = 1 n j = 1 m r x i j 2 4 , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
which gives, for the real σ i , the following value
i = 1 n j = 1 m σ i j = i = 1 n j = 1 m r x i j 2 = i = 1 n j = 1 m r 0 i j + α i j x i j 2 = i = 1 n j = 1 m r 0 i j 2 + α i j 2 x i j , n = 1 , 2 , , N ; m = 1 , 2 , , N ,
where r x i j = r 0 i j + α i j x i j , α i j and r 0 i j can be given and estimated empirically from
i = 1 n j = 1 m α i j = i = 1 n j = 1 m 0.5 ln ( z h i j / z 0 i j ) , n = 1 , 2 , , N ; m = 1 , 2 , , N
and
i = 1 n j = 1 m r 0 i j = i = 1 n j = 1 m 4 5 r 2 i j , n = 1 , 2 , , N ; m = 1 , 2 , , N
in [44], respectively. Finally, Equations (9), (23), (24) and (42) constitute the FGNDWM.

4. Comparisons and Analysis of Two Different Wake Models for OSWFs

In this section, we will give the definition of Wind Speed Deficit (WSD) and further discuss the relationship between FGWM and FGNDWM. Usually, comparing WSD is a very important approach in different Wake Models. The WSD of OSWFs is expressed as by the following equation:
Δ v i j v 0 i j = v 0 i j v x i j v 0 i j , i = 1 , 2 , , N ; j = 1 , 2 , , N .
Firstly, the WSD of FGWM is derived by Equations (7) and (45):
Δ v F G W M i j = v 0 i j v x i j v 0 i j = v 0 i j 1 1 2 1 1 2 C T i j ( r x i j / r r i j ) 2 v 0 i j v 0 i j = 1 2 1 1 2 C T i j ( r x i j / r r i j ) 2 , i = 1 , 2 , , N ; j = 1 , 2 , , N .
Meanwhile, we obtain the WSD of FGNDWM. With the condition: r i j = 0 , the WSD of FGNDWM on the axis is interpreted based on Equations (24) and (45):
Δ v F G N D W M i j = v 0 i j v x i j , r i j v 0 i j = v 0 i j 1 1 1 C T i j 2 r r i j σ i j 2 exp r i j 2 2 σ i j 2 v 0 i j v 0 i j = 1 1 C T i j 2 r r i j σ i j 2 = 1 1 C T i j 2 σ i j r r i j 2 , i = 1 , 2 , , N ; j = 1 , 2 , , N .
Substituting σ i j = r x i j 2 into Equation (47), we obtain
Δ v F G N D W M i j = v 0 i j v x i j , r i j v 0 i j = 1 1 C T i j 2 r r i j r x i j 2 2 = 1 1 2 C T i j ( r x i j / r r i j ) 2 , i = 1 , 2 , , N ; j = 1 , 2 , , N .
The ratio of the WSD from FGWM to FGNDWM is calculated by Equations (46) and (48):
Δ v F G W M i j Δ v F G N D W M i j = 1 2 1 1 2 C T i j ( r x i j / r r i j ) 2 1 1 2 C T i j ( r x i j / r r i j ) 2 = 1 2 , i = 1 , 2 , , N ; j = 1 , 2 , , N .
We can find out the WSD of FGWM is half as small as that of FGNDWM on the axis.

5. Experimental Comparisons and Analysis of Two Different Wake Models for OSWFs

In this section, we collect and use data of five cases to confirm the different characteristics between FGWM and FGNDWM. Usually, using miniature WT with the Large-Eddy Simulation (LES) data was known as the standard case in the literature.
The main data and parameters of the five cases (OSWFs in Yangshan port, Shanghai) are shown in the following Table 1, in which z h i j is the height of the Hub, and z 0 i j is the rate of surface sea or roughness. These roughness lengths shown for C a s e s ( b e ) in Table 1 are representative of different sea surface types, including very rough terrain, for instance, islands with different sizes ( z 0 i j = 0.1 m), sea surface with reefs, rocks and shoal rocks ( z 0 i j = 0.01 m), sea surface with medium waves and large waves ( z 0 i j = 0.001 m), and sea surface with small waves ( z 0 i j = 0.00001 m). d r i j is the diameter of the rotor, α i j is the axial induction factor, r 0 i j is the downstream rotor radius, C T i j is the thrust coefficient, and v 0 i j is cut-in wind speed. Here, i and j are row vector and column vector for Wind Turbine in large-scale offshore Wind Farms, respectively.
LES are applied in many fields of flow simulations. The initial conditions have a very significant influence on the LES results. x i j / d r i j = 3 and Δ v i j / v 0 i j m a x = 0.5 were chosen as the initial conditions used for LES simulation conducted in this study.
Using some given data in Table 1, we conduct the simulation experiments with FGWM and FGNDWM. By analyzing Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9, from simulations and experiments, we can find the WSD of FGWM are half times smaller then that of FGNDWM on the axis. These simulations validate the results obtained in the last sections.
If the length of the near wake region x 1 i is taken into consideration and the parameters are selected according to Table 1, the result is shown in Figure 5. The FGNDWM appeared superior to FGWM. The problem was that FGWM did not take x 1 i into consideration, whereas the FGNDWM took it into consideration. In addition, the characteristic width of the FGNDWM was obtained by fitting the LES data in the experiment of this study.
It can be seen from Figure 8 that the accuracy of the FGNDWM and the FGWM is better than that in Figure 7 and Figure 9. The maximum Wind Speed Deficit (WSD) of this FGNDWM can be proved to be twice as large as that of the FGWM; if α i j and r 0 i j are estimated by (43) and (44), respectively, the result is shown in Figure 9. As previously stated, (43) is not applicable to case a and case e. As a result, the accuracy of the FGNDWM and FGWM in Figure 6 is worse than that in Figure 5 and Figure 7, Figure 8 and Figure 9.

6. Analysis and Enlightenment of Wind Rose, Wind Weibull Probability Density Distribution and ELM Prediction

Based on the actual situation in OSWFs in Yangshan port, Shanghai, this section summarizes and describes the analysis of wind rose, wind Weibull probability density distribution and ELM prediction. We obtain the real-time and actual data from this website [45]. The WE cases are analyzed and studied based on mathematical models, and WRs are abstracted through the variable wind directions and wind speeds in OSWFs in Yangshan port, Shanghai. These data of WE resources are collected and shown in the following figures.
From 22 August 2017 to 22 August 2018, the portraits of wind direction and wind speed (m/s) in OSWFs of Yangshan port, Shanghai are shown with the wind rose in Figure 10. From the southwest direction, most of the wind mean speeds in a whole year are greater than 6 m/s, whereas, from the northwest direction, a small part of wind mean speeds are close to 6 m/s in a whole year.
From 22 August 2016 to 22 August 2017 with a whole year in OSWFs of Yangshan port, Shanghai, we keep obtaining the real-time and actual data from this website [45]. The wind speeds are collected and their Mean Wind Speed (MWS) is 3.4934 m/s. The details of them are shown in Figure 11.
From 22 August 2016 to 22 August 2017 with a whole year in OSWFs of Yangshan port, Shanghai, the rose portraits of Wind Direction (WD) and Average Wind Direction (AWD) are shown in Figure 12. Among them, the north direction is 0 degrees (North = 0 ℃). The probabilities of most of the parts of WDs from the southwest direction are greater than 5% and are close to 6%. Meanwhile, a small part of WD from the northeast direction is greater than 5% and close to 6%.
In Figure 13, the wind rose shows the portraits of wind mean speed (m/s) in OSWFs of Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 with a whole year. The majority of the wind mean speeds from the southwest are greater than 5 m/s, whereas, a small amount of wind mean speed from the northwest is approximate to 5 m/s. Therefore, WTs should adjust the direction to the southwest in OSWFs of Yangshan port, Shanghai throughout the summer, for even more periods.
Figure 14, Figure 15 and Figure 16 show all the wind situations, including the time series of air relative humidity, wind direction and wind temperature in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018, respectively. In this period, the mean air relative humidity is 82.0978 % , the mean wind direction is 169.4492 (North direction = 0 ), and the mean wind temperature is 10.0367 ℃.
ELM is a learning algorithm, initially introduced to train a Single Layer Feedforward Neural network [46]. In ELM theory, the input weights are randomly generated according to any continuous distribution function, while the output weights are analytically computed by the minimum norm solution of a linear system.
Here, as shown in Figure 17, the proposed ELM can be seen as three hidden layer neural networks, trained using the ELM algorithm. ELM is applied to wind direction and wind speed prediction in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018. The wind direction and wind speed of OSWFs are relevant to local air relative humidity and local wind temperature. Some simulation results show the corresponding rationality. Figure 18 and Figure 19 show the portraits of comparison of wind direction and wind speed forecasting results with ELM in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018, respectively. Here, the number of forecast data pieces is 2248. Figure 20 and Figure 21 show the portraits of comparison of wind direction and wind speed forecasting results with ELM in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018, respectively. At the same time, the number of forecast data pieces is 500.
Through the study of the Weibull Probability Distribution of Wind Velocity Data in OSWFs in Yangshan port, Shanghai, we obtain the following results. The cumulative distribution and linearized curve are plotted and shown in Figure 22, Linearized curve and fitted line comparison are shown in Figure 23, Weibull probability density function and Cumulative Weibull probability density function are shown in Figure 24. In light of wind tower, measuring data and plotted wind speed histograms, we inferred and estimated two parameters of the Weibull distribution by maximum likelihood estimation method, i.e., c = 3.7660 , and k = 1.7153 . The histogram of wind speed at hub height with the fitted Weibull probability density distribution are plotted and shown in Figure 25. Wind speed at hub height conforms to the Weibull distribution.

7. Conclusions

From this study, we found that the feature of FGWM is simple and intuitive, while the feature of FGNDWM is complex and precise. The accuracy of the FGNDWM is inherently better than that of the FGWM. When describing the relationship between FGNDWM and FGWM, Equation (42) can reflect the essential characteristics of FGNDWM, whereas Equation (7) reflects the basic characteristics of FGWM. Equation (42) is more accurate than (7) in expressing the characteristics of the large off-shore wind farms. By taking x 1 i into consideration, the accuracy of the FGNDWM and the FGWM can be improved. Their accuracy depends on the axial induction factor α i j . The maximum Wind Speed Deficit (WSD) of this FGNDWM can be proved to be twice as large as that of the FGWM if α i j and r 0 i j are estimated by (43) and (44), respectively.
Currently, the experiments show that the accuracy of ELM predictions needs to be improved based on the actual situation in OSWFs in Yangshan port, Shanghai. In future research, we will work on hybrid wake models for near wakes and far wakes and improve the ELM predictions accuracy of large off-shore wind farms. The data assimilation and reduced order modelling will be provided in a future paper concerning induced large off-shore wind farms dynamics.

Author Contributions

Conceptualization, M.L. and H.X.; Methodology, L.P. and C.X.; Software, L.P. and C.X.; Validation, M.L., H.X. and L.P.; Formal Analysis, M.L. and H.X.; Investigation, M.L., L.P. and C.X.; Resources, M.L. and H.X.; Data Curation, M.L. and H.X.; Writing–Original Draft Preparation, L.P. and C.X.; Writing–Review and Editing, M.L. and H.X.; Visualization, L.P. and C.X.; Supervision, H.X. and L.P.; Project Administration, H.X.; Funding Acquisition, H.X.

Funding

This research was funded by the PhD research startup foundation of Wuhan University of Technology grant number 40120291, the Fundamental Research Funds for the Central Universities grant number 173118001, the Open Project Program of National Engineering Research Center for Water Transport Safety grant number A2019006, and the Open Project Program of State Key Laboratory of Ocean Engineering grant number 1812. The APC was funded by National Key Technology Research and Development Program of the Ministry of Science and Technology of China grant number 2015BAF06B01-2.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
OSWFsOff-Shore Wind Farms
WTWind Turbine
ELMExtreme Learning Machine
WPWind Power
LESLarge-Eddy Simulation
EEAExtended Exergy Accounting
WPWind Power
HAWTsHorizontal Axis Wind Turbines
WFWind Farm
WEWind Energy
DFIGDoubly-fed Induction Generator
PMSGPermanent Magnet Synchronous Generator
FSIGFixed-Speed Induction Generator
WTPGSWind Turbine Power Generation System
FGWMFrandsen Generalized Wake Model
FGNDWMFrandsen Generalized Normal Distribution Wake Model
WSDWind Speed Deficit
WDWind Direction
AWDAverage Wind Direction

References

  1. Marden, J.R.; Ruben, S.D.; Pao, L.Y. A Model-Free Approach to Wind Farm Control Using Game Theoretic Methods. IEEE Trans. Control Syst. Technol. 2013, 21, 1207–1214. [Google Scholar] [CrossRef]
  2. Pao, L.Y.; Johnson, K. A tutorial on the dynamics and control of wind turbines and wind farms. In Proceedings of the American Control Conference, St. Louis, MO, USA, 10–12 June 2009; pp. 2076–2089. [Google Scholar]
  3. Bitar, E.; Seiler, P. Coordinated control of a wind turbine array for power maximization. In Proceedings of the American Control Conference (ACC), Washington, DC, USA, 17–19 June 2013; pp. 2898–2904. [Google Scholar]
  4. Remmers, T.; Cawkwell, F.; Desmond, C.; Murphy, J.; Politi, E. The Potential of Advanced Scatterometer (ASCAT) 12.5 km Coastal Observations for Offshore Wind Farm Site Selection in Irish Waters. Energies 2019, 12, 206. [Google Scholar] [CrossRef]
  5. Hübler, C.; Weijtjens, W.; Gebhardt, C.G.; Rolfes, R.; Devriendt, C. Validation of Improved Sampling Concepts for Offshore Wind Turbine Fatigue Design. Energies 2019, 12, 603. [Google Scholar] [CrossRef]
  6. Li, H.; Wang, J.; Lu, H.; Guo, Z. Research and application of a combined model based on variable weight for short term wind speed forecasting. Renew. Energy 2018, 116, 669–684. [Google Scholar] [CrossRef]
  7. Romanic, D.; Parvu, D.; Refan, M.; Hangan, H. Wind and tornado climatologies and wind resource modelling for a modern development situated in Tornado Alley. Renew. Energy 2018, 115, 97–112. [Google Scholar] [CrossRef]
  8. Ahmed, A.S. Wind energy characteristics and wind park installation in Shark El-Ouinat, Egypt. Renew. Sustain. Energy Rev. 2018, 82, 734–742. [Google Scholar] [CrossRef]
  9. Aghbashlo, M.; Tabatabaei, M.; Hosseini, S.S.; Dashti, B.B.; Soufiyan, M.M. Performance assessment of a wind power plant using standard exergy and extended exergy accounting (EEA) approaches. J. Clean. Prod. 2018, 171, 127–136. [Google Scholar] [CrossRef]
  10. Ahmad, T.; Basit, A.; Anwar, J.; Coupiac, O.; Kazemtabrizi, B.; Matthews, P.C. Fast Processing Intelligent Wind Farm Controller for Production Maximisation. Energies 2019, 12, 544. [Google Scholar] [CrossRef]
  11. Shakoor, R.; Hassan, M.Y.; Raheem, A.; Wu, Y.K. Wake effect modeling: A review of wind farm layout optimization using Jensen’s model. Renew. Sustain. Energy Rev. 2016, 58, 1048–1059. [Google Scholar] [CrossRef]
  12. Kuenzel, S.; Kunjumuhammed, L.; Pal, B.; Erlich, I. Impact of Wakes on Wind Farm Inertial Response. IEEE Trans. Sustain. Energy 2014, 5, 237–245. [Google Scholar] [CrossRef]
  13. Pan, L.; Voos, H.; Li, Y.; Darouach, M.; Xu, Y.; Hu, S. A wake interaction model for the coordinated control of Wind Farms. In Proceedings of the 2015 IEEE 20th Conference on Emerging Technologies Factory Automation (ETFA), Luxembourg, 8–11 September 2015; pp. 1–7. [Google Scholar]
  14. Pan, L.; Voos, H.; Pan, Y.; Darouach, M. A generalized interaction Wake Model with its variation for control in Wind Farms. In Proceedings of the 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; pp. 8759–8764. [Google Scholar]
  15. Pan, L.; Voos, H.; Li, Y.; Xu, Y.; Darouach, M.; Li, Z. A class of Improved Wake Interaction Model for the coordinated control of wind farms. In Proceedings of the Chinese Automation Congress (CAC), Wuhan, China, 27–29 November 2015; pp. 1322–1327. [Google Scholar]
  16. He, P.; Arefifar, S.A.; Li, C.; Wen, F.; Ji, Y.; Tao, Y. Enhancing Oscillation Damping in an Interconnected Power System with Integrated Wind Farms Using Unified Power Flow Controller. Energies 2019, 12, 322. [Google Scholar] [CrossRef]
  17. March, V. Key issues to define a method of lightning risk assessment for wind farms. Electr. Power Syst. Res. 2017, in press. [Google Scholar] [CrossRef]
  18. Chen, K.; Song, M.; Zhang, X. The investigation of tower height matching optimization for wind turbine positioning in the wind farm. J. Wind Eng. Ind. Aerodyn. 2013, 114, 83–95. [Google Scholar] [CrossRef]
  19. Sørensen, K.L.; Galeazzi, R.; Odgaard, P.F.; Niemann, H.; Poulsen, N.K. Adaptive Passivity Based Individual Pitch Control for Wind Turbines in the Full Load Region. In Proceedings of the 2014 American Control Conference, Portland, OR, USA, 4–6 June 2014; pp. 554–559. [Google Scholar]
  20. Barreiro-Gomez, J.; Ocampo-Martinez, C.; Bianchi, F.; Quijano, N. Model-free control for wind farms using a gradient estimation-based algorithm. In Proceedings of the 2015 European Control Conference (ECC), Linz, Austria, 15–17 July 2015; pp. 1516–1521. [Google Scholar]
  21. Long, M.; Becerra, M.; Thottappillil, R. On the attachment of dart lightning leaders to wind turbines. Electr. Power Syst. Res. 2017, 151, 432–439. [Google Scholar] [CrossRef]
  22. Thukaram, D. Accurate modeling of doubly fed induction generator based wind farms in load flow analysis. Electr. Power Syst. Res. 2018, 155, 363–371. [Google Scholar]
  23. Farajzadeh, S.; Ramezani, M.H.; Nielsen, P.; Nadimi, E.S. Statistical modeling of the power grid from a wind farm standpoint. Electr. Power Syst. Res. 2017, 144, 150–156. [Google Scholar] [CrossRef]
  24. Tian, L.; Zhu, W.; Shen, W.; Zhao, N.; Shen, Z. Development and validation of a new two-dimensional wake model for wind turbine wakes. J. Wind Eng. Ind. Aerodyn. 2015, 137, 90–99. [Google Scholar] [CrossRef] [Green Version]
  25. Park, J.; Law, K.H. A data-driven, cooperative wind farm control to maximize the total power production. Appl. Energy 2016, 165, 151–165. [Google Scholar] [CrossRef]
  26. Marseglia, G.R.; Arbasini, A.; Grassi, S.; Raubal, M.; Raimondo, D.M. Optimal placement of wind turbines on a continuous domain: An MILP-based approach. In Proceedings of the American Control Conference (ACC), Chicago, IL, USA, 1–3 July 2015; pp. 5010–5015. [Google Scholar]
  27. Van Dam, F.; Gebraad, P.; van Wingerden, J.W. A maximum power point tracking approach for wind farm control. In Proceedings of the Science of Making Torque from Wind, Oldenburg, Germany, 9–11 October 2012. [Google Scholar]
  28. Chen, J. Development of offshore wind power in China. Renew. Sustain. Energy Rev. 2011, 15, 5013–5020. [Google Scholar] [CrossRef]
  29. Zhixin, W.; Chuanwen, J.; Qian, A.; Chengmin, W. The key technology of offshore wind farm and its new development in China. Renew. Sustain. Energy Rev. 2009, 13, 216–222. [Google Scholar] [CrossRef]
  30. Ebrahimi, F.; Khayatiyan, A.; Farjah, E. A novel optimizing power control strategy for centralized wind farm control system. Renew. Energy 2016, 86, 399–408. [Google Scholar] [CrossRef]
  31. Song, Z.; Zhang, Z.; Chen, X. The decision model of 3-dimensional wind farm layout design. Renew. Energy 2016, 85, 248–258. [Google Scholar] [CrossRef]
  32. Varzaneh, S.G.; Abedi, M.; Gharehpetian, G. A new simplified model for assessment of power variation of DFIG-based wind farm participating in frequency control system. Electr. Power Syst. Res. 2017, 148, 220–229. [Google Scholar] [CrossRef]
  33. Hossain, M.E. A non-linear controller based new bridge type fault current limiter for transient stability enhancement of DFIG based Wind Farm. Electr. Power Syst. Res. 2017, 152, 466–484. [Google Scholar] [CrossRef]
  34. Yao, J.; Li, J.; Guo, L.; Liu, R.; Xu, D. Coordinated control of a hybrid wind farm with PMSG and FSIG during asymmetrical grid fault. Int. J. Electr. Power Energy Syst. 2018, 95, 287–300. [Google Scholar] [CrossRef]
  35. Li, D.Y.; Li, P.; Cai, W.C.; Song, Y.D.; Chen, H.J. Adaptive Fault Tolerant Control of Wind Turbines with Guaranteed Transient Performance Considering Active Power Control of Wind Farms. IEEE Trans. Ind. Electron. 2017, 65, 3275–3285. [Google Scholar] [CrossRef]
  36. Chaurasiya, P.K.; Ahmed, S.; Warudkar, V. Comparative analysis of Weibull parameters for wind data measured from met-mast and remote sensing techniques. Renew. Energy 2018, 115, 1153–1165. [Google Scholar] [CrossRef]
  37. Atighechi, H.; Hu, P.; Ebrahimi, S.; Lu, J.; Wang, G.; Wang, L. An effective load shedding remedial action scheme considering wind farms generation. Int. J. Electr. Power Energy Syst. 2018, 95, 353–363. [Google Scholar] [CrossRef]
  38. Suganthi, S.; Devaraj, D.; Ramar, K.; Thilagar, S.H. An Improved Differential Evolution algorithm for congestion management in the presence of wind turbine generators. Renew. Sustain. Energy Rev. 2018, 81, 635–642. [Google Scholar] [CrossRef]
  39. Wan, C.; Xu, Z.; Pinson, P.; Dong, Z.Y.; Wong, K.P. Probabilistic forecasting of wind power generation using extreme learning machine. IEEE Trans. Power Syst. 2014, 29, 1033–1044. [Google Scholar] [CrossRef]
  40. Nikolić, V.; Motamedi, S.; Shamshirband, S.; Petković, D.; Ch, S.; Arif, M. Extreme learning machine approach for sensorless wind speed estimation. Mechatronics 2016, 34, 78–83. [Google Scholar] [CrossRef]
  41. Lazarevska, E. Wind speed prediction with extreme learning machine. In Proceedings of the 2016 IEEE 8th International Conference on Intelligent Systems (IS), Sofia, Bulgaria, 4–6 September 2016; pp. 154–159. [Google Scholar]
  42. Wu, S.; Wang, Y.; Cheng, S. Extreme learning machine based wind speed estimation and sensorless control for wind turbine power generation system. Neurocomputing 2013, 102, 163–175. [Google Scholar] [CrossRef]
  43. Frandsen, S.; Barthelmie, R.; Pryor, S.; Rathmann, O.; Larsen, S.; Højstrup, J.; Thøgersen, M. Analytical modelling of wind speed deficit in large offshore wind farms. Wind Energy 2006, 9, 39–53. [Google Scholar] [CrossRef]
  44. Pookpunt, S.; Ongsakul, W. Optimal placement of wind turbines within wind farm using binary particle swarm optimization with time-varying acceleration coefficients. Renew. Energy 2013, 55, 266–276. [Google Scholar] [CrossRef]
  45. Iowa State University of Science and Technology The Iowa Environmental Mesonet (IEM). Available online: http://mesonet.agron.iastate.edu/ (accessed on 5 March 2019).
  46. Huang, G.B.; Zhu, Q.Y.; Siew, C.K. Extreme learning machine: A new learning scheme of feedforward neural networks. Neural Netw. 2006, 2, 985–990. [Google Scholar]
Figure 1. China’s first off-shore wind farm in Yangshan Port.
Figure 1. China’s first off-shore wind farm in Yangshan Port.
Energies 12 00863 g001
Figure 2. Off-shore wind farm of Yangshan Deepwater Harbor in the Port of Shanghai, China.
Figure 2. Off-shore wind farm of Yangshan Deepwater Harbor in the Port of Shanghai, China.
Energies 12 00863 g002
Figure 3. The portrait of Frandsen Generalized Wake Model (FGWM) with stream tube [43].
Figure 3. The portrait of Frandsen Generalized Wake Model (FGWM) with stream tube [43].
Energies 12 00863 g003
Figure 4. The portrait of Frandsen Generalized Normal Distribution Wake Model (FGNDWM) with stream.
Figure 4. The portrait of Frandsen Generalized Normal Distribution Wake Model (FGNDWM) with stream.
Energies 12 00863 g004
Figure 5. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e a (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Figure 5. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e a (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Energies 12 00863 g005
Figure 6. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e b (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Figure 6. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e b (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Energies 12 00863 g006
Figure 7. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e c (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Figure 7. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e c (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Energies 12 00863 g007
Figure 8. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e d (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Figure 8. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e d (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Energies 12 00863 g008
Figure 9. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e e (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Figure 9. The portrait of the maximum Wind Speed Deficit (WSD) in C a s e e (fitted by LES data, i = 1 , 2 , , N ; j = 1 , 2 , , N ).
Energies 12 00863 g009
Figure 10. Rose portrait of WD (North = 0) and speed (m/s) in off-shore wind farms (OSWFs) in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 10. Rose portrait of WD (North = 0) and speed (m/s) in off-shore wind farms (OSWFs) in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g010
Figure 11. The portraits of wind speed time series and mean speed ( v ¯ = 3.4934 m/s) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 11. The portraits of wind speed time series and mean speed ( v ¯ = 3.4934 m/s) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g011
Figure 12. Rose portrait of WD (North direction = 0 ) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 12. Rose portrait of WD (North direction = 0 ) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g012
Figure 13. Rose portrait of Wind Mean Speed (m/s) in OSWFs of Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 13. Rose portrait of Wind Mean Speed (m/s) in OSWFs of Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g013
Figure 14. The portraits of air relative humidity in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 14. The portraits of air relative humidity in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g014
Figure 15. The portraits of wind direction (north direction = 0 ) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 15. The portraits of wind direction (north direction = 0 ) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g015
Figure 16. The portraits of wind temperature in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 16. The portraits of wind temperature in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g016
Figure 17. The proposed Extreme Learning Machine (ELM) can be seen as a three hidden layer neural network trained using the ELM algorithm.
Figure 17. The proposed Extreme Learning Machine (ELM) can be seen as a three hidden layer neural network trained using the ELM algorithm.
Energies 12 00863 g017
Figure 18. The portraits of comparison of wind speed forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 2248).
Figure 18. The portraits of comparison of wind speed forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 2248).
Energies 12 00863 g018
Figure 19. The portraits of comparison of wind direction forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 2248).
Figure 19. The portraits of comparison of wind direction forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 2248).
Energies 12 00863 g019
Figure 20. The portraits of comparison of wind speed forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 500).
Figure 20. The portraits of comparison of wind speed forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 500).
Energies 12 00863 g020
Figure 21. The portraits of comparison of wind direction forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 500).
Figure 21. The portraits of comparison of wind direction forecasting results (ELM) in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018 (number of predicted data pieces is 500).
Energies 12 00863 g021
Figure 22. The portraits of distribution and cumulative distribution extracted from the time series in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 22. The portraits of distribution and cumulative distribution extracted from the time series in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g022
Figure 23. The portraits of linearized curve and fitted line comparison in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 23. The portraits of linearized curve and fitted line comparison in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g023
Figure 24. The portraits of Weibull and Cumulative Weibull probability density functions in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 24. The portraits of Weibull and Cumulative Weibull probability density functions in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g024
Figure 25. The portraits of wind speed histogram in hub height and the fitted Weibull probability density distribution in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Figure 25. The portraits of wind speed histogram in hub height and the fitted Weibull probability density distribution in OSWFs in Yangshan port, Shanghai from 22 August 2017 to 22 August 2018.
Energies 12 00863 g025
Table 1. Different experiments and Large-Eddy Simulation (LES) case to validate the Generalized model of OSWFs in Yangshan port, Shanghai.
Table 1. Different experiments and Large-Eddy Simulation (LES) case to validate the Generalized model of OSWFs in Yangshan port, Shanghai.
Cases d rij (m) z hij (m) v 0 ij (m/s) C Tij z 0 ij (m) α ij r 0 ij (m)
C a s e a 0.150.1252.20.41940.000030.1190.066
C a s e b 666560.39160.10.1138.08
C a s e c 669860.29440.010.0840.48
C a s e d 827560.22560.0010.0643.52
C a s e e 706560.22560.000010.06241.12
Note: C a s e a represents the Large-Eddy Simulation (LES) data; C a s e b represents OSWF-b in Yangshan port, Shanghai; C a s e c represents OSWF-c in Yangshan port, Shanghai; C a s e d represents OSWF-d in Yangshan port, Shanghai; C a s e e represents OSWF-e in Yangshan port, Shanghai.

Share and Cite

MDPI and ACS Style

Li, M.; Xiao, H.; Pan, L.; Xu, C. Study of Generalized Interaction Wake Models Systems with ELM Variation for Off-Shore Wind Farms. Energies 2019, 12, 863. https://doi.org/10.3390/en12050863

AMA Style

Li M, Xiao H, Pan L, Xu C. Study of Generalized Interaction Wake Models Systems with ELM Variation for Off-Shore Wind Farms. Energies. 2019; 12(5):863. https://doi.org/10.3390/en12050863

Chicago/Turabian Style

Li, Mingcan, Hanbin Xiao, Lin Pan, and Chengjun Xu. 2019. "Study of Generalized Interaction Wake Models Systems with ELM Variation for Off-Shore Wind Farms" Energies 12, no. 5: 863. https://doi.org/10.3390/en12050863

APA Style

Li, M., Xiao, H., Pan, L., & Xu, C. (2019). Study of Generalized Interaction Wake Models Systems with ELM Variation for Off-Shore Wind Farms. Energies, 12(5), 863. https://doi.org/10.3390/en12050863

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop