Analysis of Primary Energy Factors from Photovoltaic Systems for a Nearly Zero Energy Building (NZEB): A Case Study in Lithuania

: Following a new climate and energy plan, the European Union (EU) gives big attention to energy savings. The overall assessment of energy saving measures is very important. Thus, it is crucial to estimate in a proper way the primary energy factor, which is used in calculations of primary energy consumption from renewable energy (RE) sources in a Nearly Zero Energy Building (NZEB). The conduced studies of the literature and national regulations showed that di ﬀ erent methods to determine energy from photovoltaic (PV) systems are used. The aim of this paper is to evaluate the primary energy factors of energy from photovoltaics and determine the average value. To achieve this aim, the data of 30 photovoltaic systems from Lithuania were analyzed. The results show a 35% diversiﬁcation in the values of non-renewable primary energy factor, depending on the PV systems’ capacities, with the average on a level of 1.038.


Overview on Renewable Energy
At the turn of the 20th and 21st centuries, increased impacts on the environment were observed, and thus studies on sustainable development as well as actions aimed to mitigate climate change became more popular [1]. The United Nations Framework on Climate Change (UNFCCC) in 1992 set the background for the limitation of greenhouse gas emissions, as actions in this area were found to be essential for the prevention of doom-laden environmental consequences [2]. Additionally, the United Nations proposed the 17 Sustainable Development Goals in 2015, including topics related to ways of reducing greenhouse gas emissions [3]. Most attention was focused on the energy field and especially on the building sector. Existing buildings consume about 40% of the EU's total energy and cause above one third of its CO 2 emissions [4]. This is despite the progress made in implementing energy efficiency policies under national plans and actions. Thus, Member States do their best to meet the requirements of Directive 2006/32/EC of the European Parliament and of the council of 5 April 2006 on energy end-use efficiency. However, still the union has hardly succeeded in achieving its energy grids [20][21][22][23]. It means that the supplied energy from renewable energy sources to the NZEB boundary is a mixture of renewable and non-renewable energy. Figure 1 demonstrates the supply chain, starting from the primary energy from renewable sources and ending in its usage in a building. All the direct and indirect energy transformations that are part of the electricity generation process were indicated, including the energy spent on extraction, transportation, and transformation, and it takes into account the real share of renewable energy sources.
Energies 2020, 13, x FOR PEER REVIEW 3 of 15 Figure 1 demonstrates the supply chain, starting from the primary energy from renewable sources and ending in its usage in a building. All the direct and indirect energy transformations that are part of the electricity generation process were indicated, including the energy spent on extraction, transportation, and transformation, and it takes into account the real share of renewable energy sources. Renewable electricity generation and consumption differ due to climate conditions, the season of the year, and the time of day. For example, in the case of photovoltaic systems connected to the local electricity grid, the electricity produced during the day can be used immediately in the building or transmitted to the electricity grid. In the evening, when the photovoltaic system is not working, electricity is supplied from the mains. No battery is required for this type of system [24]. For these reasons, the lack of electricity available at the peak of electricity consumption can appear. In addition, conversion devices consume renewable energy during periods of low electricity consumption, or they can even discharge electricity to the grid while the needs are high. In those cases, conversion devices normally use additional non-renewable energy (for instance, any fuels) [25,26].
Directive 2018/2002/EU (indent 21) states the importance of considering all the steps in the energy chain while calculating the energy savings in order to multiply the potential of savings in electricity transmission and distribution [9]. The Primary Energy Factor (PEF), described also as the conversion factor, helps to evaluate the primary energy consumption, including the chain of energy generation based on the final energy consumption data ( Figure 1). It shows the amount of primary energy used to produce a unit of electricity or, in specific situations, the unitary useable thermal energy.

