Optimal Distribution of Renewable Energy Systems Considering Aging and Long-Term Weather E ﬀ ect in Net-Zero Energy Building Design

: Generation system interruptions in net-zero energy buildings (NZEBs) may result in missing the net-zero targets by a great margin. Consequently, it is signiﬁcant to incorporate a realistic reliability model for renewable energy systems (RESs) that considers aging and long-term weather conditions. This study proposed a robust design optimization method that deals with the selection of RES to achieve NZEB. Di ﬀ erent case studies were evaluated: 1. Deterministic approach; 2. Markov chain-based reliability without the aging e ﬀ ect; 3. Markov chain-based reliability with the aging e ﬀ ect. The results showed that the optimal sizes of RES, considering the aging e ﬀ ect, were much larger than the other two cases based on the annual energy balance. Moreover, the consideration of the aging e ﬀ ect on the reliability assessment of the generation system for NZEB opens a pathway for a more robust and economic design of RES.


Introduction
Net-zero energy buildings (NZEBs) have been widely recognized as a promising way to solve the energy crisis and environmental problems. The topic has become so significant that nowadays a range of supporting policies and a financial incentive have been provided for the development of the NZEB in different countries. In Europe, the goal of "nearly net-zero energy buildings" has been established for all the new buildings from 2020 [1]. In the USA, a net-zero energy target for all new residential buildings by 2020 and for all new commercial buildings by 2030 has been set by the California Public Utilities Commission [2]. Substantial research has been conducted on how to design a NZEB in an economic and efficient way. In general, three approaches have been identified to be key factors for achieving NZEB: (1) passive design strategy (e.g., building orientation, cooling/heating strategy) [3,4]; (2) energy efficiency technologies (e.g., HVAC, lighting, and appliances) [5,6]; (3) energy production technologies (e.g., combined cooling, heat and power, battery, photovoltaic panel, and wind turbine) [7,8]. A schematic diagram of the approaches is shown in Figure 1. The literature presents many models on NZEB, for instance, the investigation conducted by Thalfeldt et al. [3] on optimal cost solutions regarding building facade solutions (i.e., the ratio of window-to-wall, external wall insulation, window properties, and shading) for nearly zero energy buildings. Wu et al. [5] investigated HVAC technologies for a residential NZEB in order to achieve the target of energy-saving, cost-effective, and comfortable indoor environment. In the study of Lu et al. [7], the selection of RES was conducted using a single-objective method based on genetic algorithm and multi-objective optimization method based on non-dominated sorting genetic algorithm (NSGA-II) for low/zero energy buildings, respectively. Figure 1 presents a classification of the main approaches to achieve NZEB.
Considering deterministic methods of analysis may not predict the levelized cost of energy accurately, Maheri [10] proposed two algorithms to find the optimum margin of safety concerning the system reliability by the Monte Carlo simulation for standalone hybrid PV-wind diesel systems. In the study of Lu et al. [11], a penalty cost function was introduced for the NZEB design, and the potential benefits for investors were investigated by comparing with the case under no financial support. Many of the studies have considered the effect of uncertain parameters on system selection for NZEB and its performance. The uncertainties were mainly focused on building energy consumption (physical parameters, design parameters, scenario parameters) [9,12,13] and renewable sources (i.e., solar radiation, wind velocity, and biodiesel price) [11,14,15]. As indicated by Huang et al. [16], most of the present studies related to NZEB ignore component degradation and maintenance. Therefore, they proposed a robust design of NZEB systems by considering system degradation (i.e., HVAC system, generation system, and energy storage system) and maintenance. However, the effect of the reliability of generation units (considering failure and the aging effect) on the selection of RES for NZEB still needs to be investigated.
Sustainability 2020, 12, x FOR PEER REVIEW 2 of 21 presents many models on NZEB, for instance, the investigation conducted by Thalfeldt et al. [3] on optimal cost solutions regarding building facade solutions (i.e., the ratio of window-to-wall, external wall insulation, window properties, and shading) for nearly zero energy buildings. Wu et al. [5] investigated HVAC technologies for a residential NZEB in order to achieve the target of energysaving, cost-effective, and comfortable indoor environment. In the study of Lu et al. [7], the selection of RES was conducted using a single-objective method based on genetic algorithm and multiobjective optimization method based on non-dominated sorting genetic algorithm (NSGA-II) for low/zero energy buildings, respectively. Figure 1 presents a classification of the main approaches to achieve NZEB. The carbon reduction targets and integration of new technologies have shifted the focus of the research mainly on the installation of RES in NZEB. However, how to select an appropriate generation system to achieve NZEB is an interesting field that attracts many researchers. A summary of design methodologies and generation system combinations is listed in Table 1. For instance, Zhang et al. [9] proposed a response-surface-model-based system sizing method that aims to reduce the number of the Monte Carlo simulation, and it was confirmed to work satisfactorily for the energy performance of NZEB.
Considering deterministic methods of analysis may not predict the levelized cost of energy accurately, Maheri [10] proposed two algorithms to find the optimum margin of safety concerning the system reliability by the Monte Carlo simulation for standalone hybrid PV-wind diesel systems. In the study of Lu et al. [11], a penalty cost function was introduced for the NZEB design, and the The carbon reduction targets and integration of new technologies have shifted the focus of the research mainly on the installation of RES in NZEB. However, how to select an appropriate generation system to achieve NZEB is an interesting field that attracts many researchers. A summary of design methodologies and generation system combinations is listed in Table 1. For instance, Zhang et al. [9] proposed a response-surface-model-based system sizing method that aims to reduce the number of the Monte Carlo simulation, and it was confirmed to work satisfactorily for the energy performance of NZEB. Table 1. Summary of the design optimization of generation systems for NZEB. (WT-Wind turbine, BDG-Biodiesel generator, TES-Thermal energy system, EES-Electrical energy storage, SDHW-solar domestic hot water, COP-Coefficient of performance).
The consideration of generation units' aging effect and quantitative reliability assessment is considered to provide a robust system and more comprehensive information about generation systems for users' selection. A question arises naturally: how reliable will the NZEB be facing the failure and the degradation of the generation component?
The main contributions of the paper can be summarized in the following aspects: (1) an innovative reliability assessment approach is applied to evaluate the reliability of RES components in terms of the aging effect for the design of NZEB; (2) the impact of the number of years considered for net energy balance is firstly proposed and investigated, which can greatly affect the size or RES components and provide basic guidance for NZEB design principles; (3) the incorporation of the model that shows the relationship between mismatch ratio (the difference between energy demand and energy generation) and the probabilities of being NZEB is identified and updated considering both the main uncertainty parameters and generation system degradation.

