Wind Power Integration: An Experimental Investigation for Powering Local Communities

The incorporation of wind energy as a non-conventional energy source has received a lot of attention. The selection of wind turbine (WT) prototypes and their installation based on assessment and analysis is considered as a major problem. This paper focuses on addressing the aforementioned issues through a Weibull distribution technique based on five different methods. The accurate results are obtained by considering the real-time data of a particular site located in the coastal zone of Pakistan. Based on the computations, it is observed that the proposed site has most suitable wind characteristics, low turbulence intensity, wind shear exponent located in a safe region, adequate generation with the most adequate capacity factor and wind potential. The wind potential of the proposed site is explicitly evaluated with the support of wind rose diagrams at different heights. The energy generated by ten different prototypes will suggest the most optimum and implausible WT models. Correspondingly, the most capricious as well as optimal methods are also classified among the five Weibull parameters. Moreover, this study provides a meaningful course of action for the selection of a suitable site, WT prototype and parameters evaluation based on the real-time data for powering local communities.


Introduction
Owing to the escalation in energy demand to meet fundamental life necessities, researchers around the world are trying to find the best ways to employ renewable energy resources. Several techniques are proposed in the literature, different methods are introduced and numerous policies are presented to get an optimum solution. The installation of wind energy resources along the coastline is getting much attention due to the opportunity of integration of both on-shore and off-shore wind farms. In this work Pakistan is considered as an example, due to its 1600 km long coastline and the energy crisis Pakistan is facing, which has badly affected the economy and daily life of its 220 million citizens. The Government of Pakistan has tried many initiatives to reduce load shedding and provide clean, affordable and continuous electric supply. As a result there is a need to explore more available options instead of relying only on traditional fossil fuels. The Alternative Energy Distribution Board (AEDB) has taken serious steps to promote renewable energy by exploring the available wind, solar and hydro energy resource potential. The National Electric Power Regulation Authority (NEPRA) has encouraged the generation of renewable energy by giving letters of intent (LOIs) to many contractors. Figure 1 shows the share of power generation as a percentage of different resources [1]. The country is moving towards sustainable renewable energy resources as compared to past five years, but currently, a significant percentage share in country's electrical power generation is still based on fossil fuels. Energy Distribution Board (AEDB) has taken serious steps to promote renewable energy by exploring the available wind, solar and hydro energy resource potential. The National Electric Power Regulation Authority (NEPRA) has encouraged the generation of renewable energy by giving letters of intent (LOIs) to many contractors. Figure 1 shows the share of power generation as a percentage of different resources [1]. The country is moving towards sustainable renewable energy resources as compared to past five years, but currently, a significant percentage share in country's electrical power generation is still based on fossil fuels. The wind energy potential in Pakistan is the major focus for power generation as renewable energy is a most fundamental necessity to help mankind and meet energy demands, so an analysis of long term pathways for power generation is described in [2]. It is also depicted that compromising on the issues of renewable energy just for the sake of convenience is not a practical approach. Likewise, selection criteria for the most appropriate sites in Pakistan for wind energy generation based on an analytic hierarchal process is described in [3]. An analysis on energy security according to renewable the energy policy of Pakistan is elaborated in [4]. A corresponding research analysis considering the Badin and Pasni territories of Pakistan is presented in [5]. The contemporary development of wind energy, challenges and future recommendations are deliberated in [6]. The wind energy potential and power law assessment in Malaysia, economic value assessment analyzed for an optimal sizing of an energy storage system and integration of digital control techniques based on power electronics converters are described in [7][8][9].
Domestic use, as well as imports of these conventional sources, has put a massive burden on the economy of the country and also caused environmental and health issues. Pakistan has a tremendous wind energy potential available in Sindh and Baluchistan. There is a terrible need for the country to move towards renewable energies like hydro, wind, solar and bioenergy. The Ministry of Planning, Development, and Reforms in Pakistan has taken the initiative under the 11 th five-year plan, in which energy crisis is cogitated, and strategies to meet the energy demands of Pakistan have been proposed [1]. Figure 2 shows the demand and supply projections for fuel from 2013 to 2018. Oil (including LPG), gas, coal, total primary energy supply (TPES) and imports as the  The wind energy potential in Pakistan is the major focus for power generation as renewable energy is a most fundamental necessity to help mankind and meet energy demands, so an analysis of long term pathways for power generation is described in [2]. It is also depicted that compromising on the issues of renewable energy just for the sake of convenience is not a practical approach. Likewise, selection criteria for the most appropriate sites in Pakistan for wind energy generation based on an analytic hierarchal process is described in [3]. An analysis on energy security according to renewable the energy policy of Pakistan is elaborated in [4]. A corresponding research analysis considering the Badin and Pasni territories of Pakistan is presented in [5]. The contemporary development of wind energy, challenges and future recommendations are deliberated in [6]. The wind energy potential and power law assessment in Malaysia, economic value assessment analyzed for an optimal sizing of an energy storage system and integration of digital control techniques based on power electronics converters are described in [7][8][9].
Domestic use, as well as imports of these conventional sources, has put a massive burden on the economy of the country and also caused environmental and health issues. Pakistan has a tremendous wind energy potential available in Sindh and Baluchistan. There is a terrible need for the country to move towards renewable energies like hydro, wind, solar and bioenergy. The Ministry of Planning, Development, and Reforms in Pakistan has taken the initiative under the 11th five-year plan, in which energy crisis is cogitated, and strategies to meet the energy demands of Pakistan have been proposed [1]. Figure 2 shows the demand and supply projections for fuel from 2013 to 2018. Oil (including LPG), gas, coal, total primary energy supply (TPES) and imports as the percentage of TPES resulted in a million tons of oil equivalent (MTOE). In 2017-2018 total domestic fuel produced by Pakistan in MTOE was about 48.43 and fuel imports of 45.87 MTOE. As compared to the previous year 2016-2017, the domestic fuel production increased only 10.45%. Figure 3 shows the per annum demand for power (MW) by NTDC and KESC. During the period of 2017-2018, power demand by National Transmission and Dispatch Company (NTDC) and Karachi Electric Supply Company (KESC) was 26.53 GW and 4.5 GW, respectively.   Total power demand based on the performance of models is tested and evaluated to suggest more accurate Weibull distribution model. One-year data is analyzed for this purpose. Wind characteristics are estimated on the basis of daily as well as monthly average speeds, the WR, air density, WPD, energy in terms of kW/m 2 , shear exponent coefficient (∝) for boundary layer of the site, and turbulence intensity of the proposed site. Ten different wind turbine models are used to estimate wind power output and energy generated at the proposed site. More capacity factor of each wind turbine is calculated to suggest a more efficient model for this site.
MATLAB is used for the assessment of wind potential at the proposed site. Data is acquired with the support of Pakistan Metrological department (PMD). However, for more accuracy, data quality assurance test is conducted to avoid errors in results, by generating a code which organizes and inspects wind data for errors. The excellent observation data is used to estimate stated parameters. Results suggested that this site has a vast potential for wind energy that can be exploited to meet a moderate portion of energy demand in Pakistan. The proposed site is in premises of 600 kV HVDC line from Matiari to Lahore which is planned up to 2021-2022 under CPEC [10].
In this paper, five different methods based on Weibull distribution techniques have been cogitated. The real time data at various heights from a site located in coastal zone of Pakistan has been employed to acquire different parameters. Based on the calculations and analysis, most optimal wind turbine having ability to be operated on the highest wind potential and most adequate capacity factor is sorted out. The wind potential at the different heights and in proper direction on proposed site is explicitly evaluated through wind rose diagrams. The energy production capacity of ten different prototypes is calculated and the most optimum and implausible WT models are suggested along with their pros. and cons. according to the different wind classes. Congruently, the most improvident as well as finest methods are also rated among five Weibull parameters. Additionally, this paper also provides a significant strategy for the determining an appropriate site, wind turbine prototype and parameters evaluation based on the real-time data with an instance from coastal area of Pakistan for providing electrical energy in the local communities. It is expected that this study will be an imperative contribution in understanding and applying the numerous Weibull distribution techniques and also in selection of wind turbine models by considering various parameters. Therefore, the results will be undoubtedly helpful in making an efficient energy policy in the future.
The rest of the paper is organized as follows: Wind data assessment and analysis are thoroughly described in Section 2, site characteristics are appropriately addressed in Section 3, the analysis, results, and discussion are conducted in Section 4. Whereas, conclusions are drawn in Section 5.