Evaluation of Primary Energy Factor
Directive 2018/2002/EU (indent 40) provides three methods for setting the PEF value: • "physical energy content"-a method used for the generation of nuclear electricity and heat; • "technical conversion efficiency"-a method used for the production of electricity and heat from biomass and different fossil fuels; • "total primary energy"-a method that is used for non-combustible renewable energy.
The basic guidelines of the "total primary energy" method are provided in EN 15603:2014 [27]. It suggests that the PEF of electricity can be obtained when we divide the raw primary energy demand of electricity generation by the electricity produced, in cases when only the generation of energy is taken into account. The basic idea of the primary energy balance performed in EN 15603:2014 [27] is that each energy flow crossing the assessment boundary is characterized by the total associated primary energy and is given by sum of the non-renewable PEF (fP,nren) and renewable PEF (fP,ren) ( Figure 2). Renewable electricity generation and consumption differ due to climate conditions, the season of the year, and the time of day. For example, in the case of photovoltaic systems connected to the local electricity grid, the electricity produced during the day can be used immediately in the building or transmitted to the electricity grid. In the evening, when the photovoltaic system is not working, electricity is supplied from the mains. No battery is required for this type of system [24]. For these reasons, the lack of electricity available at the peak of electricity consumption can appear. In addition, conversion devices consume renewable energy during periods of low electricity consumption, or they can even discharge electricity to the grid while the needs are high. In those cases, conversion devices normally use additional non-renewable energy (for instance, any fuels) [25,26].
Directive 2018/2002/EU (indent 21) states the importance of considering all the steps in the energy chain while calculating the energy savings in order to multiply the potential of savings in electricity transmission and distribution [9]. The Primary Energy Factor (PEF), described also as the conversion factor, helps to evaluate the primary energy consumption, including the chain of energy generation based on the final energy consumption data ( Figure 1). It shows the amount of primary energy used to produce a unit of electricity or, in specific situations, the unitary useable thermal energy.

Evaluation of Primary Energy Factor
Directive 2018/2002/EU (indent 40) provides three methods for setting the PEF value: • "physical energy content"-a method used for the generation of nuclear electricity and heat; • "technical conversion efficiency"-a method used for the production of electricity and heat from biomass and different fossil fuels; • "total primary energy"-a method that is used for non-combustible renewable energy.
The basic guidelines of the "total primary energy" method are provided in EN 15603:2014 [27]. It suggests that the PEF of electricity can be obtained when we divide the raw primary energy demand of electricity generation by the electricity produced, in cases when only the generation of energy is taken into account. The basic idea of the primary energy balance performed in EN 15603:2014 [27] is that each energy flow crossing the assessment boundary is characterized by the total associated primary energy and is given by sum of the non-renewable PEF (f P,nren ) and renewable PEF (f P,ren ) ( Figure 2).
According to Directive 2018/2002/EU (indent 40), 100% conversion efficiencies are assumed in the case of non-combustible renewables, and the PEF (f P,ren ) is 1 for all energy sources for the system boundaries. However, the PEF value of non-renewable energy (f P,nren ), which includes transmission and distribution losses, is cannot easily determined. Depending on the type of primary energy and regarding the calculation method used, the primary energy factor may differ significantly for the same renewable energy source [24,28]. The analysis of the latest scientific studies [29][30][31][32][33] showed that primary energy is evaluated by different methods: • "Zero-equivalent method". This method does not evaluate the electrical and thermal energy production from renewables, thus, in this case, the total primary energy factor is considered as ƒ P,tot = 0 [29]. • "Direct equivalent". In this case, the electrical and thermal energy production from non-fossil renewable energy as well as nuclear energy sources are taken into consideration, thus ƒ P,tot = 1 [30]. • "Amount of physical energy". This method evaluates the primary form of energy obtained in the generation process and sets the factor as ƒ P,tot = 1 [31].

•
The "Alternate" method takes into account the primary form of energy. It is included in the statistical energy balance before the conversion to the secondary or tertiary form of energy, so the factor is relatively high and amounts to ƒ P,tot = 2.5 [32,33].

