scieee AI-readable full text Open interactive document viewer

Model predictive control and optimization of solar thermal energy plants

Sánchez del Pozo Fernández, Adolfo Juan

Abstract

Los objetivos de la tesis son el desarrollo de nuevas estrategias de control predictivo basadas en modelo y optimización de plantas solares parabólicas a gran escala. La tesis se centra en tres aspectos del control de plantas solares: (1) Desenfoque de los colectores, (2) control de plantas solares en presencia de limitaciones de potencia eléctrica y (3) optimización del campo solar. El dise˜no de los controladores propuestos y los resultados obtenidos en la tesis se basan en dos modelos de simulaci´on de campo solar, el campo ACUREX y una planta solar comercial de 50 MW. Es importante mantener la temperatura del fluido t´ermico dentro de los l´ımites de temperatura de seguridad proporcionados por el fabricante y evitar su degradaci´on. Para evitar que la temperatura exceda el l´ımite de degradaci´on, las plantas comerciales disponen de una estrategia de seguridad de desenfoque de colector. Esto se aplica de forma escalonada mediante desenfoque total o parcial, que se realiza modificando el ´angulo del colector. Este mecanismo de desenfoque de colector es reactivo y altamente ineficiente debido a los saltos térmicos causados por la aplicación del desenfoque total o parcial. En esta tesis se presentan nuevos algoritmos de Control Predictivo basado en Modelo para el desenfoque de los colectores de los lazos del campo solar. Se aplicarán Controladores Predictivos Generalizados basados en eventos y Controladores Predictivos basados en Modelos en el espacio de estados para el control del desenfoque del cuarto y tercer colector de los lazos del campo solar. Se presenta un Controlador Predictivo Generalizado basado en eventos para el tracking de potencia eléctrica generada. Este controlador se diseña para los casos en los que la planta se tiene que mover a un modo de operación especial, ”limitación de potencia”. Esta situación se produce cuando la planta recibe una orden del Operador del Sistema de Transporte indicando que se debe disminuir la potencia eléctrica generada. En estos casos, la planta se ve obligada a disminuir su producción eléctrica y mantener la consigna de potencia determinada por el Operador del Sistema de Transporte. En estas situaciones el objetivo es doble: cumplir con el setpoint de potencia indicado por el Operador del Sistema de Transporte y el seguimiento de la temperatura de salida de campo solar. En estas situaciones ser´a muy importante tener en cuenta la estrategia de desenfoque para evitar sobrecalentar el fluido térmico (HTF) (generalmente aceite). Además de los controladores propuestos para el desenfoque y limitación de potencia, una parte de la tesis se centrar´a en la optimización del campo solar. Se centrar´a en obtener una homogeneización del campo solar (balance térmico). Se tendrán en cuenta que factores como la reflectividad pueden ser diferentes en cada uno de los lazos provocando grandes diferencias en las temperaturas de los lazos. Esto puede provocar pérdidas de potencia, diferencias elevadas de temperatura entre los lazos y acciones innecesarias de desenfoque con el correspondiente deterioro de los actuadores, mangueras flexibles y estructuras. Se proponen algoritmos de optimización basados en modelos no lineales para el control de las válvulas de entrada de cada uno de los lazos del campo solar con el fin de homogeneizar la temperatura de salida del campo solar, o lo que es lo mismo, un balance térmico del campo solar. Este equilibrio térmico disminuir´a las grandes diferencias de temperatura entre los diferentes lazos, así como la posible pérdida de energía debido al desenfoque innecesario de parte del campo solar con mayor eficiencia óptica o mayor radiación solar. En una ´ultima parte de la tesis, se aplicar´a el algoritmo de optimización desarrollado para el balance térmico del campo en casos de nubosidad parcial sobre el campo solar. Dado que las plantas comerciales ocupan extensiones de terreno muy grandes, es posible encontrar que partes del campo estén cubiertas por nubes mientras que otras partes no lo estar´an. ´ Esto provocar´a que los lazos con menor radiación solar se encuentren a temperaturas muy inferiores a otras partes del campo solar. El controlador global de temperatura de salida de campo solar podrá estar funcionando correctamente (seguimiento de la temperatura de salida de campo solar), pero para poder cumplir con su objetivo disminuir´a el caudal haciendo que los lazos con más radiación puedan alcanzar el punto de desenfoque, produciendose pérdidas de energía.

Full text

