A Multi-Criteria Decision Maker for Grid-Connected Hybrid Renewable Energy Systems Selection Using Multi-Objective Particle Swarm Optimization

Combating climate change issues resulting from excessive use of fossil fuels comes with huge initial costs, thereby posing difficult challenges for the least developed countries in Sub-Saharan Africa (SSA) to invest in renewable energy alternatives, especially with rapid industrialization. However, designing renewable energy systems usually hinges on different economic and environmental criteria. This paper used the Multi-Objective Particle Swarm Optimization (MOPSO) technique to optimally size ten grid-connected hybrid blocks selected amongst Photo-Voltaic (PV) panels, onshore wind turbines, biomass combustion plant using sugarcane bagasse, Battery Energy Storage System (BESS), and Diesel Generation (DG) system as backup power, to reduce the supply deficit in Sierra Leone. Resource assessment using well-known methods was done for PV, wind, and biomass for proposed plant sites in Kabala District in Northern and Kenema District in Southern Sierra Leone. Long term analysis was done for the ten hybrid blocks projected over 20 years whilst ensuring the following objectives: minimizing the Deficiency of Power Supply Probability (DPSP), Diesel Energy Fraction (DEF), Life Cycle Costs (LCC), and carbon dioxide (CO2) emissions. Capacity factors of 27.41% and 31.6% obtained for PV and wind, respectively, indicate that Kabala district is the most feasible location for PV and wind farm installations. The optimum results obtained are compared across selected blocks for DPSP values of 0–50% to determine the most economical and environmentally friendly alternative that policy makers in Sierra Leone and the region could apply to similar cases.


Introduction
Energy is an essential component in the development of any society.The understanding of this fundamental principle is deeply embedded in the United Nations Sustainable Development Goal 7 (UNSDG7) [1], which clearly emphasizes universal access to affordable, reliable, and modern energy services, whilst increasing substantially the share of renewable energy in the global energy mix by 2030.Whilst some countries are taking leading roles in meeting this agenda, most countries in Sub-Saharan Africa (SSA) have been snail paced.The region still has the lowest household electrification rate worldwide, with about 588 million people living without (2016 report) [2].Though Sierra Leone, a small country located along the Atlantic Ocean in West Africa, has a small population of about 7 million, the daily power delivery is not guaranteed and is currently far below the overall demand.The obligation of the Ministry of Energy (MoE) to provide sustainable electricity for its citizens over the past decades has been met with several obstacles including the civil war (1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001) leading to vandalizing existing energy infrastructure, crippling management systems, and minimum funds to invest in renewable energy technologies.This leaves the current access to electricity to be barely 15% with about 3% rural access, which is well below the current average in SSA.Political will and economic instability have been driving factors limiting governments from investing in state of the art technologies.The unemployment rate is high (4.5%)[3] for a country with a very small population and an electricity tariff of 18.76 USc/kWh, sitting amongst the highest in region, making electricity less affordable to the poor.These measures are in complete disagreement with the UNSDG7 agenda.Irrespective of other domestic and international factors, lack of electricity access has influenced to a great extent the snail paced growth of private sector investment.
At the moment, desperate but uneconomical and unsustainable measures are being adopted by governments in trying to meet the unsuppressed national demand of about 700 MW, inclusive of existing and projected demand of mining companies and new settlements across the country.However, these measures are taken with little concern for the environment as more rental power with diesel generation is been injected to close the gap between generation and demand.According to the Sustainable Energy for All (SE4ALL) initiative launched in 2012, the target trajectory for Sierra Leone requires increasing electricity access from 15% to 92% by 2030 [4].This would require considerably improving the share of renewable energy technologies to meet the demand in order to minimize the current burden on the government to heavily subsidize the current cost of electricity at a rate of 0.15 US$/kWh for private residents and 0.26 US$/kWh for industry.Albeit the potentials of Solar (1600 kWh/m 2 /year) or an estimated 2180 h of sunshine, hydro (about 2000 MW for 27 sites), and biomass, about 2000 GWh/year, the total installed government owned generation capacity is about 130 MW, and another operational 50 MW rental power from an Independent Power Producer (IPP).Subscribing to IPPs running on diesel engines usually comes with high tariffs and this accounts for very high operational costs owing to high and fluctuating fuel costs.For the government-owned diesel plants, a huge percentage of the total revenue collected is utilized for the purchase of imported heavy fuel oil and diesel oil to run the plants, especially during the dry season.This causes the government to heavily subsidize the electricity generation.It is against this backdrop that this research has been undertaken.
Many authors have done work on hybrid energy systems ranging from residential, institutional, industrial, off-grid to grid applications.Renewable Distributed Generations (RDGs) are essential in power systems as they provide energy security and reduce power losses whilst increasing the overall efficiency and environmental protection [5].However, renewable energy sources are stochastic in nature and thus the application of hybrid energy systems is essential in overcoming the weaknesses of stochasticity with the strength of the other, or predictability of a conventional technology [6].Literature [7] analyzed four different hybrid energy options to reduce overall energy consumption and CO 2 emissions from residential sector.The authors in [8] presented the possibility of integrating a grid-connected renewable energy system to meet the energy demand of institutional buildings for an Indian scenario.Additionally, the authors employed a modern approach based on the Analytic Hierarchy Process for the optimum planning of the electric power system.Establishing standard methodologies to integrate RDGs to overcome today's challenges towards sustainable development is crucial and very complicated.Literature [9] examined the use of multi-criteria decision making methods for ranking three subcidy schemes of grid-connected and autonomous solar Photo-Voltaic (PV) technologies, for promoting PV systems in Cyprus.An optimization based on techno-economic analysis of grid-connected hybrid systems was done also in India, with comparison of three different scenarios based on per unit cost of electricity production, cost of operating diesel generator, and minimization of Green House Gases (GHG) [10].The authors in [11] optimally sized stand-alone and grid-connected PV technology for a case in Australia, whilst comparing with the optimized model of Diesel Generation (DG) system.The literature used the basis of economic feasibility and environmental impact anlysis to compare the optimized models.Resource assessment and techno-economic analysis of Hybrid Energy systems (HES) has been done with a case study approach in Saudi Arabia [12].The authors did resource assessments for wind and PV technologies to determine their technical and economic feasibility for four cities in Saudi Arabia and achieved a reduction of energy costs and GHG emissions.Reference [13] carried out resource assessments using similar approaches in Nigeria, a case study in SSA.Literature [14] carried out an assessment of grid-connected PV and Biomass energy system to provide electricity for rural communities in Egypt.The results obtained revealed that the hybrid renewable energy system proves efficient for emissions and cost reductions.However, with the exception of assessments done by Japanese International Cooperation Agency (JICA), World Bank, and other donor organizations, no academic research has been carried out on both resource assessment and techno-economic analysis of grid-connected systems in Sierra Leone.Additionally, from the database of similar work done on multi-criteria analysis for grid-connected systems, no work has been done considering resource assessments and analysis of ten different hybrid considerations for five different technologies.Therefore the novelty of this research stems from this background with focus on the Sierra Leone scenario, whilst establishing practical boundary conditions on the resources based on the results obtained from the resource assessments.In most engineering problems the established multiple objectives are often conflicting which implies any improvement in one leads to a suppression of the other.Therefore establishing a trade-off between different objectives give rise to different solutions.Several researchers have used different tools for integration of renewable resources, like the Holistic Grid Resource Integration and Development (HiGRID) tool [15] and HOMER simulator [16].The literature [17] examines dynamic operation and control strategies for a microgrid hybrid energy system and applied Particle Swarm Optimization (PSO) algorithm to analyze the performance of the PV power system.The authors in literature [5] applied a hybrid PSO and Gravitational Search Algorithm (PSOGSA) to deduce the optimum location of PV and wind systems for minimizing system power losses and operating costs whilst improving voltage profile and stability.Genetic algorithm (GA) was used to optimally size a PV hybrid renewable energy system in [18] with the objective of decreasing operation costs and emission of pollutants.MOPSO has been used in this study to optimally minimize the DEF, DPSP, LCC, and CO 2 emissions.
In general, the PSO algorithm performs very well with higher scale or as the problem becomes complex.Using evolutionary algorithms, such as Non-dominated Sorting Genetic Algorithm (NSGA-II), are usually computationally expensive in that the crowding distance works in the objective space only [19].In addition to its speed of convergence and ease of application requiring the use of simple operation of real numbers, MOPSO does not require sorting of fitness values of solutions in any process.It uses leaders or region based selection instead of individuals within the search space to guide the search.Hence, this presents a significant computational advantage as the size gets larger [20].Ting LI et al. elaborated on the methodology of the standard MOPSO algorithm and presented a comprehensive survey whilst listing references on its applications in power systems economic dispatch including Economic Environmental Dispatch (EED) and multi-area EED [21].The authors further highlighted the advantages of MOPSO over other optimization methods.However, it is worth noting that MOPSO has computational limitations of handling constraints and the static analysis of tradeoffs.WL Theo et al. [22], also presented a review on system planning and optimization techniques applicable to hybrid energy systems integration whilst comparing various mathematical programming methods.Reference [23] used MOPSO algorithm to compute sufficient capacity of solar PV that the power system can accommodate while maintaining the voltage stability of the system.The literature [24] also applied the MOPSO technique to optimally size grid-connected ground mounted and rooftop PV system in Freetown, Sierra Leone, to reduce the daily generation cost and gap between generation and demand.
The importance of this research is thus three folds; it can help policy makers to understand the magnitude of Green House Gas (GHG) emissions obtained from thermal generation in Sierra Leone, though relatively small (under 0.03% of Global emissions) at the moment; it helps in providing both technical and economic basis for the adoption of suitable renewable energy technologies and the best hybrid options for integration; finally, it helps in establishing a holistic decision making process from resource availability to conversion and its environmental effects.