•
The "Effectiveness of technical conversion" method. It considers the entire energy generation chain, with the assumption of independent estimations of both renewable and non-renewable energy, resulting in a non-renewable primary energy factor ƒ P,nren = 0.032 and a renewable primary energy factor ƒ P,ren = 2.5 [32,33]. • "Amount of physical energy". This method evaluates the primary form of energy from the generation process and recommends the values of factors ƒ P,nren = 0.032 and ƒ P,ren = 1 [31]. According to Directive 2018/2002/EU (indent 40), 100% conversion efficiencies are assumed in the case of non-combustible renewables, and the PEF (fP,ren) is 1 for all energy sources for the system boundaries. However, the PEF value of non-renewable energy (fP,nren), which includes transmission and distribution losses, is cannot easily determined. Depending on the type of primary energy and regarding the calculation method used, the primary energy factor may differ significantly for the same renewable energy source [24,28]. The analysis of the latest scientific studies [29][30][31][32][33] showed that primary energy is evaluated by different methods: • "Zero-equivalent method". This method does not evaluate the electrical and thermal energy production from renewables, thus, in this case, the total primary energy factor is considered as ƒP,tot = 0 [29]. • "Direct equivalent". In this case, the electrical and thermal energy production from non-fossil renewable energy as well as nuclear energy sources are taken into consideration, thus ƒP,tot = 1 [30]. • "Amount of physical energy". This method evaluates the primary form of energy obtained in the generation process and sets the factor as ƒP,tot = 1 [31].

•
The "Alternate" method takes into account the primary form of energy. It is included in the statistical energy balance before the conversion to the secondary or tertiary form of energy, so the factor is relatively high and amounts to ƒP,tot = 2.5 [32;33].

•
The "Effectiveness of technical conversion" method. It considers the entire energy generation chain, with the assumption of independent estimations of both renewable and non-renewable energy, resulting in a non-renewable primary energy factor ƒP,nren = 0.032 and a renewable primary energy factor ƒP,ren = 2.5 [32;33]. • "Amount of physical energy". This method evaluates the primary form of energy from the The analysis shows that different methodologies use various PEF values for the evaluation of photovoltaic systems, thus the results are not comparable or useful for general studies. The analysis of PEF calculation options for electricity performed by Esser and Sensfuss [25] showed that the actually used methodology was outdated and needed a review. The PEF set at a level of 2.5 reflected the old data for the European power system, while it does not refer to the current energy market with a significant share of renewable sources in the power generation balance. As noted, in 2015 the average renewable energy sources share for electricity in the EU was 28%, and it is increasing continuously to meet a projected 45% level by 2030, based on the EU renewable energy sources' target.
The outcomes of the literature review are that the information about the setting of PEF values is not sufficient; the PEF values varied for the same type of energy, and several diversified methodologies were used for calculations. It is crucial to set a proper PEF value, as just accurate data on photovoltaic renewable (ƒ P,ren ) and non-renewable (ƒ P,nren ) primary factors could allow the objective estimation of the amount of renewable and non-renewable primary energy consumed in NZEB buildings. Therefore, Energies 2020, 13, 4099 5 of 14 this paper's aim is to investigate and evaluate the primary factors. The proposal of this work is also to estimate the average value of the primary energy factor as the energy factor of energy from photovoltaic systems, based on the data of the produced and consumed energy, as it intends to investigate its dependence on the capacity of PV systems. A detailed case study of the Lithuanian photovoltaic energy could be later used as a template for a similar analysis to be carried out for other national markets, and the foundlings can be relevant to policy makers.
The structure of this work consists of four sections: Methodology, Results, Discussion, and Conclusion. The research object, climate data, and primary energy calculation methodology are provided in the Methodology section. The section on the results provides a concise description of the experimental results and their interpretation. In the section on the discussion, the results are provided to be interpreted in a context of other studies and this working hypothesis. The study findings are discussed in the EU law and recommendations background. The final and significant outcomes are provided in the section on the conclusion.

Research Object
The statistical data from the year 2018 suggests that the total number of photovoltaic (PV) systems in Lithuania was 3050, in which the sum of the installed capacity was 88 MW. The amount of total produced energy was 67 GWh.
Data for the investigation were collected from 30 PV systems operating in Lithuania by interviewing PV system owners/operators and by taking into account the outcomes from the reports of electricity transmission system operators in the country. The data of the produced and consumed energy were collected over the period of 2014-2018. The PV systems were divided into groups of different capacities: group A > 20 (kW); group 20 > B > 10 (kW); and group C < 10 (kW). The main characteristics of the investigated PV systems are presented in Table 1.

