Next Article in Journal
Design of a Novel Hybrid Concentrated Photovoltaic–Thermal System Equipped with Energy Storages, Optimized for Use in Residential Contexts
Previous Article in Journal
New Approaches in Finite Control Set Model Predictive Control for Grid-Connected Photovoltaic Inverters: State of the Art
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Approach to Power Optimization of Concentrating Solar Power Plants with Thermal Storage

1
German Aerospace Center (DLR), Institute of Solar Research, Langenbroich 13, 52428 Jülich, Germany
2
Solar-Institut Jülich, University of Applied Sciences FH Aachen, Heinrich-Mußmann-Str. 5, 52428 Jülich, Germany
3
Department of Fundamental and Applied Mathtematics, Zaporizhzhia National University, Zhukovsky St. 66, 69600 Zaporizhzhia, Ukraine
*
Author to whom correspondence should be addressed.
Solar 2024, 4(3), 509-525; https://doi.org/10.3390/solar4030024
Submission received: 7 May 2024 / Revised: 16 September 2024 / Accepted: 19 September 2024 / Published: 21 September 2024

Abstract

:
This paper deals with the problem of determining the optimal capacity of concentrated solar power (CSP) plants, especially in the context of hybrid solar power plants. This work presents an innovative analytical approach to optimizing the capacity of concentrated solar plants. The proposed method is based on the use of additional non-dimensional parameters, in particular, the design factor and the solar multiple factor. This paper presents a mathematical optimization model that focuses on the capacity of concentrated solar power plants where thermal storage plays a key role in the energy source. The analytical approach provides a more complete understanding of the design process for hybrid power plants. In addition, the use of additional factors and the combination of the proposed method with existing numerical methods allows for more refined optimization, which allows for the more accurate selection of the capacity for specific geographical conditions. Importantly, the proposed method significantly increases the speed of computation compared to that of traditional numerical methods. Finally, the authors present the results of the analysis of the proposed system of equations for calculating the levelized cost of electricity (LCOE) for hybrid solar power plants. The nonlinearity of the LCOE on the main calculation parameters is shown.

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 N D N I a v a i l 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
E y e a r n e e d = E y e a r C S P + E y e a r P V ,
E y e a r C S P = 0 τ P C S P d t = P C S P m a x h C S P F L H
where τ —the time period considered, usually τ = 1   y e a r ;
P C S P m a x —the necessary electric power to fulfill the demand schedule from CSP;
h C S P F L H —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 P C S P m a x , we obtain the imaginary number of hours h max that the CSP should operate to obtain E y e a r n e e d . If a CSP is designed with less than the required capacity for a given area P C S P d e s i g n , 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 P C S P d e s i g n and possibly three periods of plant shutdown. P C S P d e s i g n —the normalized design capacity of the CSP. P h y b r i d = P n e e d p e a k —the maximum required electric output power of the hybrid plant from demand schedule. For CSP without PV, P C S P m a x = P n e e d p e a k .
The full-load hours are calculated as follows:
h C S P F L H = E y e a r C S P P C S P m a x
The CSP electricity generation per year is calculated as follows:
E y e a r d e s i g n = P C S P d e s i g n h C S P F L H
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]:
L C O E = T L C C y = 1 Y E y 1 + i y
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
L C O E = C R F C A P E X + O P E X E y e a r
The capital recovery factor is as follows:
C R F = 1 + i Y i 1 + i Y 1
A weighted average cost of capital method [29] for the discount rate is used to depict the costs of equity i e q u and debt i d e b t less the income taxes:
W A C C = 1 f d e b t i e q u + f d e b t 1 r t a x , i n c i d e b t .
To account for further tax costs, the fixed charge rate FCR, according to [30], is introduced:
F C R = C R F 1 r t a x , i n c P r V d e p r 1 r t a x , i n c ,
where P r V d e p r —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:
P r V d e p r = y = 1 Y d e p r 1 Y d e p r 1 + W A C C y
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 Z C S P f a i l and costs due to penalties Z C S P p e n a l t i e s :
L C O E = F C R C A P E X + O P E X E y e a r + Z C S P f a i l + Z C S P p e n a l t i e s m i n
According to the life cycle assessment performed in reference [31], the operation phase decreases with increasing capital expenditures in CSP. So O P E X = f α C S P · N D N I a v a i l . 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:
S M = Q ˙ S F m a x Q ˙ P B d e s i g n
where Q ˙ S F m a x —the thermal power output of the solar field at the reference DNI value;
Q ˙ P B d e s i g n —the design nominal electric output power CSP according to its efficiency.
For the averaged values,
S M = P C S P m a x P C S P d e s i g n η C S P d e s i g n η C S P m a x
where η C S P m a x and η C S P d e s i g n —the effectiveness CSP for P C S P m a x and P C S P d e s i g n .
η C S P —efficiency CSP. In this study, it is assumed to be linearly dependent on the receiver efficiency. Then, according to [34], η C S P = a · l n f j + b . Here, a, b—const a > 0, b < 1. Hence, η C S P does not depend on the CSP capacity.
The amount of electricity to be produced by the CSP plant is given by
E y e a r C S P = E y e a r d e s i g n + E y e a r T E S
where E y e a r T E S —the energy from thermal storage. So
E y e a r C S P = P C S P d e s i g n h C S P F L H + E y e a r T E S
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. D F 0 , 1 is defined as
D F = P C S P m a x P n e e d p e a k
where P C S P m a x —the design nominal electric output power of the CSP according to its efficiency.
P n e e d p e a k —the maximum required electric output power of the hybrid plant from the demand schedule.
If P n e e d m a x = I P C S P m a x D F , the maximum electric output power of the CSP plant is required from the demand schedule (the design power of the entire power plant), so
E y e a r n e e d = P C S P d e s i g n h C S P F L H + E d a y T E S d a y + 1 D F P n e e d p e a k h C S P F L H
where E d a y T E S —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
C P T E S = Q d a y T E S Q ˙ P B d e s i g n = Q d a y T E S η p o w P C S P d e s i g n
where η p o w —the efficiency of the power block;
Q d a y T E S —the thermal energy storage capacity in MWh/day.
Q d a y T E S = E d a y T E S / η T E S
E y e a r T E S = E d a y T E S h C S P F L H 24
Therefore,
E y e a r n e e d = P C S P d e s i g n h h y b r i d F L H + C P T E S η p o w P C S P d e s i g n η T E S h h y b r i d F L H 24 + P C S P m a x 1 D F D F h h y b r i d F L H ;
E y e a r n e e d = P C S P d e s i g n h h y b r i d F L H 1 + η T E S η p o w C P T E S 1 24 + S M η C S P m a x η C S P d e s i g n 1 D F D F ;
E y e a r C S P + T E S = P C S P d e s i g n h C S P F L H 1 + η T E S η p o w C P T E S 1 24
h C S P F L H = h C S P a v a i b l e , if   h C S P F L H > h C S P a v a i b l e E y e a r n e e d P C S P d e s i g n 1 + η T E S η p o w C P T E S 1 24 , if   h C S P F L H < h C S P a v a a b l e
This equation shows the relationship between the required energy and the designed capacity of the CSP output unit through the parameters C P T E S ,   S M ,     α C S P and the efficiency of the plant. Assume that the dependence of the CSP efficiency on the receiver efficiency is known. Then, P C S P d e s i g n   is dependent on three parameters, C P T E S ,   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 P C S P m a x is
A C S P m a x = S M · A C S P d e s i g n = E y e a r C S P 0 n D N I n d h η C S P d e s i g n
where d h —the hours of DNI data resolution.
Assuming a constant efficiency for the normalized available operating hours, the following equation results:
h C S P F L H = 0 n D N I n d h S M A C S P d e s i g n P C S P m a x η C S P d e s i g n
where A C S P d e s i g n —the area of the solar field for reaching P C S P d e s i g n e .
To consider the minimum and maximum capacity of the power generation system during deviations from the selected average value,
D N I n = 0   f o r   D N I n   10 % D N I d e s i g n D N I n       f o r   10 % D N I d e s i g n P C S P n 120 %   D N I d e s i g n 120 %   D N I d e s i g n         f o r     D N I n   120 %   D N I d e s i g n  
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:
C P T E S = E y e a r C S P E y e a r d e s i g n η p o w P C S P d e s i g n η T E S 24 h C S P F L H
And Equation (18) is used:
C P T E S = S M A C S P d e s i g n 0 n D N I n d h η C S P d e s i g n P C S P d e s i g n h C S P F L H η p o w P C S P d e s i g n η T E S 24 h C S P F L H
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:
C P T E S = P C S P m a x η C S P d e s i g n h C S P F L H P C S P d e s i g n h C S P F L H η p o w P C S P d e s i g n η T E S 24 h C S P F L H
C P T E S = 24 P C S P m a x η C S P d e s i g n P C S P m a x S M η C S P d e s i g n η C S P m a x η p o w P C S P m a x S M η C S P d e s i g n η C S P m a x η T E S
C P T E S = 24 S M 1 η p o w η T E S
Given that the dependence A C S P d e s i g n is determined by the type and internal characteristics of the designed CSP, it is assumed that A C S P d e s i g n P C S P d e s i g n   r e f 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
h C S P a v a i l a b l e = h C S P F L H f f j , h s d , D N I r e f
where D N I r e f —the reference DNI for the normalized value for operating hours;
h s d —the sunshine duration, h/year;
f j —the relative load corresponding to the average arrival of solar radiation [36].
f j = e x p β D j D m i n D m a x D m i n j = 1 J e x p β D j D m i n D m a x D m i n
where D m a x   and   D m i n 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:
L C O E = F C R C A P E X + O P E X P C S P d e s i g n h C S P F L H 1 + η T E S η p o w C P T E S 1 24 + S M η C S P m a x η C S P d e s i g n 1 D F D F + Z C S P f a i l + Z C S P p e n a l t i e s m i n
where Z C S P f a i l —mathematical expectation of the specific impaired value due to losses caused by the flow of failures in the network;
Z C S P p e n a l t i e s —the cost and the penalty function.
The CAPEX can be divided into overnight capital costs OCC and costs for network connections C t r a n s i t as well as costs for financing the construction C f i n .
C A P E X = O C C + C t r a n s i t + C f i n
The EPC costs are divided into component costs, land costs including preparation, and service costs and are extended by a contingency factor f c o n t . The OCC cost components all depend on the size of the plant and therefore on D F . The scaling can be performed using an exponential economy of scale relationship, as described in [38].
O C C = O C C C S P + O C C P V ;
O C C C S P = ( C c o m p m a x D F s c o m p + C s e r v i c e m a x D F s s e r v i c e + C l a n d m a x     D F ) f c o n t
O C C P V = ( C c o m p m a x 1 D F + C s e r v i c e m a x 1 D F + C l a n d m a x     1 D F )   f c o n t
where Ccomp—capital investment for all components and equipment [€];
Cservice—capital investment for all engineering, procurement and construction services [€];
C l a n d —capital investment to purchase and prepare the land [€];
D F —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].
C s e r v i c e m a x = C c o m p m a x   f s e r v i c e
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:
C t r a n s i t = K m a x t r   n   η t r   L B + Z m i n t r     E y e a r n e e d
where K m a x t r —capital investments in electricity transit networks, €/m;
n—the number of parallel strands of wire 1 or 2;
η t r —the network performance indicator;
L B —the transit network length, m;
Z m i n t r —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 r f e e s , for the interest during construction i c o n s t r , assumed to be accountable for half of the construction period N c o n s t r , and additional money loaned for reserve accounts r r e s :
C f i n = f d e b t O C C + C t r a n s i t ( r f e e s + i c o n s t r N c o n s t r 2 + r r e s )
The annual operational expenditures OPEX are divided into fixed and variable costs C f i x and C v a r , both in €, and are introduced as [39]
O P E X = C f i x + C v a r
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 P C S P d e s i g n :
C f i x = c f i x P C S P d e s i g n
where c f i x —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
C v a r = c v a r E y e a r
where c v a r —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:
Z C S P f a i l = C S C P n m i r r o r 1 h C S P m a x   1 exp ω m i r r o r n m i r r o r h C S P F L H h y e a r + C P V n P V 1 h P V m a x × 1 exp ω P V n P V h P V m a x h y e a r
where C C S P   ,   C P V —the cost of repairing a mirror unit or PV panel €/(pcs);
n m i r r o r ,   n P V —the number of mirrors or PV panels;
ω m —the mirror failure rate, 1/(year);
h P V m a x —the number of hours of PV operation per year;
A m i r r o r —the single mirror area.
Or
Z C S P f a i l = S M C m A C S P d e s i g n A m i r r o r 1 h C S P F L H   1 exp ω m i r r o r n m i r r o r h C S P F L H h y e a r + C P V 1 D F E y e a r n e e d P P V p c s 1 h P V m a x × 1 exp ω P V n P V h P V m a x h y e a r
where P P V p c s —the power of one panel, €/MWt.
Using this equation, it should be noted that α C S P 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 C C S P p e n a l t i e s = 7.69   / M W h   , given the theoretical support from the government for the introduction of a green tariff at which Z C S P p e n a l t i e s < 0 , so C Z can be negative.
The target function has been calculated as
L C O E = F C R C A P E X + O P E X E y e a r + S M C m A C S P d e s i g n A m i r r o r 1 h C S P F L H 1 exp ω m i r r o r n m i r r o r   h C S P F L H h y e a r + C P V 1 D F E y e a r n e e d P P V p c s 1 h P V m a x 1 exp ω P V n P V   h P V m a x h y e a r m i n
E y e a r p r o d u c e = P C S P d e s i g n h C S P F L H 1 + η T E S η p o w C P T E S 1 24 + S M η C S P m a x η C S P d e s i g n 1 D F D F ;
h C S P F L H = h C S P a v i a b l e , i f h C S P F L H > h C S P a v i a b l e E y e a r n e e d P C S P d e s i g n 1 + η T E S η p o w C P T E S 1 24 , i f h C S P F L H < h C S P a v i a b l e ;
C P T E S = S M A C S P d e s i g n 0 n D N I n d h η C S P d e s i g n P C S P d e s i g n h C S P F L H η p o w P C S P d e s i g n η T E S 24 h C S P F L H f o r S M = c o n s t 24 1 1 S M 1 η C S P m a x η p o w 1 S M 1 η C S P m a x η T E S f o r Q ˙ P B d e s i g n = c o n s t
P C S P m a x = D F   P n e e d m a x
P C S P d e s i g n = D F   P n e e d m a x   S M η C S P d e s i g n η C S P m a x ;
C A P E X = O C C C S P + O C C P V + C t r a n s i t + C f i n ;
O C C C S P = ( C c o m p m a x D F 0.8 + C s e r v i c e m a x 0.05 D F 0.4 + C l a n d m a x     D F )   f c o n t ;
O C C P V = ( C c o m p m a x 1 D F 0.8 + C s e r v i c e m a x 0.05 1 D F 0.4 + C l a n d m a x     1 D F )   f c o n t ;
O P E X = C f i x + c v a r _ C S P D F E y e a r n e e d + c v a r P V 1 D F E y e a r n e e d .
This function has three variables, S M ,   P C S P d e s i g n   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 α C S P = 0 . 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 Q ˙ P B d e s i g n 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.

