Optimal Sizing and Techno-Economic Analysis of Minigrid Hybrid Renewable Energy System for Tourist Destination Islands of Lake Tana, Ethiopia

: Achieving universal electricity access is a challenging goal for the governments of devel-oping countries such as Ethiopia. Extending the national grid to the remotely located, scattered, and island populations demands a huge investment. This paper aims to show the techno-economic feasibility of minigrid renewable energy system to electrify Kibran Gabriel island in Ethiopia, through the execution of simulation, optimization and sensitivity analysis using Hybrid Optimization Mod-els for Energy Resources (HOMER Pro) software. The minigrid systems were compared with both diesel generation (DG) and grid extension systems. The hybrid PV/DG/battery system is more economically feasible compared with other minigrid systems, and the best cost-effective option is the one including load flow (LF) strategy with 25 kW of PV, 10 kW of DG, 40 kWh of battery, and 5 kW of bi-directional convertor. The optimal PV/DG/Battery system, having levelized cost of energy (COE) of USD 0.175/kWh, net present cost (NPC) of USD 119,139 and renewable fraction (RF) of 86.4%, reduces the pollutant emissions by 33,102 kg/yr compared with the stand-alone DG system. The optimal minigrid sensitivity to the variations in global horizontal irradiance (GHI), diesel price and load consumption were considered in the sensitivity analysis, and the result shows that the system will operate reasonably well.


Introduction
Energy poverty is a serious problem that must be addressed for the sustainable development of a community. Of different forms of energy, electrical energy is the one that has become part and parcel of our lives. Nowadays, electricity access is closely linked to aspects of human development such as healthcare, sanitation, educational access, and employment rate [1]. That is why it is given priority attention in the sustainable development goals (SDGs)-specifically, as stated in SDG-7: access to affordable, reliable, sustainable and modern energy for all [2]. According to the International Energy Agency (IEA), the number of people in sub-Saharan Africa without access to electricity was estimated to be 580 million in 2019, while the number of people without access to clean cooking facilities was estimated to be about 900 million [3]. Similar studies showed that a number of people without access to electricity reside in rural areas of sub-Saharan Africa and South Asia [4].
The rural off-grid communities rely heavily on conventional energy resources such as firewood, animal dung, kerosene, petroleum, and candles. However, these resources have negative consequences on the environment as well for human health [5].
The main challenge today is to supply electricity to people who are not yet connected to the grid in sub-Saharan African countries, such as Ethiopia. There are four possible In accordance with worldwide trends, the Government of Ethiopia has set key targets for its energy sector, with the goal of providing electric power to all citizens by the end of 2030 [25]. The government has a vision to provide high quality energy services to all sectors of the economy, as well as to export power to neighboring countries via the East African Power Pool [26]. Renewable energy sources such as hydropower, photovoltaics, wind, and geothermal energy have been identified as promising possibilities for Ethiopia's energy security and environmental management in order to mitigate the harmful consequences of climate change [26,27]. However, rural electrification in the remote areas and small islands of Ethiopia through grid extension demands a huge capital investment because of the dispersed population, and difficult geographic terrain. According to the IEA, African energy outlook 2019 special report, Ethiopia has an electricity access rate of only 45%; nearly 60 million people do not have access to electricity [28]. Only 11% of the population has access to electricity through decentralized solutions, mainly from standalone diesel and PV.
Considering the research gaps and the real problem on the ground at the specified site, the authors proposed optimal sizing and techno-economic analysis of the minigrid hybrid renewable energy system in comparison with grid-extension and diesel power generation for Kibran Gabriel island of Ethiopia. In addition to the optimal minigrid system, e.g., other alternative configurations, the pollutant emission gases that can be reduced, the optimal system sensitivity to the variations of GHI, load consumption and diesel price are addressed in the analysis.
The remainder of this paper is organized as follows: Section 2 describes materials and methods, whereas in Section 3 resources and demand of the case study are presented. Minigrid hybrid system discerptions, and component cost and financial assumptions are discussed in Sections 4 and 5, respectively. The simulation results and discussion are provided in Section 6, then a brief conclusion follows in Section 7.

Materials and Methods
This study mainly aims to show the techno-economic viability of the minigrid system for Kibran Gabriel island. Using HOMER Pro software, the feasible configurations and system sizing for the optimal system were carried out. The required data for simulation analysis were acquired from multiple sources: the weather data of the island was obtained from NASA surface meteorology and solar energy database integrated in HOMER Pro [29]. In addition, the load data, diesel fuel price, interest rate and inflation rate, grid extension cost, and e-market surveys for equipment costs were required for the analysis.

