Optimal Sizing and Scheduling of Hybrid Energy Systems: The Cases of Morona Santiago and the Galapagos Islands

Hybrid energy systems (HESs) generate electricity from multiple energy sources that complement each other. Recently, due to the reduction in costs of photovoltaic (PV) modules and wind turbines, these types of systems have become economically competitive. In this study, a mathematical programming model is applied to evaluate the techno-economic feasibility of autonomous units located in two isolated areas of Ecuador: first, the province of Galapagos (subtropical island) and second, the province of Morona Santiago (Amazonian tropical forest). The two case studies suggest that HESs are potential solutions to reduce the dependence of rural villages on fossil fuels and viable mechanisms to bring electrical power to isolated communities in Ecuador. Our results reveal that not only from the economic but also from the environmental point of view, for the case of the Galapagos province, a hybrid energy system with a PV–wind–battery configuration and a levelized cost of energy (LCOE) equal to 0.36 $/kWh is the optimal energy supply system. For the case of Morona Santiago, a hybrid energy system with a PV–diesel–battery configuration and an LCOE equal to 0.37 $/kWh is the most suitable configuration to meet the load of a typical isolated community in Ecuador. The proposed optimization model can be used as a decision-support tool for evaluating the viability of autonomous HES projects at any other location.


Introduction
Today, the world faces unprecedented environmental and economic challenges due to the unsustainable use of natural resources. The growth of the world population, the rising global energy demand, the depletion of conventional energy sources, and the volatility of fossil fuel prices have triggered the enthusiastic search for new alternative energy sources [1]. Renewable energy sources (RESs) have become an essential and promising option to face the abovementioned challenges [2][3][4]. As a result, in the last two decades, researchers, decision-makers, and industry leaders have prioritized the technological development and widespread deployment of low-carbon energy technologies [5]. Decentralized renewable units comprised of technologies such as wind turbines, photovoltaic (PV) solar panels, and micro-hydropower have been proposed as suitable solutions for a more sustainable and clean electricity supply [6][7][8]. However, despite the advantages of renewable energy technologies, the generation of electricity from RESs presents numerous challenges to their deployment, including policy barriers, technological barriers (e.g., intermittency and uncertainty due Energies 2020, 13, 3933 2 of 20 to its dependence on meteorological conditions, wind speed, or sunlight availability), and economic and financial barriers (e.g., high upfront costs, capital and operation costs) [4,9,10]. Hence, to provide a reliable and adequate electricity supply, there is a need for the hybridization of RESs with conventional generating sources [11].
Hybrid energy systems (HESs) generate electricity from multiple energy sources that complement each other. Historically, the global efforts to research and utilize these types of systems increased after the oil crisis in the early 1970s; since then, the number of HES installations has experienced a strong upsurge [12]. Generally, HESs incorporate one or more renewable energy sources (e.g., solar, wind, hydro, fuel cell), a conventional energy source (e.g., microturbines, gas turbines, and internal combustion engines), and an energy storage system, such as batteries, compressed-air energy storage (CAES) systems, flywheels, ultracapacitors, and superconducting magnetic energy storage (SMES) systems [13,14]. Furthermore, their configuration varies depending on their potential application, the availability of local resources, and the energy demand in a considered area. The basic components of an HES are shown in Figure 1.
Energies 2020, 13, x FOR PEER REVIEW 2 of 21 to its dependence on meteorological conditions, wind speed, or sunlight availability), and economic and financial barriers (e.g., high upfront costs, capital and operation costs) [4,9,10]. Hence, to provide a reliable and adequate electricity supply, there is a need for the hybridization of RESs with conventional generating sources [11]. Hybrid energy systems (HESs) generate electricity from multiple energy sources that complement each other. Historically, the global efforts to research and utilize these types of systems increased after the oil crisis in the early 1970s; since then, the number of HES installations has experienced a strong upsurge [12]. Generally, HESs incorporate one or more renewable energy sources (e.g., solar, wind, hydro, fuel cell), a conventional energy source (e.g., microturbines, gas turbines, and internal combustion engines), and an energy storage system, such as batteries, compressed-air energy storage (CAES) systems, flywheels, ultracapacitors, and superconducting magnetic energy storage (SMES) systems [13,14]. Furthermore, their configuration varies depending on their potential application, the availability of local resources, and the energy demand in a considered area. The basic components of an HES are shown in Figure 1. In the literature, HESs have been classified based on several criteria, such as types of technology, type of energy produced, capacity size, and level of integration with the grid [15]. Figure 2 shows an overview of the HES classifications. There are two main classifications of HESs on the basis of their interaction with the distribution grid. An HES that is connected to the distribution grid and that can buy (import) or sell (export) electricity is called an on-grid system. Conversely, an HES that is located in a remote area and therefore not connected to the distribution grid is called an off-grid system, also referred to as autonomous, standalone, or islandic. Another common classification of HESs is based on the criterion of installed capacity. HESs range from small-scale units with a few kilowatts to largescale systems with enough capacity to power entire communities or industrial establishments. Smallscale hybrid energy systems have an installed capacity of less than 5 kW and are used for typical loads of remote homes, telecommunication systems, or meteorological stations. Medium-scale HESs have a capacity in the range of 5 to 100 kW and are used to electrify small communities in rural areas, isles, or mid-sized utilities. Lastly, large-scale HESs have an installed capacity greater than 100 kW and, in most cases, are connected to the distribution grid or the national grid [16]. In the literature, HESs have been classified based on several criteria, such as types of technology, type of energy produced, capacity size, and level of integration with the grid [15]. Figure 2 shows an overview of the HES classifications. There are two main classifications of HESs on the basis of their interaction with the distribution grid. An HES that is connected to the distribution grid and that can buy (import) or sell (export) electricity is called an on-grid system. Conversely, an HES that is located in a remote area and therefore not connected to the distribution grid is called an off-grid system, also referred to as autonomous, standalone, or islandic. Another common classification of HESs is based on the criterion of installed capacity. HESs range from small-scale units with a few kilowatts to large-scale systems with enough capacity to power entire communities or industrial establishments. Small-scale hybrid energy systems have an installed capacity of less than 5 kW and are used for typical loads of remote homes, telecommunication systems, or meteorological stations. Medium-scale HESs have a capacity in the range of 5 to 100 kW and are used to electrify small communities in rural areas, isles, or mid-sized utilities. Lastly, large-scale HESs have an installed capacity greater than 100 kW and, in most cases, are connected to the distribution grid or the national grid [16]. A particularly interesting feature of small-and medium-scale hybrid energy systems is their potential to supply electricity to remote unelectrified villages and rural areas in developing countries [17]. In 2018, the International Energy Agency (IEA) estimated that 992 million people did not have access to electricity and projected that by the year 2030, around 674 million would remain without access [18]. Consequently, local governments, international organizations, and private firms have developed a myriad of strategies and proposed numerous solutions aimed at expanding electricity services to households in remote villages. A conventional strategy is extending the national grid to the territories without electricity coverage. Unfortunately, this solution is prohibitively expensive, since the cost of transmission increases significantly with the remoteness of the location, and in most cases, it is impractical due to the small load demand of these regions [19,20]. Thus, a viable alternative solution that has been widely implemented for the electrification of rural areas in developing countries is the use of RESs to form hybrid energy systems [21].
An increasing number of articles report on the feasibility of hybrid energy power generation systems in developing countries. Particularly, there have been considerable efforts in the development of sizing methods and optimization techniques for these types of systems. For a comprehensive literature review encompassing the domains of optimization, planning, and sizing of hybrid energy systems, see [22][23][24][25]. The authors of [26] carried out a feasibility study of hybrid diesel-RES and grid-RES energy systems in the rural village of Sheikh Abolhassan, Iran. The study consisted of simulating different configurations of standalone HESs and evaluating the economic impact of grid connection and extension to the village with the use of the Hybrid Optimization Model for Electric Renewables (HOMER) computational tool. In [27], the authors applied a linear programming model for the study of hybrid off-grid energy systems in the unelectrified villages of Narendra Nagar, India, and Titutmate, Colombia. The proposed model served as a tool to understand the impact of crucial parameters such as fuel prices, load demand, and investment costs in the optimal configuration of the HESs. In [28], the authors used a simulated annealing (SA) algorithm-based chaotic search and harmony search to find the optimal size of a hybrid energy system while considering the total life cycle costs of a wind turbine, PV collectors, and an energy storage system. Their study focused on the remote region of Kerman, Iran, and the simulation was implemented in MATLAB. The authors of [29] presented an energy management strategy based on a receding horizon optimization (RHO) approach to determine the optimal capacity size of a hybrid renewable energy system (HRES) while satisfying the load demand and minimizing the operation, maintenance, and environmental costs. In the model, photovoltaic panels and wind turbines were considered as the main power technologies, with a diesel generator and a battery storage as backup units. Their study demonstrated the benefits of combining real-time weather forecasting, an RHO procedure, and a demand response strategy to increase the overall efficiency of the HRES.
In [30], the authors used a genetic algorithm (GA) to find the optimum capacity size of an autonomous HES comprised of photovoltaic (PV) panels, a backup source (either a microturbine or diesel generator) and a battery unit. The objective function minimized the overall production and A particularly interesting feature of small-and medium-scale hybrid energy systems is their potential to supply electricity to remote unelectrified villages and rural areas in developing countries [17]. In 2018, the International Energy Agency (IEA) estimated that 992 million people did not have access to electricity and projected that by the year 2030, around 674 million would remain without access [18]. Consequently, local governments, international organizations, and private firms have developed a myriad of strategies and proposed numerous solutions aimed at expanding electricity services to households in remote villages. A conventional strategy is extending the national grid to the territories without electricity coverage. Unfortunately, this solution is prohibitively expensive, since the cost of transmission increases significantly with the remoteness of the location, and in most cases, it is impractical due to the small load demand of these regions [19,20]. Thus, a viable alternative solution that has been widely implemented for the electrification of rural areas in developing countries is the use of RESs to form hybrid energy systems [21].
An increasing number of articles report on the feasibility of hybrid energy power generation systems in developing countries. Particularly, there have been considerable efforts in the development of sizing methods and optimization techniques for these types of systems. For a comprehensive literature review encompassing the domains of optimization, planning, and sizing of hybrid energy systems, see [22][23][24][25]. The authors of [26] carried out a feasibility study of hybrid diesel-RES and grid-RES energy systems in the rural village of Sheikh Abolhassan, Iran. The study consisted of simulating different configurations of standalone HESs and evaluating the economic impact of grid connection and extension to the village with the use of the Hybrid Optimization Model for Electric Renewables (HOMER) computational tool. In [27], the authors applied a linear programming model for the study of hybrid off-grid energy systems in the unelectrified villages of Narendra Nagar, India, and Titutmate, Colombia. The proposed model served as a tool to understand the impact of crucial parameters such as fuel prices, load demand, and investment costs in the optimal configuration of the HESs. In [28], the authors used a simulated annealing (SA) algorithm-based chaotic search and harmony search to find the optimal size of a hybrid energy system while considering the total life cycle costs of a wind turbine, PV collectors, and an energy storage system. Their study focused on the remote region of Kerman, Iran, and the simulation was implemented in MATLAB. The authors of [29] presented an energy management strategy based on a receding horizon optimization (RHO) approach to determine the optimal capacity size of a hybrid renewable energy system (HRES) while satisfying the load demand and minimizing the operation, maintenance, and environmental costs. In the model, photovoltaic panels and wind turbines were considered as the main power technologies, with a diesel generator and a battery storage as backup units. Their study demonstrated the benefits of combining real-time weather forecasting, an RHO procedure, and a demand response strategy to increase the overall efficiency of the HRES.
In [30], the authors used a genetic algorithm (GA) to find the optimum capacity size of an autonomous HES comprised of photovoltaic (PV) panels, a backup source (either a microturbine or diesel generator) and a battery unit. The objective function minimized the overall production and emission costs while meeting the electricity demand of a small rural community located in Palestine. Their analysis indicated the importance of adequately designed HESs since they are capable of reducing not only energy costs but also greenhouse gas emissions. The authors of [31] investigated the potential application of a cuckoo search (CS) algorithm in finding the optimal size of a remote case study in Almora District of Uttarakhand, India. The CS algorithm performance was compared to other optimization algorithms like genetic algorithm and particle swarm optimization (PSO). Their results showed that CS algorithms outperform other well-known optimization methods due to their fast convergence. In [32], the authors focused on the feasibility study of an HES to supply electricity to a rural community in Santa Elena, Ecuador. By using the computational tool HOMER, the authors determined that the optimal HES configuration for the studied location comprised a PV array, wind turbines, a diesel generator, and a battery for electric storage.
The authors of [33] presented a PSO algorithm for the optimal sizing of a hybrid standalone power system and used empirical data from Ardebil, a province in Northwest Iran, to test the optimization model. The study proposed an approximate method for reliability evaluation and included a comprehensive reliability assessment of a hybrid wind/photovoltaic/fuel cell generation system. In [34], the authors proposed a linear programming formulation for the optimization of an autonomous energy system located in Lebanon (composed of wind, solar, diesel, and battery storage). In the study, various control policies and operation modes of diesel engines and batteries were evaluated. The authors of [35] simulated a PV/wind HES with battery storage in ARENA 12.0 software and compared its results to an optimization method using an SA algorithm. The study focused on an HES located in Turkey and provided a systematic approach for the design and evaluation of HESs. Moreover, their study suggested the use of loss of load probability (LLP) as a method to evaluate the autonomy of these types of systems. Besides the methods mentioned above, some other expert systems have been applied to determine the optimal component size, including iHOGA, Hybrid2, and HOMER [36][37][38]. Table 1 presents a summary of optimization techniques along with their area of study.
Although a large number of studies have reported on the feasibility of HESs in various regions of the world (for comprehensive reviews, see [13,39,40]), there is a gap in the literature on the potential of hybrid energy systems in Ecuador, in particular on the application of these types of systems for rural electrification. With this in mind, the purpose of this article is multifold. First, it presents a computable model that allows the determination of the capacity of each system component in an autonomous HES and, simultaneously, find the optimal operating schedule of the considered technologies. Second, on the methodological ground, the proposed model contributes to the literature by extending the linear problem formulation given in [27] and makes use of recent advances in the domain of subsystem component modeling described in [41][42][43][44]. Third, this paper extends the available literature related to the analysis of load profiles for remote areas. It considers the case of a rural community and formulates a load profile given the number of electrical appliances and usage habits of a low-income Ecuadorian household. Fourth, the proposed model is applied to investigate the potential of HESs to power off-grid, low-income communities in two regions of Ecuador, which is a relatively unexplored topic (in terms of innovative means for rural electrification).
The reminder of the paper is organized as follows: In Section 2, the problem definition is presented. In Section 3, the formulation of the computable model for optimal sizing and scheduling of autonomous HES is provided. The case studies and input data for the model are described in Section 4. Section 5 summarizes the model results of the two case studies. Finally, conclusions are presented in Section 6.

