Economic feasibility analysis and optimization of hybrid renewable energy systems for rural electrification in Peru

The majority of rural communities in developing countries (such as Peru) are not connected to the electrical grid. Hybrid energy production from available renewable resources (e.g., wind and solar) and diesel engines is considered as an economically viable and environmentally friendly alternative for electrification in these areas. Motivated by the lack of a comprehensive investigation dedicated to the techno-economic analysis of hybrid systems (PV–wind–diesel) for off-grid electrification in Peru, the present work is focused on determining the optimal configuration of these systems for remote Peruvian villages. Three small communities without access to the grid (Campo serio, El potrero, and Silicucho), which are located in different climatic zones of Peru, have been accordingly selected as case studies. Seven different configurations including single component systems (solar, wind, and diesel) and hybrid ones are considered. While taking into account the meteorological data and load characteristics of the communities along with the diesel fuel’s price and the cost of components, HOMER software is utilized to determine the optimal sizing of the system [resulting in the lowest net present cost (NPC)] considering different scenarios. The obtained configurations are then compared considering other state-of-the-art economic indices [initial capital cost, total annual operating cost, and the cost of energy (COE)], the generation fractions, and the resulting CO2 emissions. The obtained results have revealed that, for all of the investigated communities, the hybrid solar–wind–diesel system is the most economically viable scenario. Considering the latter scenario, the obtained optimal configuration leads to an NPC of USD 227,335 (COE: 0.478 USD/kWh) for Campo serio, USD 183,851 (COE: 0.460 USD/kWh) for El potrero, and USD 146,583 (COE: 0.504 USD/kWh) for Silicucho. Furthermore, employing the optimal configurations a renewable fraction (with respect to the total generation) of 94% is obtained for Campo serio and Silicucho, while the achieved renewable fraction for El potrero is 97%. Moreover, for the case of Campo serio, the resulting CO2 emission of the obtained optimal system is determined to be 6.1% of that of a diesel-only unit, while the latter ratio is determined to be 2.7% for El potrero and 9.9% for that of Silicucho. The optimal configurations that are obtained and presented in the present paper can be utilized as guideline for designing electrification systems (with a minimized cost) for the considered communities and other villages with similar characteristics (population and climatic conditions).