•
The load profile of the island was estimated and categorized into primary loads (power is required at a specific time) and deferrable loads (power is required but time is flexible).

•
Since the island is not fully electrified yet, in the simulation analysis, electrifying the island with diesel generator was considered as a reference for other comparisons. Despite the pollution, the economic feasibility of diesel power generation relative to grid-extension cost was also considered for further comparisons.

•
We considered different feasible configurations of PV/wind/battery, PV/wind/diesel/battery, PV/battery, and PV/diesel/battery systems compared with the standalone DG system, and finally the optimal system was determined. • Once the optimal system was determined, sensitivity analysis was carried out to observe the effects of the variation of important parameters such as global horizontal irradiation, diesel fuel price, and load consumption on the optimal system. • Finally, the optimal size and techno-economically feasible configuration were determined. The minimum LCOE, and NPC, the total emissions reduced (compared with stand-alone diesel system), RF, capacity shortage fraction, and BGED and other parameters were computed.

Description of Simulation Software
HOMER Pro software was developed by US National Renewable Energy Laboratory (NREL). The software is a powerful tool for designing, analyzing and planning of minigrid systems to determine the optimal size of system components through carrying out techno-economic analysis [29]. Figure 1 depicts the HOMER Pro simulation process architecture diagram.

Site Description
Kibran Gabriel is one of Lake Tana islands located 11 km North West of Bahir Dar city, Ethiopia [30] with latitude 11.6516° and longitude 37.3644°. Lake Tana is one of the UNESCO world biosphere reserve sites in Ethiopia, registered on 19 June 2015 [31]. It is the home of 37 islands, of which 20 have Ethiopian Orthodox churches and monasteries. Kibran Gabriel is endowed with historical and religious heritages that make it the most popular tourist destination among the Lake Tana islands. However, the communities living on the island do not have access to electricity. Negligible SHS and portable diesel generators are used during holidays. The major economic activities practiced are agriculture, fishery and sightseeing. The geographic location is indicated in Figure 2 below.

Global Horizontal Irradiation and Wind Speed
Providing the latitude and longitude of Kibran Gabriel island to HOMER Pro software, GHI and wind speed data of NASA can be found via Internet connection. According to the data obtained, daily average GHI is 5.94 kWh/m 2 /day. As depicted in Figure 3, the daily average GHI for each month varies from the minimum 5.05 kWh/m 2 /day in July to the maximum 6.68 kWh/m 2 /day in April. The sky is relatively clearer in winter (November to May) than summer (June to October). Similarly, the annual average wind speed data of the site taken from NASA. The annual average wind speed is 3.78 m/s at 50 m above the surface of the earth. As shown in Figure 4, the relative minimum and the maximum monthly average wind speeds are measured in October (2.89 m/s) and June (4.7 m/s), respectively.

Load Profile
The estimated load for Kibran Gabriel island is depicted in Figure 5 below. In this paper, considering the environmental conditions of the island and the community living there, loads are categorized as summer (June to October) load, winter (November to May) load, and weekend load. According to the estimation, a total of 16.00 kWp power, of which 6.00 kWp were deferrable power, 86.24 kWh/day of winter load, 77.84 kWh/day of summer load, and 62.22 kWh/day of weekend loads were considered for the analysis. However, the scaled annual average load taken for the analysis was 76.94 kWh/day. In addition, depending on the demand and its time of use, the load is categorized into primary and deferrable load. The primary loads are loads that require a certain amount of operating reserves to respond immediately to the sudden increase in the electric load and/or due to a decrease in power output [32], whereas deferrable loads permit some flexibility when power is insufficient. That is, deferrable load is the second priority after the primary load; however, its importance comes ahead of charging batteries. In this study, lighting loads, fans, sound amplifiers, refrigerators, sewing machines, flourmills, and water pumps were the major loads estimated for the analysis. However, the flourmill and sewing machines were excluded from the weekend loads considering the community religious practices in the island. Similarly, since the summer season is rainy and relatively cold, water pumps and fans were excluded from the summer load. In addition, water pumps are considered as deferrable load during the winter season. The estimated load frequency distribution is shown in Figure 6. According to the scaled frequency distribution, 75% of the distribution is 2 kW-5 kW, and the remaining 25% is 6 kW-22 kW.