Problem Identification
The current studies on NZEB are generally based on the assumption of ideal renewable energy generation, ignoring the impact of aging on the reliability of generation systems. In such cases, the obtained design option (C1) of RES has a lower probability to achieve the expected energy balance during the real application, as shown in Figure 2, where the y-axis indicates the probability to achieve net energy balance and the x-axis is the mismatch ratio. A simple consideration of the reliability assessment of RES is to give a constant availability value for them, which ignores the degradation of the generation system due to aging. Then, the obtained design option (C2) of RES may have a higher probability than that of design option C1, however, it is still lower than the expected value. A more accurate method is to consider the failures of generation systems, where the tendency of failure will increase concerning aging impact. In the third situation, a much larger size of generation systems (C3) may be required whilst it produces an accurate reliability assessment for generation adequacy. In this study, the effect of reliability of generation systems was identified and quantified on the target of net-zero energy balance design.

Problem Identification
The current studies on NZEB are generally based on the assumption of ideal renewable energy generation, ignoring the impact of aging on the reliability of generation systems. In such cases, the obtained design option (C1) of RES has a lower probability to achieve the expected energy balance during the real application, as shown in Figure 2, where the y-axis indicates the probability to achieve net energy balance and the x-axis is the mismatch ratio. A simple consideration of the reliability assessment of RES is to give a constant availability value for them, which ignores the degradation of the generation system due to aging. Then, the obtained design option (C2) of RES may have a higher probability than that of design option C1, however, it is still lower than the expected value. A more accurate method is to consider the failures of generation systems, where the tendency of failure will increase concerning aging impact. In the third situation, a much larger size of generation systems (C3) may be required whilst it produces an accurate reliability assessment for generation adequacy. In this study, the effect of reliability of generation systems was identified and quantified on the target of net-zero energy balance design. Mismatch ratio ( ) is an index used to identify the difference between on-site energy generation ( ) and building energy demand ( ) for the selected design system, and it can be calculated in the form as given in (1). The probability of being NZEB is an index that can be assessed by using the ratio between the number of satisfying years and the total years' simulation [14] as given in (2).
Optimal selection of RES for NZEB is conducted for 100 operation years, and a diagram of the design optimization procedure considering RES reliability with aging features is shown in Figure 3. In the first step, electricity generation models (i.e., RES models) and building energy consumption model (energy system models and cooling load model) were developed. By using a Markov chain, a reliability assessment of renewable energy systems in each year was conducted.
Probability to achieve net energy balance Mismatch ratio (ε) is an index used to identify the difference between on-site energy generation (P gen ) and building energy demand (P dem ) for the selected design system, and it can be calculated in the form as given in (1). The probability of being NZEB is an index that can be assessed by using the ratio between the number of satisfying years and the total years' simulation [14] as given in (2). ε = (P gen − P dem )/P dem (1) Optimal selection of RES for NZEB is conducted for 100 operation years, and a diagram of the design optimization procedure considering RES reliability with aging features is shown in Figure 3. In the first step, electricity generation models (i.e., RES models) and building energy consumption model (energy system models and cooling load model) were developed. By using a Markov chain, a reliability assessment of renewable energy systems in each year was conducted.
Subsequently, the profile of 100 years' parameters (e.g., solar radiation, wind velocity, cooling load, and building other load, etc.) was generated by the employment of Monte Carlo simulation in combination with the profile of a standard year. For each selected RES size, a mismatch ratio was calculated based on energy generation and energy consumption. The capacity of each RES combination needs to further identify whether it meets the evaluation criterion of NZEB (e.g., annual energy balance, 100-years energy balance). In the third step, the total cost (TC) and grid interaction index (GII) were obtained for those satisfied RES sizes. Exhaustive search approach is applied in this study for design optimization. In the last step, the optimal option of RES and the corresponding NZEB performance were sorted according to the normalized total cost (TC n ) or normalized grid interaction index (GII n ).  Subsequently, the profile of 100 years' parameters (e.g., solar radiation, wind velocity, cooling load, and building other load, etc.) was generated by the employment of Monte Carlo simulation in combination with the profile of a standard year. For each selected RES size, a mismatch ratio was calculated based on energy generation and energy consumption. The capacity of each RES combination needs to further identify whether it meets the evaluation criterion of NZEB (e.g., annual energy balance, 100-years energy balance). In the third step, the total cost ( ) and grid interaction index ( ) were obtained for those satisfied RES sizes. Exhaustive search approach is applied in this study for design optimization. In the last step, the optimal option of RES and the corresponding NZEB performance were sorted according to the normalized total cost ( ) or normalized grid interaction index ( ).