Climate Data
Lithuania is in the south-western sub-area of the Atlantic continental zone. Lithuania is between the parallel of the 54 • and 56 • north latitude. Thus, in June, the height of the sun is 58.5 • ; in December, it is 11.5 • ; and on the equinox of spring and autumn, it is 35 • . At noon on a sunny day, at the indicated times of the year, the Lithuanian surface receives, respectively, 86%, 21%, and 57% of solar energy, regarding the average amount of energy on Earth above the sea surface. Respectively, the length of the day varies from 17.3 h in midsummer to 12 h on the equinoxes and up to 7.2 h in midwinter.
The average annual solar radiation impinging onto the horizontal surface of Lithuania reaches approximately 1000 kWh/m 2 . Sunlight lasts the longest in the western part (1840-1900 h/year) and shortens while moving towards the east (1700 h/year) (Figure 3).

Climate Data
Lithuania is in the south-western sub-area of the Atlantic continental zone. Lithuania is between the parallel of the 54° and 56° north latitude. Thus, in June, the height of the sun is 58.5°; in December, it is 11.5°; and on the equinox of spring and autumn, it is 35°. At noon on a sunny day, at the indicated times of the year, the Lithuanian surface receives, respectively, 86%, 21%, and 57% of solar energy, regarding the average amount of energy on Earth above the sea surface. Respectively, the length of the day varies from 17.3 h in midsummer to 12 h on the equinoxes and up to 7.2 h in midwinter.
The average annual solar radiation impinging onto the horizontal surface of Lithuania reaches approximately 1000 kWh/m 2 . Sunlight lasts the longest in the western part (1840-1900 h/year) and shortens while moving towards the east (1700 h/year) (Figure 3). The generative efficiency of photovoltaic power stations is directly dependent on solar radiation. In order to use solar radiation to its maximum under the Lithuanian climate conditions, photovoltaic modules must be directed towards the south a with 30°-60° lean angle from the horizontal plane. Such an orientation of one square meter of the area of a photovoltaic power station generates 130-160 kWh electrical energy annually under the Lithuanian climate conditions.

Primary Energy Calculation Methodology
The total primary energy factor may be calculated using the Equation (1)   The generative efficiency of photovoltaic power stations is directly dependent on solar radiation. In order to use solar radiation to its maximum under the Lithuanian climate conditions, photovoltaic modules must be directed towards the south a with 30 • -60 • lean angle from the horizontal plane. Such an orientation of one square meter of the area of a photovoltaic power station generates 130-160 kWh electrical energy annually under the Lithuanian climate conditions.

Primary Energy Calculation Methodology
The total primary energy factor may be calculated using the Equation (1) where ƒ P,tot is the total primary energy factor; ƒ P,nren is the non-renewable primary energy factor; ƒ P,ren is the renewable primary energy factor. The analysis is based on the assumption that all of the energy supplied to the building is attributable to renewable energy because it is made from renewable PV energy. Accordingly, the renewable primary energy factor ƒ P,ren is equal to 1. The value of the efficiency of electricity production η el and the losses of electricity transportation e T were evaluated by the calculations. The building normative documents of Lithuania for building electricity consumption calculations indicate that e T = 0 and η el = 2.8 [35]. The value of the primary non-renewable energy factor f P,nren produced by the PV system is given by the Equation (2), which is provided in EN 15603: 2014 [27]: where E c is the amount of consumed energy in kWh/year; E p is the amount of produced energy in kWh/year; η el is the coefficient of the efficiency of electricity production (η el = 2.8); e T is the coefficient of the losses of electricity transportation (e T = 0). The calculation of the value of PEF of the electricity produced by PV systems was performed according to the data provided by the power systems ( Table 2). The descriptive statistics were applied in this study. In order to estimate reliably the mean of the range µ, the sample mean x was estimated (Equation (3)): where: n is the number of items in a sample; j is the variable (j = 1,2, . . . n).
The mean of the range µ is calculated according to the mathematical expression (Equation (4)): where: N is the number of items in the population; j is the variable (j = 1,2, . . . n).