Methodology
An in-depth understanding of the extent of availability of renewable energy resources and the suitability in terms of the appropriate location to carryout the project is critical to establishing the financial viability of any of such projects.Fossil fuel technologies do not directly fall under this scope of assessment as the fuel resources are mostly located at or very close to the generating facilities.Policy makers thus need to be well informed about the quality, proximity of the resource to demand centers and take a closer look at the potential to reduce costs through the sharing of grid infrastructure and streamlined permitting [25].At the moment there are limited comprehensive approaches that could be adopted for estimating national renewable energy generation in SSA.To help in overcoming this drawback, estimates of solar, wind, and biomass for every country in SSA has been done [26].
Amidst the several approaches adopted and mentioned in literatures [27,28], this study does not only focus on the technical and economic components of the factors that influence the viability of the project, but also throws light on the social and environmental sustainability of the project.Therefore, a summary indicating the inter dependency among the key components considered has been illustrated in the pentagonal relationship as highlighted in the bullets below and in Figure 1.

Political will
Usually in developing countries, public service providers and legal systems often lack the firmness and reliability over the medium to long term to institute and enforce policies governing private sector and supportive incentives for the integration of renewable energy.Therefore, the need for governments to give their political backing through the promulgation of energy policies with the necessary regulatory framework that would also encourage the private sector to invest in renewable energy systems for power delivery has become imminent [29].
The components of such policies include provision of long term investment security, drive upgrades of existing national grids and accelerate roll-out of mini-grids, address needs of non-electrical energy forms in productive sectors, energy efficiency improvement measures, and as much as possible, apart from operational and deployment activities, enhance domestic renewable energy technology manufacturing and assembly [30]. •