System Generation Reliability Model
The reliability modelling requires the implementation of Markov chains, which is defined as a representation of all possible states in a space state diagram [23]. Each state is connected between them by variables called transition rates. The states employed to define the reliability model of each generation unit are: "1" operating, "2" repairable failure, and "3" obsolescence. The transitions rates are defined as failure , repair , and degradation [24].

System Generation Reliability Model
The reliability modelling requires the implementation of Markov chains, which is defined as a representation of all possible states in a space state diagram [23]. Each state is connected between them by variables called transition rates. The states employed to define the reliability model of each generation unit are: "1" operating, "2" repairable failure, and "3" obsolescence. The transitions rates are defined as failure λ r , repair µ r , and degradation λ F [24].
In this paper, Case 2 is considered to be an alternating renewal process with two states, as presented in Figure 4a. This model disregards the aging effect and brings inaccuracies in a realistic reliability model of the unit generation systems. Therefore, Case 3 is proposed which employs the model presented by Alvarez-Alvarado and Jayaweera [24]. The model contemplates the aging effect by adding a third state called obsolescence, which is an absorbing state, as presented in Figure 4b. In this paper, Case 2 is considered to be an alternating renewal process with two states, as presented in Figure 4a. This model disregards the aging effect and brings inaccuracies in a realistic reliability model of the unit generation systems. Therefore, Case 3 is proposed which employs the model presented by Alvarez-Alvarado and Jayaweera [24]. The model contemplates the aging effect by adding a third state called obsolescence, which is an absorbing state, as presented in Figure 4b.  Later, the degradation rate can be obtained as a function of the , , and the lifetime of the component , that is [24]: Later, the degradation rate can be obtained as a function of the λ r , µ r , and the lifetime of the component T W , that is [24]: On the other hand, the probability of being in each state can be mathematically determined from the stochastic matrix of transition states H. This matrix is the infinitesimal generator of dimensions z × z, where z is the total number of states. The terms inside this matrix are defined as follows: d ij = λ ij (i j) and d ii = − i j λ ij . Therefore, the H matrices for the models presented in Figure 4a,b are as shown in (4) and (5), respectively.
Then, the probability vector of all possible states is determined by (6), where υ is the eigenvalues of H T , v is the eigenvectors of H T , and k is a constant given by the initial state; T indicates the transpose of the matrix [25,26].
Applying this criterion given in (4) and (5), the probabilities of being in each state are as presented in (7) and (8), respectively. case 2 : case 3 : A measure of reliability is the availability of the component, and it is defined as the sum of the probabilities of being in the operating state. For the models given in Figure 4a,b, there is just the operating state (state "1"), therefore, their availability is defined as presented in (9) and (10), respectively.
Then, considering the reliability data of PV, WT, and BDG listed in Table 2 and following the procedure stated above, the availability for each generation system is obtained as shown in Figure 5. In addition, Table 3 provides the availability ranges for the generation systems for each case study.   It is clear from Table 3 that without consideration of the aging factor, the reliability of the RES does not vary significantly. However, the availability shows a significant variation ranging from 0.40 to 0.95 for generation units. It is clear from Table 3 that without consideration of the aging factor, the reliability of the RES does not vary significantly. However, the availability shows a significant variation ranging from 0.40 to 0.95 for generation units.

