Optimal Selection of Hybrid Renewable Energy System Using Multi-Criteria Decision-Making Algorithms

Hegazy Rezk,Irik Z.Mukhametzyanov,Mujahed Al-Dhaifallahand Hamdy A.Ziedan

1College of Engineering at Wadi Addawaser,Prince Sattam Bin Abdulaziz University,Al-Kharj,11911,KSA

2Department of Electrical Engineering,Faculty of Engineering,Minia University,Minia,61517,Egypt

3Center for Strategic and Interdisciplinary Research,Ufa Federal Research Centre of the Russian Academy of Sciences,Ufa,450054,Republic of Bashkortostan,Russia

4Department of Systems Engineering,King Fahd University of Petroleum&Minerals,Dhahran,31261,KSA

5Department of Electrical Engineering,Faculty of Engineering,Assiut University,Assiut,71518,Egypt

Abstract: Several models of multi-criteria decision-making (MCDM) have identified the optimal alternative electrical energy sources to supply certain load in an isolated region in Al-Minya City,Egypt.The load demand consists of water pumping system with a water desalination unit.Various options containing three different power sources:only DG,PV-B system,and hybrid PV-DG-B, two different sizes of reverse osmosis (RO) units; RO-250 and RO-500,two strategies of energy management;load following(LF)and cycle charging(CC),and two sizes of DG;5 and 10 kW were taken into account.Eight attributes,including operating cost,renewable fraction,initial cost,the cost of energy,excess energy,unmet load,breakeven grid extension distance,and the amount of CO2,were used during the evaluation process.To estimate these parameters,HOMER?software was employed to perform both the simulation and optimization process.Four different weight estimation methods were considered; no priority of criteria, based on a pairwise comparisons matrix of the criteria,CRITIC-method,and entropy-based method.The main findings(output results)confirmed that the optimal option for the case study was hybrid PV-DG-B with the following specification:5 kW DG,RO-500,and load following control strategy.Under this condition, the annual operating cost and initial costs were$5546 and$161022,respectively,whereas the cost of energy was 0.077$/kWh.The excess energy and unmet loads were 40998 and 2371 kWh, respectively.The breakeven grid extension distance and the amount of CO2 were 3.31 km and 5171 kg per year, respectively.Compared with DG only, the amount of CO2 has been sharply reduced by 113939 kg per year.

Keywords: Al-Minya city (Egypt); energy efficiency; multi-criteria decision-making; optimization; renewable energy; reverse osmosis units

1 Introduction

Global warming is one of the greatest challenges experienced by humanity in today’s era.The best way to reduce or eliminate its effects is by limiting CO2emissions.This can be achieved by utilizing renewable energy sources (RES) to generate electrical power instead of using fossil fuel sources, which have negative effects on the environment [1].RES is environment-friendly and can replace all the conventional sources for power supply.

Figure 1:Assessment criteria for using RES technologies

There are several kinds of RES, which can be used in any project.Furthermore, these can be separate sources or hybrid systems which depend on the location of the case study (project)and availability of the natural sources, i.e., wind energy, solar radiation, falls or dams etc.The process of selecting the suitable kind of RES for any project is a strenuous task.For any successful RES project, decision-makers must study, analyze, and consider multiple factors such as economic, environmental, and technological factors.Decision-makers should resolve this problem in the framework to understand suitable source without concessions [2].The design of the RES uses several technical methodologies and algorithms based on the optimal analysis.Therefore, the process of selecting a powerful and optimum RES is a complex problem that needs evidencebased decision.Furthermore, the decision of using/installing a RES normally includes various stakeholders that may have various interests and goals related to the project.An assessment criterion for using RES technologies include technical, institutional, economic, environmental, and social/ethical criteria as is illustrated in Fig.1.

The variety of available technologies and equipment to achieve specific goals is a characteristic feature of the modern era.The ratios of parameters such as price, performance, reliability, durability, safety, environmental friendliness, ergonomics, etc.constitute the issue of optimal choice.In a situation when there are 5-7 choices with immense object attributes and competing criteria;there is no obvious choice for the optimal solution.

