Multi-Objective Decision-Making for Hybrid Renewable Energy Systems for Cities: A Case Study of Xiongan New District in China

: The application of renewable energy has become increasingly widespread worldwide because of its advantages of resource abundance and environmental friendliness. However, the deployment of hybrid renewable energy systems (HRESs) varies greatly from city to city due to large di ﬀ erences in economic endurance, social acceptance and renewable energy endowment. Urban policymakers thus face great challenges in promoting local clean renewable energy utilization. To address these issues, this paper proposes a combined multi-objective optimization method, and the speciﬁc process of this method is described as follows. The Hybrid Optimization Model for electric energy was ﬁrst used to examine ﬁve di ﬀ erent scenarios of renewable energy systems. Then, the Technique for Order Preference by Similarity to an Ideal Solution was applied using eleven comprehensive indicators to determine the best option for the target area using three di ﬀ erent weights. To verify the feasibility of this method, Xiongan New District (XND) was selected as an example to illustrate the process of selecting the optimal HRES. The empirical results of simulation tools and multi-objective decision-making show that the Photovoltaic-Diesel-Battery o ﬀ -grid energy system (option III) and PV-Diesel-Hydrogen-Battery o ﬀ -grid energy system (option V) are two highly feasible schemes for an HRES in XND. The cost of energy for these two options is 0.203 and 0.209 $ / kWh, respectively, and the carbon dioxide emissions are 14,473 t / yr and 345 t / yr, respectively. Our results provide a reference for policymakers in deploying an HRES in the XND area.


Research Background
China is the largest energy consumer worldwide due to its continued industrialization and urbanization, and the majority of the energy supply is still provided by non-renewable resources such as coal and natural gas. To cope with this situation and to promote ecological environmental improvement, the Chinese government has decided to intensely develop solar, wind, hydropower, and other renewable energy sources. Some related policies and regulations reflect the detailed staged goals of China and establish feasible ways to achieve these goals up to 2050, which provides policy support for the development of renewable energy. The Renewable Energy Law of the People's Republic of China [1] provides a new framework for the development of renewable energy from the aspects

The Potential Contribution of This Study
The priority of a candidate option in HOMER is determined by the principle of the minimum NPC. Thus, HOMER is a great tool for simulating scenarios in undeveloped areas where the economy and technology are the obvious limiting factors for local development. However, for developed areas, the limiting factors may be the society and environment, rather than the economy and technology. Although some predecessors have used the method of combining TOPSIS and HOMER to evaluate HOMER simulation results, these studies have some shortcomings, which can be divided into two categories. First, these studies did not consider the impacts of local policies on system selection [26]. For example, political willingness plays a special role in regional development in China. This method considers comprehensive indicators affecting the selection of HRES options, including economic (the initial cost, O&M cost and COE), technological (renewable fraction, ease of installation and excess electricity), environmental (CO 2 emissions, area requirements and environmental impact) and social (acceptability and willingness to pay) factors. Different from previous studies, social indicators related to local energy policies are included in the index evaluation system for the first time. In addition, the subjective weighting method represented by the AHP method has strong subjectivity, which ignores the variability and conflict within indicator elements. Second, the objective weighting method represented by the CRITIC method has strong subjectivity, which is not in line with the reality of real-world decision-making [28]. Thus, using only subjective weighting or objective weighting can introduce certain errors into the results.
Combining subjective and objective weighting to reduce errors, this paper presents an evaluation method based on the HOMER model and the TOPSIS method. This contribution unifies the indicator unit in the evaluation process of a complex system, which is helpful to the scientific ranking of alternatives. Based on this method, we further determine the optimal HRES option in the XND in China.

Methodology
This paper proposes an evaluation method combining the HOMER model and the TOPSIS method. Generally, a large amount of heat energy will be released in the process of coal-fired power production in China, which is used as a by-product of electric energy to provide heating services. Since this study focuses on the optimization of the power system, the research does not consider the heat energy supply. The energy system consists of the already-produced power to the end-users. In this text, the HOMER model is used to simulate different scenarios in the research area, and the simulation results are evaluated as options. The TOPSIS method is utilized to rank the HRES options, considering all aspects affecting HRES selection. In the decision-making process of the TOPSIS method, three different weighting methods are proposed-the AHP method, CRITIC method, and AHP-CRITIC method-which represent subjective empowerment, objective empowerment and the combination of subjective and objective empowerment. Figure 1 shows a flow chart of this method, and the details of each step are presented in the following sections.

HOMER Simulation
HOMER was developed by the NREL. The main functions of this model include system simulation, system optimization, and sensitivity analysis, which can be used for system design, energy planning, market analysis and policy making. To obtain candidate HRES options in the research area, HOMER based on the minimum net present value is applied.

HOMER Simulation
HOMER was developed by the NREL. The main functions of this model include system simulation, system optimization, and sensitivity analysis, which can be used for system design, energy planning, market analysis and policy making. To obtain candidate HRES options in the research area, HOMER based on the minimum net present value is applied.