Introduction
Power access is at the forefront of governments' preoccupations, particularly in nations in which electricity is essential for certain basic activities such as lighting, refrigeration and running of household appliances (Kanase-Patil et al. 2010). The vast majority of rustic communities in developing countries like Peru are not entirely connected to electrical grid due to geographical obstacles and small population, which make the required investments for grid extension unjustifiable. Energy production from available and sustainable sources such as wind and sun has been considered as a viable and environmentally friendly alternative (Mamaghani et al. 2016b). Due to the ever-increasing price of petroleum derivatives on a global scale and concerns regarding the emission of greenhouse gases (GHG) (Najafi et al. 2015), remarkable attention has been directed toward green renewable technologies for catering growing energy demand (Mamaghani et al. 2017).
Renewable energy sources (RES) are abundant in most parts of the world, and, unlike fossil fuels, can be harnessed without any cost for the resource. In this regard, many governments have already started to finance renewable technologies by means of direct grants, loans and tax incentives (Liu et al. 2012). Electricity generation in Peru through hydro, wind, solar, geothermal, biomass, tidal power or other RES is subjected to an annual maximum 20% depreciation regime for income tax purposes (Irena 2014).
Despite the aforementioned upsides of RES, there are a number of technical difficulties which must be resolved to make renewable energy systems reliable and self-sufficient. The most important issue with the electricity generation by RES is the lack of stability which stems from their intermittent nature. Such nature is due to the variations in the atmospheric conditions which can result in substantial fluctuations in power generation seasonally or even daily (Bekele and Boneya 2012). To provide a stable supply of electrical power, the application of power storage systems such as batteries (Fragaki and Markvart 2008) or combining RES with non-renewable technologies such as diesel, natural gas or biomass driven generators (Montuori et al. 2014) have been proposed in the previous studies. Hybrid energy systems, which usually comprise of at least two power sources, have been utilized to reach higher electrical efficiency and more uniform power supply. Another shortcoming of RES is the significantly higher capital cost of such systems, compared to the conventional diesel generators. Nevertheless, with respect to stand-alone diesel generators, hybrid systems require less fossil fuel (Kalantar and Mousavi 2010) (i.e., reducing operational cost), reduce carbon footprints (Bentouba and Bourouis 2016), while enhancing quality of service (Valente and De almeida 1998). For some configurations of RES projects including the wind turbine-based units, the cost of energy (COE) and net present cost (NPC) can be notably decreased by expanding the capacity of the installed power system and addressing a higher electrical load (Diab et al. 2016). Since several factors [such as different possible combinations of RES and non-RES sources, the demand profile, seasonal availability of RES, capital cost of components, and the fuel cost (Aminyavari et al. 2016)] must be considered while attempting to determine optimal configuration of hybrid energy systems, optimization tools should be utilized in order to obtain a comprehensive evaluation of different scenarios (Gu et al. 2017). Genetic algorithm (Najafi et al. 2011), cuckoo search algorithm (Ray et al. 2018), modified electric system analysis (Zahboune et al. 2016), and game theory (Khare et al. 2016) are some of the frequently utilized techniques that are employed for optimizing the sizing of components (Rajkumar et al. 2011) aiming at minimizing the cost (Yousefi et al. 2017) of these systems. Arceo et al. (2018) demonstrated that, by utilizing the optimal configuration of a hybrid electrification system in remote areas of Western Australia, the overall environmental impact is reduced by 16%, although it leads to increasing the total life cycle costs by 4%. Roy and Kulkarni (2016) determined the optimal configuration of a PV-diesel generator hybrid system for rural areas in India, using which 70% of the energy demand is met by the diesel generator and 30% is addressed by the PV panels. Flores et al. (2016) optimized (through minimizing the COE) the configurations of hybrid wind, PV, and biomass-based generation for rural electrification in Honduras. Hrayshat (2009) showed that utilizing the optimal configuration of hybrid wind-diesel generation units in remote Jordanian settlements leads to an annual reduction of 21.3% in the diesel consumption.
Several works have utilized hybrid optimization model for electric renewables (HOMER) software to perform techno-economic feasibility study, sensitivity analysis, and optimization (Singh and Baredar 2016) on hybrid microgrids (Dekker et al. 2012). The optimal configuration of PV and battery system, which was obtained using HOMER in a study conducted by Alsharif (2017), was demonstrated to be an energy efficient and cost-effective alternative for supplying heterogeneous cellular networks. In another study conducted using HOMER, Brandoni and Bošnjakovic (2017) demonstrated that using the obtained optimal configuration of a hybrid system (constituting of PVs, wind turbines and internal combustion engines) 33-55% of the energy demand of a wastewater facility located in Sub-Saharan Africa can be addressed, while the COE is lower than the local cost of electricity. Marneni et al. (2015) instead used HOMER to find the optimal sizing of solar photovoltaic generation to enhance the voltage profile of a rural feeder (3.06 MW peak load) in Mysuru, India.
It is generally accepted that answering the electrical demand by hybrid systems (i.e., more than one RES or non-RES) is more logical than only depending on RES. This stems from the fact that naturally reliance on a single RES necessitates over-sizing the system to be able to cater the demand considering the variations in solar irradiation/ wind speed throughout the day or seasonally. Many studies have been dedicated to performance evaluation and feasibility analysis of hybrid systems such as PV-wind units (Arribas et al. 2010), wind-diesel-battery, and wind-fuel cell systems (Khan and Iqbal 2005). Shaahid and Elhadidy (2007) performed a techno-economic feasibility analysis on a hybrid PV-diesel-battery system. In a similar study conducted for the case study of Ireland, it was found that, owing to the climatic characteristics of the area, wind is the dominant component of the majority of optimal hybrid power systems (Goodbody 2013). In a study comparing PV-diesel and PV-battery systems, it was concluded that the former is by far more cost-effective for loads higher than 13 kWh/day while the latter is more economical for 3-13 kWh/day range (Lilienthal 2015). PV-wind-diesel-battery hybrid system was observed to be the most practical option (Bekele and Palm 2010) to supply a community model living in an Ethiopian remote area.
Regarding the environmental impacts of hybrid systems, Hafez and Bhattacharya (2012) assessed the emissions of a microgrid arrangement including diesel, wind, PV, battery, and hydro and the CO 2 emission were estimated to be 1078.4 t/year. A study conducted by Ajlan et al. (2017) revealed that the PV-wind and PV-wind-diesel can reduce the CO 2 emission by 100% and 70% and the COE by 30% and 45% with respect to the diesel generators. In another study (Shezan et al. 2015), it was shown that the optimal configuration of wind-DG-battery systems results in a renewable fraction (RF) of 0.0914. In another conducted by Hossain et al. (2017), which was focused on electrification for tourist resorts in Malaysia, it was concluded that the obtained optimal configuration of wind-diesel-battery hybrid system results in a lower COE (0.279 vs 0.343 USD/kWh) and CO 2 emission (2,571,131 kg/year vs 5,432,244 kg/year) compared to the diesel-only system; which make it a promising alternative for reducing the carbon emission intensity  in these areas. Rajbongshi et al. (2017) showed that the COE for a grid-connected hybrid system is lower than an off-grid one for similar load profiles due to the fact that grid-tied systems allow export of excess electricity into the grid rather than storing it with battery. Regarding the impact of projected variations in influential parameters, Dorji et al. (2012) stated that 20% drop in the price of PV led to 4.9-8% decrease in NPC.
Despite the promising potentials of RES for power production in Peru and existence of abundant resources, feasibility studies to explore green and cost-effective technologies such as PV or wind are scarce. To the best of our knowledge, there is no thorough study on techno-economic analysis of hybrid systems (PV-Wind-Diesel) in Peru. The present work aims at finding the optimal combination of available RES to satisfy the energy demand of three offgrid villages in Peru. These territories have been selected according to geographical and population consultation centre of Peru (INEI 2012). INEI (2012) provides statistical data such as population, access to power network, and distance from large urban areas. Meteorological data of solar irradiation and wind speed were borrowed from NASA atmospheric science data center using the location of each community (NASA 2017). Solar irradiation, wind speed, and electricity demand of each community are provided as inputs to HOMER software to conduct the techno-economic analysis. Seven different possible scenarios including single component systems (diesel, solar, or wind) and hybrid ones (Solar-Wind, Solar-Diesel, Wind-Diesel, or Solar-Wind-Diesel) are considered. Each of these scenarios is modeled in HOMER software and, considering the NPC as the economic main index, their corresponding optimal sizing is determined. The obtained configuration is next evaluated based on various economic and environmental criteria including the initial capital cost, operating cost, COE, RF and pollutants emission rates. Since HOMER is utilized in this study, it incorporates all of the limitations of the corresponding models; however, it has extensively been recognized as a promising tool in the literature to achieve the defined aims of the present study.