One of the approaches to solve this problem of optimal choice is the use of various multi-criteria solution techniques i.e., Multi-criteria decision-making (MCDM/MCDA) [3-5] and integration of methods into the engineered design process.Though, MCDM models are partially formalized and for them there is no concept of an absolute optimal solution.Nevertheless, as practice shows, MCDM models allow the selection of best option among predefined alternatives.If several alternatives have some of the attributes as “strong” and approximately the same part as “week,” then the performance indicators of such alternatives will differ significantly, and the alternatives will be hardly distinguishable.This requires a comprehensive analysis using various MCDM models and solution-based sensitivity analysis, where partially formalized quantitative (or qualitative) analysis is the basis for decision-making.Tab.1 reviews the previous literature on MCDM methods showing aim of the work, modeling, solution, and main results of the literature work [6-19].

Table 1:Literature review’s summary

Table 1:Continued

Table 1:Continued

Table 1:Continued

The main contributions of this research work can be summarized as:

(1) Based on multi-criteria decision-making models, implementing the integration of MCDM rank methods into the process of engineering design of hybrid renewable energy systems.

(2) Four different methods of weight estimation are considered for usage; no priority of criteria, criteria based on a pairwise comparisons matrix, CRITIC-method, and entropybased method.

(3) A step-by-step methodology for forming various MCDM models and subsequent analysis of the results is described.

2 Description of the Case Study

The case study represents a flat 70 acres site in the Al-Minya city (Egypt) as an example of the far region location.The latitude and longitude is 28°N and 30°E, respectively.The nearest electrical grid point was 12 km from this location.The site is wealthy with solar irradiance.The mean solar irradiance level was 5.97 kWh/m2/day [20].The maximum and minimum solar radiation values were 8.056 and 3.555 kWh/m2/day, respectively, for June and December.The sunshine period was about 9 to 11 h per day all year except for a few cloudy days.The hourly solar irradiance profile for every month is presented in Fig.2.

On the site, there was a well with the following specifications:150 m depth, 40 m static level of water, 120 m3the hourly rate of discharge.The salinity of brackish water was 2500 mg/l.It had been scheduled to cultivate part of the land with crops using the raw brackish water.The reminder was cultivated with Wheat as it cannot grow with brackish water.The salinity of the water needed to be lower than 800 mg/l.The amount of desalinated water was 250 m3/day.The required amounts of brackish water were 350-500 m3/day and 250-300 m3/day in the summer and winter periods, respectively.

The required energy to extract the brackish water was around 110 kWh/day with a peak of 15 kW.For desalination of the brackish water, it was scheduled to employ a reverse osmosis(RO) unit.Two different sizes of RO units, RO-250 and RO-500, were considered.The electrical peak demand values were 15 and 29.5 kW for RO-250 and RO-500 respectively [21].To collect 250 m3, RO-250 operates 24 h every day.So, the total required energy was 360 kWh/day.While,to collect the same amount by RO-500, 12 h of operation was required with a total consumption of 354 kWh/day.

Figure 2:The solar irradiance profile for every month (kW/m2) per day, Al-Minya city (Egypt)

Figure 3:Renewable energy system (RES) schematic graph, Al-Minya city (Egypt)

Fig.3 illustrates the suggested hybrid system, which contains a Fixed PV array at a tilt angle of 28-degree, DG, power conditioning unit, and battery storage bank.The model of battery was Trojan L16P (360 Ah, 2.16 kWh).The input techno-economic specification data for different elements of a hybrid system are shown in Tab.2 [22-24].These data were used for determining the system’s best sizes using HOMER?software [25,26].Various options containing three different power sources:only DG, PV-B system, and hybrid PV-DG-B, two different sizes of reverse osmosis(RO) units; RO-250 and RO-500, two strategies of energy management; load following (LF) and cycle charging (CC), and two sizes of DG; 5 and 10 kW were taken into account.The optimal size of each option was selected based on the minim cost of energy.

To identify the best option for the case study, MCDM tools were used.Eight parameters, including operating cost, renewable fraction, initial cost, the cost of energy, excess energy,unmet load, breakeven grid extension distance, and the amount of CO2were considered during the determination of the best option.Tab.3 shows the output of eight parameters for every considered option.

Table 2:Specification of different elements of the hybrid system

Table 3:The eight output parameters for all alternatives

3 Basic Methods and Formulas of MCDM

3.1 Aggregation Models

3.1.1 SAW(Simple Additive Weighting)Method[3-5]

Performance indicatorQiof thei-th alternative was determined as the entire standardized estimations of the attributesrijwith the weightwjof thej-th criteria:

3.1.2 COPRAS(Complex Proportional Assessment)Method[27]

The aggregation method uses the construction of a performance indicator of alternatives based on the homogeneous function of the two argumentsS+iandS?i:

where;

