1. Introduction
Current global developments once again demonstrate the need for reliable, readily available and predictable access to (electrical) energy. Solar energy has the potential to be a large quantitative source of long-term, reliable energy that is not dependent on any particular market player [
1]. In order to provide electrical power according to the demand of the grid, the fluctuating availability of local solar power must be balanced. Concentrating Solar Power (CSP) provides relatively inexpensive medium-term energy storage in the form of heat that can be converted into electricity on demand. There are line-focusing and point-focusing CSP systems.
Point source CSP plants concentrate sunlight from a large area via mirrors (heliostats) onto a small receiver aperture at the top of a tower, which can easily reach temperatures above 1000 °C [
2]. Due to the potential to reach high temperatures, solar tower plants (STPs) also offer good opportunities for chemical processing to produce substances that act as storage and/or fuel for other processes (ammonia, H2, etc.) [
3]. Line-focusing CSP plants use parabolic mirrors to concentrate sunlight onto receiver tubes through which a heat transfer fluid flows. Parabolic trough plants (PTPs) typically operate in the 300–400 °C temperature range [
4].
In solar towers and parabolic troughs, the heat generated can be stored in molten salt storage tanks and converted into electricity by a power unit. The power unit consists of a heat engine using a thermodynamic cycle and a generator. The design of a CSP plant is mainly determined by the size of the solar field and receiver, the storage capacity, and the gross electrical rated output of the power unit.
The operating costs of generating electricity and meeting demand generally increase with the size of the components. At the same time, the capital cost per unit of electricity decreases. In addition, the estimation of these values is highly dependent on the daily and monthly variations in solar energy as well as the energy demand schedule. The balance between demand and generation is achieved by energy storage. Therefore, finding the economically optimal design for a given site and demand profile becomes a non-trivial task that can be solved using optimization methods. In this article, a new analytical optimization approach is proposed to investigate the fraction of available solar energy that should be used by CSP at a given site, so that the rest can be covered by other systems e.g., photovoltaic (PV).
The method presented in this research can also be used to calculate the optimal ratio of heat production from CSP to electricity production from PV in solar fuel production systems, according to the scheme presented in [
5].
The focus of this study is the projected concentrated solar power plant as a part of hybrid renewable power plants. The subject of this study is the functional dependence between the levelized cost of electricity (LCOE) and the projected CSP power. LCOE is an indicator used to compare the cost of generating electricity from different sources and is the main criterion for the economic viability of CSP construction [
6]. It includes capital, operating, maintenance, fuel, financing and decommissioning costs over the life of the plant.
The purpose of the article is to build a target function to determine the optimal integrated CSP capacity to obtain the minimum levelized cost of electricity. Two dimensionless parameters were used for this purpose: the design factor (DF) and solar multiple (SM). The SM is a factor used in CSP installations to compare the actual thermal output of the solar field with the thermal output required for the power unit. It helps to optimize the size of the solar field, which allows for efficient energy storage and maximizes electricity production [
7]. An SM greater than 1 indicates an excess capacity for storage or increased generation, while an SM equal to 1 means that the solar field produces exactly the same thermal power as the power unit [
8]. The design factor is the coefficient equal to the ratio of the CSP capacity to the maximum required capacity of the power plant.
2. Materials and Methods
The issue of power optimization in CSP plants has been raised by different authors to determine the features of a hybrid CSP plant [
9], to develop a strategy to mitigate the effects of high-power electricity shortages [
10], to investigate the impact of weather forecast uncertainties [
11] and to improve the efficiency of CSP cavity absorbers [
12]. Optimizing the design of a mirror field is a compromise between the possible thermal power of the field and the cost [
13]. One of the most effective procedures for determining the characteristics of the mirror field is to determine the radial shift [
14]. In this paper, it is proposed to consider the optimization of the solar power plant capacity embedded in a hybrid system to determine the main parameters of influence on which to focus and to find new engineering solutions to improve CSP plants.
The main criteria for the implementation of solar power plants is the determination and comparison of the following values: the levelized cost of electricity (LCOE), net present value of the project (NPV), benefit-to-cost ratio (BCR), total capital cost, and annual load satisfaction factor (ALSF) [
15]. According to the authors, the most influential of the above criteria that affects the use of CSP and the construction of new plants is the price per kWh of electricity produced: the LCOE. At a given power consumption and due to the nature of the change in solar radiation during the year (given by the schedule of solar irradiance for the duration of the year), it is necessary to determine the characteristic parameters of the CSP (solar receiver, tower height, field size and layout, power block size, thermal energy storage (TES) size) considering the minimization of the price per kWh of electricity produced.
The available solar power of the solar power plant
and its geographically advantageous location are determined on the basis of GIS (geographical information system) assessment using the average annual amounts of direct normal irradiation (DNI) [
16]. As a rule, this task is solved with the help of state-of-the-art tools to optimize the LCOE of CSP plants like SAM (NREL) or Greenius from DLR [
17,
18].
SAM is a techno-economic computer model designed to facilitate decision making for people involved in the renewable energy industry [
19]. For example, a model of a SUPCON Delingha 50 MW MS ST plant was developed with SAM and a simulation took place in order to calculate the annual electrical output performance of the commercial power plant [
20]. SAM is a very useful tool, making it possible for solar energy professionals to analyze photovoltaic systems and concentrate solar power parabolic trough systems in the same modeling platform using consistent financial assumptions. Such a PV-CSP power plant with a parabolic trough of 50 MW gross capacity has been simulated with SAM for two different locations in Southern Spain and South Africa [
21]. It has also been used for a research study which investigated the concept of providing both heat and power from a PV and CSP hybrid plant to meet the energy demand of LNG export terminals [
22].
Another approach is the use of a Hybrid Optimization and Performance Platform (HOPP), as mentioned in [
23], in order to evaluate the technological and financial performance of a CSP-PV hybrid system without detailed modeling of annual operations. A HOPP is an open-source modeling tool that uses a Python-based scripting interface to access and combine underlying single-technology performance models in NREL’s System Advisor Model (SAM) to evaluate the performance and financial viability of hybrid renewable energy systems [
24].
Greenius [
25] is a simulation tool for the annual yield assessment of solar energy systems. It can be used to calculate the operation of a solar power plant depending on irradiance data in hourly or sub-hourly fidelity using a set of technical parameters to define the plant design as well as a predefined operation strategy. It is usually used to evaluate a typical operational year (TOY) of a solar power plant with operational results describing the plant state in every time step of a typical meteorological year (TMY). With the results of the TOY and additional economic parameters, defining the costs and financing a detailed economic evaluation can be performed, calculating the cash flow and key performance indicators like the levelized cost of electricity (LCOE) and the capacity factor (CF). The modeling in Greenius is based on the energy flow between the components, assuming steady-state behavior. Some transient effects in the scope of the time step size are accounted for using simplified equations. Such effects occur for components with large thermal inertia, e.g., the heating up and cooling down process of the solar field or the power block. To optimize a plant layout for a certain site, Greenius can be used to evaluate different plant designs via grid search with varying technical parameters and compare the results using the key performance indicators.
In contrast to a numerical solution, here, the authors propose an analytical approach to optimization. The analytical approach is simple and fast to solve and allows for the minimization of the investment to obtain a profit close to the maximum. It also shows the capacity limit below which the construction of a CSP plant is not cost-effective. The analytical approach to the equation also allows the designer of a CSP plant to quickly perform a comprehensive analysis and find the conditions under which the LCOE at a given site will change its character. This will allow for the analysis of how changes in these conditions will change the optimum of the LCOE function. It is planned that the proposed analytical approach, in combination with existing LCOE calculation software, will allow for a comprehensive analysis of the optimization of several variables in the design process of CSP plants.
The proposed analytical approach is centered around the incorporation of the design factor and the solar multiple factor into a mathematical optimization model. This model specifically targets the power capacity of concentrated solar power plants with thermal storage as the primary energy source.
Hybrid CSP systems can combine the best of CSP and PV renewable energy. By seamlessly integrating solar thermal and electrical conversion methods, these systems ensure a consistent and efficient power supply. Such hybrid system solutions, including connection with a PV field, have been identified as a viable solution to reduce the LCOE of CSP plants while maintaining the flexibility and high capacity factors granted by the TES unit [
26]. The PV unit is used to directly produce cheap power during daylight, while the CSP unit offers the storage capacity, e.g., of molten salt TES, as shown in [
27], which guarantees the required flexibility and power production during evening and night, increasing the plant capacity factor. With the added benefit of energy storage integration, hybrid CSP systems promise uninterrupted power even during periods of low sunlight or high demand.
In optimization models, the input information is the load forecast data and the projected electricity demand. The amount of electricity to be generated is defined as
where
—the time period considered, usually
;
—the necessary electric power to fulfill the demand schedule from CSP;
—the number of hours of operation per year, according to which the amount of electricity produced will be equal to the yearly energy consumption. It is usually referred to as the full-load hours (FLH), or the annual full-load hours for a year.
Figure 1 shows the main considerations regarding the dependence of the required CSP power demand on the CSP operating time. The area of the graph (yellow) corresponds to the CSP power demand. If we represent the selected area (the same amount of electricity generated) as a rectangle with a height of
, we obtain the imaginary number of hours h max that the CSP should operate to obtain
. If a CSP is designed with less than the required capacity for a given area
, the lack of electricity can be compensated for by energy storage and PV. Thus, the red line shows the theoretical operation of a CSP plant with a capacity of
and possibly three periods of plant shutdown.
—the normalized design capacity of the CSP.
—the maximum required electric output power of the hybrid plant from demand schedule. For CSP without PV,
.
The full-load hours are calculated as follows:
The CSP electricity generation per year is calculated as follows:
The target function of the complex optimization of energy supply is the sum of the specific costs of electricity, €/(MWh). For this approach, the concept of the levelized costs of electricity (LCOE) to define the economically viable specific costs of electricity production in €/MWh is used [
28]:
where
TLCC—total life cycle cost;
Ey—produced energy in year y;
Y—total number of years considered for the operation of the plant;
i—considered discount rate.
The results of annual simulations for PV CSP hybrid concepts are mainly compared, using the levelized electricity cost as the figure of merit for comparing the different configurations, as shown in [
29].
Using the assumption of constant yearly operational costs defined as operational expenditures (OPEX) and the capital invested in the beginning defined as capital expenditures (CAPEX) as well as a constant yearly electricity production of
Eyear, the annuity factor to simplify the equation is introduced as
The capital recovery factor is as follows:
A weighted average cost of capital method [
29] for the discount rate is used to depict the costs of equity
and debt
less the income taxes:
To account for further tax costs, the fixed charge rate FCR, according to [
30], is introduced:
where
—the present value of the depreciation cash flows.
The depreciable money benefits can be calculated depending on the present value of the depreciated cash flows assuming a linear depreciation schedule:
To represent the optimal yearly electricity production in consideration of the production costs and the demand–coverage ratio, the LCOE is extended by further specific cost components in €/MWh representing the specific reduced cost due to the damage caused by the network failure flow
and costs due to penalties
:
According to the life cycle assessment performed in reference [
31], the operation phase decreases with increasing capital expenditures in CSP. So
. But there is still no specific model for building an optimization function for changing the cost of electricity from the power of the plant.
3. The Proposed Function of Optimization
3.1. The Mathematical Model
For CSP with thermal storage, it is necessary to design a solar field larger than the required thermal load for the power unit. Therefore, the solar multiple factor (SM) [
32,
33] is introduced as the ratio of the available thermal power from the solar field to the input thermal power of the power unit:
where
—the thermal power output of the solar field at the reference DNI value;
—the design nominal electric output power CSP according to its efficiency.
For the averaged values,
where
and
—the effectiveness CSP for
and
.
—efficiency CSP. In this study, it is assumed to be linearly dependent on the receiver efficiency. Then, according to [
34],
=
. Here,
a,
b—const
a > 0,
b < 1. Hence,
does not depend on the CSP capacity.
The amount of electricity to be produced by the CSP plant is given by
where
—the energy from thermal storage. So
To find the analytical optimization function, the concept of the design factor by analogy with a combined heat and power plant is introduced as in [
35]. Here,
DF is a coefficient equal to the ratio of the CSP capacity to the maximum required capacity of the power plant.
is defined as
where
—the design nominal electric output power of the CSP according to its efficiency.
—the maximum required electric output power of the hybrid plant from the demand schedule.
If
, the maximum electric output power of the CSP plant is required from the demand schedule (the design power of the entire power plant), so
where
—the energy from storage, MWh/day;
day—the number of days per year of CSP working time.
The storage capacity relative to the nominal thermal input power of the power block in h/day is calculated using
where
—the efficiency of the power block;
—the thermal energy storage capacity in MWh/day.
This equation shows the relationship between the required energy and the designed capacity of the CSP output unit through the parameters and the efficiency of the plant. Assume that the dependence of the CSP efficiency on the receiver efficiency is known. Then, is dependent on three parameters, SM and DF.
To reduce the number of input parameters, we can define the location where the CSP is designed. The electrical power of the CSP depends on the efficiency, the DNI value, and the mirror area.
The area of the solar field for range
is
where
—the hours of DNI data resolution.
Assuming a constant efficiency for the normalized available operating hours, the following equation results:
where
—the area of the solar field for reaching
.
To consider the minimum and maximum capacity of the power generation system during deviations from the selected average value,
Assuming that excess power is generated due to deviations between the maximum power generation and design and is fed into the energy storage, the following equation results:
And Equation (18) is used:
Or suppose that the energy for the energy storage is extracted by increasing the area of the mirror field; then, the following equation can be used:
Given that the dependence is determined by the type and internal characteristics of the designed CSP, it is assumed that is known.
In reality, not all of the solar energy can be used because the storage has a maximum capacity and the components need to warm up and cool down. Depending on the storage size, the irradiance distribution, the sunshine hours, and the solar multiplier, the real normalized operating hours can be defined as
where
—the reference DNI for the normalized value for operating hours;
—the sunshine duration, h/year;
—the relative load corresponding to the average arrival of solar radiation [
36].
where
denote the maximum and minimum distance among historical days to the current forecasting day, respectively.
β is a free parameter controlling the effective number of profiles which are averaged.
Dj is the distance between the current day and the historical day.
The accuracy of load forecasting for CSP is directly affected by input uncertainties such as the weather forecast. These are specified by black box models, or white box models, to provide a probabilistic load forecast [
37]. Using Gaussian kernel density estimation, the procedure converts the point load forecast into a probabilistic load forecast based on historical data provided by the internal and external monitoring system.
3.2. The Function of Optimization
The target function is the following:
where
—mathematical expectation of the specific impaired value due to losses caused by the flow of failures in the network;
—the cost and the penalty function.
The
CAPEX can be divided into overnight capital costs
OCC and costs for network connections
as well as costs for financing the construction
.
The EPC costs are divided into component costs, land costs including preparation, and service costs and are extended by a contingency factor
. The
OCC cost components all depend on the size of the plant and therefore on
. The scaling can be performed using an exponential economy of scale relationship, as described in [
38].
where
Ccomp—capital investment for all components and equipment [€];
Cservice—capital investment for all engineering, procurement and construction services [€];
—capital investment to purchase and prepare the land [€];
—design factor (Equation (14));
fcont—contingencies relative to the investment costs [%];
scomp—exponential scaling factor for component costs,
scomp = [0.7, 0.9] [
38];
sservice—exponential scaling factor for service costs,
sservice = [0.3, 0.5] [
38].
where
fservice—factor of service costs to component costs, which are usually around 5% [
38].
The transit costs are mainly associated with the construction of transit networks from the CSP to the distribution point:
where
—capital investments in electricity transit networks, €/m;
n—the number of parallel strands of wire 1 or 2;
—the network performance indicator;
—the transit network length, m;
—the minimum cost of the necessary equipment, €/MWh.
The construction finance costs depend on the capital investment and therefore can be introduced relative to the
OCC and transit costs. They include costs for up-front and commitment fees
, for the interest during construction
, assumed to be accountable for half of the construction period
, and additional money loaned for reserve accounts
:
The annual operational expenditures
OPEX are divided into fixed and variable costs
and
, both in €, and are introduced as [
39]
Fixed operating costs are the costs of land lease, insurance, personnel, and the operation and maintenance of components. They can be represented by a functional relationship to the electrical capacity of the plant
:
where
—the relative fixed OPEX costs per year in
€/kW-yr.Variable operational costs are costs in € depending on the yearly electricity produced, e.g., for auxiliary electric consumption or water use. This leads to
where
—the relative variable OPEX costs per year in
€/MWh-yr.The mathematical expectation of the specific reduced cost due to the damage caused by the internal network failure flow is given by the following equation:
where
—the cost of repairing a mirror unit or PV panel €/(pcs);
—the number of mirrors or PV panels;
—the mirror failure rate, 1/(year);
—the number of hours of PV operation per year;
—the single mirror area.
Or
where
—the power of one panel, €/MWt.
Using this equation, it should be noted that represents the ratio of the used plant size to the maximum plant size not only for the solar field but for all components and corresponding processes, influencing the costs as well as the yearly energy produced.
The function of costs and penalties is proposed and investigated by the authors in [
40] for a 50 MW CSP plant with molten-salt-based
, given the theoretical support from the government for the introduction of a green tariff at which
so
can be negative.
The target function has been calculated as
This function has three variables, and DF.
4. Results and Discussion
The following results are calculated according to the proposed mathematical dependencies with different technical and economic indicators of the plant. A techno-economic analysis has been developed for 15–150 MW CSP which used the function of optimization in this research. The storage capacity was from 0 to 10 h.
The type of CSP is that from a solar tower. The solar multiple is 1.5. The grid connection costs and costs of transit are not included. The location is Almeria, Spain. The efficiency of the CSP is a function of the power capacity and the efficiency ratio of the TES is set at 0.9. The prices of the equipment and service were chosen as those for 2022.
Figure 2 shows the validation of the analytical model with the above parameters and the numerical model calculated by Greenius.
As can be seen from
Figure 2, for the selected parameters, the data convergence is very high, which allows us to recommend the selected analytical approach for building an optimization model of hybrid power plants.
The same SPP parameters as before (without TES) and the same extremes of the function are obtained when the price of Battery Energy Storage Systems (BESSs) for PV storage is added (
Figure 3). From the analysis of
Figure 2 and
Figure 3, it can be seen that with the selected economic indicators, the optimization function has a minimum at the point
. So, with lower PV prices, and the higher the share of PV in the plant, the better. Therefore, the methods for finding the roots of the optimization function are not presented in this article. But this is true only for selected technical and economic indicators and may vary depending on the economic and political conditions of the area.
Figure 4 shows the percentage of electricity generated by different components of the hybrid power plant, calculated using the analytical model. In this case, it is assumed that all the excess energy from the SM goes to the energy storage, including additional energy losses. This figure shows that due to the TES, the generation of electricity from the PV and CSP is not linear.
Since, in the proposed analytical method, the storage size function is set through the solar multiple factor,
Figure 5 shows the change in the LCOE with the design factor DF. As the power of the concentrated station and the size of the mirror field increase, the size of the drive also increases, but the LCOE does not increase linearly. The nonlinearity is caused by the difference in the cost of the CP TES and the mirror field.
Despite the fact that the developed optimization model under the current economic conditions has an optimum at DF → 0, the very method of building this model allows us to analyse various indicators of the designed hybrid power plant.
Figure 6 also shows the existence of an SM optimum, which shows the importance of calculating the optimal capacity of a hybrid solar station for a specific geopolitical location in order to obtain the minimum LCOE.
5. Conclusions
This paper presents a novel analytical approach for optimizing the capacity of concentrated solar power plants in hybrid configurations. By incorporating the design factor and the solar multiple factor into a mathematical optimization model, the method provides a more detailed analysis of the design possibilities while offering faster computational speeds. This advancement contributes to ongoing efforts to improve the efficiency and feasibility of concentrated solar power plants, particularly in hybrid systems.
The developed analytical approach to the performance optimization of CSP plants has shown that it is possible to use analytical functions to obtain the same results as those usually obtained by numerical methods.
The cost of electricity in a hybrid power plant has been analyzed. Hybrid power plants show a relative LCOE with increasing storage capacity, but the cost of the PV-BESS system increases significantly due to higher storage costs. The costs of heat production from CSP without storage and PV systems with an electric heater are comparable.
As the design factor increases, the capital expenditures increase, but the operating costs decrease. Further research in the field of analytical calculations of the levelized costs of electricity for hybrid power plants could involve the calculation of the threshold level of economic indicators for which the use of CSP will be economically more feasible and the function LCOE = f(DF), DF ∈ (0;1) will be extreme.
The dependence LCOE = f (SM) is not linear and has a clear optimum (for the selected economic indicators at SM = 1.3; see
Figure 5).
From the above, there are two different ways to determine the optimal storage capacity. The first way is to take the optimal SM value as a constant of the ratio of the mirror field area to the CSP capacity according to the CSP efficiency change graph, and direct the energy that is higher than this to the storage during peak hours. The second way is to take as constant the value of the of the nominal electrical output power of the CSP according to its efficiency and to increase the SM to obtain energy for the energy storage. According to the authors, the best way to commence this research is to optimize these two functions using the regulator simplex method.