Next Article in Journal
Efficiency Enhancement of Switched Reluctance Generator Employing Optimized Control Associated with Tracking Technique
Next Article in Special Issue
Advanced Techniques for the Modeling and Simulation of Energy Networks
Previous Article in Journal
An Attempt to Apply Canonical Analysis to Investigate the Dependencies between the Level of Organic Farming Development in Poland and the Chosen Environmental Determinants
Previous Article in Special Issue
A Statistical Assessment of Blending Hydrogen into Gas Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling Long-Term Transition from Coal-Reliant to Low-Emission Power Grid and District Heating Systems in Poland

Faculty of Electrical and Control Engineering, Gdańsk University of Technology (Gdańsk Tech), ul. Gabriela Narutowicza 11/12, 80-233 Gdańsk, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(24), 8389; https://doi.org/10.3390/en14248389
Submission received: 13 November 2021 / Revised: 7 December 2021 / Accepted: 10 December 2021 / Published: 13 December 2021

Abstract

:
Energy systems require technological changes towards climate neutrality. In Poland, where the power system is dominated by outdated coal-fired power plants, efforts to minimize the environmental impact are associated with high costs. Therefore, optimal paths for the development of the energy sector should be sought in order to achieve ambitious long-term strategic goals, while minimizing the negative impact on the consumers’ home budget. A methodology and a model for the development of the electricity and heat generation structure were developed and implemented in market allocation (MARKAL) modelling framework. Two scenarios were presented, i.e., business as usual (BAU) and withdrawal from coal (WFC) scenarios. The calculations showed a significant role of nuclear energy and offshore wind power in the pursuit of climate neutrality of electricity generation. In the BAU scenario, the model proposes to stay with coal technologies using carbon capture and storage systems. Withdrawal from coal (WFC scenario) makes it necessary to replace them by gas-fired power plants with CO2 sequestration. Solar energy can be used both in electricity and district heating. In order to build on the latter technological option, appropriate energy storage techniques must be developed. Geothermal energy is expected to be the key option for district heat generation in the long-term horizon. The proposed development paths guarantee a significant reduction in greenhouse gases and industrial emissions. However, complete climate neutrality is uncertain, given the current degree and dynamics of technological development.

1. Introduction

This article is devoted to modelling the development of energy systems in terms of the structure of electricity and heat generation. The subject of the work is in line with the idea of sustainable development and helps to set and achieve the goals of the national energy and climate strategy. The need to reduce greenhouse gases and industrial emissions in order to improve air quality and mitigate climate change was reflected in the legal and strategic documents of the European Union. It is also stimulated by the functioning of the EU Emissions Trading System (EU ETS) [1,2,3,4]. Climate change and emission trade have a broader dimension as the ETS problem was recently addressed in studies, e.g., in China [5,6,7]. These policies force the use of sources that apply low- and zero-emission technologies, including those based on renewable energy resources, often with uncontrolled variability of power generation [8]. The implementation of, e.g., solar technologies is therefore a considerable challenge, as the power system operator must maintain power system stability [9] and tackle voltage issues [10]. In order to be able to balance the power systems at all times and meet energy quality standards, sources of various types of operation must be installed. A rapidly growing market of energy storage systems may constitute a promising solution, while batteries themselves may not guarantee the stable operation of the power system, in which intermittent generation sources will prevail [11]. The growing interest in emission-free modes of transport, in particular electric vehicles, will certainly contribute to a significant increase in demand for electricity. At the same time, the demand for batteries will increase, which may limit their availability for the power industry [12]. However, other solutions can be available, such as the conversion of electricity to hydrogen and its subsequent use, e.g., in fuel cells [13].
The need to use indigenous renewable energy resources and reduce losses in energy supply systems requires the transformation of these systems from a traditional centralized to decentralized structure, planned and coordinated at the local level. However, the total transformation to decentralized systems may not materialize by the middle of this century due to the limitations in technological progress in the field of energy storage and energy management systems, as well as the relatively low energy awareness and activity in the field of energy production in society.
Power systems developed dynamically in the 1960s and 1970s. Initially, they were based mainly on fossil fuels. Currently, many of these systems require modernization, especially in terms of the structure of electricity generation. A similar issue concerns district heating systems. In Poland for instance, heating plants and coal-fired combined heat and power (CHP) plants still dominate. As of the end of 2020, electric installed capacity in Poland was 49.1 GW, of which 46.01% was installed in coal-fired and 17.15% in lignite-fueled utility (public) power plants, the share of gas-based utility power plants was 5.63%, the contribution of industrial autoproducing CHP plants (of which many are still fired with hard coal) was 5.76%, the utility hydro plants share was 4.77%, while wind farms and other renewables share was 20.69% [14]. The domination of coal is even more visible when analyzing annual gross domestic electricity generation. The latter amounted to 152.3 TWh/a in 2020, which was lower than in preceding years, i.e., 158.8 TWh/a in 2019, and 165.2 TWh/a in 2018 [14]. The 2020 shares in total gross domestic annual electricity production were [14]: hard-coal-fired utility power plants—46.97%, lignite-fired utility power plants—24.93%, gas-fueled utility power plants—9.14%, industrial autoproducing CHP plants—6.43%, utility hydro plants—1.77%, wind and other renewables—10.75%. District heat production (not including heat for industrial or commercial processes) totaled 285 PJ/a in 2020, of which public thermal plants (CHP plants and turbines installed in thermal power plants) generated 60%, public heating plants 21%, heat-only boilers in public thermal plants 12%, autoproducing CHP plants 6%, and non-public heating plants 1% [15].
The necessity to change the method of generating and supplying energy to end users creates an opportunity to build an energy system that is environmentally friendly. Past studies [16] indicated that rebuilding a power system based on fossil fuels is also associated with high electricity generation costs. Relying on zero- and low-emission sources is also a financial challenge, but it is worth considering in order to reduce the external costs of energy generation, including the costs of deteriorating health and loss of human life resulting from a polluted environment.
Bearing in mind the above considerations, the authors seek answers to the following questions: (1) How to ensure electricity supply in response to the growing demand in the face of the development of the electric vehicle market? (2) What is an optimal choice of technologies to meet the renewable energy share goals? (3) Should the electricity generation structure completely change from highly centralized to decentralized? (4) What impact will greenhouse gas emission allowance trading schemes have on the choice of electricity and heat generation technologies? (5) When and under what economic conditions will it be justified to build a nuclear power plant? (6) What role should coal play in the transformation of the electricity sector? (7) Is it possible to build a climate-neutral energy system or under what conditions is it feasible?
In order to obtain answers to the above questions, authors developed a methodology that builds on existing modelling frameworks for energy systems optimization, but takes into account modern technologies and market mechanisms. The main area of interest of this work is generation expansion planning of energy systems, where the supply of carriers takes place through power and heating networks, taking into account distributed energy sources and energy microgrids. Their development will be studied in the long term, i.e., until 2060. The modelling methodology is intended to be a tool supporting decisions in the field of investment planning of energy facilities. These investments are expected to lead to the building of sustainable energy systems, so that the existing fossil fuel resources can serve the needs of future generations, not only in terms of energy supply, but also in the production of other goods.
Simultaneous consideration of all long-term energy policy objectives when planning optimal energy systems is a complex process. The databases on energy systems are impossible to handle without computer help. Therefore the use of supporting tools is recommended and the development of models is an important scientific and engineering problem. The most commonly used modelling frameworks are: market allocation (MARKAL), energy flow optimization model (EFOM), model for energy supply system alternatives and their general environmental impacts (MESSAGE), the integrated MARKAL-EFOM system (TIMES), and open source energy modelling system (OSeMOSYS).
The use of energy modelling tools to address generation expansion planning dilemmas was a subject of the previous scientific work. Dagoumas and Koltsaklis [17] discussed generation expansion planning tools with special emphasis on renewable energy sources integration. Gaur et al. [18] went further in their analyses by means of using a unit commitment module to incorporate short-term operational constraints in energy system planning. Yu et al. [19] analyzed air pollutant reduction paths for China using the integrated MARKAL EFOM system (TIMES). Sarica and Tyner [20] used the MARKAL-MACRO model to study the effect of increased US natural gas imports on the energy system. OSeMOSYS was employed to address the problem of balancing high shares of variable renewables through the use of flexible power generation [21] and to build a global energy system model [22]. One of the tool extensions is OSeMOSYS-PuLP [23] which handles large real-world data sets and incorporates Monte Carlo simulations. Most recent analyses on Polish energy system optimal development involve the contribution of Kaszyński and Kamiński [23], who studied demand for hard coal and brown coal in a long-term perspective (by 2050) in view of the current environmental regulations.
Research presented in this paper contributes in many ways to the scientific debate on energy system development, specifically in Polish conditions characterized by a high share of coal-based electricity and heat. Relative to previous work by the authors [16,24,25,26,27,28] and published literature discussed above, this article: (1) presents the method of modelling carbon capture and storage (CCS) technologies; (2) updates datasets used; (3) adds new technological options to the model, e.g., battery energy storage systems; (4) explains the model of EU ETS based on step-wise supply curve for emission allowances; (5) presents new electricity and heat demand projections based on the dynamics of macroeconomic parameters and taking into account dynamic development of the electric car market; and (6) extends the planning horizon to the year 2060.
This article is structured as follows: Section 2 demonstrates the modelling tool used, reference energy system (RES) description, and datasets used. Section 3 summarizes the results, while Section 4 provides discussion of the most important outcomes. Concise conclusions are drawn in Section 5. In Appendix A, detailed results of the analysis are presented, whereas Appendix B describes the MARKAL modelling framework.

2. Materials and Methods

This chapter describes the assumptions and methodology for modelling the development of sustainable energy systems in Poland. The model allows for the optimization of the production structure of energy carriers corresponding to these systems. It includes the production of electricity and network heat from commercial energy sources, i.e., system power plants, combined heat and power plants, and heating companies, supplemented with generation of energy from distributed generation sources, i.e., independent power plants and renewable energy sources, as well as microgeneration.
This model searches for a technology structure that meets the optimization criterion, limited by linear constraints that reflect the functioning of energy systems and the assumptions of national and European energy policy. The main results are the technological structure of net attainable power (capacity), annual net electricity and heat generation, and the amount of emissions. The tool the model is implemented in is MARKet ALlocation (MARKAL) [29], described in Appendix B. The model implemented in this tool allows, inter alia, the inclusion of the EU Emissions Trading System (ETS, EU ETS). The mapping of carbon dioxide capture and storage (CCS) technology is proposed.

2.1. Reference Energy System (RES)

The energy system is understood as a network of interconnections of energy resources in various forms with end users through a set of energy technologies. In the described model, the structure of connections between the flows of individual energy carriers and the corresponding technologies is called the reference energy system (RES). This takes into account data availability constraints and the planning horizon. In the models based on the RES concept, the network infrastructure is presented in a simplified manner. The power system and district heating systems are represented by busbars. They are assigned unit variable operating costs and unit investment expenditure, related to the installed capacity in new production technologies, the construction of which is proposed by the optimization procedure. The model user also determines the efficiency of these systems, taking into account transmission and distribution losses. Connections between power grids, e.g., power stations, are mapped. Electricity imports and exports represent trans-boundary power system connections. The RES diagram for the power infrastructure is shown in Figure 1, and for district heating systems it is presented in Figure 2.

2.2. Modelling Emission Trading Schemes (ETS)

