1. Introduction
With worldwide commitment to the ambitious goal of carbon neutrality, the penetration of renewable energy will inevitably increase [
1,
2]. However, the electrical grid system, under high penetration of renewable energy, will face many unavoidable challenges due to the strong uncertainty of renewable sources, such as security, stability and reliability [
3,
4]. In order to alleviate the impact of uncertainties on the grid system, exploiting controllable power generation plants to maintain the quality requirements of grid system is a necessity. Among them, flexible operation of conventional fossil fuel plants has gained unprecedented increasing attention. However, affordable accommodation for renewable energy of a conventional coal power plant is usually limited at 50% of the rated capacity because of various physical constraints in energy conversion processes [
5]. In addition, the fossil fuel plants produce a large amount of carbon dioxide, causing their decommission in some countries [
6]. There are many alternative power generation plants. Solar power plants, as a commercially mature technology, is one such alternative, and has even taken a large market share in the grid system. For this reason, this paper attempts to investigate the possibility of solar power plants, specifically the concentrating solar power plant with thermal energy storage (CSP-TES) system, to operate flexibly in the regulation of the electrical grid system.
The CSP-TES system has an interesting characteristic of dispatchability, essentially due to the coordination function of the thermal energy storage subsystem. This characteristic has been mentioned by several works. For instance, ref. [
7] examined the dispatchability by using two country cases of Morocco and Algeria, and found that cost advantages of CSP with TES should be highlighted when the gird system is in high renewable penetration. Ref. [
8] also concluded that CSP with TES can provide dispatchable energy and power based on the operation date of Andasol 3, and such an advantage enables the CSP to have many attractive functions. Based on the review of CSP plant development and operation market rules in eight countries, ref. [
9] concluded that CSP’s dispatchability and flexibility should be highly valued and encouraged by markets and policies. The influences of dispatchability on the CSP system are also investigated in many works [
10,
11,
12].
However, the dispatchability of CSP-TES is mostly applied to economic dispatch in existing works, where the basic idea is to firstly establish model constraints according to equipment characteristics, and then optimize a specific objective function and calculate the schedule of the CSP-TES via various operation strategies. For instance, a mixed-integer linear problem was also established by Kost et al. to maximize both the operational profit and investment annuity of a CSP plant over the planning horizon [
13]. Similarly, a mixed-integer programming model of a CSP-TES plant was proposed in ref. [
10]. In addition, a model predictive control operation strategy was also exploited to minimize the objective composed of tracking the schedule and generating the optimal schedule for the next day. Studies of optimal operation strategy of CSP-TES can also be found in many works [
14,
15,
16].
In fact, the optimal schedules of CSP-TES have to be delivered to controllers of subsystems. In existing works, the controllers are designed mainly from a subsystem level in a decoupling way, as the coupling effects among subsystems are considered in the scheduling layer. Extensive works based on the decoupling way have been reported. For instance, a nonlinear model-based adaptive control strategy has been proposed to control the outlet temperature of a parabolic trough plant [
17]. An optimal model predictive control with lump disturbance to maintain the outlet temperature has been proposed [
18]. A linear quadratic regulator was designed to regulate the output oil enthalpy of the thermal storage subsystem in ref. [
19]. This idea of controller design can be found in many studies [
20,
21].
However, whether optimal dispatch or controller design of the studies as previously mentioned, they are all carried out without considering CSP participation in the regulation of the grid system. Therein, CSP can calculate the best performance itself through optimal operation strategies, and then approach this self-optimization by tracking the optimized instructions, without considering the contribution to the performance of the grid system.
In contrast, the work here studies the flexible operation of CSP-TES and investigates its potential function in grid regulation, based on its dispatchable capability. Different from the traditional approach, when the CSP-TES participates in the grid regulation, the power instruction is not obtained from its own optimization, but delivered from the grid operator. This indicates that the role of the CSP-TES transfers from independent optimal operation to a clean regulator of the grid system. Regarding the capability of the role transition of CSP-TES, it has been claimed by [
8] that CSP can offer similar operational attributes as conventional power plants, and can also support regulation and frequency response. However, the similarities and differences between CSP and conventional plants, and how to realize the flexible regulation of CSP-TES, have not yet been studied. To this end, dynamic characteristics of the CSP-TES system is investigated based on the dynamic model in ref. [
22], and a coordinated control strategy is further proposed to realize the flexible operation of CSP-TES. To the authors’ best knowledge, research on CSP-TES from this perspective has not yet been conducted.
This work is a further extension of our previous work [
22], and organized as follows.
Section 2 examines the dynamic characteristics of CSP-TES in ref. [
22], and analyzes its control-oriented characteristics for flexible operation. Afterwards, a coordinated control strategy is designed in
Section 3.
Section 4 testifies the capacity of CSP-TES to operate flexibly in grid regulation by case studies even under wide-range operation and subjected to strong uncertain irradiation. Finally, conclusions are drawn in
Section 5.
2. Dynamic Simulation and Control Characteristic Analysis
An in-depth understanding of the controlled system is critical to controller design.
Figure 1 shows the schematic diagram of a parabolic trough CSP with TES, where it consists of solar field (SF), TES and power block (PB) subsystems. In control practice, the SF outlet temperature of heat transfer fluid (HTF) and its efficiency are dominated by the HTF flowrate of SF; electricity generation and main steam pressure of PB are controlled by HTF flowrate and the opening valve of the PB subsystem; the thermal charging/discharging of TES is manipulated by the MS flowrate of TES; and solar irradiation is the energy input of the CSP system. Consequently, typical step change in process control engineering is performed to evaluate the control characteristics of the CSP-TES system. In this work, the CSP-TES is assumed on a sunny day where the TES subsystem is under charging state.
For this system, we have established a system-level analytical model in the form of the ordinary differential equations in ref. [
22], whose accuracy has been verified by the steady-state and dynamic tests. Therefore, it can be reliably applied to design controllers of CSP plants. Then, based on the analytical model and with the given stable operating parameters in
Table 1, the system step response is analyzed as follows.
2.1. System Response to Step Change of HTF Flowrate of SF
Figure 2 shows the system response to step change of HTF flowrate of SF. As can be seen, the HTF flowrate of SF has a strong effect on the controlled variables of the system. When the HTF flowrate increases, the SF outlet temperature decreases gradually as the input solar energy remains. The decrease of SF outlet temperature further leads to the reduction in electricity and main steam pressure. Because the HTF flowrate to PB remains at 500 kg/s, the lower HTF temperature to PB weakens the heat transfer of heat exchangers, which causes less energy to be transferred to steam. Conversely, the MS temperature increases first as the HTF flowrate of SF steps. This increase is attributed to more HTF flowrate delivered to the TES subsystem. However, as the SF outlet HTF temperature further decreases, the heat transfer efficiency of the HTF-MS heat exchanger degrades, which leads to the decline of MS temperature. However, the responses are not monotonous, as all the chosen variables climb again after about 900 s. Such behaviors maybe be due to the increase in inlet HTF temperature of SF, and this increase contributes to the increases in SF outlet HTF temperature. With the SF outlet HTF temperature increasing, the MS outlet temperature, the electricity and the main steam pressure climb smoothly to new steady states.
2.2. System Response to Step Change of HTF Flowrate of PB
As the inlet HTF flowrate of PB increases, more steam is generated because the heat transfer capacity in the heat exchanger train is enhanced, thereby resulting in a monotonous increase in the electricity, together with the steam pressure. These indicate that the response characteristics of CSP with TES are similar to conventional plants. Meanwhile, the change of inlet HTF flowrate of PB also affects the SF outlet HTF temperature and MS temperature, as proven in
Figure 3. However, the MS outlet temperature responds quickly, while the HTF temperature changes later, with a time delay of about 500 s. This is because the decrease of the HTF flowrate of TES subsystem leads to the increase of the HTF flowrate of PB, under the constant total HTF flowrate. Such a decrease leads to the drop in the MS outlet temperature in the TES. For SF outlet HTF temperature, although the total HTF flowrate in the SF is the same, the flowrate redistribution changes the SF inlet HTF temperature, and ultimately affects the SF outlet HTF temperature after a transport delay.
2.3. System Response to Step Change of Main Valve Opening
Figure 4 illustrates the system response to main valve opening steps. As can be seen, an increase in the opening of the main valve will instantly lead to an increase in electricity and a decrease in steam pressure, which is similar to conventional plants. However, the influence on MS temperature and SF outlet HTF temperature are not revealed because the changes of the main valve have a minor direct influence on the SF inlet HTF temperature, and indirect impact on the outlet temperature of SF and MS; nevertheless, the influence is very limited.
2.4. System Response to Step Change of MS Flowrate of TES
The step change of MS flowrate of TES also has an impact on system dynamics, as illustrated in
Figure 5. When the MS flowrate steps, the MS temperature will drop instantly; however, the SF outlet HTF temperature decreases after a certain time delay. The reason is similar to that of system step response to HTF flowrate of PB; the two cases both reduce SF inlet HTF temperature directly, and ultimately reduce the SF outlet HTF temperature after a certain transport delay. The changes in SF outlet HTF temperature further impact on the performance of the PB subsystem, despite the impact being limited.
2.5. System Response to Step Change of Solar Irradiation
Solar irradiation is the input energy to CSP and has the features of uncertainty and uncontrollability. To this end, dynamic response of the CSP system to irradiation is examined in
Figure 6. It can be seen from
Figure 6 that step change of irradiation has a significant impact on the overall CSP system. With about 7.1% increase of solar irradiation, the SF outlet temperature and MS temperature both increase by around 8 °C, and the electricity and main steam pressure increase by 7.1% and 3.8%, respectively. Since solar energy is intermittent and fluctuates strongly, the dynamic responses indicate the importance of disturbance rejection control. Moreover, solar energy has a very broad influence on system performance, as its increase leads to an increase in the SF outlet HTF temperature, and as a result, the performance of the downstream subsystem changes greatly. The slow decrease also suggests the large inertia of the CSP system.
Based on the analysis above, it can be seen that there are many similarities between the dynamic characteristics of the CSP-TES and conventional fossil fuel plant. For instance, they are both a multivariable system with large thermal inertia and certain time delays. These similarities, to some degree, demonstrate the potential of CSP-TES to flexible operation. Nevertheless, there are also differences between CSP-TES and conventional plants.
- (1)
The CSP-TES system is a weak coupling system. Although the CSP-TES system comprises three subsystems, the dispatchability of the TES subsystem depresses the coupling degree between subsystems. The connection between SF and PB subsystems is dominantly the HTF temperature, as the HTF flowrate is determined automatically when HTF flowrate of SF and PB is controlled. Once the HTF temperature between SF and PB is controlled well, the SF and PB can be considered to be almost decoupled.
- (2)
The CSP-TES system is characterized by switch nonlinearity. The nonlinear behavior is determined by HTF flowrate of TES, which depends on the charging/discharging states of TES, while the states of TES are restricted by HTF flowrate between SF and PB. Under the assumption that TES is in charging state, the HTF flowrate of TES is equal to the HTF flowrate of SF minus the flowrate of PB, indicating that the flexible operation capacity of CSP system varies from that of conventional plants significantly.
- (3)
The CSP system is an uncertainty-dominated system. The largest uncertainty stems from the intermittence of solar energy, and the large variation in solar energy further deduces model mismatches and parameter variation, making the CSP system have strong uncertainties. Thus, from the perspective of disturbance rejection, it is necessary to attenuate the disturbance quickly from where it occurs to prevent solar energy fluctuations affecting the downstream subsystems. Consequently, disturbance rejection control of the CSP system should be a top priority.
3. Coordinated Control Strategy of CSP-TES System
According to the aforementioned control characteristics, and with the aim of flexible operation of the CSP-TES system, a coordinated control strategy was proposed, as shown in
Figure 7. Since the CSP-TES system is a weak coupling system, decentralized proportional-integral-derivative (PID) controllers are designed, providing that the controllers can maximize the potential engineering application due to their simple structure and high maturity. As for uncertainty alleviation, a disturbance observer-based (DOB) feedforward–feedback scheme is designed and applied to the SF subsystem to reject the disturbance of solar energy rapidly. For the PB subsystem, a decoupling controller is designed to decouple the electricity generation and main steam pressure, providing the similarities between the PB subsystem and traditional plants. The difference of HTF flowrate of SF and PB are fed forward to the TES, together with SF outlet HTF temperature, to compensate the switch nonlinearity. In addition, a feedback controller is designed to manipulate the MS flowrate of TES to control the MS temperature. Details of the control scheme are as follows.
3.1. DOB-Based Feedforward–Feedback Control Scheme of SF
Disturbance rejection is a common issue that exists in almost all control practices. In general, disturbance is categorized into internal disturbance, such as model mismatch and parameter perturbation, and external disturbance. Measurable external disturbance can usually be compensated by a feedforward strategy. However, the unmeasurable external and internal disturbance are challenging to fully eliminate. The CSP-TES system is a solar-driven system, which means there is unavoidable large external disturbance. Moreover, large variation in solar irradiation causes off-design operation conditions, resulting in certain internal disturbance. The large external and internal disturbances together put great challenges on the control of the CSP system. To solve the problem, a disturbance observer-based feedback–feedforward control strategy is proposed for the SF of the CSP system, of which the block diagram is illustrated in
Figure 8.
As can be seen, measurable solar radiation is compensated by the feedforward method to reduce its influence in time. However, the feedforward approach can only reject the measurable disturbance. In order to reduce the influence of unmeasurable disturbance, a disturbance observer is designed to obtain the equivalent disturbance. The disturbance estimated by the observer is then sent to the input as an additional control signal to compensate for the estimated disturbance. Finally, the outlet HTF temperature is fed back to the feedback controller to further boost the performance of disturbance attenuation.
Transfer functions GSF,ir and GSF,fl can be easily identified from the single-input single-output step response of the CSP-TES system under the open loop condition, accordingly. The issue is the design of the controllers FF1 and PI1, and the filter Q.
From the dynamic analysis, we found that
GSF,ir and
GSF,fl are relatively simple and can be identified as the first order process. In this way, the feedforward controller can be designed according to the invariance principle, namely
where
GSF,ir is the transfer function from irradiation to outlet HTF temperature, and
GSF,fl is the transfer function from HTF flowrate to outlet HTF temperature of
SF.
G−1 is the inverse of
GSF,fl.
FF1 is the feedforward controller,
PI1 is the feedback controller, and
Q is a filter.
Filter
Q is the core of the disturbance observer and determines the robustness and anti-interference of control performance. To guarantee physical realizability, its relative order must not be lower than
GSF,fl. In this work, the filter is designed in the following form:
where
β is designed according to ref. [
23].
Since the block diagram is established based on a linear system, the disturbance estimator and the feedback controller can be designed, respectively, according to the principle of superposition. Therefore, PI1 is designed based on the identified transfer function GSF,fl and tuned manually. With this control scheme, the disturbance on the SF outlet HTF temperature can be attenuated rapidly.
3.2. Decoupling Control Scheme of PB
Flexible operation of the CSP-TES system directly depends on the flexible response of PB to grid commands. As PB has similar characteristics to conventional plants, the coordination idea is stimulated where the HTF flowrate majorly contributes to the electricity tracking, while the main valve largely controls the main steam pressure. To achieve the control purpose, a decoupling control scheme is proposed here, as demonstrated in
Figure 9, where
G is the transfer function matrix,
D is the decoupling matrix, and
PI3 and
PI4 are the controllers.
Elements of the transfer function matrix can be identified from the step response of control variables, respectively, and also based on the open loop simulation of the simplified dynamic model in ref. [
22]. Note that the inlet HTF temperature of PB is set as a constant during the identification simulation. In this way, the relationship between control and controlled variables is expressed as below:
where
u1 and
u2 are the HTF flowrate and opening of the main valve, respectively; and
E and
P are the electricity generation and steam pressure, respectively.
The non-diagonal elements indicate the coupling effect between the controlled variables. To eliminate the non-diagonal effect, a compensation method is applied. This approach treats the coupling effects as the feedforward signal, and compensates it by choosing a proper decoupling compensation matrix, as shown in the block diagram of the control structure given in
Figure 9. Here, the coupling effect is fully eliminated by setting
D11 and
D22 as 1, and
Having deigned the decoupling matrix, the original coupled controlled issue is degraded to two independent single-input single-output problems, where the PI3 and PI4 can be designed conveniently.
3.3. Feedforward-Feedback Control Scheme of TES
From the CSP overall control strategy diagram, it can be seen that the HTF flowrate of the PB subsystem is responsible for power tracking, while the HTF flowrate of the SF is responsible for regulating SF outlet HTF temperature. The difference between the two flowrates are then delivered into the TES subsystem. Therefore, for the TES subsystem, the HTF flowrate is not controllable but depends on system operation conditions. Consequently, it is treated as a feedforward signal for the controller of TES. In order to ensure the efficiency of the TES subsystem, it must control the MS flowrate to compensate the influence of varying HTF flowrate. In addition, the measurable HTF temperature also serves as a feedforward signal to improve control performance of the TES subsystem. Finally, the block diagram of the feedback–feedforward control scheme of the TES subsystem is shown in
Figure 10.
According to the block diagram, disturbance of HTF flowrate and HTF temperature can be eliminated theoretically based on the invariance principle. By proper identification, the feedforward compensator can be designed as
where
GTES,hx is the transfer function from MS inlet flowrate to MS outlet temperature,
GTES,fl is the transfer function from HTF flowrate to HTF outlet temperature, and
GTES,tem is the transfer function of SF outlet HTF temperature to MS outlet temperature.
FF2 and
FF3 are feedforward compensators for HTF temperature and flowrate, respectively.
PI2 is the feedback controller.
Similar to the feedback controller of SF, PI2 is designed and tuned on the basis of GTES,hx.
5. Conclusions
Different from existing research of CSP-TES systems, investigating the potential of CSP-TES systems to flexibly operate in grid system regulation is the core purpose of this research. To achieve this goal, this work first analyzed the dynamic characteristics of a CSP-TES system; specifically, step response of the HTF flowrate, HTF temperature of PB, opening of main valve and solar irradiation were studied. The results show that CSP and conventional plants both have a multiple variable control problem with large inertial time, but the control problem of the CSP-TES system is more complex because of strong uncertainties and switch nonlinearity, despite its weak coupling. This indicates that the control tasks are not only the well-known disturbance rejection, but also the coordination of different subsystems to the flexible operation of the CSP system.
Afterwards, a decentralized coordinated control strategy was designed to realize the flexible operation of the CSP system. Specifically, a DOB-based feedforward and feedback control strategy was designed for SF, a feedforward-feedback scheme for TES, and decoupling control for PB. These control schemes were then implemented on a system-level analytical dynamic system. The simulation results demonstrate that the proposed control strategy had good disturbance rejection performance, and could realize the flexible operation of the CSP-TES system even under wide-range load variation or strong disturbance, or both.
However, with the transition of conventional economic operation from independent to flexible participation in grid regulation, the role transfer of the CSP-TES system will face many challenges. For instance, how to design advanced control strategies to realize the deep variable working operation of CSP-TES; how to optimize the configuration and fully tape the flexible potential of energy storage devices; and how the semi-controllable CSP system collaborates with uncontrollable renewable energy sources, such as wind power and photovoltaics, remains uncertain. These all require further in-depth study in the future.