Description of the XND
To focus on the non-capital functions of Beijing in China-considering the application of a new optimal development mode for densely populated areas; adjusting and optimizing the layout and spatial structure of the Beijing, Tianjin and Hebei districts; and cultivating a new engine for innovation-driven development-the government of China decided to develop the XND area.
Located in the hinterland of Beijing, Tianjin and Baoding at 38°51.8′ N, 115°52.2′ E, the XND area covers the administrative districts of Xiongxian, Rongcheng and Anxin counties, as well as the surrounding areas, which have distinct regional advantages, such as a low extent of current development, basic conditions that provide a good starting point and high development and construction standards.
The plan involves developing the XND area in advance, starting with specific areas, and the planned start-up area of the XND area is 38 km 2 , in which the planned construction land area is 26 km 2 , accounting for 68.4% of the total planning area according to the local regulations published in [29]. This area will inherit the non-capital function stock of Beijing and effectively attract part of the population of Beijing, thereby accruing highly skilled workers, and a population structure that is

Description of the XND
To focus on the non-capital functions of Beijing in China-considering the application of a new optimal development mode for densely populated areas; adjusting and optimizing the layout and spatial structure of the Beijing, Tianjin and Hebei districts; and cultivating a new engine for innovation-driven development-the government of China decided to develop the XND area.
Located in the hinterland of Beijing, Tianjin and Baoding at 38 • 51.8 N, 115 • 52.2 E, the XND area covers the administrative districts of Xiongxian, Rongcheng and Anxin counties, as well as the surrounding areas, which have distinct regional advantages, such as a low extent of current development, basic conditions that provide a good starting point and high development and construction standards.
The plan involves developing the XND area in advance, starting with specific areas, and the planned start-up area of the XND area is 38 km 2 , in which the planned construction land area is 26 km 2 , accounting for 68.4% of the total planning area according to the local regulations published in [29]. This area will inherit the non-capital function stock of Beijing and effectively attract part of the population of Beijing, thereby accruing highly skilled workers, and a population structure that is compatible with the function positioning of innovative intelligence and other new areas will gradually form to satisfy the development needs of high-end and prosperous industries. This paper selects this area as the research area, and Figure 2 shows the detailed information of the research area in this study. The sources of the pictures on the left and right of Figure 2 are HOMER and the Planning Outline of the XND area in Hebei Province [30], respectively.

Hybrid Energy System Options
According to the construction plan, the XND area will realize a 100% clean power supply and become the first urban power grid with a continuous 100% clean power supply. This plan will enable clean power consumption to account for 52% of the final energy consumption. To fulfil this purpose, it is necessary to make full use of local renewable resources to optimize the power system configuration. The policy planning outline (Hebei 2017) shows that clean energy sources include wind, solar, hydrogen, and other renewable energy with no waste. In view of the dense population, nuclear power is not feasible in the XND. This paper applied HOMER to simulate HRESs with different storage modes, including battery and hydrogen energy storage equipment, which are applied to solve the problems of excess electricity and power output volatility. In the present study, PV panels, wind turbines, batteries, fuel cells, electrolytic cells, hydrogen tanks, converters, and diesel generators constitute the candidate components of the HRES options. This study constructed five different scenarios in the XND area, including a grid-connected HRES with a battery storage mode (option I), a grid-connected HRES with a hydrogen storage mode (option II), an off-grid HRES with a battery storage mode (option III), an off-grid HRES with a hydrogen storage mode (option IV), and an off-grid HRES with modes that combine battery and hydrogen storage (V). Figure 3 shows the detailed configurations of the above systems.

Hybrid Energy System Options
According to the construction plan, the XND area will realize a 100% clean power supply and become the first urban power grid with a continuous 100% clean power supply. This plan will enable clean power consumption to account for 52% of the final energy consumption. To fulfil this purpose, it is necessary to make full use of local renewable resources to optimize the power system configuration. The policy planning outline (Hebei 2017) shows that clean energy sources include wind, solar, hydrogen, and other renewable energy with no waste. In view of the dense population, nuclear power is not feasible in the XND. This paper applied HOMER to simulate HRESs with different storage modes, including battery and hydrogen energy storage equipment, which are applied to solve the problems of excess electricity and power output volatility. In the present study, PV panels, wind turbines, batteries, fuel cells, electrolytic cells, hydrogen tanks, converters, and diesel generators constitute the candidate components of the HRES options. This study constructed five different scenarios in the XND area, including a grid-connected HRES with a battery storage mode (option I), a grid-connected HRES with a hydrogen storage mode (option II), an off-grid HRES with a battery storage mode (option III), an off-grid HRES with a hydrogen storage mode (option IV), and an off-grid HRES with modes that combine battery and hydrogen storage (V). Figure 3 shows the detailed configurations of the above systems.