Wind Data Assessment
The characteristics of wind have time a varying nature at each instant. There are various parameters that are key indicators to estimate that whether a particular site is suitable for a utility-scale wind power project [11][12][13].
These parameters are air density, ambient temperature, turbine class and the hub height of turbine at which blades usually capture the wind energy. For accuracy, it is essential to assess wind energy potential of proposed site critically by considering aforementioned wind characteristics. Wind power density is considered as another critical factor in determining how much wind power (W) is available at per unit area (m 2 ) area. This can be achieved by probability distribution function [14]. The wind map of Pakistan is provided in Figure 4, where several wind classes are classified by taking wind power generation into account [13].  These parameters are air density, ambient temperature, turbine class and the hub height of turbine at which blades usually capture the wind energy. For accuracy, it is essential to assess wind energy potential of proposed site critically by considering aforementioned wind characteristics. Wind power density is considered as another critical factor in determining how much wind power (W) is available at per unit area (m 2 ) area. This can be achieved by probability distribution function [14]. The wind map of Pakistan is provided in Figure 4, where several wind classes are classified by taking wind power generation into account [13].

Wind Characteristics
To assess the wind potential of a particular site, it should be considered that wind energy is a time-varying entity which not only changes the magnitude but also the direction [15][16]. Moreover, above characteristics are significant for adequate evaluation of wind potential available in particular site. Wind turbine performance is affected by the wind rose and frequency of wind distribution during the higher occurrence of wind speed [17][18][19][20]. An appropriate method for analyzing the wind characteristics and wind energy potential is referred in [21]. Moreover, if the data for more than two years have been considered, then measure-correlate-predict (MCP) method is clearly described and implemented in [22]. Whereas, the presence of probable errors and accuracy of wind speed data is also discussed in the same article.

