Aquila Optimization Algorithm for Wind Energy Potential Assessment Relying on Weibull Parameters Estimation

: Statistical distribution approaches have been developed to describe wind data due to the intermittent and unpredictable nature of wind speed. The Weibull distribution with two parameters is thought to be the most accurate distribution for modeling wind data. This study seeks wind energy assessment via searching for the optimal estimation of the Weibull parameters. For this target, analytical and heuristic methods are investigated. The analytical methods involve the maximum likelihood, moment, energy pattern factor, and empirical methods, while the heuristic optimization algorithms include particle warm optimization and the Aquila optimizer (AO). Both analytical and heuristic methods are assessed together to ﬁt the probability density function of wind data. In addition, nine models are submitted to ﬁnd the most appropriate model to represent wind energy production. The error between actual and estimated wind energy density is computed to the model for each study site which has less error of energy. The ﬁt test is performed with real data for the Zafarana and Shark El-Ouinate sites in Egypt for a year. Additionally, different indicators of ﬁtness properties are assessed, such as the root mean square error, determination coefﬁcient (R 2 ), mean absolute error, and wind production deviation. The simulation results declare that the proposed AO optimization algorithm offers greater accuracy than several optimization algorithms in the literature for estimating the Weibull parameters. Furthermore, statistical analysis of the compared methods demonstrates the high stability of the AO algorithm. Thus, the proposed AO has greater accuracy and more stability in the obtained outcomes for Weibull parameters and wind energy calculations.