Hybrid Energy System Options
According to the construction plan, the XND area will realize a 100% clean power supply and become the first urban power grid with a continuous 100% clean power supply. This plan will enable clean power consumption to account for 52% of the final energy consumption. To fulfil this purpose, it is necessary to make full use of local renewable resources to optimize the power system configuration. The policy planning outline (Hebei 2017) shows that clean energy sources include wind, solar, hydrogen, and other renewable energy with no waste. In view of the dense population, nuclear power is not feasible in the XND. This paper applied HOMER to simulate HRESs with different storage modes, including battery and hydrogen energy storage equipment, which are applied to solve the problems of excess electricity and power output volatility. In the present study, PV panels, wind turbines, batteries, fuel cells, electrolytic cells, hydrogen tanks, converters, and diesel generators constitute the candidate components of the HRES options. This study constructed five different scenarios in the XND area, including a grid-connected HRES with a battery storage mode (option I), a grid-connected HRES with a hydrogen storage mode (option II), an off-grid HRES with a battery storage mode (option III), an off-grid HRES with a hydrogen storage mode (option IV), and an off-grid HRES with modes that combine battery and hydrogen storage (V). Figure 3 shows the detailed configurations of the above systems.

AHP Method
The AHP method was proposed in the mid-1970s by the American operations research scientist Tomas L. Saaty [31]. This method regards complex multi-objective decision-making problems as systems, decomposes the objective into multiple objectives, and then obtains several levels of multi-indexes through further decomposition. Through the nine-scale scoring method put forward by Saaty, single and total rankings are calculated, which can be used as a systematic method of multi-index and multi-scheme decision-making [32]. The calculation method for each evaluation indicator and the steps of the AHP method are explained in Sections 2.2.1 and 2.2.2, respectively.

The Calculation Process of Each Indicator
To comprehensively evaluate the HRES options in the XND area, eleven indicators covering technology, economy, environment and society are selected, which can be divided into qualitative indicators and quantitative indicators. The detailed calculation processes of the quantitative indicators and evaluation criteria of the qualitative indicators are provided in this section.
(1) Renewable fraction The renewable fraction is a quantitative and technical index, and this indicator refers to the fraction of energy provided by renewable energy resources to the load, which is calculated by the following equation: where E nonren is the electricity generated by a non-renewable energy resource (kWh/yr), H nonren is the thermal energy generated by a non-renewable energy resource (kWh/yr), E served is the total electrical load supplied (kWh/yr), and H served is the total thermal load supplied (kWh/yr). In HRES deployment, the larger the fraction of renewable energy, the smaller the impact on the environment; therefore, the renewable fraction is a benefit indicator.
(2) Ease of installation The ease of installation is a qualitative and technical indicator. The number of system components and the ease with which the various components fit together are used to evaluate this indicator. In the deployment of renewable energy systems, an easier installation will reduce the investment cost. Thus, the ease of installation is a benefit indicator.
(3) Excess electricity The excess electricity is a quantitative and technical indicator. In the HRES, the electricity generated from renewable energy sources first satisfies the local load demand, and the rest of the electricity is divided into two parts, with one part stored in batteries or hydrogen tanks and the other part wasted. The excess electricity is applied to describe the fraction of wasted electricity, and the smaller this indicator in the HRES is, the better the HRES functions. Thus, the excess electricity is a cost indicator, and this value is calculated by the following equation: where E excess is the total excess electricity (kWh/yr) and E prod is the total electrical production (kWh/yr).
(4) Initial cost The initial cost is a quantitative and economic indicator. This indicator refers to the investment required for each component at the beginning of the HRES. The smaller this indicator is, the lower the initial capital investment is. Thus, this parameter is a cost indicator. capital investment, and the total required investment is represented by the O&M cost. This parameter is a cost indicator.
(6) Cost of energy (COE) The COE is a quantitative and economic indicator. To specify the cost of generating electricity per kilowatt hour, the concept of the COE is introduced. The value of this indicator is calculated by the following equation: where C ann,tot refers to the total annualized cost ($/yr), c boiler refers to the boiler marginal cost ($/kWh), H served refers to the total thermal load supplied (kWh/yr) and E serverd refers to the total electrical load supplied (kWh/yr). In the HRES deployment, the lower the COE is, the lower the investment is. Thus, the COE is a cost indicator. (7) CO 2 emission CO 2 emission is a quantitative and environment indicator. In the configuration of an HRES, a diesel generator is used to mitigate the harm caused by fluctuations in renewable energy generation. In addition, diesel combustion produces a certain amount of carbon dioxide, which is represented by the CO 2 emission. Although a certain amount of carbon dioxide is produced during the life cycle of renewable energy components, the amount generated is much lower than that generated by diesel generators. Thus, this paper assumes that there are no carbon dioxide emissions in the process of electricity generation using renewable energy components. The value of this indicator is calculated by the following equation: where TCE refers to the total carbon dioxide emissions (t/yr), E i refers to the carbon emission intensity of diesel combustion (t CO 2 /MW h) and Q i refers to the electricity generated from diesel (MWh). The lower the CO 2 emission is, the smaller the greenhouse effect is. Thus, this parameter is a cost indicator.
(8) Area requirement This indicator represents the total floor area of each device, including PV panels, wind turbines, electrolysers, hydrogen tanks, and generators. The land area in the city is very limited, and the smallest HRES footprint is sought. Thus, the area requirement is a cost indicator. The installed surface area of a diesel generator station is relatively small, and this study takes only the installed surface area of renewable energy components into consideration. In this study, the installed surface areas of a 1-kW PV panel and a 1.5-MW wind turbine are 7.5 m 2 and 1500 m 2 , respectively. Table 1 shows the detailed information of land surface requirement for each option.
(9) Environment impact The environment impact is a qualitative and environmental indicator. This indicator refers to the noise impact caused by an HRES. In addition to the greenhouse effect, the noise impact also causes physical discomfort to humans, which is why we select this indicator to evaluate the candidate HRES schemes in this study. The environment impact is a cost indicator.
(10) Acceptability Acceptability is a qualitative and social indicator. This indicator refers to the degree of recognition and mastery of technology from the perspective of the public. In the XND, the public is more familiar with diesel generator technology than with renewable energy components technology and is more familiar with battery storage systems than with hydrogen storage systems. Acceptability is a benefit indicator.
(11) Willingness to pay Willingness to pay is a qualitative and social indicator, and this indicator represents the degree of willingness of investors or policymakers to invest in different HRESs. Hence, willingness to pay is a benefit indicator. This section describes the detailed calculation process of the AHP method in the XND area, and the specific steps are as follows: Step 1: Establish a hierarchical model Figure 4 shows the hierarchical structure of this study. The highest level in this study is the selection of the best HRES option in the target region. The renewable fraction, ease of installation and other nine indicators are the evaluation indicators in this study.
The environment impact is a qualitative and environmental indicator. This indicator refers to the noise impact caused by an HRES. In addition to the greenhouse effect, the noise impact also causes physical discomfort to humans, which is why we select this indicator to evaluate the candidate HRES schemes in this study. The environment impact is a cost indicator.
(10) Acceptability Acceptability is a qualitative and social indicator. This indicator refers to the degree of recognition and mastery of technology from the perspective of the public. In the XND, the public is more familiar with diesel generator technology than with renewable energy components technology and is more familiar with battery storage systems than with hydrogen storage systems. Acceptability is a benefit indicator. (11) Willingness to pay Willingness to pay is a qualitative and social indicator, and this indicator represents the degree of willingness of investors or policymakers to invest in different HRESs. Hence, willingness to pay is a benefit indicator.