Author Contributions

Conceptualization, A.C.; Supervision, A.C.; Methodology, A.C.; Investigation, A.C., O.K.; Writing—original draft, A.C.; Writing—review and editing, S.A.; Formal analysis, A.P.; Validation, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Alqallaf, N.; Ghannam, R. Immersive Learning in Photovoltaic Energy Education: A Comprehensive Review of Virtual Reality Applications. Solar 2024, 4, 136–161. [Google Scholar] [CrossRef]
  2. Wieghardt, K. Concentrating Solar Thermal Technologies—Status and Perspectives. In Proceedings of the EST2015 Energy Science Technology Conference, Karlsruhe, Germany, 20–22 May 2015. [Google Scholar]
  3. Ackermann, S.; Michael, T.; Cheffe, J.A. Steinfeld Reticulated porous ceria undergoing thermochemical reduction with high-flux irradiation. Int. J. Heat Mass Transf. 2017, 107, 439–449. [Google Scholar] [CrossRef]
  4. Teleszewski, T.J.; Żukowski, M.; Krawczyk, D.A.; Rodero, A. Analysis of the Applicability of the Parabolic Trough Solar Thermal Power Plants in the Locations with a Temperate Climate. Energies 2021, 14, 3003. [Google Scholar] [CrossRef]
  5. Benghanem, M.; Mellit, A.; Almohamadi, H.; Haddad, S.; Chettibi, N.; Alanazi, A.M.; Dasalla, D.; Alzahrani, A. Hydrogen Production Methods Based on Solar and Wind Energy: A Review. Energies 2023, 16, 757. [Google Scholar] [CrossRef]
  6. Kim, Y.; Skaggs, A. A Multi-Stage Approach to Assessing the Echo-Tech Feasibility of a Hybrid SAM-CREST Model for Solar PV Power Plants in Maryland, USA. Solar 2024, 4, 246–268. [Google Scholar] [CrossRef]
  7. Ait Lahoussine Ouali, H.; Alami Merrouni, A.; Chowdhury, S.; Techato, K.; Channumsin, S.; Ullah, N. Optimization and Techno-Economic Appraisal of Parabolic Trough Solar Power Plant under Different Scenarios: A Case Study of Morocco. Energies 2022, 15, 8485. [Google Scholar] [CrossRef]
  8. Montes, M.J.; Abánades, A.; Martínez-Val, J.M.; Valdés, M. Solar multiple optimization for a solar-only thermal power plant, using oil as heat transfer fluid in the parabolic trough collectors. Solar Energy 2009, 83, 2165–2176. [Google Scholar] [CrossRef]
  9. De la Fuente, E.; Martin, M. Site specific process design for Hybrid CSP-Waste plants. Comput. Chem. Eng. 2020, 135, 106770. [Google Scholar] [CrossRef]
  10. Wang, Q.; Chang, P.; Bai, R.; Liu, W.; Dai, J.; Tang, Y. Mitigation Strategy for Duck Curve in High Photovoltaic Penetration Power System Using Concentrating Solar Power Station. Energies 2019, 12, 3521. [Google Scholar] [CrossRef]
  11. Ferretti, F.; Hogendijk, C.; Aga, V.; Ehrsam, A. Addressing Forecast Uncertainty Impact on CSP Annual Performance. AIP Conf. Proc. 2017, 1850, 030015. [Google Scholar] [CrossRef]
  12. Moghimi, M.A.; Craig, K.J.; Meyer, J.P. Optimization of a trapezoidal cavity absorber for the Linear Fresnel Reflector. Sol. Energy 2015, 119, 343–361. [Google Scholar] [CrossRef]
  13. Iñigo-Labairu, J.; Dersch, J.; Hirsch, T.; Giuliano, S.; Loevenich, M.; Córdoba, D. Techno-Economic Evaluation of CSP–PV Hybrid Plants with Heat Pump in a Temperature Booster Configuration. Energies 2024, 17, 2634. [Google Scholar] [CrossRef]
  14. Rouibah, A.; Benazzouz, D.; Kouider, R.; Al-Kassir, A.; García-Sanz-Calcedo, J.; Maghzili, K. Solar Tower Power Plants of Molten Salt External Receivers in Algeria: Analysis of Direct Normal Irradiation on Performance. Appl. Sci. 2018, 8, 1221. [Google Scholar] [CrossRef]
  15. Shakeel, M.R.; Mokheimer, E.M.A. A techno-economic evaluation of utility scale solar power generation. Energy 2022, 261, 125170. [Google Scholar] [CrossRef]
  16. Djebbar, R.; Belanger, D.; Boutin, D.; Weterings, E.; Poirier, M. Potential of concentrating solar power in Canada. Energy Procedia 2014, 49, 2303–2312. [Google Scholar] [CrossRef]
  17. Dersch, J.; Schwarzbozl, P.; Richert, T. Annual Yield Analysis of Solar Tower Power Plants with GREENIUS. J. Sol. Energy Eng.-Trans. ASME 2011, 133, 031017-1–031017-9. [Google Scholar] [CrossRef]
  18. Inigo-Labairu, J.; Dersch, J.; Schomaker, L. Integration of CSP and PV Power Plants: Investigations about Synergies by Close Coupling. Energies 2022, 15, 7103. [Google Scholar] [CrossRef]
  19. Blair, N.; DiOrio, N.; Freeman, J.; Gilman, P.; Janzou, S.; Neises, T.; Wagner, M. System Advisor Model (SAM) General Description (Version 2017.9.5); NREL/TP-6A20-70414; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2018. [Google Scholar]
  20. Alexopoulos, S. Modeling and Simulation of the 50 MW molten salt solar tower power plant in Delingha Real tower power plant simulation with SAM. In Proceedings of the SolarPACES 2023, Sydney, Australia, 10–13 October 2023. [Google Scholar]
  21. Riffelmann, K.-J.; Weinrebe, G.; Balz, M. Hybrid CSP-PV Plants with Integrated Thermal Storage. In Proceedings of the SolarPACES 2021, Online, 27 September–1 October 2021. [Google Scholar]
  22. Hasni, S.; Platzer, W.J. Case Study on Decarbonization Strategies for LNG Export Terminals Using Heat and Power from CSP/PV Hybrid Plants. In Proceedings of the SolarPACES 2022, Albuquerque, NM, USA, 26–30 September 2022. [Google Scholar]
  23. Hamilton, W.T.; Martinek, J.; Cox, J.; Newman, A. Modeling and Optimizing CSP-PV Hybrid Systems Using the Hybrid Optimization and Performance Platform (HOPP). In Proceedings of the SolarPACES 2022, Albuquerque, NM, USA, 26–30 September 2022. [Google Scholar]
  24. Tripp, C.E.; Guittet, D.; Barker, A.; King, J.; Hamilton, B. Hybrid Optimization and Performance Platform (HOPP). [Computer Software]. 14 November 2019. Available online: https://github.com/NREL/HOPP (accessed on 1 July 2024). [CrossRef]
  25. Dersch, J.; Inigo Labairu, J.; Hirsch, T. Greenius—A free Software Tool for Simulating Electricity and Heat Generation Systems with Concentrating and Non-Concentrating Solar Collectors. In Proceedings of the Cologne Solar Colloquium 2024, Cologne, Germany, 26 June 2024. [Google Scholar]
  26. Green, A.; Diep, C.; Dunn, R.; Dent, J. High Capacity Factor CSP-PV Hybrid Systems. Energy Procedia 2015, 69, 2049–2059. [Google Scholar] [CrossRef]
  27. Trevisan, S.; Guédez, R.; Laumert, B. Techno-economic Analysis of Hybrid Supercritical CO2 Based Molten Salts Tower CSP-PV. In Proceedings of the SolarPACES 2021, Online, 27 September–1 October 2021. [Google Scholar]
  28. Gasa, G.; Prieto, C.; Lopez-Roman, A.; Cabeza, L.F. Life cycle assessment (LCA) of a concentrating solar power (CSP) plant in tower configuration with different storage capacity in molten salts. J. Energy Storage 2022, 53, 105219. [Google Scholar] [CrossRef]
  29. Gedle, Y.; Schmitz, M.; Schmitz, P.; Herrmann, U.; Boura, C.; Mahdi, Z.; Chico, C.; Ricardo, A.; Merige, P.; Dersch, J. Analysis of an Integrated CSP-PV hybrid power plant. In Proceedings of the 26th SolarPACES Conference 2020, Albuquerque, NM, USA, 28 September–2 October 2020. [Google Scholar]
  30. Short, W.; Packey, D.; Holt, T. A Manual for the Economic Evaluation of Energy Efficiency and Renewable Energy Technologies; National Renewable Energy Laboratory: Golden, Colorado, 1995. [Google Scholar]
  31. Kumar, S.; Agarwal, A.; Kumar, A. Financial viability assessment of concentrated solar power technologies under Indian climatic conditions. Sustain. Energy Technol. Assess. 2021, 43, 100928. [Google Scholar] [CrossRef]
  32. Coelho, B.; Varga, S.; Oliveira, A.; Mendes, A. Optimization of an atmospheric air volumetric central receiver system: Impact of solar multiple, storage capacity and control strategy. Renew. Energy 2014, 63, 392–401. [Google Scholar] [CrossRef]
  33. Khosravi, A.; Malekan, M.; Pabon, J.J.G.; Zhao, X.; Assad, M.E.H. Design parameter modelling of solar power tower system using adaptive neuro-fuzzy inference system optimized with a combination of genetic algorithm and teaching learning-based optimization algorithm. J. Clean. Prod. 2020, 244, 118904. [Google Scholar] [CrossRef]
  34. Rosenstiel, A.; Monnerie, N.; Dersch, J.; Roeb, M.; Pitz-Paal, R.; Sattler, C. Electrochemical Hydrogen Production Powered by PV/CSP Hybrid Power Plants: A Modelling Approach for Cost Optimal System Design. Energies 2021, 14, 3437. [Google Scholar] [CrossRef]
  35. Cheilytko, A. Vykorystannya Vtorynnykh Enerhoresursiv: Navchalʹno–Metodychnyy Posibnyk Dlya Studentiv ZDIA Napryamu 6.050601 “Teploenerhetyka” Dennoyi i Zaochnoyi Formy Navchannya; ZDIA: Zaporizhzhia, Ukraine, 2014. [Google Scholar]
  36. Mohammadzadeh, N.; Truong-Ba, H.; Cholette, M.E.; Steinberg, T.A.; Manzolini, G. Model-predictive control for dispatch planning of concentrating solar power plants under real-time spot electricity prices. Sol. Energy 2022, 248, 230–250. [Google Scholar] [CrossRef]
  37. Segarra, E.L.; Ruiz, G.R.; Bandera, C.F. Probabilistic Load Forecasting Optimization for Building Energy Models via Day Characterization. Sensors 2021, 21, 3299. [Google Scholar] [CrossRef]
  38. Hirsch, T.; Mehos, M. SolarPACES Guideline for Bankable STE Yield Assessment. 2017. Available online: https://www.solarpaces.org/guideline-for-bankable-ste-yield-assessment/ (accessed on 15 March 2024).
  39. Ho, J.; Becker, J.; Brown, M.; Brown, P.; Chernyakhovskiy, I.; Cohen, S.; Cole, W.; Corcoran, S.; Eurek, K.; Frazier, W.; et al. Regional Energy Deployment System (ReEDS) Model Documentation: Version 2020; National Renewable Energy Laboratory: Golden, CO, USA, 2021; NREL/TP-6A20-78195. Available online: https://www.nrel.gov/docs/fy21osti/78195.pdf (accessed on 1 July 2024).
  40. Vasallo, M.; Bravo, J.; Marín, D.; Gegúndez, M. Economic MPC applied to generation scheduling in CSP plants. IFAC-PapersOnLine 2017, 50, 115–120. [Google Scholar] [CrossRef]