Problem Statement
As discussed in Section 1, this paper contributes to the existing literature by providing a framework capable of identifying the optimal component size and determining the optimal generation schedule of autonomous HESs. This is achieved through the formulation of an optimization model based on a linear programming (LP) approach. The objective function is defined as the minimization of the total life cycle cost (TLCC) of the HES, taking into account the technical constraints associated with each subsystem. Moreover, the model considers local meteorological and economic parameters. The model addresses the problem of the techno-economic feasibility of HES projects in isolated areas of Ecuador. To solve the abovementioned problem, the following simplifying assumptions are adopted: 1.
The hybrid energy systems are assumed to work independently and are expected to power isolated rural communities by supplying electricity in an efficient, reliable, and economical way.

2.
The optimization procedure is conducted over the period of one year with an hourly resolution.

3.
The meteorological conditions and electricity consumption values are kept constant within each one-hour time step. 4.
The hybrid energy systems may be comprised of solar-, wind-, battery-, and diesel-powered subsystems. Figure 3 shows a sample layout of a potential HES configuration. 5.
The annual meteorological conditions, load profiles, maintenance, and capital costs are considered to be representative of each analyzed case. This allows to evaluate the economic feasibility of the project through the computation of the TLCC (total life cycle cost) of the system. A generic example of the energy usage of a typical household and a daily solar energy production profile are presented in Figure 4, illustrating the mismatch between energy supply and demand in an HES.
Energies 2020, 13, 3933 6 of 20 6. The generation of renewable energy technologies is based on hourly capacity factors. This is calculated considering the specific technical properties of each technology as well as hourly meteorological conditions. subsystems. Figure 3 shows a sample layout of a potential HES configuration. 5. The annual meteorological conditions, load profiles, maintenance, and capital costs are considered to be representative of each analyzed case. This allows to evaluate the economic feasibility of the project through the computation of the TLCC (total life cycle cost) of the system. A generic example of the energy usage of a typical household and a daily solar energy production profile are presented in Figure 4, illustrating the mismatch between energy supply and demand in an HES. 6. The generation of renewable energy technologies is based on hourly capacity factors. This is calculated considering the specific technical properties of each technology as well as hourly meteorological conditions.