Evaluation Indicator
In this study, the total cost of NZEB comprises of the initial cost (i.e., investment of RES), operation cost (i.e., electricity cost from the grid and the cost of biodiesel) during the building application stage, and the imposed cost (i.e., an additional investment cost due to reliability consideration). This is given in (11).
Grid interaction index (GII) is an indicator representing the average stress of a building on the power grid, and it is defined as the standard deviation of the building-grid interaction over the period as formed in (12). The building-grid interaction is calculated based on the ratio between net exported energy and the average energy demand in the building during a given period, as shown in (13).
f grid,i,T = P ex,i − P im,i t2 t1 P dem,i dt/T The objective in this study was to minimize the total cost or grid interaction index from the selected ranges of RES for NZEB, as indicated in (14), where TC n , GII n , and X are the normalized total cost, normalized grid interaction index, and the decision variables on the selected combination sizes of each generation unit, respectively. In order to achieve the target of net-zero energy building, the selected size of each generation unit should be within the set size ranges (15), the amount of energy generation during certain periods should be larger than or equal to the building energy demand (16), and the energy/cooling balance must be met in each hour (17). As shown in Figure 3, the exhaustive searching method will be applied to solve the optimization.
Objective function : f = TC n or GII n Subject to AX ≤ a

Building Features
The first zero-carbon building (ZCB) in Hong Kong, China, was used to identify the effect of reliability of RES on design option and corresponding performance. This is a three-story building (area = 1520 m 2 ) with the air-conditioned area of approximately 990 m 2 . Two types of generation systems, i.e., photovoltaic (PV) of 1015 m 2 and a biodiesel generator (BDG) with a rated power of 100 kW, are employed to provide energy for the grid-connected building. In this study, the wind turbine (WT) system was also assumed to be one of the potential generation systems, in addition to PV and BDG. The cooling load of the building was undertaken by three electric chillers, one absorption chiller driven by BDG, and the peak cooling load was about 163 kW. A schematic diagram of the energy systems in the studied building is shown in Figure 6.
The main key parameters for the energy system are listed in Table 4. The total number of PV, WT, and BDG size combination was 945 (Table 4: 21 × 9 × 5 = 945), and the lifetime of the three systems was set to be 25 years.
Four uncertain parameters, i.e., solar radiation, wind velocity, cooling load, and building other load (electricity usage, except for cooling load), were considered to be varying over a 100-years, and it was assumed that they follow uniform distributions. It was assumed that when the component reaches the obsolescence state (end lifetime), this is immediately replaced and the lifecycle starts again. The study was conducted for the next 100-years simulation to investigate the impact of uncertainty from energy demand, and energy generation can also be reflected.     The cooling system of the building is provided by three electric and one absorption chiller, as shown in (18). The absorption chiller ( ) is driven by BDG and is given priority to handle the cooling load, whereas the rest of the cooling load is undertaken by the electric chillers ( ), as shown in (19) and (20).
The total electricity consumption of the building ( ) consists of five parts, as shown in (21), where Pec is the power consumed by electric chillers, Ppump is the power consumed by cooling water pumps and chiller water pumps, Pct is the power consumed by cooling tower fans, is the power consumed by AHU (air handling unit) fans, and Pother is the power consumed by other appliances, such as lighting and socket outlet. The power supply system consists of on-site RES (i.e., PV, WT, and BDG) and the grid, as shown in (22). The hourly electricity consumption in the building should be The cooling system of the building is provided by three electric and one absorption chiller, as shown in (18). The absorption chiller (Q ac ) is driven by BDG and is given priority to handle the cooling load, whereas the rest of the cooling load is undertaken by the electric chillers (Q ec ), as shown in (19) and (20).
The total electricity consumption of the building (P dem ) consists of five parts, as shown in (21), where P ec is the power consumed by electric chillers, P pump is the power consumed by cooling water pumps and chiller water pumps, P ct is the power consumed by cooling tower fans, P f an is the power consumed by AHU (air handling unit) fans, and P other is the power consumed by other appliances, such as lighting and socket outlet. The power supply system consists of on-site RES (i.e., PV, WT, and BDG) and the grid, as shown in (22). The hourly electricity consumption in the building should be equal to the amount of supplied power by the sources given in (22). It is noted that heating load does not exist for the building in Hong Kong since Hong Kong belongs to the hot summer and warm winter zone. P t dem = P t ec + P t pump + P t ct + P t f an + P t other (21) P t supply = P t PV + P t WT + P t BDG + P t grid (22) Electricity consumption of the electric chillers (P ec ), pumps models (P pump ), cooling tower (P ct ), and AHU fan (P f an ) can be modelled in the following form, based on the empirical model and the on-site data in ZCB [7]: Energy generation in a PV system, WT system, and BDG system can be evaluated by using (27), (28), and (29), respectively. The power generation of BDG depends on the building cooling demand (26) since it is controlled following the cooling load. Furthermore, fuel consumption of BDG (F bio ) during the operation stage can be assessed by (30).
F bio = A G P BDG + B G P rated,BDG (30)