The above equation represents the sum of the normalized attribute values with weight revenue criteria and cost criteria.The best alternative was the one with the most elevatedQiscore.

3.1.3 TOPSIS(Technique for Order of Preference by Similarity to Ideal Solution)[3]

To determine the performance indicator of thei-th alternativeQi, a homogeneous function was used:

where;

S+iandS?iwere the distancesdbetween the ideal and anti-ideal objects respectively.Whereas,the alternativeAiin then-dimension attributes space, which are defined in one of theLp-metrics(Section 3.2).The TOPSIS ranking result depends on the choice of distance metric.The best alternative was the one with the highestQiscore.

3.1.4 GRA(Gray Relation Analysis)[28]

It evaluates the effectiveness of alternatives in two groups with respect to ideal and anti-ideal objects.The sequence of calculations is as follows:

Step 1:Define two sets of attributes i.e., ideal and anti-ideal:

Step 2:Determine the matrix of deviations of normalized values from the ideal and anti-ideal:

Step 3:Determine the matrices the gray relational coefficient:

Step 4:Determination of the indicator performance for the alternativeQi:

Here, the best alternative was the one with the highestQiscore.

3.1.5 VIKOR(VIsekriterijumsko KOmpromisno Rangiranje)[29]

Step 1:Determination of “ideal” and “anti-ideal” object can be expressed as:

Step 2:Weighted normalization:

Step 3:MaximalRand the group utilitySstrategies can be expressed as:

Step 4:Calculate the values ofQi:

Here, vassumes the part of balancing factor between the general advantage (S) and the maximum individual deviation (R).Smaller estimations ofvaccentuate bunch gain, while bigger qualities increased the weight controlled by singular deviations.“Voting by majority rule” (v>0.5);or “by consensus” (forv=0.5); or “with a veto” (for v<0.5).

Step 5:The aftereffect of the system is the three-rating recordsS,R, andQ.The options were assessed by arranging the estimations ofS,R, andQmodels of the base worth.

Step 6:As a compromise arrangement, optionA1was proposed, which was best assessed byQ(minimum) if the accompanying two conditions were met:

Condition C1:“Allowable advantage”:

where;A2is an alternative in contrast to the second situation in theQranking rundown:

Condition C2:“Adequate soundness in decision-making”:AlternativeA1ought to likewise be best assessed bySor/andR.

Step 7:If one of the conditions 1 or 2 was not fulfilled, a lot of negotiating arrangements were proposed, which comprises of:

AlternativesA1andA2; if conditionC2is not met, or

AlternativesA1,A2,...,Ak; if conditionC1is not fulfilled.Where,Akis controlled by the connection:

3.1.6 PROMETHEE(Preference Ranking Organization Method for Enrichment Evaluations)[30]

Step 1:Set the inclination work for two items for every modelHj=H(dis,p,q).When in doubt, they have two boundaries:p-indifference edge, it mirrors the way that if the distinction of these estimations of two optionsiandsare immaterial, objects by standardjwere comparable.If the distinction in the limit esteempwas surpassed, an inclination connection was built up between the items.Similarly, if the distinction in edgeqwas surpassed, the inclination work, which compares the “strong preference” of variationiconcerningsvariation as for thejmeasure.With the distinction ofdisin the stretch fromptoq, the inclination work was under 1, which compares a “weak preference”.The decision of the inclination work was controlled by the leaders.A few sorts of capacities favoredH(d) were introduced in Tab.4.

Table 4:Preference functions for PROMETHEE-method

Step 2:Compute the distinction in the estimations of the models for the two items and calculate the inclination recordsV:

Step 3:Determination of the preference factors:

The best option is the one with the most elevatedQiscore.

3.1.7 ORESTE(Organization,Arrangement to Sinteze of Relational Data)[31]

Step 1:Change from network DM to ranks matrix (the columns of the matrix are supplanted by their ranks).

Step 2:Determine the ranks of criteria:

Step 3:The projections of ranks were computed:

p-one of:p=1 Average (CityBlock, TaxiCab or Manhattan) distance,

p=2 Mean Square (Euclidean) distance,

p=∞Chebyshev distance.

Step 4:Calculating Ranksdij

Step 5:Calculating RanksRi(ORESTE 1)

Step 6:Calculation of preference factorsCik

The best option was the one with the lowestQiscore.

3.2 Distance Metric

Distance metric was used to choose measurement to quantify the distance between twondimensional objectsxandy:

3.3 Normalization Methods

The study used the following 8-normalization methods [32,33] as listed in Tab.5.