Wind Characteristics
To assess the wind potential of a particular site, it should be considered that wind energy is a time-varying entity which not only changes the magnitude but also the direction [15,16]. Moreover, above characteristics are significant for adequate evaluation of wind potential available in particular site. Wind turbine performance is affected by the wind rose and frequency of wind distribution during the higher occurrence of wind speed [17][18][19][20]. An appropriate method for analyzing the wind characteristics and wind energy potential is referred in [21]. Moreover, if the data for more than two years have been considered, then measure-correlate-predict (MCP) method is clearly described and implemented in [22]. Whereas, the presence of probable errors and accuracy of wind speed data is also discussed in the same article.

Average Wind Speed, Variance and Standard Deviation
The mean or average wind speed 'v avg ' is obtained from Equation (1). Whereas, the variance 'σ 2 ' and standard deviation (SD) 'σ', for wind speed data are calculated from Equations (2) and (3) respectively: where n is a number of observations, v i is the i-th wind speed, i is the i-th observation. The air density (ρ) for the proposed site can be obtained by using Equation (4) [23][24][25], wind power density (WPD) and energy (E) from Equations (5) and (6): where, P is air pressure (Pa or N/m 2 ), R is the specific gas constant (287 J/kg), and T is air temperature at the site in Kelvin (C + 273 • ): where ρ is the air density, A T is the swept area of turbine blades (m 2 ), P w is wind power (W), v is wind velocity (m/s) and C p is the Betz limit equal to 0.593 or the maximum value of C p , the performance for the ideal WT. Furthermore, due to mechanical deficiency of a real turbine, the fraction of the power extracted from the wind will be less than that for an ideal WT. In other words, this limit states how effectively a WT converts the wind energy into electricity [26][27][28]: where, E is energy output in terms of Weibull distribution at the proposed site in (kWh/m 2 ). P(V), T, V, k and c are wind turbine's power curve, time period, wind velocity, shape and scale parameters, respectively [29].