Diesel Price
Even though the diesel fuel price market in Ethiopia is relatively fluctuating from time to time, the government subsidizes it to minimize the cost. In this study considering the local market, a diesel fuel price of USD 0.62/L is taken for the analysis [33].

Economic Assignment Criteria
In this study real interest rate, capital recovery factor (CRF), net present cost (NPC), cost of energy (COE), and break-even grid extension distance were considered as the economic evaluation criteria [32,34].

Annual Real Interest Rate
As stated in Equation (1) HOMER Pro uses the annual real interest rate to determine discount factor and annualized costs from present costs.
where i is annual real interest rate, i' is nominal interest rate, and f is expected inflation rate (%).

Capital Recovery Factor (CRF)
The CRF is the ratio that is used to compute the present value of the annuity in the project lifetime.
where i is the annual real interest rate (%), and t is the number of years.

Net Present Cost
According to HOMER Pro, total NPC value signifies the project life cycle cost, and is computed using Equation (3). The cost includes the initial capital cost and the cash flow of the t-year over the factor (capital cost, fuel cost, operation cost, maintenance cost, replacement cost, etc.). Income contains electricity selling and salvage value after project life cycle.
CF0 is initial capital cost (USD); CFt is cash flow of t-year (USD); t is number of year (yr); i is annual real interest rate (%); N is project lifetime (yr).

Cost of Energy
HOMER Pro defines COE as the average cost per kWh of useful electrical energy generated by the system [15]. It is calculated by dividing total annualized cost (TAC) by total annualized useful electrical energy generated as shown in Equation (4).
TAC is the annualized value of NPC in USD/yr, Eserved is total electrical load served (primary load, deferrable load, energy sales to the grid) in kWh/yr, CRF is capital recovery factor, i is the real interest rate (%) and N is the project lifetime (yr).

Grid-Extension Cost
In this research, the grid extension cost is determined through the computation of breakeven grid extension distance (BGED). The BGED is the distance away from the grid at which the NPC of expanding the grid equals the NPC of a minigrid (stand-alone) system. When the grid extension distance exceeds the BGED, the minigrid system is more cost-effective [19,32]. The BGED is computed using Equation (7).
where total NPC of the minigrid system is in USD, CRF is capital recovery factor, i real discount rate in %, N is project lifetime in years, Edemand is total annual electrical demand including primary and deferrable in kWh/yr, cpower is cost of power from gird in USD/kWh, ccap is capital cost of grid extension in USD/km, com is O&M cost of grid extension in USD/yr/km.

Electrical Assignment Criteria
The descriptions below are some of the electrical assignment criteria references according to HOMER Pro [34].

Renewable Fraction
Renewable fraction (RF) is the fraction of the energy delivered to the load that is originated from renewable power sources. The percentage RF is expressed as shown in Equation (8).

Excess Electricity Fraction
Excess electricity fraction (fexcess) is the ratio of total excess electricity to total electrical production as shown in Equation (9).

Capacity Shortage Fraction
Capacity shortage fraction (fcs) is the ratio of total capacity shortage and total electrical demand as described in Equation (10). HOMER Pro considers a system viable only if the capacity shortage fraction is less than or equal to the maximum annual capacity shortage [35].

Sensitivity Variables
In the rural areas, there is a low load demand which varies on a time basis [4]. On the other hand, solar radiation and wind speed vary depending on the daily environmental conditions. Sensitivity analysis monitors the effect of certain variables. That is, different values are assigned to these variables within a given range in order to assess their influence on the optimal minigrid system [19]. In this study, considering the optimal minigrid system, the sensitivity variables with the variations of GHI from 5.05 to 6.69 kWh/m 2 /day, fuel price from USD 0.62 to 1/L, and load consumptions from 76.94 to 107.72 kWh/day are considered for the analysis.

MiniGrid System Schematic Diagram
The schematic diagram of the minigrid power system is depicted in Figure 7. The minigrid system consists of wind power generator (WG), PV system (PVs), diesel generator (DG), power convertor (PC), storage system, primary and deferrable loads. PVs and WG are operated in parallel with DG, and AC coupled. PVs and WG supply power to the load to reduce the output of DG in order to reduce fuel intake [35]. Whereas the grid system is connected only for BGED analysis.

