Next Article in Journal
Biorefining of Pigeon Pea: Residue Conversion by Pyrolysis
Next Article in Special Issue
Grid Synchronization and Islanding Detection Methods for Single-Stage Photovoltaic Systems
Previous Article in Journal
Increasing the Sensitivity of the Method of Early Detection of Cyber-Attacks in Telecommunication Networks Based on Traffic Analysis by Extreme Filtering
Previous Article in Special Issue
A Coordinated Voltage and Reactive Power Control Architecture for Large PV Power Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Design of Photovoltaic Power Plant Using Hybrid Optimisation: A Case of South Algeria

by
Tekai Eddine Khalil Zidane
1,*,
Mohd Rafi Adzman
1,2,
Mohammad Faridun Naim Tajuddin
1,
Samila Mat Zali
1,
Ali Durusu
3 and
Saad Mekhilef
4,5,6
1
School of Electrical Systems Engineering, Universiti Malaysia Perlis, Arau 02600, Malaysia
2
Centre of Excellence for Renewable Energy (CERE), School of Electrical Systems Engineering, Universiti Malaysia Perlis, Arau 02600, Malaysia
3
Department of Electrical Engineering, Davutpasa Campus, Yildiz Technical University, Istanbul 34220, Turkey
4
Power Electronics and Renewable Energy Research Laboratory (PEARL), Department of Electrical Engineering, University of Malaya, Kuala Lumpur 50603, Malaysia
5
School of Software and Electrical Engineering, Swinburne University of Technology, Victoria 3122, Australia
6
Center of Research Excellence in Renewable Energy and Power Systems, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Energies 2020, 13(11), 2776; https://doi.org/10.3390/en13112776
Submission received: 5 March 2020 / Revised: 27 April 2020 / Accepted: 30 April 2020 / Published: 1 June 2020
(This article belongs to the Special Issue Grid-Connected PV Plants)

Abstract

:
Considering the recent drop (up to 86%) in photovoltaic (PV) module prices from 2010 to 2017, many countries have shown interest in investing in PV plants to meet their energy demand. In this study, a detailed design methodology is presented to achieve high benefits with low installation, maintenance and operation costs of PV plants. This procedure includes in detail the semi-hourly average time meteorological data from the location to maximise the accuracy and detailed characteristics of different PV modules and inverters. The minimum levelised cost of energy (LCOE) and maximum annual energy are the objective functions in this proposed procedure, whereas the design variables are the number of series and parallel PV modules, the number of PV module lines per row, tilt angle and orientation, inter-row space, PV module type, and inverter structure. The design problem was solved using a recent hybrid algorithm, namely, the grey wolf optimiser-sine cosine algorithm. The high performance for LCOE-based design optimisation in economic terms with lower installation, maintenance and operation costs than that resulting from the use of maximum annual energy objective function by 12%. Moreover, sensitivity analysis showed that the PV plant performance can be improved by decreasing the PV module annual reduction coefficient.

1. Introduction

Nowadays, solar photovoltaic energy is being utilised in electrical energy generation to meet the quick-growing consumption and the urgent need for power [1]. Grid-connected photovoltaic (PV) systems with a capacity of 3 kW PV modules could meet the electric demand of a 60–90 m2 for residential building [2]. By contrast, large-scale PV power plants face some major challenges for the use of vast amounts of components in relation to the cost, reliability, and efficiency, requiring an optimal design of the PV power plant. Recently, the drop in PV module prices of up to 86% from 2010 to 2017 [3] resulted in a decrement in the levelised cost of energy (LCOE) of large-scale PV power plants reaching 0.03 ($) [4].
TRNSYS software has been used to determine the optimum PV inverter sizing ratios [5]. The simulation has been carried out using three types of inverters with low, medium and high efficiency to determine the maximum total output of the PV system. Furthermore, the PV inverter sizing ratio of the grid-connected has been investigated for eight European locations. Mondol et al. suggested that the installation of a PV system with high-efficiency inverter in the sizing of PV and inverter is more flexible than that of a low-efficiency inverter. Artificial intelligence (AI) methods have also been used to optimise the grid-connected PV power plant, as presented in [6], whereas the PV plant global solution is solved through particle swarm optimisation technique (PSO) and compared with a genetic algorithm (GA), based on the total net economic benefit. However, the PSO algorithm showed better performance than the GA approach used in this study. The optimisation design of the grid-connected PV system is introduced in [7]. The decision variables of the proposed methodology are the type of PV modules, inverter, and tilt angle. The study supports the mathematical models of the PV array, inverter and solar irradiance on tilt PV modules surface. The optimisation process considered three types of inverter, four types of PV modules and seven values of tilt angle, as well as the hourly solar irradiance and ambient temperature. As a result, the optimal design of the system is selected based on maximum efficiency.
In 2012 [8], Sulaiman, S.I., et al. proposed a sizing methodology by using an evolutionary programming sizing algorithm. The optimisation procedure supports all possible combinations of PV modules and inverters considering different types of PV modules and inverters. The technical and economic aspects are included in this method, and both the maximum yield factor and the net present value of the PV system were calculated. Chen et al. have proposed an iterative method for the optimal size of inverter for PV systems with maximum savings in nine locations in the USA [9]. The optimisation procedure has selected the gainful inverter size for each location. Additionally, optimum inverter size lower than or the same as that of PV array rated size can be installed, due to the inverter intrinsic parameters, economic and weather considerations. In 2014, Perez-Gallardo proposed an optimal configuration of the grid-connected PV power plant of different PV technologies by using the GA technique, by considering economic, technical and environmental criteria [10]. This study aims to maximise annual energy generation. Another methodology was proposed in [11] to design a PV plant for the self-consumption mechanism for different capacities in the range of 450–1250 kWp for the university campus. The simulation was performed using PV*SOL software.
A study in [12] investigated the selection and configuration of inverter and PV modules for a PV system for minimising costs. The purchasing costs can be reduced by 16.45% of 10 kW by using this model. However, this evaluation model is applicable only at the lowest price and cannot be applied to achieve the highest efficiency in power production. A mathematical procedure is presented in [13,14,15] to determine the optimal number of rows and a PV module tilt angle for maximising the profit during PV plant lifetime, by considering the effect of shading on the PV module output power. A work in [16] investigated the design of PV systems grid-connected, considering the PV module degradation rate, to select the optimum inverter size for increased energy and reduced cost. Actual inverters with high efficiency offer a wider range than inverter with low-efficiency for sizing factor to increase the energy generation. Research presented in [17] proposed an eco-design for grid-connected PV systems, on the basis of the combination of multi-objective optimisation and other software. The techno-economic and environmental criteria were optimised simultaneously. The installation of thin-film PV modules in PV systems show an advantage over crystalline silicon ones. A methodology for achieving the optimal configuration of large-scale PV power plants to improve its performance is presented in [18]. The optimisation process was performed using different algorithms and is considered to minimise the LCOE by using crystalline silicon and thin-film cadmium telluride PV module technology. According to this study, the proposed technique of grey wolf optimiser showed improved results compared with the other methods in solving the optimal design of the PV power plant. The PV plant LCOE with the thin film had a lower value than crystalline silicon and is more productive. The work reported in [19] proposed a method to convert the design of PV power plants to binary linear programming to achieve an economical design. However, in this method, only the number of inverters and PV modules connected in series and parallel were considered as the design variables. Some other methods have also been employed and published by researchers in this topic to propose a suitable configuration and determine the best solution that considers the environment, economic and technical aspects of the PV system [20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36]. Additionally, references [37,38,39,40,41] reviewed the grid-connected PV system optimisation and challenges.
The average time for the input meteorological data is an essential factor in PV system design, because the monthly and daily average time of the meteorological data fails to determine an optimum design, resulting in the oversizing system and high energy losses and increasing the financial risk of the PV plant. Additionally, the geographic latitude of the PV plant installation site can lead to a significant variation of the PV module optimal tilt angle from one location to another, to convert maximum solar irradiance into electricity and make the PV system more profitable [42].
This paper intends to present a methodology for designing PV power plants by considering semi-hourly time-resolution (i.e., 30 min-average) to address the accuracy of the meteorological data variation, and thus determine the PV plant optimal design and increase its performance. The procedure considers the detailed specifications of the different alternatives of PV modules and inverters to determine the optimum component and system topology for the location under study. Three meteorological parameters of solar irradiation, wind speed, and ambient temperature were measured for 1 year at the installation field are considered. Hybrid grey wolf optimiser-sine cosine algorithm (HGWOSCA) [43] and sine cosine algorithm (SCA) [44] were applied as optimisation techniques to solve the PV plant design problem for two different objectives, including minimum levelised cost of energy (LCOE) and maximum annual energy, while considering many design variables for improving the system performance. The contributions of this article to the book of knowledge in this research field are described below.
  • The proposed methodology is suitable to be executed using semi-hourly time-resolution (i.e., 30 min-average) values of the meteorological input data in designing the PV power plant and by introducing an actual PV plant field model, by considering the shape and size of the PV power plant installation area, to arrange all the existing components properly.
  • The application of a HGWOSCA optimisation approach after the consideration of two objective functions to design the PV power plant was presented.
  • A sensitivity study was performed to investigate the effect of the annual PV module reduction coefficient on PV plant performance.
  • A review of the Algerian renewable energy target and its integration was presented.