Wind Turbulence Intensity, Shear, and Power Law
Turbulence Intensity (TI) is defined as the ten-minute standard deviation of the velocity divided by the ten-minute mean velocity of the wind as given by Equation (7) [25]: where σ is considered a ten-minute (SD) of wind velocity, and V is the ten-minute average velocity of the proposed site. The exponent of wind shear and shear of the site is calculated from Equations (8) and (9) respectively [30][31][32]. The term α is essential to estimate wind velocities at higher altitude by processing the wind velocities measured at lower or previous altitudes. The power law is used to calculate wind speed at hub height by using Equation (10) [32]: where z 1 and z 2 are heights, v 1 and v 2 are wind speeds, α is the coefficient of wind shear and h is the reference height.

Wind Power Classes
Elliot and Schwartz classified wind power into seven classes, considering the wind speed and power density of a particular site [33][34][35][36][37]. The wind power class 1 and 2 are for rural applications, and class 4 and beyond are for commercial purposes [38][39][40][41]. They defined these classes at heights  [42]. Wind speed above 5.5 m/s yields power generation that is economical and located in class 3. While classes 1 and 2 are for micro-generation purpose [43]. Several parameters regarding wind potential in Pakistan, based on the utility scale of wind class, are articulated in Table 1, which depict that 26,362 km 2 of land can be used to produce 131,820 MW of electricity through wind turbines [12]. Whereas, the international standards of wind power generation classification at various heights are categorized in Table 2 [24].

Weibull Distribution
Weibull distribution of the wind speed data collected from a particular location is a probability density function f (v) as well as the cumulative distribution function F(v), it can be calculated from Equations (11) and (12), respectively [44][45][46][47]: v avg can be calculated as follows: whereas σ can be obtained from Equation (14), as given below: where Γ is a gamma function of (y), as given by:

Different Weibull Methods
The shape (k) and scale (c) parameters for Weibull distribution function can be calculated by numerous methods suggested in the literature. In this study five methods namely, GM, MMLM, EPF, EMJ and EML are used in this study to calculate Weibull parameters. In the MMLM method, it is necessary that wind speed data should be in FD format and a number of iterations should be performed to find the shape (k) and scale (c) parameters for Weibull distribution function, which can be obtained from Equations (16) and (17) [48,49]: In the EML method, Lysen suggested that the shape (k) and scale (c) parameters are obtained from Equations (18) and (19), respectively [50]: The EMJ, also known as an empirical method of Jestus in which, the shape (k) can be calculated by Equation (18) and scale (c) parameter is calculated by Equation (20) [51,52]: In the Graphical Method (GM), the cumulative distribution function of Weibull distribution is used, in which wind data is sorted into bins due to the least squares regression. The graphical method's equation can be obtained by taking double logarithms of Equation (12) [53,54]: Comparing Equation (20) with y = ax + b, we get: where, y and x can be calculated by using measured wind speed data. The standard least regression method is applied to obtain the slope (a) and the intercept (b), and then the shape (k) and the scale (c) parameters can be determined as: Energies 2019, 12, 621 9 of 24 From EPF, the average wind speed data is considered to calculate shape (k) and scale (c) parameters of the Weibull distribution. First, the parameter for aerodynamic design of the turbine should be defined as [55]: From Equation (24), the energy pattern factor E p f can be calculated as a ratio of the average of the cubic value of wind speed data over the cubic value of average wind speed data. The shape parameter (k) obtained as [55]: The scale parameter (c) can be calculated same as Equation (20).