The Process of the AHP Method
This section describes the detailed calculation process of the AHP method in the XND area, and the specific steps are as follows: Step 1: Establish a hierarchical model Figure 4 shows the hierarchical structure of this study. The highest level in this study is the selection of the best HRES option in the target region. The renewable fraction, ease of installation and other nine indicators are the evaluation indicators in this study. Step 2: Compare the different indicators When comparing the importance of the ith and the jth elements in the previous layer, use the quantitative relative weight aij for the analysis. The aij value in the pairwise comparison matrix is assigned according to the scale proposed by Saaty, and detailed information on the Saaty scale table can be found in Appendix A [33,34].
The pairwise comparison matrix has the following properties: aij ＞ 0, and aij = 1/aji; (when i is equal to j, the value of aij is 1). Step 2: Compare the different indicators When comparing the importance of the ith and the jth elements in the previous layer, use the quantitative relative weight a ij for the analysis. The a ij value in the pairwise comparison matrix is assigned according to the scale proposed by Saaty, and detailed information on the Saaty scale table can be found in Appendix A [33,34].
The pairwise comparison matrix has the following properties: a ij > 0, and a ij = 1/a ji ; (when i is equal to j, the value of a ij is 1).
Step 3: Calculate the weight vector and evaluate the consistency As shown in Tables 2-6, five n×n judgement matrices can be obtained in Step 2. Then, the maximum eigenvalue (λ max ) and corresponding eigenvectors are calculated for each paired comparison matrix, and the consistency test is performed. Detailed information on the consistency test can be found in Appendix A.     Step 4: Ranking each indicator The vector of each indicator can be obtained by multiplying the weight vector of the sub-criterion layer by the weight vector of the criterion layer. If the value of the indicator vector is high, this indicator has a higher importance in the whole HRES and vice versa.

