Next Article in Journal
On the Need for a Paradigm Change in the Valuation of Concrete with Waste Materials Based on the Example of Concrete with Crumb Rubber
Next Article in Special Issue
Satisfaction with Media Information about Renewable Energy Investments
Previous Article in Journal
Receiving Robust Analysis of Spatial and Temporary Variation of Agricultural Water Use Efficiency While Considering Environmental Factors: On the Evaluation of Data Envelopment Analysis Technique
Previous Article in Special Issue
Deterioration of Novel Silver Coated Mirrors on Polycarbonate Used for Concentrated Solar Power
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Prediction of Wind Energy Resources Based on WOA—A Case Study in Jordan

by
Ayman Al-Quraan
1,*,
Bashar Al-Mhairat
1,
Ahmad M. A. Malkawi
2,
Ashraf Radaideh
1 and
Hussein M. K. Al-Masri
1
1
Electrical Power Engineering Department, Hijjawi Faculty for Engineering Technology, Yarmouk University, Irbid 21163, Jordan
2
Mechatronics Engineering Department, The University of Jordan, Amman 11942, Jordan
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(5), 3927; https://0-doi-org.brum.beds.ac.uk/10.3390/su15053927
Submission received: 3 January 2023 / Revised: 3 February 2023 / Accepted: 16 February 2023 / Published: 21 February 2023
(This article belongs to the Special Issue Renewable Energy Technologies for Sustainable Development)

Abstract

:
The average wind speed in a given area has a significant impact on the amount of energy that can be harvested by wind turbines. The regions with the most attractive possibilities are typically those that are close to the seaside and have open terrain inland. There is also good potential in several mountainous locations. Despite these geographical restrictions on where wind energy projects can be located, there is enough topography in most of the world’s regions to use wind energy projects to meet a significant amount of the local electricity needs. This paper presents a new method of energy prediction of wind resources in several wind sites in Jordan, which can be used to decide whether a specific wind site is suitable for wind farm installation purposes. Three distribution models, Weibull, Gamma and Rayleigh, were employed to characterize the provided wind data. Different estimation methods were used to assign the parameters associated with each distribution model and the optimal parameters were estimated using whale optimization algorithms which reduce the error between the estimated and the measured wind speed probability. The distribution models’ performance was investigated using three statistical indicators. These indicators were: root mean square error (RMSE), coefficient of determination (R2), and mean absolute error (MAE). Finally, using the superlative distribution models, the wind energy for the chosen wind sites was estimated. This estimation was based on the calculation of the wind power density (ED) and the total wind energy (ET) of the wind regime. The results show that the total wind energy ranged from slightly under 100 kWh/m2 to nearly 1250 kWh/m2. In addition, the sites recording the highest estimated wind energy had the optimum average wind speed and the most symmetrical distribution pattern.

1. Introduction

Currently, renewable energy is a crucial and essential need to mitigate the hazards caused by fossil fuels. Although the global capacity of renewable energy in 2019 was 2588 GW, the percentage of electricity generation from renewable resources was about 27.3% versus 72.7% for traditional resources [1]. Therefore, all essential steps must be taken by governments worldwide to support these strategies.
One of the most highly developed sectors in the field of renewable energy is the wind energy sector. In 2015, 63.8 GW of wind power capacity was installed in networks throughout the globe, recording the greatest yearly increase in the capacity of wind power generation. In 2019, the second largest yearly increase in the capacity of wind power generation was recorded with a value of 60 GW [1,2].
The climate of Jordan enhances the utilization of renewable energy sources. It contains various locations that are appealing for wind energy investment. However, Jordan is still 93% dependent on its energy from gas and oil imported from other countries. In 2018, renewable energy contributed 10.7% of the total national energy production, in which the wind energy contribution was around 3.52% [1].
A wind energy prediction process can be accomplished by determining the speed and direction of the wind [3]. However, these two parameters cannot be specified precisely, due to the wind’s erratic behavior. The characteristics of wind should be established based on analyzing and recognizing the natural behavior of wind at the candidate site. The parts of the wind energy conversion system (WECS) can be efficiently designed once the wind characteristics have been established.
Any wind energy project has a lifecycle that consists of several phases and several steps. Accomplishing each phase properly leads to the success of the project. The prediction of wind resources is the first stage in the lifecycle of a wind energy project. Generally, wind resource prediction involves four main steps which are individually discussed as follows [3,4].
The first step is “site identification and data collection”: this step includes identifying one or more sites across countries, states, or provinces. Wind atlases or airport wind data are used to provide a preliminary prediction of the wind potential in those sites. Several aspects should be considered when identifying sites such as: market situations, local society and government support, transmission access and capacity, site constructability, and environmental aspects. Wind data are measured and collected from the metrological stations or airports near the candidate wind sites. Wind data are usually measured in ten minutes basis at 10 m height, which is considered the standard height.
The second step is “analyzing the collected data”: in this step, the collected data (mean wind speed) are analyzed by utilizing a statistical distribution model to obtain the frequency distribution of the average wind velocity for the candidate site over a certain period. Each distribution model includes several parameters that are assigned by utilizing a variety of estimation methods. Such methods can be classified into two main categories: the numerical methods and the optimization algorithms.
The third step is “choosing the best performance indicator model”: during this stage, a comparison between the selected distribution models must be carried out to choose the optimal distribution function. This is achieved by using various types of goodness of fit indicators such as the root mean square (RMS) error, the coefficient of determination (R2), chi-square (χ2), and others.
The fourth step is “estimating the wind energy”: this step involves evaluating the extracted energy from both the wind regime and wind turbines that are installed in the candidate sites based on the best distribution models. In this stage, the distribution model can be represented as a probability distribution function (PDF). Then, the power curve model (P-V characteristic) is utilized to evaluate the extracted wind energy from a specific wind turbine installed in a candidate wind site. By the end of this phase, it will be feasible to determine whether or not the potential wind site can support the installation of a wind turbine.
Several research studies have made great efforts in the evaluation of wind energy potential, attempting to find the optimum wind distributions and estimation methods [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30]. Most of these studies focused on parametric distribution models and numerical estimation methods. However, some researchers considered the non–parametric distribution models [13]. Various parametric models were utilized to assess wind energy resources such as Weibull [4,5,6,7,9,10,11,12,13,14,15], Rayleigh [8,18,19,20,21], Lognormal [21,22,23,24,25,26], and others.
According to the literature, the Weibull distribution was the most widely utilized by researchers, especially when optimizing algorithms were used to determine their parameters. Based on wind data gathered for over a year from a measuring station built at the Atlm University campus area in Ankara, Turkey, Bilir et al. used the Weibull distribution as an assessment tool in [7]. Wang et al., in [11], showed that the Weibull model was the most suitable distribution among the five selected distribution models, which were tested according to the RMSE indicator. In [14], Li et al. conducted a comparison of onshore and offshore wind profiles and their potentials for generating wind energy in two sites in China’s southeast coastline area. The authors verified the reliability of utilizing the Weibull model in both candidate sites. The results showed that the wind energy at an offshore location was greater than onshore wind energy for a specified region. The Weibull distribution was compared with the Rayleigh, Gamma, and Lognormal distributions in [18] by Jiang et al. Afterward, considering the low wind profile, the Weibull model was utilized for evaluating the associated potentials in China.
An analysis of 46 studies released during the years of 2010–2018 was performed in [27]. This analysis demonstrated that in 44 studies, the Weibull distribution was the optimum method. The Weibull distribution has not always been thought to be the best option [22,24]. Wu et al., in [24], employed the Weibull, Logistic, and Lognormal distributions to assess wind energy at typical locations in Inner Mongolia, China. The Weibull performance was the weakest, and the Logistic and Lognormal performances were better. According to [28], a number of mixed distributions perform better than the traditional Weibull distribution, including the Gamma–Weibull function (GW), bimodal Weibull function (WW), mixture truncated normal function (NN), and truncated Normal–Weibull function (NW).
A variety of estimation techniques, including the maximum likelihood method (MLM) [4,7], moment method (MM) [10,12,13,14,15], least square estimation (LSE) [18,19], empirical method (EM) [10,11,12,15], power density method (PDM) [14,26], and energy pattern factor (EPF) [9,15], should be used to accurately value the parameters of the associated models. The modified maximum likelihood (MMLM) [25] and modified energy pattern method (MEPM) [12] are two examples of modified approaches that can also be used. The first four mentioned methods were most commonly used by many studies. In [9], authors found that the MLM was superior over the other methods according to three performance indicators, which were the root mean square error (RMSE), coefficient of determination (R2), and chi-square (χ2). The EM was the best choice for wind energy assessment in three candidate potential sites in Australia to allocate parameters of the Weibull distribution, as presented in [10]. Gugliani et al. [12] proposed a new estimation method that was compared with six common estimation methods via a Monte Carlo simulation study with sample sizes varying from 100 to 100,000. For each sample size, the MEPF achieved the best results and was superior over the other methods.
Some researchers proposed non-parametric distribution models to overcome the drawbacks of the parametric models. One of these drawbacks is the need of determining the values of the parameters appropriately by various types of estimation methods. According to [16,17,28,29,30], the Kernal density method is the well-known used method.
Wind energy exploration in Jordan began in 1979 via the Royal Scientific Society (RSS) [31], and various studies have been performed since then to analyze wind energy resources [31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48]. The technical details of eight studies published between 2007 and 2022 as well as the proposed study are listed in Table 1. Markedly, the available span for wind data is between one and 20 years; nevertheless, there was an inverse relationship between the length of the study and the precision of the data available, as shown in the table. These studies were reported to be restricted to the implementation of the Weibull or Rayleigh distribution models. Moreover, it is noticeable that numerical methods dominated the estimation of model parameters in these studies.
Different wind locations in Jordan were examined as study objects in this work to investigate various types of distribution functions for wind energy evaluation. The following were the primary objectives of this study in comparison to those presented by other studies in the same scope:
  • In this paper, the mathematical model of wind analysis based on the Gamma distribution, which was infrequently utilized in the literature, has been derived;
  • Using real examples, this paper provides technical comparisons between several types of wind distribution models: Weibull, Gamma, and Rayleigh distribution functions;
  • The wind energy has been estimated for the designated wind sites using Weibull, Gamma, and Rayleigh distribution models.