Model Formulation
This section presents the mathematical model for the optimization of autonomous hybrid energy systems. The proposed linear programming model is an extension of the one developed by [27], and 5. The annual meteorological conditions, load profiles, maintenance, and capital costs are considered to be representative of each analyzed case. This allows to evaluate the economic feasibility of the project through the computation of the TLCC (total life cycle cost) of the system. A generic example of the energy usage of a typical household and a daily solar energy production profile are presented in Figure 4, illustrating the mismatch between energy supply and demand in an HES. 6. The generation of renewable energy technologies is based on hourly capacity factors. This is calculated considering the specific technical properties of each technology as well as hourly meteorological conditions.

Model Formulation
This section presents the mathematical model for the optimization of autonomous hybrid energy systems. The proposed linear programming model is an extension of the one developed by [27], and

Model Formulation
This section presents the mathematical model for the optimization of autonomous hybrid energy systems. The proposed linear programming model is an extension of the one developed by [27], and it is built upon the formulations presented in [41][42][43][44]. A new feature in our model is the direct inclusion of the total life cycle cost (TLCC) incurred by each technological component in the objective function. With these new cost elements, the model can be used to obtain the minimum levelized cost of energy (LCOE) and define the optimal capacity of each subsystem and their optimal operation schedule. Similarly to [41], in this study, we assume that the annual electricity generated and the annual costs of production remain constant over the system's lifetime. With this assumption and the calculation of a capital recovery factor (CRF), it is possible to allocate the results obtained from the optimization Energies 2020, 13, 3933 7 of 20 procedure for one year to the entire lifetime of the hybrid energy system. The CRF is determined by [45].
where i is the interest rate, and A is the project lifetime.
The objective function of the model is presented in Equation (2). It minimizes the total life cycle costs of the HES. In this equation, TLCC is calculated as the sum of life cycle costs (in USD) of each power-generating technology (TLCC t )plus the life cycle costs of the batteries (TLCC B ).
The total life cycle costs of the photovoltaic (TLCC PV ) and wind (TLCC W ) power generation subsystems are determined by Equations (3) and (4). The TLCCs of these technologies include capital costs (IC PV and IC W in $/kW) and annual operation and maintenance costs (MC PV and MC W ; these costs are estimated as a fixed percentage of the capital costs (%)). The capital costs of the DC-to-AC inverter and other electrical and electronic equipment, as well as the inverter replacement costs, are included in IC PV [41]. The abovementioned costs are multiplied by the variables representing the installed capacity of each specific technology (C PV and C W in kW).
The total life cycle cost of the diesel generator subsystem (TLCC D ) is given in Equation (5). Analogous to the TLCC of the solar and wind-powered subsystems, TLCC D includes the capital costs (IC D ) in $/kW, maintenance costs (MC D ) as a percentage of the capital costs (%), and variable production costs [27,41]. Similar to other works [46], the variable production costs are determined by the diesel generator power output p D in kW, diesel price (DP) in $/kWh , and the efficiency of the diesel generator (η D ).
The total life cycle cost of the battery subsystem (TLCC B ) is specified by the capital costs (IC B ), maintenance costs (MC B ), and operating costs. The total operating costs are calculated by considering the annual operating cost of the battery (OC B , in $/kWh) and the sum of energy charged (b bc,h ) and discharged (b dc,h ) during the year. Moreover, it is assumed that the battery bank will need to be replaced once throughout the system's lifetime (after 10 years); hence, the TLCC B includes a one-time replacement cost (RC B ) which is equivalent to the capital cost (IC B ) discounted to its present monetary value [47]. TLCC B is defined using Equation (6) [41]: The minimization of the TLCC is subject to constraints (7)- (13). Equation (7) is the energy balance which states that the energy produce d by each technology p t,h and the energy flow from (b dc,h ) and to b dc,h the battery must satisfy the hourly energy demand (D h ): Equation (8) restricts the power output of each technology t. This limit is defined as the product of the installed capacity C t and the capacity factor CF t,h : Energies 2020, 13, 3933 8 of 20 The capacity factor CF t,h used in this study is defined as the fraction of the installed capacity of each technology t that is available at the time segment h [27]. Since the operation of the diesel generator is independent of the meteorological factors, and technical interruptions are not taken into consideration, the capacity factor for this technology CF D,h is set to 1 for every time segment t.
For photovoltaic and wind-powered technologies, the capacity factors are calculated beforehand according to (9) and (10). The capacity factor of the PV subsystem CF PV,h is calculated as the ratio of the solar irradiation incident on the PV array (I h ) to the incident irradiation at standard test conditions (I STC ) for each time segment h [42,43]. This is given by Equation (9). The effect of temperature on the PV array is neglected, and it is assumed that for a typical flat-plate module the incident irradiation at standard conditions is 1 kW/m 2 [48].
The wind subsystem capacity factor CF W,h is based on the simplifying assumption proposed by [44] and is calculated as follows: where v h is the wind speed in m/s, v in is the cut-in speed, v r is the rated speed, and v o f f is the cut-out speed.
Equations (11) and (12) relate the capacity of the electrical storage subsystem (C B ) to the capacity factor of this technology, which in this case is the depth of discharge δ. The depth of discharge is defined as the ratio of the discharged energy in kWh to the usable capacity of the battery (kWh), usually ranging 80-100% for Li-ion batteries [49]. It is modeled to reduce the degradation of the battery by increasing the number of available life cycles. The restriction on the energy flow from the battery (b dc,h ) is given by Equation (11), and the restriction on the energy flow to the battery b ch,h is given by Equation (12) [27]. The battery charge and discharge efficiencies are considered in the model and allow the battery round-trip efficiency to be determined.
The storage level at the time segment BS h+1 is equal to the energy stored at the beginning of the previous hour (BS h ) plus the net energy charged (b ch,h ) minus the energy released (b dc,h ) to cover the energy demand.
To compare the competitiveness of the analyzed HESs to other alternative solutions, the levelized cost of energy (LCOE) is calculated using Equation (14). The LCOE is usually referred to as the price at which electricity must be sold to break even [50].
A flowchart of the methodological framework proposed in this paper is provided in Figure 5. It shows the importance of the construction of the mathematical equations and the input data (economic, meteorological, and technical data of the system) on the modeling and design of hybrid energy systems. These two elements are essential for the optimization process. The mathematical model was implemented in General Algebraic Modeling System (GAMS) 27.1.0 and solved using the Energies 2020, 13, 3933 9 of 20 CPLEX 12.9.0.0 solver [51]. The model's output was exported to a spreadsheet software, and MATLAB 2019a was employed for data structuring, postprocessing, and plotting. Finally, the numerical and graphical results of the optimization process were compiled and analyzed.