Goodness of Fit Test
The performance of five Weibull distribution methods is compared by using five statistical analyses, including mean squared error (MSE) given by Equation (26), root mean squared error (RMSE) from Equation (27), mean absolute error (MAE) by using Equation (28), the coefficient of correlation (R) by using Equation (29) and coefficient of determination (R 2 ) from Equation (30) as follows: where y i represents the i-th actual wind speed, x i is the i-th predicted wind speed, y is mean of the actual wind speed and N is the number of observations.

Site Characteristics
The Sanghar site wind mast funded by the World Bank is located in Sanghar, Sindh, Pakistan. The height of the mast is 80 m, and the geographic location of the site is 25 • 48 57.26 N and 69 • 2 15.12 E. Figure 5 shows wind mast installed at the Sanghar site with the terrain and surroundings. The proposed site is flat, wide and opens having no obstruction with an elevation higher than 20 m. Figure 6 shows topographic maps for the Sanghar site. The ruggedness index (RIX) at a specific location is the percentage of the ground surface that has a slope above a given threshold (here 30%) within a certain distance (here 3.5 km). The site is located 250 km (4-h drive) away from Jinnah International Airport Karachi, and easily accessible by any type of vehicles. One year 10 min' average data is analyzed from November 2016 to October 2017.
Energies 2019, , x FOR PEER REVIEW 11 of 26 drive) away from Jinnah International Airport Karachi, and easily accessible by any type of vehicles.
One year 10 minutes' average data is analyzed from November 2016 to October 2017.

Wind Speed Measurement
In Figure

Air Density and Turbulence Intensity Measurement
Air density () is calculated by using temperature, and pressure (Pa) measured at the heights of 76 m and 4 m, respectively. It is a critical factor as it is directly proportional to the wind power density (WPD). 10-minutes average values for temperature and pressure were used in the estimation of air density to achieve accurate WPD.
In Figure 8, the relation between monthly average of air density and temperature is indicated.

Air Density and Turbulence Intensity Measurement
Air density (ρ) is calculated by using temperature, and pressure (Pa) measured at the heights of 76 m and 4 m, respectively. It is a critical factor as it is directly proportional to the wind power density (WPD). 10-min average values for temperature and pressure were used in the estimation of air density to achieve accurate WPD.
In Figure 8, the relation between monthly average of air density and temperature is indicated. length for specific areas is 0.5 m and 0 m for towns and rivers, respectively. Therefore, the site has minimum turbulence and shear. The Turbulence Intensities (TI) at 80 m, 60 m, and 40 m heights are analyzed using 10 min interval yearly data of standard deviation and average wind speed. TI at different heights are shown in Figure 9, which clearly depicts that as altitude rises above ground level from 40 m and above, turbulence intensities are reduced. Therefore, relatively at 80 m or below heights, the wind velocities are more uniform. For the design purpose of a wind turbine, International Electro technical Commission (IEC) has formed standards for TI up to 18%, for a wind speed of 15 m/s, and this site has minimum TI, as perceived in Figure 9. The overall air density (ρ) observed is 1.167 kg/m 3 , and overall temperature is 27.27 • C. During the summer season, air density is lower than that of the winter season with negligible variation. During the cool-dry winter, hot-dry spring, rainy summer and autumn seasons, the air densities are 1.205 kg/m 3 , 1.155 kg/m 3 , 1.148 kg/m 3 and 1.1675 kg/m 3 , respectively.
It can be observed from the topographical map shown in Figure 6 that the roughness length is 0.07 m, corresponding to an open field with distributed rows of trees and low buildings. Roughness length for specific areas is 0.5 m and 0 m for towns and rivers, respectively. Therefore, the site has minimum turbulence and shear. The Turbulence Intensities (TI) at 80 m, 60 m, and 40 m heights are analyzed using 10 min interval yearly data of standard deviation and average wind speed. TI at different heights are shown in Figure 9, which clearly depicts that as altitude rises above ground level from 40 m and above, turbulence intensities are reduced. Therefore, relatively at 80 m or below heights, the wind velocities are more uniform. For the design purpose of a wind turbine, International Electro technical Commission (IEC) has formed standards for TI up to 18%, for a wind speed of 15 m/s, and this site has minimum TI, as perceived in Figure 9. The overall air density () observed is 1.167 kg/m 3 , and overall temperature is 27.27 °C . During the summer season, air density is lower than that of the winter season with negligible variation. During the cool-dry winter, hot-dry spring, rainy summer and autumn seasons, the air densities are 1.205 kg/m 3 , 1.155 kg/m 3 , 1.148 kg/m 3 and 1.1675 kg/m 3 , respectively.
It can be observed from the topographical map shown in Figure 6 that the roughness length is 0.07 m, corresponding to an open field with distributed rows of trees and low buildings. Roughness length for specific areas is 0.5 m and 0 m for towns and rivers, respectively. Therefore, the site has minimum turbulence and shear. The Turbulence Intensities (TI) at 80 m, 60 m, and 40 m heights are analyzed using 10 min interval yearly data of standard deviation and average wind speed. TI at different heights are shown in Figure 9, which clearly depicts that as altitude rises above ground level from 40 m and above, turbulence intensities are reduced. Therefore, relatively at 80 m or below heights, the wind velocities are more uniform. For the design purpose of a wind turbine, International Electro technical Commission (IEC) has formed standards for TI up to 18%, for a wind speed of 15 m/s, and this site has minimum TI, as perceived in Figure 9. The use of standard (constant) air density i.e., 1.255 kg/m 3 for the analysis can cause the variation of wind power density results [27]. When the wind power density is calculated using the air pressure and temperature data of candidature site, the results are slightly different. Therefore, it is better to use air density data of the selected site instead of a standard value for air density [56,57]. Whereas in this research work, authors have utilized the 10 min average actual data of temperature and air pressure of selected site to evaluate the 10 min average air density. Moreover, the calculated 10 min average air density values were utilized to calculate the 10 min average wind power density. Finally, the monthly and yearly average wind power densities were calculated to ensure the accuracy in results. Hence, Figure 8 shows the monthly average air density and temperature.