The body of this article is structured as follows: Section 2 provides the measured wind data. The direction of the collected wind data is also presented in Section 2. The estimation of these wind data based on Weibull, Rayleigh, and Gamma distribution functions and their parameters is presented in Section 3. The results are discussed in Section 4 before concluding in Section 5.

2. Wind Data

The wind data should be assembled from the nearest meteorological stations to the potential wind site in order to perform an accurate wind resource evaluation. Such information is typically collected at 10 m above ground level and can be extrapolated to a higher altitude (usually the turbine’s installation height). This gives a precise evaluation of the wind for big wind projects. Key wind statistics, such as average wind speed and direction, may be categorized into 10 min, hourly, daily, monthly, and annual time intervals. The most typical time scale, a 10 min basis, can offer enough resolution to provide an accurate wind potential estimation [3,49,50,51,52]. Jordan is divided geographically into three main regions: Jordan Valley region, Western Highlands Range, and the Desert Regions. Jordan Valley is a narrow region along the western strip considered a low-lying agricultural land that can reach up to 400 m below sea level; thus, it is a poor region in terms of wind energy resources. Western Highlands Range is a mountain range along Jordan’s western facade representing an encouraging environment for investment in wind energy projects, as well as the Desert Regions which are vast areas including the eastern and southern parts of Jordan [31,38,39,40]. The wind data were collected and recorded for this study, at a height of 10 m above ground, from nine wind stations in Jordan: Queen Alia Airport, Civil Amman Airport, King Hussein Airport, Irbid, Mafraq, Ma’an, Ghor El-Safi, Safawi, and Irwaished. The geographical distribution of these sites across Jordan is depicted in Figure 1. In this work, the wind data were collected for a one-year period, sufficiently representing all the year’s seasons. The geographic position of the meteorological stations at these sites are shown in Table 2. Additionally, Table 1 includes the wind data time spans, data sampling rate, and the number of obtained data for each wind station.

2.1. Wind Speed Analysis

The direct statistical analysis of wind speed data to find the mean value and standard deviation are illustrated in Figure 2. This chart demonstrates that the average wind speed at all sites was roughly 5 m/s. A deeper look at this data reveals that the average wind speed for all sites ranged between 2.5 m/s and slightly above 6 m/s. The minimum standard deviation regarding wind speed data across all locations was approximately 1 m/s in Irbid, while the maximum values are around 2.7 m/s in Ghor El-Safi. The distribution patterns of wind data are illustrated in Figure 3. Wind data skewness was positive at all locations. Ghor El-Safi had the highest skewness value for wind speed data, with a value of about 12.5. This result shows that its distribution curve had a noticeable asymmetry. The smallest kurtosis values were found at King Hussein Airport and Safawi, with a value of less than one for both sites. Remarkably, Ghor El-Safi had the highest kurtosis value compared to the other sites, indicating a significant difference. The distribution curve for this site was very steep compared to the typical distribution as a result.
The wind speed was split across many divisions in order to be organized. Each division had a speed interval of 1 m/s. The observed probability of each class was calculated for all sites, as shown in Table 3. A closer inspection of this table shows that class (3) recorded the highest probability value in five wind sites, followed by class (4) with two wind sites. The highest probability values achieved by classes (2) and (3) were in the Irbid wind site with around 84%. This indicates that the prevailing wind speed in the Irbid site is within these classes, which is considered as low wind speeds. The same observation was made for the Ghor El-Safi wind site. It is clear from Table 3 that the classes of (11) and above recorded the lowest probability values in all wind sites, as they did not exceed 5% (see the Safawi wind site). Figure 4 represents the histograms of the observed probabilities for two wind sites, one was classified as richly endowed with wind resources (Safawi site) while the other was a resource-poor site, such as Ghor El Safi, according to Table 2 and Figure 2 and Figure 3. In general, the data shown in Table 2 and Table 3 show significant disparities in the features and specifications of the proposed wind sites. This was due to the differences in geography between these areas.

2.2. Wind Direction Analysis

The available data were utilized to determine the wind directions. The frequency of each wind direction at each wind station is well described in Table 4. This table clearly shows that the north-west direction showed the highest recurrence rate across all sites (26.29%). With a mere 3.87% chance of occurring, the eastward wind had the lowest frequency of any direction. Surprisingly, King Hussein Airport had the largest percentage of times that wind was coming from the north (68.55%). In general, knowing the predominant wind direction at a particular wind site gave an indicator of the wind turbine’s proper direction. In most locations, their range should be from northwest to southwest. The wind roses were drawn for both Safawi and Ghor El-Safi sites, as presented in Figure 5.

3. Methodology

In this section, a comprehensive description of the developed models and methods is provided. The work plan of the suggested wind energy prediction is represented in step-wise-step style, graphically, as shown in Figure 6.

3.1. Wind Energy Fundamentals