Location and population
Peru has one of the lowest rural electrification rates in Latin America (INEI 2012). Considerable efforts have been made to extend Peru's national grid as part of the country's plan to augment rural electrification from 71 to 88.2% between 2013 and 2017 (SE4All 2014). For this study, three regions with drastically low power coverage, Cajamarca, Puno and Loreto, have been selected (SE4All 2014). Electrification via grid network to these areas is inconceivable because of impassable routes through the rough mountains, forests and dispersion of the settlements. The vast majority of population in these areas suffer from poverty and cannot afford expenses of electrical installations and maintenance (Rehman et al. 2007). Based on the master plan of rural electrification with RES of Peruvian ministry of mine and energy (MINEM Peru 2008), three small residential communities were chosen from previously mentioned regions.
The regions with lack of access to electricity grid and the chosen communities within these regions are demonstrated in Fig. 1. Loreto is Peru's northernmost region with very small population due to its remote location in the Amazon Rainforest. Cajamarca is situated in the northern highlands of Peru with a tropical climate which is characterized by mild and sunny days. Puno is located in south-eastern of Peru at high elevation which results in harsh climate conditions compared to the other selected communities. Table 1 summarizes the geographic and demographic information and climatic characteristics of these communities.

Load estimation
Current population of each community has been extracted from the last national surveys provided by INEI (2012). Energy demands of each town were estimated by considering the fundamental needs of individuals such as lighting and communication. The maximum connected load for each house is assessed to be 81 W including lighting (36 W fluorescent lamps), communication (15 W radio receiver) and a reserve (30 W). Number of houses is determined considering the population of each community. It should be highlighted that since people in these areas are mostly involved in farming and livestock breeding, they spend most of their day outside of the house, and, thereby, most of the energy consumption takes place in the early night. Two essential services, a health center and a school with loads of 1.4 KW and 1.3 KW, were taken into account during analysis. For the school eight 36 W fluorescent lamps are considered as well as computers, printers and a TV. It was assumed that students attend the classes from 8:00 a.m. to 13:00 p.m. For the case of health center, appliances such as refrigerator for storing medicines, radios and lighting were considered. Public lighting of 400 W for each group of 5 houses was considered in the calculation of load profile for these areas. These approximations are made based on the guideline presented by Bekele and Boneya (2012), which introduced an approach to obtain the load distribution during the day. Figure 2 displays the hourly distribution of the overall load during workdays and weekends in Campo Serio. For the case of brevity only the load profile of this community is presented. It can be observed that the electricity consumption declines in January, February, and March in Peru schools are closed from January to March as a result of the termination of the academic year.
The determined overall load profile of each community is then provided to HOMER. Next, in order to add randomness to the load data to make it more realistic, timestep-to-time-step and day-to-day variability values of 20% and 15% are applied to the provided raw data. The peak load in the generated profile for Campo Serio, El Potrero, and Silicucho is determined to be 18.9 kW, 15.1 kW and   Weekday from Apr to Dec Weekday from Jan to Mar Weekend 9.8 kW. Figure 3 depicts the average load, based on the generated load profiles, for each months of the year for the considered communities. The determined average yearly load for Campo serio, El potrero, and Silicucho is 4.3 kW, 3.6 kW, and 2.6 kW.