CRITIC Method
The CRITIC method is an objective weighting method proposed by Diakoulaki in 1995 [28]. In this method, the objective values of each indicator element are determined based on the variability and conflict within the indicator element. The measurement of the variability and conflict within the indicator element is realized in the form of the standard deviation in the indicator element and the correlation coefficient between the indicator elements. The process of determining the objective weight of indicator elements by the CRITIC method can be divided into five parts, namely, the dimensionless treatment of the elements of the decision matrix, the calculation of the standard deviation, the calculation of the correlation coefficient, the calculation of the information value and the calculation of the weight of each indicator element. The detailed calculation process can be found in Appendix B. Step 1: The dimensionless treatment of the elements of the decision matrix The units of each indicator element in the initial decision matrix are different, and the differences between the indicator elements bring errors to the subsequent results. Before objective weight assignments using the CRITIC method, the first step is the dimensionless processing of the decision matrix.
Step 2: The calculation of the standard deviation of the indicator elements To quantitatively measure the variability within the evaluation indicator element, the standard deviation is selected as the tool with which to evaluate variability in this paper.
Step 3: The calculation of the correlation coefficient between indicator elements In order to quantitatively measure the conflict between indicator elements, this paper selects the correlation coefficient as the tool to evaluate the conflict.
Step 4: The calculation of the information value of each indicator element To comprehensively evaluate the variability and conflict within in each indicator element, this paper selects the information value of the indicator element to evaluate the variability and conflict within the indicator element.
Step 5: The calculation of the weight of each indicator element In the CRITIC method, the weight value of each indicator element is calculated according to the ratio of the information value of each indicator element to the sum of the information value of each indicator.

AHP-CRITIC Method
The AHP method determines the weight of each indicator element based on expert consultation and scoring and does not consider the objective information of the indicator elements themselves. In the CRITIC method, the objective information of the indicator element itself is the basis for determining the weight of the indicator element, without considering the subjective intention in the real world. To make full use of the advantages of the two methods, this paper proposes an AHP-CRITIC combination weighting method. The process of weighting can be achieved by the following equation.

TOPSIS Method
TOPSIS is an acronym for the Technique for Order Preference by Similarity to an Ideal Solution. This method was developed by C.L. Hwang and K. Yoon in 1981 [35]. The main idea of this method is to rank all candidate options according to their distance from idealized goals. If the evaluation option is the closest to the positive ideal solution and the farthest from the negative ideal solution, it is the optimal option, and if not, the evaluation option is not the best option. In this paper, TOPSIS is used to select the optimal HRES option in the XND area, and the calculation process can be found in Appendix B. TOPSIS is applied in this study according to the following process: Step 1: Formulate the decision matrix According to the AHP method, eleven indicators are used to evaluate the HRES options in the target area. These indicators can be divided into two types, namely, cost and benefit indicators. The quantitative indicator values in this study are calculated based on the HOMER simulation results, and the qualitative indicator values are determined according to expert opinions by the AHP method. Table 7 shows the detailed information of the decision matrix. The analysis process of each indicator can be found in Section 4. Step 2: Normalize the decision matrix To eliminate errors due to the order of magnitude, this paper normalizes the decision matrix. Detailed information on the normalization can be found in Appendix B.1, and Table 8 shows the detailed information of the decision matrix after normalization. Step 3: Calculate the weighted normalized matrix The weighted normalized matrix is obtained after the normalization process by Equation (A11) in Appendix C. There are three weighting systems in this paper, which are the weights determined by the AHP method, the CRITIC method, and the AHP-CRITIC method, and the weight determination process is shown in Sections 2.2-2.4.
Step 4: Determine the ideal solutions among all the candidate options In this step, the ideal solution has the shortest distance from the positive solution indicator and the longest distance from the negative solution indicator and is considered the optimal reference option.
Step 5: Calculate the relative distances of the option to the best and worst solutions In this step, the relative distances are calculated from the ideal and non-ideal solutions. The longer the distance of the option to the ideal solution, the more feasible the option is.
Step 6: Calculate the closeness coefficient The closeness coefficient is the relative closeness of each option to the ideal solution, which is used to determine the optimal option out of all the candidate options.

Parameter Setting
Detailed information on the data input is provided in this section. The data input is divided into two types: The HOMER model and the AHP method. In the section on HOMER, the power load estimation approach, available renewable energy resources and cost parameters are introduced. In addition, the establishment of the matrix weights and the consistency evaluation process of the AHP method are described in Section 3.2.

Electricity Consumption Estimation
To predict the population density of the XND area, some researchers use the grey system prediction model. According to the simulation results, the population density of the XND area is 685 people per km 2 [36]. The construction floor area is 26 km 2 in the start-up area of the XND, and the population of the research area in the future scenario is 17,810 people.
According to the technical provisions for urban planning management in Hebei Province, the average annual electricity consumption ranges from 6000 to 8000 kWh; that is, the estimated yearly electricity demand estimation ranges from 106.86 to 142.48 GWh. To ensure that the simulation results can meet the electricity demand of the region to the greatest extent, the annual electricity consumption considered for this region is 142.48 GWh, which means that the daily electricity consumption is 390 MWh. In addition, the energy system is subject to random variability, and the daily and hourly random variability in this study are 10% and 20%, respectively, according to previous research [37].