Resource assessment and suitability
The availability and suitability of any resource ensures that there is sufficient generating capacity to support the demand for targeted level of reliability with minimal cost.Hence carrying out critical assessment to determine the capacity factor of variable, uncertain and spatially diverse renewable energy resource using well-known methods, such as the Weibull distribution for wind assessment [31], is key to establishing the project viability.Data collection and identification of potential constraints that includes assessing of climatic factors, biographical factors, and geographical features such as altitude, temperature, densed urban areas, protected areas, large water bodies, etc. are also essential components in determining the suitability of the project [32].

Social and environmental sustainability assessment
This aspect involves the identification of mitigation measures that will not affect the livelihood of the project area negatively; keeping renewable harvest rates within regeneration rates; keeping emissions to a bare minimum; ensuring that sources of raw materials needed by humanity should not be exhausted for electricity generation especially with regards biomass and hydro power generation; ensuring job creation [33]; and ensuring that livelihood of the inhabitants of the project area is improved.

Technology assessment
Renewable Energy Technologies (RETs) exist in a wide range of options for electricity generation; wind power generation can either be onshore or offshore; solar energy has solar heating, Concentrated Solar Power (CSP) and PV; hydro power has reservoir and run-of-the-river; and biomass power has combustion plants, biogas technology, grate technology, Bubbling Fluidized Bed (BFB) and Circulating Fluidized Bed (CFB), etc.This implies that selection of the preferred technology based on the resources available requires careful considerations ranging from type, fuel flexibility, load ramping capability, investment cost, and plant size [34] and thus the adoption of Clean Energy Technology Assessment Methodology (CETAM) [35].

Study Area
The scope of this research is nationwide.Table 1 lists all the existing generating facilities across the nation.However, not all generating facilities are in full operation due to aging and persistent operational problems.Hence their availability is very low and cannot be fully relied upon.Figure 2 shows the capacity of the existing RETs in the country for a period of 11 years to account for the slow development of RETs in the country.At the moment there is no grid connected PV and wind energy systems.The off grid PV installations mainly constitute solar street lights in selected parts of the country and few mini grids within the range of 6-36 kW.There is also no off-grid installation of wind turbines.Of the 15% electricity access, about 80% is dedicated to supply the capital city, Freetown, located in the Western area.This implies that for the remaining Eastern, Northern, and Southern regions where the bulk of the population resides, access to electricity is practically non-existent.These regions also host the major mining companies, as shown in Table 1, which rely entirely on captive generation to meet their energy requirements since the continuity and quality of energy supply are crucial to their energy requirements.Less than 2% of these companies are directly connected to the existing grid.Figure 3 shows the wide gap existing between generation and demand taking into consideration the utilization of existing energy generation, the operations and production patterns of industries using captive generation and their proposed expansions [37,38].Since the existing grids do not cover the entire country, the accuracy of the profile is limited to ±25%.The existing transmission and distribution infrastructure is one of the major constraints to expanding electricity access covering only about 40% of the residents in the Western area [39].Hence, the industrial and commercial sectors have to greatly depend on running expensive diesel generators, giving rise to over 35,000 off-grid diesel generators.It is against this background that the government has partnered with donor organizations to urgently address the energy needs of the country.One of such partnerships have therefore caused the approval of the Coted'Ivoire, Liberia, Sierra Leone, and Guinea (CLSG) electricity networks interconnection project funded by the African Development Bank (AfDB), which involves the construction of a 1357 km long double high voltage (225 kV) line to connect the national networks of the four countries.The construction of this line is part of the priority projects of the West African Power Pool (WAPP) master plan [40,41] and completion of the construction work will bring about a huge boost in the generation sector as their will be enough transport capacity to cater for the energy needs of the country.The transmission line, as shown in Figure 4, almost evenly divides the entire country into halves, traversing seven districts.Hence in the subsections following, careful resource assessments have been done for wind, solar, and biomass along the proposed and existing grids in an attempt to preparing informed database for policy makers to access in their decision making.

Resource Assessment
Central to the development of variable renewable energy technologies, among other considerations such as demand, the nature of the network and required energy storage, a holistic planning process involves provision of an understanding of where the resources are strongest in order to enable reliability, increased generation, reduction of life cycle costs, helping policy makers develop policy incentives for renewable energy sources with the strongest potential, and, by extension, enable regional energy planners to appropriately reflect renewable energy contributions in their energy master plans [42].Hence in this research well known methods have been used to make an approximate assessment of solar and wind energy potential and characteristics.

Wind Resource Assessment
Making accurate assessments of wind energy potential and characteristics requires long term meteorological observations [43].For this study, one year analysis has been carried out on two locations; Kabala District in the Nothern region and Kenema District in Southern Region, denoted for simplicity as (D1) and (D2) respectively, along the proposed 225 kV as shown in Figure 4. Figure 5a,b shows the wind speed distribution for the two locations at 100 m hub height.Deployment of the renewable energy technology closer to the proposed infrastructure, denoted as (B) in Figure 4, will greatly reduce overall grid connection costs.Point (C) in the figure, denotes the proposed distance (within 30 km) from the existing and proposed grids to install the technologies.Among the many methods known for wind assessment, Rayleigh and Weibull distributions have been widely used to describe the wind speed distribution [44,45].The literature [44] also highlighted on the versatility of the 2-parameter Weibull distribution method and referenced the work of many researchers that have used it for wind profile characterization of any site.Hence, in this study, we have employed this method to carry out sensitivity analysis to determine the (c) and (k) values of the sites (D1) and (D2).The Weibull distribution function or Probability Density Function (PDF) is expressed mathematically as: where f(v) is the probability of occurence of the observed wind speed v, k is the unitless Weibull shape parameter and c is the Weibull scale parameter.