To reflect European Union Emission Trading Scheme (EU ETS) in the MARKAL model, sets of elements representing the type of emissions (in MARKAL marked as ENV—environmental indicator) were created, which distinguish CO2 emissions from installations participating in the EU ETS from those outside the ETS (non-ETS). Then, subsets were created reflecting the emissions from various types of installations participating in the EU ETS, mainly due to the division of installations in the plans of free allocation of emission allowances in the first three trading periods.
The EU ETS group includes technologies for which the typically thermal power supplied in the fuel is higher than 20 MW [2]. At the same time, installations using renewable fuels and those with thermal power in fuel below 20 MW belong to the non-ETS group. The MARKAL model, based on information entered by the user (including emission indicators) and the values of variables representing energy production, calculates the amount of emissions from particular types of technologies. The emissions are constrained in two ways. Firstly, by the number of allowances granted for free (only until 2015 of the model, covering the years 2015–2019), and secondly, by the number of allowances available on the EU ETS market—from the model year 2010 (2010–2014) onwards. The second constraint, however, should additionally take into account the share of allowances available to energy entities covered by the MARKAL model—i.e., power plants, combined heat and power plants, and heating plants in Poland, which in turn would require the inclusion of other EU ETS participants in the model. In the 2010 modelling period, the emission limit was set at the national level, while in the periods 2015–2060—at the level of the entire EU ETS. If the emissions of a given group of installations exceed the limit of free allowances, that group has to buy them at auctions or from other participants in the EU ETS. If the national or European limit of available allowances is exceeded, participants who fail to surrender sufficient number of EUA are fined.
The EU ETS model works as follows. Each electricity- or heat- or both-generating technology in which installations covered by the EU ETS are located, were assigned, apart from the standard CO2 indicator, an emission equivalent (TC1-TC8), depending on whether they belong to the group of producers, described in Figure 3. One EU emission allowance (1 EUA) is required for each emission unit (1 t). Some of these allowances were allocated to power plants for free, but only until the end of the third trading period in the EU ETS (years 2013–2020). The remaining part must be purchased at auctions that are mapped in the model using the IMPT(X)C(Y) purchase option, where X is replaced by a number from 1 to 8 representing a group of producers, and Y is a number from 1 to 5 representing the allowance price level in a given planning subperiod. The sale of allowances is mapped by the export of goods EXP(X)TC1 and only one price level is assumed here. The selling price must be lower than the lowest purchase price of the allowance in order to avoid the problem of finding the optimal value of the objective function. The minimization of the system costs would then lead to the maximization of exports, which, in the absence of a boundary on the quantity of exported goods, would make the optimization problem unsolvable. The number of allowances granted for free to a given group of producers and the number of allowances available at a given price are limited by equations using the ad hoc relationships (ADRATIO) IMPT(1–8)C(1–5) and EXPT(1–8) C1, which are constraints imposed on the variables representing the quantities of commodities being CO2 equivalents. Such a mapping of the auction mechanism is a simplification aimed at taking into account the fact that, apart from Polish power plants, combined heat and power plants, and heating plants, there are other participants of the EU ETS, not only in other sectors of the economy, but also in other countries participating in the system.
Allocations of free allowances, made possible in the first three trading periods of the EU ETS, were taken into account as emission limits; if exceeded, emission allowances must be purchased at auctions. In the first two settlement periods, the sum of the allowances granted was the national emission limit. From the third settlement period, the emission limit was established for the entire EU ETS (ETS cap). As its value after the fourth trading period (years 2021–2030) is not known, it can only be forecast on the basis of the reduction targets proposed in the EU’s energy policy until 2030 [30] and the decarbonization strategy included in the Energy Roadmap to 2050 [3]. It is a separate issue as to what part of the pool of allowances will be available for purchase by Polish entities. Allowance price forecasts have been linked to the number of allowances available in the EU ETS. Five price levels were proposed for each planning horizon and it was determined what share of the entire allowance pool would be available at that price.
Another method, much more labor-intensive, but allowing for more detailed analysis, is the construction of a multi-regional model, in which Poland would be one of them. This task, however, goes beyond the research presented in this article. The reference energy system for the EU ETS is presented in Figure 3.

2.3. Modelling Carbon Capture and Storage (CCS)

In the face of the decarbonization of the energy sector in the EU, CCS technologies are the only chance to preserve coal technologies. They were left in this model to support the decision on whether to leave coal-fired power plants in the power system, e.g., in the case of an absence of other affordable options. The CCS systems have been included in the model by means of a set of technologies and energy carriers (commodities) described in the RES fragment, presented in Figure 4.
The names of energy carriers that allow one to track the sequestration process are written vertically, and the names of emission factors horizontally (The MARKAL model enables processes in which only energy carriers are the input commodities, hence the use of the so-called tracking carriers). Technologies using CCS have been defined in the RES of the MARKAL model, taking into account additional capital expenditure and operating costs as well as lowering the efficiency of electricity generation due to CO2 capture from flue gas or its separation before gasification [31,32,33]. Each of the technologies belonging to the group of power plants equipped with a system to separate carbon dioxide (CCS_ELE) was assigned the CCS emission factor—which is the equivalent of the captured amount of CO2 emissions, the TC2 emission factor—which is the equivalent of the emitted CO2 (assuming separation efficiency of 90% [33,34]), and the carrier energy CCSOUT tracks the amount of CO2 captured from the power plant. Subsequently, the CCSOUT tracking medium is the input of the TRNCCS process, which maps the liquefaction (condensation) and transport of liquid CO2 through pipelines to the place of permanent deposit. The output carrier is CCSSTR which is the equivalent of CO2 delivered to the storage site. The next step is to map the injection using the SALAQU, HYDFLD, and CBMCOA processes, the output carriers of which are CO2 equivalents stored in saline aquifers (CCSAQU), hydrocarbon deposits (CCSHYD), and carbon seams, or by injection with methane recovery from coal beds (CCSCBM). The respective processes EXPCCSAQU1, EXPCCSHYD1, and EXPCCSCBM1 represent the final deposition of CO2 and have a cumulative limit corresponding to the storage potential of CO2 determined on the basis of [35].

2.4. Datasets

The time horizon of the analysis was divided into five-year periods. The introduced assumptions and the results of the solution to the optimization problem, e.g., energy production and emissions, are appropriate for the period of one year, but at the same time are equal in the entire five-year period of time. This is a kind of simplification that is allowed in long-term planning analyses due to the uncertainty of the input data and the need to shorten the calculation time of the model. The first period of the model analysis, called the base year, lasts from 2010 to 2014, for which the data on the condition of the energy sector from 2014 are representative. In the MARKAL model, this period is designated as its first year (2010). This is a period of time that is entirely relevant to the past and for which historical data have been provided. The period 2015 (2015–2019) is the calibrating period—the model computes variable values and proposes new capacity additions, where actual investment has taken place, while capacities and annual generations are constrained by historic data from this period. The planning time horizon ends in 2060 (2060–2064).
The division of the year was made by assigning the full months of the seasons to particular model seasons. Winter covers the months of January–March, spring and autumn (transition season)—April–June and October–December, while summer—July–September. The times of the day (day and night) were selected so that they were equal in all seasons. The day covers the hours from 8.00 am to 7.59 pm, and the night—from 8.00 pm to 7.59 am. The discount rate in the model was set for the entire system at 10%, but energy technologies have individually set discount rates.
Representations of existing power plants and new electricity generation technologies have been drawn up separately. The forecast of changes in the net capacity of the existing power plants and combined heat and power plants was made on the basis of the power transmission system development plan until 2025 [36] and the data contained in the catalogs of public power plants and combined heat and power plants [37] and industrial combined heat and power plants [38]. Data on investments in 2014–2021 were verified on the basis of information on generation resources in the National Power System (NPS) [39]. The operating period of the power unit was determined on the basis of the year of installation and information on the modernization completed. The forecast of changes in electric power is available in the existing utility thermal power plants as well as utility and industrial heat and power plants.
Table 1 presents the technical and economic characteristics of the electricity generation technologies. The forecast of trends in changes in cost indicators was made on the basis of [40]. Discount rates for individual generation technologies were selected on the basis of available studies [41,42].
Figure 5 illustrates the current and forecast fuel prices, including delivery costs, calculated on the basis of the data from [40,43,47].
Figure 6 presents weighted average price of emission allowances. It was obtained on the basis of step-wise supply curves of EUA.
Forecasts of electricity and heat demand at end users were also prepared. Analysts use various approaches to this problem, from the simplest—based on linear increment—to the more complex one, programmed in models specially adapted for this purpose, i.e., MAED. The authors used an approach that is partly based on the methodology of the MAED model. It was aimed at generating electricity and heat demand paths in various sectors of the economy, taking into account different types of sources in sectors, e.g., industrial heat and power plants in the industrial sector, also due to the need to make sectoral demand dependent on various macroeconomic factors or taking into account sectoral trends in changes in demand. The demand or the dynamics of the variability of the demand for final electric energy in the relevant sector of the economy was calculated as a dependent variable of multiple regression. In the case of heat demand, due to difficulties in determining the correlation between demand and economic factors, multiple regression was not applied. The dynamics of changes in demand were based on the extrapolation of trends of changes from the past (2006–2016) in individual subsectors. The independent variables of the regression equation were chosen by the authors of the article. They are quantities having a potential impact on the volatility of demand for final electricity or heat in the economy sector. Forecasts of the volatility of the values of these variables in the considered time horizon were prepared by the authors on the basis of functions describing the volatility trend based on data from the past, taken from the databases of Statistics Poland (GUS) [50,51]. It was assumed that by 2060 the share of electric passenger cars would be 100% and trucks 60%. Electricity demand forecasts are presented in Figure 7, and heat demand forecasts in Figure 8. Despite the measures leading to the efficiency improvement of the final electricity use, the value of the latter increases to 350 TWh/a by 2060, which is mainly the result of the growing market of electric vehicles of all kinds, specifically passenger cars and trucks. For comparison, Energy Policy of Poland to 2040 (PEP2040) [52] assumes that in 2040 the final demand for electricity will amount to 192 TWh/a, and the final demand for district heat will be 214,870 TJ/a.

3. Results

3.1. Electricity Generation and Installed Capacity

3.1.1. BAU (Business as Usual) Variant

BAU variant assumes a dynamic increase in demand caused by the growing share of electric vehicles and slow changes in energy efficiency. Therefore, it is a variant that requires large capacity additions in the electricity generation sector. It also chooses the least-cost technology option, taking into account emission cost and boundaries on capacity additions in offshore wind farms and nuclear reactors. Figure 9 and Figure 10 present the electricity generation and net electric attainable power (capacity) for Poland to 2060 by technology groups. More detailed data are presented in Appendix A. Table A1 contains the historical and planned values of electricity production in generation technologies in Poland until 2060. Table A2 shows the historical and planned values of net attainable electric power in Poland until 2060, and Table A3—the utilization factors of these capacities.
The model analysis showed that in the case of significant electrification of transport, assumed in the BAU variant, annual electricity production should be planned at the level of at least 350 TWh/a in 2060, which is more than twice the value recorded in the period 2015–2019. In addition, the model indicates the need to import an additional 20 TWh/a of electricity in 2060. Considering the above conditions and model assumptions, the total net achievable power in 2060 should be over 103 GW, although it should be emphasized that the peak power of solar photovoltaic and microgeneration plants is 39 GW, with the average annual factor of net attainable electric power utilization being 0.10. The restrictions imposed in the model, reflecting the strategic goals of decarbonizing the energy sector, were implemented in the form of assumptions regarding: (1) high prices of CO2 emission allowances; (2) a limited number of allowances available to Polish generation sources; (3) the indicative target in terms of the share of renewable energy sources in the annual final electricity consumption; (4) tightening environmental standards imposed on power plants and combined heat and power plants.
The implementation of these assumptions naturally forces the search for low- or zero-emission sources. On the other hand, the developed hard coal and lignite mining industry and experience in the use of coal technologies have allowed Poland to satisfy domestic energy needs for many years, maintaining a high level of security of energy fuel supplies in the long term and electricity security in the short term. The energy transformation in this variant is gradual, hence the presence of coal technologies throughout the entire time horizon, taking into account new solutions for the capture and storage of CO2 emissions. Changes in the power system will lead to an increase in the degree of technological diversification, but the whole process is a great challenge in terms of financial expenses for the modernization of the power sector and new solutions in the field of energy storage and management.
Utility thermal power plants built and commissioned by the end of 2014, fired mainly with hard coal and lignite, will be gradually decommissioned until the end of the model time horizon, when they will completely disappear from the electricity generation structure. Constraints imposed on renewable energy sources, in particular wind farms and nuclear power plants, mean that, instead of aging coal-fired power plants, new ones, equipped with CO2 capture and storage (CCS) systems, should be foreseen, even with high assumed emission allowance prices for CO2 and constraints on their number. However, the role of such sources will be limited, and the share in the total national net electricity production will be gradually reduced. At the same time, the amount of captured CO2 in 2060 will amount to 96 million t/a, and in total in the years 2040–2060 it will be 544 million t, which will make use of the assumed CO2 storage potential at the level of 3.5%.
Hard-coal-fired CHP plants, both utility and industrial, will be partially replaced by gas-fired CHP plants. On the other hand, no new investments in biomass combined heat and power plants are expected after 2024, i.e., after the planned investments are completed. The role of cogeneration in electricity production may decrease despite promoting energy efficiency and supporting the technology itself. High electricity needs, with the simultaneous stabilization of the level of heat demand, will result in the fact that energy resources will be used mainly for the production of electricity, and renewable energy resources will be used for the production of heat in heating plants, which are more difficult to use in the power industry in national conditions, e.g., geothermal energy.
According to the model proposal, assuming a unit capacity of 1500 MW, the first block of a nuclear power plant should be commissioned in the period 2030–2034; the next ones are to be built according to the plan implemented by the author, assuming that in one planning period it is possible to build units with a total capacity of up to 3000 MW. Construction of nuclear power plants in the considered planning horizon requires capital expenditure in excess of (2018) EUR 56 billion in the years 2040–2060.
In order to achieve the assumed indicator targets for the share of electricity in renewable sources (50% by 2060), investments in wind and biomass power plants should be provided for, while the production volume in hydro power plants does not change much in the time horizon under consideration. In line with the applicable law [53], the model has limited investments in onshore wind farms. After 2040, the total available capacity of onshore wind farms is close to zero, and for offshore wind farms it would have to increase to 11 GW in the coming years (by 2024), and then to 15 GW in 2060. The total share of renewable sources in net electricity production in 2060 would amount to 31% in this scenario. Solar power plants, despite their high share in the total net attainable electric power (37.6% in 2060), have a 9.8% share in electricity production in 2060, due to the low average annual capacity utilization factor (0.10—Table A3). Biomass power plants play an important role in the production of electricity (6.6% in 2060), as they are treated as a renewable source with zero net CO2 emissions. The high share of sources with unpredictable electric power generation characteristics makes it necessary to maintain intervention sources based on gas turbine technology, with a total net attainable electric power of 1.3 GW in 2060. Their capacity utilization is at zero level in the periods 2020–2034 and 2050–2064, which may mean that they will only constitute a reserve of power in the power system during this period.
Although the average annual factor of utilization of net attainable electric power in the entire power system decreases with time, mainly due to the increase in the share of weather-dependent renewable energy sources, this parameter for basic power plants, i.e., nuclear or coal, lignite and biomass-fired power plants, is most often at a level close to the limit values. The capacity utilization in currently operating coal-fired power plants is decreasing year by year, due to the increasing costs of maintaining these sources, in particular the costs of purchasing CO2 emission allowances. The investment outlays planned for the period 2020–2060 on electricity generation sources amount to almost (2018) EUR 290 billion, of which most will be needed to build offshore wind farms, photovoltaic power plants, and nuclear power plants. Capital expenditures will increase significantly from the planning period 2035–2039, when they will exceed (2018) EUR 25 billion, and in the next five-year periods, these values will further increase, and may even reach (2018) EUR 59 billion.