Available Renewable Energy Resources in the Research Area
To verify the optimal HRES option in the XND, the potential renewable energy resources in the research area must be evaluated. This study utilized meteorological data consisting of wind speed and solar energy data for the evaluation. Figure 5 shows detailed information on the local renewable resources. As shown in the left panel of Figure 5, the average annual solar global horizontal irradiance and clearness index are 4.07 kWh/m 2 /day and 0.508, respectively. The highest daily solar irradiation occurs in May and July, which means that there is abundant solar energy in these months. Moreover, the wind speed information is displayed in the right panel of Figure 5. According to this figure, the average annual wind speed is 3.52 m/s, which is measured by an anemometer at 50 m above sea level, and the highest wind speed is observed in April. The above information is retrieved from NASA meteorology and solar energy data resources.
Energies 2020, 13, x FOR PEER REVIEW 13 of 26 estimation approach, available renewable energy resources and cost parameters are introduced. In addition, the establishment of the matrix weights and the consistency evaluation process of the AHP method are described in Section 3.2.

Electricity Consumption Estimation
To predict the population density of the XND area, some researchers use the grey system prediction model. According to the simulation results, the population density of the XND area is 685 people per km 2 [36]. The construction floor area is 26 km 2 in the start-up area of the XND, and the population of the research area in the future scenario is 17,810 people.
According to the technical provisions for urban planning management in Hebei Province, the average annual electricity consumption ranges from 6000 to 8000 kWh; that is, the estimated yearly electricity demand estimation ranges from 106.86 to 142.48 GWh. To ensure that the simulation results can meet the electricity demand of the region to the greatest extent, the annual electricity consumption considered for this region is 142.48 GWh, which means that the daily electricity consumption is 390 MWh. In addition, the energy system is subject to random variability, and the daily and hourly random variability in this study are 10% and 20%, respectively, according to previous research [37].

Available Renewable Energy Resources in the Research Area
To verify the optimal HRES option in the XND, the potential renewable energy resources in the research area must be evaluated. This study utilized meteorological data consisting of wind speed and solar energy data for the evaluation. Figure 5 shows detailed information on the local renewable resources. As shown in the left panel of Figure 5, the average annual solar global horizontal irradiance and clearness index are 4.07 kWh/m 2 /day and 0.508, respectively. The highest daily solar irradiation occurs in May and July, which means that there is abundant solar energy in these months. Moreover, the wind speed information is displayed in the right panel of Figure 5. According to this figure, the average annual wind speed is 3.52 m/s, which is measured by an anemometer at 50 m above sea level, and the highest wind speed is observed in April. The above information is retrieved from NASA meteorology and solar energy data resources.

Cost parameters of the HRES options in the research area
In this paper, a techno-economic analysis was conducted with HOMER to determine the optimal HRES option in the XND area among five different scenarios. Table 9 lists the detailed cost information of each component [37,38].

Cost parameters of the HRES options in the research area
In this paper, a techno-economic analysis was conducted with HOMER to determine the optimal HRES option in the XND area among five different scenarios. Table 9 lists the detailed cost information of each component [37,38].

AHP Method
In this study, each candidate HRES option was evaluated from the perspectives of the economy, technology, environment, and society. The AHP method is applied in this study to determine the most influential indicator in HRES deployment. Table 1 summarizes the different weights of the indicators of the criterion layer, and most of these types are influenced by social indicators. Furthermore, the weight of each indicator in the corresponding type (sub-criterion or indicator criterion) is presented in Tables 5-9. From the evaluation results, the renewable energy fraction (with a weight of 0.63), COE (with a weight of 0.63), CO 2 emission (with a weight of 0.63), and acceptability (with a weight of 0.75) are the most influential indicators among the different types. The overall weight of each indicator is provided in the following part. The weight of each indicator is determined by the sum of the elements in a single row relative to the sum of the elements in each row.

Results and Discussion
The purpose of this study is to determine the optimal HRES configuration in the XND area by considering all factors, including technology, economy, environment and society. The HOMER simulation and the TOPSIS evaluation results are presented in this section.

HOMER Simulation Results
The HOMER model is applied to determine the candidate HRES options in the research area based on the simulation results. In this section, the different technical analyses of hybrid energy are divided into four parts: System configuration, economic analysis, technology analysis, and environmental analysis.

System Configuration
To determine the optimal configuration of each option, five different scenarios are built via HOMER simulations. According to the simulation results, each optimal HRES option is distinguished. Table 10 contains detailed information on each configuration in the research area and the various configurations identified as candidate options to be assessed by the TOPSIS method. It is worth noting that diesel generators are used in each optimal system because diesel generators can prevent system damage caused by the intermittent nature of renewable energy power generation.