PV System
The PV module used for this study is CS6K-295MS mono-crystalline silicon solar cell module, which has the parameters shown in Table 1. The PV module is expected to serve for 25 years and the power output of the module (Ppv) is calculated according to Equation (11) [36,37].
where Ppv,STC is PV rated power power (kWp), fpv is PV derating factor (%), GT is solar irradiance on the surface of the PV module (kW/m 2 ), GT,STC is standard solar irradiance (1 kW/m 2 ), KP is the PV module temperature coefficient of power (−0.40%/°C), TC is temperature of the PV module (°C), TSTC is temperature of the PV module under standard test conditions (STCs) (25 °C). In this study a derating factor of 85% was considered.

Wind Generator System
A Bergey (XL10)) wind turbine of 10 kWp AC power output at 12 m/s, a minimum operating wind speed of 2.5 m/s, and 20 years of lifetime was considered for this analysis [39]. The power curve of the wind turbine, and the wind speed at the hub (Uhab) are shown in Figure 8, and Equation (12), respectively.
where Uref is reference wind speed (m/s), given by NASA measurements at a reference height of 50 m from the surface (Href). Hhub is the hub height (30 m in this case) where the wind speed at the hub is calculated. α is surface roughness-the value is set to 0.14 [40,41].

Diesel Generator
Due to the fact that the PV and wind modules have intermittent natures, the diesel generator will ensure the off-grid hybrid power system's stability and reliability. In addition, in this study the diesel generator is considered as the base case reference for relative comparison using LCOE and NPC. In HOMER Pro, we input the load profile and select ''Autosize Genset'' for the DG specification, and set the maximum annual capacity shortage fraction of the system operational constraint to 0%. The generic 10 kW diesel generator with efficiency of 30%, and fuel curve slope of 0.286 L/h/kW were considered for the initial case from HOMER Pro catalogue. Then, the software automatically adjusted to the proper capacity and fulfilled the system's operational constraints.

Power Convertor and Storage System
A bi-directional DC-AC converter was employed for power conversion. When surplus power is generated by PV, diesel, and/or WG, the excess AC power is converted to DC and stored in a battery, to supply it back when necessary. A conversion efficiency of 95% and 15 years of lifetime are assumed in the simulation [33]. In this study, a Li-ion battery with a round-trip efficiency of 90%, an initial state of charge of 100%, a minimum SOC of 20%, and 10 years of lifetime with negligible O&M cost is considered [42].

Grid System
Kibran Gabriel island is located about 11 km from Bahir Dar city; however, extending the grid from the city to Kibran has not yet been realized due to Lake Tana. In this study, we consider the grid system to compute BGED using Equation (7) so as to compare the proposed off-grid systems cost with the grid extension cost.

Dispatch Strategy
Excess energy from renewables is always expected to charge the battery, and one of the most important decisions that the dispatch system must make is when to charge the batteries with energy from the generator. In this research, we look at a Combined Dispatch (CD) approach that combines Load Following (LF) and Cycle Charging (CC) dispatch strategies [43]. In the CC strategy, whenever the generator needs to be on to meet the net load (load minus intermittent renewable power), it operates at full capacity and the surplus power charges the battery bank. In LF strategy, the generator is dispatched just to fulfill the net load and any operational reserve needs. The CD strategy determines whether to charge the battery with the generator based on the current net load, employing the CC strategy when the net load is low and the LF strategy when the net load is high [20,32].

MiniGrid Component Cost
Taking into account the local market and the reduction in system component costs from time to time, the summary of component costs depicted in Table 2 were considered for the simulation analysis. As shown in the table, the cost includes capital cost, replacement cost, and operation and maintenance (O&M) cost. The capital cost of the PV (USD 900/kW) includes the PV inverter cost, and the capital cost of the WG (USD 3500/kW) includes wind turbine, 30 m tower and wind invertor costs. The estimate costs of all the components also include transportation and installation costs.

Grid Extension Component Cost
In this study, in order to verify the feasibility of minigrid power supply to Kibran Gabriel island, BGED was determined using Equation (7). For the BGED analysis, the minimum grid extension capital cost of USD 23,000/km and O&M cost of USD 200/km were considered [46]. With the government subsidies, a grid power price of USD 0.06/kWh is taken for the analysis according to Ethiopian Electric Utility revised tariff [47].