=
• ∑ A flowchart of the methodological framework proposed in this paper is provided in Figure 5. It shows the importance of the construction of the mathematical equations and the input data (economic, meteorological, and technical data of the system) on the modeling and design of hybrid energy systems. These two elements are essential for the optimization process. The mathematical model was implemented in General Algebraic Modeling System (GAMS) 27.1.0 and solved using the CPLEX 12.9.0.0 solver [51]. The model's output was exported to a spreadsheet software, and MATLAB 2019a was employed for data structuring, postprocessing, and plotting. Finally, the numerical and graphical results of the optimization process were compiled and analyzed.

Selected Sites in Ecuador
Compared to other countries of Latin America, there is a limited body of literature that explores the potential as well as the implementation of HESs in Ecuador [52]. Therefore, the proposed optimization framework is used to evaluate the techno-economic feasibility of autonomous hybrid energy systems located in two isolated areas of this country, one in the Amazonian tropical forest and one on a subtropical island. The selected areas are in the provinces of Morona Santiago and Galapagos. While the former has one of the lowest national electrification rates and is inhabited by small, isolated, dispersed indigenous communities [53], the latter has only recently introduced local environmental policies aimed at promoting the implementation of an HES based on renewables, seeking to achieve a zero-emission fossil fuels objective [54]. Grid extension in these two provinces has been insufficient due to the difficult geographical conditions, inadequate infrastructure, and high investment costs. Nonetheless, in recent years, the Ecuadorian government, in coordination with international agencies, has proposed the use of distributed renewable energy technologies as mechanisms to expand the electricity access to isolated communities in the Andean and Amazonian regions [55].
Moreover, in the last decade, Ecuador has experienced profound changes in its energy sector. The country is transitioning from a fossil-fuel-based energy system to a sustainable system where renewable electricity sources account for approximately 60% of the total installed capacity. The Ecuadorian government invested nearly USD 6 billion in eight different projects comprised of runof-river and flexible reservoir systems [56][57][58]. While hydropower constitutes the bulk of Ecuador's installed capacity, resources like solar and wind have had little penetration in the Ecuadorian energy mix, mainly due to policy barriers and limited funding for emerging technologies. Despite the few solar and wind energy systems in the country, Ecuador has a high solar resource potential, with an average solar insolation of 4575 Wh/m 2 /day, and an available wind potential of 1.6 GW [58]. Hence, the deployment of renewable energy systems may be an effective instrument to increase the electrification of isolated regions in the country. Figure 6 shows the (a) solar and (b) wind potential maps of Ecuador. In this study, the meteorological data, including global horizontal irradiation (GHI) and wind speed of the two selected sites were obtained from the National Solar Radiation Database (NSRDB) [59]. The data sets are for one full year with an hourly resolution.

