Multi-criteria planning of microgrids for rural electrification

Aim: The objective of this study is to facilitate decision makers’ selection of groups of components for isolated microgrid which should not only be economically attractive and efficient but also acceptable from various aspects. The microgrid is designed to serve a rural remote region of Tanzania with an approximate energy demand of 1000 kWh/day. Methods: A comprehensive comparison has been presented using (Hybrid Optimization Model for Electric Renewable) HOMER Pro software and fuzzy analytic hierarchy process (AHP)-based multi-criteria decisionmaking (MCDM) (fuzzy-AHP-MCDM) technique for planning of isolated microgrid systems which is formed by using a group of sources, like photo-voltaic (PV) generator, wind-turbine, micro-hydro, bio-gas based generation and battery storage system. In this study, ten various alternatives comprising a blend of the aforementioned sources have been considered. A group of households having an approximate load demand of 1000 kWh/day located near Arusha in Tanzania has been selected for this study. Results: A levelized cost of energy of 0.0694 $/kWh has been found for the optimized case using HOMER Pro comprising all five aforesaid components. Conclusion: The fuzzy-AHP-MCDM technique also exhibits highest priority for a similar combination when only the economic criteria are considered and all other criteria are being suppressed significantly. On the other side, a different group consisting only PV and wind-based generation is found prioritized when all criteria are taken Page 121 Bohra et al . J Smart Environ Green Comput 2021;1:120-34 I http://dx.doi.org/10.20517/jsegc.2021.06 into account with different weights, either based on experts’ score or equal importance. Moreover, a sensitivity analysis has also been conducted to observe the impact of alternative selections.