Table 5:The 8-normalization methods used in the study

Table 5:Continued

3.4 Normalization for Cost Criteria

Two-Step Res-algorithm for inversion of cost credits into advantage ascribes [33]:

where; the linear normalization methodNorm(aij) in the first step was applied to both, the benefit attributes and the cost attributes, and indexj* meets the cost criteria.

3.5 Methods for Weight Estimation

3.5.1 No Priority of Criteria

3.5.2 Estimation the Weights Based on a Pairwise Comparisons Matrix of the Criteria[34]

Step 1:Determine pairwise comparison matrixPinT,Saaty scale.

Step 2:Determine eigenvector(v)for max eigenvalueλmaxof matrixP(the calculations use the MATLAB functioneigs()):

Step 3:Calculate consistency index (C.I.) and compare with Tab.5 of random consistency index (R.I.(n)).

Step 4:Check the consistency of the pairwise comparison matrix.Compare the consistency index (C.I.) with the values of (R.I.(n)):

if C.I./R.I.(n)<0, 1‘Pairwise comparison matrix P consistent’else go to Step 1 end

Step 5:Calculate weights of criteria:

In this case study, the values of thePwere listed in Tab.6.

Table 6:Values of the P-matrix in the case of study

C.I./R.I.(n)=0,07.

3.5.3 CRITIC-Method (Criteria Importance Through Intercriteria Correlation) for Weight Estimation[12]

Step 1:Determine ‘best’(b) and ‘worst’(t) solution ([1×n]-vector) for all attributes.

Step 2:Determine relative deviation matrixV=(vij)[m×n]-matrix

Step 3:Determine standard deviation (St) ([1×n]?vector) forcollsofV.Step 4:Determine correlation matrix (Cr) ([n×n]) forcollsofV.

Step 5:Determine vector (c) and calculate the weight of criteriawk.

In this case study, the values of theB,T,and wwere listed as:

3.5.4 Entropy-Based Method for Weight Estimation[13]

The step-by-step algorithm for estimating the weights of the criteria using the entropy method was presented as follows:

Step 1:Standardized decision matrix (Max-Min method) for benefit criteria:

For cost criteria:

Step 2:Calculate the equity contribution of thei-th attribute for each criterion:

Step 3:Calculate the entropy of each criterion:

Step 4:Calculate the weight of each criterion:

The value (1?ej) is the internal intensity of the contrast of each criterion or is the degree of divergence of the internal information of each criterion [35,36].The smaller value of the entropy,the larger the entropy-based weight.In this case study, the values of thewwere listed as:

w=[0,093 0,082 0,106 0,023 0,180 0,032 0,086 0,398].

4 MCDM Techniques

4.1 MCDM Rank Model

The ranking-based MCDM model for every electiveAidecides a specific exhibition level of the choicesQibased on the ranking of the other options, and the ensuing decision-making was obtained [3-5]:

The MCDM rank model incorporates the decision of a lot of alternatives (A) and a set of criteria (C), an evaluation of the estimations of the characteristics of choices with respect to every criterion-a decision-making matrix (aij), a method for assessing the weight or priority of criteria(w), a choice of a normalization method (‘nm’) decision-making matrix, a choice of metric for calculating distances inn-dimensional space of criteria (‘dm’), a choice of preference functions(‘pr’), and the definition of aggregation function of alternatives’attributes (F) to calculate efficiency indicator (Q) of each alternative.Based on the calculation of the aggregate performance indicator of alternativesQ, alternatives were ranked, i.e., SAW ranking model is simplified as:

whererijare the standardized estimations of the regular estimations of the qualitiesaij, acquired utilizing one of the standardization techniques.None of the arguments to F were unambiguous.The choice of A and C was not formalized, the estimatesaijwere not accurate, the choice of the method for evaluating the weights of the criteria, the method of the standardized method, the method of aggregation, and the choice of the distance metric were not formalized, as there were no selection criteria.Therefore, different combinations of 8 basic parameters in Eqs.(44)-(46)define different MCDM models.

4.2 MCDM Formalization of the Problem of Choosing Hybrid Renewable Energy Systems

In the present study, different models were defined by combining the‘nm’normalization method, the aggregation methodF, the choice of different distance metrics, and different preference functions.Thus, for integration into the engineering design process of hybrid renewable energy systems, 55 models or variations of the basic ranked MCDM methods have been identified in Tab.7.

Table 7:Constructor of alternative ranking models