Selected Sites in Ecuador
Compared to other countries of Latin America, there is a limited body of literature that explores the potential as well as the implementation of HESs in Ecuador [52]. Therefore, the proposed optimization framework is used to evaluate the techno-economic feasibility of autonomous hybrid energy systems located in two isolated areas of this country, one in the Amazonian tropical forest and one on a subtropical island. The selected areas are in the provinces of Morona Santiago and Galapagos. While the former has one of the lowest national electrification rates and is inhabited by small, isolated, dispersed indigenous communities [53], the latter has only recently introduced local environmental policies aimed at promoting the implementation of an HES based on renewables, seeking to achieve a zero-emission fossil fuels objective [54]. Grid extension in these two provinces has been insufficient due to the difficult geographical conditions, inadequate infrastructure, and high investment costs. Nonetheless, in recent years, the Ecuadorian government, in coordination with international agencies, has proposed the use of distributed renewable energy technologies as mechanisms to expand the electricity access to isolated communities in the Andean and Amazonian regions [55].
Moreover, in the last decade, Ecuador has experienced profound changes in its energy sector. The country is transitioning from a fossil-fuel-based energy system to a sustainable system where renewable electricity sources account for approximately 60% of the total installed capacity. The Ecuadorian government invested nearly USD 6 billion in eight different projects comprised of run-of-river and flexible reservoir systems [56][57][58]. While hydropower constitutes the bulk of Ecuador's installed capacity, resources like solar and wind have had little penetration in the Ecuadorian energy mix, mainly due to policy barriers and limited funding for emerging technologies. Despite the few solar and wind energy systems in the country, Ecuador has a high solar resource potential, with an average solar insolation of 4575 Wh/m 2 /day, and an available wind potential of 1.6 GW [58]. Hence, the deployment of renewable energy systems may be an effective instrument to increase the electrification of isolated regions in the country. Figure 6 shows the (a) solar and (b) wind potential maps of Ecuador. In this study, the meteorological data, including global horizontal irradiation (GHI) and wind speed of the two selected sites were obtained from the National Solar Radiation Database (NSRDB) [59]. The data sets are for one full year with an hourly resolution.
The geographical locations of the selected areas in the provinces of Morona Santiago and the Galapagos Islands are characterized by high solar resource potential [56] and therefore suitable for the deployment of off-grid photovoltaic systems. In terms of wind power conditions, the selected site in Morona Santiago presents significantly lower mean wind speeds than those recorded in the Galapagos Islands. This is due to the contrasting climatic regions of Ecuador and of the two analyzed areas. The coastal and marine environment in the Galapagos Islands favors the deployment of small, low-wind-speed turbines [62], whereas, in the province of Morona Santiago, the low wind velocities make the application of wind-powered technologies less effective. The average monthly wind speeds and global horizontal irradiation of the considered site in the Galapagos province are shown in Figure 7. The geographical locations of the selected areas in the provinces of Morona Santiago and the Galapagos Islands are characterized by high solar resource potential [56] and therefore suitable for the deployment of off-grid photovoltaic systems. In terms of wind power conditions, the selected site in Morona Santiago presents significantly lower mean wind speeds than those recorded in the Galapagos Islands. This is due to the contrasting climatic regions of Ecuador and of the two analyzed areas. The coastal and marine environment in the Galapagos Islands favors the deployment of small, low-wind-speed turbines [62], whereas, in the province of Morona Santiago, the low wind velocities make the application of wind-powered technologies less effective. The average monthly wind speeds and global horizontal irradiation of the considered site in the Galapagos province are shown in Figure  7. Numerous works have studied the energy use and daily demand patterns of rural households across Latin America [63]. For both case studies, the daily load profile of a low-income rural consumer is generated with the use of LoadProGen, a novel MATLAB tool that uses a bottom-up stochastic approach for the development of load profiles of off-grid systems for rural electrification [64]. Information on the number of electrical appliances and usage habits-including appliance functioning times and functioning windows-of low-income households living in isolated communities in Ecuador were obtained from [65,66] to compute the daily load profile. Furthermore, based on [66], it is assumed that the unelectrified community is composed of 30 houses, with a typical  The geographical locations of the selected areas in the provinces of Morona Santiago and the Galapagos Islands are characterized by high solar resource potential [56] and therefore suitable for the deployment of off-grid photovoltaic systems. In terms of wind power conditions, the selected site in Morona Santiago presents significantly lower mean wind speeds than those recorded in the Galapagos Islands. This is due to the contrasting climatic regions of Ecuador and of the two analyzed areas. The coastal and marine environment in the Galapagos Islands favors the deployment of small, low-wind-speed turbines [62], whereas, in the province of Morona Santiago, the low wind velocities make the application of wind-powered technologies less effective. The average monthly wind speeds and global horizontal irradiation of the considered site in the Galapagos province are shown in Figure  7. Numerous works have studied the energy use and daily demand patterns of rural households across Latin America [63]. For both case studies, the daily load profile of a low-income rural consumer is generated with the use of LoadProGen, a novel MATLAB tool that uses a bottom-up stochastic approach for the development of load profiles of off-grid systems for rural electrification [64]. Information on the number of electrical appliances and usage habits-including appliance functioning times and functioning windows-of low-income households living in isolated communities in Ecuador were obtained from [65,66] to compute the daily load profile. Furthermore, based on [66], it is assumed that the unelectrified community is composed of 30 houses, with a typical Numerous works have studied the energy use and daily demand patterns of rural households across Latin America [63]. For both case studies, the daily load profile of a low-income rural consumer is generated with the use of LoadProGen, a novel MATLAB tool that uses a bottom-up stochastic approach for the development of load profiles of off-grid systems for rural electrification [64]. Information on the number of electrical appliances and usage habits-including appliance functioning times and functioning windows-of low-income households living in isolated communities in Ecuador were obtained from [65,66] to compute the daily load profile. Furthermore, based on [66], it is assumed that the unelectrified community is composed of 30 houses, with a typical household consisting of five people. Table 2 summarizes the number of appliances and their average total daily time of use.
The techno-economic parameters of the constituent technologies are presented in Table 3, and the technological specifications and assumptions used in the model are shown in Table 4.  Table 4. Technological specifications used in the model.