This paper is organised as follows: Section 2 presents an overview of the renewable energy potential of Algeria, and Section 3 presents the work methodology, including the formulation of the design problem, the PV system description and meteorological data and, the proposed design optimization. In Section 4, the HGWOSCA algorithm is described. Section 5 presents the obtained results with the sensitivity study. Finally, Section 6 presents the conclusions of the paper.

2. The Renewable Energy Potential of Algeria

Algeria has an important potential for electricity generation from renewable energy sources, as performed in several recent studies. However, according to reference [45], approximately 0.415% of electricity in Algeria is generated from renewable energy sources in 2014. The diesel generator is the dominant energy source in rural and Saharan regions in Algeria [46].

2.1. Solar Energy

The potential of solar energy in southern Algeria is the largest in all Mediterranean basins, with 1,787,000 km2 of Sahara desert, according to the German Aerospace Centre (DLR). The insolation time of almost all the national territory exceeds 2000 h annually and reaches 3900, as shown in Figure 1 (high plains and Sahara) [47]. Over most of the country and during the day, the energy obtained on a horizontal surface of 1 m2 is nearly 5 kWh or about 2263 kWh/m2/year in the south and 1700 kWh/m2/year in the north [48]. This great potential in solar energy compels Algeria to go towards the exploitation of solar energy for power generation, rather than oil and gas. Table 1 shows the rate of sunshine for each region of Algeria.

2.2. Other Renewable Energy Potential

Algeria has a considerable average of wind speed that can reach 6 m/s in approximately 50% of the country’s surface. The government plans to promote this energy source. Hydroelectricity potential is modest and has little benefit to the Algerian economy. For geothermal energy, the Renewable Energy Development Centre listed more than 200 hot springs, thereby presenting a favourable outlook to exploit this resource [48]. Finally, the biomass potential in Algeria is considerable. It reaches 5 Mtoe/year, but it is still not consumed. This source of energy offers great promises and should be enhanced [51].

2.3. Algerian Renewable Energy Program

Currently, Algeria focuses on the production of electrical energy from renewable sources, especially on solar energy and grid-connected photovoltaic power plants, for a capacity of several megawatts to reduce its dependence on oil rents, thereby representing approximately 96% of export and 17.36% of GDP in 2014 [52]. It also aims to reduce the use of diesel generators, which is the dominant energy source in rural and Saharan regions in Algeria [46].
The development of green energy sources is a key priority for the Algerian government. Nowadays, enhancing the exploitation of renewable energy is a necessity to decrease the amount of CO2 emissions, considering that Algeria has shown great interest in signing the historic Paris agreement in 2015 [53]. Furthermore, with such a measure, the government can save conventional resources, which are used to generate electricity. The renewable sources showed a poor share in the total energy compared with the conventional sources [54]. The residential electricity sector reached approximately 42% of the total energy consumption [55]. However, the ambitious national renewable energy program allowed one to reach 27% of renewable energy in the national energy mix [56].
Renewable energy sources are the focus of the 2011–2030 development program adopted by the Algerian government to achieve the installation of 22,000 MW of renewables by 2030, including 10,000 MW for export and 12,000 MW for meeting the national market demand [57]. Notably, photovoltaic energy is the dominant renewable source and is expected to reach a capacity of 13,575 MW, representing 62% of the total power installation, as shown in Table 2. In this ambitious program, the government strategy focuses on the development of photovoltaics on a large scale and prepares for the future of Algeria. The adopted program includes the development of wind and CSP energy, biomass, cogeneration and geothermal sources.

2.4. Photovoltaic Power Plants Installed in Algeria

In August 2019, the Algerian company for electricity and gas (Sonelgaz) signed an agreement with five companies to construct nine PV power plants in the southern big Sahara region, with a total capacity of 50 MW [59]. This project aims to make a hybrid energy system with an existing gas turbine and diesel generator. The PV plants are installed in different locations, as shown in the following Table 3:
Several large-scale PV power plants have been installed and connected to the electric grid in different locations across the country, and most of the PV plants are located in the big Sahara in the South, thereby indicating that Algeria benefits from this source of energy, as shown in Table 4. In addition to the absence of batteries that reduce the total capital cost, this system allows power generation surplus in terms of consumption of the load to be automatically injected into the electric grid.

3. Methodology

3.1. Formulation of the Design Problem

The single objective optimisation function is used to find the optimum solution corresponding to the minimum or maximum value defined by the objective function. In contrast, multi-objective optimisation combines two or more individual objective functions to determine a set of trade-off solutions, which allow decision makers to select the most suitable solution based on the problem requirements [60]. In this study, in sizing optimisation methodology depending on the requirements of the power plant designer, each of the two objectives can be used to produce an optimal design for the PV power plant. In addition, for comparison purposes, the optimum values are calculated by using each objective function individually to evaluate the PV power plant performance.
Furthermore, multi-objective optimisation can be used in the design of PV systems with a small capacity in the range of kW, with a small number of PV modules and inverters or in hybrid renewable energy systems for example (PV-wind) or (PV-diesel denerator-battery). However, in large scale PV power plants (i.e., >200 kW nominal power rating—the largest plants reaching several tens of MW of capacity), with a considerable number of components required in PV plant installation, it is well-known that the levelised cost of energy (LCOE) is applied to enable the reduction of the PV plant cost per watt of nominal power that is installed [61,62], for this reason, single objective optimisation is used. Additionally, a recent study is presented in [63] to investigate the LCOE of large scale PV power plants at 8 PV plants ranging from 1 to 46 MWp and many similar studies can be found in the literature.
In this section, two objective functions are considered to evaluate the PV power plant performance and to solve its complex design problem. The design variables and constraints of the proposed methodology are also explained.

3.1.1. Objective Function

In this work, the LCOE and maximum annual energy were set as objective functions to determine the optimal solution of the PV plant design. These two objective functions can be combined to form a single optimisation function.
The first part presents the LCOE which is calculated on the basis of the sum of maintenance, operation and installation costs of the plant divided by the total energy generation of the plant during its lifetime. The LCOE method is generally applied to compare power plants with different energy generation sources, by considering the appropriate cost structures. However, the best LCOE for power plants presents the lowest possible investment with high annual energy production. The second part presents the maximum amount of annual energy that can be captured by the PV modules during the PV plant in its lifetime, which is 25 years. The single optimisation function is expressed by the following equation:
X ¯ m i n [ ( C c ( X ) + C M ( X ) E t o t ( X ) ) · a ( ( 1 a ) · ( P p l a n t ( X ) · n s · E A F ) ) ]
where n s is equal to 1 year.
The optimum values are calculated by using each objective function individually. In other words, in the objective function, a is a binary number; if a is equal to 0, the target of the objective function is maximum energy and, if a is equal to 1, the objective function target is minimum LCOE.

3.1.2. Design Variables

The proposed optimisation algorithm was used for the calculation of all the decision variables, to determine the optimum design of the PV power plant. The chosen optimisation algorithm should have high performance in determining the best design variables and solving the design problem. In this methodology, the proposed decision variables, including the number of PV modules connected in series ( N s ) and parallel ( N p ) , number of PV module lines per row ( N r ) , the distance between two adjacent rows ( F y ) , the tilt angle of the PV module ( β ) , the orientation of PV modules ( P V o r i e n ) , that can be installed vertically or horizontally, optimum PV module ( P V i ) , and inverter ( I N i ) , can be selected on the basis of several alternatives from a list of possible candidates.
The vector of the decision variables are summarized as given by the following expression:
X = [ N s   N p     N r   β   F y P V o r i e n     P V i     I N i ]

3.1.3. Constraints

During the design of the PV power plant, many constraints are considered to account for the limits of the different parameters of the whole system. The following expression shows the limitation of some variables:
N s , m i n N s N s , m a x
1 N p N p , m a x
1 N r N r , m a x
0 ° β 90 °
S o c u p i e d A v a l a i b l e a r e a
The following equality constraint expressions were used to select the PV module and inverter from the list of candidates:
P V 1 + P V 2 + = 1
I N V 1 + I N V 2 + = 1

3.2. System Description and Meteorological Data

This paper focuses on the optimal design of large-scale PV power plants connected to the electric grid. The optimisation process was carried out using HGWOSCA and SCA techniques. The proposed methodology supports PV plant configurations for both central and string topologies and the arrangement of the components within the installation area. Moreover, it considers a list with different types of PV modules and inverters and their specifications as candidates to design the PV power plant, as illustrated in Table 5 and Table 6. The actual meteorological data, such as solar irradiance, wind speed, and ambient temperature, are also considered. The optimisation process in this methodology aims to determine the number of PV modules connected in series and parallel, the number of PV module lines per row, the distance between two adjacent rows, the tilt angle of the PV module, the orientation of PV modules that can be installed vertically or horizontally, the optimum PV module and inverter, from a list of possible candidates. Two objective functions, namely the minimum LCOE and maximum annual energy, were considered. Furthermore, the effect of the annual PV module reduction coefficient on PV plant performance was determined.

3.2.1. System Description

The architecture of PV plants is composed of several hundreds of PV modules that produce DC power depending on the meteorological parameters (solar radiation and temperature) [64] and inverters that permit the conversion from DC to AC power and ensure the maximum power extraction from the PV modules [65]. The generated power is injected directly into the electric network at the point of common coupling by using step-up transformers [66].
In PV plants, PV modules are connected in series ( N s ) to form a string. These strings are rationalised and connected in parallel ( N p ≥ 1) to form a PV array. In the case of central inverters topology, several hundreds of PV modules are connected to one inverter, and junction boxes need to be used through the DC main cable before leaving for the inverter. In string inverter topology, one string is connected directly to one inverter, and junction boxes need not be placed in the installation field.
The PV arrays within the available area are arranged in multiple rows. Each row consists of multiple lines of PV modules. The number of lines per row is equal to N r . Considering the shading effect, the inclined adjacent rows are installed with space in between. The tilt angle (β) of PV modules is constant during the PV plant lifetime.
The universal transverse Mercator coordinate (UTM, X: east, Y: north) system was used to model the PV power plant area. Several PV modules vary from one row to another, as the row length varies on the basis of the installation area shape.