The time-rate of the flow of wind energy across certain space is called wind power, which is a function of air’s volume, speed and mass. The fundamental equation for wind power can be written as follows [53,54].
P W =   1 2   ρ   A   v 3
where v is the wind speed, ρ is the air density and A is the cross-sectional area of wind flow. This equation clearly shows that wind speed is the most important factor in wind energy calculations. When the wind speed doubled, the available power increased eightfold. Furthermore, wind power is directly related to air density and implicitly related to radius squared.

3.2. Wind Distribution Models

Three distinct distribution functions were employed in this study to illustrate the wind behavior of the specified sites. Many wind researchers are familiar with and employ Weibull and Rayleigh. The Gamma function, which is infrequently applied in the field of wind, was chosen as the third function in this study.

3.2.1. Weibull Model

The Weibull model’s probability distribution function is described as follows: [55,56,57].
f W ( v ) = k c ( v c ) k 1 e ( v c ) k
where k and c are the Weibull shape and scale parameters, respectively. Integrating Equation (2), the cumulative distribution of the Weibull function, Fw(v), is given by:
F W ( v ) = 1 e ( v c ) k

3.2.2. Rayleigh Model

The probability distribution function fR(v) of the Rayleigh model is stated as follows [3,58].
f R ( v ) = v c 2 e ( v 2 2 c 2 )

3.2.3. Gamma Model

Gamma distribution, fG(v), can be used to evaluate wind energy. Gamma’s probability distribution function is written as follows [59].
f G ( v ) = ( v c ) k e ( v c ) Γ ( k )   v
where Γ ( a ) is the gamma function.

3.3. Wind Distribution Models Parameters Estimation

The distribution models were used to represent the probability distribution functions of a set of wind speed data for the candidate wind sites. Each distribution model had several parameters which were estimated by numerical and artificial intelligent (AI) techniques. The numerical techniques used in this study were the maximum likelihood method and the moment method, while the artificial intelligent technique was the whale optimization. The mathematical formulas of the numerical methods utilized in this study are presented in [59]. The whale optimization algorithm will be discussed in the following subsection.

Whale Optimization Algorithm