Results
As previously mentioned, in this paper, a mathematical programming model was employed to determine the optimal configuration and generation schedule of two autonomous hybrid energy systems. The systems were assumed to serve over a lifetime of 20 years, with an 8.6% annual interest rate. Two case studies were carried out using data from sites located in the provinces of (i) Morona Santiago and (ii) Galapagos, Ecuador. System configurations such as diesel, diesel-PV, diesel-wind, diesel-battery, diesel-wind-battery, diesel-PV-battery, PV-battery, PV-wind-battery, and diesel-PV-wind-battery were assessed for each site. The subsections below show the results of the model computed for both cases. Table 5 summarizes the outputs of the model for the hybrid energy project located in the Galapagos province. The results show that a diesel-PV-wind-battery system leads to the lowest value of LCOE (0.30 $/kWh) and a total life cycle cost of approximately $36,352. The system consists of a 1.27 kW diesel generator, a 5.46 kW PV array, a 2.84 kW wind turbine, and a 29.72 kWh battery. In this system, renewable energy sources contribute nearly 95% (63% PV and 32% wind turbine) of the total power generation, while the diesel generator meets the remaining 5%. It is worth noting that in 2007, the Ecuadorian government launched the "Zero Fossil Fuels" initiative in the Galapagos Islands [54], a program aimed at decreasing the use of fossil fuels in electricity generation. Therefore, even though a PV-wind-battery configuration does not lead to the lowest LCOE (0.36 $/kWh) due to its high investment costs, it is reasonable to assume that this configuration should be considered as the optimal case for the Galapagos Islands. Table 5 also shows the annual CO 2 emissions from the combustion of diesel fuel for all the configurations analyzed. As can be seen from the table, the hybridization of energy sources has a substantial impact on the TLCC and their respective emissions. In comparison to Energies 2020, 13, 3933 12 of 20 the diesel-only system, the diesel-PV-wind-battery configuration emitted nearly 13 times less CO 2 . Figure 8 shows the component capacities of the different system configurations. combination of a wind turbine with a diesel generator reduces the LCOE and TLCC by nearly 33%. For this configuration, 84% of electricity is generated by the wind subsystem, and the diesel generator serves as a secondary source of electricity.   The results also reveal that a system configuration comprised of a standalone diesel or a diesel-battery system leads to an LCOE of 0.69 $/kWh. For both cases, a 6 kW diesel generator would be required to meet the load demand of a typical low-income unelectrified village. However, the combination of a wind turbine with a diesel generator reduces the LCOE and TLCC by nearly 33%. For this configuration, 84% of electricity is generated by the wind subsystem, and the diesel generator serves as a secondary source of electricity. Figure 9 presents the estimates of the TLCCs for each system configuration, including capital and O&M costs. In the case of a standalone diesel generator, capital costs represent 2.6% of the TLCC, and the remaining 97.4% is related to fuel and maintenance costs. For the configuration comprised of a diesel generator and PV-wind-battery subsystems, the majority of the capital costs are split between the battery, PV module, and wind turbine (49%, 29.2%, and 19.9%), while the fueling costs of the diesel generator represent the largest share of O&M costs (62.7%). Such results highlight the importance of RESs and their significant impact on the TLCC of a hybrid energy system. and the remaining 97.4% is related to fuel and maintenance costs. For the configuration comprised of a diesel generator and PV-wind-battery subsystems, the majority of the capital costs are split between the battery, PV module, and wind turbine (49%, 29.2%, and 19.9%), while the fueling costs of the diesel generator represent the largest share of O&M costs (62.7%). Such results highlight the importance of RESs and their significant impact on the TLCC of a hybrid energy system. The model proposed in this study is capable of providing an optimal operating schedule for the considered technologies in addition to determining the appropriate capacity of each component. Figure 10 illustrates the optimal power generation of the Diesel-PV-Wind-Battery (D-P-W-B) system configuration for 24 h. It can be observed that the solar-powered subsystem can generate electricity between 6 AM and 6 PM, with its highest output in the middle of this period. Most of the surplus energy is stored in the battery, which is later discharged and backed up by the diesel generator. The figure shows the battery storage depth of discharge, which, in the case of this study, was set to be 90%. The model proposed in this study is capable of providing an optimal operating schedule for the considered technologies in addition to determining the appropriate capacity of each component. Figure 10 illustrates the optimal power generation of the Diesel-PV-Wind-Battery (D-P-W-B) system configuration for 24 h. It can be observed that the solar-powered subsystem can generate electricity between 6 AM and 6 PM, with its highest output in the middle of this period. Most of the surplus energy is stored in the battery, which is later discharged and backed up by the diesel generator. The figure shows the battery storage depth of discharge, which, in the case of this study, was set to be 90%. between the battery, PV module, and wind turbine (49%, 29.2%, and 19.9%), while the fueling costs of the diesel generator represent the largest share of O&M costs (62.7%). Such results highlight the importance of RESs and their significant impact on the TLCC of a hybrid energy system. The model proposed in this study is capable of providing an optimal operating schedule for the considered technologies in addition to determining the appropriate capacity of each component. Figure 10 illustrates the optimal power generation of the Diesel-PV-Wind-Battery (D-P-W-B) system configuration for 24 h. It can be observed that the solar-powered subsystem can generate electricity between 6 AM and 6 PM, with its highest output in the middle of this period. Most of the surplus energy is stored in the battery, which is later discharged and backed up by the diesel generator. The figure shows the battery storage depth of discharge, which, in the case of this study, was set to be 90%.