Solar irradiation
Utilizing the climatic information from NASA surface meteorology and solar energy (NASA 2017), the average monthly solar irradiations are determined for each location (Fig. 4a). Among selected communities, El potrero has the highest average yearly irradiation (3.145 kWh/m 2 ) and Silicucho has the lowest one (2.326 kWh/m 2 ). As can be noticed in Fig. 4a, during winter (months of May-July), solar irradiation reaches its lowest values, leading to small power production from the PV source.

Wind speed
Average monthly data for wind speed have been extracted from NASA surface meteorology and solar energy (NASA 2017). As displayed in Fig. 4b, Silicucho is the windiest region among all of selected communities with average wind speed of 7.16 m/s. Due to the fact that HOMER calculates the generated power based on the wind speed at the hub height, the altitude of each area should be given as input to the software. HOMER uses logarithmic law to assess the speed of wind at the hub height, as described in Eq. (1).
where z hub is the hub height of the wind turbine in m , z anem is the anemometer height in m , z 0 is the surface roughness length in m , V z hub is wind speed at the hub height of the wind turbine in m/s and V z anem is wind speed at anemometer height in m/s.

Diesel fuel
After high growth rates for several consecutive years and reaching around USD100/barrel in 2013 (Kojima 2016) oil price has experienced a sharp fall in the recent years and in 2017 it reached ca. USD 50/barrel. Following these changes, diesel price estimation for remote areas in Peru has decreased from 1.41 USD/L in 2013 (GIZ 2013) to 0.8 USD/L in 2017 (Global Petrol 2017).