3.1.2. Withdrawal from Coal (WFC) Variant

The WFC variant assumes withdrawal from coal-based electricity generation by 2050. Consideration of such a move is the result of the Paris Agreement [54] and the conclusions of subsequent Conferences of the Parties (COPs). Figure 11 and Figure 12 present net electric energy production and net attainable electric power planned for Poland by 2060. More detailed data are presented in Appendix A. Table A4 contains the historical and planned values of electricity production in generation technologies in Poland until 2060. Table A5 shows the historical and planned values of net attainable electric power in Poland until 2060, and Table A6—the utilization factors of these capacities.
The calculations show that coal-based electricity generation is mainly replaced by energy from natural-gas-fueled systems.

3.2. Heat Production

3.2.1. BAU Variant

The structure of district heat production for heating purposes and process heat is presented in Table 2. The growing demand for electricity means that energy resources are mainly used for its production, which forces the search for primary energy sources, the use of which in the power industry is less likely than in heating. Such sources may include geothermal energy. However, building a large number of heat sources of this type in the time schedule proposed in the model may prove to be a considerable challenge. In the absence of the appropriate potential, it would be necessary to lean towards solar collectors, but these would require a backup source that would compensate for production losses due to low solar levels or significant storage capacities. An alternative solution would be to obtain heat by converting surplus electricity in the power system. With the shutdown of the existing coal and lignite-fired power plants, the production of district heat from these power plants drops to zero in the considered time horizon. A downward trend can be observed in commercial heating plants, but after 2050 it is reversed due to the aforementioned necessity to use a large amount of energy fuels for the production of electricity, as well as due to the high costs of emission allowances.

3.2.2. WFC Variant

The structure of district heat production for heating purposes and process heat is presented in Table 3. The structures of heat generation are similar in both variants.

3.3. Emissions

Figure 13, Figure 14, Figure 15 and Figure 16 compare emission levels and reduction paths in relation to the base year (2010) that was arbitrarily chosen to show the effectiveness of the modelling tool in designing low-emission energy systems. Detailed results on emissions of sulfur dioxide, nitrogen oxides, and carbon dioxide in the BAU variant are presented in Table A7, Table A8 and Table A9 in Appendix A, respectively. Table A10, Table A11 and Table A12 show the emissions of sulfur dioxide, nitrogen oxides, and carbon dioxide, respectively, in the WFC variant.
There is an observable convergence of emission paths in both the BAU and WFC variants. Emissions in the WFC variant are generally lower than in the BAU scenario, as natural gas replaces coal in power generation.
SO2 emissions from large combustion plants (LCP) mapped in the model amounted to 308 kt/a in 2010, and corresponding NOx emissions in the same year—190 kt/a. CO2 emissions from utility power plants and combined heat and power plants totaled 164 million t/a in 1988, the reference year for the emission reduction targets set out in the Kyoto Protocol.
In relation to the year 2010, the SO2 emissions from electricity and heat generators belonging to the large combustion plants (LCP) group decrease by 2060 by 78% in the BAU variant and 99% in the WFC variant. NOx reduction in the same period was 61% and 91%, respectively, in the BAU and WFC variant.
Emissions outside LCP also decline due to the slowdown in the development of biomass power plants and CHP plants and a small share of sources emitting these compounds, the power of which does not qualify them for the LCP group. It may also be partly related to the functioning of the EU ETS, although this applies only to the group with thermal power supplied in fuel above 20 MW. Thus, the reduction of SO2 and NOX emissions may be a secondary effect of the functioning of the greenhouse gas emission reduction system.
In the period 2010–2060, the reduction of CO2 emissions from energy plants mapped in the model and participating in EU ETS was 89% in the BAU variant and 92% in the WFC variant. The corresponding reduction in utility power plants and CHP plants, belonging to both the ETS and non-ETS sector, was 75% and 79% in the BAU and WFC variant, respectively. If the emission level from the year 1988 (base year for Kyoto Protocol commitments) is used as a reference value, instead of the one from 2010, CO2 emissions drop by 78% in the BAU and 82% in the WFC variant.

4. Discussion

Referring to the research hypotheses posed in the work, after the model analysis, the authors came to the following conclusions:
  • The electrification of passenger and truck transport will cause a dynamic increase in the demand for electricity, and its pace will depend on the scale of projects aimed at improving the efficiency of end-use energy. It is estimated that gross electric energy demand in 2060 may reach values exceeding 350 TWh/a in the worst case (because it requires the largest investments). New baseload power plants will have to meet demanding emission standards, therefore the model indicates the need to build nuclear power plants, and, if possible, those fired with biomass. In the case of power plants based on fossil fuels, investments in CO2 capture systems and infrastructure for its transmission and storage will be necessary to keep them in the fuel mix. The power flexibility requirements will promote solutions based on natural gas combustion.
  • Offshore wind farms will be key investments to meet the benchmark target for the share of electricity from renewable sources. However, in the modelling perspective, it was predicted that their total installed capacity would not exceed 15 GW. Another important option will be photovoltaic systems, although due to the low level of peak power utilization, it will be necessary to build a large number of them, which will allow one to obtain the appropriate level of available power. The costs of their implementation may, however, be higher than envisaged in the model, because in order to maintain the safe operation of the power system, it will be necessary to equip them with energy storage systems. The development of “green hydrogen” technology, although it may result in a further increase in costs, is a promising technological option, especially when one uses fuel cells with trace emissions of industrial pollutants and greenhouse gases.
  • The model promotes solutions in the form of centralized power plants, but if it is necessary to completely withdraw from the combustion of coal and natural gas, it is necessary to look for hybrid system solutions, including solar, wind, and fuel cells. Returning to investing in onshore wind farms will be beneficial from the point of view of the power system if one is able to store and manage electricity from these sources. Small nuclear reactors, which are characterized by lower investment risk and create opportunities for cogeneration, could also contribute to the transformation of the structure. The transition to a completely decentralized generation structure does not seem realistic, however.
  • The CO2 Emissions Trading Scheme is and will be a key tool of the decarbonization policy. The increase in the prices of emission allowances will encourage investors to choose low- and zero-emission technologies. Therefore, in the conditions of a dynamic increase in the price of allowances, the model proposes the construction of nuclear power plants, and in the case of options based on fossil fuels, it is necessary to equip the units with CCS systems. In the district heating sector, technologies based on the use of renewable energy resources, predominantly geothermal energy, should be developed.
  • Taking into account the current pace of technological development, it can be assumed that within 40 years the energy system will be able to function as a set of climate neutral objects. However, looking at the current state of the energy system in general (including transport), the available technological options (including slowly developing fuel cells) and potentially new options (e.g., batteries) are unlikely.
Thus, the model study confirmed the correctness of the research hypotheses. At the same time, the presented vision of development indicates the need to undertake a number of important projects that are also consistent with the assumptions of the Polish energy policy:
  • Commissioning of the first offshore wind farms by 2025;
  • Preparation of industry and human resources for the construction of nuclear power plants by 2025–2030;
  • Preparation of workforce for the operation of nuclear power plants in the 2035–2040 perspective;
  • Conducting research and development works in the field of energy storage systems and energy management in power transmission and distribution systems, in particular in the face of increasing capacity in photovoltaic systems and offshore wind farms;
  • Detailed studies of the domestic storage potential and possibilities for the development of liquefied CO2 transport infrastructure;
  • Detailed research on the potential of geothermal energy to produce district heat for large urban agglomerations.

5. Conclusions

Modelling tools for the development of energy systems are an element of supporting investment decisions in electricity and heat generation technologies. It is true that it is difficult to consider them as an oracle and to rely solely on the results of modelling, but one can set directions for the development of generation infrastructure. With this knowledge, one can plan construction projects, capital expenditures, research and development programs, etc. The elaboration of the model results is only one stage of the entire planning process for the development of energy systems. The next stage includes detailed analyses for individual technological options and projects for the construction of specific power plants. Each project must be preceded by a feasibility study, including a detailed technical and economic analysis, which will help to determine the conditions for the profitability of the project. Each investment is burdened with investment risk and the higher the risk, the higher the investment expenditure on a given project.
Despite covering many of the aspects of generation expansion planning in this study, still there are some issues that require further consideration. Long-term models have limited capabilities to address the intermittency of solar and wind installations, which limits the deployment of energy storages in the modelling solution. It is recommended to include in the model the hybrid systems, consisting of intermittent energy sources and energy storages, which should be preceded by detailed techno-economic analyses aimed at sizing storage capacity and determining cost characteristics of a technology. Furthermore, the storage of large amounts of CO2 may prove to be infeasible over several decades’ worth of perspectives. Inclusion of the options that facilitate utilization of CO2 should be considered in future studies. As the technology is still immature, the techno-economic data, including technology learning rates, are burdened with uncertainty.
Currently, Poland aims at implementing programs that are a part of the national policy document (PEP2040). Construction of the offshore wind farms is at the beginning stage of the investment process, i.e., contracts and agreements are being signed with a perspective of a commissioning by the middle of the present decade. The nuclear program is at the stage of studies of potential localizations, limited to a short list including two sites in Northern Poland. The year 2022 should bring the choice of technology vendor and a final decision regarding the installation of the first nuclear block, which is expected to be completed by the year 2034. The government also commissioned a detailed study on the potential and techno-economic feasibility of carbon capture, utilization, and storage (CCUS) options to examine the possibility of keeping coal-based energy plants in operation. Renewable energy sources are supported by the auctioning system, i.e., specific capacity is commissioned by the president of the Energy Regulatory Office and the winning bid benefits from a premium electricity price (higher than the average wholesale market price) over the 15-year period or until the year 2035. Prosumers (mostly PV system owners) will benefit from net metering until the end of March 2022. After that, excess electricity (not used by the prosumer) will be cleared using the average wholesale market price. The latter move is expected to decelerate the pace of capacity growth in prosumer installation capacity and encourage them to size PV systems to cover mostly their needs. The major concern is the stability of the network and insufficient capacity in electricity storage systems.
Building sustainable energy systems, although difficult, is feasible, and the condition of the current infrastructure, requiring the reconstruction of a significant part of the power plants built in the second half of the twentieth century, is a chance to realize the dream of low-emission energy production. At stake in this game are the life, health, and wellbeing of society, and these will certainly improve when the negative impact of energy on the environment is reduced.