Morona Santiago Province
The meteorological conditions in the province of Morona Santiago differ substantially from those in the Galapagos Islands. While the annual average wind speed in the province of Galapagos is around 4 m/s, the annual average wind speed in the Amazonian tropical forest is approximately 1 m/s. The results in Table 6 and Figure 11 show that the significantly lower wind speeds in the province of Morona Santiago excluded the wind power sources from the optimal HES configuration.
In the system configuration powered only by a diesel generator, the LCOE (0.69 $/kWh) and TLCC ($84,221.42) are roughly the same as in the case of the Galapagos Islands. However, the optimization results show that the usage of PV and battery in the optimal design of the HES leads to a lower LCOE. A diesel-PV-battery system with a TLCC of $45,035.31 and comprised of a 0.8 kW diesel generator, an 11.38 kW PV array, and a 39.95 kWh battery results in the lowest LCOE, roughly 0.37 $/kWh. This is 15.6% higher than the LCOE of the diesel-PV-battery system in the Galapagos Islands. In the province of Morona Santiago, the solar irradiation is lower than in the Galapagos Islands, so the evident cost difference is due to the higher PV capacity of the system.   Figure 12 presents the capital and O&M cost of each system configuration. Fuel expenses constitute the majority of the O&M costs in most cases, yet the O&M costs drastically decrease for system configurations with solar power and battery storage subsystems. On the other hand, battery storage has a relatively high investment cost, and it accounts for nearly half of the TLCC in the HES powered by PV. Analogous to the case of the Galapagos Islands, the hybridization of energy sources resulted in a significant reduction of CO 2 emissions. Although the emissions of a diesel-PV-wind-battery standalone power system in Morona Santiago are slightly higher than a similar system in the Galapagos Islands (because of the low wind speeds in the area), the results presented in Table 6 show that the application of these technologies in Morona Santiago can help reduce the reliance of remote communities on diesel for power generation. Figure 12 presents the capital and O&M cost of each system configuration. Fuel expenses constitute the majority of the O&M costs in most cases, yet the O&M costs drastically decrease for system configurations with solar power and battery storage subsystems. On the other hand, battery storage has a relatively high investment cost, and it accounts for nearly half of the TLCC in the HES powered by PV. Figure 12 presents the capital and O&M cost of each system configuration. Fuel expenses constitute the majority of the O&M costs in most cases, yet the O&M costs drastically decrease for system configurations with solar power and battery storage subsystems. On the other hand, battery storage has a relatively high investment cost, and it accounts for nearly half of the TLCC in the HES powered by PV.