Combining 4 different methods for evaluating criteria weights, gives 220 options for ranking alternatives.Besides, the following model notation was used in the form of a cortege:#={‘F’,‘w’,‘nm’,‘dm’,‘pr’}.For example, model #18 = {TOPSIS, (ii), Sum,L1} uses the TOPSIS attribute aggregation method, (ii)-a method for evaluating criteria weights, the Sum normalization method and the CityBlock-metric [3-5,27-31].

How much the ranking results differ depends on many factors.First, a ranking of alternatives was determined by the partial preference of various alternatives among themselves according to individual attributes.Suppose one of the alternatives has a preference over the other alternative according to several criteria, and vice versa, the other alternative dominates over the first according to another group of criteria.In that case, the performance indicators of these alternatives differ slightly.Although the aggregation methods, normalization methods, and the choice of parameters for preference functions affect the ranking result insignificantly, their small variations,together with a weak distinction of alternatives determine the ranking results [37].Another parameter determining the ranking is the criterion weight.The weights directly determine the preference of alternatives over each other according to certain criteria.Therefore, the assessment of the criteria weights requires a justification of the chosen method and subsequent comparative analysis and correction of the weights of various criteria.Following this, one of the tasks of the study is to determine several best alternatives based on the analysis of the ranking results when varying the methods and parameters in the MCDM models.

4.3 Description of Alternatives and Attributes of Hybrid Renewable Energy Systems;Decision Matrix

Tab.8 presents a matrix of decisions for the selected list of alternatives and their attributes.

Table 8:Matrix of decision D [12×8]

All criteria except the secondwere “cost” criteria (C?j).Therefore, to aggregate the attributes of benefit and cost jointly, the inversion of the normalized values was used through the ReS-algorithm [33].The ReS-algorithm allows the same normalization method to be applied to both benefit and cost attributes and is effective for all normalization methods.

To determine the priority of alternatives, it was not enough to compare the absolute values of the efficiency indicatorQi.Attribute values may not be accurate due to many factors.For example, an attribute can be measured where the data source may be unreliable, there was error in measurement, the measurements for various alternatives were carried out using different methods,some attributes may be random values or determined by the values of intervals, etc.Thus,the value of the performance indicator was determined with an error ofQi±ΔQi, and the distinguishability of alternatives was determined by the errorΔQi.

In many cases, it was not possible to estimate the error.Then use the “a priori” or expert estimate, expressing it as a percentage.For example, as follows:the error in assessing the indicator of the alternative’s effectiveness was 5% of its value.Considering that alternatives were ranked according to their place in the ordered list of performance indicators, it was advisable to determine a relative indicator to assess the distinguishability of alternatives:

where;Qpis the value of the performance indicator corresponding to the p-rank alternative,rng(Q)=QI?Qm.Following Eq.(48), dQprepresents the relative (given in theQscale) increase or decrease in the efficiency indicator for an ordered list of alternatives.Afterward, two alternatives:the relative increase in dQof which differ less than the value of the given a priori error, should be considered indistinguishable.The dQindicator was used to assess the distinguishability of alternatives and to compare the results of aggregation performed by different methods.

4.4 Estimation of Weights of Criteria

Tab.9 presents the criterion weights obtained using the 4-methods of estimation [12,13,34,36].

Table 9:Values of the weighting coefficients of the criteria obtained by using various methods

In the second part of Tab.9, values of the relative difference (%) were given between the weights obtained by different methods.For the highlighted cells, the criteria weights differ significantly more than 70%.

For method (iv), the weight values forC8were greatly overestimated (by 4-5 times) and the weight ofC4andC6was greatly underestimated (by 4-5 times) in comparison with the weights determined by other methods.This overestimates the contribution of attribute 8 to the performance indicator of alternatives.It was expected that the priority 8 attribute alternatives will receive priority in the performance indicator.These are alternatives toA2,A8,A9, andA10.

The weights obtained using methods (ii) and (iii) differ on average by 30%.Both methods consider the relationship between different criteria in general, rather than the difference in attributes like the entropy-based method.

5 Results and Discussion

Calculations for various models were performed using the MCDM_tools software (version 2020), developed in the MATLAB system.MCDM_tools (version 2018) were posted in the public domain in a MathWorks File Exchange service on the website of the developer company Math-Works [38].For each MCDM model, the performance indicator of each alternativeQi, the relative intensity iQ, the relative increment dQwere calculated and the ranks of the alternativesAiwere determined.An example of calculated indicators was presented in Tab.9.