Author Contributions

Conceptualization, M.J. and P.B.; data curation, P.B.; formal analysis, P.B.; investigation, M.J.; methodology, M.J. and P.B.; resources, M.J.; software, M.J.; supervision, P.B.; validation, M.J. and P.B.; visualization, M.J.; writing—original draft, M.J.; writing—review and editing, P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The APC was funded by The Faculty of Electrical and Control Engineering, Gdańsk University of Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Net electricity production in Poland to 2060 by technology (TWh/a)—MARKAL-PL—BAU variant (own study).
Table A1. Net electricity production in Poland to 2060 by technology (TWh/a)—MARKAL-PL—BAU variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)63.444.737.324.26.364.342.540.29--
Thermal PP—lignite, biomass (existing)52.132.230.918.312.95.932.40---
CHP district heating—many fuels (existing)21.117.914.39.929.499.935.780.460.140.12
CHP industrial autoproduction—many fuels (existing)7.294.794.212.492.692.622.161.871.030.48
Hydro PP (existing)2.812.812.942.942.942.943.073.072.472.10
Renewable sources (existing)8.898.898.898.758.897.510.330.330.330.33
Thermal PP—hard coal-29.929.929.929.9-----
Thermal PP—lignite-3.783.783.78------
Thermal PP—hard coal + CCS----20.464.868.570.870.870.8
Thermal PP—lignite + CCS--13.826.239.847.051.053.854.254.6
Nuclear PP---10.921.843.665.476.387.298.2
Gas PP and intervention units-------2.011.751.51
Gas PP with CCS------10.410.410.410.4
CHP district heating—hard coal0.050.190.190.190.190.190.050.100.07-
CHP district heating—natural gas--2.511.111.111.111.111.961.961.96
CHP district heating—biomass0.030.120.120.120.120.150.11---
CHP district heating—waste0.180.390.650.971.291.612.552.873.193.19
CHP industrial autoproduction—natural gas---0.131.151.641.642.302.541.70
CHP industrial autoproduction—biomass-----0.781.562.132.923.24
Wind farms—inland6.616.616.616.616.61-----
Wind farms—offshore-10.810.810.810.810.821.638.648.948.9
Solar farms and PV microgeneration 0.271.381.383.279.2316.025.827.632.134.8
Thermal PP—biomass-0.730.7326.026.026.026.026.024.123.6
Thermal PP—waste1.371.371.371.371.371.37----
Electrical energy storages----------
Total net electricity production164167170188213248292321344356
Import3.9010.516.417.218.119.019.920.922.023.1
Export3.293.293.293.293.293.293.293.293.293.29
Total net electricity supply165174184202228264309338363376
Table A2. Net electric generation capacity of Poland to 2060 by technology (GW)—MARKAL-PL—BAU variant (own study).
Table A2. Net electric generation capacity of Poland to 2060 by technology (GW)—MARKAL-PL—BAU variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)15.213.29.085.422.842.461.770.25--
Thermal PP—lignite, biomass (existing)9.447.575.265.264.502.561.32---
CHP district heating—many fuels (existing)5.965.695.153.923.232.831.770.200.050.05
CHP industrial autoproduction—many fuels (existing)1.771.361.100.990.770.640.540.490.280.16
Hydro PP (existing)2.232.232.232.232.232.232.232.232.232.22
Renewable sources (existing)4.174.174.174.174.173.430.090.090.090.09
Thermal PP—hard coal-3.803.803.803.803.803.803.803.803.80
Thermal PP—lignite-0.480.480.480.480.480.480.480.480.48
Thermal PP—hard coal + CCS----2.598.228.698.988.988.98
Thermal PP—lignite + CCS--1.753.325.055.966.476.836.876.92
Nuclear PP---1.503.006.009.0010.512.013.5
Gas PP and intervention units-------1.271.311.31
Gas PP with CCS------1.431.431.431.43
CHP district heating—hard coal0.030.100.100.100.100.100.100.100.08-
CHP district heating—natural gas0.151.531.531.531.531.531.532.702.702.70
CHP district heating—biomass0.010.050.050.050.050.050.04---
CHP district heating—waste0.050.150.200.200.200.230.400.510.851.00
CHP industrial autoproduction—natural gas---0.181.592.262.262.782.782.78
CHP industrial autoproduction—biomass-----0.260.430.490.720.72
Wind farms—inland3.133.133.133.133.13-----
Wind farms—offshore-3.283.283.283.283.286.5711.714.814.8
Solar farms and PV microgeneration 0.301.541.543.6510.317.928.730.835.838.8
Thermal PP—biomass-0.100.103.573.573.573.573.573.473.47
Thermal PP—waste0.240.240.240.240.240.24----
Electrical energy storages----------
Total net electric capacity42.748.643.247.056.668.081.289.298.7103
Table A3. Utilization factor of the net electric generation capacity of Poland to 2060 by technology (-)—MARKAL-PL—BAU variant (own study).
Table A3. Utilization factor of the net electric generation capacity of Poland to 2060 by technology (-)—MARKAL-PL—BAU variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)0.480.390.470.510.260.200.160.13--
Thermal PP—lignite, biomass (existing)0.630.490.670.400.330.260.21---
CHP district heating—many fuels (existing)0.400.360.320.290.340.400.370.270.300.28
CHP industrial autoproduction—many fuels (existing)0.470.400.440.290.400.470.450.440.420.35
Hydro PP (existing)0.140.140.150.150.150.150.160.160.130.11
Renewable sources (existing)0.240.240.240.240.240.250.420.420.420.42
Thermal PP—hard coal-0.900.900.900.90-----
Thermal PP—lignite-0.900.900.90------
Thermal PP—hard coal + CCS----0.900.900.900.900.900.90
Thermal PP—lignite + CCS--0.900.900.900.900.900.900.900.90
Nuclear PP---0.830.830.830.830.830.830.83
Gas PP and intervention units-------0.180.150.13
Gas PP with CCS------0.830.830.830.83
CHP district heating—hard coal0.220.220.220.220.220.220.050.110.11-
CHP district heating—natural gas--0.190.080.080.080.080.080.080.08
CHP district heating—biomass0.270.270.270.270.270.340.34---
CHP district heating—waste0.400.300.370.550.730.800.730.640.430.36
CHP industrial autoproduction—natural gas---0.080.080.080.080.090.100.07
CHP industrial autoproduction—biomass-----0.340.410.500.460.52
Wind farms—inland0.240.240.240.240.24-----
Wind farms—offshore-0.380.380.380.380.380.380.380.380.38
Solar farms and PV microgeneration 0.100.100.100.100.100.100.100.100.100.10
Thermal PP—biomass-0.830.830.830.830.830.830.830.790.78
Thermal PP—waste0.650.650.650.650.650.65----
Electrical energy storages----------
Average net electric capacity utilization factor0.440.390.450.460.430.420.410.410.400.39
Table A4. Net electricity production in Poland to 2060 by technology (TWh/a)—MARKAL-PL—WFC variant (own study).
Table A4. Net electricity production in Poland to 2060 by technology (TWh/a)—MARKAL-PL—WFC variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)63.444.737.324.212.704.342.540.29--
Thermal PP—lignite, biomass (existing)52.132.230.930.412.95.932.40---
CHP district heating—many fuels (existing)21.117.914.99.809.4910.16.140.470.150.12
CHP industrial autoproduction—many fuels (existing)7.294.794.212.482.602.672.161.871.050.47
Hydro PP (existing)2.812.812.942.942.942.943.073.072.492.32
Renewable sources (existing)8.898.898.898.898.897.510.330.330.330.33
Thermal PP—hard coal-29.929.929.929.929.90.590.690.380.13
Thermal PP—lignite-3.783.783.783.8-----
Nuclear PP---10.921.843.665.476.387.298.2
Gas PP and intervention units--7.47.47.4--1.831.831.83
Gas PP with CCS---1.734.569.7115115115115
CHP district heating—hard coal0.050.190.390.190.190.190.100.100.07-
CHP district heating—natural gas--1.961.111.111.111.111.081.131.13
CHP district heating—biomass0.030.120.240.120.120.240.11---
CHP district heating—waste0.180.390.650.971.291.612.552.873.193.19
CHP industrial autoproduction—natural gas---0.180.991.801.982.092.271.80
CHP industrial autoproduction—biomass-----0.231.091.692.492.62
Wind farms—inland6.616.616.616.616.61-----
Wind farms—offshore-10.810.810.810.810.820.848.859.759.7
Solar farms and PV microgeneration 0.271.383.213.219.3416.826.128.132.735.0
Thermal PP—biomass-0.730.7326.827.127.127.127.125.225.2
Thermal PP—biogas--3.773.773.773.773.77---
Thermal PP—waste1.371.371.371.371.371.37----
Electrical energy storages----------
Total net electricity production164167170188210242283312335347
Import3.8910.516.417.218.119.019.920.922.023.1
Export−3.29−3.29−3.29−3.29−3.29−3.29−3.29−3.29−3.29−3.29
Total net electricity supply165174183201224258299329354367
Table A5. Net electric generation capacity of Poland to 2060 by technology (GW)—MARKAL-PL—WFC variant (own study).
Table A5. Net electric generation capacity of Poland to 2060 by technology (GW)—MARKAL-PL—WFC variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)15.213.29.085.422.842.461.770.25--
Thermal PP—lignite, biomass (existing)9.447.575.265.264.502.561.32---
CHP district heating—many fuels (existing)5.965.695.153.923.232.831.770.200.050.05
CHP industrial autoproduction—many fuels (existing)1.771.361.100.990.770.640.540.490.280.16
Hydro PP (existing)2.232.232.232.232.232.232.232.232.232.22
Renewable sources (existing)4.174.174.174.174.173.430.090.090.090.09
Thermal PP—hard coal-3.803.803.803.803.803.803.803.803.80
Thermal PP—lignite-0.480.480.480.480.480.480.480.480.48
Nuclear PP---1.503.006.009.0010.512.013.5
Gas PP and intervention units--1.021.021.021.021.021.021.021.02
Gas PP with CCS---0.234.759.5915.8415.8415.8415.84
CHP district heating—hard coal0.030.100.100.100.100.100.100.100.08-
CHP district heating—natural gas0.151.531.531.531.531.531.531.491.571.57
CHP district heating—biomass0.010.050.050.050.050.050.04---
CHP district heating—waste0.050.150.200.200.200.230.390.630.891.04
CHP industrial autoproduction—natural gas---0.241.372.492.492.492.492.49
CHP industrial autoproduction—biomass-----0.050.260.420.620.62
Wind farms—inland3.133.133.133.133.13-----
Wind farms—offshore-3.283.283.283.283.286.3114.818.118.1
Solar farms and PV microgeneration 0.301.543.583.5810.418.829.131.336.539.0
Thermal PP—biomass-0.100.103.683.733.733.733.733.633.63
Thermal PP—biogas--0.760.760.760.760.76---
Thermal PP—waste0.240.240.240.240.240.24----
Electrical energy storages----------
Total net electrical capacity42.748.645.345.855.666.282.589.999.7104
Table A6. Utilization factor of the net electric generation capacity of Poland to 2060 by technology (-)—MARKAL-PL—WFC variant (own study).
Table A6. Utilization factor of the net electric generation capacity of Poland to 2060 by technology (-)—MARKAL-PL—WFC variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)0.480.390.470.510.510.200.160.13--
Thermal PP—lignite, biomass (existing)0.630.490.670.660.330.260.21---
CHP district heating—many fuels (existing)0.400.360.330.290.340.410.400.270.320.30
CHP industrial autoproduction—many fuels (existing)0.470.400.440.290.390.480.450.440.420.34
Hydro PP (existing)0.140.140.150.150.150.150.160.160.130.12
Renewable sources (existing)0.240.240.240.240.240.250.420.420.420.42
Thermal PP—hard coal-0.900.900.900.900.900.020.020.010.00
Thermal PP—lignite-0.900.900.900.90-----
Nuclear PP---0.830.830.830.830.830.830.83
Gas PP and intervention units--0.830.830.83--0.200.200.20
Gas PP with CCS---0.830.830.830.830.830.830.83
CHP district heating—hard coal0.220.220.440.220.220.220.110.110.11-
CHP district heating—natural gas--0.150.080.080.080.080.080.080.08
CHP district heating—biomass0.270.270.550.270.270.550.34---
CHP district heating—waste0.400.300.370.550.730.800.750.520.410.35
CHP industrial autoproduction—natural gas---0.080.080.080.090.100.100.08
CHP industrial autoproduction—biomass-----0.550.480.460.460.48
Wind farms—inland0.240.240.240.240.24-----
Wind farms—offshore-0.380.380.380.380.380.380.380.380.38
Solar farms and PV microgeneration 0.100.100.100.100.100.100.100.100.100.10
Thermal PP—biomass-0.830.830.830.830.830.830.830.790.79
Thermal PP—biogas--0.570.570.570.570.57---
Thermal PP—waste0.650.650.650.650.650.65----
Electrical energy storages----------
Table A7. SO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—BAU variant (own study).
Table A7. SO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—BAU variant (own study).
Emission Category2015202020252030203520402045205020552060
Emissions from PP and CHP (non-LCP)22.015.615.38.939.677.737.025.954.013.18
Emissions from DHP (non-LCP)12.210.68.496.015.242.551.670.140.040.04
Emissions from LCP29822820914412512410276.372.069.3
Total emissions from PP, CHP, and DHP33225523315914013411182.476.072.5
Table A8. NOx emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—BAU variant (own study).
Table A8. NOx emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—BAU variant (own study).
Emission Category2015202020252030203520402045205020552060
Emissions from PP and CHP (non-LCP)11.98.527.5214.715.115.614.914.913.412.9
Emissions from DHP (non-LCP)0.000.000.000.000.000.000.000.000.000.00
Emissions from LCP18414714210596.310292.379.877.474.6
Total emissions from PP, CHP, and DHP19615614912011111810794.790.887.5
Table A9. CO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (Mt/a)—MARKAL-PL—BAU variant (own study).
Table A9. CO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (Mt/a)—MARKAL-PL—BAU variant (own study).
Emission Category2015202020252030203520402045205020552060
PP and CHP (non-EU ETS)5.664.864.6221.421.623.023.624.323.623.4
Utility PP existing (EU ETS)10770.663.339.119.09.864.610.230.000.00
Utility PP new (EU ETS)0.0025.726.326.826.09.4910.311.111.110.6
Utility CHP existing (EU ETS)25.221.317.011.711.011.76.580.250.030.02
Utility CHP new (EU ETS)0.070.282.951.471.471.471.252.222.192.08
Industrial autoproduction CHP existing (EU ETS)3.642.402.061.201.231.161.090.940.490.31
Industrial autoproduction CHP new (EU ETS)0.000.000.000.141.231.751.752.452.711.81
Utility DHP existing (EU ETS)6.225.554.613.462.150.750.661.921.992.01
Utility DHP new (EU ETS)0.000.000.000.000.000.000.000.000.000.00
Total emissions (EU ETS)14212611683.962.036.226.319.218.516.8
Total emission over all categories14813112110583.759.249.943.442.140.2
Total emissions from utility PP and CHP13812311410179.155.546.438.136.936.1
Table A10. SO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—WFC variant (own study).
Table A10. SO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—WFC variant (own study).
Emission Category2015202020252030203520402045205020552060
Emissions from PP and CHP (non-LCP)22.015.615.58.989.878.147.195.954.093.12
Emissions from DHP (non-LCP)12.210.68.496.015.242.551.670.240.052.32
Emissions from LCP29822820515610474.736.07.654.643.24
Total emissions from PP, CHP, and DHP33225522917111985.444.913.88.788.68
Table A11. NOx emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—WFC variant (own study).
Table A11. NOx emissions from electricity and heat-generating plants in Poland to 2060 by technology (kt/a)—MARKAL-PL—WFC variant (own study).
Emission Category2015202020252030203520402045205020552060
Emissions from PP and CHP (non-LCP)11.98.5211.2718.319.119.318.514.913.512.8
Emissions from DHP (non-LCP)0.000.000.000.000.000.000.000.000.000.00
Emissions from LCP18414713410474.057.135.019.118.016.9
Total emissions from PP, CHP, and DHP19615614612293.176.353.534.031.429.7
Table A12. CO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (Mt/a)—MARKAL-PL—WFC variant (own study).
Table A12. CO2 emissions from electricity and heat-generating plants in Poland to 2060 by technology (Mt/a)—MARKAL-PL—WFC variant (own study).
Specification2015202020252030203520402045205020552060
PP and CHP (non-EU ETS)5.664.8610.026.927.528.228.724.323.623.3
Utility PP existing (EU ETS)10770.663.352.324.09.864.610.230.000.00
Utility PP new (EU ETS)0.0025.727.627.227.823.04.545.194.984.82
Utility CHP existing (EU ETS)25.221.317.711.711.011.97.030.270.040.03
Utility CHP new (EU ETS)0.070.282.651.471.471.471.321.291.311.21
Industrial autoproduction CHP existing (EU ETS)3.642.402.061.211.231.201.090.940.510.30
Industrial autoproduction CHP new (EU ETS)0.000.000.000.191.051.922.112.232.421.92
Utility DHP existing (EU ETS)6.225.554.613.462.150.750.802.222.352.95
Utility DHP new (EU ETS)0.000.000.000.000.000.000.000.000.000.00
Total emissions (EU ETS)14212611897.468.750.021.512.411.611.2
Total emissions over all categories14813112812496.278.250.236.635.234.6
Total emissions from utility PP and CHP13812312111991.774.446.231.229.929.4