Conclusions
This paper develops a techno-economic model based on a linear programming approach to optimize autonomous hybrid energy systems. The model has been implemented in the modeling language GAMS, and it is employed to determine the optimal component size and generation schedule of two independent HESs. The model considers investment, O&M, replacement, and fueling costs. Technical data such as component efficiencies, rated irradiation, rated wind speed, cut-in wind speed, and cut-out wind speed, as well as meteorological data (wind speed and irradiation), are used as input parameters. The model evaluation demonstrated that meteorological conditions are crucial during the design phase of an off-grid HES, but economic parameters (fuel prices and investment costs, among others) also have a significant impact on the optimal solution.

Conclusions
This paper develops a techno-economic model based on a linear programming approach to optimize autonomous hybrid energy systems. The model has been implemented in the modeling language GAMS, and it is employed to determine the optimal component size and generation schedule of two independent HESs. The model considers investment, O&M, replacement, and fueling costs. Technical data such as component efficiencies, rated irradiation, rated wind speed, cut-in wind speed, and cut-out wind speed, as well as meteorological data (wind speed and irradiation), are used as input parameters. The model evaluation demonstrated that meteorological conditions are crucial during the design phase of an off-grid HES, but economic parameters (fuel prices and investment costs, among others) also have a significant impact on the optimal solution.
The model was applied to assess two different hybrid energy projects located in isolated areas of Ecuador, one in the Amazonian tropical forest and one on a subtropical island. The energy projects are expected to supply power to low-income households living in rural communities. Considering the economic and environmental aspects of the HES project located at the Galapagos province, the most cost-efficient system is comprised of a 10.51 kW PV array, 1.52 kW wind turbine, and 49.18 kWh battery. The computations showed that this system capable of providing electricity at an LCOE of approximately 0.36 $/kWh and a TLCC of $44,153. Although the entirety of the electricity is produced by RES, the lithium-ion battery storage provides flexibility and continuous service in times of unfavorable meteorological conditions. For the second case study, the HES project located in the province of Morona Santiago, the optimization model excluded the wind-powered subsystem due to the insufficient wind speeds. The optimized system configuration is comprised of a 0.8 kW diesel generator, 11.38 kW PV array, and 39.95 kWh battery; this configuration has an LCOE of 0.37 $/kWh and a TLCC of $45,035.
The results obtained from the application of the mathematical model developed in this study are in accordance with other sizing and optimization techniques commonly used for the feasibility analysis of renewable-energy-based hybrid systems. For example, [68,78,79] reported LCOE estimates ranging from 0.196 to 0.410 $/kWh. Furthermore, these studies also showed that the hybridization of energy sources in autonomous systems provides a substantial reduction of the total annualized costs and CO 2 emissions.
To conclude, the case studies show that the proposed model adequately considers technical and economic constraints while minimizing the life cycle cost of the HES. Such a model is a step towards the development of computable models to foster global electrification, inclusive communities, and sustainable development. Currently, it constitutes an efficient and valuable decision-making tool in the sizing and scheduling of HESs.
Our results also demonstrate that the model presented in this study can be used as a building block for other optimization methods and simulation tools. For future research, the optimization model presented in this research will be extended to consider other techno-economic factors such as payment and tariff methods, taxes, distribution charges, and energy conservation methods. Moreover, the computable model will be compared with other optimization tools.
Funding: This research received no external funding.
Acknowledgments: This study was conducted within the statutory research of the Mineral and Energy Economy Research Institute of the Polish Academy of Sciences.

Conflicts of Interest:
The authors declare no conflict of interest.