The Cumulative distribution function (CDF) of the Weibull distribution is given by
where v, c, and k are the same as expressed in Equation ( 1).The maximum wind speed v max and maximum frequency of the wind speed v mod are expressed as: Usually, wind power is expressed as: where ρ (kg/m 3 ) is the air density taken to be 1.225 kg/m 3 , A (m 2 ) is swept area.The Weibull wind power density PD, per unit area is given as: where γ is a gamma function.Since wind direction is much useful for positioning of wind turbines at the site, the assessment in this study is limited only to wind speed.Table 2 shows the scale and shape parameters and Power Density (PD) values for the sites.The scale parameter, c, indicates the intensity of the wind in the locations, whilst the shape parameter, k, indicates how peaked the wind distribution is [46].Figure 6a-c shows the mean daily wind speed, PDF, and CDF for both locations.Figure 6d shows a simple model that simulates the electrical output power of the wind turbine generator.V ci , V r , and V co are the cut-in, rated, and cut-off speeds, respectively, where the cut-in speed is necessary to start generating energy, the rated speed being the maximum yield point of the output power and the cut-off speed is the speed at which the generator stops in order to protect its mechanism.The behaviour of this system is usually expressed mathematically as shown in Equations ( 7) and ( 8) [44,47].
The average output power, P w,avg , of the wind turbine is modeled as: The capacity factor, CF, is expressed as: The CF values for both locations have been listed in Table 2.These values are a measure of the energy production efficiency of that facility over a period of time based on the solar resource potential of the site.According to NREL's new transparent cost data base [48], of the 35 projects evaluated, onshore wind has a minimum capacity factor of 26%, a median value of 38%, and a maximum value of 52%.Therefore, from the values obtained, it is clearly seen that D1 location in the Northern province is preferred for wind farm installation as the capacity factor is relatively higher than that of D2 and above the minimum value of 26%.It shows that wind energy generation is possible especially in mountainous regions in the Northern part of Sierra Leone albeit, this value is obtained at 100 m hub height.

Solar Resource Assessment
In the same manner for wind resource assessment, obtaining solar resource information is critical in the design stages of a solar project.Figure 7 shows the global horizontal irradiation (GHI) and suitability areas identified along the existing and proposed grids.Figure 8a,b shows the D1 and D2 solar radiations respectively.Analysis of monthly and annual values of GHI was carried out by the authors in [49,50].For this research, we first obtained the CFs for both sites and were recorded as 27.41% and 25.13% for D1 and D2, respectively.The daily irradiation data for the best and worst days were then analysed to obtain the mean hourly distribution.The frequency count and probability distribution plots were obtained over a period of 1 year as shown in Figure 9a.The averaged hourly irradiance, G t , is obtained from the Gaussian distribution function, f (G t ), as expressed in Equation (12).From the statistical and probabilistic results obtained, it can be observed that the results for both sites are closely marched.However, site D1 proves to be preferable since the CF obtained is slightly higher than the average CF for SSA (27%) [51].Site D2 also shows good potential as the CF is higher than the mean value obtained for Solar projects at the NRELs new transparent cost data base (20.30%).Figure 9b shows the normalized PV generator power curve.
P pv (t) is the generated power at time t and P pvr is the PV generator capacity assumed.
where µ G t is the mean value and σ G t is the standard deviation of the solar radiation data set.

Biomass Resource Assessment
There has been growing interest in the adoption of biomass in energy master plans worldwide mainly due to advances in biomass conversion technologies that bring out low GHG emissions.However, carrying out an assessment to map the location of the feedstock supply base with excess amount of land requiring least amount of resettlement, and quantify the minimum amounts of biomass needed is crucial to confirming the technical and financial feasibility of a biomass project in SSA.Determination of the security of supply, including supplier risks, seasonal variability, and the probable need for supplementary biomass feedstocks in the event of shortage of supply is all encompassing.The authors have highlighted several models for assessing biomass feedstock availability [34,53,54].For this research, we focused on the assessment for the availability of sugarcane bagasse, a residue from sugarcane, to be used as fuel for the Circulating Fluidized Bed (CFB) biomass plant.Sugarcane bagasse has gained notoriety for use in eclectric power generation especially in Brazil, where its continuous usage to supply energy is increasing yearly [55].Kambia district in Northern Sierra Leone has been considered the most suitable location for growing sugarcane because of its proximity to Freetown port (138 K) and Kolente river which provides access to irrigation [56].Characterization of the biomass feedstock, which includes the physical and chemical properties, has been done in Table 3. Annual sugarcane production was obtained over 26 years as shown in Figure 10.The following calculation methodology was developed [34] to quantify the available sugarcane baggase, also known as supply potential, suitable for energy production annually: where Q Bagg is the amount of sugarcane bagasse (tons) available for energy production, S ann is the annual production of sugarcane(tons), K RB is the ratio of sugarcane baggase to primary sugarcane (30%), K AC is the accessibility coefficient (95%), K HC is the harvest coefficient (100%), and K UF is the unused fraction (80%) [34].The current annual sugarcane production has been assumed to be the average production over the 26 years data obtained (70,000 tons/year).This implies that from Equation ( 13), the yearly supply potential of sugarcane bagasse is about 16,000 tons for energy production.Dedicated land would thus be required to increase the supply potential to 25,000 tons in order to meet the maximum desired capacity of the plant (45,000 kW).The selection of the appropriate biomass combustion power plant capacity is discussed in the Biomass System Model Section.

Configuration and Scheme of Hybrid Systems
Ten combinations of hybrid systems have been introduced according to the scheme shown in Figure 11.The letters A, B, C, D, and E are used to designate the PV plant, wind turbines, biomass plant, battery bank, and diesel generator set, respectively, for ease of referencing the components of the blocks (Blocks 1-10) considered.The specifications of each component of the hybrid system for the respective blocks are fed into the optimization algorithm (MOPSO) as clearly illustrated in Figure 11.This process is explained in the next section (Problem Formulation).The configuration of the hybrid systems is representative of block 1 in Figure 12.