Wind Rose and Wind Shear Measurement
The wind rose obtained by using 10 min. average wind direction, and wind speed data at different heights is shown in Figure 10 . This is ideal for a wind turbine to capture more wind to work effectively and produce more power, instead of changing direction. This also means at hub height, the wind turbine can harness more available wind power at this site. In Figure 10  The monthly mean wind shear exponent () and shear profile of this site are shown in Figure  11a

Wind Speed Distribution and Methods
Moreover, the lowest values were observed in April, May, June, July, August, and September. In the shear profile, the wind velocities of all sensors at heights of 80 m, 60 m, and 40 m were fitted to estimate hub height velocity at 100 m. This is estimated to be 6.75 m/s as shown in Figure 12. The

Wind Speed Distribution and Methods
Moreover, the lowest values were observed in April, May, June, July, August, and September. In the shear profile, the wind velocities of all sensors at heights of 80 m, 60 m, and 40 m were fitted to estimate hub height velocity at 100 m. This is estimated to be 6.75 m/s as shown in Figure 12. The comparison among five Weibull methods fitting over wind speed data at each height including V hub is shown in Figure 13. For wind speed distribution at heights of 80 m and 60 m, it can be observed that all methods fitted appropriately; however, the Graphical Method (GM) showed poor performance over the measured data. At the 40 m height, GM and MMLM almost fit the data, but EMJ, EML, and EPF have slightly differences in fitting, so MMLM is considered as the finest method for estimating wind power potential at this site. Moreover, details on the monthly and overall values of mean wind speed, shape (k), scale (c), wind power density and energy obtained by using MMLM method is shown in Table 3, where the estimated monthly mean and yearly values of wind speed, shape (k), scale (c), wind power density (W/m 2 ) and energy density (kWh/m 2 ) at different heights are shown. These parameters are obtained using the MMLE method for a Weibull distribution. A one-year period of data is analyzed. At hub height (100 m), the maximum mean wind speed is 9.35 m/s in June, shape parameter (k)