UNIVERSIDAD DE SEVILLA ESCUELA TECNICA SUPERIOR DE INGENIERIA PhD THESIS MODEL PREDICTIVE CONTROL AND OPTIMIZATION OF SOLAR THERMAL ENERGY PLANTS Adolfo Juan S´ anchez del Pozo Fern´ andez 2019 UNIVERSIDAD DE SEVILLA ESCUELA TECNICA SUPERIOR DE INGENIERIA PhD THESIS MODEL PREDICTIVE CONTROL AND OPTIMIZATION OF SOLAR THERMAL ENERGY PLANTS A thesis submitted to the Departamento de Ingenier´ ıa de Sistemas y Autom´ atica, Escuela T´ ecnica Superior de Ingenier´ ıa, in fulfilment of the requirements for the degree of Doctor of Philosophy, Universidad de Sevilla, Espa˜ na Author: Adolfo Juan S´ anchez del Pozo Fern´ andez Supervisors: Prof. Eduardo Fern´ andez Camacho, PhD. Juan Manuel Esca˜ no Gonzalez, PhD. Sevilla, September 2019 A mi familia To my family Acknowledgments First of all I would like to deeply thank my supervisors Eduardo Fern ´ andez Camacho and Juan Manuel Esca ˜ no Gonz ´ alez. Their dedication, trust, support, patience and great experience as well as their contribution of ideas have been vital for the development of this thesis. I want to thank Eduardo Fern ´ andez Camacho for introducing me to the topic of solar energy as well as all the knowledge I have acquired thanks to his vast experience in control engineering, especially in model predictive control and control of solar plants. I want to thank Juan Manuel Esca ˜ no Gonz ´ alez, for having given me valuable knowledge due to his great experience in automation, model predictive control, fuzzy logic control, machine learning techniques as well as his advice and moral support. But most importantly, thank you both for your friendship. Thank you for believing in me, for your continuous advice and for supporting me in the good and bad moments. I am proud to have had you as supervisors. It would not have been possible to carry out this thesis without my family. To my parents, Adolfo and Esther, and my sisters, Esther and Mar ´ ıa, to whom I owe everything, education, studies, unconditional love, infinite patience and understanding and help in all senses and moments of my life. Thank you, I would be nothing without you. After several years of effort, there have been many people who have contributed, in one way or another, to the development of this work. I am very grateful to all the members of the predictive control group of the Department of Systems and Automation Engineering, University of Seville. During these years I have been able to attend group meetings, seminars, talks and presentations where I have learned a lot from the work done by its members. Thank you all for your support, advice and encouragement. To thank the OCONTSOLAR project for financing this thesis as well as its members. Part of this thesis has been developed in the Nimbus Research Center, Cork Institute of Technology, Ireland. Thanks to all the members, colleagues and friends of the Nimbus Research Center, for having enjoyed very good times during my stay in Ireland. Thanks in particular to my colleagues and friends of the Advanced Control Group of Nimbus Centre, Niel Canty, Kritchai Witheephanich, Jean Michel Rubillon, Samira Roshany-Yamchi and Louis Cronin with whom I had the privilege of learning and working on various research and development projects. I want to thank my friends from Jerez, Seville, Madrid and Tarragona, for their great friendship, unconditional support, understanding and for their infinite patience since during these years of thesis I have stopped giving signs of life in certain occasions for some or other reasons and they have always been there. Thank you all for being part of my life. I would like to especially thank my colleague and friend Antonio Javier Gallego Len, with whom I started working in the control of solar plants, under the direction of Eduardo Fern ´ andez Camacho, in several solar energy research projects. Thank you Antonio, since your great experience in this field has been vital to be able to start and finish this thesis. Thank you for your advice, contributions, ideas and support at all times, good and bad. Even in your worst moments you always had time to give me advice or to just have a chat with a coffee. Thank you for your friendship. Try not to become a man of success, but rather try to become a man of value. Albert Einstein Agradecimientos En primer lugar me gustar ´ ıa agradecer profundamente a mis directores de tesis Eduardo Fern ´ andez Camacho y Juan Manuel Esca ˜ no Gonz ´ alez. Su dedicaci ´ on, confianza, apoyo, paciencia y gran experiencia as ´ ı como su aportaci ´ on de ideas han sido vitales para el desarrollo de esta tesis. Agradecer a Eduardo Fern ´ andez Camacho por introducirme en el campo de la energ ´ ıa solar as ´ ı como todos los conocimientos que he adquirido gracias a sus grandes conocimientos y experiencia en ingenier ´ ıa de control, especialmente en control predictivo basado en modelo y en control de plantas solares. Agradecer a Juan Manuel Esca ˜ no, por haberme proporcionado valiosos conocimientos por su gran experiencia en control autom ´ atico, predictivo, borroso, t ´ ecnicas de aprendizaje y machine learning as ´ ı como sus consejos y apoyo moral continuo. Pero lo m ´ as importante, gracias a ambos por vuestra amistad. Gracias por creer en mi y apoyarme en los buenos y malos momentos. Es un orgullo haberos tenido como directores. Esta tesis nunca habr ´ ıa sido posible sin mi familia. A mis padres, Adolfo y Esther, y a mis hermanas, Esther y Mar ´ ıa, a los que les debo todo, educaci ´ on, estudios, amor incondicional, infinita paciencia y compresi ´ on y ayuda en todos los sentidos y momentos de mi vida. Gracias, no ser ´ ıa nada sin vosotros. Tras varios a ˜ nos de esfuerzo, han sido muchas las personas que han contribuido, de una forma u otra, en el desarrollo de este trabajo. Estoy muy agradecido a todos los compa ˜ neros y amigos del Departamento de Ingenier ´ ıa de Sistemas y Autom ´ atica de la Escuela T ´ ecnica Superior de Ingenier ´ ıa y del grupo de control predictivo. Durante estos a ˜ nos he podido adquirir valiosos conocimientos en seminarios, charlas, reuniones y presentaciones de los trabajos realizados por todos ellos. Gracias a todos por vuestro apoyo, consejo y animo. Agradecer al proyecto OCONTSOLAR por financiar esta tesis as´ ı como a sus miembros. Parte de esta tesis ha sido desarrollada en el Nimbus Research Centre, Cork Institute of Technology, Irlanda. Gracias a todos los miembros, compa ˜ neros y amigos del Nimbus Centre, por haber disfrutado de muy buenos momentos durante mi estancia en Irlanda. Gracias en particular a mis compa ˜ neros y amigos del grupo de Control Avanzado del Nimbus, Niel Canty, Kritchai Witheephanich, Jean Michel Rubillon, Samira Roshany- Yamchi y Louis Cronin con los que tuve el privilegio de aprender y trabajar en diversos proyectos de investigaci´ on y desarrollo. Agradecer a mis amigos de Jerez, Sevilla, Madrid y Tarragona, por su gran amistad, apoyo incondicional, comprensi ´ on y su paciencia ya que durante estos a ˜ nos de tesis he dejado de dar se ˜ nales de vida en ciertas ocasiones por unos u otros motivos y siempre han estado ah ´ ı. Gracias a todos por vuestra amistad y por ser parte de mi vida. He de agradecer muy especialmente a mi compa ˜ nero de trabajo y amigo Antonio Javier Gallego Len, con quien empec ´ e a trabajar en el control de plantas solares, bajo la direcci ´ on de Eduardo Fern ´ andez Camacho, en varios proyectos de investigaci ´ on en energ ´ ıa solar. Gracias Antonio, porque tu gran experiencia en este campo ha sido vital para poder empezar y terminar esta tesis. Gracias por tus consejos, contribuciones, aportaciones de ideas y apoyo en todo momento, buenos y malos. Incluso en tus peores momentos siempre tuviste tiempo para darme consejos o simplemente charlar con un caf´ e. Gracias por tu amistad. Intenta no convertirte en un hombre de ´ exito, sino m´ as bien intenta convertirte en un hombre de valor. Albert Einstein Abstract Model Predictive Control and Optimization of Solar Thermal Energy Plants Adolfo Juan S´ anchez del Pozo Fern´ andez The goals of the thesis are the development of new model predictive control strategies and optimization of large-scale parabolic solar plants. The thesis focuses on three aspects of the control of solar plants: (1) defocusing of the collectors, (2) control of solar plants in the presence of limitations of electric power and (3) optimization of the solar field. The design of the proposed controllers and the results obtained in the thesis are based on two solar field simulation models, the ACUREX field and a 50 MW commercial solar plant. It is important to keep the oil temperature within the safe temperature limits provided by the manufacturer and to avoid its degradation. To prevent the temperature from exceeding the degradation limit commercial plants have a safety collector defocus strategy. This is applied in a staggered manner as total or partial defocusing, which is done by modifying the angle of the collector. This mechanism of defocusing the collector is reactive and highly inefficient due to the thermal jumps caused by the application of full or partial defocusing. In this thesis new Model Predictive Control (MPC) algorithms for the defocusing of the solar field loops collectors are presented. Event-based Generalized Predictive Control and State Space Model Predictive Control strategies will be applied for the defocusing of the fourth and third collector of the solar field loops. An Event-based Generalized Predictive Control strategy for the generated electrical power tracking is presented. This controller is designed for cases in which the plant has to move to a special operation mode, ”power limitation”. This situation arises when the plant receives an order from the Transmission System Operator (TSO) indicating that the power generation must be reduced. In these cases, the plant is forced to decrease its electric production and maintain the power set-point determined by the TSO. In these situations the objective is double: fulfilling the TSO power set-point and temperature tracking. In these situations it will be very important to take into account the defocus strategy to avoid overheating the Heat Transfer Fluid (HTF) (usually oil). iv CONTENTS 7 Thermal balance of a solar trough field under cloud coverage 141 7.1 Introduction ......................... 141 7.2 Nonlinear model based optimization for valve control . . . 143 7.3 Simulated plant scenario .................. 146 7.4 Results ............................ 149 7.5 Conclusion ......................... 153 8 Conclusion 154 References 159 List of Figures 1.1 Solar energy power comparative .............. 1 1.2 Global primary energy consumption. Source: Our World in Data ............................. 2 1.3 Solar energy reaching earth ................. 3 1.4 A wall painting from the Uffizi Gallery, Stanzino delle Matematiche, in Florence, Italy, shows the Greek mathematician Archimedes’ mirror burning Roman military ships. Painted in 1600 by Gieulio Parigi .............. 4 1.5 Solar powered motor at the Cawston’s Ostrich Farm. Source: DeGolyer Library ...................... 5 1.6 PS20 PS10 solar tower plants in Sanlucar la Mayor, Seville, Spain. Source: Koza1983, Wikimedia ........... 8 1.7 Crescent Dunes Solar Energy Facility. Source: Amble, Wikimedia .......................... 8 1.8 Ivanpah Solar Power Facility (all towers online). Source: Sbharris, Wikimedia .................... 9 1.9 Linear Fresnel Collector solar cooling plant located at the Engineering School (ESI) of Seville ............ 10 1.10 Novatec Fresnel Power Station Puerto Errado 2. Source: Novatec Solar, Wikimedia ................. 11 1.11 Parabolic Solar Collector at Mojave Desert, Mojave Solar Project. Source: Z22, Wikimedia .............. 12 1.12 Parabolic collector schematic ................ 13 1.13 Top view of the new TCP-100 parabolic trough facilities in Almer´ ıa (PSA-CIEMAT) (Spain) .............. 15 v vi LIST OF FIGURES 1.14 Solar Energy Generating Systems solar power plants IIIVII at Mojave Desert, California. Source: Alan Radecki Akradecki, Wikimedia ................... 18 1.15 Andasol Solar Power Station, Guadix (Spain). Source: BSMPS, Wikimedia ..................... 18 2.1 Acurex solar field schematic ................ 30 2.2 Parabolic trough plant general schematic .......... 32 2.3 EUROTrough Prototype at the Plataforma Solar de Almer ´ ıa. Source: European Commission technical report ...... 33 3.1 Solar field step response of linear models. Output increments when unitary steps are applied to the input ..... 44 3.2 GS-GPC + FeedForward control scheme .......... 45 3.3 GS-GPC + FF tracking results (medium irradiance) . . . . 46 3.4 GS-GPC + FF. Flow-Rate at maximum due to high irradiance leading to temperature limit violation ......... 47 3.5 Efficiency - Defocus angle curve .............. 48 3.6 Collector step response of linear models (low and lowmedium flow-rates) ..................... 49 3.7 Collector step response of linear models (medium-high and high flow-rates) ....................... 50 3.8 Flow GS-GPC + FF + C4 EGS-GPC loop temperature tracking results. Top plot: field fluid temperatures. Bottom plot: Irradiance and flow-rate ................... 52 3.9 Flow GS-GPC + FF + C4 EGS-GPC loop temperature tracking results. Top plot: loop fluid temperatures. Bottom plot: Fourth Collector defocus angle control actions ....... 53 3.10 Flow GS-GPC + FF + C4 EGS-GPC control scheme . . . 54 3.11 Power cycle step response of linear models. Output increments when unitary steps are applied to the input ..... 55 3.12 Power Limitation Event block state graph ......... 56 LIST OF FIGURES vii 3.13 Flow GS-GPC + FF + C4 EGS-GPC outlet temperature limit violation. Fully defocused fourth collector. Top plot: field fluid temperatures. Bottom plot: Irradiance, flow-rate and electric power ........................ 58 3.14 Flow GS-GPC + FF + C4 EGS-GPC outlet temperature limit violation. Fully defocused fourth collector. Top plot: loop fluid temperatures. Bottom plot: Fourth collector defocus angle ............................. 59 3.15 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC control scheme ....................... 60 3.16 Flow GS-GPC + FF + C4 EGS-GPC + Power EGS-GPC, 40 MW TSO limitation on at 11:40am and off at 15:45pm (60 min ramp). Field HTF temperature, irradiance, flow-rate and electric power ...................... 61 3.17 Flow GS-GPC + FF + C4 EGS-GPC + Power EGS-GPC, 40 MW TSO limitation at 11:40am and TSO limitation off at 15:45pm (60 min ramp). Top plot: loop fluid temperatures. Bottom plot: Fourth Collector defocus angle control actions 62 3.18 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 30 MW TSO limitation at 11:40am and TSO limitation off at 15:45pm (60 min ramp). Top plot: field fluid temperatures. Bottom plot: Irradiance, flow-rate and electric power . . . 63 3.19 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 30 MW TSO limitation at 11:40am and TSO limitation off at 15:45pm (60 min ramp). Top plot: loop fluid temperatures. Bottom plot: Fourth collector defocus angle ........ 64 3.20 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 30 MW TSO limitation at 11:40am and TSO limitation off at 15:45pm (60 min ramp). Top plot: loop fluid temperatures. Bottom plot: Third collector defocus angle ......... 65 3.21 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 40 MW TSO limitation at 11:40am (60 min ramp) and TSO limitation off at 15:45pm (30 min ramp). Top plot: field fluid temperatures. Bottom plot: Irradiance, flow-rate and electric power ........................ 66 viii LIST OF FIGURES 3.22 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 40 MW TSO limitation at 11:40am (60 min ramp) and TSO limitation off at 15:45pm (30 min ramp). Top plot: loop fluid temperatures. Bottom plot: Fourth collector defocus angle control actions .................... 67 3.23 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 40 MW TSO limitation at 11:40am (60 min ramp) and TSO limitation off at 15:45pm (30 min ramp). Top plot: loop fluid temperatures. Bottom plot: Third collector defocus angle control actions .................... 68 3.24 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 30 MW TSO limitation at 11:40am (60 min ramp) and TSO limitation off at 15:45pm (30 min ramp). Top plot: field fluid temperatures. Bottom plot: Irradiance, flow-rate and electric power ........................ 69 3.25 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 30 MW TSO limitation at 11:40am (60 min ramp) and TSO limitation off at 15:45pm (30 min ramp). Top plot: loop fluid temperatures. Bottom plot: Fourth collector defocus angle control actions .................... 70 3.26 Flow GS-GPC + FF + C4/C3 EGS-GPC + Power EGS-GPC, 30 MW TSO limitation at 11:40am (60 min ramp) and TSO limitation off at 15:45pm (30 min ramp). Top plot: loop fluid temperatures. Bottom plot: Third collector defocus angle control actions .................... 71 4.1 Collector fluid temperatures estimation ........... 82 4.2 Collector efficiency estimation using the concentrated parameter model ........................ 83 4.3 Flow/Power EGS-GPC + C4/C3 MPC control scheme . . . 84 4.4 Efficiency - Defocus angle curve with modeling errors . . 85 4.5 GS-GPC vs MPC results. Real Defocus Curve 1. 40 MW limitation received (60 min ramp). Top plot: Fluid temperature. Bottom plot: Collector defocus angle control action . 86 LIST OF FIGURES ix 4.6 GS-GPC vs MPC results. Real Defocus Curve 2. 40 MW limitation received (60 min ramp). Top plot: Fluid temperature. Bottom plot: Collector defocus angle control actions . 87 4.7 GS-GPC vs MPC results. Real Defocus Curve 1. 30 MW limitation received (60 min ramp). Top plot: Fluid temperature. Bottom plot: Collector defocus angle control actions . 88 4.8 GS-GPC vs MPC results. Real Defocus Curve 2. 30 MW limitation received (60 min ramp). Top plot: Fluid temperature. Bottom plot: Collector defocus angle control actions . 89 4.9 GS-GPC vs MPC results. Real Defocus Curve 2. GS-GPC control exciting high frequency dynamics. Top plot: Fluid temperature. Bottom plot: Collector defocus angle control actions ............................ 90 5.1 Acurex GS-GPC + FeedForward control scheme ...... 94 5.2 Acurex Solar Field step response of linear models ..... 95 5.3 GS-GPC+FeedForward outlet temperature tracking results 96 5.4 GS-GPC+Feedforward outlet temperature tracking and Loops outlet temperatures. Loops reflectivities = [0.65, 0.67, 0.68, 0.7, 0.715, 0.73, 0.75, 0.77, 0.79, 0.8] ........... 97 5.5 Acurex field loops, valves and pump schematic ...... 98 5.6 Classification and Regression Tree node splitting ..... 99 5.7 CART based temperature estimates are compare to the nonlinear model. Temperature profiles at the output of the collectors (top plot) and estimation errors (bottom plot) . . 101 5.8 August case. Loop efficiency estimation .......... 103 5.9 December case. Loop efficiency estimation ........ 103 5.10 Proposed control scheme .................. 106 5.11 Case 1. GS-GPC+FeedForward results (no valve control, all at 100%) ........................... 107 5.12 Case 1. GS-GPC+FeedForward results with reflectivity estimation and valve control. Loops temperatures (top plot) and Valves apertures (bottom plot) ............. 108 5.13 Case 2. GS-GPC+FeedForward results (no valve control, all at 100%). Defocus control is active in some loops . . . . . 109 x LIST OF FIGURES 5.14 Case 2. GS-GPC+FeedForward results with valve control. Loops temperatures (top plot) and Valves apertures (bottom plot) ............................. 109 5.15 Case 3. GS-GPC+FeedForward results (no valve control, all at 100%). Defocus control is active in some loops . . . . . 110 5.16 Case 3. GS-GPC+FeedForward results with valve control. Loops temperatures (top plot) and Valves apertures (bottom plot) ............................. 111 5.17 Case 2. Flow-rate comparison. GS-GPC+FeedForward results with valve control .................. 111 5.18 Case 3. Flow-rate comparison. GS-GPC+FeedForward results with valve control .................. 112 6.1 Loops temperatures (top) and defocus angles (bottom) when there is no valve control under stable radiation ....... 120 6.2 Field outlet temperature (top) and field flow-rate/irradiance (bottom) (no valve control applied) ............. 121 6.3 Optimization computation times depending on the number of the decision variables .................. 126 6.4 Example 1: Clustering result (10 groups) for 90 loops when loops are not defocusing .................. 129 6.5 Example 2: Clustering result (10 groups) for 90 loops when there are loops being defocused ............... 130 6.6 Loops temperatures (top) and defocus angles (bottom) results of the valve control strategy for thermal balance with 10 clusters for a stable radiation .............. 131 6.7 Loops temperatures (top) and defocuse angles (bottom) results of the valve control strategy for thermal balance with 20 clusters for stable radiation ............... 132 6.8 Field outlet temperature (top) and flow-rate (bottom) comparison with stable radiation ................ 133 6.9 Valve control actions for the 10 (top) and 20 (bottom) clusters simulations with stable radiation ............ 134 6.10 Loops traveled defocus angles comparison with stable radiation ............................. 135 LIST OF FIGURES xi 6.11 Loops temperatures (top) and defocus angles (bottom) when there is no valve control under transient radiation . . . . . 136 6.12 Loops temperatures for 10 variables (left) and 20 variables (right) results of the valve control strategy for thermal balance with transient radiation ................ 137 6.13 Field outlet temperature (top) and flow-rate (bottom) comparison with transient radiation ............... 138 6.14 Valve control actions for 10 clusters simulations with transient radiation ........................ 138 6.15 Valve control actions for 20 clusters simulations with transient radiation ........................ 139 6.16 Loops traveled defocus angles comparison with transient radiation ........................... 139 7.1 Loops rows and numbering for each solar field (upper and lower) ............................ 144 7.2 Upper solar field loops block division ........... 144 7.3 Lower solar field loops block division ........... 145 7.4 Solar field partial coverage by passing clouds (collectors 2 and 3 being covered) .................... 147 7.5 Solar field partial coverage by passing clouds (collectors 1 and 4 being covered) .................... 148 7.6 Loops temperatures. Top plot: Loop blocks outlet temperatures, no valve control. Bottom plot: Loop blocks outlet temperatures, NLVC .................... 149 7.7 Field state. Top plot: Field outlet temperature. Bottom plot: Field flow-rate ........................ 150 7.8 Loops defocus angles comparison. Top plot: Loop blocks defocus actions, no valve control. Bottom plot: Loop blocks defocus actions, NLVC ................... 150 7.9 NLVC valves control actions. Top plot: Loop blocks B1-B5 (upper field) valves apertures. Bottom plot: Loop blocks B6-B10 (lower field) valves apertures ........... 151 List of Tables 2.1 EuroTrough ET150 parameters. .............. 33 3.1 Simulated Loop Efficiencies. ................ 70 4.1 ITAE comparison. ...................... 90 4.2 ISE comparison. ....................... 91 5.1 Simulated scenarios. Reflectivities. ............. 108 5.2 Mean Power Production (MW) (12am - 18pm). ...... 113 6.1 Reflectivity ranges and Loops ................ 119 6.2 Clusters and Loop members ................ 129 6.3 Number of defocus and valve actions, case 1 (stable radiation) 135 6.4 Number of defocus and valve actions, case 2 (transient radiation) ............................ 140 7.1 Defocus Reduction (10am - 16pm). ............. 151 7.2 Travelled Angles Reduction (10am - 16pm). ........ 152 7.3 Power (MW) Improvement (10am - 16pm). ........ 152 xii Acronyms BESS Battery Energy Storage System CARIMA Controlled Auto Regressive Integrated Moving Average CART Classification and Regression Tree CIEMAT Centro de Investigaciones Energ´ eticas, Medioambientales y Tecnol´ ogicas CPM Concentrated (Lumped) Parameter Model CSP Concentrated Solar Plant DNI Direct Normal Irradiance DPM Distributed Parameter Model EGS-GPC Event Based Gain Scheduling Generalized Predictive Control EKF Extended Kalman Filter FDMC Filter Dynamic Matrix Control FF Feed Forward GA Genetic Algorithm GPC Generalized Predictive Control GS Gain Scheduling GS-GPC Gain Scheduling Generalized Predictive Control HCE Heat Collector Element HTF Heat Transfer Fluid ISE Integral Square Error ITAE Integral of Time Absolute Error LFC Linear Fresnel Collector LFR Linear Fresnel Reflector MLD Mixed Logical Dynamical (MLD) model MPC Model-based Predictive Control NLVC Nonlinear Valve Control NMPC Nonlinear Model Predictive Control xiii 4Introduction Syracuse (III century B.C.). Fig. 1.4 shows an illustration of Archimedes mirrors. However, this is part of a legend and it is not known how much of this is true in the writings. Centuries later, mirrors were used for metal smelting and burning of trees, among others [4]. Figure 1.4: A wall painting from the Uffizi Gallery, Stanzino delle Matematiche, in Florence, Italy, shows the Greek mathematician Archimedes’ mirror burning Roman military ships. Painted in 1600 by Gieulio Parigi In 1515, Leonardo started one of his many projects, although it would never be finished. His project was to build a concave mirrors based solar concentrator for industrial steam and heat production. In 1615, the French engineer Salomon de Caux engineered an engine, whose source was solar energy through glass lenses, which through the expansion of air could pump water [5]. It was the first recorded mechanical application of the Sun’s energy. It must be highlighted the work of Horace-B ´ en ´ edict de Saussure in 1767. He wanted to see what temperature could be reached in a glass box thanks to the greenhouse effect. He managed to reach 109 C and without knowing it, he created the so-called solar collector. Antoine-Laurent de Lavoisier, French chemist, built in 1792 a solar oven consisting of two powerful lenses that concentrated the solar radiation in a focus point allowing him to reach high temperatures to melt metals. The second half of the 19th century was especially interesting in the development of new techniques for capturing solar energy. Several were the 1.1. Solar Energy: History 5 pioneers during this stage who developed systems to heat domestic water, motors powered by the use of solar energy, as well as the first thermosolar plants. Augustin Mouchot’s conviction was that coal, which fueled the Industrial Revolution, was not going to last forever and led him to notice the energy that came from the Sun. In 1860 he built a solar kitchen and later he designed a water-filled boiler enclosed in glass, which would be exposed to the heat of the sun until the water boiled. A small steam engine motor would use the steam produced by the sun boiler. By August 1866, Mouchot had developed the first parabolic trough solar collector [6]. John Ericsson, Swedish engineer and inventor, designed in 1870 a solar collector attached to a 373 W hot air motor [7]. Charles Wilson, Swedish engineer, designed and built a desalination plant using solar energy in the desert of Atacama (Chile). This plant produced around 20,000 litres per day of fresh water [8]. In 1891, Clarence Kemp designed and patented the first water heater with solar energy thanks to the use of a collector. Registered the patent as ’Climax’, promoting the commercial development of this type of energy [4]. Figure 1.5: Solar powered motor at the Cawston’s Ostrich Farm. Source: DeGolyer Library [9] In 1901, Aubrey Eneas built the first successful experiment of a solar 6Introduction powered motor for commercial used. He installed it at the Cawston’s Ostrich Farm in South Pasadena, Fig. 1.5. The parabolic dish mirrors concentrated solar energy on a boiler to produce steam to power an engine to pump underground water to the ostrich farm and founded the first solar company, The Solar Motor Co. [10]. In 1909 William Bailey introduced a substantial improvement in the development of solar water heater. Bailey proposed and patented a solar water heater by dividing it into two stages. These stages would separate the water storage system from the solar system that would heat it by using a circulation and heat conservation system. Its design was the basis of solar heaters that are currently used [4]. One of the most ambitious was the American Frank Shuman. Shuman founded his company Sun Power Co. in 1911. In 1912-1913, he built the first parabolic trough solar plant in Maadi, Egypt. The plant, called ”Solar Engine One”, was composed of five 60 m length parabolic reflectors, north-south oriented, with a mechanical tracker mechanism to track the sun along the day to concentrate the sun rays in order to generate steam to pump 22,700 l/m approximately, to irrigate the desert near the Nile River. But his ambition had no limits, so he struggled to get a solar plant capable of generating all the energy consumed by the whole world. Shuman calculated that a 52,600 square kilometers parabolic trough plant, built in the Sahara, could produce the equivalent energy to all the fuel burnt in the world in 1909 [11]. Before the beginning of the war, Shuman reached agreements with the British and German governments for his proposal of solar energy in the tropics. The German Reichstag agreed to give an advance of 200.000 marks to start the project [12], but the start of World War I changed everything. Shuman’s project was automatically paralyzed and all workers at the Maadi solar plant returned to Germany to fight on the German side. The installations at Maadi were destroyed. Schuman died during the course of the war in 1918. At the end of the war, with the defeat of Germany, the loss of all African colonies and the discovery of oil and gas reserves, which were sold at very low prices, the interest in solar energy was lost and Shuman’s project fell into oblivion [13]. The oil crisis during the second half of the 70s caused an interest in renewable energy unprecedented, although this interest was not due to the 1.2. State of the Art 7 renewable nature but to economic factors because oil prices soared. However, when oil prices decreased again, interest in renewables was lost. However, this interest has resurged, although not due to economic factors caused by the price of oil, but due to environmental awareness. Due to global warming and with the objective of reducing harmful pollutant gases emissions from conventional fossil power plants such as SO2 , CO , NOx , HC and CO2 , the interest in renewable energies has, once again, resurged [14 – 16]. Currently renewable energies with the greatest impact on society are solar, wind and hydraulic. Solar energy being, by far, the most abundant. 1.2 State of the Art The first thermosolar plants for experimentation purposes were plants of reduced size, e.g., the old solar trough plant ACUREX and the new TCP-100 field in the Plataforma Solar de Almer ´ ıa (PSA) [17] or the Juelich solar tower research facility [18]. However, nowadays, commercial solar thermal power plants cover vast areas of land. There are several current technologies to produce electricity from the direct solar radiation. Among them, the following stand out: 1. Solar Tower Power Plants 2. Linear Fresnel Reflector Power Plants 3. Solar Parabolic Trough Power Plants 1.2.1 Tower and Fresnel Power Plants The solar thermal tower power plants are point-focus systems that concentrate the solar energy in a central receiver located at the top of a tower. A field of mirrors or ”heliostats” are in charge of reflecting the direct solar radiation in said receiver. This receiver will convert the concentrated heat into steam which will then be converted into electrical energy by means of a turbine. In general, a solar tower power plant is composed of the aforementioned tower and the field of heliostats that may surround it partially or completely depending on the number of receivers of the tower. For example, the PS10 8Introduction (10 MW) and PS20 (20MW) towers, owned by Atlantica Yield [19], are made up of partial fields of heliostats as can be seen in Fig. 1.6, while in Fig. 1.7 shows the image of the Crescend solar plant in which it can be seen how the heliostat field completely surrounds the central tower. The heliostats have a solar tracking system for concentrating the sun’s rays throughout the day [20]. This type of plants can also include thermal storage systems through the use of molten salts, generally nitrates, for power generation in periods without solar radiation. Figure 1.6: PS20 PS10 solar tower plants in Sanlucar la Mayor, Seville, Spain. Source: Koza1983, Wikimedia [21] Figure 1.7: Crescent Dunes Solar Energy Facility. Source: Amble, Wikimedia [22] 1.2. State of the Art 9 Solar tower plants are one of the most promising technologies for the production of electrical energy from solar energy. One of the aspects that continue to be part of the core of the research in this type of plants is the optimal design of the heliostat field. Its complexity is due to the annual performance of the field since it depends strongly on the location, the annual radiation periods and the positions between the heliostats that will cause shadows on other neighboring heliostats thus decreasing the solar energy collection [23]. The solar tower plants have maintained a constant development and have gone from experimental plants to plants of few MW and in the last years to hundreds of MW of production [24]. Among the tower plants that are currently in production, it can be found the already mentioned PS10 (10 MW) and PS20 (20 MW) in Seville (Spain) [25,26], the Supcon Solar Project (50MW) in Delingha (China) [27], Crescent Dunes Solar Energy Facility in Nevada (United States) with a power production of 100 MW and an heliostat field covering 6,474,970 m2 approximately [28], and the Ivanpah Solar Electric Generating System which is a set of three plants with a total gross production of 392 MW (126 MW, 133 MW and 133 MW) over 2,600,000 m2 [29], see Fig. 1.8. Figure 1.8: Ivanpah Solar Power Facility (all towers online). Source: Sbharris, Wikimedia [30] On the other hand, the Fresnel linear collectors or reflectors (LFC, LFR) differs both in the arrangement of the field and in the point and form of solar 10 Introduction energy collection. The mirrors are arranged horizontally and linearly. A linear receiver located over the mirrors is responsible for receiving the energy and transferring it to the internal fluid. Unlike the tower and parabolic plants, the horizontal arrangement implies a lower relation between the aperture and optical efficiency. Although it is important that these plants occupy relatively little space, there must be a compensation for the loss of optical efficiency with a decrease in the cost per m2 of aperture compared to parabolic solar plants. One of the advantages of LFC or LFR technology is the reduced cost due to the use of simple structures and flat mirrors whose manufacture is cheaper than the curved mirrors of the PTC. The use of flat mirrors makes the design of these plants differ, in general, in size, arrangement and dimensions of the mirrors [31,32]. Companies such as SUNCINM or SOLTIGUA [33,34], among others, compete in the development and installation of LFC plants. Typically, the tower and parabolic trough plants have been the technologies used for the generation of electric power from solar energy in large scale. Fresnel plants began their development as solar thermal plants for low and medium scale such as Puerto Errado 1 Thermosolar Power Plant [35] in Murcia (Spain) with a power production of 1.4 MW. Given the small size of these plants it is interesting to use them on the roofs of the buildings for various services. As an example it can be highlighted the experimental Linear Fresnel Collector solar cooling plant located at the Engineering School (ESI) of Seville [36], see Fig. 1.9. Figure 1.9: Linear Fresnel Collector solar cooling plant located at the Engineering School (ESI) of Seville 1.2. State of the Art 11 However, they are also beginning their large-scale development. Proof of these are the 30 MW Puerto Errado 2 Thermosolar Power Plant [37], covering 302,000 m2 of land, sited in Murcia (Spain), see Fig. 1.10, and Dhursar Fresnel plant in Rajasthan (India) with a production of 125 MW and 3,000,000 m2of land area [38]. Figure 1.10: Novatec Fresnel Power Station Puerto Errado 2. Source: Novatec Solar, Wikimedia [39] As an additional comment, it can be found another type of plants that is quickly spreading, photovoltaic plants (PV). Although these types of plants also use direct solar radiation (not exclusively), they are not considered solar concentration systems since the solar cells convert sunlight directly into electricity, therefore, they are outside the scope of this document. However, it should be noted that these types of plants can be used at any scale, being possible to find them both for large-scale use for the generation of large quantities of electrical energy as well as domestic and are becoming in a strong commitment to the integration of renewable sources in power distribution networks [40]. This thesis is mainly focused on cylinder-parabolic technology, although some of the concepts that are exposed are valid also for other of the described technologies. 12 Introduction 1.2.2 Solar Parabolic Trough Power Plants Parabolic Trough Power Plants (PTPP) also named as Parabolic Trough Collector (PTC) Concentrated Solar Plant (CSP) are formed by parallel rows of parabolic mirrors whose focal line is used to set a special receiving pipe through which a synthetic fluid circulates. Each of these rows are commonly called solar field loops and current commercial plants contain from several tens to several hundred loops. Figure 1.11: Parabolic Solar Collector at Mojave Desert, Mojave Solar Project. Source: Z22, Wikimedia [41] Each of the loops that make up the solar field is composed of independent modules called Solar Collector Assemblies (SCA). In general, each of these SCAs are made up of the following elements: 1. Metal structure. 2. Parabolic reflectors (mirrors). 3. Absorber pipe (receiver tube). 4. Flexible hoses. 5. Tracking system, sensors and control. 1.2. State of the Art 13 Figure 1.12: Parabolic collector schematic The mirrors concentrate the direct solar energy in the focus line, where the receiver tube is, and the heat is transferred to the fluid that circulates through the pipe. The heated fluid is transported to the electric generation phase. By means of a heat exchanger the heat of the fluid is used to obtain steam and by means of a turbine electrical energy is generated. Although most generation systems in parabolic trough plants use heat exchangers for steam generation, there are other systems in which the steam is produced directly in the solar field [42]. The receiver tube or heat collection element (HCE), installed along the focal line of the loops, is composed of a thin metal pipe surrounded by a glass cover. In the space left between the metal pipe, where the HTF will flow, and the glass envelope, a vacuum is created to reduce thermal losses to the environment due to convection losses. The solar tracking system is a fundamental element in PTC solar plants. It is important to know the angle between the solar rays and the surface or aperture of the collector. If the collector is not directly pointing to the sun, some of the energy is lost since it is not reflected on the focal line, which is where the receiver pipe with the thermal fluid is located [43]. The thermal fluid or most commonly known as Heat Transfer Fluid (HTF) is, generally, a synthetic thermal oil that supports high temperatures before its degradation. The vast majority of commercial parabolic cylinder plants currently work in temperatures in the range of 300 to 400 ◦C and use this type of oil. The nominal operating point of the plants that use synthetic 20 Introduction produced by the conventional electric generation industry. The solar thermal renewable energy industry is facing many problems in the massive production of electric power. One of the most significant issues is the cost per kWh or MWh. The costs of producing electricity in solar thermal power plants are not yet competitive in comparison with conventional fossil fuel plants. In fact, one of the great challenges of the century identified by the National Academy and the European Commission is to make solar energy economical and competitive. [102,103]. Another aspect that commercial solar plants have to face is the temporary availability of the radiation resource [104]. This means that the production of the solar plants may not be as constant as in conventional fossil fuels plants with the corresponding economic losses and stability in the electricity grid. Therefore, it is important to emphasize the great effort that is being made both at the research and industrial level in the development of improved thermal energy storage units that will contribute to a more stable production even when there is no solar radiation. The application of advanced control and optimization algorithms can play an important role in improving the overall efficiency of solar energy thermal plants and thus improving the penetration of these kind of plants into the global market [1]. If in addition, the research and development of new and improved thermal energy storage systems, as well as the optimization of their loading and unloading by means of advanced control techniques, could considerably increase the number of hours of electrical production, which would imply a reduction in the production costs per MWh making the thermosolar technology [105] more competitive and attractive. Actually, both the research groups and the solar industry focus their efforts on solving the following problems, among others: • Improve the efficiency of solar collectors. Whether thermal or photovoltaic technology. • Improve the efficiency of the conversion cycles of the radiation captured in both thermal and photovoltaic. • The development of improved and new energy storage systems. In this case, in solar thermal is Thermal Energy Storage (TES), and in photovoltaic it is Battery Energy Storage System (BESS). 1.2. State of the Art 21 • The optimization of the operation of the electric generation processes in solar plants. •The reduction of design, manufacturing and installation costs. In PTC CSPs, The efficiency of the solar collector is related both to the tracking system and to the construction materials of the parabolic mirrors, which have a form factor that indicates the level of perfection with which it was built. This factor is usually found in a value close to unity (it is taken as efficiency) but never equal to one due to imperfections, so that part of the efficiency in the solar rays collection is lost due to this imperfection in the mirror manufacturing. On the other hand, the reduction of installation costs is closely related to the use of wider aperture trough collectors [90]. The optimization of the performance of the generation processes can be obtained partially through the use of optimal control strategies of the solar field. It is said partially since the solar plant is composed of two main parts, as already mentioned, the solar phase and the steam phase. In the steam phase, the performance improvement refers to the optimal use of the different elements that make up this system, such as the turbine, the different heat exchangers, condensation phase, expansion stage, steam overheating, etc. Classical control systems such as PIDs are not optimal strategies for maximizing the performance of the solar field due to its simplicity and reactivity in the control. Solar plants contain highly nonlinear dynamics and variable time delays that depend on the HTF flow. One of the main problems in the control of the solar fields is the appearance of antiresonant modes that are problematic when it comes to obtaining fast and stable behavior of the solar field [106,107]. It is also essential to bear in mind that current commercial solar trough cover vast extensions of land. In the case of the two solar trough plants of Mojave cover 700 hectares and they are composed of 282 loops each [96]. The solar trough plant SOLANA is even larger. It is composed of 808 loops covering 780 hectares [45]. These large-scale solar plants show up new challenges for the application of advanced control strategies. New advanced control techniques have to be devised and developed to address these issues. In order to develop new advanced control strategies for large scale solar plants, The Advanced Grant Optimal Control of Thermal Solar Energy 22 Introduction Systems (OCONTSOLAR) funded by the European Research Council, is being conducted by the ”Model Predictive Group” at the Systems and Automatics Engineering Department, University of Seville, Spain. This thesis is framed and financed by the OCONTSOLAR project. One of the main objectives of this project is to develop radically new model predictive control (MPC) algorithms which use mobile solar sensor to obtain estimations and predictions of the solar radiation mapping [108]. In general, the control strategies proposed in the literature use the direct solar radiation provided by pyrheliometers. They consider that this measurement is the same for the whole solar field. If the solar field is small, such as the ACUREX field, this assumption can be considered reasonable, but in large scale plants different levels of solar radiation affect the solar field due to passing clouds. Furthermore, the efficiency of the loops can be substantially different because a group of them has been cleaned whereas another has not been [70]. The most efficient loops have to be defocused to avoid excessive temperatures. Paradoxically, the most efficient loops will have the higher energy losses. To avoid this energy loss, the valves of the most efficient loops would have to be opened to increase the HTF flow. However, any movement of the valve in one of the loops will influence the flow of the rest of the loops. Loop valves are only used in current plants for steady state flow balancing. 1.3 Thesis Objectives The objectives of the thesis are mainly focused on the development of new model predictive control strategies and optimization of large-scale parabolic solar plants. In the first part novel control algorithms are proposed for the control of the defocusing of the parabolic trough collectors as well as a control strategy for an operation mode which has not previously been taken into account in the literature, power limitation, which restricts the possibility of maximizing the production of the plant since the plant is forced to reduce the efficiency of the solar field in certain seasons of the year. It is important to keep the oil temperature within the safe temperature limits provided by the manufacturer and to avoid its degradation. To prevent the temperature from exceeding the degradation limit, commercial plants 1.3. Thesis Objectives 23 have a safety strategy: collector defocus. This is applied in a staggered manner as total or partial defocusing, which is done by modifying the angle of the collector. Given that, in commercial plants, this is considered to be a safety mechanism, it is usually carried out only on the basis of thresholds. This mechanism of defocusing the collector is reactive and highly inefficient due to the thermal jumps caused by the application of full or partial defocusing and may cause oscillations in the outlet temperature of the loops [109,110]. Commercial plants loops are, generally, formed by four solar collectors. Generally, the defocusing action of the fourth collector is, in normal situations, sufficient to control the outlet temperature around a set-point, avoiding exceeding the established thermal limit. MPC controllers are proposed for the defocusing of the fourth collector, avoiding the use of techniques based on partial or total defocuses. The proposed controllers have been simulated using a nonlinear model of a large-scale 50 MW plant composed of 90 loops. However, when power limitations appear, defocusing the fourth collector is not enough to avoid exceeding the safe thermal limit of the HTF, and the third collector defocus has to be applied. This is due to the fact that to reduce power it is necessary to reduce the HTF flow-rate of the plant as it is directly related to the electrical power generation [109,110]. Generally, power limitations appear when the electrical grid is saturated and they are commanded by the Transmission System Operator (TSO). In these cases, the plant is forced to decrease its electric production and maintain the power set-point determined by the TSO. Therefore, maximum power production no longer makes sense. In this situation the objective is double: fulfilling the TSO power set-point and temperature tracking. Therefore, controllers are also proposed for the defocusing of the third collector and for the control of the electrical power generation in cases of power limitation. The second part of the thesis focuses on the optimization of parabolic trough plants. The optimization of large-scale solar trough plants operation poses important challenges which requires new advanced control techniques to address them: • The optical efficiency of different groups of loops may be substantially different in large scale solar plants. The most efficient loops will probably have to be defocused to avoid excessive temperatures. 24 Introduction Paradoxically, the most efficient loops will have the higher energy losses because of defocusing. To avoid this energy loss, the valves of the most efficient loops would have to be opened to increase the HTF flow. However, any movement of the valve in one of the loops will influence the flow of the rest of the loops. Loop valves are only used in current plants for steady state flow balancing. • Scattered clouds may only affect the locations where the sensors are placed, while the rest of the plant may be under the effect of intense DNI, or vice versa. Sudden changes in DNI produced by scattered clouds induce oscillations so severe that the solar field may have to be defocused or shutdown. This fact tends to cause not only energy losses but plant deterioration. A spatially distributed DNI nowcasting can be used to improve the plant operation and optimizing the production. Optimization algorithms are proposed for the control of the inlet valves of each of the loops in order to homogenize the solar field outlet temperature, or what is the same, a thermal balance of the solar field. This balance will decrease the high temperature gradients between the different loops as well as the possible loss of energy due to the defocusing of part of the solar field with higher optical efficiency or higher solar radiation. The proposed algorithms have been simulated using the ACUREX model in a preliminary phase, and adapted to a 50 MW plant in a second phase as a case study. It is shown how the proposed strategy is valid both for when the thermal unbalance is produced by the difference in the overall efficiency of the loops and when there is cloudiness affecting only part of the solar field (spatially distributed DNI). In addition to the controllers, different techniques have been used to perform estimations of the temperature and reflectivity or efficiency of the loops. Model based and data driven or machine learning observers have been applied such as the Lumped Parameter Model, the Unscented Kalman Filter, Classification and Regression Trees and Fuzzy Logic. 1.4 Thesis Organization This section details the organization of the Thesis. 1.4. Thesis Organization 25 In Chapter 2, the models of parabolic trough solar plants used as the main platform for the simulation and design of controllers and estimators are described: (1) the ACUREX solar field sited in Almer ´ ıa (Spain) (PSACIEMAT), and (2) a large-scale 50 MW solar plants composed of 90 loops. The mathematical description of concentrated and distributed parameters models that are used throughout the thesis are described. Chapter 3presents Gain-Scheduling Generalized Predictive Controllers (GS-GPC) based on events for the control of the defocusing of the third and fourth collectors of the loops of a 50 MW field as well as for the control of the generated electrical power when the TSO commands power limitations. Chapter 4 proposes an improvement in the control of the defocus control proposed in chapter 3through the use of an adaptive state space MPC control scheme. Among the improvements, an estimator of the reflectivity of the collectors by means of the concentrated parameters model to improve the adaptation of the model of the space of states at each moment as well as the estimation of the states of the plant (temperatures) using an Unscented Kalman Filter (UKF) are proposed. In Chapter 5a preliminary study of the proposal of a nonlinear model based optimization algorithm for the homogenization of the outlet temperature of the solar field loops of ACUREX. This is done by controlling the inlet valves taking into account that the global optical efficiency of the loops do not have to be the same. The estimation of the temperatures is obtained by means of a Classification and Regression Trees which is previously trained with data obtained from the simulation of the distributed parameter model. In Chapter 6modifications are made and the optimization algorithm is adapted for use in large-scale solar trough plants. It is applied to a 50 MW plant simulation model and the advantages of applying a model-based optimization technique for the control of the inlet valves of the loops are presented, as well as the problems that arise in plants with a large number of loops. Chapter 7presents the problem of solar radiation in large and very large-scale fields. Due to the large size of solar plants, parts of the field may be covered by clouds (temporarily) while other areas will not. Similarly, it is possible to apply online optimization algorithms for the control of the inlet valves as the cloudiness penetrates certain areas of the field, avoiding energy losses, high temperature gradients and decreasing the defocus control actions. Finally, Chapter 8presents the 26 Introduction conclusions of the work carried out in this Thesis and future lines of work. 1.5 Main Contributions This thesis focuses on the following contributions: 1. Development of Gain-Scheduling Generalized Predictive Controllers (GS-GPC) based on events for the defocus control of the third and fourth collectors of the loops of a 50 MW commercial solar plant as well as for the control of the generated electrical power when the TSO commands power limitations. 2. Improvement of the GS-GPC defocus control strategy through the use of an adaptive state space MPC control scheme. 3. An estimator of the reflectivity of the collectors by means of the concentrated parameters model to improve the adaptation of the model of the space of states at each moment as well as an estimator of the states of the loops temperatures using an Unscented Kalman Filter (UKF) are proposed. 4. A nonlinear model based optimization algorithm for the homogenization of the outlet temperature of the solar field loops. This is done by controlling the inlet valves taking into account that the global optical efficiency of the loops do not have to be the same. 5. An estimator of the temperatures by means of a Classification and Regression Trees, previously trained with data obtained from the simulation of the distributed parameter model, is proposed. 6. A model-based optimization technique for the control of the inlet valves of the loops to achieve a solar field thermal balance, in large scale solar trough plants, is presented. A clustering technique is applied to reduce the number of decision variables in the optimization problem. 7. An online optimization algorithm for the control of the loops inlet valves is designed, for large scale solar trough plants, to thermally 1.6. Publications List 27 balance the solar field as the cloudiness penetrates certain areas of the solar field, avoiding energy losses, high temperature gradients and decreasing the defocus control actions. 1.6 Publications List The following articles have been published or sent for publication during the development of this thesis: Journal articles: 1. A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Temperature homogenization of a solar trough field for performance improvement. Solar Energy., 165C:1–9, May 2018. 2. A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Event based mpc for defocusing and power production of a parabolic trough plant under power limitation. Solar Energy, 174:570–581, November 2018. 3. A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Adaptive incremental state space mpc for collector defocusing of a parabolic trough plant. Solar Energy, 184:105–114, May 2019. 4. A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Thermal balance of large scale parabolic trough plants: A case study. Solar Energy, (Submitted to Solar Energy, Elsevier, June 2019). Conference articles: 1. A. J. S ´ anchez, J. M. Esca ˜ no, N. Canty, A. J. Gallego and E. F. Camacho. Solar radiation estimator and fault tolerant model predictive control of a parabolic-trough field. 26th Irish Signals and Systems Conference (ISSC), pp:1-7, June 2015. 2. A. J. Gallego, L. Yebra, E. F. Camacho and A. J. S ´ anchez. Mathematical modeling of the parabolic trough collector field of the tcp-100 28 Introduction research plant. In Conference: Modelling and Simulation - 9th EUROSIM 2016, at Oulu, Finland, September 2016. 3. A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Fault tolerant MPC of a solar trough field based on classification and regression trees. 3rd Conference on Control and Fault-Tolerant Systems (SysTol), pp:152-157, September 2016. 4. A. J. S ´ anchez, J. M. Esca ˜ no, C. Bordons and E. F. Camacho. Estimador borroso de una planta solar cilindro-parab ´ olica. XXXIX Jornadas de Autom ´ atica, Universidad de Extremadura, Badajoz, Espa ˜ na, Septiembre de 2018. Award: Best paper in intelligent control. Book chapters: 1. E. F. Camacho, A. J. S ´ anchez and A. J. Gallego. Chapter title: Model predictive control of large scale solar trough plants. Book: Solar Energy Systems: Progress and Future Directions. Nova Science Publishers, Inc, New York, USA, (Submitted and accepted for publication). Chapter 2 Solar trough field models 1. A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Event based mpc for defocusing and power production of a parabolic trough plant under power limitation. Solar Energy, 174:570–581, November 2018. 29 155 temperature thresholds. However, applying the defocus control in this way can lead to oscillations in the outlet temperature of the loops. In this thesis two types of controllers have been presented for the defocusing of the collectors for loop outlet temperature set-point tracking instead of applying control signals based on thresholds. An Event-based GS-GPC and an Adaptive Incremental State Space MPC. Both controllers have been simulated in the 50 MW model. The two controllers have shown a good behavior in the tracking of the loop outlet temperature with a smooth control, avoiding reaching the security limit. The state space MPC has shown a similar level of control performance to that of the EGS-GPC although with certain improvements. The EGS-GPC controller is somewhat inferior (ITAE and ISE indexes values are higher) to the MPC in the state space since, although it is based on global linear models at different operating points (Gain Scheduling), its performance may vary due, among others, to modeling errors and other factors that are considered in the EGS-GPC as the efficiency of the loops. In general, in the literature, the efficiency of the loops has been assumed constant and equal in all loops. This can be a valid assumption for smallsize solar fields. However, throughout the thesis it has been taken into account that the efficiency can be different in each loop due to cleaning factors (reflectivity) and condition of the mirrors. It is convenient to take into account this type of differences in large scale solar fields since they may affect the control performance of the solar plant. For this, different estimators of the efficiency of the loops (or collectors) have been applied depending on the needs of each control strategy. In some cases, an efficiency estimation has been applied by means of an optimization based on a non-linear model due to the fact that sufficient time was available for its computation. In other cases, the concentrated parameter model has been used due to the demands of the sampling times of the controllers. Both have proven to provide good results in the estimation and have been included into the control strategies in order to have an adaptive control model. Throughout the thesis it has been commented that on some occasions, a commercial solar plant may receive a power limitation from the TSO. These commands make the plant move into an operation mode in which the objective is not maximum production. In these cases, the plant is forced to decrease its electric production and maintain the power set-point commanded by the TSO. Therefore a double objective problem arises: (1) temperature tracking and (2) generated power set-point tracking. The plant will have a time period to reduce its generated electric power to the set-point determined 156 Conclusion by the TSO. In this mode, the flow-rate should be decreased until the power set-point is achieved but at the cost of increasing the outlet temperature. These orders must be met, otherwise the plant will face economic sanctions. For the control of electrical power, an Event-based GS-GPC control scheme has been proposed in this thesis. In this case, given that the plant has a time to reach the established limit, it is possible to apply to the predictive controller a ramp of future references in each sampling time which allows a better control of the power thanks to the use of the sliding horizon of the model predictive control scheme. Small solar fields such us the ACUREX field can be considered as one single equivalent collector loop in order to obtain an overall dynamic model and design control strategies in which the temperature to be controlled is the weighted average temperature of the loops. This approach is reasonable when the solar field is relatively small. Parameters such as efficiency of the receiver tubes, reflectivity or the shape factor of the collectors can be considered quite similar for all loops. Another important point is that the DNI is measured locally by pyrheliometers. Considering that all the solar field is affected by similar levels of solar radiation is reasonable in small plants but not in large solar fields. In large scale solar plants the flow-rate controller, usually design as it has been stated above, may cause energy losses or unnecessary overreaction due to the temperature discrepancy between the different loops. This effect occurs not only due to the difference between the loops but also due to clouds. Controlling the average temperature of all loops presents a drawback: if the solar field does not have a good thermal balance, due to differences in the loops optical efficiencies or partial clouds over the solar field, some loops may reach much higher temperatures than the others. The optical efficiency of the loops can be different due to aspects such as reflectivity, efficiency of the tube, shape factor and structural state of the collectors/loops. These imbalances in the optical efficiency may lead to the most efficiency loops to be defocused to avoid overheating problems. In this thesis, control strategies based on non-linear optimizations have been presented for the control of the flow inlet valves to each of the loops. In commercial solar trough plants, the input valves are generally used to achieve a hydraulic or static thermal balance. Although this is normally done on rare occasions. It has been shown, by simulations, that a thermal balance of the solar field can be achieved (homogenizing loop temperatures) by applying the proposed optimization algorithm based on a simplified nonlinear model to control the inlet valves aperture. In addition to avoid large temperature gradients between the loops, it is also possible to avoid energy 157 losses (in certain occasions) as well as to reduce or completely eliminate the defocusing actions applied to each of the collectors, which will imply a life extension of the actuator. For the design of these optimization strategies, a clustering algorithm (K-Means) was used to group the loops according to their similarity in certain variables and parameters. Decreasing the number of decision variables in a non-linear optimization problem reduces computation times and can be adjusted to the desired control sampling times. This strategy has been applied for cases of differences in the efficiency of the loops as well as for cases in which partial clouds covers part of the solar field leaving another part of the field free of transients. In each of the cases, different sampling times have to be used since the needs of the problem to be solved are different. Future Works Some of the research lines to follow as a continuation of this thesis are: 1. Complete analysis of the effects of defocusing. Although two MPC strategies have been developed for defocus control, it is a system that has barely been touched in the literature. A deeper analysis of the problem is necessary. In this thesis has been applied to the third and fourth collector, which is the normal procedure in commercial plants. However, it may be interesting to analyze the first and second collectors. It will be studied if by applying fixed defocusing signals in the first two collectors, similar results can be obtained with lower control signals from the actuators. The temperature reference setpoints will be studied more carefully to find the optimal temperature references that may vary according to the operating modes. 2. Distributed optimization strategies for thermal balance of the solar field (OCONTSOLAR project). Valve control would be carried out through a distributed optimization based on nonlinear models. In this way, it is expected to be able to further reduce the computation time of the optimization problems in order to increase the number of decision variables that can be used. It is an important point because there are plants with up to 808 loops which will require high computation times. 3. Optimization strategies for thermal balance using game theory (OCONTSOLAR project). Similar to how it was done with the clustering of the valves, coalitions could be applied to the solar field loops. The 158 Conclusion number of agents in the coalitions would change dynamically depending on factors such as reflectivity, defocus and cloudiness among others. 4. Integration of the proposed strategies with the radiation estimates generated by UAVs (OCONTSOLAR project). By making use of the spatial estimation of the solar radiation of the field and/or having a short term forecasting, much more anticipatory controllers could be obtained to avoid energy losses loss, more stable plant operation and to avoid undesired effects in the temperatures of the solar field. 5. Another point to develop is in the topic of temperature estimates. In this thesis data-based techniques (CARTs) and model-based (UKF) observers have been used. It is true that the performance of the UKF is indisputable due, among others, to the use of a simplified nonlinear model and its adaptation to the state of the plant at all times. However, it has the disadvantage of the computation time. Given that the OCONTSOLAR project is focused on the development of new control algorithms for plants with large extensions, the number of loops to control will be very high. When applying new controllers the estimates of the temperatures of the loops segments may exceed the sampling times of the controllers. Although the observer based on CART has shown good results, the size that it can occupy in memory is a disadvantage. Although fast-track observers have already begun to be develop, as a fuzzy observer in [161], it remains an observer based on static data and therefore research on data driven observers with reduced computation times that can include plant dynamics will be continued. References [1] E. F. Camacho and A. J. Gallego. Optimal operation in solar trough plants: A case study. Solar Energy, 95:106 – 117, 2013. [2] Stanford University. Global Climate and Energy Project. http://gcep.stanford.edu/research/exergy/resourcechart.html, 2018. (Accessed June 1, 2019). [3] H. Ritchie and M. Roser. Our World in Data. Energy Production & Changing Energy Sources. https://ourworldindata.org/energyproduction-and-changing-energy-sources. (Accessed June 1, 2019). [4] K. Butti and J. Perlin. A Golden Thread. 2500 Years of Solar Architecture and Technology. Cheshire Books, Palo Alto. Van Nostrand Reinhold Company., 1980. [5] European Commission. European Research on Concentrated Solar Thermal Energy. https://ec.europa.eu/research/energy/pdf/cst en.pdf, 2004. (Accessed June 1, 2019). [6] J. M. Gordon, editor. Solar Energy: The State of the Art. Earthscan, 2013. [7] J. Ericsson. Contributions to the Centennial Exhibition. Printed by John Ross, 1876. [8] J. R. Hirschmann. Solar distillation in chile. Desalination, 17(1):31 – 67, 1975. [9] DeGolyer Library. Sola-motor at the Ostrich Farm, Southern Methodist University (SMU), 1902. (Accessed June 1, 2019). [10] R. Thomas. The Quarterly Magazine. South Pasadena Local History: One of the Strangest Sights in America. 159 160 References https://www.southpasadenaca.gov/home/showdocument?id=10169, 2001. (Accessed June 1, 2019). [11] G. Jones. Profits and Sustainability: A History of Green Entrepreneurship. Oxford University Press, 2017. [12] The New York Times. American inventor uses egypt’s sun for power, 1916. https://timesmachine.nytimes.com/timesmachine/1916/07/02/ 104680095.pdf. (Accessed June 1, 2019). [13] J. Perlin. From Space to Earth. The Story of Solar Electricity. Harvard University Press, 2002. [14] D. Y. Goswami, F. Kreith, and J. F. Kreider. Principles of Solar Engineering. Taylor & Francis, 2th edition, 2000. [15] M. Romero and J. Gonz ´ alez-Aguilar. Solar thermal csp technology. Wiley Interdisciplinary Reviews: Energy and Environment, 3(1):42– 59, 2014. [16] M. J. Blanco and L. R. Santigosa. Advances in Concentrating Solar Thermal Research and Technology. Woodhead Publishing, 1 edition, 2017. [17] A. J. Gallego, L. Yebra, E. F. Camacho, and A. J. S ´ anchez. Mathematical modeling of the parabolic trough collector field of the tcp-100 research plant. In Conference: Modelling and Simulation - 9th EUROSIM 2016, At Oulu, Finland, pages 912–918, 09 2016. [18] DLR. Institute of Solar Research. The Experimental Solar Power Plant Juelich. https://www.dlr.de/sf/en/desktopdefault.aspx/tabid- 8560/15527 read-44868/. (Accessed June 1, 2019). [19] Atlantica Yield. Atlantica Yield Sustainable Infraestructure. Our Assets. https://www.atlanticayield.com/web/en/company-overview/ourassets/, 2017. (Accessed June 1, 2019). [20] P. Breeze. Chapter 5 - Solar Towers. In Paul Breeze, editor, Solar Power Generation, pages 35 – 40. Academic Press, 2016. [21] Koza1983, Wikimedia. Ps20 and ps10 solar tower plants. creative commons attribution 3.0 unported license. https://commons.wikimedia.org/wiki/File:PS20andPS10.jpg, September 2007. (Accessed June 1, 2019). References 161 [22] Amble, Wikimedia. Crescent Dunes Solar Energy Facility. Attribution-Share Alike 4.0 International license., 2014. https://commons.wikimedia.org/wiki/File:Crescent Dunes Solar December 2014.JPG (Accessed June 1, 2019). [23] F. J. Collado and J. Guallar. A review of optimized design layouts for solar power tower plants with campo code. Renewable and Sustainable Energy Reviews, 20:142 – 154, 2013. [24] J. Ballestr ´ ın, E. Carra, R. Monterreal, R. Enrique, J. Polo, J. Fern ´ andez-Reche, J. Barbero, A. Marzo, J. Alonso-Montesinos, G. L ´ opez, and F. J. Batlles. One year of solar extinction measurements at plataforma solar de almer ´ ıa. application to solar tower plants. Renewable Energy, 136:1002 – 1011, 2019. [25] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Planta Solar 10. https://solarpaces.nrel.gov/plantasolar-10, 2017. (Accessed June 1, 2019). [26] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Planta Solar 20. https://solarpaces.nrel.gov/plantasolar-20, 2017. (Accessed June 1, 2019). [27] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Supcon Solar Project. https://solarpaces.nrel.gov/supcon-solar-project, 2019. (Accessed June 1, 2019). [28] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Crescent Dunes Solar Energy Project. https://solarpaces.nrel.gov/crescent-dunes-solar-energy-project, 2016. (Accessed June 1, 2019). [29] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Ivanpah Solar Electric Generating System. https://solarpaces.nrel.gov/ivanpah-solar-electric-generating- system, 2014. (Accessed June 1, 2019). [30] Sbharris, Wikimedia. Ivanpah Solar Electric Generating System (ISEGS) (California). Creative Commons Attribution-Share Alike 3.0 Unported license., 2014. 162 References https://en.wikipedia.org/wiki/File:Ivanpah Solar Power Facility Online.jpg. (Accessed June 1, 2019). [31] G. Morin, J. Dersch, W. Platzer, M. Eck, and A. H ¨ aberle. Comparison of linear fresnel and parabolic trough collector power plants. Solar Energy, 86(1):1 – 12, 2012. [32] G. Zhu, T. Wendelin, M. J. Wagner, and C. Kutscher. History, current state, and future of linear fresnel concentrating solar collectors. Solar Energy, 103:639 – 652, 2014. [33] Soltigua. Linear Fresnel Collector – Manufacturer. http://www.soltigua.com/flt-introducion/. (Accessed June 1, 2019). [34] SUNCINM. Solar Steam Generators. https://www.suncnim.com/en. (Accessed June 1, 2019). [35] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Puerto Errado 1 Thermosolar Power Plant. https://solarpaces.nrel.gov/puerto-errado-1-thermosolar-power- plant, 2011. (Accessed June 1, 2019). [36] E. F. Camacho, A. J. Gallego, A. J. S ´ anchez, and M. Berenguel. Incremental state-space model predictive control of a fresnel solar collector field. Energies, 12:3, 12 2018. [37] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Puerto Errado 2 Thermosolar Power Plant. https://solarpaces.nrel.gov/puerto-errado-2-thermosolar-power- plant, 2013. (Accessed June 1, 2019). [38] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Dhursar. https://solarpaces.nrel.gov/dhursar, 2014. (Accessed June 1, 2019). [39] Novatec Solar, Wikimedia. Novatec Fresnel Power Station Puerto Errado 2. Creative Commons Attribution-Share Alike 3.0 Unported license, 2012. https://commons.wikimedia.org/wiki/File:Novatec Solar Puerto Errado 2 BoP PI.jpg (Accessed June 1, 2019). References 163 [40] E. Rakhshani, K. Rouzbehi, A. J. S ´ anchez, A. Cabrera-Tobar, and E. Pouresmaeil. Integration of large scale pv-based generation into power systems: A survey. Energies, 12:1425, 04 2019. [41] Z22, Wikimedia. Parabolic trough at Harper Lake in California (Mojave Solar Project). Creative Commons Attribution-Share Alike 3.0 Unported license., 2013. https://commons.wikimedia.org/wiki/File:Parabolic trough at Harper Lake in California.jpg (Accessed June 1, 2019). [42] E. Zarza, L. Valenzuela, J. Le ´ on, K. Hennecke, M. Eck, H.-Dieter Weyers, and M. Eickhoff. Direct steam generation in parabolic troughs: Final results and conclusions of the diss project. Energy, 29(5):635 – 644, 2004. SolarPACES 2002. [43] E. F. Camacho, M. Berenguel, F. R. Rubio, and D. Mart ´ ınez. Control of Solar Energy Systems. Springer-Verlag London, 1st edition, 2012. [44] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Helioenergy 1. https://solarpaces.nrel.gov/helioenergy-1, 2017. (Accessed June 1, 2019). [45] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Solana Generating Station. https://solarpaces.nrel.gov/solana-generating-station, 2015. (Accessed June 1, 2019). [46] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Solaben 3. https://solarpaces.nrel.gov/solaben- 3, 2017. (Accessed June 1, 2019). [47] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. NOOR II. https://solarpaces.nrel.gov/noor-ii, 2018. (Accessed June 1, 2019). [48] D. Kearney, U. Herrmann, P. Nava, B. Kelly, R. Mahoney, J. Pacheco, R. Cable, N. Potrovitza, D. Blake, and H. Price. Assessment of a molten salt heattransfer fluid in a parabolictrough solar field. Journal of Solar Energy Engineering, 125(2):170–176, May 2003. 164 References [49] D. Kearney, B. Kelly, U. Herrmann, R. Cable, J. Pacheco, R. Mahoney, H. Price, D. Blake, P. Nava, and N. Potrovitza. Engineering aspects of a molten salt heat transfer fluid in a trough solar field. Energy, 29(5):861 – 870, 2004. SolarPACES 2002. [50] M. Liu, N. H. S. Tay, S. Bell, M. Belusko, R. Jacob, G. Will, W. Saman, and F. Bruno. Review on concentrating solar power plants and new developments in high temperature thermal energy storage technologies. Renewable and Sustainable Energy Reviews, 53:1411 – 1432, 2016. [51] G. Alva, L. Liu, X. Huang, and G. Fang. Thermal energy storage materials and systems for solar energy applications. Renewable and Sustainable Energy Reviews, 68:693 – 706, 2017. [52] U. Pelay, L. Luo, Y. Fan, D. Stitou, and M. Rood. Thermal energy storage systems for concentrated solar power plants. Renewable and Sustainable Energy Reviews, 79:82 – 100, 2017. [53] I. Sarbu and C. Sebarchievici. A comprehensive review of thermal energy storage. Sustainability, 10(1):1–32, January 2018. [54] C. Prieto, A. Rodr ´ ıguez, D. Pati ˜ no, and L. F. Cabeza. Thermal energy storage evaluation in direct steam generation solar plants. Solar Energy, 159:501 – 509, 2018. [55] L. Roca, J. Bonilla, M. M. Rodr ´ ıguez-Garc ´ ıa, P. Palenzuela, A. de la Calle, and L. Valenzuela. Control strategies in a thermal oil – molten salt heat exchanger. AIP Conference Proceedings, 1734(1):130017, 2016. [56] G. Peir ´ o, C. Prieto, J. Gasia, A. Jov ´ e, L. Mir ´ o, and L. F. Cabeza. Two-tank molten salts thermal energy storage system for solar power plants at pilot plant scale: Lessons learnt and recommendations for its design, start-up and operation. Renewable Energy, 121:236 – 248, 2018. [57] E. Gonz ´ alez-Roubaud, D. P ´ erez-Osorio, and C. Prieto. Review of commercial thermal energy storage in concentrated solar power plants: Steam vs. molten salts. Renewable and Sustainable Energy Reviews, 80:133 – 148, 2017. References 171 [117] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Ibersol Ciudad Real (Puertollano). https://solarpaces.nrel.gov/ibersol-ciudad-real-puertollano, 2017. (Accessed June 1, 2019). [118] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Helios I. https://solarpaces.nrel.gov/helios-i, 2013. (Accessed June 1, 2019). [119] National Renewable Energy Laboratory (NREL). Concentrating Solar Power Projects. Palma del R ´ ıo I. https://solarpaces.nrel.gov/palmadel-rio-i, 2013. (Accessed June 1, 2019). [120] S. Rohani, T. P. Fluri, F. Dinter, and P. Nitz. Modelling and simulation of parabolic trough plants based on real operating data. Solar Energy, 158:845 – 860, 2017. [121] National Renewable Energy Laboratory (NREL). System Advisor Model (SAM). https://sam.nrel.gov/, 2017. (Accessed June 1, 2019). [122] M. Geyer, E. L ¨ upfert, R. Osuna, A. Esteban, W. Schiel, A. Schweitzer, E. Zarza, P. Nava, J. Langenkamp, and E. Mandelberg. Eurotrough - parabolic trough collector developed for cost efficient solar power generation. In 11th SolarPACES International Symposium on Concentrated Solar Power and Chemical Energy Technologies, September 2002. [123] D. W. Kearney. Parabolic trough collector overview. Parabolic trough work-shop, NREL., 2007. [124] CIEMAT, CRES, DLR, FICHTNER, FLABEG, INABENSA, and SBP. Development of a low cost european parabolic trough collector - eurotrough. Technical report, European Commision. CORDIS. FP4- NNE-JOULE C. Grant agreement ID: JOR3980231., 2001. [125] SCHOTT Solar CSP GmbH. Schott ptr R  70 receivers. https://www.us.schott.com/csp/english/schott-solar-ptr-70- receivers.html, September 2018. (Accessed June 1, 2019). [126] Therminol VP1 HTF. https://www.therminol.com/products/Therminol- VP1, 2019. (Accessed June 1, 2019). 172 References [127] H. Schenk, J. Dersch, T. Hirsch, and T. Polklas. Transient simulation of the power block in a parabolic trough power plant. In The 11th International Modelica Conference Versailles, France, pages 605– 614. Link ¨ oping University Electronic Press, Link ¨ opings universitet, September 2015. [128] R. M. Monta ˜ n ´ es, J. Windahl, J. Palsson, and M. Thern. Dynamic modeling of a parabolic trough solar thermal power plant with thermal storage using modelica. Heat Transfer Engineering, 39(3):277–292, 2018. [129] F. Burkholder, M. Brandemuehl, H. Price, J. Netter, C. Kutscher, and E. Wolfrum. Parabolic trough receiver thermal testing. In Energy Sustainability, ASME 2007 Energy Sustainability Conference, pages 961–970, 2007. [130] E. L ¨ upfert, K. J. Riffelmann, H. Price, F. Burkholder, and T. Moss. Experimental analysis of overall thermal properties of parabolic trough receivers. Journal of Solar Energy Engineering, 130(2), May 2008. [131] M. Beschi, S. Dormido, J. Sanchez, A. Visioli, and L. J. Yebra. Eventbased pi plus feedforward control strategies for a distributed solar collector field. IEEE Transactions on Control Systems Technology, 22(4):1615 – 1622, July 2014. [132] B. Khoukhi, M. Tadjine, and M. S. Boucherit. Nonlinear continuoustime generalized predictive control of solar power plant. Int. J. Simul. Multisci. Des. Optim., 6, May 2015. [133] M. J. Blanco and S. Miller. 1 - introduction to concentrating solar thermal (cst) technologies. In Manuel J. Blanco and Lourdes Ramirez Santigosa, editors, Advances in Concentrating Solar Thermal Research and Technology, Woodhead Publishing Series in Energy, pages 3 – 25. Woodhead Publishing, 2017. [134] M. J. Montes, A. Ab ´ anades, J. M. Mart ´ ınez-Val, and M. Vald ´ es. Solar multiple optimization for a solar-only thermal power plant, using oil as heat transfer fluid in the parabolic trough collectors. Solar Energy, 83(12):2165 – 2176, 2009. [135] N. B. Desai and S. Bandyopadhyay. Optimization of concentrating solar thermal power plant based on parabolic trough collector. Journal of Cleaner Production, 89:262 – 271, 2015. References 173 [136] A. J. S ´ anchez, A. J. Gallego, J. M. Esca ˜ no, and E. F. Camacho. Temperature homogenization of a solar trough field for performance improvement. Solar Energy., 165C:1–9, May 2018. [137] E. F. Camacho and C. Bordons. Model Predictive control. Springer- Verlag London, 2th edition, 2007. [138] E. F. Camacho, M. Berenguel, and F. R. Rubio. Application of a gain scheduling generalized predictive controller to a solar power plant. Control Engineering Practice, 2(2):227–238, 1994. [139] M. J. Vasallo and J. M. Bravo. A mpc approach for optimal generation scheduling in csp plants. Applied Energy, 165:357 – 370, 2016. [140] T. A. Elias, P. C. Mendes, J. E. Normey-Rico, and C. B. Alba. Optimal solar collectors defocusing based on maximum temperature. In 9th International Renewable Energy Congress (IREC), Hammamet, Tunisia., 2018. [141] U. Maeder, F. Borrelli, and M. Morari. Linear offset-free model predictive control. Automatica, 45:2214–2222, 2009. [142] D. Limon, I. Alvarado, T. Alamo, and E. F. Camacho. Robust tubebased mpc for tracking of constrained linear systems with additive disturbances. Journal of Process Control, 20:248–260, 2010. [143] D. di Ruscio. Model predictive control with integral algorithm: A simple mpc algorithm. Modeling, Identification and Control, 34 (3):119–129, 2013. [144] R. E. Kalman. A new approach to linear filtering and predictions problem. Trans. ASME Jounal of Basic Engineering, 82:35–60, 1960. [145] S. Haykin. Kalman Filtering and Neural Networks. A Wiley- Interscience Publication, 2001. [146] A. Romanenko and J. A. Castro. The unscented Kalman filter as an alternative to the EKF for nonlinear state estimation: a simulation case study. Computers and Chemical Engineering, 28:347–355, 2004. [147] M. St-Pierre and D. Gringras. Comparison between the unscented kalman filter and the extended kalman filter for the position estimation module of an integrated navigation information system. IEEE Intelligent Vehicles Symposium, 0-7803-8310-9/04:831–835, 2004. 174 References [148] L. Li, C. Hua, and H. Yang. A new adaptive unscented kalman filter based on covariance matching technique. In International Conference on Mechatronics and Control (ICMC), 2014. [149] Y. Zhou, C. Zhang, Y. Zhang, and J. Zhang. A new adaptive squareroot unscented kalman filter for nonlinear systems with additive noise. International Journal of Aerospace Engineering, pages 1–10, 2015. [150] Y. Wang, Z. Qiu, and X. Qu. An improved unscented kalman filter for discrete nonlinear systems with random parameters. Discrete Dynamics in Nature and Society, pages 1–10, 2017. [151] C. M. Cirre, M. Berenguel, L. Valenzuela, and R. Klempous. Reference governor optimization and control of a distributed solar collector field. European Journal of Operational Research, 193:709–717, 2009. [152] L. Breiman, J. Friedman, C. J. Stone, and R. A. Olshen. Classification and Regression Trees. The Wadsworth and Brooks-Cole statisticsprobability series. Taylor & Francis, 1984. [153] A. Malinowska. Classification and regression tree theory application for assessment of building damage caused by surface deformation. Natural Hazards, 73(2):317–334, 2014. [154] R. Siciliano and F. Mola. Multivariate data analysis and modeling through classification and regression trees. Computational Statistics & Data Analysis, 32(3–4):285 – 301, 2000. [155] A. Nai Oleari, J. R. D. Frejo, E. F. Camacho, and A. Ferrara. A model predictive control scheme for freeway traffic systems based on the classification and regression trees methodology. In 2015 European Control Conference (ECC), pages 3459–3464, July 2015. [156] A. J. S ´ anchez, J. M. Esca ˜ no, A. J. Gallego, and E. F. Camacho. Fault tolerant mpc of a solar trough field based on classification and regression trees. In 3rd Conference on Control and Fault-Tolerant Systems (SysTol), pages 152–157, September 2016. [157] G. J. VanWylen, R. E. Sonntag, and C. Borgnakke. Thermodynamics. McGraw-Hill, New York, 2th edition, 2002. References 175 [158] Y. A. C¸ engel and M. A. Boles. Thermodynamics: An engineering approach. McGraw-Hill, New York, 8th edition, 2015. [159] Dan Simon. Optimal state estimation: Kalman, H infinity, and nonlinear approaches. John Wiley & Sons, Inc, Hoboken, New Jersey., 2006. [160] S. P. Lloyd. Least squares quantization in pcm. IEEE Transactions on Information Theory, 28(2):129–137, 1982. [161] A. J. S ´ anchez, J. M. Esca ˜ no, C. Bordons, and E. F. Camacho. Estimador borroso de una planta solar cilindro-parab ´ olica. In XXXIX Jornadas de Autom ´ atica, Universidad de Extremadura, Badajoz, Espa ˜ na, September 2018.