Technology Analysis
This study selects the renewable fraction, excess electricity and grid power fraction as the evaluation criteria; thus, our technology analysis considers these three indicators. Detailed information on these two indicators can be found in Table 11. For the grid-connected HRES options, an amount of electricity is supplied by the grid, leading to a low renewable fraction. In this study, the grid power fractions of options I and II are 48.5% and 16.3%, respectively. In addition, for energy storage purposes, the effect of battery energy storage is better than that of hydrogen energy storage, which can be evaluated according to the excess electricity. In this study, the excess electricity of the battery storage mode option is lower than that of the hydrogen storage mode option, regardless of whether the option is a grid-connected or off-grid HRES option. Additionally, this paper constructs three different hydrogen energy storage HRESs to be evaluated because of the policy requirements of the XND area.
According to the development plan of the research area, the XND area will have more hydrogen energy vehicles and refuelling stations in the near future. Therefore, this study regards the hydrogen energy storage HRES as a plan option to be evaluated.

Economic Analysis
For each candidate HRES option to be evaluated, economic data including the initial capital cost, O&M costs and COE can be obtained from the HOMER simulation results. Table 12 summarizes the economic data of all candidate HRES options. The grid power and sellback prices are 0.1 $/kWh and 0.2 $/kWh, respectively. The grid electricity portion is large in the grid-connected HRES option; thus, the electricity price of the off-grid HRES is higher than that of the grid-connected HRES. In this study, the COEs of options I and II are lower than those of the other options. Moreover, the cost of hydrogen storage equipment is much higher than that of battery storage equipment, so the COE and initial cost of option II are higher than those of option I in the grid-connected HRES, and the COE and initial cost of options IV and V are higher than those of option III in the off-grid HRES. Although this study does not take the cost of the disposal process of waste batteries into consideration, the cost is too high to be ignored, which will bring obstacles into the deployment of HRES to some extent [39]. Carbon dioxide, as an indicator of climate change, cannot be ignored in the process of HRES selection. According to the simulation results of the HOMER model, the optimal configuration of carbon dioxide emissions in each scenario is obtained, which are 51,322 (option I), 24,405 (option II), 14,473 (option III), 58,995 (option IV), and 345 (option V) t/yr. From these results, we can conclude that options I and V have the highest and lowest CO 2 emissions because of the different renewable fractions, respectively. In this study, 632 g carbon dioxide is emitted into the atmosphere by the grid per 1 kWh of electricity provided.

AHP Results
The AHP method was applied in this study to calculate the relative importance of each indicator of the candidate HRES options. This section is divided into two parts, namely, the qualitative indicator evaluation and the indicator weight calculation by the AHP method.

Evaluation of the Qualitative Indicators
The indicators in this study can be divided into quantitative and qualitative indicators. The quantitative indicators are required for the HOMER simulation results, and detailed information on these indicators can be found in Section 3.1. In this section, we evaluate the qualitative indicators by their Saaty scale values, and Table 13 contains detailed information on these indicators in each candidate HRES option. From Table 13, we can obtain the priority order for each indicator in all candidate HRES options. At the technical level, we choose the high slope of the learning curve as our preferred criterion. The hydrogen storage mode equipment includes an electrolyser, hydrogen tank, dispenser, and fuel cell, which are more complex than battery storage equipment [40]. Therefore, options I and III are easier to install than the other options, and if grid power is considered, the installation will become complex, which is why the score of option III is higher than that of option IV. Furthermore, the environmental impact mainly refers to the noise impact on the environment, and wind turbines generate a large amount of noise in the process of generating electricity; therefore, options IV and V have the lowest scores for this cost indicator. In addition, the mastery of hydrogen energy storage technology in the research area is worse than that of battery energy storage technology, so the scores of acceptability of option I and III in this study are higher than those of other options. Lastly, since the disposal of waste batteries will have an impact on the environment, local policymakers increase their investment in hydrogen energy storage technology, so the willingness-to-pay scores of option I and III in this study are lower than those of other options.

The Overall Weight of Each Option
In this section, the results of the AHP method are summarized in Table 14. From this table, we can obtain the local and overall weights of each indicator. From the perspective of the local index weight, the renewable fraction, COE, CO 2 emission and acceptability are the most important indicators in each class, and the weights of these indicators in the classes are 0.63, 0.63, 0.63, and 0.75, respectively. Furthermore, CO 2 emission is the most important indicator among all indicators, and the overall weight of this indicator is 0.353, which means that policymakers should pay more attention to CO 2 emissions in the HRES deployment process, which is basically consistent with the local construction planning, and a 100% clean power supply will be realized in the XND, thereby allowing the area to become the first urban power grid to provide a continuous 100% clean power supply in China according to the local development planning.

CRITIC Results
According to the differences and conflicts between the indicator elements, the CRITIC method is used to assign weights to the indicator elements. In this study, this method is used to calculate the weights of 11 indicators in the indicator system.
As shown in Table 15, the standard deviation, information value and weight of each indicator are calculated. In terms of the standard deviation, the environmental impact, acceptability and willingness-to-pay indicators have the highest scores, which means that these three indicators have more variability within the indicator element. Furthermore, the willingness-to-pay indicator has the highest score, which means that this indicator is highly weighted, and the weight of this indicator is 0.19.