INTRODUCTION
The availability of low-priced electricity has completely transformed the modern world. And this has been achieved along with the use of petroleum products to form the energy backbone of modern life. According to the World Bank's World Development Indicators, about 85% of the world's households now have access to electricity, and the rest 15% belong to the world's absolutely poorest households [1] . Access to high grade electrical energy should be an absolute priority. Modern energy is now an essential requirement and none of the Sustainable Development Goals will be accomplished if this aspect is not complied with. Sub-Saharan Africa (SSA) has the lowest energy access rates in the world [2] . The electricity reaches only about half of its people, while clean cooking only to one-third; and roughly 600 million people lack electricity [3,4] . Thirteen countries in SSA have less than 25% access, compared to only one in developing Asia [5] .
A few years ago, only economical parameters were of prime importance for power system planners and operators. But now decision makers have realized that the planning should require not only financial aspects but also other facets of the project as far as facilitating electricity to underprivileged households is concerned [6] . This includes nearby resources of fuel and at reasonable price, local skills required to keep running the installed system, use of renewable energy sources to meet carbon emission norms, autonomy, self-supply, distance from load-center, acceptable noise level, and so on [7] .
An attempt was made to review quantitative and qualitative electricity planning and related implementation aspects by considering 49 SSA countries including four regional power pools and the sub-continent as a whole [8] . The formalized multi-criteria decision making (MCDM) tools have been employed in analysis. A microgrid is defined as a smaller power system employing one or more conventional and/or renewable distributed energy sources based generating units. This may be grid supportive or independently operating to feed for remote areas or villages where installation of transmission lines is economically not justifiable. Moreover, microgrid planning is a complex problem which involves technical, economic, social, and environmental considerations. Conventionally, the design and planning of microgrid targets best utilization of energy sources, minimum overall cost, better ancillary services, and resiliency and reliability. In order to plan a microgrid for a given location, approximate load demand should be known beforehand. Such assessment has been done for Tanzania based on interview and measurement, though it does not reflect the true load demand due to the lack of instruments and methodology adopted [9] . One of the IET's reports entitled "Power Network Joint Vision" elaborates on significant changes required in view of current practice in the assessment, planning and development of the power network facilities for future system operation. The report also discussed key hurdles, their solutions, planning techniques [10] . Similar, study was conducted in a scenario of power system planning under a high-penetration of solar PV [11] . The microgrid being emerging technology in the field of power system, the various technical aspects, key drivers and challenges are investigated [12] . The optimal output planning for each individual unit of microgrid comprising wind power, solar PV, gas turbines, and energy storage has been modeled using optimization algorithm with an objective of significant reduction in system operation cost by Wang et al. [13] . A microgrid planning model ensuring adequate supply to critical loads under the uncertain condition is discussed in Wu et al. [14] with minimizing the total cost during the whole project life, construction cost, and operation cost for the distributed energy resources. The overview of microgrid and its associated challenges, like planning, operation, and control, as well as protection aspects are briefly discussed in [15] . It also helps planners and decision makers involved in setting up microgrids for rural households/villages becoming electrified. The concept of optimal design, planning, techno-economic analysis, sizing, and operation of a hybrid renewable energy-based grid-connected microgrid is explained in Hafez et al. [16] with an objective to minimize the lifecycle cost. A comparison and evaluation are done mainly for economics, operational performance, and environmental emissions by using the HOMER software. An extensive survey of mathematical methods applied to microgrid planning based on economic feasibility from a capacity sizing standpoint has been conducted in two prior studies [17,18] and can be useful for microgrid planners. MCDM is a methodology of operational research discipline to evaluate projects involving multiple and mutually conflicting features. Numbers of MCDM techniques are advocated suitable for providing optimal solutions to a various class of energy management and planning problems [19][20][21] . An analytic hierarchy process (AHP)-entropy method has been introduced for microgrid planning considering economy, reliability, basic structure of planning and environmental protection by He et al. [22] .
The applications of HOMER software in optimal planning of hybrid renewable energy sources have been surveyed by Bahramara et al. [23] for wide range of electrical load and stand-alone and grid-connected mode as well. On the similar line of work, the mix of diesel and renewable energy are found in [16] to have minimum net present value and less carbon footprint while pure renewable exhibits almost zero carbon emission but higher net present value for grid-connected microgrid analysis. The technical and economic analyses of hybrid renewable energy systems are considered to be essential ones for the efficient utilization of renewable energy resources for microgrid applications [24] . The optimum planning of renewable energy system comprising hybrid PV-Wind-BESS, which has been conducted for a small village in India having peak demand of approximately 27 kW, has been investigated by Fulzele and Dutt [25] . Various important aspects addressing the challenges for rural microgrid design with an objective of optimal investment and planning using different techniques and tool are discussed by Heleno et at. [26] . Framework for microgrid has been proposed using multidisciplinary design optimization under uncertainty approach [27] . Jiménez-Estévez et al. [28] also incorporated life-cycle-assessment module, so that the configuration can take care of economic and environmental impacts of the power system. The planning and design of isolated microgrids is discussed by Vafaei et al. [29] focusing on their technical and investment costs. In the same work, they have employed traditional optimization technique for selecting sources and power generation technologies and storage devices to ensure minimal operating cost for microgrid. Logenthiran et al. [30] has explained differences between classical integer minimization problem and evolutionary strategy technique for sizing power equipment for isolated microgrid with an objective of minimizing capital, operational, and maintenance cost of energy sources. The applications of heuristic optimization techniques such as Genetic Algorithm, Particles Swarm Optimization, Simulated Annealing, Ant Colony (AC), and Artificial Immune Systems for renewable energy systems sizing have been discussed by Erdinc [31] . Chen et al. [32] proposed a method based on the cost-benefit analysis for optimal sizing of an energy storage system in a microgrid. A designing of stand-alone microgrid to achieve optimum cost and sizing of storage for the Island of Singapore, and selection and capacity sizing of generators for a rural microgrid have been explained by Fan et al. [33] and Kumaravel et al. [34] , respectively and by other researchers as well.
However, though economic factors are greatly influential in deciding the type and size of sources of microgrid, their significance in terms of types and capacity as well may change if other non-tangible factors are given appropriate priority and this fact is being highlighted in this paper. Moreover, the Cost of Energy (CoE) is one of the most sensitive and governing factors to decide the acceptance of a project by local consumers which has been looked at in this study. The aim of this study is to facilitate decision makers to achieve economic benefits along with other factors to be satisfied for planning of an isolated microgrid system catering a group of households, but it can be equally applicable to any size of the systems. The remaining of the paper is organized as follows. In section 2, the selected sample reference system is described. Section 3 explains the methodology adopted for analysis. The results are discussed in section 4. Finally, section 5 concludes the paper.