Results
In order to determine the PEF, first of all the overall balance of the generated and consumed electric power in the analyzed PV systems was investigated. The results of the analysis of the average monthly distribution of the produced and consumed electric power in A, B, and C groups of the studied PV systems are presented in Figure 4.
The average annual consumption percentage value from the total produced electricity amount of the PV systems of A group equals 0.7%. In the case of the PV systems of group B, the average annual consumption percentage value from the total produced electricity amounts to 0.9%. Finally, the average annual consumption percentage value for the C group of PV systems reaches 2.2%. Thus, the obtained results suggest that there is no relation between the produced and consumed electrical energy in the PV systems. Figure 5 shows the probability density function of the consumed electrical energy of PV systems, which has a normal Gauss distribution. That means that the electricity consumption is almost constant each month. The fluctuations are small and can be explained by hardware factors. The obtained analysis results allow us to reliably determine the PEF values.    As was mentioned before, the PEF or fP,tot (they are the same) values consist of photovoltaic renewable (ƒP,ren) and non-renewable (ƒP,nren) primary factors. The renewable primary energy factor ƒP,ren is equal to 1, according to the directive [18]. Hence, the non-renewable primary energy factor fP,nren was determined after the assignment of analyzed PV systems to groups of different capacities (A, B, and C). The results of the calculated fP,nren factors of the studied PV systems are presented in Figure 6.    As was mentioned before, the PEF or fP,tot (they are the same) values consist of photovoltaic renewable (ƒP,ren) and non-renewable (ƒP,nren) primary factors. The renewable primary energy factor ƒP,ren is equal to 1, according to the directive [18]. Hence, the non-renewable primary energy factor fP,nren was determined after the assignment of analyzed PV systems to groups of different capacities (A, B, and C). The results of the calculated fP,nren factors of the studied PV systems are presented in Figure 6. As was mentioned before, the PEF or f P,tot (they are the same) values consist of photovoltaic renewable (ƒ P,ren ) and non-renewable (ƒ P,nren ) primary factors. The renewable primary energy factor ƒ P,ren is equal to 1, according to the directive [18]. Hence, the non-renewable primary energy factor f P,nren was determined after the assignment of analyzed PV systems to groups of different capacities (A, B, and C). The results of the calculated f P,nren factors of the studied PV systems are presented in Figure 6.
The data presented in Figure 6 show that the average numeric indicator of the f P,nren factor in PV systems is 0.038 (the dotted line). The average numerical value of the f P,nren indicator for PV systems of A group, where the capacity of PV systems is more than 20 (kW), equals 0.021. A similar result was found in B group, where the capacity of PV systems is 20-10 (kW), and the average of f P,nren equals 0.028. However, the average of the f P,nren value for C group (capacity is less than 10 (kW)) reaches 0.065. The data presented in Figure 6 show that the average numeric indicator of the fP,nren factor in PV systems is 0.038 (the dotted line). The average numerical value of the fP,nren indicator for PV systems of A group, where the capacity of PV systems is more than 20 (kW), equals 0.021. A similar result was found in B group, where the capacity of PV systems is 20-10 (kW), and the average of fP,nren equals 0.028. However, the average of the fP,nren value for C group (capacity is less than 10 (kW)) reaches 0.065.
The lowest fP,nren indicator (the value 0.01) was estimated for the 1A, 5A, and 8A power stations, and the highest fP,nren value 0.09 was for the 3C and 5C investigated PV systems. The obtained results lead to the conclusion that the capacity of the PV system may have an impact on the fP,nren indicator.
The frequency distribution (Figure 7) meets a normal Gauss distribution, but it is slightly shifted to the left side. This result shows that the values of fP,nren for C group (capacity is less than 10 (kW)) is bigger 32-43% compared with the A and B groups. The result of the primary energy factor calculation is presented in Table 2. The lowest f P,nren indicator (the value 0.01) was estimated for the 1A, 5A, and 8A power stations, and the highest f P,nren value 0.09 was for the 3C and 5C investigated PV systems. The obtained results lead to the conclusion that the capacity of the PV system may have an impact on the f P,nren indicator.
The frequency distribution ( Figure 7) meets a normal Gauss distribution, but it is slightly shifted to the left side. This result shows that the values of f P,nren for C group (capacity is less than 10 (kW)) is bigger 32-43% compared with the A and B groups. The data presented in Figure 6 show that the average numeric indicator of the fP,nren factor in PV systems is 0.038 (the dotted line). The average numerical value of the fP,nren indicator for PV systems of A group, where the capacity of PV systems is more than 20 (kW), equals 0.021. A similar result was found in B group, where the capacity of PV systems is 20-10 (kW), and the average of fP,nren equals 0.028. However, the average of the fP,nren value for C group (capacity is less than 10 (kW)) reaches 0.065.
The lowest fP,nren indicator (the value 0.01) was estimated for the 1A, 5A, and 8A power stations, and the highest fP,nren value 0.09 was for the 3C and 5C investigated PV systems. The obtained results lead to the conclusion that the capacity of the PV system may have an impact on the fP,nren indicator.
The frequency distribution ( Figure 7) meets a normal Gauss distribution, but it is slightly shifted to the left side. This result shows that the values of fP,nren for C group (capacity is less than 10 (kW)) is bigger 32-43% compared with the A and B groups. The result of the primary energy factor calculation is presented in Table 2. The result of the primary energy factor calculation is presented in Table 2. The investigation has shown that the PEF value of the PV systems of group C with less than 10 (kW) of capacity is bigger (f P,tot = 1.065) than that of the PV systems of the A and B groups with more than 10 (kW) of capacity (respectively, f P,tot = 1.021 and 1.028). This difference consists of 32-43%.