Results and Discussion
RES design was optimized for all three cases. The effect of reliability of generation systems can be identified by comparing the results of the three cases. The required RES size in Case 3, when the aging effect is considered, was expected to be larger than that in the other two cases. Furthermore, the number of years concerned for net energy balance evaluation may have a great impact on system selection and NZEB performance since the degeneration of RES due to the aging effect was considered. Therefore, both scenarios are investigated in the following section.

Performance Evaluation: Annual Energy Balance
In terms of the annual energy balance, the optimal design results based on a 100-years simulation are presented in Table 5. It was found that the optimal sizes for Case 3 were much larger than that for the ones in Case 1 and Case 2; for example, the combination of optimal RES consisted of WT of 160 kW, BDG of 20 kW, and PV of 200 m 2 in Case 1, whilst the PV size increased to be 400 m 2 in Case 2, and it increased to be 3,600 m 2 in Case 3. The main reason was that a slight reduction of availability was found in each generation unit (it was more than 0.95) in Case 2, whilst a large variation range of availability (between 0.4 and 0.95) was seen in Case 3, as shown in Table 3. In addition, the total cost (T Cn ) and grid interaction index (GII n ) in Case 3 were identified to be two times and three times, respectively, that of Case 1.
In order to have a better understanding of actual NZEB performance under the selected RES sizes in each year, a comparison of annual energy generation and annual mismatch ratio for the optimal options in the three cases is presented in Figure 7 (for T Cn ) and Figure 8 (for GII n ), respectively.
In terms of the total cost, the average values of annual energy generation and annual mismatch ratio were 209,775 kWh and −8.2% (Case 1), 223,214 kWh and −2.3% (Case 2), and 438,226 kWh and 91.9% (Case 3), respectively.      Figures 9 and 10 further present the average annual energy generation of the selected option and the corresponding probability to be NZEB when the aging effect is taken into consideration. It can be seen that the average annual energy generation was less than its energy consumption in Case 1 and Case 2, whilst it was about two times that of energy consumption in Case 3. It was also interesting to find that the probability to be NZEB during practical application was only 0.35 for Case 1 and 0.43 for Case 2 when the total cost was considered, while it was 0.2 for Case 1 and 0.17 for Case 2 when the grid interaction index was considered. Thus, the aging effect had a significant impact on system selection and expected target.   Figures 9 and 10 further present the average annual energy generation of the selected option and the corresponding probability to be NZEB when the aging effect is taken into consideration. It can be seen that the average annual energy generation was less than its energy consumption in Case 1 and Case 2, whilst it was about two times that of energy consumption in Case 3. It was also interesting to find that the probability to be NZEB during practical application was only 0.35 for Case 1 and 0.43 for Case 2 when the total cost was considered, while it was 0.2 for Case 1 and 0.17 for Case 2 when the grid interaction index was considered. Thus, the aging effect had a significant impact on system selection and expected target. Figures 9 and 10 further present the average annual energy generation of the selected option and the corresponding probability to be NZEB when the aging effect is taken into consideration. It can be seen that the average annual energy generation was less than its energy consumption in Case 1 and Case 2, whilst it was about two times that of energy consumption in Case 3. It was also interesting to find that the probability to be NZEB during practical application was only 0.35 for Case 1 and 0.43 for Case 2 when the total cost was considered, while it was 0.2 for Case 1 and 0.17 for Case 2 when the grid interaction index was considered. Thus, the aging effect had a significant impact on system selection and expected target.