Problem Formulation
As seen in Figure 3, the gap between generation and demand is very wide and the grid power is dominated by diesel generation.In order to meet the electricity demand in the most economical and environmentally friendly manner, optimum configuration of wind, PV, biomass, DG, and BESS has been designed for 10 different hybrid blocks within certain levels of system reliability.The MOPSO algorithm is used to determine the optimal values of each block at a time.In addition to the brief highlight given in the Introduction Section, during the optimal operation each particle achieves the best fitness value and maintains that value.Then the particle with the best fitness value is calculated and updated during the iterative process.Each particle represents the configurations of the technologies to be considered (A, B, C, D, and E) for which the search space dimensions and the decision variables shown in Figure 11, will be 5.Each particle will then be evaluated by the multi-objective functions corresponding to LCC, DPSP, DEF, and CO 2 emissions.Hence in summary [58], the first step, which is the initialization involves requesting the input data of hourly radiation, wind speed, national demand requirements, all data to evaluate the technical and economic functions, and finally, data to access constraints for which upper and lower bounds have been stated in Equation (48).The energy output of the technologies are then calculated from Equations ( 14)- (21).Random selection of the position and velocity of particles is done to generate the initial population and applied to the objective functions for each iteration.The fitness value with minimum objective function is then calculated and updated.The updated value is then applied to find the optimal solution of the objective functions.If the objective functions for each combination is satisfied, then the optimal configuration for each combination is selected.Otherwise, the process starts again from evaluating and classifying the objective functions as shown in Figure 11.The optimality of the process is then confirmed after several runs of the optimization process.The optimal values of DPSP (0%, 5%, 10%, 20%, 30%, 40%, and 50%) were deduced as indicators for policy makers to decide on the specific percentage of deficiency they want to achieve whilst carefully selecting which technologies would be appropriate to consider within the scope of the available resources.Therefore ten combinations were designed and evaluated whilst the LCC, DPSP, CO 2 emissions, and DEF were minimized.The LCC constitutes the capital cost (C C ), the Net Present Values (NPV) of the replacement cost (C R ), the resale cost (C S ), and the operations and maintenance costs (C OM ) of the technologies for each hybrid system.

Physical System Criterion
For each of the combinations evaluated, the decision variables of the objective functions are land area occupied by PV panels (A PV ), area occupied by wind turbine's blades (A W ), quantity of baggase available for energy production (Q Bagg ), number of DG sets to be injected as additional back up power (Q DG ), and the battery bank capacity (Q Bat ).The performance of the PV panels and wind turbines depends on the meteorological conditions of the site.The amount of energy generated hourly by the systems is measured in kWh.

Wind System Model
The system behavior of the wind turbine is expressed mathematically as mentioned earlier in Equations ( 7) and (8).The decision variable is the total area swept by the wind turbine blades A w and the rated power from Equation ( 7) is modelled as shown in Equation ( 14) below.
where, ρ a is the air density C p is the power coefficient and η w is the efficiency of the wind turbine.
The rated power, P r , is measured in kW.The adjustment of the wind profile from the reference hub height H R (10 m), to the new hub height H (100 m), is given as [59]: V H represent the wind speed (m/s) at 100 m hub height, V H R is the wind speed at 10 m hub height, and α is the roughness length index taken to be 0.14 for an area with low mature crops.

PV System Model
The decision variable for the PV system is the total area A pv covered by the PV panels at the project site.The power supplied by the PV panels can be expressed as: where G(t) is the hourly irradiation measured in watts per meter squared (W/m 2 ), η pv is the efficiency of the PV module.P pv is measured in kW.

Biomass System Model
For the CFB steam based power plant, the capacity per unit time was constrained to a maximum plant capacity of 45,000 kW, for which the maximum supply potential of baggase has been chosen to be 25,000 tons yearly.The dispatch scheme, illustrated in the reliability flow chart in Figure 13, is used to determine the yearly operating time T, of the plant.The decision variable, chosen to be the yearly supply potential of the plant Q Bagg , is a function of this value T.An increase in the operating time would require a corresponding increase in Q Bagg to maintain the desired output of the plant.This phenomenon is explained in the mathematical model expressed below [60]: where P BM is the plant capacity (kW), Q Bagg is determined from Equation ( 13) to be obtained from the optimization, NCV Bagg is the net calorific value of baggase, η CFB is the efficiency of the biomass plant, and T is the total yearly operating time.The CFB combustion plant was selected because it can be used in the mid range of 10-40 MWe and its excellent emission performance.Additionally, its fuel handling flexibility for feedstocks with low alkali content and typical moisture content values less than 60% [61] gives it advantage for this research scope over stoker grate boilers, which also uses baggase as one of its fuels.

DG System Model
A minimum plant unit, N DG (10,000 kW), was selected for the diesel generator to be injected.The decision variable for this system is the total number of DG units, Q DG , to be obtained from the optimization.The total generated output power of the DG is then formulated as:

BESS Model
The battery energy storage is sized to meet the demand during periods of unavailability in accordance with the hybrid systems reliability flow chart.We have chosen the lithium ion battery bank for this research because of its high-cycle efficiency and fast response times.Recent advances in the technology development has rendered it cost competitive with that of the traditional lead acid battery.Additionally, lithium-ion represented more than 80% of the installed power and energy capacity of large-scale battery storage in operation in the United States at the end of 2016 in comparison with Lead acid, which represent 2%-3%, even though it is one of the oldest forms of battery storage [62].The charging of the battery bank, represented by the state of charge (SOC) at any hour t, is calculated as follows: where P Gen (t) = P grid (t) + P PV (t) + P W (t). The battery bank discharge at any time is given by: The SOC at any given time is subject to the following constraint: where the minimum SOC is a function of the depth of discharge DOD, and nominal capacity of the battery bank C B , given by SOC min = (1 − DOD) × C B with the DOD taken to be 50%.

