# Assessment of Wind Characteristics and Wind Power Potential of Gharo, Pakistan

Authors: Zahid Hussain Hulio
Journal: Journal of Renewable Energy
 (2021)

---

## Abstract

The objective of this research work is to assess the wind characteristics and wind power potential of Gharo site. The wind parameters of the site have been used to calculate the wind power density, annual energy yield, and capacity factors at 10, 30, and 50 m. The wind frequency distribution including seasonal as well as percentage of seasonal frequency distribution has been investigated to determine accurately the wind power of the site. The coefficient of variation is calculated at three different heights. Also, economic assessment per kWh of energy has been carried out. The site-specific annual mean wind speeds were 6.89, 5.85, and 3.85 m/s at 50, 30, and 10 m heights with corresponding standard deviations of 2.946, 2.489, and 2.040. The mean values of the Weibullk parameter are estimated as 2.946, 2.489, and 2.040 while those of scale parameter are estimated as 7.634, 6.465, and 4.180 m/s at 50, 30, and 10 m, respectively. The respective mean wind power and energy density values are found to be 118.3, 92.20, and 46.10 W/m2 and 1036.6, 807.90, and 402.60 kWh/m2. As per cost estimation of wind turbines, the wind turbine WT-C has the lowest cost of US$ Cents 0.0346/kWh and highest capacity factors of 0.3278 (32.78%). Wind turbine WT-C is recommended for this site for the wind farm deployment due to high energy generation and minimum price of energy. The results show the appropriateness of the methodology for assessing the wind speed and economic assessment at the lowest price of energy.


---

## Body

## 1. Introduction

The United Nations General Assembly (UNGA) in the year 2015 provided a powerful mechanism for an international corporation on sustainable development goals (SGD’S) of Earth. The agenda 2030 comprises three important facets including decreasing poverty, right full distribution of resources and justice, and sustainable energy and environment of the planet. On the face of it, the green energy was central to an agenda that focuses on reliable, affordable, and accessible energy to all. The overall objective of UNGA to meet the future energy demands via renewable energy is achievable while minimizing the impact of dangerous environment.In recent years, the rising trend of electricity generation by renewable energy is taking a positive way all over the world. Owing to growth of renewable energy resource, the wind energy generation is ahead of all natural resources. There are some noteworthy reasons that played constructive role in development of wind energy generation including atmosphere, state-of-the-art design of wind turbines, and low cost/kWh. Wind is a natural and irregular source of energy that somewhat makes it less efficient than the conventional sources of energy. However, the conventional sources of energy have some issues too. The changing price of oil in world markets and climate contamination is making it less favourable for generation of energy. More specifically, the less developing world is facing severe socioeconomic problems. Pakistan is also facing similar problems of energy generation with fossil fuels.Energy generation in Pakistan is mainly achieved through the conventional resources. This is a costly business for Pakistan owing to import of the costly oil. The rise in oil prices in international market can put a huge pressure on the country to increase the prices of oil. The oil import statistics showed a rising trend at 3.8% since last two decades. During 1990 to 2014, the oil consumption increased from 28.6–67 MTOE (million tons of oil equivalent) [1]. The overall consumption of oil and gas accounted for 72%. Around 65% of energy generation is being materialized by burning the crude oil. Owing to increased use, natural gas has started to deplete. Also, energy shortages are a major cause of loss of GDP (2%) during the last decade [1]. During the year 2015-16, 111,300 GWh electricity was generated against the demand of 106,966 GWh in 2014-15 [2]. However, the consumption increased from 85,818 GWh to 90,431 GWh during last year. Pakistan has around 9.06% appropriate land area of the total which can be utilized for the development of wind farms [3].The Government of Pakistan in the year 2004 constituted the Alternate Energy Development Board (AEDB) to assess the natural resource and install the renewable energy power plants. The major cause behind the formation of AEDB was the rapid rise of prices of oil. This led to the foundation of renewable energy policy formulation and assessment of wind energy, biomass, and small hydro and solar resource assessment and later led to the installation of the wind mast throughout the country. In this regard, the external expert institutions including National Renewable Energy Laboratories helped to prepare wind mapping of Pakistan. Also, indigenous institutions like Pakistan Meteorological Department (PMD) installed wind masts at different locations of Pakistan. The geographical wind mapping and intensity at 80 m is given in Figure1. The surface roughness length of Pakistan is given in Figure 2.Figure 1
Pakistan geographical wind resource assessments at 80 m height [4].Figure 2
Geographical surface roughness length factor of Pakistan [5].The wind characteristics including speed, direction, and temperature of Kiribati were investigated in [6] in which the researchers used the Weibull parameters and found the moment method (MM) to be the accurate method. In another research work, the wind characteristics and power potential of Nooriabad (Pakistan) were investigated in [7], in which the researchers applied five methods to achieve the accuracy of Weibull parameters. The wind power potential assessment of Lithuania was examined by [8]. The authors used Weibull probability distribution function to assess the wind resource. For Chad located in North-Central Africa, the wind resource assessment has been carried out at 10 m in which the authors used the Weibull probability function [9].The continuous published research works are showing the importance of regional development of wind energy potential. Here are some of the mentioned regional studies that are showing the progress in the field of wind energy. The wind potential of Shaharbak city in Kaman province of Iran was assessed by [10]. In another resource study by [11], the author assessed the wind power potential at 10, 20, and 40 m of Yazd province, Iran. Also, another resource assessment of Iran was carried out to determine the energy potential. The considered region was the capital of Iran, Tehran. The energy potential was based on the eleven-year wind speed records calculated in [12]. The wind energy potential of Zarinah city of Iran was examined in [13], the Binalood wind resource study was investigated at 10, 30, and 40 m in [14], and Semnan province wind resource assessment study was conducted by [15].Teimourian et al. [16] studied the wind power potential of Lotak and Shandol, Iran. The authors considered the wind measurements at 10, 30, and 40 m above the ground level based on 10-minute interval of wind measurement and concluded that sites are suitable for power generation at the lowest cost. The wind resource assessment of Turkmenistan has been assessed by Bahrami, Arian. The author assessed wind power potential of 18 different locations and concluded that energy can be generated at the lowest cost [17]. In another research work, Brahami Arian et al. studied the wind power potential of 17 different locations of Uzbekistan. The author evaluated annual wind power density, energy density, and capacity factor and concluded that Nukus, Kungrad, Ak Bajtal, and Buhara are the best sites for generation of energy at lowest cost [18].The wind assessment of Jubail city, Saudi Arabia, using the 24 hourly data wind measurements at three different heights was investigated in [19]. The authors used the Weibull distribution function. Similarly, in another work carried out in the aforementioned region of Jubail, the authors took wind measurements of seven locations and investigated the k and c Weibull parameters by means of maximum likelihood, least square regression, and WAsP algorithms [20]. For selection of best possible site to install wind power plant, the authors used multicriteria decision approach (MCDA) and geographic information system (GIS) in [21]. The wind characteristics of Jeddah, Saudi Arabia, were investigated in [22]. Similarly, the wind resource assessment study for an industrial city of Yanbu, Saudi Arabia, is carried out in [23], as well as for seven different stations of Eastern province of Saudi Arabia, it is carried out in [24], and Lidar-based wind measurements are carried out in [25]. There are also some more wind potential studies of Saudi Arabia has been carried out to assess the wind energy potential considering different wind turbines [26, 27].The existing rate of potential wind speed of a site is a key factor for realizing the wind energy. Apparently, it is supposed that the wind speeds between 3–25 m/s are conducive for the conversion of wind to energy. Furthermore, the occurrences of the effective wind speeds exhibit the available rate of energy resource. The valuation of wind power using the seasonal and diurnal numbers of Waterloo of Canada was carried out in [28]. Similarly, the wind power potential assessment of Borj Cedria in Tunisia was assessed by [29] at 10, 20, and 30 m heights in which the authors projected the wind energy density by assessing the seasonal wind speeds. Also, the wind power potential of five cities including Tangier, Tetuan, Al Hoceima, Nador, and Larach of Northern Morocco was investigated by [30]. The wind characteristics of Port Said in Egypt [31] were used to investigate the energy arena. The authors used the energy flux method. The measured wind speed data were used for Tindouf in [32], which considered the eight-year data, and for Timimoun region of Algeria in [33].The wind parameters of South Banat constituency of Serbia were investigated by [34]. The authors considered the wind speed data at 10, 40, 50, and 60 m measurement heights and analysed the measured wind speed, direction, and energy density of the site. Furthermore, there are number of wind resource assessment studies available in the literature carried out at diverse parts of world including Korea [35], China [36], Malaysia [37], India [38], Egypt [39], Pakistan [7, 40, 41], and Columbia [42].In this paper, an analysis of wind characteristics is carried out to assess the wind energy potential of Gharo site. The measured wind speeds at 10, 30, and 50 m are analysed to assess the wind power prospective of the mentioned site. The wind frequency distribution including seasonal wind frequency and percentage of wind frequency has been examined to determine the accurate wind power potential of the site. In this paper, the applications of Weibull distribution function for the estimation of wind energy resource assessment of Gharo have been carried out. The annual wind power density, wind energy, and capacity factors are calculated at the three measured heights. Also, the economic assessment of the site has been assessed to check the viability of energy yield from suggested wind turbines at the lowest cost (US$/kWh).