Performance Evaluation: 100-years Energy Balance
In terms of the 100-years energy balance, the optimal design results are presented in Table 6. For the total cost, the selected options in Case 3 were larger than those in Case 1 and Case 2, but it was reduced from 3600 m2 under annual year balance concerned (Table 5) to 600 m2 under 100-years energy balance. The total cost did not show much reduction in Case 1 and Case 2, however, in Case 3, compared to the annual year balance, the cost was reduced to half. In addition, a remarkable reduction from 92% to 3% was seen in the mismatch ratio in Case 3, indicating a much smaller size was required for the 100-years energy balance. A similar tendency was found for the performance of the grid interaction index. A further comparison of annual energy generation and annual mismatch

Performance Evaluation: 100-Years Energy Balance
In terms of the 100-years energy balance, the optimal design results are presented in Table 6. For the total cost, the selected options in Case 3 were larger than those in Case 1 and Case 2, but it was reduced from 3600 m 2 under annual year balance concerned (Table 5) to 600 m 2 under 100-years energy balance. The total cost did not show much reduction in Case 1 and Case 2, however, in Case 3, compared to the annual year balance, the cost was reduced to half. In addition, a remarkable reduction from 92% to 3% was seen in the mismatch ratio in Case 3, indicating a much smaller size was required for the 100-years energy balance. A similar tendency was found for the performance of the grid interaction index. A further comparison of annual energy generation and annual mismatch ratio in the three cases under 100-years energy balance is presented, when the aging effect was taken into consideration, in Figure 11 (regarding TC n ) and Figure 12 (regarding GII n ).  In terms of the total cost, the average values of annual energy generation and annual mismatch ratio were 196,337 kWh and −14.1% (Case 1), 196,377 kWh and −14.1% (Case 2), and 236,652 kWh and 3.6% (Case 3), respectively. In terms of the grid interaction index, the average values of annual energy generation and annual mismatch ratio were 129,719 kWh and −39.3% (Case 1), 143,157 kWh and −33.0% (Case 2), and 218,019 kWh and 2.2% (Case 3), respectively. Following a similar pattern of the annual energy balance, the fluctuations of both parameters showed approximately periodic variations. The average annual energy generation of the selected option and the corresponding probability to be NZEB in the three cases were further compared when the aging effect was taken into consideration, as shown in Figures 13 and 14. Different from that under annual energy balance criteria, the required annual energy generation in Case 3 was much less than that in Case 1 and Case 2.  Year ɛ Figure 11. Annual performance of NZEB under optimal options with the aging effect (considering TC n ): (a) Annual energy generation; (b) Annual mismatch ratio.
The average annual energy generation of the selected option and the corresponding probability to be NZEB in the three cases were further compared when the aging effect was taken into consideration, as shown in Figures 13 and 14. Different from that under annual energy balance criteria, the required annual energy generation in Case 3 was much less than that in Case 1 and Case 2.
In terms of the total cost, the probability to be NZEB in Case 1 (0.3) was the same as that in Case 2 since they had the same design option, whilst it was about 0.5 in Case 3 during practical application. In terms of the grid interaction index, the probability to be NZEB during practical application was only 0.0, 0.06, and 0.46 for Case 1, Case 2, and Case 3, respectively. Thus, the number of years for net energy balance evaluation also had a significant impact on system selection and expected target. In terms of the total cost, the probability to be NZEB in Case 1 (0.3) was the same as that in Case 2 since they had the same design option, whilst it was about 0.5 in Case 3 during practical application. In terms of the grid interaction index, the probability to be NZEB during practical application was only 0.0, 0.06, and 0.46 for Case 1, Case 2, and Case 3, respectively. Thus, the number of years for net energy balance evaluation also had a significant impact on system selection and expected target.