Economic System Criterion
Evaluating the economic criterion for a project is key in optimizing the size of a hybrid energy system.For this research we have considered the project life N of all blocks evaluated to be 20 years.All economic parameters have been listed in Table A2.As mentioned earlier, the cost component considered is LCC.The constituent components were then evaluated in the next sections.

Wind System Model
Considering an initial cost of installation of one turbine to be α w ($/m 2 ), the capital cost of investment is expressed as: The lifespan of a wind turbine is considered to be 20 years.Therefore the replacement cost is zero (R NPV w = 0).
The annual operation and maintenance (OM w ) cost is expressed as: α OM w × A w ($/m 2 /year).Since the annual growth rate of the cost is considered as µ w and the interest rate is taken to be i, the NPV of the total O&M cost is given as: The resale price is expressed as s w ($/m 2 ).The total cost recovered from resale is given as S w Tot = s w × A w .Hence, the resale price at the end of the project life is given as: The LCC of the wind power system is obtained as follows:

PV System Model
Since the lifespan of the PV system is the same as the project life, the cost analysis follows the same procedure done for the wind system in Equations ( 22)-( 25), whilst replacing the subscript w with the subscript pv.

Biomass System Model
As mentioned earlier, the decision variable for the biomass model is the quantity of baggase, Q Bagg (tons).The total cost of supplying the baggase at the plant site is expressed as: where C baggase , C storage , and C loading are the cost of baggase, storage, and loading per ton of baggase respectively.C transport is the cost per ton of baggase per kilometer [34].The capital cost of investment of the biomass plant is given as: where max(P BM ) (kW) is the maximum output power of the biomass plant obtained from the optimization.The annual operation and maintenance cost, OM BM , is expressed as: (α OM BM × sum(P BM )) + C Bagg Tot ($/kWh/year).With the annual growth rate of the cost considered to be µ BM , the NPV of the total O&M cost is given as: If the resale price is considered as s bm ($/kW), then the total cost recovered from resale is given as S BM Tot = s bm × max(P BM ).Hence, the resale price at the end of the project life is given by: The replacement cost, R NPV BM , is considered to be zero.Therefore the LCC of the biomass power plant is obtained as follows:

DG System Model
For the cost analysis of the DG system considered to operate in dispatch mode as backup power, we employed the unit cost of $1000/kW for 10,000 kW oil fired diesel plant.This cost was an approximate value obtained from cost estimates done for engines with nominal capacities of 40, 50, 100, and 160 MWe [63].The choice of model for this study can be dual fuel MAN Engine, Nigatta, or Wartsila thermal engines that the workforce already has operational experiences with.The capital cost of the DG power plant, C DG , to be installed for this scope of work is thus: Since the DG plant is a dual fuel system, fuel consumption and cost analysis was done for Heavy Fuel Oil (HFO) and Diesel Oil (DO) operation based on the total yearly running hours, n, obtained from the optimization constrained by the dispatch schedule illustrated in Figure 13.The total fuel cost is therefore denoted as τ DG .The annual O&M cost, (OM DG ), is given as α OM DG + τ DG .The replacement cost is zero (R NPV DG = 0).Therefore, for a project life of 20 years, the total O&M (OM NPV DG ) and resale (S NPV DG ) costs are given by: The LCC of the DG power plant is thus obtained as follows:

BESS System Model
The lifetime of the lithium ion battery bank is taken to be 10 years.The resale cost is zero (S NPV Bat = 0) and the capital cost of investment is given by: where C b is the nominal capacity of the battery bank.Since the project life is 20 years, the battery bank will be replaced once.Therefore the replacement cost after 10 years is expressed as: In this case, j = 10.The LCC of the BESS is therefore given by:

Environmental Criterion
In this study, we have considered the Diesel Energy Fraction (DEF) and CO 2 emissions as the environmental performance indicators for each block.The DEF is the fraction of energy supplied from diesel generation.This fraction can be reduced by increasing the share of renewable energy and limiting the diesel operation to backup mode.The CO 2 emissions are therefore taken to be a function of the quantity of fuel burnt annually, Q F , from the existing and injected diesel generation (HFO and DO) and the amount of baggase burnt from biomass generation.CO 2 emissions have been evaluated firstly, for 100% injection of diesel generation at certain levels of reliability from 0% to 50% and then compared against the hybrid blocks for the same levels of reliability as indicated in the case study and discussion of results section.Table A2 of Appendix A shows the Net Calorific Value (NCV) and Emissions Factors (EF s ) for HFO, DO, and Baggase.The CO 2 emissions coefficient, EF s , of the fuel type, i, has been estimated in accordance with the default values given in the IPCC guidelines [64,65].The fuel consumption for the existing dual fuel (HFO and DO) and light fuel oil (DO) engine models, as well as the injected engine model, have been listed in Table A1.Therefore, the DEF and amount of CO 2 emissions (Q CO 2 ) produced yearly for all blocks have been expressed as: where P Renewable(t) is the total renewable generation, P Diesel(t) is the total generation from DG s , Q F,i is the quantity of the ith fuel (HFO, DO, and Baggase) burnt annually, and WEF F,i is the weighted average emission factors of the ith fuel.WEF F,i is calculated as NCV i is the weighted average net calorific value (mmBtu/gal) of the ith fuel and EF CO 2,i is the weighted average CO 2 emissions factor (KgCO 2 /mmBtu) of the ith fuel.The annual CO 2 emissions avoided (Ann CO 2 ,Avoided ) by each block is determined from the difference between the reference emission (R emss ) and the block emissions (B emss ), as shown in the equation below: This value gives policy makers a broad overview of the environmental friendliness of each combination as discussed in Section 5.The reference emission is the annual additional CO 2 emissions by the base case (DG Only) and the block emissions are the annual additional emissions caused by the injection of each block.