Interest Rate and Inflation Rate
As depicted in Figure 9a,b the nominal interest rate and inflation rate, respectively, announced by National Bank of Ethiopia [48] for the past seven years (2013-2019) were considered for the simulation analysis. The minimum nominal interest rate of 11.88%, the maximum nominal interest rate of 13.5%, and the average of the nominal interest rate for the considered period was 12.59%. Similarly, the minimum, maximum, and average inflation rate for the considered period were 7.4%, 14.6% and 10.51%, respectively. In this paper, the project lifetime of 25 years, the average values of nominal interest rate (12.59%) and the inflation rate (10.51%) were used to calculate the annual real interest rate using Equation (1).

Simulation Results and Discussion
The proposed minigrid system depicted in Figure 7 is optimized using HOMER Pro based on LCOE and NPC. Once the optimal minigrid system was determined, various sensitivity variables were considered to check the effect on the optimal system. According to the simulation result, the possible configurations for the selected site are discussed as follows:

Stand-Alone Diesel System
As depicted in Table 3, when only diesel power generation is considered the LCOE of USD 0.346/kWh, NPC of USD 235,734 and 38.3 t/yr of pollutant emission gases were the optimization results. However, since the simulation result of only diesel system is the reference for other configurations, to make it cost effective excesses energy of not more than 3% and a maximum capacity shortage of 0.09% were considered. Once the optimization result of the stand-alone diesel system was determined, its feasibility relative to the grid extension cost was determined using BGED analysis. According to BGED analysis, a stand-alone diesel system is relatively feasible for the selected site compared to grid extension cost when the electrification distance is farther than 3.11 km, as shown in Figure  10a. However, due to the high level of polluting gases emitted, the high operating fuel cost, and the tourist-attracting nature of the island, the stand-alone diesel system is not the favored solution to electrify the specified island.

PV/Wind/Battery System
The PV/wind/battery system is the other alternative to electrify the specified site. According to the simulation analysis this system was optimized to the LCOE of USD 0.323/kWh and NPC of USD 219,736. The PV/wind/battery system is free of pollutant gas emissions, and economically feasible compared with the stand-alone diesel system. However, the PV/wind/battery system demands a relatively high initial investment cost (USD 129,986) due to its high component cost compared with the other configurations considered in this study.

PV/Wind/Diesel/Battery System
Of the hybrid systems considered in this study, the PV/wind/diesel/battery system, having LCOE of USD 0.304/kWh, NPC of USD 207,388, and emissions of 18,324 kg/yr, was composed of a better energy mix than the other systems discussed in Sections 6.1.2, 6.1.4 and 6.1.5. Even though the PV/wind/diesel/battery system had a better in energy mix than the other configurations considered in this study, the higher magnitudes of NPC, LCOE Grid extension Stand-alone diesel Grid extension PV/Diesel/Battery and emissions, as shown in Table 3, could lead us to search for another cost effectivealternative to electrify Kibran Gabriel island.

PV/Battery System
The PV/battery system is the second-lowest-cost alternative for electrifying the specified island. The simulation result indicated that the LCOE of 0.258 USD/kWh, NPC of USD 175,506 with zero pollutant emission gasses. In this system, in order to deliver the desired power to the community with 0% capacity shortage, a PV size of 60 kWp, and a battery size of 75 kWh were required. However, due to the intermittent nature of solar energy, the PV/Battery system must be of a relatively big size to cope up the day and night loads. Having a big size leads the system to have high excess electricity (66.3%). On the other hand, when we calculated a capacity shortage of 5%, the excess electricity was reduced to 35.6% and the unmet load was 4.5%. This assured that in the PV/battery system, when 0% capacity shortage is considered, the system size will be higher. When a capacity shortage of greater than zero is considered, the system is incapable of handling peak loads. Considering these limitations, the authors continued to search for another optimal solution to electrify the island.