3.2.2. Selected Site

The Djanet village is located in the South East of Algerian Sahara in the province of Illizi, and it is characterised by a hot desert climate according to longitude 9.28° E, latitude 24.15° N and an altitude of 1030 m. Djanet village electrification is still dominated by diesel generation (DG). The delivery of fuel leads to an increase in the cost and the maintenance of the units. However, a PV plant grid connected with a capacity of 3 MW was installed in 2015, and its extension is expected to reach 7 MW to supply the village and decrease the units of DG.

3.2.3. Meteorological Data

The performances of PV modules depend on solar irradiation, ambient temperature, and wind speed. These data have been recorded to step times of semi-hourly and hourly. Solar irradiation is expressed in Watts per square meter (W/m2). Ambient temperature is in degrees Celsius (°C). Wind speed is in m/s. Figure 2 and Figure 3 show the daily assessment of the meteorological characteristics. The semi-hourly data of solar irradiance and ambient temperature are observed as more accurate than the hourly data. That is because the peak values of the solar irradiance and ambient temperature may not be recorded in the hourly data. This happens because the data only records point values for every hour, while within a period of time the meteorological data may have significant fluctuation. Thus, semi-hourly data is more precise and accurate compared to hourly data measurements.
To determine the energy generated by the PV power plant, knowledge of the solar irradiance profile during the year is required. The solar irradiation for hourly and semi-hourly (i.e., 30 min-average) data of the selected location are plotted in Figure 4 and Figure 5. As we can see, the radiation intensity remained high over the year. It is observed that the maximum value of solar irradiance is reached in March. As for the minimum value of irradiance, this is recorded in June. Each vector represents the same size for 1 year, where the one for semi-hourly step is equal to 17,520 data and the one for hourly step is equal to 8760 data. The climatic conditions of the location are as follows: high solar irradiance potential; ambient temperature with a maximum average of (29.7 °C) in summer; and the sky is mostly clear during the whole year.

3.3. Proposed Design Methodology

The proposed methodology was applied to solve the PV power plant design and aimed to determine the optimal sizing and configuration of the PV plant, as shown in Figure 6. In the following section, the PV plant design parameters were calculated step by step by considering the measured meteorological data of the location, PV modules, inverter specifications, area coordinates, and cost units.

3.3.1. Irradiance Model