Hybrid Systems Reliability Criterion
We have used the DPSP to define the reliability of the system whilst biomass and DG plants are operated in dispatch mode, with priority given to biomass injection when PV and wind generation cannot meet the load demand.The DPSP approach is considered as a technical criterion for sizing a hybrid system that employs a battery bank [66].Figure 13 (reliability flow chart) illustrates the reliability model with the operational sequence summarized below:

•
If the generated power, P Gen , from the injected renewables (PV and Wind) is greater than the load demand, P L , for a given time (t), the excess power is used to charge the battery bank as shown in Equation ( 19), until the SOC(t) = SOC max .The new excess power, EPG, is unused.

•
If P L is greater than the energy generated from PV and wind injection, the deficiency of power, DPS, is calculated and biomass power, P BM , is used to meet the deficit until the maximum intended capacity of P BM is obtained.If the deficit is less than the intended minimum injection of biomass P BM min , then the capacity of P BM is taken to be P BM min .However, if the deficit is greater than the intended maximum biomass power, P BM , then the capacity of biomass power is taken to be P BM max .The new deficit is calculated and diesel generation, P DES , is used to meet the deficit.The minimum and maximum diesel generation, P DES min and P DES max , respectively, are used to constrain the diesel generation in the same manner as described for the biomass injection.

•
When the calculated P DES ≥ P DES max , this implies there is still deficit.In this case the battery bank is discharged according to Equation (20) to satisfy the demand.If the SOC(t) < SOC max the new deficit, DPS(t), is calculated.Otherwise, EPG(t) is calculated.The excess and deficient power for any given hour, t, are thus expressed as follows: where P Gen (t) = P grid (t) + P PV (t) + P W (t).For this research we have calculated the EPG to determine the unused power generated by the system.The Relative Excess Power (REPG) and DPSP are evaluated as shown in the reliability flow chart.Values of DPSP obtained from the optimization that are equal to or closer to 0 implies the load will be always or nearly satisfied.Hence, policy makers can have an insight into the cost and environmental implications simultaneously, when choosing any block for a desired DPSP.

Summary of Objectives
The objective functions considered for this study are thus summarized below: where OF ECO , OF DEF , OF DPSP , and OF ENV are the economic objective or LCC, DEF, DPSP, and environmental objective or CO 2 emissions respectively.

Case Study and Discussion of Results
The technical, environmental, and economic parameters of the PV system, wind farm, biomass plant, DG plant, and battery bank applied in the different hybrid blocks considered are listed in Table A2 of Appendix A. Upper bounds were established on the decision variables based on resource availability and suitability.Therefore, the total area swept by the wind turbines' blades, PV panel project area, quantity of baggase available for energy production, number of DG units, and the battery bank capacity have been varied accordingly below: The total project areas allocated for PV system and Wind farm installations occupies less than 0.03% of the total land area of Sierra Leone (71,740 km 2 ), thereby posing no threat to the living area of the inhabitants.
Two-parameter Weibull distribution and Gaussian distribution methods were used to assess the potential of wind and PV where the capacity factors revealed that wind generation is feasible in Kabala (Site D1) with a capacity factor of 31.6% for 100 m hub height as compared to capacity factor of 9.95% in Kenema (Site D2).Capacity factors of 27.41% and 25.13% also revealed the feasibility of PV injection at sites D1 and D2, respectively.Biomass resource assessment revealed that in order to achieve a maximum injection of 45 MW additional 9000 tons of yearly baggase supply is needed to complement the 16,000 tons of baggase obtained.This would require a dedicated land at or near the plant site as mentioned in the Section 2.2.3.Multi-objective optimization technique using MOPSO was applied to follow the reliability model with the objective of minimizing DPSP, DEF, LCC, and CO 2 emissions.In order to justify the accuracy of the considered algorithm, several runs were made for each block from which the best solutions were deduced.The optimal results were evaluated and compared for different values of DPSP (0-50%).It was observed that due to maintaining the upper bounds of the resources consistently across all blocks, DPSP values ranging from 0-50% were obtained for 5 blocks whilst DPSP values from 30-50% were obtained for the remaining 5 blocks.In order to avoid ambiguity in the analysis, comparison on blocks with DPSP values within the range of 0-50% has been done in Tables 4 and 5. Figure 14 shows the optimal pareto front plots obtained within this range, whilst Figure A1 in Appendix A shows plots for the remaining blocks that has DPSP values in the range 30-50% respectively.

Annual Avoided CO 2 Emissions on Varying Reliability Options for Each Combination
The annual avoided CO 2 emissions for each block have been evaluated against the set of DPSP values as shown in Table 4 and Figure 15. Figure 16 illustrates the actual emissions intensity for the various levels of reliability.DG only block has been considered as base case for comparing the performance of the other blocks.From the results obtained it can be observed that blocks 3 (ABD), 8 (AD), and 9 (BD) are the most environmentally friendly since they have the highest avoided emissions for all levels of DPSP.High levels of CO 2 emissions are realized in the base case scenario (DG Only) due to the daily operation of DGs to meet demand.Figures 15 and 16 clearly shows that Block 1 (ABCDE) has moderately low emission values as compared to blocks 7 (ADE) and the base case.This is clearly seen in the amount of tCO 2 emissions avoided yearly by Block 1 for each DPSP level.This shows that increasing the share of renewable technologies with DG operated only in backup mode can considerably reduce CO 2 emissions.Block 7 is seen to be pretty high for 10% and 20% DPSP.However, in order to reach a conclusion on the choice blocks appropriate for the level of reliability that could be decided by policy makers, LCC values were evaluated as well in the next subsection.