PV/Diesel/Battery System
The initial simulation conditions considered for this study were: annual average solar radiation of 5.94 kWh/m 2 /day; annual average wind speed of 3.78 m/s; nominal interest rate of 12.59%; primary load demand of 76.94 kWh/day; and a project lifetime of 25 years. From the feasible system configurations shown in Table 3, the PV/diesel/battery system offers the optimal solution for the proposed system. The detailed is presented in Table 4. As depicted in Table 4, the PV/diesel/battery system could deliver the required power to the island with LCOE of 0.175 USD/kWh, NPC of USD 119,139 and emissions of 5159 kg/yr. When these parameters compared to the simulation outputs of the stand-alone diesel system, the pollutant emission gasses, NPC, and LCOE were reduced by 86.52%, 49.46% and 49.42%, respectively. The BGED, and a detailed cost summary of components, are shown in Figures 10b and 11, respectively.  On the other hand, the technical performance of the PV/Diesel/battery system was simulated considering the winter and summer environmental conditions with the corresponding loads. The result shown in Figure 12 displays the detailed time series analysis of the power outputs of the PV and diesel system, battery charging and discharging, and total loads (including primary and deferrable) for the winter's day of January 3. In the winter time, due to high solar radiation (up to 6.68 kWh/m 2 /day) and clearance index (up to 0.69) in the country, specifically at the study site, the PV system generates significant amounts of power from 06:00 to 18:00. The deferrable loads (6 kW) include water pumps were served from 06:00 to 11:00 when there is excess power. However, during the peak loading condition (10:00 to 12:00) the diesel generator supported the PV system to fulfill the power demand and the battery charging power started to decline at 11:00. During the nighttime, from 18:00 to 06:00 the load demand receives power from battery and/or diesel generator. As shown in Figure 12, starting from 17:00 to 04:00 the battery discharged to deliver the required power to the load demand. However, from 04:00 to 06:00, when the battery was incapable of fulfilling the load demand, the diesel generator started to produce the required electricity. Similarly, the time series analysis of power outputs of the PV and diesel system, battery charging and discharging, and total loads for the summer's day of July 6 is depicted as follows in Figure 13.   During the summer season in the study area experiences a relatively high amount of rain, with lower solar radiation (5.05 kWh/m 2 /day) and lower clearance index (0.48). As shown in Figure 13, the total load and the primary load are the same (overlapping) because the deferrable loads (water pumps) were neglected since the summer is rainy season in the study area. From the time series analysis, we observed that the main power source during the daytime (07:00 to 17:00) was the PV system, whereas for the remaining time, from 17:00 to 03:00 the battery system, and from 03:00 to 07:00 the diesel generator contributed the power demanded by the load. As we can observe from Figures 12-14a the, PV power output of the summer season was relatively lower than the PV power output of the winter season. The yearly power contributions of the PV and diesel generator, and the state of charge of the battery for the PV/diesel/battery hybrid system were depicted in Figure 14. As shown in Figure 14a during, the daytime, especially from 8:00 to 16:00, the PV system contributed a relatively high power share, whereas the remaining time and the peak load conditions were handled by the diesel generator and the battery system, as depicted in Figure 14b,c.

Sensitivity Analysis
The optimization process that enabled us to find out the optimal solution (PV/diesel/battery system) used average values such as the annual average solar radiation of 5.94 kWh/m 2 /day, the annual average wind speed of 3.78 m/s, a primary load demand of 76.94 kWh/day, an ad fuel cost of USD 0.6/L. However, in the real system the solar radiation and wind speed vary depending on the daily environmental conditions, the load demand varies dynamically, and the fuel cost fluctuates depending on the world's fuel market. Having this in mind, the effects of the variations of GHI, fuel price, and load consumption on the optimal PV/diesel/battery system were considered in the sensitivity analysis.

Global Horizontal Irradiance (GHI)
Based on the data taken from National Renewable Energy Laboratory Database, the minimum GHI (5.05 kWh/m 2 /day) and the maximum GHI (6.69 kWh/m 2 /day) were the ranges on which the sensitivity variables varied. Considering ±7% changes of the average value (5.94 kWh/m 2 /day), the sensitivity variables of GHI were 5.05, 5.46, 5.94, 6.35, 6.68 kWh/m 2 /day. As we can observe from the optimal plot (Figure 15), maintaining the scale average load (76.96 kWh/day) when diesel price is minimum (USD 0.6/L) and GHI is maximum (6.68 kWh/m 2 /day), the minimum LCOE is USD 0.170/kWh. On the contrary, when the diesel price is maximum (USD 1/L) and GHI is minimum (5.05 kWh/m 2 /day) the maximum LCOE is USD 0.198/kWh. The effect of GHI variation on LCOE, NPC and RF can also be observed from Figures 16 and 17. When we set the load to 100.01 kWh/day and the diesel price USD 0.62L, the LCOE was lowered to USD 0.168/kWh and RF increased to 79.4%, as depicted in Figure 16. Keeping the minimum diesel price and load but increasing GHI (from 5.05 to 6.68 kWh/m 2 /day) will reduce the NPC from USD 124,086 to USD 116,059 (Figure 17), the LCOE from USD 0.182 to USD 0.170/kWh (Figure 15), and increase the RF from 82.7% to 87.3%.