Wind Speed Distribution and Methods
Moreover, the lowest values were observed in April, May, June, July, August, and September. In the shear profile, the wind velocities of all sensors at heights of 80 m, 60 m, and 40 m were fitted to estimate hub height velocity at 100 m. This is estimated to be 6.75 m/s as shown in Figure 12. The comparison among five Weibull methods fitting over wind speed data at each height including Vhub is shown in Figure 13. For wind speed distribution at heights of 80 m and 60 m, it can be observed that all methods fitted appropriately; however, the Graphical Method (GM) showed poor performance over the measured data. At the 40 m height, GM and MMLM almost fit the data, but EMJ, EML, and EPF have slightly differences in fitting, so MMLM is considered as the finest method for estimating wind power potential at this site. Moreover, details on the monthly and overall values of mean wind speed, shape () k , scale () c , wind power density and energy obtained by using MMLM method is shown in Table 3, where the estimated monthly mean and yearly values of wind speed, shape () k , scale () c , wind power density (W/m 2 ) and energy density (kWh/m 2 ) at different heights are shown. These parameters are obtained using the MMLE method for a Weibull distribution. A one-year period of data is analyzed. At hub height (100 m), the maximum mean wind speed is 9.     The seasonal average values of the shape (k) and scale (c) parameters calculated using MMLM, EMJ, EML, EPF and GM methods at 100 m hub height for the Sanghar site, are shown in Figure 14. The shape (k) parameter calculated using MMLM, EMJ, EML, and EPF shows a small variation for all seasons, whereas the shape (k) parameter estimated by the GM method shows lower values in winter, spring and summer. The scale (c) parameter for all five methods shows closeness for all seasons with the slight increase for GM in spring and summer.   The seasonal average values of the shape (k) and scale (c) parameters calculated using MMLM, EMJ, EML, EPF and GM methods at 100 m hub height for the Sanghar site, are shown in Figure 14. The shape (k) parameter calculated using MMLM, EMJ, EML, and EPF shows a small variation for all seasons, whereas the shape (k) parameter estimated by the GM method shows lower values in winter, spring and summer. The scale (c) parameter for all five methods shows closeness for all seasons with the slight increase for GM in spring and summer.
In this research work, the authors followed the international standard, IEC 61400-12-1 in order to select the WTs to estimate the power performance of a wind turbine. Through the IEC standards, the power curve of the wind turbine is considered for calculating the annual energy production [58]. The performance of a wind turbine is indicated by the power curve. There are various methods available in the literature to develop the power curve model. In which the turbine specifications from the manufacturers and the wind speed data of a particular site are used [59,60]. Moreover, in this study, the discrete model is used. The estimation of the output power of a wind turbine is done by arranging the wind speed data of the studied site into bins and the reliable hours for every bin are collected. According to the IEC international standards, the discrete model is a simple method that does not require any mathematical functions for describing the curve [58]. The power curve of each turbine was modeled into tabular form and then the number of probable hours for each bin and corresponding power output were calculated.
The average energy generated in GWh/year by ten different turbines at four hub heights is estimated and shown in Figure 15a. The detail of wind turbine models with their characteristics is provided in Table 4 The capacity factor (%) of ten different turbines at four hub heights on the proposed site is calculated and revealed in Figure 15b. The capacity factor for all ten turbines is located in the appropriate range with maximum values of 63.5%, 59.9%, 54.8% and 49.5% for Goldwind GW121/2500 at four hub heights. Followed by Vestas V126/3300 and Gamesa G97/2.0, with a slight difference between both, estimated C.F is 58.3, 54.4, 48.9, 43.4 and 58.1, 54.1, 48.7, 43.3%, respectively. The lowest C.F 37.5, 33.4, 28.6 and 24.3% was estimated for the Nordex n80/2500 wind turbine at four hub heights, which is an economical range of C.F. For all turbine models, C.F. lies in the excellent range, which shows that this site has a tremendous capacity to generate electrical power for most of the available wind turbines. More details are specified in Table 5, whereas, Table 6 shows the annual values of shape (k) and scale (c) parameters obtained by EMJ, EML, EPF GM and MMLM methods at the heights of 80 m, 60 m, and 40 m. The overall results of performance evaluation of different Weibull models based on the parameters explained in Section 2.5 are compared in Table 7. In this research work, the authors followed the international standard, IEC 61400-12-1 in order to select the WTs to estimate the power performance of a wind turbine. Through the IEC standards, the power curve of the wind turbine is considered for calculating the annual energy production [58]. The performance of a wind turbine is indicated by the power curve. There are various methods available in the literature to develop the power curve model. In which the turbine specifications from the manufacturers and the wind speed data of a particular site are used [59][60]. Moreover, in this study, the discrete model is used. The estimation of the output power of a wind turbine is done by arranging the wind speed data of the studied site into bins and the reliable hours for every bin are collected. According to the IEC international standards, the discrete model is a simple method that does not require any mathematical functions for describing the curve [58]. The power curve of each turbine was modeled into tabular form and then the number of probable hours for each bin and corresponding power output were calculated.
The average energy generated in GWh/year by ten different turbines at four hub heights is estimated and shown in Figure 15(a). The detail of wind turbine models with their characteristics is provided in Table 4. The performance of Vestas V126/3300 and Goldwind GW121/2500 outperformed the rest of turbine models at four hub heights (120 m, 100 m, 80 m, and 60 m). The In additions, Figure 15, Tables 4 and 5 show the units generated (GWh/year) and CF by particular wind turbines. However, the CF was estimated by the ratio of the average power output to the rated power of the generator [61]. It is used to estimate the average energy production (units generated) of a wind turbine to estimate the cost of the system [61][62][63].  Table 5. Power generated, the energy produced and capacity factor of ten different turbine models.