AHP-CRITIC Results
The weight determined by the AHP method is the product of expert scoring or consultation, which has certain subjectivity. Furthermore, this paper uses the method of combining subjective weighting and objective weighting to reduce the error in the decision result caused by subjectivity. Table 16 shows the detailed information of each indicator's weight.

TOPSIS Results
The TOPSIS method is applied to rank all the candidate HRES options in the research area. In the process of the TOPSIS method, once the decision matrix and weight sets are established, the next step is to evaluate all the candidate options using the TOPSIS method. To express the results more accurately, Benedetto and Kleme used the method of drawing a performance strategy map to directly express the advantages and disadvantages of the results [41]. In this paper, Figure 6 shows the score results using the TOPSIS method with three different weighting systems. In all three weighting methods, the least feasible scheme is option IV, and the reason behind this result is that the excess electricity in option IV is the highest among all the candidate options, which means that the majority of the electricity produced by renewable energy components is wasted. For option I and option II, part of the electricity is supplied by the grid, and this study considers that the electricity of the grid is provided by thermal power generation, which means that certain amounts of carbon dioxide emissions and higher grid transmission costs exist, so the feasibility of these two options decreases. In all three weighting methods, the most preferred option is option III, followed by option V. In this paper, the score of option III is higher than that of option V, and the reasoning behind this result can be divided into four factors, including technology, economy, environment, and society. From the perspective of the technical indicators, the proportion of excess electricity is the lowest and the ease of installation is the highest, which means that option III has a better technical performance and simpler construction techniques. From the perspective of the economy, option III's initial cost and COE are lower than those of other distributed renewable energy system options because the cost of hydrogen storage systems is higher than that of battery storage systems. From the perspective of the environment, the area requirement and environmental impact of option III perform better than do those of option V. On one hand, the requirement of land surface area of PV panels is high, and the number of PV panels in option III is higher than that in option V. On the other hand, option V's configuration includes wind turbines, which have greater a noise effect on the environment than option III. From the perspective of society, option III has lower technical requirements and is easier to implement because people in the research area have a better understanding of battery energy storage technology than hydrogen energy storage.

Conclusions
Promoting renewable energy utilization is challenging because it is related not only to local renewable resources but also to local policies and economic development levels. The XND area, as

Conclusions
Promoting renewable energy utilization is challenging because it is related not only to local renewable resources but also to local policies and economic development levels. The XND area, as an area designated for green and low-carbon development, needs to fully utilize local renewable energy resources. To achieve this target in a short time, a method that combines the HOMER model and the TOPSIS method was proposed in this study, which is used to quantitatively evaluate the feasible renewable energy allocation strategy in the XND area.
The optimal configuration for each scenario, which are regarded as the candidate options to be assessed, can be determined through the simulation and optimization results of the HOMER model. Based on the TOPSIS method results, the optimal HRES options are option III and option V, which are composed of a 134.40 MW PV array, a 24 MW diesel generator, a 38.43 MW converter, and 19.6 MW batteries and a 168.43 MW PV array, 28 kW wind turbines, 15 MW diesel generators, a 39.20 MW converter, 75.6 MW batteries, a 20 MW fuel cell, a 10 MW electrolyser, and a 10 t hydrogen tank, respectively. To perform a more comprehensive assessment of the scenarios, we conducted technical, economic, environmental and social analyses of these two optimal scenarios. First, the renewable fractions of options III and V are 85.8% and 99.7%, respectively, which means that most of the electricity is generated from renewable energy. In addition, the excess electricity portions in option III and V are 6.77% and 25.1%, respectively. The excess electricity in option V can be improved by the development of hydrogen storage technology. Apart from the system configuration and technical analysis, this study also presents economic analysis results based on three indicators, including the initial cost, O&M cost and COE. The initial cost, O&M cost and COE of options III and V are 253 and 355 M$, 9.35 and 2.24 M$ and 0.203 and 0.209 $/kWh, respectively. Moreover, we also conduct an environmental analysis of each option. The CO 2 emissions in option V are 345 t/yr, and this value is the minimum value among all the candidate options. The carbon dioxide emission reduction in option III ranks second, at 14,473 t/yr. In general, these two options have good environmental benefits and play important roles in improving the ecological environment and promoting carbon emission reduction.
Optimal HRES configuration selection in the XND area is realized by the quantitative evaluation method in combination with the HOMER model and the TOPSIS method. However, there are still limitations and uncertainties. On one hand, the limitations lie in the low applicability level of this method because different regions have different policies. This method is suitable for economically developed regions in which distributed renewable and hydrogen energy resources urgently need to be developed. On the other hand, uncertainties are embodied in the strong subjectivity of this approach because different experts have different opinions, which lead to different weights being applied in the TOPSIS method. Although this study adopts the AHP-CRITIC weighting method, there is still strong subjectivity because there are not enough samples of expert grading tables and there are few qualified experts. To alleviate the influence of subjectivity, other multi-objective evaluation methods will be applied in a future study. Despite these shortcomings, the research results provide a new approach for decision-makers and provide a useful reference for local renewable energy planning.