SYSTEM DESCRIPTION
The system under study is located in Tanzania and the objective is to make an electrification plan leveraging the multifold benefits of microgrids. Tanzania is the sixth-most populous country in SSA; it connects six land-locked countries to the Indian Ocean, and with the discovery of natural gas reserves, it opens opportunities for investment. Tanzania has abundant and world-class wind and solar resources [35] . But it is a low-income country, with a population of about 59.6 million in 2020. About 66% of the population lives in rural areas. Tanzania has outlined its medium-term objective of becoming a middle-income country by 2025, and the country also participates in seven regional projects with a total commitment value of $550.8 million. The identified sectors in the national projects include: transport (22%), urban development (20.5%), energy (13.5%), education (10%), and governance (10%) [36] . Arusha is a city in East Tanzania, located at the base of volcanic Mt. Meru. It is a gateway to safari destinations and to Africa's highest peak, Mt. Kilimanjaro, and to the west lie Serengeti National Park, which is home to wildlife.
Tanzania has immense solar and wind power potential, a global horizontal irradiation of 4-7 kWh/m 2 per day, and average wind speed of 8.5-10 m/s. Moreover, biomass is a renewable energy source developed from living or recently living plant and animal materials. Methane rich bio-gas is produced when the biomass is anaerobically degraded by micro-organisms using anaerobic digestion [37] . According to a report by Prof. Siza Tumbo, the director general of the Arusha-based Centre for Agricultural Mechanisation and Rural Technology (Camartec) and implementing the Tanzania Domestic Biogas Programme, Tanzania has the third largest livestock population on the African continent and hence has enormous potential for biogas plants [38] .
To model the examined system, that is, configuring the microgrid and its components, carrying out simulations and determining different energy planning scenarios considering different constraints, the HOMER Pro [39] software is used. The microgrid is comprised of electric load, PV system, group of windturbines, bio-gas generator set (genset), micro-hydro power plant, battery energy storage system (BESS), and associated AC-DC converter as shown in Figure 1. The principal characteristics of the microgrid system are: 25 years project lifetime, 3% expected inflation rate, and 5% nominal discount rate. The various components for simulation study of the project including the number of units and capacity of each unit considered are shown in Table 1.  [39] .
Many combinations could be possible depending on different energy sources and their units. In order to limit these combinations, the system has been simulated into HOMER Pro to identify a group of energy sources that offer the most economic performance. Therefore, out of several combinations of sources integrated into HOMER Pro forming microgrid, the one which has least CoE has been considered as a reference or base system for subsequent analysis. Thus, it infers that this system analysis takes mainly economic constraints into account. The respective economic parameters are now supplied to fuzzy-AHPbased MCDM method. Moreover, these base values are also used for sensitivity analysis to explore the impact of tangible parameters, like CoE and net present cost (NPC). For such projects, the lowest value of CoE acts as one of the most sensitive and deciding parameters for acceptance by households.