LCC Outcomes on Varying Reliability Options for Each Combination
LCC values were compared for the base case and the blocks considered in this literature in order to determine the most economical grid scenario.Table 5 and Figures 17 and 18 fully describe this scenario.The capital costs of the technologies and on-site capacity factors for the wind and solar PV obtained earlier from the assessments conducted significantly influence the LCC values for each block.Usually, technologies with lower investment costs and higher on-site capacity factors result in lower LCC.We observed that the base case is the most expensive grid scenario due to the significantly high penetration of DGs running for a longer term as seen from the high DEF values, high capital costs, and operations and maintenance costs.For example, considering a 90% electricity access (DPSP of 10%), running DGs only would cost the government an LCC of 6.14 Bn USD over 20 years.Block 1 (ABCDE) was also observed to have comparatively high LCC value for DPSP of 5% (2.23 Bn USD).This outcome resulted from high penetration of low capacity factor solar PVs (27.41%) for the site selected and biomass plant with high capital investment and total operations cost.However, as seen from Table 5 and Figure 17, block 1 has the lowest LCC values from 20% (1.0 Bn USD) to 50% (380 M USD) DPSP for the project life of 20 years, which makes it very competitive with the other blocks.Block 9 (BD) has the lowest LCC values for 0%, 5%, and 10% DPSP.Block 3 (ABD) also has lower LCC values than block 9 for reliability levels of 40% and 50%.

Conclusions
A multi-criteria decision maker for grid-connected renewable energy systems in Sierra Leone was considered as the Sub-Saharan African case.A pentagonal relationship was used to summarize the inter dependency of the key components that should be considered to determine the feasibility of any renewable energy project.A fair resource assessment of solar PV, wind, and biomass power potential was done for the district towns of Kabala, Kenema, and Kambia, where the PV plant, wind farm, and biomass plant will be constructed to be connected along the proposed WAPP 225 kV line and closer to the existing grid.The results from the assessment revealed that Site D1 at Kabala district has the highest capacity factor for PV (27.41%) and Wind (31.6%), indicating that the Northern region is most feasible for installing solar PV panels and wind turbines.The biomass assessment also revealed the need for a dedicated land for an additional 9000 tons of yearly bagasse supply to achieve the maximum capacity of 45 MW.Ten hybrid combinations were studied in order to establish the most economical and environmentally friendly hybrid combinations for the selected technologies.The input data set for the MOPSO algorithm consists of hourly solar global irradiation, wind speed, load requirements for the year 2017, and specifications of the system devices.The optimum results revealed that Block 1, a combination of PV, wind, biomass, BESS, and DG, has low LCC values, 1.0 Bn USD to 380 M USD for a wider range of DPSP values of 20% to 50%, respectively.At 10% DPSP, Bolck 1 has an LCC value of 1.5 Bn USD, which is also very competitive with the LCC of the other blocks for the same level of reliability.The heat tiles (Figure 16 and Table 4) show moderately low CO 2 emissions for Block 1. From the perspective of job creation, the fraction of each renewable energy technology considered in the respective blocks influences the average employment gained in construction, installation, O&M and fuel processing jobs.This study, however, did not deal with the quantitative analysis of job creation due to limited data obtained.However, from the methodologies adopted by [67,68], considering a rich energy mix and creating jobs in the agricultural sector as well, block 1 is the most ideal for a wider scope of sectors, from Agriculture to Construction, for the locations considered.Therefore, we have recommended Block 1 to achieve the long term goals of the government, which is 92% electricity access (a DPSP of approximately 10%) by 2030, as mentioned in Section 1.
In our future work, we intend to include quantitative analysis to accurately determine the job creation potentials of each block for all the reliability levels, whilst evaluating the power trading potential within the Cote d'ivoire, Liberia, Sierra Leone, and Guinea block in the West African Power Pool.As we look forward to making more contributions we hope this information will add to the database of similar work to help policy makers in the region.

Acknowledgments:
The authors acknowledge the unwavering support of the co-wokers of D.A Konneh, operations staff of the Electricity Generation and Transmission Company and the Distribution Supply Authority of Sierra Leone, for helping in providing relevant data with regards the daily operations of the power plants across the country.

Conflicts of Interest:
The authors delcare no conflict of interests

Abbreviations
The following abbreviations are used in this manuscript:

Figure 6 .
Figure 6.Assessment of wind farm locations continued.(a) Mean daily wind speeds of both locations.(b) Weibull probability density function.(c) Cumulative probability density function.(d) Turbine power curve for both locations.

Figure 8 .
Figure 8. Assessment of solar farm locations.(a) D1 solar irradiation in the Northern region.(b) D2 solar irradiation in the Southern region.

Figure 9 .
Figure 9. Assessment of solar farm locations.

2 Figure 18 .
Figure 18.Weighted LCC and CO 2 emissions for each block.

Renewable Energy Project Feasible 1.Political Will 2.Resource & Technical Assessment/ Suitability 5. Economic Assessment (LCCvsDPSP) 3. Social and Environmental Sustainability Assessment 4. Technology Assessment Figure 1. Pentagonal
decision criteria for renewable energy integration.

Table 2 .
Weibull parameters and capacity factor.

Table 3 .
Characterization of the biomass feedstock.

Table A1 .
OM NPV w net present value of annual operation and maintenance cost of wind turbine.NPV of the replacement cost of wind turbine.C Bagg Tot total cost of supplying the baggase at the plant site.BM annual growth rate of the BM cost.OM BM annual operation and maintenance cost of BM. S NPV BM NPV of the resale price of biomass plant.S BM Tot total cost recovered from resale.δ BM initial cost of biomass plant.LCC BM life cylce cost of biomass power plant.OM NPV BM NPV of the total operation and maintenance cost of biomass plant.R NPV BM NPV of the replacement cost of biomass plant.C DG capital cost of the DG power plant.α DG initial cost of DG.OM NPV DG NPV of the total operation and maintenance cost of DG.OM DG operation and maintenance cost of DG. µ DG annual growth rate of the DG cost.S NPV DG NPV of the resale price of DG. S DG Tot total resale price of DG at the end of the project life.Cont.
µ δ Bat initial cost of battery.

Table A2 .
Physical, Environmental, and Economic Parameters.