Solar irradiance on tilted PV modules surface is a very important factor in the optimal design of PV power plants. Installation areas with PV modules facing south are suitable for PV power plants [35]. Additionally, a good agreement was shown by PV plants oriented towards the south by using the isotropic model [67].
Several models, which are classified as isotropic and anisotropic, can be used to estimate the solar irradiance on a tilted plane [67,68]. The treatment of diffuse radiation was the only difference between these two models, while the rest is treated the same. However, the isotropic model developed by Liu and Jordan [69] was applied to this methodology to estimate the solar irradiance on the tilted PV module surface [67].
The total radiation received on a horizontal surface (global radiation: I) can be divided into two components: beam and diffused radiation. The estimation of solar radiation on the tilted surface is calculated on the basis of these two components. The total radiation received on the horizontal surface is given by the following equation:
I ( t ) = I b ( t ) + I d ( t )
The index of transparency of the atmosphere or the clearness index k T of the sky is an essential factor. The clearness index is the function of the ratio between the extraterrestrial and horizontal radiation, as expressed by the following equation:
k T ( t ) = I ( t ) I 0 ( t )
The diffuse fraction of total horizontal radiation depends on the clearness index of the sky [67] and is expressed by the following equation:
I d ( t ) I ( t ) = { 1.0 0.09 k T ( t ) , k T ( t ) 0.22 0.9511 0.1604 k T ( t ) + 4.388 ( t ) k T 2 16.638 k T 3 ( t ) + 12.336 k T 4 ( t ) ,   0.22 < k T ( t ) 0.8 0.165 , k T ( t ) 0.80
manipulating Equation (10), the beam radiation is given by the following expression:
I b ( t ) = I ( t ) I d ( t )
The total incident solar radiation on tilted surface is the sum of three components, namely, beam radiation from direct radiation of the inclined surface, diffuse radiation and reflected radiation.
I T ( t , β ) = I B ( t ) + I D ( t ) + I R ( t )
The beam irradiance on an inclined surface can be calculated on the basis of multiplication between beam horizontal irradiance and beam ratio factor R b , as shown in the following expression:
I B ( t ) = I b ( t ) R b ( t , β )
where the beam ratio factor R b is a function of the ratio between beam irradiance on the inclined surface and horizontal irradiance, as expressed in the Equation (18).
The first component is the incidence angle cos ( t , β ) , which can be derived as follows:
cos ( t , β ) = sin δ ( t ) sin φ cos β sin δ ( t ) cos φ sin β cos γ + cos δ ( t ) cos φ cos β cos ω ( t ) + cos δ ( t ) sin φ sin β cos γ cos ω ( t ) + cos δ ( t ) sin β sin γ sin ω ( t )
where δ is the solar declination angle, φ is the location latitude, γ is the surface azimuth angle, and ω is the hour angle. The global radiation on the inclined surface calculation model’s error was lower than 3% [70].
The second component deals with solar zenith angle cos θ z and can be calculated using the following equation:
cos θ z ( t , β ) = cos γ cos δ ( t ) cos ω ( t ) + sin φ sin δ ( t )
R b ( t , β ) = cos ( t , β ) cos θ z ( t , β )
Diffuse irradiance on an inclined surface is computed on the basis of the isotropic sky model. A well-known isotropic model was introduced by Liu and Jordan (1963). This model is simple, and the diffuse radiation has a uniform distribution over the skydome. The diffuse radiation on the inclined surface increases with an increasing amount of seen by the inclined surface, as expressed in Equation (19).
I D ( t ) = I d ( t ) ( 1 + cos β 2   )
where β is the surface tilt angle and considered as a design variable. Its optimal values are computed by the optimisation algorithm.
The reflected irradiance on an inclined surface is expressed by Equation (20) and depends on the transposition factor for ground reflection R r given by Equation (21) and the reflectivity of the ground ρ that is equal to 0.2 [68].
I R ( t ) = I ( t ) ϱ R r
R r ( t ) = 1 cos β 2

3.3.2. Area Calculation Model

In actual cases, the PV power plant installation area is limited in surface and does not have a uniform shape. However, this proposed methodology can be applied to all actual area shapes to determine the optimal size and configuration of large-scale PV power plants. This methodology supports actual area shapes by using the coordinates of the location under study. Additionally, the universal transverse Mercator coordinate (UTM, X: east, Y: north) system is used to model the PV power plant area. Furthermore, the PV plant occupying the surface, length, and width of each row, junction boxes, and cable length are computed on the basis of the coordinates. As mentioned in the previous section, PV nodules are oriented towards the south in the installation field, as illustrated in Figure 7.
The Y-axis was used to calculate the total number of rows supported by the PV plant area and the width ( W T ) . of each row. The space between two adjacent rows ( F y ) , which is a design variable, and its optimum value were calculated by the algorithm process. The following equation expresses the total number of rows:
N r o w = f l o o r ( max ( Y ) m i n ( Y ) W T + F y   )
After calculating the total number of rows ( N r o w ) , it can be used to calculate the north coordinates of each row ( Y N i ) in the PV plant, as expressed in the following equation:
Y N i = min ( Y ) + ( W T + F y ) N i
The parameters ( N i ) and ( Y N i ) are subject to a constraint in this methodology, as provided in Equations (24) and (25):
1 N i N r o w
m i n ( Y ) < Y N i < max ( Y )
The X-axis presents the east coordinates and is used for calculating the length of each row in the PV plant X N i , as expressed in the following straight-line equation:
X N i = ( X 2 X 1 ) ( Y N i Y 1 ) Y 2 Y 1 + X 2
where ( X 1   , Y 1 ) . and ( X 2   , Y 2 ) correspond to the coordinates of two consecutive points. The parameter X N i is a constraint in this methodology, as provided in the following expression:
m i n ( X ) X N i max ( X )
The row length ( M r o w i ) is obtained after the calculation of the east coordinates ( X N i ) of each row, considering the difference between these coordinates, and is expressed by using Equation (28):
M row i = X N i 1 + X N i 2
The PV power plant area calculation process considers other important parameters, such as row height ( H T ) , row width ( W T ) and the space between two adjacent rows ( F y ) . These parameters can be calculated on the basis of the following equations:
W T = N r L p v , 2 cos β
H T = N r L p v , 2 sin β
F y = d H T
where ( N r ) and ( F y ) are considered as design variables, and their optimal values are calculated via optimisation. Notably, ( N r ) is the number of PV module lines in each row, and ( F y ) is the distance between two adjacent rows. In this methodology, all rows in the installation area have the same lines of PV modules. The arrangement of rows and PV modules in a row within the installation area is shown in Figure 7.

3.3.3. Components Arrangement

The arrangement of the components within the installation area is an essential part of the PV plant design process in the presence of several parameters, such as the location characteristics and the device’s specifications. In addition, component arrangement depends on the optimal topology selected by the optimisation algorithm. Furthermore, the distribution of a large amount of the components among the PV power plant is computed in terms of several constraints.
However, PV modules and inverters are the two main devices considered in the PV power plant arrangements. Additionally, in case of the optimisation algorithm select central topology, the junction box arrangement is considered, and its distribution among the PV modules and the inverters is calculated on the basis of its rating power.
Finally, the PV power plant device arrangement is influenced by the amount of solar irradiance, ambient temperature, wind speed, and the geographic location. These parameters affect the tilt angle of PV modules and increase or decrease the PV module energy output, leading to the installation of varying numbers of inverters in the PV plant. Moreover, in this methodology, the aforementioned parameters are considered to control the total cost.
Dependent on PV inverter size, the number of series PV modules in each string ( N s ) and parallel PV modules ( N p ) should be computed by the algorithm to meet a specific voltage and current requirement of inverters. On the one hand, to avoid the inverter damage that can be caused by overvoltage in case of low temperature in some locations, in every string, the number of PV modules connected in series has to be optimally computed. On the other hand, the number of parallel-connected PV modules ( N p ) multiplied by its current is equal to the input current of the inverter. To avoid the inverter damage created by the overcurrent locations with high solar irradiance, a limited number of PV modules connected in parallel ( N p ) should be addressed.
The first part handles PV modules distribution among the inverters and their arrangement within the PV plant area. The number of series ( N s ) and parallel ( N p ) PV modules are computed in accordance with the optimum selected inverter by the optimisation process. In this proposed methodology, the number of PV modules connected in series ( N s ) and parallel ( N p ) were considered as the design variables, and their optimum values were calculated using the optimisation algorithm. The ( N s ) design variable involves a number of minimum ( N s , m i n ) and maximum ( N s , m a x ) PV modules, and these limitations can be calculated on the basis of the inverter input voltage range in [11,22], as expressed in the following equations:
N s , m i n = V i , m i n V m p p , m i n
N s , m i n = V i , m a x V o c , m a x
N s m , 2 = V i , m p p t m a x V m p p , m a x
N s , m a x = { N s m , 1 , N s m , 1 N s m , 2 N s m , 2 ,   N s m , 2 < N s m , 1
The maximum number of PV modules connected in parallel ( N p ) was calculated according to the selected inverter by using the nominal power ( P i ), and the PV module maximum output power ( P m p p , m a x ) was selected with respect to the optimum number of PV modules connected in series ( N s ) [22], as provided in the following expression:
N p , m a x = P i N s P m p p , m a x  
As mentioned in the previous section, the arrangement of PV modules in the PV plant area requires the use of the length of each row in the PV plant to determine the optimum number of PV modules installed in each line ( N c i )   and the total number in each row ( N r o w i ,   p v ). The total number of PV modules installed in each line ( N c i ) of rows, which are described as the function ratio between the length of each row ( M r o w i ) and the length of the optimum PV modules ( L p v , 1 ) , is given in the following equation:
N c i = M r o w i L p v , 1
The total number of PV modules installed in each row ( N r o w , p v ) depends on the number of PV module lines ( N r ) , which is a design variable in this methodology, and its optimum value is computed by the optimisation algorithm.
N r o w i , p v = N r N c i
The sum of PV modules in each row of the PV plant results in their total number in the installation area as expressed in the following equation:
N I = 1 i N r o w i , p v
The number of series ( N p ) and parallel ( N p ) PV modules are the main parameters in the inverter calculation process. These design variables determine the number of blocks, and x i n v represents the pieces of inverters in blocks, and each piece is composed of N b l o c k [22], as given in the following equations:
N b l o c k = N s N p
y = ( N i   , N b l o c k )
x inv = N i y N b l o c k
Finally, the total number of inverters is calculated on the basis of the following expression:
N i = { x i n v , ( y x i n v ) P p v , s t c 0.1 P i x i n v + 1 , ( y x i n v ) P p v , s t c > 0.1 P i

3.3.4. PV Plant Total Energy

The proposed methodology offers many alternatives for PV modules with different specifications. Additionally, the optimisation algorithm was applied to determine the best candidate for the design of the PV plant and the optimum configuration of the PV plant as a global solution. However, the PV module output power depends on the amount of solar radiation, ambient temperature, wind speed, and electrical characteristics. Moreover, a recent review [71] has covered approximately 70 important papers on PV cell modelling, and the equations used in this proposed methodology have been applied in several papers, as shown in this review. The equations have been used in a recent paper [72], and the obtained results by the proposed procedure are more accurate than the [73] model, which involves the use of the same equations. Accordingly, these equations are suitable for calculating the performance of PV modules in our proposed design procedure.
The PV power plant consists of a large number of PV modules. Additionally, the output power is assumed to be the same for all PV modules in the PV plant, except for the southernmost row, which is considered never shaded. More importantly, the degradation of PV modules is inevitable regardless of the size of a PV power plant [74,75]. However, this research considered the PV module output power derating factor ( d f ) due to soiling effect on the PV module surface, which is equal to d f = 0.069, and the annual reduction coefficient r of PV module [34], which is equal to 0.5%. Finally, PV modules output power can be calculated using the following expression:
P pv ( t , β ) = ( 1 r ) ( 1 d f ) P m p p ( t , β )
where P m p p presents the produced power by each PV module in the PV plant.
The produced energy can be affected by the shadow area on PV modules and is related to the shade impact factor (SIF) [76], and its value is equal to 2 [35]. This parameter can be obtained using the following equation:
A S i ( t ) = ξ i ( t ) S I F
where ( ξ i ( t ) ) presents the ratio of the shadow area.
The total energy of the PV power plant can be calculated according to the optimum inverter topology selected by the optimisation algorithm. Furthermore, the PV power plant produced energy and the total cost can be influenced by the selected inverter topology. For string inverter topology, the following equation is applied to calculate the PV plant output power:
P plant ( t , β ) = n t r ( 1 η cac ) ( 1 η cic ) P o i ( t , β ) N i
where ( P o i ) is the inverter output power, ( N i ) represents the total number of inverters, ( n t r ) is the transformer efficiency, ( η cac ) presents the AC cable losses and ( η cic ) is the interconnection cable losses.
In the case of central inverter topology, PV plant output power can be obtained using the following equation:
P plant ( t , β ) = ( 1 η c d c ) n m p p t n i n v n t r ( 1 η cac ) ( 1 η cic ) 1 r o w i P r o w i ( t , β )
where P r o w i ( t , β ) presents the PV row output power, n m p p t , n i n v and n t r are the efficiencies of the PV module, inverter and transformer, respectively, and η c d c and η cac , are the DC and AC cable losses, respectively.
However, in this methodology, the PV plant energy generation was directly injected to the electric network over its operational lifetime, and it was calculated using Equation (48):
E t o t = P plant ( t , β ) n s E A F
where EAF is the energy availability factor, and ( n s ) is the PV plant operational lifetime.

3.3.5. PV Plant Total Cost

The PV power plant consists of two types of costs, as expressed by Equation (49):
C t o t = C c + C M
The installation cost ( C c ) deals with the cost of the device, such as C p v , C i n v which represents the unit cost of the PV modules and inverters, respectively. In addition, C B is the PV module mounting structure cost. Moreover, C c b , C t r , C p d , and C c m represent the costs of the cable, transformer, protection devices and monitoring system, respectively. Finally, C L represents the cost of the plant area. The installation cost is expressed in Equation (50):
C c = N I C p v + N i C i n v + C L + C B + C c b + C t r + C p d + C c m
The operation and maintenance costs of the PV plant during its lifetime depend on the annual inflation rate (g), the nominal annual interest rate ( i r ) . and the operation and maintenance costs per watt ( M o p ) , as given in the following expression:
C M = N I P p v , s t c M o p ( 1 + g ) [ 1 ( 1 + g 1 i r ) n s i r + g ]

4. Hybrid Grey Wolf Optimizer-Sine Cosine Algorithm (HGWOSCA)

SCA and grey wolf optimiser (GWO) are meta-heuristic optimisation algorithms recently developed by Mirjalili et al. [44,77]. Both SCA and GWO approaches show high performance compared with other well-known meta-heuristic algorithms [44,77]. The hybrid GWO-SCA technique was introduced by N. Singh et al. [43] for combining the advantages of both approaches. In the GWO-SCA hybrid approach, GWO presents the main part, whereas the implementation of SCA assists in the optimisation of GWO. An improvement in the position, speed, and convergence of the best grey wolf individual alpha (α) by using the original equation expressed in [77], is achieved by applying the position updating equations of the SCA approach, as illustrated in [44].
The position of the current space agent is updated on the basis of the following equation:
x 2 = x β a 2 · ( d β ) , x 3 = x δ a 3 · ( d δ )
where a is random value in the gap [−2a, 2a].
The position of X β , X δ . and X α . is updated using the following equation:
x 1 + x 2 + x 3 3
Details and description of the HGWOSCA approach can be found in reference [43]. Furthermore, the computational procedure of the HGWOSCA approach is illustrated in Figure 8.
Although the GWO and SCA are able to expose an efficient accuracy in comparison with other well-known swarm intelligence optimisation techniques, it is not fitting for highly complex functions and may still face the difficulty of getting trapped in local optima [43]. Thus, a new hybrid variant based on GWO and SCA is used to solve recent real-life problems.

5. Results and Discussion

The proposed methodology has been implemented in MATLAB software and applied to the development of the optimal design of a PV plant connected to the electric grid. Solar irradiance, ambient temperature, and wind speed data for 1 year from the installation field are required. The effect of minimum LCOE and maximum annual energy objective functions on the PV plant design was determined. The HGWOSCA optimisation technique and a single SCA algorithm were applied with 400 search agents and 30 iterations to solve the design problem.
According to the results presented in Table 7, the PV plant optimal design variables depend on the selected objective function. The minimum LCOE and maximum annual energy result in two completely different optimal PV plant structures. PV power plant results are presented in Table 8. The optimisation process applying HGWOSCA outperforms the single SCA for minimum and maximum objective functions.
For both objectives using HGWOSCA, the optimisation process has selected mono-crystalline PV module type 3 (PV3) from the list of candidates. This module uses 295 W, and inverter type 3 (INV3) was selected from a list of three inverters. This inverter uses 500 kW and presents the central topology of the PV power plant. With the objective of maximum annual energy, the suggested number of PV modules is 1394 and 2 inverter to have 786.5035 (MWh). In this case, the LCOE was 32.1174 ($/MWh), and the PV plant total cost was the highest at 0.6315 (M$). With the objective of minimum LCOE, the number of PV modules is 1376 and only 1 inverter is required to have 28.6283 ($/MWh) of LCOE. In this case, the annual energy generation is equal to 786.5035 (MWh), and the total cost was reduced to 0.5557 (M$) compared with the first case. The use of LCOE’s objective function to optimise the design of PV plants can reduce the financial risks, as proven in this case study. The total cost of using minimum LCOE decreased by 12% with a benefit of 71,800 ($) in terms of installation cost, maintenance and operation costs. Figure 9 illustrates the maintenance and operational costs and the installation cost throughout the life of the PV plant for minimum LCOE and maximum annual energy generation.
The area occupied by the PV power plant can be calculated based on the summation of the occupied area by all PV rows, according to the length of each row and the inter-row area of all adjacent rows. The total available area of the installation field is equal to 3131 m2 and the installed PV modules occupied 3094 m2 of the installation site, which is nearly the same as the total area of the field. Therefore, the percentage of the occupied area by PV modules in the two cases presents 99% of the available area. The arrangement of PV modules in rows within the installation area is illustrated in Figure 10 using the LCOE objective function. The length of each row changed from one row to another according to the shape of the PV plant. Furthermore, this configuration has been designed in terms of the shape of the installation area, reflecting the actual situation. The difference obtained on the energy production using LCOE and maximum energy objective functions is due to the configuration and the arrangement of the PV modules within the available installation area. On the one hand, the optimal design of the PV plant under the maximum annual energy resulted in the minimum number of lines N r installed in each row, which is equal to 1. Additionally, this arrangement allowed the PV modules to capture more reflected radiation from the ground. Furthermore, at N r = 1 , a small distance between two adjacent rows in terms of shading effect is required, thereby increasing the total number of rows in the installation area to D r o w = 36 with one PV module line in each row and increasing the reflected radiation on PV modules. Moreover, the total number of PV modules for maximum energy is equal to 1394 and distributed among two central inverters. However, the number of PV modules for LCOE is less, leads to 1379 and arranged among only one inverter. PV modules are installed in multiple lines in case of LCOE objective function. In this configuration, the number of lines N r for each row is equal to 4 and leads only to 14 rows. Moreover, this configuration decreases the reflected radiation from the ground to be captured by PV modules and cannot be absorbed by the rest of the lines ( N r > 1 ) . The PV modules are installed horizontally for minimum LCOE ( P V o r i e n = 1 ) and vertically ( P V o r i e n = 2 ) for maximum annual energy.
Figure 11 illustrates the monthly energy generation by the PV power plant for the LCOE objective function. The PV plant energy generation remained high over the year, with an energy average of 65 (MWh) per month. The highest value of the energy generated by the PV power plant is obtained in March, because this condition is due to the high solar irradiance in this period.
For comparison, the semi-hourly average time was compared with the hourly average time meteorological data to examine the step time effect on the PV plant performance. The peaks of the meteorological data can influence the design solution. Therefore, the usage of annual semi-hourly average time rather than monthly, daily and hourly is recommended, as semi-hourly data contain the troughs and peaks of solar irradiation, ambient temperature, and wind speed. According to the results presented in Table 9 and Table 10, the step time data can affect the objective functions. The LCOE for semi-hourly average time is 28.6283 ($/MWh), and that obtained for hourly average time is higher and equal to 28.637 ($/MWh). The use of semi-hourly average time meteorological data in designing the PV plant can increase the financial benefits.
In all resulting cases, the proposed HGWOSCA optimisation approach was applied successfully and showed higher efficiency than that of a single SCA technique, with high performance in determining the optimal solution and solving the PV plant complex design problem. The convergence optimisation of annual energy and LCOE is illustrated in Figure 12 and Figure 13.

Effect of PV Module Reduction Coefficient

A sensitivity analysis was applied to evaluate the PV power plant performance. Accordingly, the variations in the PV module annual reduction coefficient were investigated. The optimisation results were obtained for different annual reduction coefficient values, from 0.3% to 0.7% per year. The annual reduction coefficient used in this study was 0.5%, as mentioned in Equation (44).
The optimum results for five different values for the annual reduction coefficient of the PV module are presented in Figure 14 and Figure 15. According to the results, by increasing the PV module reduction coefficient, the PV plant energy production is reduced throughout its lifetime period. The LCOE of the PV plant increases by increasing the PV module reduction coefficient. By contrast, the total cost of the PV power plant is not affected and has the same value for all reduction coefficient values.
In economic terms, an improved PV module annual reduction coefficient leads to the recovery of capital investment of the PV plant within a smaller time period, making the PV plant economically profitable. Moreover, the sensitivity of the PV power plant improved by the decrement of the PV module annual reduction coefficient and vice versa.

6. Conclusions

The proposed methodology was executed using semi-hourly time-resolution (i.e., 30 min-average) values of meteorological input data, including solar irradiance, ambient temperature, and wind speed. The procedure considers PV modules and inverter specifications, including a list of different commercially available PV modules and inverter technologies as candidates. The optimisation process selects only one PV module and inverter from a list of several alternatives, presenting the optimum combination. The proposed PV plant area model considers the shape and size of the installation field to properly arrange all the existing components.
The minimum LCOE and maximum annual energy objective functions were used to design the PV power plant. On the basis of the optimal results, the total cost of using the minimum LCOE objective function decreased by 12% with a benefit of 71,800 ($), including installation cost and maintenance and operation costs compared with the maximum annual energy. In this methodology, the HGWOSCA optimisation technique and a single SCA algorithm were applied. The optimum design solution shows that the proposed HGWOSCA is more efficient. Additionally, the PV plant optimal design variables depend on the selected objective function. The minimum LCOE and maximum annual energy result in two different optimal PV plant structures. LCOE improved with the use of semi-hourly average time meteorological data for designing the PV plant and can increase the financial benefits. Moreover, the sensitivity analysis shows that the PV power plant can be improved by the decrement of the PV module annual reduction coefficient and makes the PV plant economically more profitable.

Author Contributions

T.E.K.Z. contributed theoretical approaches, simulation, and preparing the article; M.R.A., M.F.N.T., S.M.Z. and A.D. contributed to supervision; M.R.A., A.D. and S.M. contributed to article editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the School of Electrical System Engineering Research Fund (SESERF), UniMAP.

Acknowledgments

We gratefully acknowledge the support of the Algerian company of electricity SONELGAZ, for providing the measurement data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohammadi, K.; Naderi, M.; Saghafifar, M. Economic feasibility of developing grid-connected photovoltaic plants in the southern coast of Iran. Energy 2018, 156, 17–31. [Google Scholar] [CrossRef]
  2. Zou, H.; Du, H.; Brown, M.A.; Mao, G. Large-scale PV power generation in China: A grid parity and techno-economic analysis. Energy 2017, 134, 256–268. [Google Scholar] [CrossRef]
  3. Feldman, D.; Margolis, R. Q2/Q3 2018 Solar Industry Update 2018; Technical Report; National Renewable Energy Lab. (NREL): Golden, CO, USA, 2018.
  4. Comello, S.; Reichelstein, S.; Sahoo, A. The road ahead for solar PV power. Renew. Sustain. Energy Rev. 2018, 92, 744–756. [Google Scholar] [CrossRef]
  5. Mondol, J.D.; Yohanis, Y.G.; Norton, B. Optimal sizing of array and inverter for grid-connected photovoltaic systems. Sol. Energy 2006, 80, 1517–1539. [Google Scholar] [CrossRef]
  6. Kornelakis, A.; Marinakis, Y. Contribution for optimal sizing of grid-connected PV-systems using PSO. Renew. Energy 2010, 35, 1333–1341. [Google Scholar] [CrossRef]
  7. Notton, G.; Lazarov, V.; Stoyanov, L. Optimal sizing of a grid-connected PV system for various PV module technologies and inclinations, inverter efficiency characteristics and locations. Renew. Energy 2010, 35, 541–554. [Google Scholar] [CrossRef]
  8. Sulaiman, S.I.; Rahman, T.K.A.; Musirin, I.; Shaari, S.; Sopian, K. An intelligent method for sizing optimization in grid-connected photovoltaic system. Sol. Energy 2012, 86, 2067–2082. [Google Scholar] [CrossRef]
  9. Chen, S.; Li, P.; Brady, D.; Lehman, B. Determining the optimum grid-connected photovoltaic inverter size. Sol. Energy 2013, 87, 96–116. [Google Scholar] [CrossRef]
  10. Perez-Gallardo, J.R.; Azzaro-Pantel, C.; Astier, S.; Domenech, S.; Aguilar-Lasserre, A. Ecodesign of photovoltaic grid-connected systems. Renew. Energy 2014, 64, 82–97. [Google Scholar] [CrossRef] [Green Version]
  11. Senol, M.; Abbasoğlu, S.; Kukrer, O.; Babatunde, A. A guide in installing large-scale PV power plant for self consumption mechanism. Sol. Energy 2016, 132, 518–537. [Google Scholar] [CrossRef]
  12. Báez-Fernández, H.; Ramirez-Beltran, N.D.; Méndez-Piñero, M.I. Selection and configuration of inverters and modules for a photovoltaic system to minimize costs. Renew. Sustain. Energy Rev. 2016, 58, 16–22. [Google Scholar] [CrossRef]
  13. Topić, D.; Knežević, G.; Fekete, K. The mathematical model for finding an optimal PV system configuration for the given installation area providing a maximal lifetime profit. Sol. Energy 2017, 144, 750–757. [Google Scholar] [CrossRef]
  14. Durusu, A.; Erduman, A. An improved methodology to design large-scale photovoltaic power plant. J. Sol. Energy Eng. 2017, 140, 011007. [Google Scholar] [CrossRef]
  15. Zidane, T.E.K.; Adzman, M.R.; Zali, S.M.; Mekhilef, S.; Durusu, A.; Tajuddin, M.F.N. Cost-effective topology for photovoltaic power plants using optimization design. In Proceedings of the 2019 IEEE 7th Conference on Systems, Process and Control (ICSPC), Malacca, Malaysia, 13–14 December 2019; pp. 210–215. [Google Scholar]
  16. Wang, H.; Muñoz-García, M.; Moreda, G.; Alonso-Garcia, C. Optimum inverter sizing of grid-connected photovoltaic systems based on energetic and economic considerations. Renew. Energy 2018, 118, 709–717. [Google Scholar] [CrossRef]
  17. Perez-Gallardo, J.R.; Azzaro-Pantel, C.; Astier, S. Combining multi-objective optimization, principal component analysis and multiple criteria decision making for ecodesign of photovoltaic grid-connected systems. Sustain. Energy Technol. Assess. 2018, 27, 94–101. [Google Scholar] [CrossRef] [Green Version]
  18. Zidane, T.E.K.; Adzman, M.R.; Tajuddin, M.F.N.; Zali, S.M.; Durusu, A. Optimal configuration of photovoltaic power plant using grey wolf optimizer: A comparative analysis considering CdTe and c-Si PV modules. Sol. Energy 2019, 188, 247–257. [Google Scholar] [CrossRef]
  19. Bakhshi-Jafarabadi, R.; Sadeh, J.; Soheili, A. Global optimum economic designing of grid-connected photovoltaic systems with multiple inverters using binary linear programming. Sol. Energy 2019, 183, 842–850. [Google Scholar] [CrossRef]
  20. Sulaiman, S.I.; Rahman, T.K.A.; Musirin, I.; Shaari, S. Sizing grid-connected photovoltaic system using genetic algorithm. In Proceedings of the 2011 IEEE Symposium on Industrial Electronics and Applications, Langkawi, Malaysia, 25–28 September 2011; pp. 505–509. [Google Scholar] [CrossRef]
  21. Al-Sabounchi, A.M.; Yalyali, S.A.; Al-Thani, H.A. Design and performance evaluation of a photovoltaic grid-connected system in hot weather conditions. Renew. Energy 2013, 53, 71–78. [Google Scholar] [CrossRef]
  22. Kornelakis, A.; Koutroulis, E. Methodology for the design optimisation and the economic analysis of grid-connected photovoltaic systems. IET Renew. Power Gener. 2009, 3, 476. [Google Scholar] [CrossRef]
  23. Koutroulis, E.; Yang, Y.; Blaabjerg, F. Co-Design of the PV array and DC/AC inverter for maximizing the energy production in grid-connected applications. IEEE Trans. Energy Convers. 2018, 34, 509–519. [Google Scholar] [CrossRef] [Green Version]
  24. Moradi-Shahrbabak, Z.; Tabesh, A.; Yousefi, G.R. Economical design of utility-scale photovoltaic power plants with optimum availability. IEEE Trans. Ind. Electron. 2013, 61, 3399–3406. [Google Scholar] [CrossRef]
  25. Okido, S.; Takeda, A. Economic and environmental analysis of photovoltaic energy systems via robust optimization. Energy Syst. 2013, 4, 239–266. [Google Scholar] [CrossRef]
  26. Paravalos, C.; Koutroulis, E.; Samoladas, V.; Kerekes, T.; Sera, D.; Teodorescu, R. Optimal design of photovoltaic systems using high time-resolution meteorological data. IEEE Trans. Ind. Inform. 2014, 10, 2270–2279. [Google Scholar] [CrossRef]
  27. Ramli, M.; Hiendro, A.; Sedraoui, K.; Twaha, S. Optimal sizing of grid-connected photovoltaic energy system in Saudi Arabia. Renew. Energy 2015, 75, 489–495. [Google Scholar] [CrossRef]
  28. Rodrigo, P.M. Improving the profitability of grid-connected photovoltaic systems by sizing optimization. In Proceedings of the 2017 IEEE Mexican Humanitarian Technology Conference (MHTC), Puebla, Mexico, 29–31 March 2017; pp. 1–6. [Google Scholar]
  29. Arefifar, S.A.; Paz, F.; Ordonez, M. Improving solar power PV plants using multivariate design optimization. IEEE J. Emerg. Sel. Top. Power Electron. 2017, 5, 638–650. [Google Scholar] [CrossRef]
  30. Camps, X.; Velasco, G.; De La Hoz, J.; Martín, H. Contribution to the PV-to-inverter sizing ratio determination using a custom flexible experimental setup. Appl. Energy 2015, 149, 35–45. [Google Scholar] [CrossRef]
  31. Demoulias, C. A new simple analytical method for calculating the optimum inverter size in grid-connected PV plants. Electr. Power Syst. Res. 2010, 80, 1197–1204. [Google Scholar] [CrossRef]
  32. Fernández-Infantes, A.; Contreras, J.; Bernal-Agustín, J.L. Design of grid connected PV systems considering electrical, economical and environmental aspects: A practical case. Renew. Energy 2006, 31, 2042–2062. [Google Scholar] [CrossRef]
  33. Gemignani, M.; Rostegui, G.; Salles, M.B.C.; Kagan, N. Methodology for the sizing of a solar PV plant: A comparative case. In Proceedings of the 2017 6th International Conference on Clean Electrical Power (ICCEP), Santa Margherita Ligure, Italy, 27–29 June 2017; pp. 7–13. [Google Scholar]
  34. Kerekes, T.; Koutroulis, E.; Sera, D.; Teodorescu, R.; Katsanevakis, M. An optimization method for designing large PV Plants. IEEE J. Photovoltaics 2012, 3, 814–822. [Google Scholar] [CrossRef]
  35. Kerekes, T.; Koutroulis, E.; Eyigun, S.; Teodorescu, R.; Katsanevakis, M.; Sera, D. A practical optimization method for designing large PV plants. In Proceedings of the 2011 IEEE International Symposium on Industrial Electronics, Gdansk, Poland, 27–30 June 2011; pp. 2051–2056. [Google Scholar] [CrossRef]
  36. Kornelakis, A. Multiobjective particle swarm optimization for the optimal design of photovoltaic grid-connected systems. Sol. Energy 2010, 84, 2022–2033. [Google Scholar] [CrossRef] [Green Version]
  37. Alsadi, S.; Khatib, T. Photovoltaic power systems optimization research status: A review of criteria, constrains, models, techniques, and software tools. Appl. Sci. 2018, 8, 1761. [Google Scholar] [CrossRef] [Green Version]
  38. Baños, R.; Manzano-Agugliaro, F.; Montoya, F.; Gil, C.; Alcayde, A.; Gómez, J. Optimization methods applied to renewable and sustainable energy: A review. Renew. Sustain. Energy Rev. 2011, 15, 1753–1766. [Google Scholar] [CrossRef]
  39. Eltawil, M.A.; Zhao, Z. Grid-connected photovoltaic power systems: Technical and potential problems—A review. Renew. Sustain. Energy Rev. 2010, 14, 112–129. [Google Scholar] [CrossRef]
  40. Khatib, T.; Mohamed, A.; Sopian, K. A review of photovoltaic systems size optimization techniques. Renew. Sustain. Energy Rev. 2013, 22, 454–465. [Google Scholar] [CrossRef]
  41. Mellit, A.; Kalogirou, S.; Hontoria, L.; Shaari, S. Artificial intelligence techniques for sizing photovoltaic systems: A review. Renew. Sustain. Energy Rev. 2009, 13, 406–419. [Google Scholar] [CrossRef]
  42. Maatallah, T.; El Alimi, S.; Ben Nassrallah, S. Performance modeling and investigation of fixed, single and dual-axis tracking photovoltaic panel in Monastir city, Tunisia. Renew. Sustain. Energy Rev. 2011, 15, 4053–4066. [Google Scholar] [CrossRef]
  43. Singh, N.; Singh, S. A novel hybrid GWO-SCA approach for optimization problems. Eng. Sci. Technol. Int. J. 2017, 20, 1586–1601. [Google Scholar] [CrossRef]
  44. Mirjalili, S. SCA: A sine cosine algorithm for solving optimization problems. Knowl.-Based Syst. 2016, 96, 120–133. [Google Scholar] [CrossRef]
  45. Energy Information Administration (EIA) USA. Total Electricity Net Generation. 2014. Available online: https://www.eia.gov/beta/international/data/browser (accessed on 30 March 2015).
  46. Fodhil, F.; Hamidat, A.; Nadjemi, O. Potential, optimization and sensitivity analysis of photovoltaic-diesel-battery hybrid energy system for rural electrification in Algeria. Energy 2019, 169, 613–624. [Google Scholar] [CrossRef]
  47. Saiah, S.B.D.; Stambouli, A.B. Prospective analysis for a long-term optimal energy mix planning in Algeria: Towards high electricity generation security in 2062. Renew. Sustain. Energy Rev. 2017, 73, 26–43. [Google Scholar] [CrossRef]
  48. Alegria CDER. 2010. Available online: http://www.cder.dz (accessed on 12 April 2010).
  49. Abada, Z.; Bouharkat, M. Study of management strategy of energy resources in Algeria. Energy Rep. 2018, 4, 1–7. [Google Scholar] [CrossRef]
  50. Solar Resource Maps and GIS Data for 200+ Countries. Available online: https://solargis.com/maps-and-gis-data/download/algerian.d. (accessed on 7 June 2017).
  51. Stambouli, A.B. Promotion of renewable energies in Algeria: Strategies and perspectives. Renew. Sustain. Energy Rev. 2011, 15, 1169–1181. [Google Scholar] [CrossRef]
  52. World Bank. World Development Indicators. 2017. Available online: http://data.worldbank.org/data-catalog/world-development-indicators (accessed on 1 March 2017).
  53. Bouznit, M.; Pablo-Romero, M.D.P. CO2 emission and economic growth in Algeria. Energy Policy 2016, 96, 93–104. [Google Scholar] [CrossRef]
  54. Boukhelkhal, A.; Bengana, I. Cointegration and causality among electricity consumption, economic, climatic and environmental factors: Evidence from North-Africa region. Energy 2018, 163, 1193–1206. [Google Scholar] [CrossRef]
  55. Laib, I.; Hamidat, A.; Haddadi, M.; Ramzan, N.; Olabi, A. Study and simulation of the energy performances of a grid-connected PV system supplying a residential house in north of Algeria. Energy 2018, 152, 445–454. [Google Scholar] [CrossRef] [Green Version]
  56. Mihoub, S.; Chermiti, A.; Beltagy, H. Methodology of determining the optimum performances of future concentrating solar thermal power plants in Algeria. Energy 2017, 122, 801–810. [Google Scholar] [CrossRef]
  57. Ministry of Energy Algeria. Algerian Energy Review N4. May 2015; ISSN 2437-0479. Available online: www.mem-algeria.org (accessed on 15 July 2015).
  58. Ministry of Energy Algeria. Objectives New Renewable Energies Program in Algeria (2015-2020-2030). 2015. Available online: www.mem-algeria.org (accessed on 5 September 2015).
  59. The Weekend Read: Sustainable Finance—Reassessing Risk, Purpose. Available online: https://www.pv-magazine.com/2019/08/16/algerian-diesel-solar-tender-concludes-with-lowest-bid-of-11-cents/n.d. (accessed on 16 August 2019).
  60. Fadaee, M.; Radzi, M.A.M. Multi-objective optimization of a stand-alone hybrid renewable energy system by using evolutionary algorithms: A review. Renew. Sustain. Energy Rev. 2012, 16, 3364–3369. [Google Scholar] [CrossRef]
  61. Panchula, A.F.; Hayes, W.; Bilash, J.; Kimber, A.; Graham, C. First year performance of a 20 MWac PV power plant. In Proceedings of the 2011 37th IEEE Photovoltaic Specialists Conference, Washington, DC, USA, 19–24 June 2011; p. 001993. [Google Scholar]
  62. Krishnamoorthy, H.S.; Essakiappan, S.; Enjeti, P.N.; Balog, R.; Ahmed, S. A new multilevel converter for Megawatt scale solar photovoltaic utility integration. In Proceedings of the 2012 Twenty-Seventh Annual IEEE Applied Power Electronics Conference and Exposition (APEC), Orlando, FL, USA, 5–9 February 2012; pp. 1431–1438. [Google Scholar]
  63. Marcos, J.; De La Parra, I.; García, M.; Marroyo, L.; Cirés, E. The potential of forecasting in reducing the LCOE in PV plants under ramp-rate restrictions. Energy 2019, 188, 116053. [Google Scholar] [CrossRef]
  64. Sharma, V.; Kumar, A.; Sastry, O.S.; Chandel, S. Performance assessment of different solar photovoltaic technologies under similar outdoor conditions. Energy 2013, 58, 511–518. [Google Scholar] [CrossRef]
  65. Cabrera-Tobar, A.; Bullich-Massagué, E.; Aragüés-Peñalba, M.; Gomis-Bellmunt, O. Topologies for large scale photovoltaic power plants. Renew. Sustain. Energy Rev. 2016, 59, 309–319. [Google Scholar] [CrossRef] [Green Version]
  66. Testa, A.; De Caro, S.; La Torre, R.; Scimone, T. A probabilistic approach to size step-up transformers for grid connected PV plants. Renew. Energy 2012, 48, 42–51. [Google Scholar] [CrossRef]
  67. Noorian, A.M.; Moradi, I.; Kamali, G.A. Evaluation of 12 models to estimate hourly diffuse irradiation on inclined surfaces. Renew. Energy 2008, 33, 1406–1412. [Google Scholar] [CrossRef]
  68. Demain, C.; Journée, M.; Bertrand, C. Evaluation of different models to estimate the global solar radiation on inclined surfaces. Renew. Energy 2013, 50, 710–721. [Google Scholar] [CrossRef]
  69. Liu, B.Y.H.; Jordan, R.C. A rational procedure for predicting the long-term average performance of flat-plate solar-energy collectors with design data for the U. S., its outlying possessions and Canada. Sol. Energy 1963, 7, 53–74. [Google Scholar] [CrossRef]
  70. Ayaz, R.; Durusu, A.; Akca, H. Determination of optimum tilt angle for different photovoltaic technologies considering ambient conditions: A case study for Burdur, Turkey. J. Sol. Energy Eng. 2017, 139, 041001. [Google Scholar] [CrossRef]
  71. Chin, V.J.; Chin, V.J.; Ishaque, K. Cell modelling and model parameters estimation techniques for photovoltaic simulator application: A review. Appl. Energy 2015, 154, 500–519. [Google Scholar] [CrossRef]
  72. Orioli, A.; Di Gangi, A. A procedure to evaluate the seven parameters of the two-diode model for photovoltaic modules. Renew. Energy 2019, 139, 582–599. [Google Scholar] [CrossRef]
  73. Elbaset, A.A.; Ali, H.; Sattar, M.A.-E.; Elbaset, A.A. Novel seven-parameter model for photovoltaic modules. Sol. Energy Mater. Sol. Cells 2014, 130, 442–455. [Google Scholar] [CrossRef]
  74. Dhoke, A.; Sharma, R.; Saha, T.K. PV module degradation analysis and impact on settings of overcurrent protection devices. Sol. Energy 2018, 160, 360–367. [Google Scholar] [CrossRef]
  75. Rosillo, F.G.; Alonso-Garcia, C. Evaluation of color changes in PV modules using reflectance measurements. Sol. Energy 2019, 177, 531–537. [Google Scholar] [CrossRef]
  76. Deline, C.; Deline, C. Partially shaded operation of a grid-tied PV system. In Proceedings of the 2009 34th IEEE Photovoltaic Specialists Conference (PVSC), Philadelphia, PA, USA, 7–12 June 2009; pp. 001268–001273. [Google Scholar] [CrossRef] [Green Version]
  77. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Horizontal Irradiation of Algeria [50].
Figure 1. Horizontal Irradiation of Algeria [50].
Energies 13 02776 g001
Figure 2. Irradiance data. (a) Hourly average time. (b) 30 min average time.
Figure 2. Irradiance data. (a) Hourly average time. (b) 30 min average time.
Energies 13 02776 g002
Figure 3. Temperature data. (a) Hourly average time. (b) 30 min average time.
Figure 3. Temperature data. (a) Hourly average time. (b) 30 min average time.
Energies 13 02776 g003
Figure 4. Hourly solar irradiance over the year (W/m2).
Figure 4. Hourly solar irradiance over the year (W/m2).
Energies 13 02776 g004
Figure 5. Solar irradiance over the year (W/m2).
Figure 5. Solar irradiance over the year (W/m2).
Energies 13 02776 g005
Figure 6. Flowchart of the proposed design methodology.
Figure 6. Flowchart of the proposed design methodology.
Energies 13 02776 g006
Figure 7. Arrangement of PV modules according to PV plant shape.
Figure 7. Arrangement of PV modules according to PV plant shape.
Energies 13 02776 g007
Figure 8. Pseudo-code of the hybrid grey wolf optimiser-sine cosine algorithm (HGWOSCA).
Figure 8. Pseudo-code of the hybrid grey wolf optimiser-sine cosine algorithm (HGWOSCA).
Energies 13 02776 g008
Figure 9. Throughout the life of the PV plant optimised by HGWOSCA.
Figure 9. Throughout the life of the PV plant optimised by HGWOSCA.
Energies 13 02776 g009
Figure 10. Rows arrangement for minimum LCOE using HGWOSCA.
Figure 10. Rows arrangement for minimum LCOE using HGWOSCA.
Energies 13 02776 g010
Figure 11. PV plant energy generation (MWh).
Figure 11. PV plant energy generation (MWh).
Energies 13 02776 g011
Figure 12. Convergence of the optimisation of annual energy using HGWOSCA algorithm for semi-hourly data.
Figure 12. Convergence of the optimisation of annual energy using HGWOSCA algorithm for semi-hourly data.
Energies 13 02776 g012
Figure 13. The convergence of the optimisation of LCOE using HGWOSCA algorithm for semi-hourly data.
Figure 13. The convergence of the optimisation of LCOE using HGWOSCA algorithm for semi-hourly data.
Energies 13 02776 g013
Figure 14. Total energy for reduction coefficient variations.
Figure 14. Total energy for reduction coefficient variations.
Energies 13 02776 g014
Figure 15. LCOE for reduction coefficient variations.
Figure 15. LCOE for reduction coefficient variations.
Energies 13 02776 g015
Table 1. Radiation in Algeria [49].
Table 1. Radiation in Algeria [49].
RegionCoastalHighlandsSahara
Area (%)41086
Average sun hours per year265030003500
Energy received kWh/m2/year170019002650
Table 2. Phases of the Algerian renewable energy program [58].
Table 2. Phases of the Algerian renewable energy program [58].
Energy Type1st Phase 2015–2020 (MW)2nd Phase 2021–2030 (MW)Total (MW)
Photovoltaic300010,57513,575
Wind101040005010
CSP-20002000
Cogeneration150250400
Biomass3606401000
Geothermal51015
Total452517,47522,000
Table 3. Projects of Photovoltaic power plant installation in Algeria.
Table 3. Projects of Photovoltaic power plant installation in Algeria.
NPV Power PlantProvincePower (MWp)
1In GuezzemTamanrasset6
2TinzaouatineTamanrasset3
3DjanetIllizi 4
4Bordj Omar DrisIllizi3
5Bordj Badji MokhtarAdrar10
6TimiaouineAdrar2
7TalmineAdrar8
8TabelbalaBechar3
9TindoufTindouf11
Table 4. Photovoltaic power plants installed in Algeria.
Table 4. Photovoltaic power plants installed in Algeria.
NPV Power PlantProvincePower (MWp)Energy Production (GWh)-June 2017Commissining DateArea (km2)Distance to Electric Grid (km)Voltage Level at Point of Coupling (kV)PV Modules Type Topology
1Adrar Adrar2059.585Oct-150.42.830PolyC
2Kabertène Adrar39.584Oct-150.060.230PolyC
3In Salah Tamanrasset 512.328Feb-160.10.530PolyC
4Timimoune Adrar923.822Feb-160.18930PolyC
5Regguen Adrar512.221Jan-160.10.2230PolyC
6Zaouiat Kounta Adrar615.213Jan-160.120.2430PolyC
7AoulefAdrar512.557Mar-160.10.430PolyC
8Tamanrasset Tamanrasset 1336.41Nov-150.268.830PolyC
9Djanet Ilizi310.729Feb-150.060.3530PolyC
10Tindouf Tindouf96.376Dec-150.180.330PolyC
11Oued NechouGhardaia1.14.593Jul-140.10.430m-Si,c-Si,Thin film, si-amorpheC
12Sedret LeghzelNaâma2040.751May-160.421.360PolyC
13Oued El KebritSouk Ahras1528.9Apr-160.326.530PolyC
14Ain SkhounaSaida3014.213May-160.421260PolyC
15Ain El Bel 1Djelfa2025.134Apr-160.43.960PolyC
16Ain El Bel 2Djelfa330.860PolyC
17Lekhneg 1Laghouat2053.576Apr-160.41260PolyC
18Lekhneg 2Laghouat400.860PolyC
19TelaghSidi-Bel-Abbes127.417Sep-160.3660PolyC
20Labiodh Sidi ChikhEl-Bayadh2319.146Nov-180.3920.960PolyC
21El HdjiraOuargla309.738Jul-160.6160PolyC
22Ain-El-MelhM’Sila2016.473Sep-170.40.960PolyC
23Oued El MaBatna2-Mar-170.30.530PolyC
Table 5. Photovoltaic (PV) modules specifications at standard test condition.
Table 5. Photovoltaic (PV) modules specifications at standard test condition.
Specification at STCUnitPV1PV2PV3
Nominal maximum power ( P m p p , s t c )W280285295
Optimum operating current ( I m p p , s t c )A8.849.029.08
Optimum operating voltage ( V m p p , s t c )V31.731.632.5
Current temperature coefficient ( K i )(%/C)0.050.050.05
Voltage temperature coefficient ( K v )(%/C)−0.29−0.32−0.29
Open circuit voltage ( V o c , s t c )V38.438.339.6
Wind speed temperature coefficient ( K r )-1.5091.46841.509
Length ( L p v , 1 )m1.651.651.65
Width ( L p v , 2 )m0.9920.9920.992
Efficiency%17.117.418
Type-MonoPolyMono
Table 6. Inverters specifications at standard test condition.
Table 6. Inverters specifications at standard test condition.
SpecificationUnitINV1INV2INV3
Nominal power ( P i )kW507500
Minimum input voltage ( V i , m i n )V250335450
Maximum input voltage ( V i , m a x )V950560880
Maximum power point tracking voltage ( V i , m p p t , m a x )V850540820
Power loss ( P i , s c )W1.51200
Efficiency ( η i n v )%97.595.398.7
Table 7. Optimal design variables using semi-hourly measurement data.
Table 7. Optimal design variables using semi-hourly measurement data.
Design VariablesMinimum Levelised Cost of Energy (LCOE)Maximum Energy
Sine Cosine Algorithm (SCA)HGWOSCASCAHGWOSCA
N s 20182116
N p 70744962
N r 5411
β 15151515
F y 1.9251.5420.6400.643
P V i PV1PV3PV3PV3
I N V i INV3INV3INV3INV3
P V o r i e n 1122
Table 8. Results of optimal design algorithms using semi-hourly measurement data.
Table 8. Results of optimal design algorithms using semi-hourly measurement data.
PV Plant ParametersMinimum LCOEMaximum Energy
SCAHGWOSCASCAHGWOSCA
LCOE ($/MWh)29.182928.628332.098332.1174
Yearly total energy (MWh)731.4199776.4012785.8698786.5035
Total energy (GWh)18.285519.41019.646719.6626
Total cost (M$)0.53360.55570.63050.6315
Installation Cost (M$)0.44650.46320.53680.5378
Maintenance cost (M$)0.08710.09250.09360.0937
Junction Boxes671818
PV modules ( N I )1365137613931394
Inverters ( N i )1122
Rows11143636
Table 9. Optimal design variables using hourly measurement data.
Table 9. Optimal design variables using hourly measurement data.
Design VariablesMinimum LCOEMaximum Energy
SCAHGWOSCASCAHGWOSCA
N s 19182214
N p 70744247
N r 5411
β 15151515
F y 1.9251.5400.640.64
P V i PV3PV3PV3PV3
I N V i INV3INV3INV3INV3
P V o r i e n 1122
Table 10. Results of optimal design algorithms using hourly measurement data.
Table 10. Results of optimal design algorithms using hourly measurement data.
PV Plant ParametersMinimum LCOEMaximum Energy
SCAHGWOSCASCAHGWOSCA
LCOE ($/MWh)28.662228.637032.098332.1552
Yearly total energy (MWh)770.0032776.1651785.6774786.2915
Total energy (GWh)19.250119.404119.641919.6573
Total cost (M$)0.55170.55570.63050.6321
Installation Cost (M$)0.46000.46320.53680.5384
Maintenance cost (M$)0.09170.09250.09360.0937
Junction Boxes671818
PV modules ( N I )1365137613931394
Inverters ( N i )1122
Rows11143636

Share and Cite

MDPI and ACS Style

Zidane, T.E.K.; Adzman, M.R.; Tajuddin, M.F.N.; Zali, S.M.; Durusu, A.; Mekhilef, S. Optimal Design of Photovoltaic Power Plant Using Hybrid Optimisation: A Case of South Algeria. Energies 2020, 13, 2776. https://doi.org/10.3390/en13112776

AMA Style

Zidane TEK, Adzman MR, Tajuddin MFN, Zali SM, Durusu A, Mekhilef S. Optimal Design of Photovoltaic Power Plant Using Hybrid Optimisation: A Case of South Algeria. Energies. 2020; 13(11):2776. https://doi.org/10.3390/en13112776

Chicago/Turabian Style

Zidane, Tekai Eddine Khalil, Mohd Rafi Adzman, Mohammad Faridun Naim Tajuddin, Samila Mat Zali, Ali Durusu, and Saad Mekhilef. 2020. "Optimal Design of Photovoltaic Power Plant Using Hybrid Optimisation: A Case of South Algeria" Energies 13, no. 11: 2776. https://doi.org/10.3390/en13112776

APA Style

Zidane, T. E. K., Adzman, M. R., Tajuddin, M. F. N., Zali, S. M., Durusu, A., & Mekhilef, S. (2020). Optimal Design of Photovoltaic Power Plant Using Hybrid Optimisation: A Case of South Algeria. Energies, 13(11), 2776. https://doi.org/10.3390/en13112776

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

Article Metrics

Back to TopTop