METHODOLOGY
The projects aiming improving living standards for underprivileged and/or people living in underdeveloped Asian or African countries should have holistic approach fulfilling multiple needs. And it is mainly measured by quantifiable things, such as income, employment opportunities, cost of goods and services, and poverty. The accessibility of electrical energy should meet certain bare minimum demand, leaving the environment minimally polluted and preserving the natural beauty of the surrounding and habitat. Furthermore, the project should also facilitate generating of employment opportunities for locals' and hence impacts the living standards directly or indirectly. The objective of this study is to highlight the fact that economic indicators alone are not sufficient to decide the overall success of such project. Other aspects also, which seem to be less influential at a first glance, should be given due importance. A techno-economic analysis for an approximate demand of 175 kWh/day by a group of houses located remotely catered by photovoltaic mini-grid systems has been projected by Daniel for Nigeria [44] . A composite nature of a load profile for a typical village-having a population of around 2500-3000 people and therefore demanding an approximate electrical demand of 1000 kWh/day has been considered. The composite load may include lights, fans, other typical domestic loads per household, water purifying system, mini/small industrial loads like cutting, sewing, and grinding machines, and a health care and school premise and its associated load [9] . In order to define a reference system and capacity of components, the HOMER Pro software has been widely used by researchers and planners in most analysis. The information of various sizes of sources with their associated parameters, such as its prevailing capital cost per kW, operation and maintenance cost, lifetime, replacement cost, fuel cost (bio-gas), available potential at chosen location for solar irradiation (PV farm) and wind velocity (wind turbine) and details of power and energy BESS are supplied.
The daily irradiation of Arusha, Tanzania is estimated to be 6.17 kWh/m 2 by the National Renewable Energy Lab database [45] and this was utilized in simulation. The annual average wind velocity for the selected location is 5.35 m/s from NASA Surface meteorology and Solar Energy database measured at 50 m above the surface of the earth [46] . Furthermore, it is estimated that under ideal conditions, 10 kg of biomass can produce 3 sq.mt of biogas. A set of simulations have been run to decide the most optimized system architecture based on CoE offered. In a later part of analysis, a comprehensive list of criteria has been gathered [47,48] as depicted in Figure 2. Mainly five criteria such as economic, operational, structural, technical, and others as well as associated sub-criteria have been identified. The "others" criteria include environmental, appearance, inclination of stakeholders, emission level, and more. The impact of weights Canadian Solar MaxPower CS6X-325P [40] 30 × 1 Storage system Gildemeister 10 kW-40 kWh CELLCUBE® FB 10-40 [41] 5 × 5 Wind turbine AWSHC 5.1kW Wind Turbine [42] 4 × 4 Hydroelectric Natel FreeJet FJ-7A 49kW [43] has been analyzed in three ways: firstly, only economic parameters are given more than 70% importance and leaving all other with equal weights; secondly, all primary criteria are assigned equal weights of 20%; and thirdly, expert opinions for the primary criteria are taken into account. The opinion of experts for set of primary criteria has been collected in the form of scores based on comparison of each other. In order to simplify the process of creating pair-wise comparison matrix with minimal influence of fuzziness involved in assigning score to criteria by experts, a method devised [49] has been incorporated.
A group consisting of 15 members was asked to score five primary criteria such as economical, structure, technical, operation, and others; equal weights by assigning scores to each based on a numerical scale of 0-10 are done. These members are representing experts/stakeholders in respective groups and holding key positions for various vital project planning matters, like financing, designing, operating and maintenance, environment protector and social welfare workers, and end-users of electricity.
The collected opinions are real numbers and not the fuzzy sets. In this method a fuzzy number is constructed from a set of scores assigned by experts. The weight determination technique is utilized to take the distribution of the scores into account when estimating the mode and spreads of a fuzzy number. It is assumed that the criterion has been scored using a scale from 0 to 9. Let s 1 , s 2 ..., s n represent the scores assigned to by various experts to same criterion. Due to preconception and expertise in particular areas, different experts are likely to score differently for the same criterion. Here, the approach developed by Cheng [49] has been adopted to identify the membership function of a fuzzy number depending on s 1 , s 2 ..., s n .
In the estimation process, experts' scores closer to the center are considered to have more importance. The estimation of this center is a weightage average of the s i . The distance between each s i is computed to find relative distance matrix defined as follows: where, d ij = |s i -s j |, which implies that d ij ≠ 0 and d ij = d ji . The average of relative distance i d , for each s i is given by where, p ij is the relative importance of s i compared to s j , which implies that p ij = 1 and p ji = 1/p ij , since, P is derived from comparison distance, and it is absolutely consistent Let w i , the real weight of s i , and 0 ≤ w i ≤ 1. Pw = nw (5) which implies that n is an Eigen value of P and w is the corresponding Eigen vector The importance degree w i serves as the weight to be associated with s i when estimating the mode of the fuzzy number.
After getting a pair-wise comparison matrix of primary criteria, this information is fed into a computing device wherein pair-wise comparison for sub-criteria of each primary criterion is done based on quantitative values either obtained from HOMER Pro or from the methodology described by Bohra et al. [48] . The detailed specifications from manufacturer of the components have been used to achieve realistic results. The normalized pair-wise matrix A norm is obtained using the following: where, a ij represents the importance of the i th criterion relative to the j th one and m indicates number of evaluated criteria/sub-criteria to be compared. Then after, criteria weight vector w, which is k-dimensional column vector is computed taking mean of each row of A norm as per Equation (8), Let the j th entry of the i th criterion is being compared with the l th of the same criterion and assuming values for the same criterion lie in the interval of [I j,max , I j,min ]. The relative importance between criteria is ranked between 1-9; the scale was suggested by Saaty [50] .
In AHP, the consistency of matrices can be measured by consistency index (CI) [51] as given in following: where, λ max is the sum of product of associated weight and sum of respective column of pair-wise comparison. The competency of CI can be validated by determining randomized (consistency) index (RI), the average CI for randomly filled matrices. The consistency ratio (CR) and RI can be calculated as below: where, λ max represents maximum Eigen value.
Based on the various systems installed across the world, ten alternatives have been chosen to form a group of sources as shown in Table 2.
The CoE and NPC for microgrid system are obtained from following: Life time cost in given life time CoE kWh = (12) (13) where, t is time of cash flow and n is total number of periods.