Discussion
The global problem regarding climate change is related to the excessive usage of fossil fuels and greenhouse gas emissions. Taking into account their negative environmental impact as well as the global temperature rise of below 2 • C [36], the EU designates a lot of attention to energy saving. Due to today's energy trends, it is recommended to reduce the energy consumption of the building sector. Currently, about 3% of Europe's building stock meets the top "A+" energy class, called "Nearly Zero Energy Buildings" (NZEBs). This means that the other 97% of buildings are marked by a low energy performance [37]. National activity in the area of building renovation depends on the Member State, while the average data indicate that only 0.4-1.2% of the stock is retrofitted and modernized every year [38]. Renewable energy, for example, solar energy, can be used in NZEBs to provide a reduction in pollutants from the power generation process [39]. As shown in the literature [40], the application of solar water heaters can result in a 10-15% reduction in energy consumption, and solar heating systems are able to make a 45% energy saving in buildings [40]. All in all, renewable energy sources should undoubtedly be used in modern buildings, and they could be an important contributor to the design and development of NZEBs [41].
The PEF is mentioned in several EU Directives for translating the primary utilization into final energy. For instance, EED [19] or EPBD [18]. The EPBD Directive aims to decrease the primary energy demand for buildings. There are possible several improvements in the building envelope, like the isolation of external walls or the replacement of old windows by more airtight ones, leading to savings in final energy. Then, the PEF is applied to convert these savings into primary energy. The EED is the EU-wide Energy Efficiency target that has been expressed in both primary energy and final energy [2]. During the estimation of Member States' savings in primary energy, the PEF is used for a conversion of the final energy savings into primary energy. It means that savings in the area of electricity can be multiplied by 2.5. However, the analysis of the PEF calculation options for electricity in 2016 [25] showed that the PEF incentivizes savings in electricity over direct fossil savings, as necessary to meet the targets of EEP. Nowadays, a significant renewable generation capacity is added to the electricity system, creating a linkage between the electricity and heat sectors. Such actions could successfully help to meet expectations and tackle challenges; however, the current PEF for electricity can be found as an obstacle during the process of decarbonization in the heating sector.
Currently, Member States have a choice to use their own PEF value according to EED and EPBD. The results of the investigation of the review sources are given in Table 3, which shows that EU countries have provided various PEF values. The national laws of many EU member states governing the construction sector do not include or sometimes do not define the PEF values, so the problem with clear values for defining photovoltaic energy appears. It is crucial to develop consistent and normalized values that could be used for both technical and scientific applications, as well as financial analyses. Thus, the PEF is applied to convert the savings in buildings into the value of primary energy, with regard to the energy efficiency of the products and types of fuels used. The PEF is applied to convert electricity (final energy) into primary energy. Such analysis should us allow to compare the energy damages or savings between different countries. Because the energy savings depend on the locations, climate conditions, political aspects, etc., of different countries, directive [22] provides an evaluation of the PEF, including the national factors. Unfortunately, the analysis of researchers' reports shows that not all the EU-28 countries evaluate the PEF. Maybe the reason for that is a negative approach to PEF. For instance, Firlag and Piasecki [52] state that the use of PEFs to estimate the energy consumption in buildings is often incorrect or results in opposite effects. Authors have analyzed systems with dual fuel boilers that are used for burning both biomass and coal that were found to be very common in some Member States-e.g., Poland. The results of their analysis showed that energy characteristics often did not meet reality, as by not using the proper proportion of coal and biomass and their PEFs in calculations, the real energy consumption and emission of air pollutants such as carbon dioxide were significantly greater than the calculated ones. Thus, the estimated ecological targets were not met in reality.
However, in another opinion presented by Esser and Sensfuss [25], PEF is based on a scientific approach and makes the factor well aligned with the actual of power generation. Authors claim that it also reflects well the EU energy mix and the increasing share of renewable sources in electricity production.
So, Lithuania chose to set the PEF for different renewable energy sources. This study focused on the evaluation of the PEFs of PV systems. The values of the total PEF of the PV systems operated in Lithuania, estimated as the average ƒ P,tot = 1.038, are similar to factors used in Germany and France. However, as presented in Table 1, in some of the EU countries the obtained results for PEF differ significantly-e.g., in Ireland or Denmark. The methodology for PEF determination in the countries mentioned is not very clear, and thus it could be related to the capacity of PV systems. This primary outcome highlights the importance of further research.
Overall, this paper proposes guidelines for PEF determination; however, it is necessary to note that each national case should be carefully analyzed. It is recommended to consider the individual NZEB data, the local climate, and regulations, especially when there is a lack of methodology to determinate the PEF values. Undoubtedly, the proper PEF values are crucial for setting a precise primary energy level. Thus, its appropriate estimation influences energy policymaking analysis, energy consumption efficiency in international and national energy statistics, and environmental impact assessments. Thus, all EU members should define the primary energy in PV systems.

Conclusions
The investigation has shown that the total average primary energy factor ƒ P,tot is 1.038. This could depend on the capacity of PV systems. The value of the non-renewable energy factor f P,nren of PV systems which have capacities of less than 10 kW is 0.065; it is 0.028 for PV systems with capacities ranging from 10 to 20 kW, and it is 0.021 for PV systems which have capacities of more than 20 kW. This difference consists of 32-43%. Hence, the total value of the primary energy factor depends on the capacity of PV systems. The study recommends using different values of the non-renewable energy factor f P,nren regarding the capacity of PV systems for the evaluation of the amount of primary energy consumed in NZEB buildings.
In order to achieve the aforementioned goals, all the EU member states must use the same or similar methodology to calculate the primary energy factor of renewable and non-renewable energy sources.
Author Contributions: All the authors contributed equally to this work. All the authors designed the calculations, discussed the results and implications, and commented on the manuscript at all stages. R.T. calculated various primary factor values of non-renewable energy of PV systems of different capacities and interviewed PV system owners/operators. D.A.K. and V.M. analyzed the climate data and prepared the report of electricity transmission system operators. J.Š. wrote the paper. All authors have read and agreed to the published version of the manuscript.