Tabs.10, 11 present the synthetic results of ranking alternatives (based on Tab.7) for various options (i)-(iv) estimates of the criterion weights obtained for 55 different MCDM models.

Table 10:Fragment of calculation results for the MCDM model #17-24 = {TOPSIS, (i), (Max,Sum, ..., Sgm), L1,}

Table 11:Statistics of alternatives of I-III ranks based on the results of calculations of 55 MCDM models

The numerical values in the Tab.11 indicate how many times each of the alternatives A was ranked as I, II, III when ranks were based on 55 variants of MCDM models.

First:Let us consider the clearly “weak” alternatives that, according to the results of calculations, did not have high ranks.According to Tab.8, it isA11andA12.

Second:The assumption (Section 3.4) that, for the method (iv), all alternatives with a priority on attribute #8 will also receive priority in the performance indicator.

Indeed, according to Tab.8 alternatives,A9andA8have I and II ranks in most models, and alternativeA2has II and III ranks in one and two cases, respectively, (and for other methods of estimating weights, I-III ranks are never achieved).The final ranks of the alternatives were presented in Tab.12.

Table 12:Final rank of alternatives

The unconditional leader was alternativeA9.However, the alternativesA10,A7,A4, andA3for some models (about 30% of variants) also had the first rank.Determining the leader by majority of votes cannot be a correct method.Additional information consists of assessing the distinguishability of alternatives using the relative performance indicators iQand dQ.

Tab.13 shows the ranks of alternatives (fragment) based on the results of calculations for 55 models in the case of determining the weights of the criteria by method (iii).

Table 13:Ranks of alternatives based on the results of 55 models in the case of determining the weights of criteria by method (iii) (Fragment)

Tab.13 and the data in Tab.7 make it possible to select models for subsequent refinement of the leader.The specificity of MCDM models shows that for some models (more precisely,an unsuccessful combination of model parameters), a result is possible in which an alternative“weak” in terms of characteristics has a high rating (rank).For example, the alternativeA1has shown II rank in the model #38={TOPSIS,(iii),St4,L∞}, Tab.11.In the absence of formal criteria for the selection of models, acceptance or rejection was possible if there were additional arguments (reasons).

Fig.4 shows various histograms of the ranks of alternatives based on the results of calculations in 55 models (option (ii)).Data were collected in separate histograms considering the distinguishability of rank I-III alternatives.

Figure 4:Distinguishability of the alternatives of rank I-III for (ii)-method

Fig.5 shows the intensity of the performance indicator of alternatives of I-III ranks (points),considering their distinguishability for various models of aggregation of attributes.The results were collected in sequential groups corresponding to the eight different normalization methods as indicated in Tab.2 (Model Builder)-{‘Max,’‘Sum,’‘Vec,’‘Max-Min,’‘Dea,’‘St4,’‘Spl,’‘Sgm.’}Colored markers illustrate the distinguishability of rank I-III alternatives.

Similar Figs.4 and 5 were obtained for all the variants for evaluating the weights of criteria(i)-(iv).Tab.14 presents a summary of the results.

The distinguishability of alternatives of I-III ranks was no more than 61.8%, and the indifference of alternatives I, II, and III ranks above 30% cannot be made unambiguously.AlternativesA9,A10,A8, andA4were recommended as suboptimal.The final decision was made by the decision-maker.

Table 14:Distinguishability of rank I-III alternatives; statistics for 55 models

Figure 5:Distinguishability of the alternatives of rank I-III for the ii-method of weight estimation and linear normalization method

6 Conclusion

Based on various multi-criteria decision-making (MCDM) models, the best electrical energy source has been identified to feed both the water pumping system and the water desalination unit, respectively.The electrical energy source alternatives were suggested to consider different sizes of water desalination units, different energy management strategies, different sizes of diesel generators, and different system configurations.Four different methods of the weight estimation were considered; no priority of criteria, based on a pairwise comparisons matrix of the criteria,CRITIC-method, and entropy-based method.The results revealed that the best/optimal alternative of hybrid PV-DG-B consists of 5 kW DG, RO-500, and load following control strategy.The yearly operating cost and initial cost for such a system were $ 5546 and $ 161022, respectively, while the cost of energy was 0.077 $/kWh.The excess energy and unmet loads were 40998 and 2371 kWh,respectively.The breakeven grid extension distance and the amount of CO2were 3.31 km and 5171 kg per year.Compared with DG only, the amount of CO2has been sharply reduced by 113939 kg per year.

Funding Statement:The authors received no specific funding for this study.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.