DISCUSSION
Since the microgrid is planned for remote and rural applications, it needs to be independent and not be interconnected with any utility grid system. In such scenario, the system CoE solely depends on various costs involved in setting up the system and its components. The base case system details are shown in Table  3. The dispatch control strategy adopted for simulation is load-following, meaning the operating generator produces only enough power to meet the primary load [39] . Moreover, no conventional fossil fuel is being used, while the renewable fraction for system is absolutely 100%. From HOMER simulation results, the optimal alternative is found to be A10. The other system relevant parameters are also depicted in Table 3. The NPC is computed as $327,231 and an initial investment cost of $262,563 for alternative A10 as shown in Table 2. The other key financial indicators like return of investment (ROI), internal rate of return and simple payback period are 76.6%, 61.1% and around 1 year 8 months, respectively.  The scenario when economic criteria are given the most importance is similar to the system simulated in HOMER Pro wherein other parameters like energy density, capacity utilization factor, modularity, footprint, performance degradation, appearance, autonomy, and others are often being overlooked. In this situation, the weights of each are shown in Figure 3A and all five sources (alternative A10) are part of the system and carries the highest priority followed by A9, A8, A7, A1, and so on as shown in Figure 3B. With various capacities of sources, an optimized case has been identified in the HOMER Pro software, which later has been used as base case for some of inputs for fuzzy-AHP-MCDM technique.   The third part of study is the sensitivity analysis. The base case composition, which has been derived during the first part of simulation analysis done in HOMER Pro software and sensitivity analysis have been carried out with the following marginal variations in sensitivity parameters mentioned below. The results are shown in Figure 4: 1. Wind-turbine cost: ± 10% 2. Biomass price ($/tonne): ± 10% 3. Bio-genset cost: ± 10% 4. Bio-genset lifetime: ± 10% 5. Flow rate for hydropower: -10% 6. Hydro plant cost: ± 10% If the price of wind-turbine deviates from the base one, then it has direct effect on CoE, NPC, and initial investment of project [ Figure 5A]. Similarly, the CoE, NPC, and operating cost get affected by the fluctuation on the bio-mass price [ Figure 5B]. The bio-mass system cost has same influence for windturbine cost [ Figure 5C] and likewise bio-gas system plant life [ Figure 5D]. Since the hydro-power plant incurs the highest capital investment, it will change CoE, NPC, and overall capital cost substantially compared to the others explained before [ Figure 5E].
Finally, as the contribution of hydropower is more in terms of energy sharing, the change in water head/ level influences not only CoE and NPC but also fuel requirement/consumption and operating hours of the bio-gas plant as shown in Figure 5F.

CONCLUSIONS
A methodology for designing and planning of isolated microgrid for a rural area of Tanzania was discussed in this paper. The population of the selected region was estimated to be 2500-3000 people consuming approximately 1000 kWh of energy per day. The ranking of alternatives was decided using HOMER Pro and MCDM technique. In order to support the work, a detailed comparison between the result obtained from HOMER Pro software and fuzzy-AHP-MCDM technique was presented. The AHP-MCDM was treated to be sitting above the results generated by HOMER Pro. In this study, the subjective answers collected from experts in terms of scores were handled by fuzzy approach in preparing pair-wise comparison matrix of key criteria and quantitative or measurable criteria using AHP-MCDM technique. The parameters from optimal case obtained in HOMER Pro were used as reference for further analysis. From a sensitivity analysis, it was identified that the water head of a micro-hydro power plant is more impactful than the variation of other parameters. It affects not only CoE and NPC but also generated kWh and consumption of bio-gas significantly. It infers that from fuzzy-AHP-MCDM analysis, when only economic parameters are made emphatic, suppressing other criteria, it results in the use of all sources together. But contradictorily, when all are given equal importance and also when experts' scores are considered, in both cases it results into amalgamation of only two renewable sources, PV and wind energy-based generation. The result may vary significantly, as it depends strongly on the scores collected from experts which may otherwise will lead to unexpected outcome. Moreover, a sensitivity analysis was carried out to examine the impact of alternative selections on various significant indicators of the project. This analysis provided a framework to check feasibility for other such remote places or villages where electricity is unavailable or laying transmission lines is expensive, but rich in renewable energy source potential. Furthermore, the similar analysis could be extended with due modification for grid connected microgrid planning.