Diesel Price
Considering the local market conditions, we set the fuel price to vary from USD 0.62-1/L in steps of 15% change from the minimum value (USD 0.62/L). Keeping the minimum load (76.94 kWh/day) and minimum GHI (5.05 kWh/m 2 /day), when the diesel price varied from USD 0.62 to 1/L the LCOE increased from USD 0.182 to 0.198/kWh ( Figure 15) and NPC also increased to USD 134,656. Similarly, the line plot of diesel price variation effect on the LCOE and NPC was presented in Figure 18 above. Maintaining the average GHI of 5.94 kWh/m 2 /day and load of 76.94 kWh/day, when the diesel price varied from USD 0.62 to 1/L, relative increments of LCOE and NPC from the optimal value were observed in Figure 18.

Load Consumption
As shown in Figure 19, when the load varied from 76.94 to 107.72 kWh/day in steeps of 10% change from the minimum (keeping the current diesel price of USD 0.62/L and average GHI of 5.94 kWh/m 2 /day) the LCOE showed a small reduction from USD 0.175 to 0.171/kWh, while the fuel intake increased from 1895 L/yr to 4473 L/yr. Similarly, the NPC increased from USD 119,160 to USD 154,780. As electricity demand rises, DG will increase the power generation to meet the demand. This causes the RF to decline, the fuel consumption and NPC to increase, and COE to become relatively decreased. As we have seen from the sensitivity analysis, the variations in the GHI, diesel fuel price and load have been considered to check the capabilities of the PV/diesel/battery system to withstand those fluctuations. When these sensitivity parameters were considered, the new optimization results achieved were LCOE of USD 0.179/kWh, NPC of USD 151,468, and RF of 69.1%. These results prove that the PV/diesel/battery system can operate reasonably well for the Kibran Gabriel island during the specified project lifetime.

Conclusions
Rural electrification in the remote regions of Ethiopia through grid extension is not economically feasible because of the dispersed population, difficult geographic terrain and low load demand in the existing conditions. On the other hand, delivering electricity from SHS and diesel generators only for lighting for a short period of time will not bring the expected economic development. Having these controversies in mind, in this study, in addition to lighting loads, sewing machines, flourmills, water pumps, and refrigerators for preserving fishery products were considered in the load estimation of Kibran Gabriel island.
Based on the estimated load demand and the environmental conditions of the island, stand-alone diesel, PV/battery, PV/wind/battery, PV/diesel/battery, and PV/wind/diesel/battery systems were considered in the analysis. The following concluding remarks are drawn.

•
Of all the minigrid systems considered in this study, the PV/diesel/battery system, having a size of 25 kW of PV, 10 kW of DG, and 40 batteries of 1 kWh each is the optimal minigrid system, with LCOE of USD 0.175/kWh, NPC of USD 119,139, and RF of 86.4%. Relatively, it is the lowest cost system among the considered configurations.

•
Compared to the base case stand-alone diesel system (LCOE of USD 0.346/kWh, NPC of USD 235,734), the optimal PV/diesel/battery system will reduce pollutant emissions by 33,102 kg/yr.

•
In addition to PV/diesel/battery system, the PV/battery, PV/wind/battery, and PV/wind/diesel/battery systems are economically feasible compared to the standalone diesel system. However, the initial investment cost of the stand-alone diesel system is relatively the lowest.

•
According to the BGED analysis, grid-extension to the island within the specified project lifetime is extremely costly compared to all the proposed options considered in this study.
• Once the optimal configuration is determined, the optimal minigrid sensitivity for the variations of GHI, diesel price, and load consumption was considered in the sensitivity analysis. The new optimization results due to the variations of the sensitivity parameters were modified to LCOE of USD 0.179/kWh, NPC of USD 151,468 and RF of 69.1 %. These optimization results showed that the system will operate reasonably well with the variations of GHI, diesel price, and load consumption. However, the impacts of renewable energy resources on the reliability and quality of power, and refining the estimated loads for practical implementation will be the subjects of future works. Informed Consent Statement: Not applicable.

Data Availability Statement:
The freely available data can be accessed through the indicated references whereas the remaining parts can be obtained from the corresponding author upon the reasonable request.