Major components
The major components of hybrid systems are photovoltaic panels, wind turbines, diesel generators, batteries, and power converters. Table 2 summarizes the capital, replacement and maintenance costs of all component of the hybrid system used in HOMER. (

Months
Campo serio El potrero Silicucho

Diesel generator
Reliability of power systems which solely consist of renewable energy technologies is a matter of concern due to intermittent characteristics of RES (Arribas et al. 2010). To provide energy with high reliability while taking into account the sustainability, diesel generator system can be integrated with renewable energy systems. Diesel generators (Generac) with 15 kW, 20 kW, 30 kW or 48 kW rated power have been selected for this study and the corresponding costs are extracted from the corresponding website (Generac 2017). The operating lifetime is considered to be 20,000 h (HOMER Energy 2017a; Generac 2017). The replacement cost is considered to be 90% of the capital cost and yearly operation and maintenance cost is assumed to be 10% of the capital cost (Mamaghani et al. 2016a). To avoid wet stacking, which is a phenomenon that happens when diesel generator is idle for a long interval, minimum load ratio for diesel generator has been chosen to be 30% (HOMER Energy 2017b).

Photovoltaic array
Solar panel electricity generation systems, also known as solar photovoltaic, harness solar energy using photovoltaic cells and convert it into DC electricity energy in direct where f PV is the PV derating factor; Y PV the rated capacity of the PV array in kW; I t the global solar irradiation incident on the surface of the PV array, in kW/m 2 ; and I S the standard amount of irradiation used to rate the capacity of the PV array which is 1 kW/m 2 . Derating factor originates from the PV panel power loss due to dust on the surface of the panels, shading and aging of panels, and elevated temperature. Following HOMER software's recommendation (HOMER Energy 2017a), the derating factor is considered to be 0.9. Two low-cost 260 W and 305 W solar panels manufactured by Astronergy, which are based on Polycrystalline cells, have been chosen as the basis and the corresponding prices are obtained from the associated website (Astronergy 2017 (2) P PV = f PV Y PV I t I s

Wind energy system
Wind turbine is a device that converts kinetic energy of wind into electrical power. The power output for the wind turbine can be calculated employing Eq. (3) (Aeolos 2017).
where v is wind speed in m/s, s is rotor swept area in m 2 , η 1 is generator efficiency, c p is maximum power coefficient and ρ is air density in kg/m 3 . Wind turbines used in this analysis are Aeolos-H 10 kW and Aeolos-H 20 kW horizontal turbines with DC output and rated capacities of 10 kW and 20 kW. For the case of Aeolos-H 10 kW, cut-in wind speed is 2.5 m/s, rated wind speed is 10 m/s, cut-out wind speed is 18 m/s and the sweep area is 44.2 m 2 . For the case of Aeolos-H 20 kW instead, cut-in wind speed is 3 m/s, rated wind speed is 10 m/s, cut-out wind speed is 20 m/s and the sweep area is 78.5 m 2 (Aeolos 2017). The capital cost of the wind energy system, obtained from Aeolos wind turbine catalogue (Aeolos 2017), accounts for the wind turbine, tower, inverter, wiring, painting, corrosion package, foundation, road construction and installation (Aeolos 2017). The replacement cost is considered to be 90% of the capital cost, and the considered lifetime for both models, extracted from the manufacturer's catalogue, is 20 years. The estimation of operating and maintenance (O&M) costs appears to be strongly correlated with turbine's age. In the first few years of operation, considering the manufacturer's warranty, a low O&M expenses can be expected; however, after 10 years, greater O&M costs are considered owing to the needed repairing and the corresponding required reinvestments (Morthorst 2017). Based on the recommendations of (3) P wt = 0.5 v 3 s 1 c p

Battery
Battery stores electricity in a chemical form and subsequently this stored energy can be recharged and reused. Considering the fact that power production from renewables can fluctuate abruptly with weather conditions, to guarantee a reliable and constant power supply, battery storage has been utilized extensively. Battery storage can also be used to reduce the number of start/stop cycles of the backup diesel generator in hybrid systems (Nfah and Ngundam 2009). Two different types of lead acid batteries, Surrette (Surrete-Ltd 2017) and Hoppecke (Hoppecke-Batterien-GmbH 2017) have been chosen as HOMER equivalent batteries for this analysis. The prices of each type of battery are extracted from the corresponding manufacturer's websites (Hoppecke-Batterien-GmbH 2017; Surrete-Ltd 2017); according to which, the lifetime of Surrette and Hoppecke can be considered to be 17 and 20 years. The replacement cost is assumed to be 90% of the capital cost. The yearly operation and maintenance cost is assumed to be 10% of the corresponding capital cost.

Power converter
In systems in which both AC and DC generating elements are dealt with, converter is an essential part of the system. HOMER software uses a converter, which is equivalent to both inverter and rectifier. The inverter is one of key components of a PV system which converts DC power from PV and battery output into AC one. In case of excess energy generation from wind turbines, rectifier converts AC power to DC to be stored in the battery storage system. SolarEdge converter with capacities of 5 kW, 6 kW, 7.6 kW and 10 kW have been selected for this analysis, the corresponding costs are obtained from the manufacturer's website (SolarEdge-Technologies-Inc. 2017). The corresponding replacement cost is assumed to be 90% of the capital cost (Mamaghani et al. 2016a) and the yearly operation and maintenance cost is considered to be 10% of the capital cost (SolarEdge-Technologies-Inc. 2017). Converter and rectifier efficiencies are considered to be 98% and 85%.

Scenarios
To conduct a comprehensive investigation on each location which covers all possible combinations of RES and diesel generator, seven different case scenarios are considered, among which the most favorable one for microgrid planning is determined. HOMER software has been employed to perform the economic optimization for each scenario and community. As depicted in Fig. 5, the configurations analyzed in HOMER can be listed as: • Diesel generator (Case 1) • Photovoltaic (Case 2) • Wind turbine (Case 3) • Solar-wind hybrid (Case 4) • Solar-wind-diesel hybrid (Case 5) • Solar-diesel hybrid (Case 6) • Wind-diesel hybrid (Case 7)

Materials and methods
While taking into account the meteorological data and load characteristics of the communities along with the diesel fuel's price and the cost of components, the optimal (most economical) sizing of the system utilizing different scenarios is determined. In the optimization procedure, the NPC is considered as the key economic index. The obtained configurations are then compared considering the other state-of-the-art economic indices together with the environmental metrics and the generation fractions. The indices which are employed in this study, including NPC, COE, the initial capital cost, and the total annual operating cost along with the generation fractions are presented in this section.

Net present cost (NPC)
NPC represents the cost of a system during its lifetime. The considered costs include capital, replacement, O&M (operation and maintenance), fuel, salvage, and penalties due to emissions. NPC can be calculated using Eq. (4): where TAC is the total annualized cost, CRF is the factor of returning capital given by Eq. (5): where N is the defined life time of the project (number of years) and i is the annual interest. The replacement cost is utilized in order to take into account the economic impact of employing the components with a lower life time than that of the project.

Levelized cost of energy (COE)
COE is the average cost per kilowatt hour (USD/kWh) of useful electrical energy produced by the system which is calculated using Eq. (6).
where E ls is the electrical energy that the microgrid is supplied with and E gr is the amount of electricity sold to the grid by microgrid and C TANN is total annualized cost of the system (USD/year).

Initial capital cost
The initial capital cost is the total cost of all of the components at the beginning of the operation. In a system composed of batteries, wind turbines, solar panels, diesel generator and converter, the initial capital cost can be calculated utilizing Eq. (7).
where N b is number of batteries, C b is the price of battery energy system per each battery, P st is the power produced by all the solar panels in kW, P sb is the power produced by one solar panel in kW, C s is the price of solar energy system per each solar panel, P st is the power produced by all of the wind turbines in kW. P wb is instead the power produced by one wind turbine in kW, C w is the price of wind energy system per each wind turbine C w is the price of wind energy system per each wind turbine, P gt is the total power produced by the diesel generator in kW, C gref is the price of 1 kW generator. P gt is the power of the converter, P gt P cb is the power of base converter and C c is price of the base converter.

Total annual operating cost (TOC)
The operating costs are the expenses corresponding to the operation of all of the components of the system in 1 year. Total operating cost can be calculated using Eq. (8): Fig. 5 The off-grid configurations analyzed in HOMER where n is the number of components of the system, C OM,i is the annual operation and maintenance (O&M) cost for ith component of the system, C f is annual total fuel cost, C R,i is annualized replacement cost for ith component of the system and C S,i is the salvage value of component i. In HOMER, the components are considered to depreciate linearly; the salvage value of a component is directly proportional to its remaining life.

Fractions of generation
The RF, as specified in Eq. (9), is the ratio of the electrical production originating from renewable power sources to the total electrical production: where E ren is renewable electrical production (kWh) and E tot is the total electrical production (kWh). Based on the approach proposed by Celik (2002), the fraction of energy that each source is generating can be subsequently calculated. As expressed in Eqs. 10 and 11, photovoltaic fraction E pv F pv and the wind generation fraction E pv F pv F wg are expressed by the ratio of energy produced by solar (E pv ) and wind (E wg ) generation units to the total generated electrical energy E tot .
The fraction of energy produced by the diesel generator can be similarly defined (Eq. 12).
where E pv E dg is the energy produced by the diesel generator.

Results and discussion
For each considered community, the achieved optimal sizing of the system, considering the cases 1-7, is represented in Tables 3, 4 and 5. As can be observed in these tables, to compare the proposed configurations from the economic point of view, four indicators, namely initial capital of the system (in terms of USD), operating cost (USD/year), total NPC (USD) and COE (USD/kWh), were selected. The latter indices provided an in depth information about the economic perspective of each system and facilitate the determination of the most economically convenient case for each location. To assess these configurations from an environmental standpoint, the yearly CO 2 production for each case has been calculated. The overall yearly production of PV, wind and diesel generation units are also determined for each case.

Optimized designs for Campo serio in Loreto region
The optimization results and the corresponding costs of investigated scenarios for Campo serio are represented in Fig. 6 and Table 3. As can be observed in Fig. 6, considering the NPC, case 5 (solar and wind and diesel) is the most economical scenario. It can be noticed though that for case 1, in which only the diesel generator is used to cater the demand, the required initial capital cost is much lower than the other cases; indeed, its corresponding initial cost (USD 12,149) is almost one-ninth of that of case 5. Despite the fact that the low initial capital investment for diesel generators seems desirable to investors, a long-term analysis reveals that case 5 is the most cost-effective configuration resulting in the lowest NPC. By comparing the operating hour of the diesel generator for cases 5 and 1, indicated in Table 3, it is observed that the corresponding value is diminished by 21 times through utilizing PV and wind technology, which is desirable considering the corresponding emission reduction (Asrari et al. 2012). Considering the operating hour of diesel generator, case 5 is indeed demonstrated to be the most promising configuration among the case which include a diesel generator (cases 1, 5, 6, and 7) resulting in the highest RF. Figure 6 also demonstrates that case 2 (stand-alone solar) is the most costly one resulting in a significantly higher NPC with respect to the other cases. The key difference between the cases 2 and 5 is the power provided by the diesel generator during the winter in which the lack of solar irradiation tremendously decreases the power output of solar panels. Table 3 also demonstrates that three case scenarios with the lowest NPC are case 4, case 5, and case 7. In case 4 and case 7, the majority of energy is produced by wind turbines, while the majority of supplied energy in case 5 is provided by PV panels. It can also be observed that case 4,  with the RF of 100%, is the greenest configuration among the scenarios with the least NPCs. Figure 7 shows the productions of solar, wind and diesel generation units for the optimal configuration of case 5 throughout the year. It can be observed that during the months with a notable solar irradiation (October to February), the PV generation constitutes the major role in the overall electricity production. On the other hand, in the winter months, and particularly during the months of May, June, and July, the PV generation is almost negligible and the corresponding lack of production is compensated by the diesel generator. Table 4 and Fig. 8 represent the optimized results and the corresponding costs of all investigated scenarios for El potrero. As observed in Table 4, the three case scenarios with the lowest NPC are case 4 (solar, wind), case 5 (solar, wind, diesel) and case 6 (solar, diesel). Solar panels supply 87% of overall generation in case 6, while wind turbines produce 81% and 50% of total generation in case 4 and case 5. Among the three cases with the lowest NPC, case 4 is the most environment-friendly scenario leading to the renewable fraction of 100%. As can be observed in Fig. 8, similar to the obtained results for Campo serio, case 5, with the lowest NPC, is the most economical configuration and case 1 is the configuration requiring the lowest initial cost. The latter is due to the low cost of diesel generators compared to renewable energy system components, though the elevated cost of this plan in the long run, and particularly its emissions, devalue the benefit of utilizing this scheme. As demonstrated in Table 4, diesel generator operating hour in case 5 is approximately one fifties of that of case 1.

Optimized designs for Silicucho in Puno region
The obtained optimal sizing and the resulting costs, employing different scenarios, for Silicucho are demonstrated in Table 5 and Fig. 9. As indicated in Fig. 9, three case scenarios with the lowest NPC are case 4, case 5 and case 7 and, similar to the previous regions, case 5 is first most economically viable configuration among the hybrid systems, case 1 still has the lowest initial cost, and case 2 is the scenario leading to the highest NPC. Case 6, as shown in Table 5, is the most environmentfriendly configuration among the scenarios running partly with diesel generator. The resulting carbon dioxide emission employing the configuration of case 6 is far less than the one corresponding to case 1.
Although, as demonstrated in Fig. 4b, the average wind speed in Silicucho is notably higher than the other two areas, the wind generation fraction of the corresponding optimal configuration is considerably lower than the one of Campo serio and El potrero. The latter is due to the fact that, as shown in Fig. 3, the average electricity consumption and peak load of Silicucho is lower than the other two villages. As demonstrated in Table 6, while the economic optimization results in choosing a 20 kW wind turbine for Campo serio and El potrero, a 10 kW wind turbine is selected for Silicucho. The corresponding wind-based power production is considerably lower which leads to a notably lower wind generation fraction. The latter fact demonstrates the necessity of conducting a comprehensive optimization in which the available renewable resources, load characteristics, and the economic indices are taken into account.

Conclusions
Techno-economic performance of stand-alone electricity generation systems for off-grid communities located in different climatic areas of Peru was investigated. Seven scenarios, including different combinations of diesel generators, wind turbine units, and solar panels, were assessed. Optimal sizing of each configuration, which minimizes the corresponding NPC, was determined and the achieved optimal systems were also evaluated considering other economic indices and their environmental performance. The analysis demonstrated that, for all of the investigated communities, the hybrid solar-wind-diesel system is the most economically viable configuration. For the case of Campo serio, although the initial capital cost of the dieselonly configuration is almost one-ninth of that of the mentioned hybrid system, the latter results in the lowest NPC in a long-term analysis. For the case of Campo serio, the optimal configuration requires an initial capital of USD 116,092, and corresponds to an operating cost of 6359 USD/y, a total NPC of USD 227,335 and a COE of 0.4738 USD/kWh. The optimal hybrid configuration for Campo serio resulted in renewable fraction of 94% leading to 4578 kg/y of CO 2 emissions which is 6.1% of the emissions of the diesel-only configuration (23,858 kg/y of CO 2 emissions). The latter ratio is determined to be 2.7% for the case of El  potrero and 9.9% for Silicucho. The variance in the obtained ratios is due to differences in the availability of renewable sources and the load characteristics of the considered areas. The Peruvian authority can play a notable role in facilitating the utilization of such technologies in the rural areas. A depreciation regime for the income tax is the only support which is presently provided to the RES-based electricity generation plant in Peru. In case adequate incentive policies would be provided, the COE of the proposed system will be notably reduced which will aid the mentioned communities to install the proposed systems.