Introduction
In recent years, the use of renewable energy resources has continuously increased in many countries worldwide. This is due to the growing popular consciousness of using clean energy, and the electrical sectors have many initiatives in place to diversify energy resources and reduce reliance on traditional fuels [1,2].
Wind energy can be considered one of the potential energy resources which can be used in a commercial way. Wind speed is the main source of wind energy by converting it into electrical power at the rotor of the wind turbine. Because wind speed has an intermittent and variable condition, the study of speed distribution is a vital task in determining the wind energy capacity of the site [3]. Nowadays, wind resources are integrated with distribution power systems [4], cogeneration systems [5], AC-DC power networks [6], and multi-area power systems [7].
Due to the intermittent nature of wind velocity, wind energy is unstable and fluctuating. The uncertainties make it difficult to make certain judgments, such as operational Wind 2022, 2 618 and investment decisions. To overcome operational challenges, wind energy should be combined with other conventional technologies. Furthermore, making more precise use of wind speed distribution might aid in making profitable wind energy investment selections [8]. In this context, many PDFs have been provided in the literature to represent the distribution of wind speed, and they have been compared to identify the most suited one.
The nature of wind speed introduces two major uncertainties, one of which is wind speed forecast uncertainty, which may be handled by an efficient PDF, as demonstrated in several studies. The second uncertainty pertains to wind speed to power conversion, which may be resolved by employing several models to represent the nonlinear portion of the power curve [3,9]. Because wind speed is stochastic, statistical analysis is critical for dealing with uncertainties in two issues: wind speed forecasting and wind output power. The Weibull distribution is a statistical technique that is widely used and recommended for calculating the PDF of wind speed [10].
Aside from the Weibull distribution, various probabilistic distributions, such as gamma, normal, and log-normal distributions, have been used to describe wind speed data and transform it from discrete values to continuous functions [11]. In addition, detailed comparisons of the Weibull and various distribution functions are presented in ref. [12]. The characterization of the wind regime may be determined from a probability distribution of wind speed, which is regarded as the main stage in the construction of wind power projects as well as making investment decisions.
The two major parameters of the Weibull distribution function are utilized to fit measured data at the wind farm location. Many researchers have presented analytical methods for estimating parameters process such as in refs. [10,[13][14][15]. In ref. [10], seven analytical methods were used for determining the best Weibull parameters for wind energy generation in the northeast region of Brazil, which are the graphical method (GM), maximum likelihood method (MLM), moment method (MM), energy pattern factor method (EPFM), empirical method (EM), modified maximum likelihood method (MMLM), and equivalent energy method (EEM). The Weibull distribution of the Zafarana projects in the Suez Gulf, Egypt, were determined based on monthly averaged data for one year as well as every 10 min for two days and compared with other analytical methods to achieve the best Weibull parameters [13]. Additionally, the power density method (PD) was presented to estimate Weibull distribution parameters and was compared with other analytical methods [14]. Furthermore, a developed method called the wind energy intensification method (WEIM) was presented to calculate the Weibull parameters of 60 locations in Pakistan based on three years (2014-2017) of hourly measured wind data at 50 m height [16]; furthermore, a comparative assessment of the analytical methods was carried out and an energy variance method (EVM) was performed [17].
The Weibull and Rayleigh, which is a special case of the Weibull distribution, have been compared to fit the measured wind speed data at Iskenderun, located in Turkey, and the wind energy potential has been evaluated based on 1-year measured hourly time-series wind speed data [18]. Additionally, the least square method (LSM) has been utilized to calculate the parameters of the Weibull distribution at Kayathar, Tamil Nadu, India [19,20]. In ref. [21], probability-weighted moments based on the power density method were executed as another analytical method and compared to other common methods. In this study, the MM, MLM, LSM, and maximum entropy principal Method (MEPM) were presented and compared as parametric models to estimate Weibull parameters. Additionally, the modified method of moment (MMOM) was provided as a novel methodology and compared to established methods, but not to any heuristic optimization algorithms [22]. There are different statistical tests which have been utilized to compare the performance of each method, such as the determination coefficient (R 2 ), root mean square error (RMSE), mean absolute error (MAE), correlation coefficient (R), standard deviation, wind energy error (WEE), and Monte Carlo simulation [15,23,24]. Where iterative features have been used, the analytical approaches still exhibit their limitations. The wind speed distribution cannot always be represented efficiently due to analytical approaches. As a result, the optimum Weibull distribution parameters may be determined using heuristic optimization methods that employ the goal function to obtain the best fitting values. In this context, the mean squared, RMSE, MAE, and MSE can be used as objective functions to obtain optimal values of Weibull parameters.
On the level of optimization algorithms, one of the heuristic optimization algorithms called social spider optimization has been presented to estimate Weibull parameters for wind speed distribution and compared to MLM, MM, and LSM as analytical methods and particle swarm optimization (PSO), genetic algorithm (GA), and cuckoo optimization algorithm (COA) as heuristic algorithms [8]. Additionally, the single and combined parameter PDFs, which are Weibull, Rayleigh, log-normal, and gamma, were developed to overcome the shortcoming exit in a single PDF. The comparison between some common analytical methods such as MMLM, EM, EPF, MM, and GM and heuristic optimization algorithms such as the grey wolf algorithm (GWO), bat optimization algorithm (BOA), GA, and bee colony algorithm (BCA) have been presented to guarantee the optimal fitting of Weibull parameters. The results of Weibull fitness tests show that heuristic optimization algorithms offer more accuracy than the numerical methods [25,26]. Additionally, the results of a harmony search (HS), COA, PSO, and ant colony optimization (ACO) have been compared to evaluate the performance of each method with those obtained by MLM, MM, EM, and the equivalent energy method (EEM) [27]. An investigated comparison between the analytical methods and heuristic estimation methods was implied [26], and the results show that heuristic estimation outperforms parametric methods in terms of fitting accuracy and operational simplicity. Additionally, the parameters of the Weibull distribution estimated using a simulated annealing algorithm were conducted in ref. [28].
Another uncertainty source is the wind power curve [9] which is produced from the nonlinearity part of the power curve. In the literature, there are several studies that presented many models for the nonlinearity part of the wind power curve. The polynomial, exponential, and cubic models have been addressed to achieve an optimal model for wind energy sites [18,21,[29][30][31].
Analytical methods are difficult to resolve as numerical iterations are needed to determine the parameters for the Weibull distribution. Moreover, the error in adapting analytical methods based on statistical analysis is high when estimating Weibull parameters. Accordingly, achieving the optimal wind PDF speed model is still considered a challenging task. Therefore, in this paper, different optimization algorithms of PSO and AO are developed, employed, and compared with other analytical techniques for handling this task. PSO, similar to GA, is a population-based algorithm. Its population and individuals are called swarm and particles, respectively. Each particle keeps track of its coordinates in the problem space, which is usually associated with the best solution obtained so far [8,11,32]. Additionally, the Aquila optimizer (AO) is a novel meta-heuristic algorithm that was presented in 2020 and implemented for solving several benchmarks [33]. Moreover, the simple structure, easy implementation, searching force, and faster convergence motivated it to be proposed to solve one of the important optimization problems.
The contribution of this work can be summarized as follows: • Analysis of the distribution of real data obtained for Zafarana and Shark El-Ouinate sites, Egypt. • Implementation of MLM, EM, MM, and EPFM as analytical methods for optimal fitting based on R 2 , RMSE, MAE, and wind production deviation (WPD). • Analysis of novel intelligence optimization methods called PSO and AO for measuring the accuracy of fitting by statistical tests.
Wind 2022, 2 620 The paper is described in an organized format with six sections. Section 1 specifies the introduction to the motivation of the research and the existing research that is related to the present. Section 2 demonstrates the analytical and heuristic methods that are used for PDF of wind speed distribution. Section 3 presents the measurement of fitting for all methods. Section 4 provides wind speed data and wind turbine data for Zafarana and Shark El-Ouinate sites and illustrates the comprehensive indexes intended to validate the Weibull parameter, statistical analysis, and convergence characteristics of PSO and AO. Section 5 depicts the different models for wind power curves to find the most appropriate model for each site. Section 6 derives the conclusion of the research.

Wind Speed Distribution Functions
The vital step before starting any wind energy project or making any investment decisions is the analysis of wind speed distribution at the selected site. The wind energy in any site is affected by any fluctuations in wind speed distribution. Therefore, the accuracy of the wind speed distribution model is especially important in assessing the site's wind energy potential. In the literature, there are several density functions, such as Weibull, Rayleigh, log-normal, and gamma, which can be used to describe the wind speed distribution [8]. The Weibull distribution is the most used which has two parameters, while the Rayleigh distribution is a subset of the Weibull distribution, which has only one parameter. The Weibull probability distribution function f (v) for speed (v) is expressed as follows [22,27]: and the Weibull cumulative distribution function F(v) is given by, (2)

Analytical Methods
There are several analytical methods already used to estimate Weibull curves estimating parameters for wind speed analysis, such as MLM, EPFM, MM, and EM [10,13,15,17,27]. The parameters of the Weibull distribution function can be estimated using analytical methods, but numerical iterations are required. Additionally, when the mean and standard deviation of the wind speed are known, the MM can be used to calculate the shape factor and scaling factors without the requirement for numerical iterations.

Heuristic Optimization Methods
Other methods have been used for managing optimization problems, which use artificial intelligence algorithms. Most of these algorithms have been utilized to obtain the best objective function in a problem's search space by utilizing intelligence. However, in order to identify optimal values of shape and scale factors, many heuristic techniques rely on probabilistic judgments rather than analytical methods based on a pure random search [25,27]. The compared algorithms, including PSO and the proposed AO, are managed in this study. The main objective function represents minimizing the sum of square error between the measured frequency distribution in the histogram of the data and forecasted values by the adjusted PDF. Hence, the objective function can be represented by the following equation: where f w (v i ) and f m (v i ) are the wind speed probability from the Weibull distribution and that observed in the histogram of measured values, respectively.

Aquila Optimizer Model
The Aquila optimizer (AO) is a novel meta-heuristic optimization approach inspired by the natural behavior of the Aquila during the prey-capturing process. Aquila's speed and agility, along with robust feet and broad, pointed talons, allow it to capture a variety of prey, mostly rabbits, hares, deeps, marmots, squirrels, and other ground animals. As a result, there are many ways of hunting that rely on the natural behavior of prey. The choice of hunting tactics is also influenced by the hunting circumstances and Aquila's ability to switch back and forth between hunting strategies deftly and rapidly [33].
The following steps express the model of hunting processes of the AO algorithm: 1.
Initialization process: the solution of any optimization process starts with random solutions of candidate solution (X ij ) in the range between upper (up j ) and lower limits (lb j ) which is expressed as follows: where rand is a random number, ith is the number of populations, and jth is the dimension size of the problem.

2.
Firstly, the Aquila explores any prey in the search space. This exploration process is accomplished at high soar, which is called expanded exploration in the search space.
Once the Aquila explores prey, it drops with a vertical stoop to grab the prey. This behavior is mathematically as follows: where X 1 (t + 1) is the position of the Aquila in the next iteration of t, and X b (t) is the best-obtained solution until tth iteration. 1 − t T is used to control the exploration process in search space. X M (t) refers to the mean value of the solutions at the previous iteration, which is calculated using Equation (6).
where N is the population size.

3.
The second process is the most used by Aquila, and the hunting process is at contour flight with a short glide attack. Therefore, the Aquila is close to the hounded prey, which leads to narrow exploration in the search space. This process is modeled as: where X R (t) is a random solution taken at the ith iteration. Levy(D) is the levy flight distribution function for dimension space (D), which is calculated as follows: where s is a constant value equal to 0.01, and u and υ are random values between 0 and 1. ω is a dynamic adaptive coefficient that is calculated as follows: where β is a constant value equal to 1.5. The exploration process has a spiral shape, which is controlled by values of x and y that can be evaluated as follows: where r 1 is a value between 1 and 20, U is equal to 0.00565, D 1 is integer numbers from 1 to the dimension of search space, and W is equal to 0.005.

4.
The third process is the expanded exploitation (X 3 ) that belongs to any flooring prey that have a slow escape response. Aquila is flying at a low level, and once it selects its prey, it grabs progressively on the prey with a slow decent attack. This step is mathematically presented as in the following equation: where α and δ are the exploitation adjustment parameters which have range values from 0.1 to 0.9 and are fixed at 0.1 according to tests for several benchmarks. 5.
The fourth process (X 4 ) is when the Aquila walks on the land and captures prey by pulling it. It is used for dealing with large prey, which is described as a narrowed exploitation step. This step is mathematically modeled as: where QF is defined as a quality function used to equilibrium the search strategies, which is calculated using Equation (17). G 1 denotes various motions of the AO which is generated as Equation (18). G 2 is the flight slope of the AO, which has a value decreasing from 2 to 0 during the prey elope from first to last position and is described using Equation (29).
The steps of the AO optimizer are developed to solve the considered optimization problem for Weibull parameters estimation in the flowchart as in Figure 1. The main drawback of AO implies that it has several parameters and constants which affect the results and the convergence characteristics and make it difficult to improve.

Assessment of Fit Performance
As mentioned, there are several methods for estimating the Weibull parameters to calculate the PDF of the wind speed. In order to compare these methods and extract the most accurate one, several statistical measures are applied [8,15,22,27]. The tested methods are evaluated by the following measures:

Root Mean Square Error (RMSE)
The accuracy of the wind speed model relative to the actual distribution can be measured by RMSE, which is calculated as follows: This value is an indication of the deviation of the predicted from actual probability values. Therefore, the method which has a smaller value is more accurate.

Determination Coefficient (R 2 )
Another measure of performance is the determination coefficient (R 2 ), which measures the variance or the efficiency of the method. R 2 is defined by:

Mean Absolute Error (MAE)
The accuracy of the wind speed distribution model relative to the actual distribution

Assessment of Fit Performance
As mentioned, there are several methods for estimating the Weibull parameters to calculate the PDF of the wind speed. In order to compare these methods and extract the most accurate one, several statistical measures are applied [8,15,22,27]. The tested methods are evaluated by the following measures:

Root Mean Square Error (RMSE)
The accuracy of the wind speed model relative to the actual distribution can be measured by RMSE, which is calculated as follows: This value is an indication of the deviation of the predicted from actual probability values. Therefore, the method which has a smaller value is more accurate.

Determination Coefficient (R 2 )
Another measure of performance is the determination coefficient (R 2 ), which measures the variance or the efficiency of the method. R 2 is defined by: Wind 2022, 2 624

Mean Absolute Error (MAE)
The accuracy of the wind speed distribution model relative to the actual distribution can be determined by the mean absolute error, which is calculated as follows:

Wind Production Deviation (WPD)
The wind energy production deviation obtained between the predicted probability distribution curve and the actual date is evaluated as follows:

Wind Speed Data
The candidate sites in this study are the Zafarana and Shark El-Ouinat wind farm sites. The first site, the Zafarana wind farm, is located in the Suez Gulf, which is approximately 200 km southeast of Cairo. The second site, Shark El-Ouinat, is in the depth of the southern Egyptian desert. The operational data of wind turbines considered in this study for the two sites, are illustrated in Table 1. Added to that, the time series data of the wind speed for the year 2019 on an hourly basis of Zafarana and Shark El-Ouinat wind farm sites are taken from ref. [34]. The statistics for this data are presented in Table 2.

Goodness of Fit
In this section, various computations and analyses are introduced to compare the accuracy of the analytical and heuristic methods that are presented in previous sections. MLM, EM, MM, and EPFM as analytical methods and PSO and AO algorithms as heuristic methods are applied for optimal fitting in estimating the Weibull parameters. For the first site, the Zafarana wind farm, Table 3 describes their statistical comparisons. For the second site, the Shark El-Ouinate wind farm, Table 4 describes the statistical comparisons of the MLM, EM, MM, EPFM, PSO, and AO algorithms. From Table 4, the novel meta-heuristic AO algorithm provides the highest value of R 2 of 0.942592. On the other side, the other heuristic algorithm, PSO, finds an R 2 of 0.941294. Nevertheless, the analytical methods such as EPF, EMP, MLM, and MM find R 2 of 0.896893, 0.892958, 0.888387, and 0.892456, respectively. The AO algorithm provides the least value of RMSE of 0.010755. On the other side, the other heuristic algorithm, PSO, finds an RMSE of 0.010875. Nevertheless, analytical methods such as EPF, EMP, MLM and MM find RMSE of 0.014413, 0.014685, 0.014996, and 0.014719, respectively. Figures 2 and 3 present the obtained Weibull PDF based on the MLM, EM, MM, EPFM, PSO, and AO algorithms in comparison to the provided probability function f(v), for hourly time series in the year 2019 for the two studied sites. As shown, the proposed AO algorithm derives the greatest correctness with the most suitable adequacy to estimate the Weibull parameters. These findings conclude that the statistical tests for heuristic methods are better than those obtained by analytical methods. Moreover, it was further observed that the novel meta-heuristic AO algorithm was completely adequate for estimating the Weibull parameters with the minimum square of errors since it provides the highest value.

Statistical Results for PSO and AO for Estimating Weibull Parameters
The results of the same optimization algorithm are different for each run. To investigate the performance of the optimization algorithms of PSO and AO, they are run for 50 separate runs, and the statistical results of best, mean, worst, standard deviation (STD), and standard error (STE) are compared as tabulated in Tables 5 and 6 for the Zafarana and Shark El-Ouinate sites, respectively. For all algorithms, the same number of fitness evaluations is considered, 5000. For PSO and AO methods, the total number of iterations is 100, whereas the population size is taken as 50.     Table 5, the acquired statistical tests based on the AO algorithm are always the least indicator compared to PSO. The AO algorithm achieves the worst objective of 0.000717186, whereas PSO achieves the worst objective of 0.022794503. Added to that, the mean objective that is obtained by the AO algorithm is less than the best objectives that are attained by the other. Additionally, the AO algorithm achieves the least standard deviation and standard error of 1.7788 × 10 −5 and 2.51561 × 10 −6 , respectively. Similar findings are obtained from Table 6; the acquired statistical tests based on the AO algorithm are always the least indicator compared to PSO. The AO algorithm achieves the best, worst, and mean objectives with the same value of 0.001735. Furthermore, the AO algorithm achieves the least standard deviation, and standard error of 1.69 × 10 −5 and 2.4 × 10 −6 , respectively.
So, Figures 4 and 5 describe the histogram of the obtained objectives by means of PSO and AO algorithms for Zafarana and Shark El-Ouinate sites, respectively. As shown, the AO algorithm has greater stability and outperforms over PSO as it is always capable of finding the optimal value.  Moreover, Figures 6 and 7 describe the convergence curves of the PSO and AO algorithms for the Zafarana and Shark El-Ouinate sites, respectively. As shown, the error convergence sequence is compared to heuristic optimization methods at selected sites. It can be observed that AO has the fastest convergence compared to the other since it provides its best solution through less than 40% of the total number of iterations.  Moreover, Figures 6 and 7 describe the convergence curves of the PSO and AO algorithms for the Zafarana and Shark El-Ouinate sites, respectively. As shown, the error convergence sequence is compared to heuristic optimization methods at selected sites. It can be observed that AO has the fastest convergence compared to the other since it provides its best solution through less than 40% of the total number of iterations. Moreover, Figures 6 and 7 describe the convergence curves of the PSO and AO algorithms for the Zafarana and Shark El-Ouinate sites, respectively. As shown, the error convergence sequence is compared to heuristic optimization methods at selected sites. It can be observed that AO has the fastest convergence compared to the other since it provides its best solution through less than 40% of the total number of iterations.

Power Curve Modeling
There is no doubt that the determination of wind output power is a vital issue in making investment decisions and other issues related to the planning and operation of the power system. Therefore, the wind power forecasted error should be reduced as much as possible using efficient models. In addition to wind speed uncertainty, there is uncertainty in wind speed to power conversion. Although the most precise power curve model may

Power Curve Modeling
There is no doubt that the determination of wind output power is a vital issue in making investment decisions and other issues related to the planning and operation of the power system. Therefore, the wind power forecasted error should be reduced as much as possible using efficient models. In addition to wind speed uncertainty, there is uncertainty in wind speed to power conversion. Although the most precise power curve model may

Power Curve Modeling
There is no doubt that the determination of wind output power is a vital issue in making investment decisions and other issues related to the planning and operation of the power system. Therefore, the wind power forecasted error should be reduced as much as possible using efficient models. In addition to wind speed uncertainty, there is uncertainty in wind speed to power conversion. Although the most precise power curve model may be created, the projected inaccuracy from speed to power conversion may persist. Because of the existence of power curve uncertainty, each wind speed does not have a single value of output power, but several values of output power can be produced at the same wind speed [3]. The parameters of the forecasting model are not fixed for the same site because of the nonlinearity and randomness of the power curve [3,35]. Thus, traditional methods cannot represent the power uncertainty in a perfect manner. Figure 8 depicts the power curve that is produced in terms of wind speed. The error of energy density results for two sites is calculated where the distributions of energy density error for the nine models are shown in Figures 9 and 10 for the Zafarana and Shark El-Ouinate sites, respectively. From Figure 9, at the Zafarana site, the second model (g2) provides the minimum error of energy density of 5.728%. After that, the fourth (g4) and seventh (g7) models come with an error of energy density of 6.005 and 6.167% respectively. On the other side, all the models illustrate close errors where there is no great divergence in the results, as the worst model provides an error of energy density of 6.593%. From Figure 10, at the Shark El-Ouinate site, the eigth model (g8), the cubic form provides the minimum error of energy density of 13.2%. After that, the ninth (g4) and seventh (g7) models come with an error of energy density of 13.6 and 16.6%, respectively.   From Figure 8, the power curve is divided into four major parts. The first part occurs at wind speed less than cut-in speed, which leads to zero output power. In the second part, the wind speed is between cut-in (v in ) and rated speed (v r ), while the output power has nonlinear characteristics. The wind turbine output power has a constant value of rated output power for all models when the wind speed value is located between the rated and cut-out speed, as shown in the third part. To protect the wind turbine from damage, the wind turbine must be stopped when the wind speed becomes greater than the cut-out speed (v out ) [30]. In the literature, various models have been used to describe this nonlinear relationship, such as the third-degree polynomial curve in ref. [21], as well as fitted polynomial function as in refs. [35,36]. The general form can be represented by: where g(v) represents the relationship between the wind speed and the wind turbine output power that can be formed in several models as follows [29,30]: where a 1 , a 2 , a 3 , a 4 , and a 5 are polynominal coefficients calculated in MATLAB by the polyfit function. According to the previous comparison between analytical and heuristic methods for estimation of the Weibull distribution parameters, the proposed AO method provides the best fit for wind speed distribution for the two sites. Therefore, its output shape and scaling factors of the Weibull parameters are utilized, and then the detailed models (g 1 -g 9 ) are considered to find the best model for each site. In order to compare the performance of these nine models in estimating the energy from a wind turbine, the error (ε) between the actual energy density (E actual ) and estimated energy (E f itted ) is calculated by the following expression [21]: The error of energy density results for two sites is calculated where the distributions of energy density error for the nine models are shown in Figures 9 and 10 for the Zafarana and Shark El-Ouinate sites, respectively. From Figure 9, at the Zafarana site, the second model (g 2 ) provides the minimum error of energy density of 5.728%. After that, the fourth (g 4 ) and seventh (g 7 ) models come with an error of energy density of 6.005 and 6.167%, respectively. On the other side, all the models illustrate close errors where there is no great divergence in the results, as the worst model provides an error of energy density of 6.593%. From Figure 10, at the Shark El-Ouinate site, the eigth model (g 8 ), the cubic form, provides the minimum error of energy density of 13.2%. After that, the ninth (g 4 ) and seventh (g 7 ) models come with an error of energy density of 13.6 and 16.6%, respectively.
(g4) and seventh (g7) models come with an error of energy density of 6.005 and 6.167%, respectively. On the other side, all the models illustrate close errors where there is no great divergence in the results, as the worst model provides an error of energy density of 6.593%. From Figure 10, at the Shark El-Ouinate site, the eigth model (g8), the cubic form, provides the minimum error of energy density of 13.2%. After that, the ninth (g4) and seventh (g7) models come with an error of energy density of 13.6 and 16.6%, respectively.   energy density error % models Figure 9. Energy density error of the nine models for Zafarana site.
Wind 2022, 2, FOR PEER REVIEW 16 Figure 10. Energy density error of the nine models for Shark El-Ouinate site.

Conclusions
This paper investigates different analytical and heuristic techniques for optimal parameter estimation of the Weibull distribution parameters to analyze the wind speed profile for the Zafarana and Shark El-Ouinate sites, Egypt. In this study, four analytical methods of MLM, MM, EM, and EPFM, and other heuristic methods of PSO and AO algorithms are introduced. Several fitness metrics are evaluated for comparative assessment between the competitive methods. The novel meta-heuristic AO algorithm shows the best performance in acquiring the highest value of R 2 , the least value of RMSE, and the least value of MAE for tuning Weibull PDF parameters. Additionally, the mean objective obtained by the AO algorithm is better than the others. Further, the AO algorithm is the most robust technique as it achieves the least standard deviation and standard error compared to the others. Furthermore, the AO algorithm demonstrates the greatest stability characteristics as its histogram of the obtained objectives is better than the PSO algorithm. Moreover, the AO algorithm presents the fastest convergence compared to the other since it provides its best solution through the first and second site parameters calculation, respectively. Furthermore, different models are analyzed to deduce the nonlinear relationship between the wind output power and the speed, where the error of wind energy density between actual and estimated is greatly minimized using the cubic model. A suggestion for future work is to obtain a rank for each method, where the mean wind speed and its standard deviation are estimated for a greater number of periods. energy density error % models Figure 10. Energy density error of the nine models for Shark El-Ouinate site.

Conclusions
This paper investigates different analytical and heuristic techniques for optimal parameter estimation of the Weibull distribution parameters to analyze the wind speed profile for the Zafarana and Shark El-Ouinate sites, Egypt. In this study, four analytical methods of MLM, MM, EM, and EPFM, and other heuristic methods of PSO and AO algorithms are introduced. Several fitness metrics are evaluated for comparative assessment between the competitive methods. The novel meta-heuristic AO algorithm shows the best performance in acquiring the highest value of R 2 , the least value of RMSE, and the least value of MAE for tuning Weibull PDF parameters. Additionally, the mean objective obtained by the AO algorithm is better than the others. Further, the AO algorithm is the most robust technique as it achieves the least standard deviation and standard error compared to the others. Furthermore, the AO algorithm demonstrates the greatest stability characteristics as its histogram of the obtained objectives is better than the PSO algorithm. Moreover, the AO algorithm presents the fastest convergence compared to the other since it provides its best solution through the first and second site parameters calculation, respectively. Furthermore, different models are analyzed to deduce the nonlinear relationship between the wind output power and the speed, where the error of wind energy density between actual and estimated is greatly minimized using the cubic model. A suggestion for future work is to obtain a rank for each method, where the mean wind speed and its standard deviation are estimated for a greater number of periods.