Figure 1. Schematic explanation of the calculation of h C S P F L H   by different types of hybrid power plant elements: yellow is the planned electricity produced by CSP, blue is the energy produced for thermal storage. The red line is the planned mode of operation. No color area—CSP+PV.
Figure 1. Schematic explanation of the calculation of h C S P F L H   by different types of hybrid power plant elements: yellow is the planned electricity produced by CSP, blue is the energy produced for thermal storage. The red line is the planned mode of operation. No color area—CSP+PV.
Solar 04 00024 g001
Figure 2. Comparison of analytical and numerical approach. The numerical approach overestimates the LCOE by about 6% but reflects the cost increase for higher design factors.
Figure 2. Comparison of analytical and numerical approach. The numerical approach overestimates the LCOE by about 6% but reflects the cost increase for higher design factors.
Solar 04 00024 g002
Figure 3. Levelized cost of electricity as a function of design factor with BESSs for PV shows higher LCOE for all DF values.
Figure 3. Levelized cost of electricity as a function of design factor with BESSs for PV shows higher LCOE for all DF values.
Solar 04 00024 g003
Figure 4. Share of electricity generated. Solar multiple (SM) = 1.3 with BESS for PV.
Figure 4. Share of electricity generated. Solar multiple (SM) = 1.3 with BESS for PV.
Solar 04 00024 g004
Figure 5. Dependence of levelized cost of electricity (LCOE) on design factor for different values of solar multiple (SM). The minimum LCOE is achieved at SM =1.3.
Figure 5. Dependence of levelized cost of electricity (LCOE) on design factor for different values of solar multiple (SM). The minimum LCOE is achieved at SM =1.3.
Solar 04 00024 g005
Figure 6. Dependence of the levelized cost of electricity from the solar multiple (SM) with design factors (DF) = 0.5. A clear extreme of the function is expressed.
Figure 6. Dependence of the levelized cost of electricity from the solar multiple (SM) with design factors (DF) = 0.5. A clear extreme of the function is expressed.
Solar 04 00024 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cheilytko, A.; Alexopoulos, S.; Pozhuyev, A.; Kaufhold, O. An Analytical Approach to Power Optimization of Concentrating Solar Power Plants with Thermal Storage. Solar 2024, 4, 509-525. https://doi.org/10.3390/solar4030024

AMA Style

Cheilytko A, Alexopoulos S, Pozhuyev A, Kaufhold O. An Analytical Approach to Power Optimization of Concentrating Solar Power Plants with Thermal Storage. Solar. 2024; 4(3):509-525. https://doi.org/10.3390/solar4030024

Chicago/Turabian Style

Cheilytko, Andrii, Spiros Alexopoulos, Andriy Pozhuyev, and Oliver Kaufhold. 2024. "An Analytical Approach to Power Optimization of Concentrating Solar Power Plants with Thermal Storage" Solar 4, no. 3: 509-525. https://doi.org/10.3390/solar4030024

APA Style

Cheilytko, A., Alexopoulos, S., Pozhuyev, A., & Kaufhold, O. (2024). An Analytical Approach to Power Optimization of Concentrating Solar Power Plants with Thermal Storage. Solar, 4(3), 509-525. https://doi.org/10.3390/solar4030024

Article Metrics

Back to TopTop