Relationship Between Probabilities and Mismatch ratio
The mismatch ratio of the selected design option has a significant impact on the probability of being NZEB, as indicated in [8] and [12]. Therefore, the selected mismatch ratio can reflect its ability on the NZEB level achieved, which can provide a general perspective for designers. All of the size

Relationship Between Probabilities and Mismatch ratio
The mismatch ratio of the selected design option has a significant impact on the probability of being NZEB, as indicated in [8] and [12]. Therefore, the selected mismatch ratio can reflect its ability on the NZEB level achieved, which can provide a general perspective for designers. All of the size

Relationship between Probabilities and Mismatch Ratio
The mismatch ratio of the selected design option has a significant impact on the probability of being NZEB, as indicated in [8] and [12]. Therefore, the selected mismatch ratio can reflect its ability on the NZEB level achieved, which can provide a general perspective for designers. All of the size combinations (945) were investigated in a 100-years operation, and it was interesting to obtain three fitting formulas for the studied three cases, as shown in Figure 15. The coefficients of determination R 2 of the three fitting formula were all above 0.98, indicating a good relationship between the mismatch ratio and the probabilities.
A detailed comparison of the probability of the three cases is listed in Table 7. For all the cases, the probability to be NZEB was 0% and around 45% under a mismatch ratio below -0.4 and about 0.0, respectively. In general, a similar trend of probability was found in Case 1 and Case 2. When the mismatch ratio was more than 0.6, the probability in Case 1 and Case 2 could hit the target of 100%, whilst it was a mismatch ratio of over 1.0 that could ensure the probability of 100% in Case 3. In addition, it seems that the probability in Case 3 was higher than that in Case 1 and Case 2 when the mismatch ratio was between -0.4 and 0.0. However, it was reversed when the mismatch ratio was between 0.0 and 0.8. combinations (945) were investigated in a 100-years operation, and it was interesting to obtain three fitting formulas for the studied three cases, as shown in Figure 15. The coefficients of determination R 2 of the three fitting formula were all above 0.98, indicating a good relationship between the mismatch ratio and the probabilities. Case 1: y1 = 36.914x 5 − 0.8983x 4 − 16.554x 3 + 0.3929x 2 + 2.9041x + 0.4636; −0.5 < x < 0.5 (R² = 0.988); Case 2: y2 = 37.906x 5 − 1.9154x 4 − 16.716x 3 + 0.6118x 2 + 2.9024x + 0.4563; −0.5 < x < 0.5 (R² = 0.988); Case 3: y3 = -0.5495x 5 + 1.5527x 4 − 1.2653x 3 − 0.4425x 2 + 1.2506x + 0.4552; −0.5 < x < 1.0 (R² = 0.997). A detailed comparison of the probability of the three cases is listed in Table 7. For all the cases, the probability to be NZEB was 0% and around 45% under a mismatch ratio below -0.4 and about 0.0, respectively. In general, a similar trend of probability was found in Case 1 and Case 2. When the mismatch ratio was more than 0.6, the probability in Case 1 and Case 2 could hit the target of 100%, whilst it was a mismatch ratio of over 1.0 that could ensure the probability of 100% in Case 3. In addition, it seems that the probability in Case 3 was higher than that in Case 1 and Case 2 when the mismatch ratio was between -0.4 and 0.0. However, it was reversed when the mismatch ratio was between 0.0 and 0.8.