Mirjilalili et al. proposed the whale optimization algorithm in 2016 [60]. The bubble-net feeding approach, a distinguishing social behavior of humpback whales when foraging, served as the model for this algorithm. By forming several bubbles along a circular or 9-shaped path, this approach works well. The statistical foraging behavior of the humpback whale is demonstrated using the procedures below.
Step one: Encircling the prey
The humpback whale can properly identify and encircle its prey. The location of the prey is considered to be the best option. The other whales (search agents) must, therefore, change their positions in accordance with the position of the ideal whale. As shown below, encircling behavior can be represented by:
D ¯ = | C ¯ · X * ¯ ( t ) X ¯ ( t ) |
X ¯ ( t + 1 ) = X * ¯ ( t ) A ¯ · D ¯
where X * ¯ refers to the best position vector of the optimal solution, X ¯ represents the position vector of a search agent, and the current iteration is defined by t. A ¯ and C ¯ are both coefficient vectors that can be expressed by the following equations:
A ¯ = 2 a ¯ · r ¯ a ¯
C ¯ = 2 r ¯
where a ¯ is a component that decreases linearly from 2 to 0 throughout the iterations, r ¯ is a random vector defined as r   ϵ   [ 0 , 1 ] .
Step two: Method of bubble-net attacking
Two approaches represent this method, as follows.
A. Shrinking encircling mechanism
The mechanism of this approach lowers the value of a ¯ . Consequently, this results in a reduction in the A ¯ coefficient’s fluctuation range, i.e., A ¯   [ a , a ] . Setting the range of A ¯ = [ 1 , 1 ] or | A | < 1 leads to locate the new position of a search agent anywhere in-between the original position and the optimal current position.
B. Spiral updating position
In this approach, the whale attacks the prey in a helix-shaped path. This can be expressed as follows.
D ¯ = | X * ¯ ( t ) X ¯ ( t ) |
X ¯ ( t + 1 ) = D ¯ · e b l · cos ( 2 π l ) + X * ¯ ( t )
where D ¯ represents the distance of the i t h whale to the prey (optimal solution obtained), l is a random number defined as l   ϵ   [ 1 , 1 ] , b is constant value related to the shape of the logarithmic spiral.
It is noteworthy that the humpback whale employs both strategies simultaneously when attacking and encircling its prey. The algorithm, therefore, considers there to be an equal chance of selecting either of the two strategies. As a result, the updated position of the whales can be represented using the following model.
X ¯ ( t + 1 ) = { E q u a t i o n   ( 6 ) ,                                       p < 0.5 E q u a t i o n   ( 11 ) ,                                       p 0.5
where p represents the probability and is a random number defined as p   ϵ   [ 0 , 1 ] .
Step three: Searching for prey
Actually, humpback whales randomly search for prey based on where they are in relation to one another. To encourage the whales (search agents) to separate in pursuit of better prey, a coefficient might be utilized. Making the vector A’s norm bigger than 1 or less than −1 might be used to model this, i.e., | A | > 1 . The new position is obtained at this stage by selecting a random search agent rather than the optimal agent selected. This can be stated as follows.
D ¯ = | C ¯ · X r a n d ¯ ( t ) X ¯ ( t ) |
X ¯ ( t + 1 ) = X r a n d ¯ ( t ) A ¯ · D ¯
where X r a n d ¯ is the random position vector.
The following steps summarize the whale optimization algorithm.
  • Establish the necessary parameters (N, Population size, Itermax) and then, initialize the population Xi (i = 1, 2..., N), as well as the coefficients a, A, C, l, and p;
  • Assess the fitness of each search agent, and then choose X* as the ideal candidate;
  • Update the following coefficients a , A , C , l and p ;
  • Determine the p value. (I) If p < 0.5, then determine the | A | value. (i) If | A | < 1 , update the position by (7). (ii) Otherwise, if | A | 1 , select a random search agent Xrand and then, update the position by (14). (II) Otherwise, if p > 0.5, then update the position by (11);
  • Verify that all whales (search agents) are taken into account. If not, move on to the next search agent; if yes, determine which search agents go over the search space and make the appropriate adjustments;
  • Calculate the fitness for all search agents;
  • Save the best solution X*.;
  • Verify that the stopping criteria are met. If not, move on to step three; if yes, provide the optimal solution X* and its corresponding fitness score.

3.4. Performance Indicators

The efficient distribution model was identified using three statistical performance tests. These tests were: RMSE, R2, and MAE. Each test is discussed in the following subsections.

3.4.1. Root Mean Square Error

The discrepancy between the estimated values of the distribution model and the actual values was computed using the root mean square error (RMSE). Therefore, RMSE is defined as follows [12]:
RMSE = 1 n i = 1 n ( y i x i ) 2
where the y i   a n d   x i are ith measured and estimated value, respectively, while n is the total number of observations. The minimum RMSE value shows that the estimating technique utilized achieved the best results.

3.4.2. Coefficient of Determination

This test assesses the level of agreement between observable and theoretical values predicted by distribution models. R2 is expressed as follows [61,62,63]:
R 2 = 1 i = 1 n ( y i x i ) 2 i = 1 n ( y i y ¯ ) 2
where the y ¯ represents the average value of the observations. The greater the R2 value, the stronger the performance of the estimating approach.

3.4.3. Mean Absolute Error

This test evaluates the absolute difference between the actual and expected values. Corresponding to RMSE, the lower the MAE value, the better the result. MAE is expressed as follows [64,65,66]:
MAE = 1 n i = 1 n | y i x i |

3.5. Wind Energy Estimation

In this section, the Weibull, Gamma, and Rayleigh distribution models were used to evaluate the energy potential of wind regimes. This evaluation was based on an investigation of the two primary components of any wind regime: total annual wind energy (ET) and wind power density (ED).
E D = 0 P v f ( v )     d v
where Pv represents the available wind power per unit area of the wind regime expressed in W/m2, while f(v) is the distribution function. The wind power can be expressed as [66,67]:
P v =   1 2   ρ     v 3
where ρ is the air density in kg/m3. In terms of kWh/m2, the total energy for a certain wind regime at a definite time (T) can be calculated as follows:
E T = T E D
The wind power density of Weibull and Rayleigh functions, (EDW)and (EDR,), can be expressed by (21) and (22), respectively [44].
E D W = ρ c 3 2 Γ ( 1 + 3 k )
E D R = ρ c 3 2 π 2
By substituting (5) and (19) in (18), the wind power density of the gamma function was calculated.
E D G = ρ 2 Γ ( k ) 0 v 2 ( v c ) k e ( v c ) d v
Let x = v/c, therefore dv = cdx. Substitute in (23), yields:
E D G = ρ c 3 2 Γ ( k ) 0 x k + 2 e x d x
The standard Gamma integral is stated as follows:
Γ n = 0 x n 1 e x d x
Comparing (24) with (25), hence, using the Gamma distribution function, the wind power density can be written as follows:
E D G = ρ c 3 Γ ( k + 3 ) 2 Γ ( k )

3.6. Objective Function and Constraint

Each distribution model included several parameters, which can be determined using different estimation methods. Two numerical methods and one Metaheuristic algorithm were applied to determine the parameters of these distribution functions. By reducing the error between the measured wind speed and the predicted values, the objective function utilized in this work aimed to maximize the estimated energy of the wind regime. The objective function in this study and its constraint can be written in a compact format by the following representations [19]:
Objective function:
Maximize   ( E T = T E D )
Subject to:
Error ( v i ) = 1 2 i = 1 n [ f m ( v i ) f c ( v i , i ) ] 2 < ε ο
where fm(vi) represents the measured frequency distribution of ith wind data, fc (vi, фi) is the predicted values obtained by the distribution functions which include the parameters фi, and n is the number of wind speed data.

4. Results and Discussion

In this study, two numerical estimating methods (MM and MLM) and one artificial intelligent optimization (WOA) were used to estimate the parameter values for all indicated models, i.e., Weibull, Rayleigh, and Gamma, for the prospective wind sites. As was already indicated, the goal of this study was to minimize the error between the likelihood of the observed wind speed and the estimated values while also maximizing the estimated energy of the wind regime. The accuracy of the anticipated wind speed in comparison to the measured values offered by Taqs Alarab Weather Forecasting Company were tested using three different performance indicators. The estimated parameters and accuracy indicators for each approach are listed in Table 5. Generally, the WOA method outperformed the traditional methods (MM and MLM) in all sites based on all indicators results; therefore, this supports the fact that artificial intelligent methods are better than traditional methods. Specifically, the Gamma model employing the whale optimization algorithm (G-WOA) achieved the best results among other models (W-WOA and R-WOA), according to RMSE, R2, and MAE outcomes. Amman Civil Airport had the lowest RMSE value, which was 0.01013, while the worst value recorded by Ma’an with 0.03222. The findings of the RMSE test, which demonstrated the dominance of G-WOA over the other approaches, were supported by the MAE test values. A closer look at the accuracy analysis results revealed that for all wind sites, values greater than 0.9 for the results generated by R2 techniques imply that the observations and the expected PDFs were very consistent. In contrast to other wind locations, where the accuracy of the Weibull model employing the whale optimization algorithm (W-WOA) was the superlative model and the indicator values of the RMSE, R2, and MAE tests were, respectively, 0.0176, 0.934, and 0.0093, the results for King Hussein Airport were different. In Figure 7, the accuracy indicator results of the best model (G-WOA) are represented graphically. Figure 8 represents the probability distribution function (PDF) and the observation histograms for each site. These graphs confirm that G-WOA offered the best fit. In comparison to the other methods, the figure generally showed that the PDFs associated with the G-WOA method matched the observed data histograms efficiently.
Table 6 summarizes the entire scenarios of wind energy outcomes for wind regime in all sites. It is noteworthy that the total wind energy calculated by WOA models were greater than the MM and MLM related models; this was because of the optimal performance of WOA according to accuracy indicators. Figure 9 illustrates the total available wind energy based on the G-WOA model. As can be seen from the figure, the Safawi and King Hussein Airport wind sites both produced high total wind energy values with 1250.82 kWh/m2 and 1140.65 kWh/m2, respectively. On the other hand, Irbid and Ghor El-Safi obtained the lowest energy values with 98.51 kWh/m2 and 168.71 kWh/m2, respectively. The mean value and the distribution of wind speeds, effectively, affected the outcomes of energy assessment. It is clear that if the mean value was high and the distribution of wind was closer to symmetry (i.e., skewness and kurtosis approach to zero), the amount of extracted energy from wind regime was high and vice versa. Both the Safawi and Irbid cases support this claim, evidently.

5. Conclusions

In this work, a novel method for evaluating the wind energy potential of nine wind sites in Jordan—Queen Alia Airport, Civil Amman Airport, King Hussein Airport, Irbid, Mafraq, Ma’an, Ghor El-Safi, Safawi, and Irwaished—was provided. The Weibull, Rayleigh, and Gamma distribution functions were used to analyze the wind energy potential. The estimation of wind power density and total wind energy—two essential wind regime parameters—formed the basis of this evaluation. It is possible to determine whether or not a wind farm may be built at a particular wind site by applying this assessment to any wind site in the world.
Three distribution models—the Weibull, Rayleigh, and Gamma distributions—were used to implement the available wind speed data in order to evaluate the created method’s correctness. Estimating the parameters for each distribution model involved using the maximum likelihood and moment methods. The root mean square error, coefficient of determination, and mean absolute error were three performance measures that were examined in order to select the best distribution model. The total wind energy varied from slightly less than 100 kWh/m2 to approximately 1250 kWh/m2. Particularly, the King Hussein Airport and Safawi wind sites both attained the maximum wind power density, whereas Irbid and Ghor El-Safi recorded the lowest wind power densities. Furthermore, the sites with the greatest estimated wind energy showed the highest average wind speed and the most symmetrical distribution pattern.

Author Contributions

Conceptualization, A.A.-Q. and B.A.-M.; methodology, A.A.-Q. and B.A.-M.; software, A.A.-Q. and B.A.-M.; validation, A.A.-Q. and B.A.-M.; formal analysis, A.A.-Q. and B.A.-M.; investigation, A.A.-Q. and B.A.-M.; resources, A.A.-Q. and B.A.-M.; writing—original draft preparation, A.A; writing—review and editing, A.A and B.A; visualization, A.A.-Q., B.A.-M., A.M.A.M., A.R. and H.M.K.A.-M.; supervision, A.A.-Q.; project administration, A.A; funding acquisition, A.A.-Q., A.M.A.M., A.R. and H.M.K.A.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The required data is included in the paper.

Acknowledgments

The authors acknowledge Yarmouk University and Taqs Alarab Weather Forecasting Company for their support in this research study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoekman, S.K.; Broch, A.; Robbins, C.; Ceniceros, E.; Natarajan, M. Review of biodiesel composition, properties, and specifications. Renew. Sustain. Energy Rev. 2012, 16, 143–169. [Google Scholar] [CrossRef]
  2. REN21. Renewables 2020 Global Status Report; REN21 Secretariat: Paris, France, 2020; ISBN 978-3-948393-00-7. Available online: http://www.ren21.net/gsr-2020/ (accessed on 15 January 2021).
  3. Sathyajith, M. Wind Energy: Fundamentals, Resource Analysis and Economics, HAR/CDR ed.; Springer Medizin: Heidelberg, Germany, 2006. [Google Scholar]
  4. Boudia, S.M.; Benmansour, A.; Hellal, M.A.T. Wind resource assessment in Algeria. Sustain. Cities Soc. 2016, 22, 171–183. [Google Scholar] [CrossRef]
  5. Bidaoui, H.; El Abbassi, I.; El Bouardi, A.; Darcherif, A. Wind Speed Data Analysis Using Weibull and Rayleigh Distribution Functions, Case Study: Five Cities Northern Morocco. Procedia Manuf. 2019, 32, 786–793. [Google Scholar] [CrossRef]
  6. Salvação, N.; Soares, C.G. Wind resource assessment offshore the Atlantic Iberian coast with the WRF model. Energy 2017, 145, 276–287. [Google Scholar] [CrossRef]
  7. Bilir, L.; Imir, M.; Devrim, Y.; Albostan, A. Seasonal and yearly wind speed distribution and wind power density analysis based on Weibull distribution function. Int. J. Hydrogen Energy 2015, 40, 15301–15310. [Google Scholar] [CrossRef]
  8. Paraschiv, L.-S.; Paraschiv, S.; Ion, I.V. Investigation of wind power density distribution using Rayleigh probability density function. Energy Procedia 2019, 157, 1546–1552. [Google Scholar] [CrossRef]
  9. Shoaib, M.; Siddiqui, I.; Rehman, S.; Khan, S.; Alhems, L.M. Assessment of wind energy potential using wind energy conversion system. J. Clean. Prod. 2019, 216, 346–360. [Google Scholar] [CrossRef]
  10. Azad, K.; Rasul, M.; Halder, P.; Sutariya, J. Assessment of Wind Energy Prospect by Weibull Distribution for Prospective Wind Sites in Australia. Energy Procedia 2019, 160, 348–355. [Google Scholar] [CrossRef]
  11. Wang, J.; Huang, X.; Li, Q.; Ma, X. Comparison of seven methods for determining the optimal statistical distribution parameters: A case study of wind energy assessment in the large-scale wind farms of China. Energy 2018, 164, 432–448. [Google Scholar] [CrossRef]
  12. Gugliani, G.; Sarkar, A.; Ley, C.; Mandal, S. New methods to assess wind resources in terms of wind speed, load, power and direction. Renew. Energy 2018, 129, 168–182. [Google Scholar] [CrossRef]
  13. Wang, J.; Hu, J.; Ma, K. Wind speed probability distribution estimation and wind energy assessment. Renew. Sustain. Energy Rev. 2016, 60, 881–899. [Google Scholar] [CrossRef]
  14. Li, Y.; Huang, X.; Tee, K.F.; Li, Q.; Wu, X.-P. Comparative study of onshore and offshore wind characteristics and wind energy potentials: A case study for southeast coastal region of China. Sustain. Energy Technol. Assessments 2020, 39, 100711. [Google Scholar] [CrossRef]
  15. Saeed, M.A.; Ahmed, Z.; Yang, J.; Zhang, W. An optimal approach of wind power assessment using Chebyshev metric for determining the Weibull distribution parameters. Sustain. Energy Technol. Assessments 2019, 37, 100612. [Google Scholar] [CrossRef]
  16. Wahbah, M.; Feng, S.F.; El-Fouly, T.H.; Zahawi, B. Wind speed probability density estimation using root-transformed local linear regression. Energy Convers. Manag. 2019, 199, 111889. [Google Scholar] [CrossRef]
  17. Boro, D.; Thierry, K.; Kieno, F.P.; Bathiebo, J. Assessing the Best Fit Probability Distribution Model for Wind Speed Data for Different Sites of Burkina Faso. Curr. J. Appl. Sci. Technol. 2020, 71–83. [Google Scholar] [CrossRef]
  18. Jiang, H.; Wang, J.; Wu, J.; Geng, W. Comparison of numerical methods and metaheuristic optimization algorithms for estimating parameters for wind energy potential assessment in low wind regions. Renew. Sustain. Energy Rev. 2017, 69, 1199–1217. [Google Scholar] [CrossRef]
  19. Alrashidi, M.; Rahman, S.; Pipattanasomporn, M. Metaheuristic optimization algorithms to estimate statistical distribution parameters for characterizing wind speeds. Renew. Energy 2019, 149, 664–681. [Google Scholar] [CrossRef]
  20. Pishgar-Komleh, S.; Keyhani, A.; Sefeedpari, P. Wind speed and power density analysis based on Weibull and Rayleigh distributions (a case study: Firouzkooh county of Iran). Renew. Sustain. Energy Rev. 2015, 42, 313–322. [Google Scholar] [CrossRef]
  21. Serban, A.; Paraschiv, L.S.; Paraschiv, S. Assessment of wind energy potential based on Weibull and Rayleigh distribution models. Energy Rep. 2020, 6, 250–267. [Google Scholar] [CrossRef]
  22. Hu, Q.; Wang, Y.; Xie, Z.; Zhu, P.; Yu, D. On estimating uncertainty of wind energy with mixture of distributions. Energy 2016, 112, 935–962. [Google Scholar] [CrossRef]
  23. Dong, Y.; Wang, J.; Jiang, H.; Shi, X. Intelligent optimized wind resource assessment and wind turbines selection in Huitengxile of Inner Mongolia, China. Appl. Energy 2013, 109, 239–253. [Google Scholar] [CrossRef]
  24. Wu, J.; Wang, J.; Chi, D. Wind energy potential assessment for the site of Inner Mongolia in China. Renew. Sustain. Energy Rev. 2013, 21, 215–228. [Google Scholar] [CrossRef]
  25. Akgül, F.G.; Şenoğlu, B.; Arslan, T. An alternative distribution to Weibull for modeling the wind speed data: Inverse Weibull distribution. Energy Convers. Manag. 2016, 114, 234–240. [Google Scholar] [CrossRef]
  26. Aries, N.; Boudia, S.M.; Ounis, H. Deep assessment of wind speed distribution models: A case study of four sites in Algeria. Energy Convers. Manag. 2018, 155, 78–90. [Google Scholar] [CrossRef]
  27. Jung, C.; Schindler, D. Wind speed distribution selection—A review of recent development and progress. Renew. Sustain. Energy Rev. 2019, 114. [Google Scholar] [CrossRef]
  28. Chang, T.P. Estimation of wind energy potential using different probability density functions. Appl. Energy 2011, 88, 1848–1856. [Google Scholar] [CrossRef]
  29. Qin, Z.; Li, W.; Xiong, X. Estimating wind speed probability distribution using kernel density method. Electr. Power Syst. Res. 2011, 81, 2139–2146. [Google Scholar] [CrossRef]
  30. Xu, X.; Yan, Z.; Xu, S. Estimating wind speed probability distribution by diffusion-based kernel density method. Electr. Power Syst. Res. 2015, 121, 28–37. [Google Scholar] [CrossRef]
  31. Ammari, H.D.; Al-Rwashdeh, S.S.; Al-Najideen, M.I. Evaluation of wind energy potential and electricity generation at five locations in Jordan. Sustain. Cities Soc. 2015, 15, 135–143. [Google Scholar] [CrossRef]
  32. Habali, S.; Hamdan, M.; Jubran, B.; Zaid, A.I. Wind speed and wind energy potential of Jordan. Sol. Energy 1987, 38, 59–70. [Google Scholar] [CrossRef]
  33. Anani, A.; Zuamot, S.; Abu-Allan, F.; Jibril, Z. Evaluation of wind energy as a power generation source in a selected site in Jordan. Sol. Wind. Technol. 1988, 5, 67–74. [Google Scholar] [CrossRef]
  34. Amr, M.; Petersen, H.; Habali, S. Assessment of wind farm economics in relation to site wind resources applied to sites in Jordan. Sol. Energy 1990, 45, 167–175. [Google Scholar] [CrossRef]
  35. Habali, S.; Amr, M.; Saleh, I.; Ta’Ani, R. Wind as an alternative source of energy in Jordan. Energy Convers. Manag. 2001, 42, 339–357. [Google Scholar] [CrossRef]
  36. Abderrazzaq, M. Energy production assessment of small wind farms. Renew. Energy 2004, 29, 2261–2272. [Google Scholar] [CrossRef]
  37. Alghoul, M.A.; Sulaiman, M.Y.; Azmi, B.Z.; Wahab, M.A. Wind energy potential of Jordan. Int. Energy J. 2007, 8, 71–78. [Google Scholar]
  38. Bataineh, K.M.; Dalalah, D. Assessment of wind energy potential for selected areas in Jordan. Renew. Energy 2013, 59, 75–81. [Google Scholar] [CrossRef]
  39. Alsaad, M. Wind energy potential in selected areas in Jordan. Energy Convers. Manag. 2013, 65, 704–708. [Google Scholar] [CrossRef]
  40. Dalabeeh, A.S.K. Techno-economic analysis of wind power generation for selected locations in Jordan. Renew. Energy 2017, 101, 1369–1378. [Google Scholar] [CrossRef]
  41. Al-Ghriybah, M.; Fadhli, Z.M.; Hissein, D.D.; Mohd, S. Wind energy assessment for the capital city of Jordan, Amman. J. Appl. Eng. Sci. 2019, 17, 311–320. [Google Scholar] [CrossRef] [Green Version]
  42. Marashli, A.; Alburdaini, M.; Al-Rawashdeh, H.; Shalby, M. Statistical Analysis of Wind Speed Distribution Based on Five Weibull Methods for Wind Power Evaluation in Maan, Jordan. J. Energy Technol. Policy 2021, 11, 55–70. [Google Scholar] [CrossRef]
  43. Alsaqoor, S.; Marashli, A.; At-Tawarah, R.; Borowski, G.; Alahmer, A.; Aljabarin, N.; Beithou, N. Evaluation of Wind Energy Potential in View of the Wind Speed Parameters—A Case Study for the Southern Jordan. Adv. Sci. Technol. Res. J. 2022, 16, 275–285. [Google Scholar] [CrossRef]
  44. Al-Mhairat, B.; Al-Quraan, A. Assessment of Wind Energy Resources in Jordan Using Different Optimization Techniques. Processes 2022, 10, 105. [Google Scholar] [CrossRef]
  45. Al-Quraan, A.; Al-Mhairat, B. Intelligent Optimized Wind Turbine Cost Analysis for Different Wind Sites in Jordan. Sustainability 2022, 14, 3075. [Google Scholar] [CrossRef]
  46. Al-Quraan, A.; Al-Masri, H.; Al-Mahmodi, M.; Radaideh, A. Power curve modelling of wind turbines-A comparison study. IET Renew. Power Gener. 2021, 16, 362–374. [Google Scholar] [CrossRef]
  47. Al-Quraan, A.; Al-Masri, H.; Radaideh, A. New Method for Assessing the Energy Potential of Wind Sites-A Case Study in Jordan. Univers. J. Electr. Electron. Eng. 2020, 7, 209–218. [Google Scholar] [CrossRef]
  48. Al-Quraan, A.; Al-Mahmodi, M.; Radaideh, A.; Al-Masri, H.M.K. Comparative study between measured and estimated wind energy yield. Turk. J. Electr. Eng. Comput. Sci. 2020, 28, 2926–2939. [Google Scholar] [CrossRef]
  49. Elsner, P. Continental-scale assessment of the African offshore wind energy potential: Spatial analysis of an under-appreciated renewable energy resource. Renew. Sustain. Energy Rev. 2019, 104, 394–407. [Google Scholar] [CrossRef] [Green Version]
  50. Han, Q.; Ma, S.; Wang, T.; Chu, F. Kernel density estimation model for wind speed probability distribution with applicability to wind energy assessment in China. Renew. Sustain. Energy Rev. 2019, 115, 109387. [Google Scholar] [CrossRef]
  51. Kaplan, Y.A. Overview of wind energy in the world and assessment of current wind energy policies in Turkey. Renew. Sustain. Energy Rev. 2015, 43, 562–568. [Google Scholar] [CrossRef]
  52. Kazimierczuk, A.H. Wind energy in Kenya: A status and policy framework review. Renew. Sustain. Energy Rev. 2019, 107, 434–445. [Google Scholar] [CrossRef]
  53. Al-Quraan, A.; Al-Qaisi, M. Modelling, Design and Control of a Standalone Hybrid PV-Wind Micro-Grid System. Energies 2021, 14, 4849. [Google Scholar] [CrossRef]
  54. Radaideh, A.; Bodoor, M.; Al-Quraan, A. Active and Reactive Power Control for Wind Turbines Based DFIG Using LQR Controller with Optimal Gain-Scheduling. J. Electr. Comput. Eng. 2021, 2021, 1–19. [Google Scholar] [CrossRef]
  55. Ko, D.H.; Jeong, S.T.; Kim, Y.C. Assessment of wind energy for small-scale wind power in Chuuk State, Micronesia. Renew. Sustain. Energy Rev. 2015, 52, 613–622. [Google Scholar] [CrossRef]
  56. Ladenburg, J.; Hevia-Koch, P.; Petrović, S.; Knapp, L. The offshore-onshore conundrum: Preferences for wind energy considering spatial data in Denmark. Renew. Sustain. Energy Rev. 2020, 121, 109711. [Google Scholar] [CrossRef]
  57. Peters, J.L.; Remmers, T.; Wheeler, A.J.; Murphy, J.; Cummins, V. A systematic review and meta-analysis of GIS use to reveal trends in offshore wind energy research and offer insights on best practices. Renew. Sustain. Energy Rev. 2020, 128, 109916. [Google Scholar] [CrossRef]
  58. Everitt, B.; Skrondal, A. The Cambridge Dictionary of Statistics; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  59. Forbes, C.; Evans, M.; Hastings, N.; Peacock, B. Statistical Distributions, 4th ed.; Wiley: Hoboken, NJ, USA, 2010. [Google Scholar]
  60. Mirjalili, S.; Lewis, A. The Whale Optimization Algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  61. Rosales-Asensio, E.; Borge-Diez, D.; Blanes-Peiró, J.-J.; Pérez-Hoyos, A.; Comenar-Santos, A. Review of wind energy technology and associated market and economic conditions in Spain. Renew. Sustain. Energy Rev. 2018, 101, 415–427. [Google Scholar] [CrossRef]
  62. Al-Quraan, A.; Alrawashdeh, H. Correlated Capacity Factor Strategy for Yield Maximization of Wind Turbine Energy. In Proceedings of the IEEE 5th International Conference on Renewable Energy Generation and Applications (ICREGA), Al-Ain, United Arab Emirates, 26–28 February 2018. [Google Scholar]
  63. Al-Quraan, A.; Stathopoulos, T.; Pillay, P. Comparison of wind tunnel and on site measurements for urban wind energy estimation of potential yield. J. Wind. Eng. Ind. Aerodyn. 2016, 158, 1–10. [Google Scholar] [CrossRef] [Green Version]
  64. Brano, V.L.; Orioli, A.; Ciulla, G.; Culotta, S. Quality of wind speed fitting distributions for the urban area of Palermo, Italy. Renew. Energy 2011, 36, 1026–1039. [Google Scholar] [CrossRef]
  65. Al-Quraan, A.; Al-Mahmodi, M.; Al-Asemi, T.; Bafleh, A.; Bdour, M.; Muhsen, H.; Malkawi, A. A New Configuration of Roof Photovoltaic System for Limited Area Applications—A Case Study in KSA. Buildings 2022, 12, 92. [Google Scholar] [CrossRef]
  66. Wimhurst, J.J.; Greene, J.S. Oklahoma’s future wind energy resources and their relationship with the Central Plains low-level jet. Renew. Sustain. Energy Rev. 2019, 115, 109374. [Google Scholar] [CrossRef]
  67. Al-Addous, M.; Saidan, M.; Bdour, M.; Dalala, Z.; Albatayneh, A.; Class, C.B. Key aspects and feasibility assessment of a proposed wind farm in Jordan. Int. J. Low-Carbon Technol. 2019, 15, 97–105. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of wind sites under test.
Figure 1. Location of wind sites under test.
Sustainability 15 03927 g001
Figure 2. Statistical analysis values of wind data including the mean wind speed and standard deviation.
Figure 2. Statistical analysis values of wind data including the mean wind speed and standard deviation.
Sustainability 15 03927 g002
Figure 3. The pattern distribution values of wind speed including the skewness and kurtosis.
Figure 3. The pattern distribution values of wind speed including the skewness and kurtosis.
Sustainability 15 03927 g003
Figure 4. Probability histogram of wind speed distribution for both Safawi (left) and Ghor El-Safi sites (right).
Figure 4. Probability histogram of wind speed distribution for both Safawi (left) and Ghor El-Safi sites (right).
Sustainability 15 03927 g004
Figure 5. Wind rose for both Safawi (top) and Ghor El-Safi (bottom) wind sites.
Figure 5. Wind rose for both Safawi (top) and Ghor El-Safi (bottom) wind sites.
Sustainability 15 03927 g005aSustainability 15 03927 g005b
Figure 6. Flowchart for the procedure of wind energy prediction.
Figure 6. Flowchart for the procedure of wind energy prediction.
Sustainability 15 03927 g006
Figure 7. RMSE and MAE (left) and R2 (right) values of best model (G-WOA).
Figure 7. RMSE and MAE (left) and R2 (right) values of best model (G-WOA).
Sustainability 15 03927 g007
Figure 8. Probability distribution functions for candidate wind sites.
Figure 8. Probability distribution functions for candidate wind sites.
Sustainability 15 03927 g008
Figure 9. The total wind energy for all sites based on best model G-WOA.
Figure 9. The total wind energy for all sites based on best model G-WOA.
Sustainability 15 03927 g009
Table 1. The recent studies of wind energy assessment in Jordan.
Table 1. The recent studies of wind energy assessment in Jordan.
StudyCandidate Site(s)Study PeriodData ResolutionDistribution Model(s)Estimation Method(s)Performance Indicator(s)Study Objective(s)
Study in Ref. [37]AmmanAqaba9 yearsDaily WeibullGraphical methodKolmogorov–Smirnov For wind regime:
  • Annual power density.
IrbidDer Alla
Ras Monief
Study in Ref. [38]HofaTafila2–4 years
9 years (Fujaij)
Monthly RayleighDirect computing of parameter (c) based on mean speed N/AFor proposed 4 models of WTs:
  • Annual energy production;
  • Capacity factor;
  • Cost of energy.
IbrahimyaFujaij
ZabdaAqaba
Ras Monief
Study in Ref. [39]Alhassan Industrial City20 yearsMonthly N/AN/AN/AFor proposed wind farm:
  • Annual energy production;
  • Capacity factor;
  • Cost of energy;
  • Payback period;
  • Annual net profit.
Fujaij
Safawi
Ras Monief
Study in Ref. [31]Azraq Q. A. Airport *5 years6 hWeibullStandard deviation methodN/AFor proposed 5 models of WTs:
  • Annual and seasonal energy pro-duction;
  • Capacity factor.
SafawiK. H. Airport **
Ras Monief
Study in Ref. [40]Amman9 yearsDaily WeibullGraphical methodN/AFor proposed 4 models of WTs:
  • Annual energy production;
  • Capacity factor.
  • No. of WT units;
  • Cost of energy.
Irbid
Aqaba
Der Alla
Ras Monief
Study in Ref. [41]Amman 7 yearsN/AWeibull Trend of k vs. V ¯ N/AFor wind regime:
  • Annual power density;
  • Annual energy density.
Study in Ref. [42]Ma’an 1 yearDailyWeibull
Rayleigh
Graphical method
Empirical method
Moment method
Energy pattern factor method
RMSE
X2 test
R2
MPE
MAPE
For wind regime:
  • Monthly power density.
Study in Ref. [43]Ma’an1 year10 min Weibull
Rayleigh
Empirical method
Moment method
Energy pattern factor method
N/AFor wind regime:
  • Annual power density;
  • Annual energy density. For proposed 4 models of WTs:
  • Annual energy production.
Aqaba
Batn Elghol
Proposed StudyQ. A. Airport *Mafraq1 year1 h, 3 h and 6 hGamma
Weibull
Rayleigh
Artificial intelligent method (WOA)
Moment method
Maximum likelihood method
RMSE
R2
MAE
For wind regime:
  • Annual power density;
  • Annual energy density.
K. H. Airport **Ma’an
A. C. Airport ***Safawi
IrbidIrwaished
Ghor El Safai
* Queen Alia Airport; ** King Hussein Airport; *** Amman Civil Airport.
Table 2. Topographical information and technical details of wind data for each designated wind site.
Table 2. Topographical information and technical details of wind data for each designated wind site.
LatitudeLongitudeElevationNumber of DataSampling RatePeriod
Queen Alia Airport31.43° N35.59° E722 m71731 h01/2019–12/2019
Amman Civil Airport31.59° N35.59° E767 m69101 h09/2018–08/2019
King Hussein Airport29.33° N35.00° E51 m67111 h01/2018–12/2018
Irbid32.33° N35.51° E618 m8996 h03/2018–02/2019
Mafraq32.22° N36.15° E686 m19093 h09/2018–08/2019
Ma’an30.10° N35.47° E1069 m10146 h03/2018–02/2019
Safawi32.09° N37.12° E647 m17033 h01/2018–12/2018
Irwaished32.30° N38.12° E686 m7646 h09/2017–08/2018
Ghor El-Safi31.02° N35.28° E−350 m7906 h01/2018–12/2018
Table 3. The observed probabilities of wind speeds for all sites.
Table 3. The observed probabilities of wind speeds for all sites.
No. of ClassSpeed Class
(m/s)
Wind Speed Observations (%)
Queen Alia AirportAmman Civil AirportKing Hussein AirportIrbidMafraqMa’anSafawiIrwaishedGhor El-Safi
1 0 V i < 1 0.040.040.000.110.000.000.000.130.13
2 1 V i < 2 3.1414.200.5843.605.081.680.296.1534.68
3 2 V i < 3 23.6723.0012.7140.0431.5318.848.6926.0541.14
4 3 V i < 4 16.9418.4716.3611.1220.5332.6424.7820.2913.42
5 4 V i < 5 12.4613.9816.413.6715.6118.9319.5514.276.20
6 5 V i < 6 16.2312.5916.381.2212.0510.0615.9713.093.04
7 6 V i < 7 10.257.7116.870.117.863.559.518.641.01
8 7 V i < 8 8.434.6013.800.004.664.445.583.530.25
9 8 V i < 9 2.892.143.890.111.105.926.463.530.00
10 9 V i < 10 2.271.192.210.000.791.184.171.440.00
11 10 V i < 11 1.740.900.520.000.260.793.410.260.00
12 11 V i < 12 0.700.640.180.000.100.790.700.650.00
13 12 V i < 13 0.500.410.010.000.260.390.530.920.00
14 13 V i < 14 0.250.070.000.000.050.490.120.650.00
15 14 V i 0.490.070.090.000.100.300.230.390.13
Table 4. The occurrence rate of each wind direction for each site.
Table 4. The occurrence rate of each wind direction for each site.
SiteOccurrence Rate (%)Overall
(%)
NNWWSWSSEENE
Queen Alia Airport1.8726.2027.4217.502.178.027.179.66100
Amman Civil Airport4.1529.3221.4926.660.874.675.707.13100
King Hussein Airport68.5516.691.004.043.462.090.363.82100
Irbid3.6714.2433.9325.142.0015.684.121.22100
Mafraq5.6151.237.4416.662.8314.351.200.68100
Ma’an7.5041.3215.3815.988.196.312.273.06100
Safawi14.1522.4918.0324.137.117.693.353.05100
Irwaished7.4620.5515.7116.3612.3011.658.907.07100
Ghor El-Safi17.2214.566.2032.533.424.431.7719.87100
Average14.4626.2916.2919.894.708.323.876.17100
Table 5. Distribution function parameters and performance indicators results.
Table 5. Distribution function parameters and performance indicators results.
SiteParameter and IndicatorEstimation Methods
W-MMW-MLMW-WOAR-MMR-MLMR-WOAG-MMG-MLMG-WOA
Queen Alia AirportK2.3262.0752.0132.0002.0002.0004.7714.1724.051
C6.5645.2725.0124.6413.6973.3691.2191.1141.013
RMSE0.039190.028460.02540.037380.029040.026900.036200.024880.02223
R20.683130.832870.882720.711690.825990.875270.729650.872320.90230
MAE0.021570.015020.013120.022010.015460.014460.017680.014170.01140
Amman Civil AirportK3.6094.2771.5032.0002.0002.0001.9402.1910.541
C5.7504.4924.2314.0633.2014.9031.1901.0981.034
RMSE0.037960.018480.015480.035370.018090.017110.034540.010820.01013
R20.735230.937270.957210.770060.939850.949120.780820.978490.98836
MAE0.025140.009000.007000.023680.009440.008210.021610.007250.00425
King Hussein AirportK3.0112.8062.5152.0002.0002.0007.6366.3395.323
C6.3975.6403.5504.5593.8002.985000.7480.7910.6820
RMSE0.020120.017600.015880.028860.030120.027460.021470.022570.01734
R20.913780.934020.943620.822650.806740.846120.901860.891490.92164
MAE0.010920.009930.007930.021920.016630.015610.012720.012010.01002
IrbidK2.7952.4622.1222.0002.0002.0006.6646.5886.211
C2.7982.3792.1391.9881.6181.5210.3740.3200.282
RMSE0.068940.042260.032260.081080.074090.054090.045040.020860.01883
R20.831090.936530.946310.766380.804910.841230.927900.984530.99123
MAE0.039910.027890.022890.060060.042310.023850.024570.014900.01291
MafraqK2.3612.0901.8902.0002.0002.0004.9044.7574.451
C5.3854.2694.1663.8082.9872.6120.9730.7920.721
RMSE0.040080.030730.026610.039330.032140.030140.034640.023100.02011
R20.727360.839730.869520.737350.824620.844550.796290.909380.92912
MAE0.017180.011690.010690.017510.012460.011460.013110.010620.00805
Ma’anK2.2792.0332.0122.0002.0002.0004.5954.5924.582
C6.1834.8684.5614.3703.4283.1221.1920.9340.832
RMSE0.060740.048120.033220.059820.048690.035110.052560.035120.03222
R20.540470.711680.794510.554390.704780.784220.656020.846390.91539
MAE0.037740.028480.021220.038060.028730.022730.030760.018890.01711
SafawiK2.5662.3122.2122.0002.0002.0005.7045.2765.076
C6.9175.7835.6454.9003.9743.7341.0770.9670.912
RMSE0.045500.033410.031340.044810.038320.028230.037430.023010.02014
R20.630020.800520.820230.641090.737500.857120.749630.905370.91533
MAE0.028770.021520.018120.030460.025360.015010.021980.015150.01112
IrwaishedK2.1501.9081.8152.0002.0002.0004.1353.7973.297
C6.2164.7644.2644.3933.4123.1121.3321.1071.007
RMSE0.048210.032680.030340.046920.031350.031340.042430.022860.01823
R20.621090.825880.835120.640980.839770.839120.706530.914820.93401
MAE0.028690.019410.017010.028540.018590.018010.025170.013880.01312
Ghor El-SafiK1.7291.6761.4762.0002.0002.0002.7794.7574.351
C5.0872.6832.4833.6172.0652.0121.6310.5030.493
RMSE0.061490.041520.038340.064660.032210.028210.058020.014570.01150
R20.485830.765580.735230.431450.858950.918950.542150.971140.98110
MAE0.025010.012030.010120.025920.010910.009910.023630.004810.00388
Table 6. Total Wind Energy for the Selected Wind Sites based MM, MLM, and WOA.
Table 6. Total Wind Energy for the Selected Wind Sites based MM, MLM, and WOA.
ET (kWh/m2)
MethodNumericalAINumericalAINumericalAI
MMMLMWOAMMMLMWOAMMMLMWOA
SiteQueen Alia AirportAmman Civil AirportKing Hussein Airport
Weibull1002.321006.721150.85634.64668.45809.34910.87992.541066.74
Rayleigh1005.821019.061319.76653.56661.49792.661080.981106.811228.51
Gamma930.56988.211070.26622.21662.62862.921002.211029.351140.65
SiteIrbidMa’anMafraq
Weibull76.5580.38109.58488.60530.80630.93780.64809.05915.65
Rayleigh81.2185.47113.77502.45537.65577.55802.88812.81890.41
Gamma70.375.798.51420.22493.92603.73705.23739.93825.12
SiteSafawiIrwaishedGhor El-Safi
Weibull1190.991209.171369.33164.89172.17194.25165.74172.17192.84
Rayleigh1202.651265.831295.95158.92177.7201.22155.56177.70189.52
Gamma1121.431170.331250.82116.22126.36192.71122.76126.36168.71
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al-Quraan, A.; Al-Mhairat, B.; Malkawi, A.M.A.; Radaideh, A.; Al-Masri, H.M.K. Optimal Prediction of Wind Energy Resources Based on WOA—A Case Study in Jordan. Sustainability 2023, 15, 3927. https://0-doi-org.brum.beds.ac.uk/10.3390/su15053927

AMA Style

Al-Quraan A, Al-Mhairat B, Malkawi AMA, Radaideh A, Al-Masri HMK. Optimal Prediction of Wind Energy Resources Based on WOA—A Case Study in Jordan. Sustainability. 2023; 15(5):3927. https://0-doi-org.brum.beds.ac.uk/10.3390/su15053927

Chicago/Turabian Style

Al-Quraan, Ayman, Bashar Al-Mhairat, Ahmad M. A. Malkawi, Ashraf Radaideh, and Hussein M. K. Al-Masri. 2023. "Optimal Prediction of Wind Energy Resources Based on WOA—A Case Study in Jordan" Sustainability 15, no. 5: 3927. https://0-doi-org.brum.beds.ac.uk/10.3390/su15053927

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