Conclusions
In this study, different wind turbine prototype models are considered based on the wind characteristics for evaluating and investigating the wind power potential along the coastline. For this purpose, a site located in coastal zone of Pakistan is evaluated owing to its vast potential for installation of commercial wind farms. Real-time one year (Nov. 2016 to Oct. 2017) wind speed data was analyzed to estimate the twelve months mean values as well as the yearly mean at the proposed site where a wind mast had been installed and funded by the World Bank. The data was observed on the daily basis (24 h) in order to predict the behavior of wind more accurately. Based on the results, the following conclusions are obtained: i.
Monthly and diurnal wind characteristics of the proposed site indicate that it has class three or higher wind potential. Therefore, according to the international standards, the proposed site is most suitable as depicted in Table 2 for wind turbine installation to supply the electricity to local communities. ii.
Basically, in the whole year, it is observed from the wind rose diagram, the wind directions at different heights show that most of the wind gusts blow from the west. iii.
The proposed site is most suitable with advantages on air density, wind shear exponent, low turbulence intensity, adequate wind speed distribution, and reliable capacity factor. iv.
In this research, the energy generated by ten different suggested wind turbine models indicates that the performance of two WTs such as; Vestas V126/3300 and Goldwind GW 121/2500 is better as compared to other prototypes for this proposed site. v.
Also, from five Weibull distribution techniques, GM shows the inadequate performance, whereas, MMLM presents the optimum performance as compared to its counter techniques. vi.
The maximum values of Weibull shape and scale parameters were also calculated to understand the seasonal wind characteristics for this site. During the spring-summer seasons, the shape values are 2.60 (spring) and 3.25 (summer) respectively. Whereas, the scale values are more than 9.0 (summer) and approximately 8.2 (spring), and even better for whole of the years.
It can be noted that this site has comparatively better wind characteristics. Additionally, this study also aims to provide a policy implication for optimal selection of a suitable site, wind turbine prototype and parameters evaluation based on the real-time data assessment and analysis. It is recommended that GoP should focus on the proposed site for future economic development and electricity generation for powering local communities and underprivileged regions.

Conflicts of Interest:
The authors declare no conflict of interest.