## 2. Methodology and Data

### 2.1. Geographical Features of Site

The small town of Gharo is geographically located in Thatta district, Sindh province, Pakistan. The Gharo wind mast is situated in the peripheries of the town. The geographical view of site is shown in Figure3. The measured wind speeds were taken at three heights: 10, 30, and 50 m. The meteorology mast is equipped with NRG data acquisition system. Table 1 shows the specifications of the atmospheric sensors. The topographical location of the site is 24°35’48”N and 67°26’39”E. Table 2 shows the site-specific surface roughness values. This site can easily host the future wind farm projects. The Government of Pakistan is taking keen interest in the development of wind energy, and it established the Alternate Energy Development Board (AEDB). Also, some prescribed standards have been set by the international wind energy forum. The standards set by NREL provide the essential insights of wind energy generation and classification that is given in Table 3, and the wind turbine design standards by International Electrotechnical Commission (IEC) corresponding to the site are given in Table 3.Figure 3
Geographical location of wind site Gharo.Table 1
Technical parameters of the installed atmospheric sensors.
Wind speed sensorTemperature sensorSensorsCup type (M# 40)6 - ICT radiation plate (M # 110)Operative assortment1–90 m/s−40∼52.5 °CCorrectness±0.8%±1.1°CTemperature assortment−55∼60°C−40∼52.5°CDistance constant3.0 m—Display assortment0–120 Hz0∼2.5 V DCWeight0.2 kg0.5 kgTable 2
Typical values of surface roughness length and power law exponent.
TerrainZ0ΑIce, mud flat10−3∼3∗10−3—Calm sea2∗10−4∼3∗10−4—Sand2∗10−4∼10−30.01Mown grass0.001∼0.01—Low grass0.01∼0.040.13Fallow field0.02∼0.03—High grass0.04∼0.10.19Forest and wood land0.1∼1—Suburb or built area1∼20.32Town, city1∼4—Table 3
IEC-specified wind turbine classes [43].
Wind turbine classIIIIIISVref (m/s)5042.537.5Values specified by designerAIref—0.16—BIref—0.14—CIref—0.12—


### 2.2. Wind Power Classification and IEC Turbine Classes

To simplify the wind power density, the wind resource is divided into seven wind classes. The wind class division is basically based on energy generated from a particular wind speed. Table4 shows the classification of wind energy generation by National Renewable Energy Laboratories. Apparently, the wind classes above 4 are generally considered as feasible for installing the wind power plants. It provides two essentials including flow of wind speed and economics of energy being generated from wind speed. Wind classes 1 and 2 are generally not considered feasible for the wind power generation. The wind energy classes also provide necessary framework for economic viability of wind turbines. The wind turbine classes are given in Table 3 under which the wind turbines are manufactured by different makers. Table 3 provides the significant manufacturing framework for the wind turbines as suggested by the International Electrotechnical Commission (IEC-61400). Similarly, the typical values of surface roughness length and power law exponent are given in Table 2.Table 4
Wind energy generation classification by National Renewable Energy Laboratories [43].
Wind classResourceHeight 30 m (AGL)Height 50 m (AGL)Wind speed (m/s)Wind power density (W/m2)Wind speed (m/s)Wind power density (W/m2)1Poor00–5.10 0–16000–5.400–2002Marginal5.1–5.916–2405.4–6.2200–3003Moderate5.9–6.5240–3206.2–6.9300–4004Good6.5–7.0320–4006.9–7.4400–5005Excellent7.0–7.4400–4807.4–7.8500–6006Excellent7.4–8.2480–6407.8–8.6600–8007Excellent8.2–11.0640–1600>8.6>800


### 2.3. Assessment of Wind Data

Wind is a referred as development of air in atmosphere. It is a highly changing atmospheric parameter that changes with respect to time. It is generally accepted that the wind speed variation is better calculated using probability density function The methodology mapping of wind resource assessment of Gharo site is given in Figure4.Figure 4
Methodology mapping of wind resource assessment of Gharo site.


### 2.4. Assessment of Wind Shear

#### 2.4.1. Log Law

While considering the atmosphere, the turbulent mix can be taken in a similar way to molecular mixing. Also, this is known ask theory. Let suppose the turbulent mixing comprising of shear forces can be derived from relationship of wind speed which is given in the following equation:(1)u=U∗klnZ−DZ0,where U is the friction, k is the von Karman constant, Z0 is the roughness length, and D is the displacement height.The wind speed is computed for a reference height is expressed in the following equation:(2)UUR=lnZ/Z0lnZR/Z0,where UR refers to wind speed at the reference height ZR.


#### 2.4.2. Power Law

Generally, power law refers to increase in wind speed with height owing easier evaluation. It can be expressed as(3)UUR=Z−DZRα,where α is the power law exponent.The power law exponent can be between 0.1 to 0.32 depending upon the landscape of the site. The exponent can be calculated from the roughness length.(4)α=lnlnZ/Z0/ZR/Z0lnZ/ZR≈1lnZ⋅ZR/Z0.


#### 2.4.3. Weibull Probability Distribution Function

The Weibull distribution function is used to achieve the effectiveness of wind potential. The probability density functionf (V) and cumulative distribution function are given as follows [7, 44]:(5)fV=kcVck−1exp−Vck,FV=1−exp−Vck.Rayleigh distribution is a diverse practice of Weibull distribution function. In such instance, the value ofk parameter is considered as 2. The probability and cumulative distribution functions are given as follows [45, 46]:(6)fV=π2VVavgexp−π4VVavg2,FV=1−exp−π4VVavg2.The average wind speedVavg is mathematically expressed as follows [44]:(7)Vavg=1N∑i=1NVi.The variance and standard deviation are expressed as follows:(8)σ2=1N−1∑i=1NVi−Vavg2,σ=1N−1∑i=1NVi−Vavg2.By using Weibull parameters, the mean and variance of wind speed can expressed as [12](9)Vavg=cΓ1+1k,σ2=c2Γ1+2k−Γ21+1k,where Γ refers to gamma function and can be represented as(10)Γx=∫0∞e−uux−1du.


### 2.5. Assessment of Coefficient of Variation

The coefficient of variation of wind speed is calculated by the following equation [47]:(11)CV=σx¯.


### 2.6. Wind Power Density

The wind power densityWP is mathematically described as [48](12)WP=12ρATV3,where V is the wind speed; ρ is the air density; and AT is the rotor area of the wind turbine. The Betz number is denoted by Cp and can be referred to as [49](13)WP=12ρCpATV3.The wind power density (Wpd) obtained from equation (11) is given as follows:(14)Wpd=PA=12ρCpV3.Wpd with Weibull pdf is expressed as [50](15)Wpd=PAT=12ρV3Γ1+3k.The energy is obtained by the wind machine and is calculated as follows. While substituting equation (2) in equation below, equation (16) can be read as(16)E=T∫0∞PVfVdV,E=T∫0∞kcVck−1exp−VckPVdV.


### 2.7. Capacity Factor

The capacity factor is stated as(17)Cf=ActualPowerApRatedPowerRp.


## 2.1. Geographical Features of Site

The small town of Gharo is geographically located in Thatta district, Sindh province, Pakistan. The Gharo wind mast is situated in the peripheries of the town. The geographical view of site is shown in Figure3. The measured wind speeds were taken at three heights: 10, 30, and 50 m. The meteorology mast is equipped with NRG data acquisition system. Table 1 shows the specifications of the atmospheric sensors. The topographical location of the site is 24°35’48”N and 67°26’39”E. Table 2 shows the site-specific surface roughness values. This site can easily host the future wind farm projects. The Government of Pakistan is taking keen interest in the development of wind energy, and it established the Alternate Energy Development Board (AEDB). Also, some prescribed standards have been set by the international wind energy forum. The standards set by NREL provide the essential insights of wind energy generation and classification that is given in Table 3, and the wind turbine design standards by International Electrotechnical Commission (IEC) corresponding to the site are given in Table 3.Figure 3
Geographical location of wind site Gharo.Table 1
Technical parameters of the installed atmospheric sensors.
Wind speed sensorTemperature sensorSensorsCup type (M# 40)6 - ICT radiation plate (M # 110)Operative assortment1–90 m/s−40∼52.5 °CCorrectness±0.8%±1.1°CTemperature assortment−55∼60°C−40∼52.5°CDistance constant3.0 m—Display assortment0–120 Hz0∼2.5 V DCWeight0.2 kg0.5 kgTable 2
Typical values of surface roughness length and power law exponent.
TerrainZ0ΑIce, mud flat10−3∼3∗10−3—Calm sea2∗10−4∼3∗10−4—Sand2∗10−4∼10−30.01Mown grass0.001∼0.01—Low grass0.01∼0.040.13Fallow field0.02∼0.03—High grass0.04∼0.10.19Forest and wood land0.1∼1—Suburb or built area1∼20.32Town, city1∼4—Table 3
IEC-specified wind turbine classes [43].
Wind turbine classIIIIIISVref (m/s)5042.537.5Values specified by designerAIref—0.16—BIref—0.14—CIref—0.12—


## 2.2. Wind Power Classification and IEC Turbine Classes

To simplify the wind power density, the wind resource is divided into seven wind classes. The wind class division is basically based on energy generated from a particular wind speed. Table4 shows the classification of wind energy generation by National Renewable Energy Laboratories. Apparently, the wind classes above 4 are generally considered as feasible for installing the wind power plants. It provides two essentials including flow of wind speed and economics of energy being generated from wind speed. Wind classes 1 and 2 are generally not considered feasible for the wind power generation. The wind energy classes also provide necessary framework for economic viability of wind turbines. The wind turbine classes are given in Table 3 under which the wind turbines are manufactured by different makers. Table 3 provides the significant manufacturing framework for the wind turbines as suggested by the International Electrotechnical Commission (IEC-61400). Similarly, the typical values of surface roughness length and power law exponent are given in Table 2.Table 4
Wind energy generation classification by National Renewable Energy Laboratories [43].
Wind classResourceHeight 30 m (AGL)Height 50 m (AGL)Wind speed (m/s)Wind power density (W/m2)Wind speed (m/s)Wind power density (W/m2)1Poor00–5.10 0–16000–5.400–2002Marginal5.1–5.916–2405.4–6.2200–3003Moderate5.9–6.5240–3206.2–6.9300–4004Good6.5–7.0320–4006.9–7.4400–5005Excellent7.0–7.4400–4807.4–7.8500–6006Excellent7.4–8.2480–6407.8–8.6600–8007Excellent8.2–11.0640–1600>8.6>800


## 2.3. Assessment of Wind Data

Wind is a referred as development of air in atmosphere. It is a highly changing atmospheric parameter that changes with respect to time. It is generally accepted that the wind speed variation is better calculated using probability density function The methodology mapping of wind resource assessment of Gharo site is given in Figure4.Figure 4
Methodology mapping of wind resource assessment of Gharo site.


## 2.4. Assessment of Wind Shear

### 2.4.1. Log Law

While considering the atmosphere, the turbulent mix can be taken in a similar way to molecular mixing. Also, this is known ask theory. Let suppose the turbulent mixing comprising of shear forces can be derived from relationship of wind speed which is given in the following equation:(1)u=U∗klnZ−DZ0,where U is the friction, k is the von Karman constant, Z0 is the roughness length, and D is the displacement height.The wind speed is computed for a reference height is expressed in the following equation:(2)UUR=lnZ/Z0lnZR/Z0,where UR refers to wind speed at the reference height ZR.


### 2.4.2. Power Law

Generally, power law refers to increase in wind speed with height owing easier evaluation. It can be expressed as(3)UUR=Z−DZRα,where α is the power law exponent.The power law exponent can be between 0.1 to 0.32 depending upon the landscape of the site. The exponent can be calculated from the roughness length.(4)α=lnlnZ/Z0/ZR/Z0lnZ/ZR≈1lnZ⋅ZR/Z0.


### 2.4.3. Weibull Probability Distribution Function

The Weibull distribution function is used to achieve the effectiveness of wind potential. The probability density functionf (V) and cumulative distribution function are given as follows [7, 44]:(5)fV=kcVck−1exp−Vck,FV=1−exp−Vck.Rayleigh distribution is a diverse practice of Weibull distribution function. In such instance, the value ofk parameter is considered as 2. The probability and cumulative distribution functions are given as follows [45, 46]:(6)fV=π2VVavgexp−π4VVavg2,FV=1−exp−π4VVavg2.The average wind speedVavg is mathematically expressed as follows [44]:(7)Vavg=1N∑i=1NVi.The variance and standard deviation are expressed as follows:(8)σ2=1N−1∑i=1NVi−Vavg2,σ=1N−1∑i=1NVi−Vavg2.By using Weibull parameters, the mean and variance of wind speed can expressed as [12](9)Vavg=cΓ1+1k,σ2=c2Γ1+2k−Γ21+1k,where Γ refers to gamma function and can be represented as(10)Γx=∫0∞e−uux−1du.


## 2.4.1. Log Law

While considering the atmosphere, the turbulent mix can be taken in a similar way to molecular mixing. Also, this is known ask theory. Let suppose the turbulent mixing comprising of shear forces can be derived from relationship of wind speed which is given in the following equation:(1)u=U∗klnZ−DZ0,where U is the friction, k is the von Karman constant, Z0 is the roughness length, and D is the displacement height.The wind speed is computed for a reference height is expressed in the following equation:(2)UUR=lnZ/Z0lnZR/Z0,where UR refers to wind speed at the reference height ZR.


## 2.4.2. Power Law

Generally, power law refers to increase in wind speed with height owing easier evaluation. It can be expressed as(3)UUR=Z−DZRα,where α is the power law exponent.The power law exponent can be between 0.1 to 0.32 depending upon the landscape of the site. The exponent can be calculated from the roughness length.(4)α=lnlnZ/Z0/ZR/Z0lnZ/ZR≈1lnZ⋅ZR/Z0.


## 2.4.3. Weibull Probability Distribution Function

The Weibull distribution function is used to achieve the effectiveness of wind potential. The probability density functionf (V) and cumulative distribution function are given as follows [7, 44]:(5)fV=kcVck−1exp−Vck,FV=1−exp−Vck.Rayleigh distribution is a diverse practice of Weibull distribution function. In such instance, the value ofk parameter is considered as 2. The probability and cumulative distribution functions are given as follows [45, 46]:(6)fV=π2VVavgexp−π4VVavg2,FV=1−exp−π4VVavg2.The average wind speedVavg is mathematically expressed as follows [44]:(7)Vavg=1N∑i=1NVi.The variance and standard deviation are expressed as follows:(8)σ2=1N−1∑i=1NVi−Vavg2,σ=1N−1∑i=1NVi−Vavg2.By using Weibull parameters, the mean and variance of wind speed can expressed as [12](9)Vavg=cΓ1+1k,σ2=c2Γ1+2k−Γ21+1k,where Γ refers to gamma function and can be represented as(10)Γx=∫0∞e−uux−1du.


## 2.5. Assessment of Coefficient of Variation

The coefficient of variation of wind speed is calculated by the following equation [47]:(11)CV=σx¯.


## 2.6. Wind Power Density

The wind power densityWP is mathematically described as [48](12)WP=12ρATV3,where V is the wind speed; ρ is the air density; and AT is the rotor area of the wind turbine. The Betz number is denoted by Cp and can be referred to as [49](13)WP=12ρCpATV3.The wind power density (Wpd) obtained from equation (11) is given as follows:(14)Wpd=PA=12ρCpV3.Wpd with Weibull pdf is expressed as [50](15)Wpd=PAT=12ρV3Γ1+3k.The energy is obtained by the wind machine and is calculated as follows. While substituting equation (2) in equation below, equation (16) can be read as(16)E=T∫0∞PVfVdV,E=T∫0∞kcVck−1exp−VckPVdV.


## 2.7. Capacity Factor

The capacity factor is stated as(17)Cf=ActualPowerApRatedPowerRp.


## 3. Results and Discussion

In this paper, the wind characteristics and power assessment of Gharo site have been studied. The site is located in wind corridor (Gharo-Ketti-Bunder) of Pakistan. The small town is located in Thatta district of province of Sindh, Pakistan. According to census department of Pakistan, there were 2,000 families living in Gharo in 2017. The wind measurements taken at 10, 30, and 50 m are examined in this paper for a period of one year. Further, the results of wind resource assessment and power assessment are discussed below in the coming sections.


### 3.1. Assessment of Surface Roughness FactorZ0 and Power Law

The power law exponent varies between the two measurement heights. The power law is checked carefully because it does not have physical representation of surface layer and does not explain the flow near to the ground. Both log law and power law generally explain the wind profile of the site. The maximum and minimum values of Gharo site are given in Table5.Table 5
Maximum and minimum values ofZ0.
MaxMinMeanZ04.90.442.298


### 3.2. Diurnal Wind Variation of Site

Figures5(a)–5(c) show the seasonal diurnal wind speed variations at Gharo site over a period of year at three heights 10, 30, and 50 m. The objective of assessment of hourly mean diurnal wind variation is to determine the seasonal variation of the site at 10, 30, and 50 m. It is observed from the results that wind speed is lower during the night, whereas wind speed after sunrise picks up and levels maximum at the evening time during 4:00 pm to 5:00 pm at 30 and 50 m. It again starts to decrease as time passes to night. The visible maximum number was found in the months of March, May, June, and August at the height of 10 m. Similarly, the nature of trends is also visible in the months of March, May, June, and August at the height of 30 m. At 50 m, the maximum values are found to be in March, May, June, and August. The seasonal diurnal wind variation has been found maximum from 00:00–05:00 hours and again from 11:00–18:00 during 24 hours of a day for a period of one year. Similar trends of diurnal wind speed variation can be seen at 10, 30, and 50 m. Figures 5(a)–5(c) show the diurnal wind speed variation of site at 10, 30, and 50 m for a period of one year.Figure 5
Diurnal wind speed variation at 10, 30, and 50 m above the ground level (AGL).
(a)(b)(c)


### 3.3. Monthly Temperature of Site

Figure6 shows the monthly mean temperature of Gharo site at 50, 30, and 10 m heights. Higher values of temperature were observed to be in the months of summer, whereas lower values were observed in winter. The mean values of temperature are found to be 24.80, 24.84, and 24.98 at 50, 30, and 10 m heights.Figure 6
Monthly mean temperature of Gharo site at 50, 30, and 10 m heights.


### 3.4. Wind Speed Frequency Distribution

In other words, the wind speed frequency distribution can be plotted by different wind bins against their relative frequencies. There are two important methods including data bins of site and frequency of wind speed of site for obtaining the necessary frequency distribution.


#### 3.4.1. Data Bins

The assortment of the data into narrower level of wind speed bands is termed as the binning of the data. In this research, the width of bin is 1 m/s. For example, a measured wind speed of 3.5 m/s is placed in 3 <X < = 4 m/s bin. However, the value of each bin is 0.5 m/s, and 1.5 m/s is used in the calculations and in frequency distribution.


#### 3.4.2. Relative Frequency

The relative frequency is proportional to wind speed in each bin. The relative frequency can be viewed as the estimate of probability of given wind speed in the bin. The relative frequency can be described as follows:relative frequency =P (Vi) = frequency of given wind speed/total periodFigure7 shows the annual cumulative frequency distribution of wind speed at three heights including 10, 30, and 50 m. It is clear from Figure 7 that the values of wind speed are greater than 5 m/s for a period of 6280 hours at 30 m. Similarly, the values of wind speed are greater than 5 m/s at 50 m for a period of 6852 hours. The site assessment showed that there is enough wind to generate the energy.Figure 7
Annual cumulative frequency distribution of wind speed at Gharo site.Figure8 shows the annual wind frequency distribution of Gharo site at 10, 30, and 50 m heights. It is clear from Figure 8 that during 730 hours, wind speed is 5 m/s, during 710 hours, wind speed is 6 m/s, during 1015 hours, wind speed is 7 m/s, during 1025 hours, wind speed is 8 m/s, during 790 hours, wind speed is 9 m/s, during 740 hours, wind speed is 10 m/s, and so forth at 50 m. This analysis of measured site shows that the wind potential can be used for the installation of utility wind turbines.Figure 8
Annual wind frequency distribution of Gharo site.Figure9 shows the percentage (%) of frequency of wind speed of Gharo site. At maximum height (50 m) of measured site, we observed that the wind speed at 85% accounts for 5 m/s, that at 8% accounts for 6 m/s, that at 12% accounts for 8 m/s, and that at 9% accounts for 9 m/s. However, the measured site observation is different at three heights. The wind speed 5 m/s is observed to be 16%, 7 m/s is 8%, and 8 m/s is 15%. This indicates that the wind energy can be generated from the prevailing wind speed of Gharo site.Figure 9
Annual % of wind frequency distribution.


### 3.5. Seasonal Analysis of Wind Frequency Distribution

Figures10(a)–10(d) present the seasonal wind frequency distribution for a period of one year. Figure 10(a) shows the wind distribution for a spring season that starts from March and lasts up to May. It can be observed from Figure 10(a) that the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, and 11 m/s for a period of 330 hours, 213 hours, 212 hours, 371 hours, 185 hours, 137 hours, and 142 hours, respectively, at 30 m. Similarly, the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, 11 m/s, 12 m/s, and 13 m/s for a period of 138 hours, 187 hours, 261 hours, 296 hours, 228 hours, 230 hours, 191 hours, 99 hours, and 110 hours, respectively, at 50 m. The observed wind speed measurement shows the increase with respect to increase in height.Figure 10
Seasonal wind frequency distribution of Gharo site. (a) March–May. (b) June–August. (c) September–November. (d) December–February.
(a)(b)(c)(d)Figure10(b) shows the wind speed for a summer season. The summer season starts from the month of June and lasts up to August.The seasonal wind speed frequency showed that the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, 11 m/s, and 12 m/s for a period of 170 hours, 162 hours, 210 hours, 496 hours, 255 hours, 238 hours, 206 hours, and 230 hours, respectively, at 30 m. The observed wind speed measurement showed the strong presence of wind at 30 m. Similarly, the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, 11 m/s, 12 m/s, 13 m/s, and 14 m/s for a period of 21 hours, 64 hours, 112 hours, 198 hours, 265 hours, 277 hours, 280 hours, 305 hours, 198 hours, and 200 hours, respectively, at 50 m.Figure10(c) shows the wind frequency distribution during the period for autumn season that starts from September and lasts up to November. The measured observations showed that the 5 m/s wind speed lasts up to 435 hours, 6 m/s wind speed for 170 hours, 7 m/s wind speed for 160 hours, 8 m/s wind speed for 280 hours, and for 9 m/s is reaches to 110 hours at the height of 30 m. Similarly, for the measurement height of 50 m, the wind speed for 5 m/s is rests up to 61 hours, 6 m/s wind speed for 125 hours, 7 m/s is for 250 hours, 8 m/s is for 210 hours, 9 m/s is for 290 hours, 10 m/s is for 255 hours, and 11 m/s wind speed for 175 hours.Figure10(d) is showing the seasonal wind speed for winter season. The wind speed reached 5 m/s, 6 m/s, 7 m/s, and 8 m/s for a period of 475 hours, 170 hours, 135 hours, and 160 hours, respectively, at 30 m. Similarly, the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, and 10 m/s for a period of 66 hours, 76 hours, 155 hours, 280 hours, 208 hours, and 280 hours, respectively, at 50 m.Similarly, the wind frequency distribution seasonal percentage (%) is given in Figures10(a)–10(d) for further detailed analysis of the wind of Gharo site (Figure 11).Figure 11
Seasonal percentage of wind frequency distribution of Gharo site. (a) March–May. (b) June–August. (c) September–November. (d) December–February.
(a)(b)(c)(d)


### 3.6. Wind Speed Variation Analysis

The wind data were collected over a period of one year in 2009 and used in the analysis. The wind analysis was carried out at 10, 30, and 50 m heights. The overall mean values of wind speeds along with corresponding standard deviations and k and c parameters at the measurement heights for a period of one year are enumerated in Tables 6–8, respectively. The mean wind speeds are found to be 3.85, 5.85, and 6.89 m/s at 10, 30, and 50 m, respectively, as listed in Table 9. It can be concluded from the measured wind data that wind speed (m/s) increases with increasing altitude. Figure 12(a) shows the monthly mean variation of wind speed at the three measurement altitudes. Figure 12(a) shows that the higher values of wind speed are found to be in July while lower values of wind speed are found to be in November at three measurement heights. However, Figure 12(b) shows the ten-minute wind measurements of Gharo site. The average seasonal (winter, spring, summer, and autumn) value at 10, 30, and 50 m heights is given in Tables 9–11, respectively. The higher number of wind speeds is observed during the summer season and the lower number of wind speeds is observed during the winter season at all respective measurement heights. The concept of most likely wind speed is used to analyse the wind. The concept of most likely wind speed is used to refer the most occurring values of wind speed in the dataset. Accordingly, the highest values of 7.19, 9.96, and 11.4 m/s of most likely winds are found in July at 10, 30, and 50 m, respectively.Table 6
Mean wind speed, standard deviation, and Weibullk and c parameters for entire dataset at 50 m.
50 m heightMonthMeanΣkc (m/s)Jan4.513.161.485.01Feb5.093.301.625.68Mar5.103.121.725.74Apr7.153.272.358.06May10.63.023.9111.7Jun8.962.923.429.98Jul11.42.555.0912.4Aug7.952.603.388.80Sep8.762.304.259.64Oct4.752.661.875.40Nov3.973.111.294.30Dec4.483.351.404.90Table 7
Mean wind speed, standard deviation, and Weibullk and c parameters for entire dataset at 30 m.
30 m heightMonthMeanσkc (m/s)Jan3.552.501.483.900Feb4.102.811.604.490Mar4.222.651.684.680Apr6.452.402.897.310May9.202.803.6510.10Jun7.802.593.308.710Jul9.962.304.9610.82Aug6.882.353.197.680Sep7.602.143.988.400Oct3.822.111.894.200Nov3.122.501.303.400Dec3.502.721.293.900Table 8
Mean wind speed, standard deviation, and Weibullk and c parameters for entire dataset at 10 m.
10 m heightMonthMeanσkc (m/s)Jan1.801.881.011.79Feb2.302.281.042.30Mar2.492.201.162.71Apr4.362.351.984.90May6.592.313.107.36Jun5.762.103.066.45Jul7.191.804.667.89Aug5.101.962.795.70Sep5.491.893.206.11Oct1.901.821.021.90Nov1.491.790.891.35Dec1.792.100.891.70Table 9
Mean and maximum wind speed, standard deviation, and Weibullk and c parameters for entire dataset.
Average annual wind speed (m/s)HeightMeanMaximumΣkc (m/s)At 50 m6.8911.42.9462.6487.634At 30 m5.859.962.4892.6006.465At 10 m3.857.192.0402.0664.180Figure 12
Monthly mean wind speed (m/s) of Gharo site at 10, 30, and 50 m height and ten-minute wind measurements of Gharo site.
(a)(b)Table 10
Seasonal wind speed, standard deviation, Weibullk and c parameters, and wind power density at 50 m.
At 50MeanΣkc (m/s)W/m2Winter4.6933.271.55.19642.30Spring7.6163.1362.668.50090.40Summer9.5032.693.96310.39207.7Autumn5.8262.692.4706.466132.7Table 11
Seasonal wind speed, standard deviation, Weibullk and c parameters, and wind power density at 30 m.
At 30MeanΣkc (m/s)W/m2Winter3.7162.6761.5464.09619.00Spring6.6232.6162.7407.36363.00Summer8.2132.4133.8169.070184.0Autumn4.8462.252.3905.333102.7


### 3.7. Wind Rose Graph

For harnessing maximum possible wind energy, the appropriate knowledge of wind direction is significant. This can help to install wind turbine at the best position. The hourly wind direction is used to find out the predominant wind direction of Gharo site. The rose diagram for a period of one year is given in Figure13. However, this is obvious from the wind speed direction that the wind blows from almost all directions with larger contribution from southwest and northwest directions.Figure 13
Wind rose graph of Gharo site.


### 3.8. Coefficient of Variation

The stability of the wind speed is effective for determining the wind energy of specific site. The average coefficient of variation was found to be 0.49, 0.50, and 0.70 at 50, 30, and 10 m heights, respectively. The maximum variation was found to be 0.78 in December and the minimum variation was 0.22 in July at 50 m. At 30 m, the maximum values were found to be 0.79 in November and the minimum values were found to be 0.23 in July. At 10 m, the maximum values were found to be 1.20 in November and the minimum values were found to be 0.25 in July. The detailed coefficient of variation of wind speed at 50, 30, and 10 m is given in Table12.Table 12
Coefficient of variation of wind speed at 50, 30, and 10 m for a period of year.
Time (months)At 50 mAt 30 mAt 10 mJan0.700.701.04Feb0.640.680.99Mar0.610.620.88Apr0.460.370.54May0.280.300.35Jun0.320.330.36Jul0.220.230.25Aug0.320.340.38Sep0.260.280.34Oct0.560.550.95Nov0.780.791.20Dec0.740.771.17Thek and c parameters are used to determine the wind power prospective of Gharo. The values of the c parameter change with time and site similar to the wind speed. A number between 1 and 2 of k parameter is indicative of low winds while its increasing trend is indicative of skewness. The obtained average values of k and c parameters were 2.066 and 4.18 m/s at 10 m, 2.600 and 6.465 m/s at 30 m, and 2.648 and 7.634 m/s at 50 m. The maximum values and average values of Weibull k and c parameters at 10, 30 and 50 m are given in Tables 4, 5 and 6; respectively. The mean values of c and k parameters are observed to be increasing with increasing altitude. Significantly high monthly mean values of shape parameters are seen during July to September, but very low monthly mean values of shape parameters are seen during October to December as shown in Figure 14. The detailed seasonal values of k parameter are given in Tables 10, 11, and 13 corresponding to 50, 30, and 10 heights. However, the higher and lower values of k parameter were observed in summer and winter, respectively.Figure 14
Weibullk parameters for a period of one year at 10, 30, and 50 m.Table 13
Seasonal wind speed, standard deviation, Weibullk and c parameters, and wind power density at 10 m.
At 10MeanΣkc (m/s)W/m2Winter1.9632.0860.9801.931.000Spring4.4802.2862.0804.9918.00Summer6.0161.9533.5036.68115.7Autumn2.9601.8331.7033.1249.70With reference to scale parameter, the higher values are observed during summer and lower values are observed during winter; the details are given in Tables10, 11, and 13. A lower value of 5.196, 4.096, and 1.930 m/s is observed in winter season, whereas a higher value of 10.39, 9.070, and 3.503 m/s is observed in summer season at 50, 30, and 10 m, respectively. The mean value of c parameter is found to be 7.634, 6.465, and 4.180 m/s at 50, 30, and 10 m, respectively. Figure 15 shows the monthly detailed c parameter at 50, 30, and 10 m measurement heights.Figure 15
Weibullc parameter for a period of one year at 10, 30, and 50 m.


### 3.9. Calculations of Wind Power Density and Energy

The calculated values of wind power density and energy density are given in Table14 at the studied measurement heights. The lower and higher amount of wind power density oscillated from 223–32 W/m2, 205–14 W/m2, and 138–1 W/m2 at 50, 30, and 10 m, respectively. The average amount of wind power density at 50, 30, and 10 m is estimated to be 118.3, 92.2, and 46.08 W/m2, respectively. The detailed mean seasonal values of wind power density corresponding to winter, spring, summer, and autumn are given in Tables 10, 11, and 13 at 50, 30, and 10 m. The maximum value of the wind power density is observed during the summer and the minimum value of the wind power density is observed during winter at all measurement heights.Table 14
Wind power density (W/m2) and annual energy density (kWh/m2) for a period of one year.
Height (m)Power densityEnergy densityMaxMinMeanMaxMinMean5022332118.319552811036.6302051492.201797123807.9010138146.1012109402.60Similarly, the lowest and highest values of energy density ranged from 281 to 1955 kWh/m2, 123 to 1797 kWh/m2, and 09 to 1210 kWh/m2 at 50, 30, and 10 m heights, whereas the mean values of energy density are found to be 1036.9, 807.90, and 402.60 kWh/m2 at 50, 30, and 10 m respectively. In general, the higher values of energy density are found in summer and lower values of energy density are found in winter season at all wind measurement heights. The monthly mean variation of wind power density and energy density is depicted in Figures 16 and 17, respectively. Evidently, higher values of both the parameters are seen in summer season and relatively lower values of those are seen in winter season. This means that more energy can be generated in summer time which is also the peak load time due to increased cooling load. Furthermore, the wind power so generated can help in supplementing the peak load demands during summer time and will result in reducing the power generation from fossil fuel.Figure 16
Monthly wind power densities (WPD) (W/m2).Figure 17
Monthly energy density (ED) (kWh/m2).To estimate the monthly and annual wind energy output at the site, six wind machines of rated capacities of around 2000 to 2750 kW are chosen. Figure18 shows the power size curves of the turbines. The technical specifications are listed in Table 15. The cut-in speed of a machine is important, and its lower value is indicative of the better performance, in general. In this case, the wind turbines have the minimum cut-in speed which is 3 m/s, and hence it is expected to perform better compared to others. The total annual and monthly energy yields and maximum and minimum capacity factors (CF) obtained using the wind power curves are provided in Tables 16 and 17, respectively. The wind turbine c produced the maximum energy of 85133205 kWh with mean CF of 32.75%. However, the wind turbine E produced the minimum energy yield of 64860327 kWh during the year with an average CF of 30.91%. Comparatively, wind turbine C seems to be more appropriate for the considered site for wind power deployment. Furthermore, the seasonal trend of energy generation is suitable for meeting the local load demand which is expected to be more in summer time.Figure 18
Wind turbine power size (kW).Table 15
Characteristics of wind turbines.
Wind turbinesCut-in speed (m/s)Cut-out speed (m/s)Rated power (kW)Rated speed (m/s)Rotor diameter (m)Hub heightWT-A3252300156060WT-B3252750148060WT-C3252300158080WT-D425250011.56060WT-E4252000158080WT-F3252750159270Table 16
Wind turbine energy output for a period of one year.
Wind turbine (kWh)WT-A (2300/82.4)WT-B (2750/80)WT-C (2300/90)WT-D (2500/80)WT-E (2000/80)WT-F (2750/92)MonthskWhkWhkWhkWhkWhkWh1201008418946973676810180657220709482738829228982612731889373638126437632952202384601532945007277595137363812687826295220239042884645096864772207640063621284462128448216486510985343118528721193411311147869100022381346105169208990964973610372640916504685922311147977171154629712601938124360151180880910442865141020548775986079753538978468757878873584759789856989285139341297100938068856607837191711130133102384053220313681229452159073242345032050121113556381145631178454011896931410007180646112201008418506342621043176250920709482680556Annual684830987050035485133205670193996486032786360512Table 17
Wind turbine capacity factor for a period of one year.
Wind turbineWT-A (2300/82.4)WT-B (2750/80)WT-C (2300/90)WT-D (2500/80)WT-E (2000/80)WT-F (2750/92)MonthsCFCFCFCFCFCF11081381211214111912171631512191217164322738283534554495951575664640514249487575262546058838334535424194439504048461012915101413117509587121081381211Average (%)28.2524.4132.7525.4130.9129.75


## 3.1. Assessment of Surface Roughness FactorZ0 and Power Law

The power law exponent varies between the two measurement heights. The power law is checked carefully because it does not have physical representation of surface layer and does not explain the flow near to the ground. Both log law and power law generally explain the wind profile of the site. The maximum and minimum values of Gharo site are given in Table5.Table 5
Maximum and minimum values ofZ0.
MaxMinMeanZ04.90.442.298


## 3.2. Diurnal Wind Variation of Site

Figures5(a)–5(c) show the seasonal diurnal wind speed variations at Gharo site over a period of year at three heights 10, 30, and 50 m. The objective of assessment of hourly mean diurnal wind variation is to determine the seasonal variation of the site at 10, 30, and 50 m. It is observed from the results that wind speed is lower during the night, whereas wind speed after sunrise picks up and levels maximum at the evening time during 4:00 pm to 5:00 pm at 30 and 50 m. It again starts to decrease as time passes to night. The visible maximum number was found in the months of March, May, June, and August at the height of 10 m. Similarly, the nature of trends is also visible in the months of March, May, June, and August at the height of 30 m. At 50 m, the maximum values are found to be in March, May, June, and August. The seasonal diurnal wind variation has been found maximum from 00:00–05:00 hours and again from 11:00–18:00 during 24 hours of a day for a period of one year. Similar trends of diurnal wind speed variation can be seen at 10, 30, and 50 m. Figures 5(a)–5(c) show the diurnal wind speed variation of site at 10, 30, and 50 m for a period of one year.Figure 5
Diurnal wind speed variation at 10, 30, and 50 m above the ground level (AGL).
(a)(b)(c)


## 3.3. Monthly Temperature of Site

Figure6 shows the monthly mean temperature of Gharo site at 50, 30, and 10 m heights. Higher values of temperature were observed to be in the months of summer, whereas lower values were observed in winter. The mean values of temperature are found to be 24.80, 24.84, and 24.98 at 50, 30, and 10 m heights.Figure 6
Monthly mean temperature of Gharo site at 50, 30, and 10 m heights.


## 3.4. Wind Speed Frequency Distribution

In other words, the wind speed frequency distribution can be plotted by different wind bins against their relative frequencies. There are two important methods including data bins of site and frequency of wind speed of site for obtaining the necessary frequency distribution.


### 3.4.1. Data Bins

The assortment of the data into narrower level of wind speed bands is termed as the binning of the data. In this research, the width of bin is 1 m/s. For example, a measured wind speed of 3.5 m/s is placed in 3 <X < = 4 m/s bin. However, the value of each bin is 0.5 m/s, and 1.5 m/s is used in the calculations and in frequency distribution.


### 3.4.2. Relative Frequency

The relative frequency is proportional to wind speed in each bin. The relative frequency can be viewed as the estimate of probability of given wind speed in the bin. The relative frequency can be described as follows:relative frequency =P (Vi) = frequency of given wind speed/total periodFigure7 shows the annual cumulative frequency distribution of wind speed at three heights including 10, 30, and 50 m. It is clear from Figure 7 that the values of wind speed are greater than 5 m/s for a period of 6280 hours at 30 m. Similarly, the values of wind speed are greater than 5 m/s at 50 m for a period of 6852 hours. The site assessment showed that there is enough wind to generate the energy.Figure 7
Annual cumulative frequency distribution of wind speed at Gharo site.Figure8 shows the annual wind frequency distribution of Gharo site at 10, 30, and 50 m heights. It is clear from Figure 8 that during 730 hours, wind speed is 5 m/s, during 710 hours, wind speed is 6 m/s, during 1015 hours, wind speed is 7 m/s, during 1025 hours, wind speed is 8 m/s, during 790 hours, wind speed is 9 m/s, during 740 hours, wind speed is 10 m/s, and so forth at 50 m. This analysis of measured site shows that the wind potential can be used for the installation of utility wind turbines.Figure 8
Annual wind frequency distribution of Gharo site.Figure9 shows the percentage (%) of frequency of wind speed of Gharo site. At maximum height (50 m) of measured site, we observed that the wind speed at 85% accounts for 5 m/s, that at 8% accounts for 6 m/s, that at 12% accounts for 8 m/s, and that at 9% accounts for 9 m/s. However, the measured site observation is different at three heights. The wind speed 5 m/s is observed to be 16%, 7 m/s is 8%, and 8 m/s is 15%. This indicates that the wind energy can be generated from the prevailing wind speed of Gharo site.Figure 9
Annual % of wind frequency distribution.


## 3.4.1. Data Bins

The assortment of the data into narrower level of wind speed bands is termed as the binning of the data. In this research, the width of bin is 1 m/s. For example, a measured wind speed of 3.5 m/s is placed in 3 <X < = 4 m/s bin. However, the value of each bin is 0.5 m/s, and 1.5 m/s is used in the calculations and in frequency distribution.


## 3.4.2. Relative Frequency

The relative frequency is proportional to wind speed in each bin. The relative frequency can be viewed as the estimate of probability of given wind speed in the bin. The relative frequency can be described as follows:relative frequency =P (Vi) = frequency of given wind speed/total periodFigure7 shows the annual cumulative frequency distribution of wind speed at three heights including 10, 30, and 50 m. It is clear from Figure 7 that the values of wind speed are greater than 5 m/s for a period of 6280 hours at 30 m. Similarly, the values of wind speed are greater than 5 m/s at 50 m for a period of 6852 hours. The site assessment showed that there is enough wind to generate the energy.Figure 7
Annual cumulative frequency distribution of wind speed at Gharo site.Figure8 shows the annual wind frequency distribution of Gharo site at 10, 30, and 50 m heights. It is clear from Figure 8 that during 730 hours, wind speed is 5 m/s, during 710 hours, wind speed is 6 m/s, during 1015 hours, wind speed is 7 m/s, during 1025 hours, wind speed is 8 m/s, during 790 hours, wind speed is 9 m/s, during 740 hours, wind speed is 10 m/s, and so forth at 50 m. This analysis of measured site shows that the wind potential can be used for the installation of utility wind turbines.Figure 8
Annual wind frequency distribution of Gharo site.Figure9 shows the percentage (%) of frequency of wind speed of Gharo site. At maximum height (50 m) of measured site, we observed that the wind speed at 85% accounts for 5 m/s, that at 8% accounts for 6 m/s, that at 12% accounts for 8 m/s, and that at 9% accounts for 9 m/s. However, the measured site observation is different at three heights. The wind speed 5 m/s is observed to be 16%, 7 m/s is 8%, and 8 m/s is 15%. This indicates that the wind energy can be generated from the prevailing wind speed of Gharo site.Figure 9
Annual % of wind frequency distribution.


## 3.5. Seasonal Analysis of Wind Frequency Distribution

Figures10(a)–10(d) present the seasonal wind frequency distribution for a period of one year. Figure 10(a) shows the wind distribution for a spring season that starts from March and lasts up to May. It can be observed from Figure 10(a) that the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, and 11 m/s for a period of 330 hours, 213 hours, 212 hours, 371 hours, 185 hours, 137 hours, and 142 hours, respectively, at 30 m. Similarly, the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, 11 m/s, 12 m/s, and 13 m/s for a period of 138 hours, 187 hours, 261 hours, 296 hours, 228 hours, 230 hours, 191 hours, 99 hours, and 110 hours, respectively, at 50 m. The observed wind speed measurement shows the increase with respect to increase in height.Figure 10
Seasonal wind frequency distribution of Gharo site. (a) March–May. (b) June–August. (c) September–November. (d) December–February.
(a)(b)(c)(d)Figure10(b) shows the wind speed for a summer season. The summer season starts from the month of June and lasts up to August.The seasonal wind speed frequency showed that the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, 11 m/s, and 12 m/s for a period of 170 hours, 162 hours, 210 hours, 496 hours, 255 hours, 238 hours, 206 hours, and 230 hours, respectively, at 30 m. The observed wind speed measurement showed the strong presence of wind at 30 m. Similarly, the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, 10 m/s, 11 m/s, 12 m/s, 13 m/s, and 14 m/s for a period of 21 hours, 64 hours, 112 hours, 198 hours, 265 hours, 277 hours, 280 hours, 305 hours, 198 hours, and 200 hours, respectively, at 50 m.Figure10(c) shows the wind frequency distribution during the period for autumn season that starts from September and lasts up to November. The measured observations showed that the 5 m/s wind speed lasts up to 435 hours, 6 m/s wind speed for 170 hours, 7 m/s wind speed for 160 hours, 8 m/s wind speed for 280 hours, and for 9 m/s is reaches to 110 hours at the height of 30 m. Similarly, for the measurement height of 50 m, the wind speed for 5 m/s is rests up to 61 hours, 6 m/s wind speed for 125 hours, 7 m/s is for 250 hours, 8 m/s is for 210 hours, 9 m/s is for 290 hours, 10 m/s is for 255 hours, and 11 m/s wind speed for 175 hours.Figure10(d) is showing the seasonal wind speed for winter season. The wind speed reached 5 m/s, 6 m/s, 7 m/s, and 8 m/s for a period of 475 hours, 170 hours, 135 hours, and 160 hours, respectively, at 30 m. Similarly, the wind speed reached 5 m/s, 6 m/s, 7 m/s, 8 m/s, 9 m/s, and 10 m/s for a period of 66 hours, 76 hours, 155 hours, 280 hours, 208 hours, and 280 hours, respectively, at 50 m.Similarly, the wind frequency distribution seasonal percentage (%) is given in Figures10(a)–10(d) for further detailed analysis of the wind of Gharo site (Figure 11).Figure 11
Seasonal percentage of wind frequency distribution of Gharo site. (a) March–May. (b) June–August. (c) September–November. (d) December–February.
(a)(b)(c)(d)


## 3.6. Wind Speed Variation Analysis

The wind data were collected over a period of one year in 2009 and used in the analysis. The wind analysis was carried out at 10, 30, and 50 m heights. The overall mean values of wind speeds along with corresponding standard deviations and k and c parameters at the measurement heights for a period of one year are enumerated in Tables 6–8, respectively. The mean wind speeds are found to be 3.85, 5.85, and 6.89 m/s at 10, 30, and 50 m, respectively, as listed in Table 9. It can be concluded from the measured wind data that wind speed (m/s) increases with increasing altitude. Figure 12(a) shows the monthly mean variation of wind speed at the three measurement altitudes. Figure 12(a) shows that the higher values of wind speed are found to be in July while lower values of wind speed are found to be in November at three measurement heights. However, Figure 12(b) shows the ten-minute wind measurements of Gharo site. The average seasonal (winter, spring, summer, and autumn) value at 10, 30, and 50 m heights is given in Tables 9–11, respectively. The higher number of wind speeds is observed during the summer season and the lower number of wind speeds is observed during the winter season at all respective measurement heights. The concept of most likely wind speed is used to analyse the wind. The concept of most likely wind speed is used to refer the most occurring values of wind speed in the dataset. Accordingly, the highest values of 7.19, 9.96, and 11.4 m/s of most likely winds are found in July at 10, 30, and 50 m, respectively.Table 6
Mean wind speed, standard deviation, and Weibullk and c parameters for entire dataset at 50 m.
50 m heightMonthMeanΣkc (m/s)Jan4.513.161.485.01Feb5.093.301.625.68Mar5.103.121.725.74Apr7.153.272.358.06May10.63.023.9111.7Jun8.962.923.429.98Jul11.42.555.0912.4Aug7.952.603.388.80Sep8.762.304.259.64Oct4.752.661.875.40Nov3.973.111.294.30Dec4.483.351.404.90Table 7
Mean wind speed, standard deviation, and Weibullk and c parameters for entire dataset at 30 m.
30 m heightMonthMeanσkc (m/s)Jan3.552.501.483.900Feb4.102.811.604.490Mar4.222.651.684.680Apr6.452.402.897.310May9.202.803.6510.10Jun7.802.593.308.710Jul9.962.304.9610.82Aug6.882.353.197.680Sep7.602.143.988.400Oct3.822.111.894.200Nov3.122.501.303.400Dec3.502.721.293.900Table 8
Mean wind speed, standard deviation, and Weibullk and c parameters for entire dataset at 10 m.
10 m heightMonthMeanσkc (m/s)Jan1.801.881.011.79Feb2.302.281.042.30Mar2.492.201.162.71Apr4.362.351.984.90May6.592.313.107.36Jun5.762.103.066.45Jul7.191.804.667.89Aug5.101.962.795.70Sep5.491.893.206.11Oct1.901.821.021.90Nov1.491.790.891.35Dec1.792.100.891.70Table 9
Mean and maximum wind speed, standard deviation, and Weibullk and c parameters for entire dataset.
Average annual wind speed (m/s)HeightMeanMaximumΣkc (m/s)At 50 m6.8911.42.9462.6487.634At 30 m5.859.962.4892.6006.465At 10 m3.857.192.0402.0664.180Figure 12
Monthly mean wind speed (m/s) of Gharo site at 10, 30, and 50 m height and ten-minute wind measurements of Gharo site.
(a)(b)Table 10
Seasonal wind speed, standard deviation, Weibullk and c parameters, and wind power density at 50 m.
At 50MeanΣkc (m/s)W/m2Winter4.6933.271.55.19642.30Spring7.6163.1362.668.50090.40Summer9.5032.693.96310.39207.7Autumn5.8262.692.4706.466132.7Table 11
Seasonal wind speed, standard deviation, Weibullk and c parameters, and wind power density at 30 m.
At 30MeanΣkc (m/s)W/m2Winter3.7162.6761.5464.09619.00Spring6.6232.6162.7407.36363.00Summer8.2132.4133.8169.070184.0Autumn4.8462.252.3905.333102.7


## 3.7. Wind Rose Graph

For harnessing maximum possible wind energy, the appropriate knowledge of wind direction is significant. This can help to install wind turbine at the best position. The hourly wind direction is used to find out the predominant wind direction of Gharo site. The rose diagram for a period of one year is given in Figure13. However, this is obvious from the wind speed direction that the wind blows from almost all directions with larger contribution from southwest and northwest directions.Figure 13
Wind rose graph of Gharo site.


## 3.8. Coefficient of Variation

The stability of the wind speed is effective for determining the wind energy of specific site. The average coefficient of variation was found to be 0.49, 0.50, and 0.70 at 50, 30, and 10 m heights, respectively. The maximum variation was found to be 0.78 in December and the minimum variation was 0.22 in July at 50 m. At 30 m, the maximum values were found to be 0.79 in November and the minimum values were found to be 0.23 in July. At 10 m, the maximum values were found to be 1.20 in November and the minimum values were found to be 0.25 in July. The detailed coefficient of variation of wind speed at 50, 30, and 10 m is given in Table12.Table 12
Coefficient of variation of wind speed at 50, 30, and 10 m for a period of year.
Time (months)At 50 mAt 30 mAt 10 mJan0.700.701.04Feb0.640.680.99Mar0.610.620.88Apr0.460.370.54May0.280.300.35Jun0.320.330.36Jul0.220.230.25Aug0.320.340.38Sep0.260.280.34Oct0.560.550.95Nov0.780.791.20Dec0.740.771.17Thek and c parameters are used to determine the wind power prospective of Gharo. The values of the c parameter change with time and site similar to the wind speed. A number between 1 and 2 of k parameter is indicative of low winds while its increasing trend is indicative of skewness. The obtained average values of k and c parameters were 2.066 and 4.18 m/s at 10 m, 2.600 and 6.465 m/s at 30 m, and 2.648 and 7.634 m/s at 50 m. The maximum values and average values of Weibull k and c parameters at 10, 30 and 50 m are given in Tables 4, 5 and 6; respectively. The mean values of c and k parameters are observed to be increasing with increasing altitude. Significantly high monthly mean values of shape parameters are seen during July to September, but very low monthly mean values of shape parameters are seen during October to December as shown in Figure 14. The detailed seasonal values of k parameter are given in Tables 10, 11, and 13 corresponding to 50, 30, and 10 heights. However, the higher and lower values of k parameter were observed in summer and winter, respectively.Figure 14
Weibullk parameters for a period of one year at 10, 30, and 50 m.Table 13
Seasonal wind speed, standard deviation, Weibullk and c parameters, and wind power density at 10 m.
At 10MeanΣkc (m/s)W/m2Winter1.9632.0860.9801.931.000Spring4.4802.2862.0804.9918.00Summer6.0161.9533.5036.68115.7Autumn2.9601.8331.7033.1249.70With reference to scale parameter, the higher values are observed during summer and lower values are observed during winter; the details are given in Tables10, 11, and 13. A lower value of 5.196, 4.096, and 1.930 m/s is observed in winter season, whereas a higher value of 10.39, 9.070, and 3.503 m/s is observed in summer season at 50, 30, and 10 m, respectively. The mean value of c parameter is found to be 7.634, 6.465, and 4.180 m/s at 50, 30, and 10 m, respectively. Figure 15 shows the monthly detailed c parameter at 50, 30, and 10 m measurement heights.Figure 15
Weibullc parameter for a period of one year at 10, 30, and 50 m.


## 3.9. Calculations of Wind Power Density and Energy

The calculated values of wind power density and energy density are given in Table14 at the studied measurement heights. The lower and higher amount of wind power density oscillated from 223–32 W/m2, 205–14 W/m2, and 138–1 W/m2 at 50, 30, and 10 m, respectively. The average amount of wind power density at 50, 30, and 10 m is estimated to be 118.3, 92.2, and 46.08 W/m2, respectively. The detailed mean seasonal values of wind power density corresponding to winter, spring, summer, and autumn are given in Tables 10, 11, and 13 at 50, 30, and 10 m. The maximum value of the wind power density is observed during the summer and the minimum value of the wind power density is observed during winter at all measurement heights.Table 14
Wind power density (W/m2) and annual energy density (kWh/m2) for a period of one year.
Height (m)Power densityEnergy densityMaxMinMeanMaxMinMean5022332118.319552811036.6302051492.201797123807.9010138146.1012109402.60Similarly, the lowest and highest values of energy density ranged from 281 to 1955 kWh/m2, 123 to 1797 kWh/m2, and 09 to 1210 kWh/m2 at 50, 30, and 10 m heights, whereas the mean values of energy density are found to be 1036.9, 807.90, and 402.60 kWh/m2 at 50, 30, and 10 m respectively. In general, the higher values of energy density are found in summer and lower values of energy density are found in winter season at all wind measurement heights. The monthly mean variation of wind power density and energy density is depicted in Figures 16 and 17, respectively. Evidently, higher values of both the parameters are seen in summer season and relatively lower values of those are seen in winter season. This means that more energy can be generated in summer time which is also the peak load time due to increased cooling load. Furthermore, the wind power so generated can help in supplementing the peak load demands during summer time and will result in reducing the power generation from fossil fuel.Figure 16
Monthly wind power densities (WPD) (W/m2).Figure 17
Monthly energy density (ED) (kWh/m2).To estimate the monthly and annual wind energy output at the site, six wind machines of rated capacities of around 2000 to 2750 kW are chosen. Figure18 shows the power size curves of the turbines. The technical specifications are listed in Table 15. The cut-in speed of a machine is important, and its lower value is indicative of the better performance, in general. In this case, the wind turbines have the minimum cut-in speed which is 3 m/s, and hence it is expected to perform better compared to others. The total annual and monthly energy yields and maximum and minimum capacity factors (CF) obtained using the wind power curves are provided in Tables 16 and 17, respectively. The wind turbine c produced the maximum energy of 85133205 kWh with mean CF of 32.75%. However, the wind turbine E produced the minimum energy yield of 64860327 kWh during the year with an average CF of 30.91%. Comparatively, wind turbine C seems to be more appropriate for the considered site for wind power deployment. Furthermore, the seasonal trend of energy generation is suitable for meeting the local load demand which is expected to be more in summer time.Figure 18
Wind turbine power size (kW).Table 15
Characteristics of wind turbines.
Wind turbinesCut-in speed (m/s)Cut-out speed (m/s)Rated power (kW)Rated speed (m/s)Rotor diameter (m)Hub heightWT-A3252300156060WT-B3252750148060WT-C3252300158080WT-D425250011.56060WT-E4252000158080WT-F3252750159270Table 16
Wind turbine energy output for a period of one year.
Wind turbine (kWh)WT-A (2300/82.4)WT-B (2750/80)WT-C (2300/90)WT-D (2500/80)WT-E (2000/80)WT-F (2750/92)MonthskWhkWhkWhkWhkWhkWh1201008418946973676810180657220709482738829228982612731889373638126437632952202384601532945007277595137363812687826295220239042884645096864772207640063621284462128448216486510985343118528721193411311147869100022381346105169208990964973610372640916504685922311147977171154629712601938124360151180880910442865141020548775986079753538978468757878873584759789856989285139341297100938068856607837191711130133102384053220313681229452159073242345032050121113556381145631178454011896931410007180646112201008418506342621043176250920709482680556Annual684830987050035485133205670193996486032786360512Table 17
Wind turbine capacity factor for a period of one year.
Wind turbineWT-A (2300/82.4)WT-B (2750/80)WT-C (2300/90)WT-D (2500/80)WT-E (2000/80)WT-F (2750/92)MonthsCFCFCFCFCFCF11081381211214111912171631512191217164322738283534554495951575664640514249487575262546058838334535424194439504048461012915101413117509587121081381211Average (%)28.2524.4132.7525.4130.9129.75


## 4. Economic Assessment of Wind Turbines

This section focuses on the assessment of cost of energy generated from wind resources using different types of wind turbines. The objective of this section is to quantify the performance-based price of energy from each type of turbine. The economics of energy can be expressed as [11](18)Pcom=I1+ir−1i1+ir,PWI−t=I1+n1−irt−1ir1+irt,NPW=PWI−tt=1t1+n1−irt−1ir1+irt,TC=PWE,where E is energy production and the annual energy yield is computed from equation (19). In equation (19), Tah is total hours in a year and RP is rated power of turbine [11]:(19)E=Tah×Rp×Cf,E=1Tah1RpCf1+n1−irt−1ir1+irt.To calculate the energy output, the considered investment cost of wind turbine is US $1000/kW. The installation cost, operation cost, and interest rate is 20 %, 3 %, and 5 % of the investment cost of wind turbine, respectively. The investment cost, installation cost, operation cost, and interest rate are taken from [48, 49][51, 52]. As per cost estimation of selected wind turbines, the wind turbine WT-C has lowest cost US$ Cent 0.0346/kWh with lowest capacity factors of 0.3275 (32. 75%). The economic assessment results show that the wind turbine WT-C has an edge over the other wind turbines in terms of price per kilowatt hour. Further, the details are given in Table18.Table 18
Economic assessment/kWh (US$) energy generated by wind turbines.
Wind turbine (US$) (Cost/kWh)WT-A (2300/82.4)WT-B (2750/80)WT-C (2300/90)WT-D (2500/80)WT-E (2000/80)WT-F (2750/92)0.04300.0500.03460.04800.03940.0400


## 5. Conclusions

In this paper, the wind resource and power assessment of Gharo site is assessed at three different measurement heights including 10, 30, and 50 m above the ground level. The wind characteristics showed the strong presence at different heights. The wind power density and energy density are estimated. The present study is dedicated to the probability distribution function and frequency of the wind speed of candidate site. Additionally, the wind characteristics and power potential of Gharo site are studied while using the wind measurements at 10, 30, and 50 m heights. The overall mean wind speeds are found to be 6.89, 5.85, and 3.85 m/s at 50, 30, and 10 m, respectively, with corresponding standard deviation values of 2.946, 2.489, and 2.040. The mean values of the Weibullk parameter are estimated as 2.946, 2.489, and 2.040 while those of scale parameter are estimated as 7.634, 6.465, and 4.180 m/s at 50, 30, and 10 m, respectively. The respective mean wind power and energy density values are found to be 118.3, 92.20, and 46.10 W/m2 and 1036.6, 807.90, and 402.60 kWh/m2. As per cost estimation of wind turbines, the wind turbine WT-C has the lowest cost US$ Cents of 0.0346/kWh and highest capacity factors of 0.3278 (32.78%). Wind turbine WT-C is recommended for this site for the wind farm deployment due to high energy generation and minimum price of energy. The resource analysis shows that Gharo site has a good potential for the deployment of wind farms.



---

*Source: 8960190-2021-02-16.xml*