Appendix B

MARKet ALlocation (MARKAL) is a computer program for the development of mathematical optimization models that reflect the current and future technological structure of energy systems. It is aimed at an integrated assessment of the development of the energy sector in one area, e.g., a single country or in many geographically or economically separated and interconnected regions [29]. The model joined the family of tools supporting integrated planning of the development of energy systems in the early 1980s [55,56], as a result of the international cooperation within the ETSAP program financed by the International Energy Agency (IEA).
MARKAL makes it possible to build a model of economic “partial equilibrium” and is characterized by balancing the production and consumption of energy carriers and other goods and fixing their prices. The cost of producing an energy carrier or other goods affects the demand for these goods, and at the same time the demand for an energy carrier or other goods affects their price. Market equilibrium is achieved at a certain price, at which the consumer does not want to buy less than what they need, and no producer wants to produce more than their current capabilities. If the market equilibrium is achieved, either profit maximization or cost minimization is sought, and this approach was used in the MARKAL program [29].
The mathematical structure of the standard MARKAL model takes the form of an optimization problem based on linear programming. As a result of solving the optimization problem, the values of decision variables are obtained, which ensure the minimum cost of the energy system, while meeting the constraints on the set of solutions imposed on it in the form of equations containing the variables and parameters of the model. The key decision variables include: the amount of net attainable power in a given group of technologies, as well as the amount of annual energy production, characteristic for a given technology of production or processing of various energy carriers. The model solution, including the proposed investments in new power plants, combined heat and power plants, and heating plants, is considered optimal for the entire time perspective under consideration.
MARKAL is based on the idea of the reference energy system (RES), which connects primary energy resources with a final or useful energy demand through a network of energy technologies as well as primary, secondary, and final energy carriers. RES is the work of the analyst’s creative invention, and its structure depends on the purpose it sets for the model and the level of detail of the available data. The shape of the RES is, however, subject to certain limitations resulting from the specificity of the tool. Technologies reflected in the MARKAL model are divided into three groups: (1) “PRC—processes”; (2) “conversion technologies” (CON); (3) “demand devices” (DMDs). “Processes” are technologies involving energy transformation that output neither electricity nor heat, which in turn are produced by “conversion technologies”, consisting of subcategories, e.g., power plants, combined heat and power plants, and heating plants. Demand devices convert secondary or final energy into useful energy, i.e., heat for space heating, domestic hot water, district cooling for room air conditioning, etc. On the primary energy side, there are options for the supply or exchange of energy or other goods—Source ENergy Carrier Price (SRCENCP).
The optimization procedure balances the RES so that the demand is met in each considered planning period and in each subperiod of the year (season, time of day). The constraining equations and inequalities make it possible to take into account the availability of power in a defined time period, including forced and planned power outages, the ability of various technologies to meet the peak power demand, the effect of “aging” of the power generation structure, the need to maintain basic power plants, and their share in covering the power demand at night.
The optimization criterion used in the MARKAL model consists in minimizing the net present value of the costs of the energy system operation, reflected in the model. Discounting takes place within a single planning period t, lasting J years, where the cash flows are brought back to the beginning of that period, and then this sum is discounted to the first year of analysis (base year) [29,57]:
z = r = 1 R t = 1 T [ ( 1 + d ) J ( 1 t ) j = 1 j = J ( 1 + d ) j K a n n ( r , t ) ]
where: z—value of the objective function equal to the net present value (NPV) of the energy system costs, mapped in the model [thous. €]; r, R—index and number of regions (geographic areas) mapped in the model; t, T—index and number of planning periods; j, J—index of the year and number of years in one planning period t; Kann(r,t)—annual costs of the operation of the energy system in the region r, in the year included in the planning period t [thous. €/a]; d—general discount rate for the energy system, mapped in the model [1/a].
The method of calculating the annual costs of the energy system operation, used in the MARKAL model, is based on the classical theory of costs in the energy sector. Fuel costs are not included in the variable operating costs because fuels are part of the energy system and their supply chain to energy facilities can be described in terms of energy resources and energy conversion technologies. Therefore, these costs are recorded in items including the extraction, acquisition (e.g., biomass) and import of fuels, energy conversion (technologies), and the supply of fuel to a specific type of technology. In addition, environmental costs are included, which may include environmental charges or emission allowances, or both. Revenues from the sale of energy commodities or other goods outside the considered area (export) reduce the annual costs of the energy system in the MARKAL model, written using the following relationship:
K a n n ( r , t ) = p = 1 P ( K c a p ( r , t , p ) + K f i x o m ( r , t , p ) + K v a r o m ( r , t , p ) + e = 1 E K d e l i v ( r , t , p , e ) ) + s e x p K s u p p ( r , t , s )          s e x p K s u p p ( r , t , s ) + v = 1 V K e n v ( r , t , v )
where: p—technology set (1, 2, …, p); e—energy carrier set (1, 2, …, E); s—supply option set, characterized by three main features: source, energy carrier, and price level (imp—imports, exp—exports, min—mining, rnw—renewables extraction); v—environmental indicators set (pollutants, greenhouse gases, emission equivalents) (1, 2, …, V); Kcap(r,t,p)—capital costs related to energy technology p; Kfixom(r,t,p)—fixed operation and maintenance costs related to energy technology p; Kvarom(r,t,p)—fixed operation and maintenance costs related to energy technology p; Kdeliv(r,t,p,e)—delivery cost of energy carrier e to energy plant built in technology p; Ksupp(r,ts)—costs associated with supply option s; Kenv(r,t,v)—environmental (emission) costs.
The expansion of the equation describing the stream of annual costs of the energy system operation was formulated as follows:
K r ( r , t ) = p = 1 P [ r c a p ( r , p ) k i n v ( r , t , p ) Δ P n ( r , t , p ) + k f i x o m ( r , t , p ) P n ( r , t , p ) + ( k v a r o m ( r , t , p )          + e = 1 E q r , t , p , e k d e l i v ( r , t , p , e ) ) w = 1 W E n ( r , t , p , w )          + s e x p l = 1 L c s u p p ( r , t , s , l ) Q r , t , s , l          s e x p l = 1 L c s u p p ( r , t , s , l ) Q r , t , s , l + v = 1 V k e n v ( r , t , v ) G r , t , v
where:
r c a p r , p = d r , p 1 + d r , p λ r , p 1 + d r , p λ r , p 1  
where: w—set of time slices (subdivisions of the year) representing season and the time of the day (1,2, …, W); l—set of price levels of a commodity (energy carrier) (1, 2, …, L); rcap(r,p)—capital recovery factor [1/a]; kinv(r,t,p)—specific investment cost related to installed capacity [thous. €/MW]; ∆Pn—net capacity addition [MW]; kfixom(r,t,p)—specific fixed operation and maintenance costs [thous. €/(MW ∙ a)]; Pn(r,t,p)—net capacity [MW]; kvarom(r,t,p)—specific variable operation and maintenance costs [thous. €/TJ]; En(r,t,p,w)—annual net energy production [TJ/a]; kdeliv(r,t,p,e)—specific delivery cost of energy carrier e to energy plant built in technology p [thous. €/TJ]; qr,t,p,e –consumption of energy carrier e related to main technology output (electricity—power plants and CHP plants, heat—heat-only plants) [–]; csupp(r,t,s,l)—price level associated with supply option s [thous. €/TJ]; Qr,t,s,l—annual amount of energy (commodity) associated with supply option s; kenv(r,t,v)—cost of emission of pollutant v (or emission allowance) [thous. €/t]; Gr,t,v—annual amount of emitted pollutant v [kt/a]; dr,p—technology-specific discount rate (if not specified, general discount rate of the energy system is applied) [1/a]; λr,p—technical lifetime of a plant built in technology p [a].
The optimal solution, assuming the determination of decision variables at the minimum value of the objective function (Equation (A1)), is limited by typical constraints imposed on these variables:
Δ P r , t , p 0 ,   P r , t , p 0 ,   E r , t , p , w 0 ,   Q r , t , s , l 0 ,   G r , t , v 0
In addition, the set of solutions is limited by linear dependencies describing, inter alia, the balance of energy carriers and goods as well as the power balance of energy facilities. The mathematical structure of the MARKAL model is complex and its discussion can be found in the documentation for this tool [29]. Additionally, the analyst can create their own equations, using the model variables and entering parameters and the value of the left hand side of the equation, in the MARKAL code defined as ADRATIO (ad hoc relationships). It is also possible to modify the source code of the program by creating new variables and equations, if necessary.

References

  1. European Commission. Directive 2010/75/EU of the European Parliament and of the Council of 24 November 2010 on industrial emissions (integrated pollution prevention and control) (Recast). Off. J. Eur. Union 2010, L334, 17–119. [Google Scholar]
  2. European Commission. Directive 2009/29/EC of the European Parliament and of the Council of 23 April 2009 amending Directive 2003/87/EC so as to improve and extend the greenhouse gas emission allowance trading scheme of the Community. Off. J. Eur. Union 2009, L140, 63–87. [Google Scholar]
  3. European Commission. Communication from the Commission: A Roadmap for Moving to a Competitive Low Carbon Economy in 2050 COM(2011) 112 Final; European Commission: Brussels, Belgium, 2011. [Google Scholar]
  4. European Commission. Energy Roadmap 2050 COM(2011) 885 Final of 15 December 2011; European Commission: Brussels, Belgium, 2012. [Google Scholar]
  5. Lin, B.; Jia, Z. The impact of Emission Trading Scheme (ETS) and the choice of coverage industry in ETS: A case study in China. Appl. Energy 2017, 205, 1512–1527. [Google Scholar] [CrossRef]
  6. Lin, B.; Jia, Z. Impact of quota decline scheme of emission trading in China: A dynamic recursive CGE model. Energy 2018, 149, 190–203. [Google Scholar] [CrossRef]
  7. Zhang, L.; Li, Y.; Jia, Z. Impact of carbon allowance allocation on power industry in China’s carbon trading market: Computable general equilibrium based analysis. Appl. Energy 2018, 229, 814–827. [Google Scholar] [CrossRef]
  8. Malkowski, R.; Bucko, P.; Jaskolski, M.; Pawlicki, W.; Stoltmann, A. Simulation of the Dynamics of Renewable Energy Sources with Energy Storage Systems. In Proceedings of the 2018 15th International Conference on the European Energy Market (EEM), Lodz, Poland, 27–29 June 2018; IEEE: Łódź, Poland, 2018; pp. 1–5. [Google Scholar]
  9. Małkowski, R.; Nieznański, J. Underfrequency Load Shedding: An Innovative Algorithm Based on Fuzzy Logic. Energies 2020, 13, 1456. [Google Scholar] [CrossRef] [Green Version]
  10. Szultka, A.; Szultka, S.; Czapp, S.; Lubosny, Z.; Malkowski, R. Integrated Algorithm for Selecting the Location and Control of Energy Storage Units to Improve the Voltage Level in Distribution Grids. Energies 2020, 13, 6720. [Google Scholar] [CrossRef]
  11. Małkowski, R.; Jaskólski, M.; Pawlicki, W. Operation of the Hybrid Photovoltaic-Battery System on the Electricity Market—Simulation, Real-Time Tests and Cost Analysis. Energies 2020, 13, 1402. [Google Scholar] [CrossRef] [Green Version]
  12. Han, X.; Liang, Y.; Ai, Y.; Li, J. Economic Evaluation of a PV Combined Energy Storage Charging Station Based on Cost Estimation of Second-Use Batteries. Energy 2018, 165, 326–339. [Google Scholar] [CrossRef]
  13. Lux, B.; Pfluger, B. A supply curve of electricity-based hydrogen in a decarbonized European energy system in 2050. Appl. Energy 2020, 269, 115011. [Google Scholar] [CrossRef]
  14. PSE SA Polish Transmission System Operator Report on Polish Power System (In Polish: Raport 2020 KSE). Available online: https://www.pse.pl/dane-systemowe/funkcjonowanie-kse/raporty-roczne-z-funkcjonowania-kse-za-rok/raporty-za-rok-2020 (accessed on 1 December 2021).
  15. Walkowska, K.; Berent-Kowalska, G.; Peryt, S.; Dziedzina, K.; Jurgaś, A.; Kacprowska, J.; Gilecki, R.; Parciński, G.; Boczek-Gizińska, R.; Szymańska, M.; et al. Gospodarka Paliwowo-Energetyczna w Latach 2019 i 2020 (Energy Statistics in 2019 and 2020); Szymańska, K., Luty, P., Eds.; Statistics Poland (GUS): Warsaw, Poland, 2020. [Google Scholar]
  16. Jaskólski, M. Modelling long-term technological transition of Polish power system using MARKAL: Emission trade impact. Energy Policy 2016, 97, 365–377. [Google Scholar] [CrossRef]
  17. Dagoumas, A.S.; Koltsaklis, N.E. Review of models for integrating renewable energy in the generation expansion planning. Appl. Energy 2019, 242, 1573–1587. [Google Scholar] [CrossRef]
  18. Gaur, A.S.; Das, P.; Jain, A.; Bhakar, R.; Mathur, J. Long-term energy system planning considering short-term operational constraints. Energy Strateg. Rev. 2019, 26, 100383. [Google Scholar] [CrossRef]
  19. Yu, Y.; Jin, Z.; Li, J.; Wu, Y.; Jia, L. An Air Pollutant Emission Reduction Path of China’s Power Industry. Atmosphere 2020, 11, 852. [Google Scholar] [CrossRef]
  20. Sarica, K.; Tyner, W.E. Economic impacts of increased U.S. exports of natural gas: An energy system perspective. Energies 2016, 9, 401. [Google Scholar] [CrossRef] [Green Version]
  21. Gardumi, F.; Welsch, M.; Howells, M.; Colombo, E. Representation of Balancing Options for Variable Renewables in Long-Term Energy System Models: An Application to OSeMOSYS. Energies 2019, 12, 2366. [Google Scholar] [CrossRef] [Green Version]
  22. Löffler, K.; Hainsch, K.; Burandt, T.; Oei, P.-Y.; Kemfert, C.; von Hirschhausen, C. Designing a Model for the Global Energy System—GENeSYS-MOD: An Application of the Open-Source Energy Modeling System (OSeMOSYS). Energies 2017, 10, 1468. [Google Scholar] [CrossRef] [Green Version]
  23. Dreier, D.; Howells, M. OSeMOSYS-PuLP: A Stochastic Modeling Framework for Long-Term Energy Systems Modeling. Energies 2019, 12, 1382. [Google Scholar] [CrossRef] [Green Version]
  24. Jaskólski, M. Application of MARKAL model to optimisation of electricity generation structure in Poland in the long-term time horizon Part I–concept of the model. Acta Energ. 2012, 12, 15–25. [Google Scholar]
  25. Jaskólski, M. Application of MARKAL model to optimisation of electricity generation structure in Poland in the long term time horizon. Part II. Model and forecasts assumptions. Acta Energ. 2012, 13, 4–13. [Google Scholar]
  26. Jaskólski, M.; Bućko, P. MARKAL long-term power generation scenarios for Poland–Concept of the model. Rynek Energii 2013, 107, 108–114. [Google Scholar]
  27. Jaskólski, M. Optimization of Nuclear Power Share in the Structure of Electricity Production in Poland in Time Perspective by 2060. Acta Energ. 2014, 20, 55–60. [Google Scholar] [CrossRef]
  28. Jaskólski, M.; Bućko, P. MARKAL long-term power generation scenarios for Poland: Increasing the share of renewable energy sources by 2040. Rynek Energii 2015, 117, 117–125. [Google Scholar]
  29. Loulou, R.; Goldstein, G.; Noble, K. Documentation for the MARKAL Family of Models. Available online: http://www.iea-etsap.org/web/MrklDoc-I_StdMARKAL.pdf (accessed on 31 August 2015).
  30. European Commission. A Policy Framework for Climate and Energy in the Period from 2020 to 2030 COM(2014) 15 Final; European Commission: Brussels, Belgium, 2014. [Google Scholar]
  31. Kannan, R.; Strachan, N.; Pye, S.; Anandarajah, G.; Balta-Ozkan, N. UK MARKAL Model Documentation. Available online: https://www.ucl.ac.uk/energy-models/models/uk-markal/uk-markal-documentation (accessed on 12 October 2015).
  32. International Energy Agency IEA-Investment Costs. Available online: http://www.worldenergyoutlook.org/weomodel/investmentcosts/ (accessed on 4 April 2016).
  33. Gerbelová, H.; Versteeg, P.; Ioakimidis, C.S.; Ferrão, P. The effect of retrofitting Portuguese fossil fuel power plants with CCS. Appl. Energy 2013, 101, 280–287. [Google Scholar] [CrossRef]
  34. Wang, X.; Du, L. Study on carbon capture and storage (CCS) investment decision-making based on real options for China’s coal-fired power plants. J. Clean Prod. 2016, 112, 4123–4131. [Google Scholar] [CrossRef]
  35. Wójcicki, A. CO2 Geological Storage Potential in Poland. Available online: https://www.iea.org/media/workshops/2009/ccswarsaw/wójcicki.pdf (accessed on 17 March 2016).
  36. Polskie Sieci Elektroenergetyczne, S.A. Plan Rozwoju w Zakresie Zaspokojenia Obecnego i Przyszłego Zapotrzebowania na Energię Elektryczną na Lata 2010–2025 (Wyciąg). Available online: http://www.pse.pl/uploads/kontener/Plan_Rozwoju_2010_2025.pdf (accessed on 31 August 2015).
  37. Agencja Rynku Energii. Katalog Elektrowni i Elektrociepłowni Zawodowych. Stan na 31.12.2014; Biczyński -Bicz, W., Wrońska, I., Eds.; Agencja Rynku Energii: Warszawa, Poland, 2015. [Google Scholar]
  38. Agencja Rynku Energii. Katalog Elektrociepłowni Przemysłowych. Stan na 31.12.2014; Biczyński-Bicz, W., Wrońska, I., Eds.; Agencja Rynku Energii: Warszawa, Poland, 2015. [Google Scholar]
  39. PSE SA Informacje o Zasobach Wytwórczych KSE (wg Stanu na 30.11.2020). Available online: https://www.pse.pl/documents/20182/89a41169-b3e1-4ff4-8bff-ce89019297cd?safeargs=646f776e6c6f61643d74727565 (accessed on 1 March 2021).
  40. Capros, P.; De Vita, A.; Tasios, N.; Papadopoulos, D.; Siskos, P.; Apostolaki, E.; Zampara, M.; Paroussos, L.; Fragiadakis, K.; Kouvartakis, N.; et al. EU Energy, Transport. and GHG Emissions. Trends to 2050. Reference Scenario 2013; Publications Office of the European Union: Luxembourg, 2014. [Google Scholar]
  41. Oxera. Discount Rates for Low-Carbon and Renewable Generation Technologies; Oxera Consulting Ltd.: Oxford, UK; Brussels, Belgium, 2011. [Google Scholar]
  42. García-Gusano, D.; Espegren, K.; Lind, A.; Kirkengen, M. The role of the discount rates in energy systems optimisation models. Renew. Sustain. Energy Rev. 2016, 59, 56–72. [Google Scholar] [CrossRef]
  43. International Energy Agency; Nuclear Energy Agency. Projected Costs of Generating Electricity 2010; OECD Publishing: Paris, France, 2010. [Google Scholar]
  44. Agencja Rynku Energii Budowane i Planowane Elektrownie. Available online: http://www.rynek-energii-elektrycznej.cire.pl/st,33,335,tr,145,0,0,0,0,0,budowane-i-planowane-elektrownie.html (accessed on 6 March 2018).
  45. Zakeri, B.; Syri, S. Electrical energy storage systems: A comparative life cycle cost analysis. Renew. Sustain. Energy Rev. 2015, 42, 569–596. [Google Scholar] [CrossRef]
  46. Ahmad, S.; Tahar, R.M. Selection of renewable energy sources for sustainable development of electricity generation system using analytic hierarchy process: A case of Malaysia. Renew. Energy 2014, 63, 458–466. [Google Scholar] [CrossRef]
  47. European Commission. Commission Staff Working Paper. Impact Assessment. Energy Roadmap 2050; COM(2011) 885 final {SEC(2011) 1566 final; SEC(2011) 1569 final}; European Commission: Brussels, Belgium, 2011. [Google Scholar]
  48. Ministerstwo Gospodarki. Projekt Polityki Energetycznej Polski do 2050 Roku; Warszawa, Poland, 2015. Available online: https://www.civitas.edu.pl/wp-content/uploads/2015/03/Securitologia_2_2016_105_118.pdf (accessed on 6 December 2021).
  49. DECC. Updated Short-Term Traded Carbon Values Used for UK Public Policy Appraisal; Department of Energy & Climate Change: London, UK, 2015. [Google Scholar]
  50. Główny Urząd Statystyczny Bank Danych Lokalnych. Available online: https://bdl.stat.gov.pl/BDL/start (accessed on 30 March 2018).
  51. Główny Urząd Statystyczny Bank Danych Makroekonomicznych. Available online: https://bdm.stat.gov.pl/ (accessed on 30 March 2018).
  52. Ministerstwo Klimatu i Środowiska. Polityka Energetyczna Polski do 2040 r. Załącznik do Uchwały nr 22/2021 z Dnia 2 Lutego 2021 r; Ministerstwo Klimatu i Środowiska: Warszawa, Poland, 2021. [Google Scholar]
  53. Sejm Rzeczypospolitej Polskiej. Ustawa z Dnia 20 Maja 2016 r. o Inwestycjach w Zakresie Elektrowni Wiatrowych; Dz. U. z 2016 r. poz. 961; Kancelaria Sejmu RP: Warszawa, Poland, 2016. [Google Scholar]
  54. United Nations. Paris Agreement; United Nations, 2015. Available online: https://unfccc.int/sites/default/files/english_paris_agreement.pdf (accessed on 6 December 2021).
  55. Egberts, G. MARKAL-Ein LP-Modell für die Internationale Energieagentur. OR Spektrum 1981, 3, 95–100. [Google Scholar] [CrossRef]
  56. Fishbone, L.G.; Giesen, G.; Goldstein, G.A.; Hymmen, H.A.; Stocks, K.J.; Vos, H.; Wilde, D.; Zölcher, R.; Balzer, C.; Abilock, H. User’s Guide for MARKAL (BNL/KFA Version 2.0); Brookhaven National Laboratory: Upton, NY, USA, 1983. [Google Scholar]
  57. Jaskólski, M. Modelowanie Rozwoju Regionalnych Systemów Energetycznych ze Szczególnym Uwzględnieniem Bioenergii; Politechnika Gdańska, Wydział Elektrotechniki i Automatyki: Gdańsk, Poland, 2006. [Google Scholar]
Figure 1. Reference energy system representation of the power grid (own study). Note: demand categories (each category has at least one corresponding demand technology, e.g., IE1)—final electricity in respective economy sectors (and subcategories): IE—industry, AE—agriculture, CE—commercial, RE—residential, TE—transport (excluding electric cars), TA—transport (electric cars); energy carriers: HCO—hard coal, LIG—lignite, NGA—natural gas, OTH—other fuels, BIO—biomass, BGS—biogas, BBG—biomass and biogas (existing plants), HYD—hydro energy, WIN—wind energy, LQD—liquid, URN—nuclear fuel, SOL—solar energy, GEO—geothermal energy, SMW—solid municipal waste; electrical grids: ELC—transmission public, ELI—industrial autoproduction (separated from public), EMG—microgrids, EHV—distribution public; conversion technologies: power substations: PTE—highest voltages (400 kV, 220 kV)/high voltages (110 kV) substations, PIN/PNI—substations’ connection of industrial autoproduction plants with public grid (in both directions); generator technology categories: E10—public coal- and lignite-fired thermal plants, E20—public nuclear plants, E30—independent renewable, E40—public gas plants, E50—public hydro and pumped-storage plants, E60—public and independent combined heat and power (CHP) plants, E80—industrial autoproduction CHP plants, storages: E9A, E9B, E9C—battery storages; resources/transboundary flows: IMPELC1, IMPELC2, EXPELC1—imports and exports of electricity, MINHCO1, IMPHCO1—domestic mining and imports of hard coal, MINLIG1—lignite mining, MINNGA1, IMPNGA1—domestic extraction and imports of natural gas, MINOTH1—extraction of other fuels (“dummy” category for existing plants), RNWBIO1, RNWBGS1, RNWBBG1—extraction of biomass and biogas (the latter is a “dummy” category), IMPURN1—nuclear fuel imports, RNWSMW1—domestic acquisition of solid municipal waste.
Figure 1. Reference energy system representation of the power grid (own study). Note: demand categories (each category has at least one corresponding demand technology, e.g., IE1)—final electricity in respective economy sectors (and subcategories): IE—industry, AE—agriculture, CE—commercial, RE—residential, TE—transport (excluding electric cars), TA—transport (electric cars); energy carriers: HCO—hard coal, LIG—lignite, NGA—natural gas, OTH—other fuels, BIO—biomass, BGS—biogas, BBG—biomass and biogas (existing plants), HYD—hydro energy, WIN—wind energy, LQD—liquid, URN—nuclear fuel, SOL—solar energy, GEO—geothermal energy, SMW—solid municipal waste; electrical grids: ELC—transmission public, ELI—industrial autoproduction (separated from public), EMG—microgrids, EHV—distribution public; conversion technologies: power substations: PTE—highest voltages (400 kV, 220 kV)/high voltages (110 kV) substations, PIN/PNI—substations’ connection of industrial autoproduction plants with public grid (in both directions); generator technology categories: E10—public coal- and lignite-fired thermal plants, E20—public nuclear plants, E30—independent renewable, E40—public gas plants, E50—public hydro and pumped-storage plants, E60—public and independent combined heat and power (CHP) plants, E80—industrial autoproduction CHP plants, storages: E9A, E9B, E9C—battery storages; resources/transboundary flows: IMPELC1, IMPELC2, EXPELC1—imports and exports of electricity, MINHCO1, IMPHCO1—domestic mining and imports of hard coal, MINLIG1—lignite mining, MINNGA1, IMPNGA1—domestic extraction and imports of natural gas, MINOTH1—extraction of other fuels (“dummy” category for existing plants), RNWBIO1, RNWBGS1, RNWBBG1—extraction of biomass and biogas (the latter is a “dummy” category), IMPURN1—nuclear fuel imports, RNWSMW1—domestic acquisition of solid municipal waste.
Energies 14 08389 g001
Figure 2. Diagram of the reference energy system for district heating systems (own elaboration). Note: demand category (each having at least one corresponding demand technology denoted, e.g., IH1): low-temperature heat demands in respective economy sectors, i.e., IH—industrial, AH—agriculture, CH—commercial, RH—residential; process heat demand in economy sector i.e., IP—industry; heating plants: H10—public heating plants, H20—private/independent heating plants, H60—heating boilers in power plants and combined heat and power plants, public, independent, and owned by heating companies, H80—heating boilers in industrial autoproduction CHP plants; remaining denotations as in Figure 1.
Figure 2. Diagram of the reference energy system for district heating systems (own elaboration). Note: demand category (each having at least one corresponding demand technology denoted, e.g., IH1): low-temperature heat demands in respective economy sectors, i.e., IH—industrial, AH—agriculture, CH—commercial, RH—residential; process heat demand in economy sector i.e., IP—industry; heating plants: H10—public heating plants, H20—private/independent heating plants, H60—heating boilers in power plants and combined heat and power plants, public, independent, and owned by heating companies, H80—heating boilers in industrial autoproduction CHP plants; remaining denotations as in Figure 1.
Energies 14 08389 g002
Figure 3. Model structure for European Union Emission Trading Scheme (EU ETS) (own elaboration). Note: TC1-TC8—emission equivalents (environmental indicators): TC1—utility (public) PP—existing, TC2—utility (public) PP—new, TC3—utility (district heating) CHP—existing, TC4—utility (district heating) CHP—new, TC5—industrial autoproduction CHP—existing, TC6—industrial autoproduction CHP—new, TC7—utility DHP—existing, TC8—utility DHP—new; T1C-T8C—corresponding EUA equivalents (energy carriers); IMPT(X)C(Y)—purchase option of EUA, EXP(X)TC1—sales option of EUA.
Figure 3. Model structure for European Union Emission Trading Scheme (EU ETS) (own elaboration). Note: TC1-TC8—emission equivalents (environmental indicators): TC1—utility (public) PP—existing, TC2—utility (public) PP—new, TC3—utility (district heating) CHP—existing, TC4—utility (district heating) CHP—new, TC5—industrial autoproduction CHP—existing, TC6—industrial autoproduction CHP—new, TC7—utility DHP—existing, TC8—utility DHP—new; T1C-T8C—corresponding EUA equivalents (energy carriers); IMPT(X)C(Y)—purchase option of EUA, EXP(X)TC1—sales option of EUA.
Energies 14 08389 g003
Figure 4. Model structure for carbon capture and storage (CCS) (own study). Note: E1E-E1H—coal-based technologies with CCS, E3H—biomass-based technology with CCS, E4C—gas-fueled technology with CCS, CCSOUT—energy carrier tracking CO2 captured in power plants, TC2—CO2 emissions from power plants covered by the EU ETS, CCSTRN—CO2 tracking medium transported to the permanent deposition site, CCSAQU, CCSHYD, CCSCBM—CO2 tracking media stored in saline aquifers (AQU), hydrocarbon fields (HYD) and coal seams or by injection with coal bed methane recovery (CBM), TRNCCS—CO2 transport process from a power plant to a deposit site, SALAQU, HYDFLD, CBMCOA—processes mapping CO2 storage, EXPCCSAQU1, EXPCCSHYD1, EXPCCSCBM1—technologies denoting the permanent deposition of CO2 in appropriate deposits.
Figure 4. Model structure for carbon capture and storage (CCS) (own study). Note: E1E-E1H—coal-based technologies with CCS, E3H—biomass-based technology with CCS, E4C—gas-fueled technology with CCS, CCSOUT—energy carrier tracking CO2 captured in power plants, TC2—CO2 emissions from power plants covered by the EU ETS, CCSTRN—CO2 tracking medium transported to the permanent deposition site, CCSAQU, CCSHYD, CCSCBM—CO2 tracking media stored in saline aquifers (AQU), hydrocarbon fields (HYD) and coal seams or by injection with coal bed methane recovery (CBM), TRNCCS—CO2 transport process from a power plant to a deposit site, SALAQU, HYDFLD, CBMCOA—processes mapping CO2 storage, EXPCCSAQU1, EXPCCSHYD1, EXPCCSCBM1—technologies denoting the permanent deposition of CO2 in appropriate deposits.
Energies 14 08389 g004
Figure 5. Fuel price forecasts until 2060, own study based on the data from [40,43,47].
Figure 5. Fuel price forecasts until 2060, own study based on the data from [40,43,47].
Energies 14 08389 g005
Figure 6. Weighted average of forecast EUA prices (own study based on the data from [40,48,49]).
Figure 6. Weighted average of forecast EUA prices (own study based on the data from [40,48,49]).
Energies 14 08389 g006
Figure 7. Forecasts of the final electricity demand (own study based on the data from [50,51]). Note: IE—industry, AE—agriculture, TE—transport (excl. electric cars), RE—residential, CE—commercial, TA—transport (electric cars).
Figure 7. Forecasts of the final electricity demand (own study based on the data from [50,51]). Note: IE—industry, AE—agriculture, TE—transport (excl. electric cars), RE—residential, CE—commercial, TA—transport (electric cars).
Energies 14 08389 g007
Figure 8. Forecasts of the final demand for district heat (own study based on the data from [50,51]). Note: IH—industrial, AH—agriculture, CH—commercial, RH- residential.
Figure 8. Forecasts of the final demand for district heat (own study based on the data from [50,51]). Note: IH—industrial, AH—agriculture, CH—commercial, RH- residential.
Energies 14 08389 g008
Figure 9. Net electrical energy production planned for Poland by 2060—BAU variant (own study).
Figure 9. Net electrical energy production planned for Poland by 2060—BAU variant (own study).
Energies 14 08389 g009
Figure 10. Net attainable electric power planned for Poland by 2060—BAU variant (own study).
Figure 10. Net attainable electric power planned for Poland by 2060—BAU variant (own study).
Energies 14 08389 g010
Figure 11. Net electrical energy production planned for Poland by 2060—WFC variant (own study).
Figure 11. Net electrical energy production planned for Poland by 2060—WFC variant (own study).
Energies 14 08389 g011
Figure 12. Net electric power (capacity) planned for Poland by 2060—WFC variant (own study).
Figure 12. Net electric power (capacity) planned for Poland by 2060—WFC variant (own study).
Energies 14 08389 g012
Figure 13. SO2 emission levels and reduction relative to 2010—power plants, combined heat and power plants, and heating plants falling into the category large combustion plants (LCP).
Figure 13. SO2 emission levels and reduction relative to 2010—power plants, combined heat and power plants, and heating plants falling into the category large combustion plants (LCP).
Energies 14 08389 g013
Figure 14. NOX emission levels and reduction relative to 2010—power plants, combined heat and power plants, and heating plants falling into the category large combustion plants (LCP).
Figure 14. NOX emission levels and reduction relative to 2010—power plants, combined heat and power plants, and heating plants falling into the category large combustion plants (LCP).
Energies 14 08389 g014
Figure 15. CO2 emission levels and reduction relative to 2010—power plants, combined heat and power plants, and heating plants participating in European Union Emission Trading Scheme (EU ETS).
Figure 15. CO2 emission levels and reduction relative to 2010—power plants, combined heat and power plants, and heating plants participating in European Union Emission Trading Scheme (EU ETS).
Energies 14 08389 g015
Figure 16. CO2 emission levels and reduction relative to 1988 and 2010—utility power plants and combined heat and power plants.
Figure 16. CO2 emission levels and reduction relative to 1988 and 2010—utility power plants and combined heat and power plants.
Energies 14 08389 g016
Table 1. Technical and economic characteristics for new electricity generation technologies used in the MARKAL model. Own study based on the data from: [31,32,40,41,42,43,44,45,46]. Note: r, t, p—sets: region, time period, and technology, respectively; tstart—first model year technology becomes available; λr,p—technical lifetime; ns—availability factor; ηe(r,t,p)—net electrical efficiency; kn(r,t,p)—specific investment cost; kes(r,t,p)—specific operation and maintenance (O&M) cost; us—capacity utilization at peak power system load; ϭP—power to heat ratio; dr,p—technology-specific discount rate; eCO2(r,t,p)emission of CO2 per unit of generated electricity; PP—power plant, CHP—combined heat and power plant; SC—supercritical; PCC—pulverized coal combustion; IGCC—integrated gasification combined cycle; CCS—carbon capture and storage; FBC—fluidized bed combustion; CCGT—combined cycle gas turbine, PWR—pressurized water reactor; ORC—organic Rankine cycle; NA—not applicable.
Table 1. Technical and economic characteristics for new electricity generation technologies used in the MARKAL model. Own study based on the data from: [31,32,40,41,42,43,44,45,46]. Note: r, t, p—sets: region, time period, and technology, respectively; tstart—first model year technology becomes available; λr,p—technical lifetime; ns—availability factor; ηe(r,t,p)—net electrical efficiency; kn(r,t,p)—specific investment cost; kes(r,t,p)—specific operation and maintenance (O&M) cost; us—capacity utilization at peak power system load; ϭP—power to heat ratio; dr,p—technology-specific discount rate; eCO2(r,t,p)emission of CO2 per unit of generated electricity; PP—power plant, CHP—combined heat and power plant; SC—supercritical; PCC—pulverized coal combustion; IGCC—integrated gasification combined cycle; CCS—carbon capture and storage; FBC—fluidized bed combustion; CCGT—combined cycle gas turbine, PWR—pressurized water reactor; ORC—organic Rankine cycle; NA—not applicable.
Technology Name and Symboltstartλr,pnsηe(r,t,p)kn(r,t,p)kes(r,t,p)usϭPdr,peCO2(r,t,p)
-a%%€/kW
/kW/a
---kg/GJ
201020602010206020102060
Thermal PP—hard coal (SC PCC/FBC)E1A20155090%43%53%1589158946.81.00NA0.10216175
Thermal PP—lignite (SC PCC/FBC)E1B20105090%43%53%2221222146.81.00NA0.10260211
Thermal PP—hard coal (SC IGCC CCS)E1D20253590%32%42%3780268199.11.00NA0.177658
Thermal PP—hard coal (SC PCC CCS)E1E20255090%32%42%3600255384.21.00NA0.172922
Thermal PP—lignite (SC PCC CCS)E1F20255090%32%42%3291285784.21.00NA0.172621
Thermal PP—lignite (SC FBC CCS)E1G20255090%32%42%5687493784.21.00NA0.172621
Thermal PP—lignite (SC IGCC CCS)E1H20253590%32%42%5867509399.11.00NA0.172621
Nuclear PWR generation III/III+E2A20306083%33%38%45003949117.01.00NA0.1300
Wind farms—inlandE3A201025*100%100%1300115033.50.35NA0.1000
Wind farms—offshoreE3B202025**100%100%4500282990.50.47NA0.1400
Solar farms—PVE3C201525***100%100%200078817.90.00NA0.0900
Microgeneration –PVE3D201520***100%100%3240259862.41.00NA0.1300
Thermal PP—biomass (IGCC GTCC)E3E20203583%58%58%3240311826.51.00NA0.1300
Thermal PP—biomass (IGCC CCS)E3F20303583%34%34%3888311899.11.00NA0.1300
Thermal PP—biomass (IGCC)E3G20253583%44%44%3240259835.11.00NA0.1300
Distributed generation—gas engine—biogasE3H20152557%30%30%2340234085.01.00NA0.1000
PP—CCGT—municipal wasteE3J20153065%50%50%66304630241.81.00NA0.1000
Intervention units—gas turbinesE4A201535100%40%40%39039015.61.00NA0.09138138
Thermal PP—natural gas (GTCC)E4B20154083%55%61%89877819.51.00NA0.0910090
Thermal PP—natural gas (GTCC CCS)E4C20253583%53%53%2200181135.11.00NA0.171010
Distributed generation—fuel cells—natural gasE4E20202550%40%22%4680195087.41.00NA0.15143143
CHP district heating—hard coalE6A20154044%23%40%2317231733.50.500.400.10407407
CHP district heating—natural gasE6B20153567%27%27%1014101430.40.500.510.09296296
CHP district heating—biomass ORCE6C20153055%22%22%31512894118.60.500.400.1300
CHP district heating—biogasE6D20152546%35%35%7742625588.90.500.900.1000
Small-scale CHP—fuel cells—natural gasE6E20202090%49%49%4000372887.41.002.460.15116116
Microcogeneration—gas microturbine (p < 120 kW)E6F20152570%33%33%4000311819.50.500.700.09242242
CHP district heating—municipal wasteE6G20152580%23%23%999960005001.000.450.0600
Industrial CHP—hard coalE8A20154044%23%23%2317231733.50.500.400.10407407
Industrial CHP—natural gasE8B20154067%27%27%1014101430.40.500.510.10296296
Industrial CHP—biomassE8C20154055%22%22%31512894118.60.500.400.1000
Industrial CHP—biogasE8D20154046%35%35%7742588588.90.500.900.1000
Energy storage (battery)—transmission gridE9A20201525%70%70%1650165016.51.00NA0.1000
Energy storage (battery)—distribution gridE9B20201025%75%75%2360236016.11.00NA0.1000
Energy storage (battery)—microgridsE9C20201025%80%80%2750275011.21.00NA0.1000
* Capacity utilization variability for inland wind farms: spring/autumn-day: 0.24; spring/autumn-night: 0.22; summer-day: 0.18; summer-night: 0.15; winter-day: 0.34; winter-night: 0.33; own computation based on [46]. ** Capacity utilization variability for offshore wind farms: spring/autumn-day: 0.38; spring/autumn-night: 0.37; summer-day: 0.31; summer-night: 0.30; winter-day: 0.45; winter-night: 0.45; own computation based on [46]. *** Capacity utilization variability for solar PV systems: spring/autumn-day: 0.18; spring/autumn-night: 0.03; summer-day: 0.25; summer-night: 0.04; winter-day: 0.12; winter-night: 0.01; own computation based on the data from [46].
Table 2. Net heat production from utility and industrial autoproduction plants in Poland to 2060 by technology (PJ/a)—MARKAL-PL—BAU variant (own study).
Table 2. Net heat production from utility and industrial autoproduction plants in Poland to 2060 by technology (PJ/a)—MARKAL-PL—BAU variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)36.425.721.513.93.662.501.460.17--
Thermal PP—lignite, biomass (existing)16.710.49.935.894.141.900.77---
CHP district heating—many fuels (existing)47039731621821822514214.44.383.64
CHP industrial autoproduction—many fuels (existing)13186.576.445.348.747.438.933.918.88.82
CHP district heating—hard coal0.431.731.731.731.731.730.430.860.64-
CHP district heating—natural gas--17.77.87.87.87.813.813.813.8
CHP district heating—biomass0.271.081.081.081.081.351.01---
CHP district heating—waste1.403.115.197.7410.312.920.423.025.525.5
CHP industrial autoproduction—natural gas---0.918.1411.611.616.217.912.0
CHP industrial autoproduction—biomass-----7.0014.019.226.229.1
District heating plants (existing)61.655.045.634.221.37.424.26---
District heating plants—geothermal/solar------84.2184201222
Heating plants—industry and other sectors------3.9132.133.133.5
Heat boilers in CHP district heating14.310.77.113.0711.38.416.080.640.040.03
Heat boilers in CHP industrial autoproduction1.130.850.560.220.880.790.610.510.430.16
Total net heat production733592503340337336338339342348
Table 3. Net heat production from utility and industrial autoproduction plants in Poland to 2060 by technology (PJ/a)—MARKAL-PL—WFC variant (own study).
Table 3. Net heat production from utility and industrial autoproduction plants in Poland to 2060 by technology (PJ/a)—MARKAL-PL—WFC variant (own study).
Technology Group Name2015202020252030203520402045205020552060
Thermal PP—hard coal, natural gas, biomass (existing)36.425.721.513.97.302.501.460.17--
Thermal PP—lignite, biomass (existing)16.710.49.939.784.141.900.77---
CHP district heating—many fuels (existing)47039733121421822814914.74.553.80
CHP industrial autoproduction—many fuels (existing)13186.576.445.047.048.338.933.919.28.56
CHP district heating—hard coal0.431.733.471.731.731.730.860.860.64-
CHP district heating—natural gas--13.87.87.87.87.87.68.08.0
CHP district heating—biomass0.271.082.171.081.082.171.01---
CHP district heating—waste1.403.115.197.7410.312.920.423.025.525.5
CHP industrial autoproduction—natural gas---1.256.9912.714.014.816.012.7
CHP industrial autoproduction—biomass-----2.069.815.222.523.6
District heating plants (existing)61.655.045.634.221.37.424.26--14.1
District heating plants—geothermal/solar------77.0191208227
Heating plants—industry and other sectors------6.1137.039.225.4
Heat boilers in CHP district heating14.310.77.113.0711.38.416.081.280.080.13
Heat boilers in CHP industrial autoproduction1.130.850.560.220.880.790.610.510.430.39
Total net heat production733592517340338337338340344350
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jaskólski, M.; Bućko, P. Modelling Long-Term Transition from Coal-Reliant to Low-Emission Power Grid and District Heating Systems in Poland. Energies 2021, 14, 8389. https://doi.org/10.3390/en14248389

AMA Style

Jaskólski M, Bućko P. Modelling Long-Term Transition from Coal-Reliant to Low-Emission Power Grid and District Heating Systems in Poland. Energies. 2021; 14(24):8389. https://doi.org/10.3390/en14248389

Chicago/Turabian Style

Jaskólski, Marcin, and Paweł Bućko. 2021. "Modelling Long-Term Transition from Coal-Reliant to Low-Emission Power Grid and District Heating Systems in Poland" Energies 14, no. 24: 8389. https://doi.org/10.3390/en14248389

APA Style

Jaskólski, M., & Bućko, P. (2021). Modelling Long-Term Transition from Coal-Reliant to Low-Emission Power Grid and District Heating Systems in Poland. Energies, 14(24), 8389. https://doi.org/10.3390/en14248389

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop