Abstract

The boiler-turbine unit is really a complex system in thermal power engineering due to its large-scale nonlinearity, unmeasured state, unknown disturbances, and constraints imposed on both controls and outputs. To design a controller with appropriate performance in above synthetical cases, this paper intends to propose an adaptively receding Galerkin optimal controller design method, in which, the mathematical dynamics of unit can be directly used as a predictive model without any linearization, and the unmeasured state in the predictive model is adaptively estimated using a predesigned state observer. With the help of a mathematical predictive model, optimal control law is then obtained based on a Galerkin optimization algorithm. Due to the application of the useful information measured at every sampling time instant, the proposed method can deal with the tracking problem with constraints rather than the stabilization problem that can be only done by the traditional Galerkin optimal control. Furthermore, it can also be easily extended to estimate and thus eliminate constant disturbances in an output channel using an independent model strategy. Some simulations suggest that satisfactory tracking performance can be achieved even when the unit experiences wide-range load change.

1. Introduction

The boiler-turbine unit plays a critical role in a thermal power plant. Due to its genuine nonlinearity, serious couplings among state variables, and physical constraints, it is difficult to design a controller with appropriate transient performance for the boiler-turbine unit [1]. In particular, some key unmeasured state variables as well as unknown disturbances bring much more difficulties in controlling the boiler-turbine unit.

Conventionally, the boiler-turbine unit was usually operated in a local load range. In this way, the proportional-integral-differential (PID) controller can achieve acceptable performance [2]. However, the unit should be run now in a large-scale load range, which leads the unit’s dynamics to be inherently nonlinear [3]. In this case, it is challenging to design a PID controller with appropriate performance, especially when some physical constraints are required for a safe and correct functioning of the unit. To achieve better performance, recent years witness a surge of interests on designing an advanced controller for the popular oil-fired drum-type boiler-turbine unit [1], for example, see [319] and the literature therein. As for the boiler-turbine unit, it is interesting and imperative to derive new control strategy that can maximize or minimize a specified control performance index while honoring the constraints imposed on the unit. Until now, several optimal controllers such as model predictive control have been investigated and achieved better performances for the unit [4, 1013, 15, 16]. Nevertheless, most of these controllers were designed on the basis of either the black-box nonlinear model identified from running data of the unit or linear models obtained by linearizing the unit’s mathematical model. More attentions recently have been paid to optimal controller design, for example, see [20]. How to design an optimal controller directly based on the unit’s mathematical model with appropriate tracking performance is still open.

When designing an optimal controller for either the unit or other nonlinear systems, it has been well recognized as an extremely challenging problem to analytically solve a state- and control-constrained nonlinear optimal control problem in particular in real-time applications. As a matter of fact, due to the genuine large-scale nonlinearity of the unit, it is widely considered to be difficult to solve an optimal controller for the unit even in nonreal time. The main difficulty arises in seeking a closed-form solution to the Hamilton-Jacobi equation or in solving the canonical Hamiltonian equations resulting from an application of the minimum principle [21]. Alternatively, a Galerkin pseudospectral method [22, 23] is one of the most efficient computational approaches, intending to solve the optimal state and control sequences through transforming the state- and control-constrained nonlinear optimal control problems into a nonlinear programming problem [18, 24]. There are some obstacles on the path to design a Galerkin optimal controller to make the unit track large-scale load demand/reference. Firstly, either the Galerkin method or other pseudospectral methods usually pay much more attentions to the stabilization problem rather than the tracking problem. To solve the tracking problem, the Galerkin method should be receding or rolling in some sense by making use of the useful information measured at every sampling time instant. Secondly, to make use of information at every sampling time instant, it requires all state variables to be measurable. It will be seen that the key unmeasured state fluid density in the drum of the unit is unmeasurable. Finally, some unknown (constant) disturbances in the output channels bring much more difficulties. How to compensate uncertainties/disturbances so as to enhance control performance can be referred, for example, to [2527].

Motivated by above statements, this paper aims to propose an adaptively receding Galerkin optimal control method for an oil-fired drum-type boiler-turbine unit with some unmeasured states as well as some unknown constant disturbances in the output channels. More precisely, a state observer is first designed to adaptively estimate the key unmeasured state so as to make the information available at every sampling time instant; then, a receding Galerkin optimal controller is constructed by sufficiently taking into account information observed at each sampling time, through borrowing the basic idea from the model predictive control method; after this, an independent model strategy is embedded into the receding Galerkin controller structure to estimate and thus eliminate the constant disturbances in the output channels. Evidently, the main contributions of this paper are in twofold: an adaptively receding Galerkin optimal control strategy with estimations of unmeasured state and unknown constant output disturbances and its application for a boiler-turbine unit.

The rest of the paper is organized as follows. Section 3 briefly recalls the Galerkin method. In Section 5, the receding Galerkin optimal control strategy is proposed after introducing the boiler-turbine unit, including the state observer and the independent model strategy. Simulation results are presented in Section 4. The last section concludes this paper.

2. Galerkin Method

The main purpose of the optimal control is to solve out the admissible control sequences minimizing a cost function based on the mathematical model of an object. The mathematical description of an optimal control problem can be described as follows: determine the state-control function pair, minimizing the following cost functional (or called a performance index) subject to the dynamics initial conditions endpoint conditions and path constraints where the running (or Lagrange) cost , the endpoint (or Mayer) cost , , , and are all Lipschitz continuous.

A Galerkin method transforms the above problem into a nonlinear programming problem through the following four steps: approximating state and control variables, discretizing the system dynamics, integrating the cost function, and discretizing other constraints.

To realize approximation or discretization, it needs to introduce the concept of the node. By converting the real-time domain into a closed interval according to a series of Legendre-Gauss-Lobatto (LGL) nodes can be calculated as the roots of where is the th order Legendre polynomial defined by . Obviously, there are LGL nodes in -space such that which correspond to in the real-time domain.

In this way, one has , , or more specifically , , . With the help of LGL nodes, a Galerkin method approximates the state and control by the th order Lagrange interpolation polynomial defined on LGL nodes as follows: where is the th order Lagrange interpolation basis function defined by

Using (9), the differential equation (2) can be approximated using the following integral formulation: with test functions . When defining test function as equal to the basis function , (10) can be rewritten as

For simplicity, define and . The following two approximate equalities can be induced: where is the Legendre differentiation matrix calculated by for , for , and for , otherwise , and are the quadrature weights, and the LGL version of quadrature weights can be calculated as

With the help of and , (11) can thus be finally simplified as

In a relatively easy way, the cost function (1) can then be approximated according to the Gauss-Lobatto integration rule as follows:

Finally, together with the following approximations the Galerkin method transforms the continuous optimal control problems (2), (3), (4), and (5) into the following discrete nonlinear programming problem: where is a constant tolerance used to guarantee feasibility of the nonlinear programming problem [21].

To solve a discrete nonlinear programming problem (17), a nonlinear programming solver such as SNOPT and IPOPT is usually used [21].

3. Main Results

3.1. Problem Formulations and Adaptively Receding Galerkin Strategy

This paper considers a 160 MW oil-fired drum-type boiler-turbine unit [1], whose flow diagram is summarized in Figure 1. The mathematical model of this unit has been established in a form of where is drum pressure (kg/cm2); is electrical output (MW); is fluid density in the drum (kg/cm3); outputs , , and are, respectively, the drum steam pressure, electrical output, and drum water level; , , and are, respectively, the normalized fuel flow rate, control valve position, and feedwater flow rate; and the coefficient and evaporation rate of steam (kg/s) are defined, respectively, as

For safety consideration, models (18), (19), (20), (21), (22), and (23) should satisfy the following constraints [1, 5]:

Remark 1. Models (18), (19), (20), (21), (22), and (23) indicate that the unit’s behaviors are genuinely nonlinear and state variables are seriously coupled. In particular, the state is unmeasurable. Furthermore, as will be seen, (constant) output disturbances will be considered. All these facts put some obstacles on the path to design a receding Galerkin optimal controller for the unit straightforwardly.

Galerkin optimal control for boiler-turbine unit: define such that , , and with expansion coefficients , the optimal control problem for the unit can be formulated as (utilizing (18), (19), (20), (21), (22), and (23)) where the outputs , output references , and ; is the controllers’ derivatives; and are positive definite weight matrices. By introducing , the states and change of controls can now satisfy and the path constraints and .

It is evident that the current Galerkin method interpreted in Section 3 is only feasible for the stabilization problem rather than the tracking problem. To deal with the tracking problem, it should take into account the useful information measured at each sampling time instant, including the information of states, outputs, and references. To solve this problem, a receding version of Galerkin optimal control strategy is proposed straightforward by borrowing the basic idea from model predictive control as below. (i)At current time instant , let the current state and control be the initial conditions, that is, , ; then, obtain the optimal discrete state and control sequences and by minimizing through the Galerkin method over the prediction horizon , where is the length of horizon. More precisely, we have with and . Note that the superscript in variables such as and is just used to distinguish optimal solutions at different sampling time instants.(ii)Apply the optimal control law on the unit and repeat the above operations in step (i) at the coming time instant .

3.2. Receding Galerkin Optimal Control with a State Observer

In order to implement the receding Galerkin optimal controller, all the states should be known in advance. However, the state variable , that is, the fluid density in the drum, cannot be measured online. Therefore, we design a state observer to estimate the unmeasured state as follows.

Proposition 1. The following observer can render the unmeasured state in the unit (18), (19), (20), (21), (22), and (23) asymptotically to its true value: where with , and constant is observer gain.

Proof. Define the error , whose derivative results in with .
Note that , then one can conclude (when , takes the maximum value), which together with leads to If , like the real range of the state , one has It means that as long as , we can achieve with a positive coefficient . In fact, it is reasonable to preset any initial value for state in its practical range. To this end, the estimated value of state can asymptotically converge to its true value.

Remark 2. To show the performance of the above state observer, a simple simulation is conducted here with  kg/cm3 and  kg/cm2 under constant control inputs . The performance is shown in Figure 2, where the true trajectory of is obtained according to (20) with initial condition  kg/cm3 and other conditions as the same as the observer. It suggests that the larger the constant is, the observer approaches to the true trajectory faster.

With the help of Proposition 1, the estimated state at the thsampling time instant will be taken as the initial condition for real to implement the receding Galerkin method. Figure 3 presents the block diagram of the receding Galerkin optimal control strategy with the state observer for the unit.

3.3. Receding Galerkin Optimal Control with an Independent Model

There are always different versions of disturbances in practice for the unit. By considering the fact that the unit is usually slowly time-varying, we just consider constant disturbance in the output channels in this paper. More precisely, we hold (18), (19), (20), and the following

In order to estimate the constant output disturbances, an independent model strategy can be introduced into the receding Galerkin optimal control strategy with the state observer, as shown in Figure 4.

It can be seen from Figure 4 that the constant output disturbance can be estimated according to where is the practical output and is the output of an independent model. Here, the independent model is simply defined as the same as models (18), (19), (20), (21), (22), and (23). To this end, the estimated disturbances are fed back instead of in (33) to implement the receding Galerkin method.

Remark 3. The existence of constant output disturbance does not affect the estimation of unmeasured state . This is because the error in will disappear at the coming sampling time instant once the constant output disturbance appears. In practice, the mathematic model of the unit is not accurate, which means (34) cannot estimate perfectly. In this case, the error cannot be zero but be bounded. In this way, the estimate error will just approach zero.

Remark 4. For either the unit (18), (19), (20), (21), (22), and (23) or the one with constant output disturbances like (18), (19), (20), (21), (22), (23), (24), (25), (27), (28), (29), (30), (31), (32), and (33), the solutions of problem (27) do in fact exist by selecting appropriate feasibility tolerance and the order of approximation , as remarked in [21, 23]. Precise bounds for can be found experimentally using a recursive refinement process through increasing the order of approximation until all the constraints in a nonlinear programming problem like (27) are satisfied.

4. Simulation

In this section, some simulations are conducted to validate the performance of receding Galerkin optimal control strategy for the oil-fired drum-type boiler-turbine unit.

To implement the receding Galerkin optimal control strategy, a nonlinear programming solver SNOPT is adopted and some controller parameters should be preset in advance such as , and . Suppose the sampling time for the unit is . Due to the fact that only the second component in the optimal control sequence should be applied on the unit, the first two approximate time points in -time domain should satisfy . Therefore, with a given node number , according to (6), we can determine the terminal time as where at the th sampling time instant, and can be calculated by (7). During our simulations, we suggest

We choose , in what follows and will discuss the influences of , and on the closed-loop system in Section 4.1. In Section 4.2, we continuously validate the controller performance through several study cases.

4.1. Influences of Controller Parameters

For simplicity in this section, the states, outputs, and control inputs are initialized for all study cases as follows: , , , and . The outputs aim to track constant references ,, and , respectively.

Firstly, we focus on the influence of the node number on the performance of the closed-loop system by considering . Given , where indicates a diagonal matrix. The simulation results are shown in Figure 5, from which we can see that too small value of (e.g., here) usually leads to unappropriate performance and that a larger value of leads to faster responses as well as control inputs. While raises up to a certain limit, the responses and control inputs will change slightly but with an increasing heavy burden of computation. In what follows, we prefer to select . Correspondingly, the length of prediction horizon can be calculated as

With fixing , , and , we can now test the influence of expansion coefficients on the performance of the closed-loop system. Here, we just consider identical by taking values in for all . Simulation results are illustrated in Figure 6, which indicates that larger usually result in faster responses with drastic changes of control inputs. To guarantee the constraints imposed on the changes of controls, we prefer smaller , that is, we suggest for all .

Besides and , matrixes and in the performance index also play a critical role. To see an insight, by fixing and , we consider the following three cases: (1) , (2) , and (3) .

The results are shown in Figure 7. It can be seen from Figure 7 that the responses of the unit will arrive at their static setpoints faster by increasing or relatively decreasing and that smaller control laws can be achieved when taking a larger value of . One can select appropriate and by considering the trade-off between the constraints on the output and control inputs.

From above discussions, we select the controller parameters as , , , and for the consequent performance validations in what follows.

4.2. Performance Validation

This subsection presents the following three different study cases so as to further validate the performance of the receding Galerkin optimal controller.

Case 1. Receding Galerkin optimal control versus PID.
In this case, we compare the performances of receding Galerkin optimal control strategy and PID. We aim to drive the outputs to from initial condition , , , and . The 2-freedom PID controller saturated in bound is designed as with the following parameters

The results are shown in Figures 8, which demonstrates that the receding Galerkin method outperforms the 2-freedom PID controller. As can be seen from Figure 8, the outputs of the receding Galerkin method are much more smooth than that of PID and the control inputs’ rate of change can be well guaranteed by the receding Galerkin method. Note that the receding Galerkin method can guarantee , whereas the PID does not. Furthermore, it is difficult to guarantee the constraints imposed on the control inputs when applying PID.

Case 2. Wide-range load tracking.
In the second case, we intend to validate the tracking performance of the receding Galerkin method when the unit experiences in a wide-range load change. More precisely, we suggest here that the electrical output tries to track load demand from a static condition with a rate of and then back to with the same rate at time . Correspondingly, the drum steam pressure (kg/cm2) rises from 110 (kg/cm2) to 120 (kg/cm2) with a rate of (kg/cm2/s) and descents back to 100 (kg/cm2) with a rate of (kg/cm2). The drum water level has to vary in range .

The simulation results are shown in Figure 9. It shows that the electrical output can track the load demands/references well and the drum water level satisfies output constraint . As well, the drum steam pressure can also track the predefined reference well during the changes of load demands.

Case 3. Eliminating constant output disturbances using an independent model.
In the last case, we intend to show the performance of the receding Galerkin method with an independent model when constant disturbances exist in output channels. On the basis of Case 1, we suggest now that some constant disturbances appear in output channels at different times. More precisely, we define such that

Figure 10 shows that the constant output disturbances can be eliminated and the outputs of unit can track back to their original reference points. However, the control inputs settle at their new steady points. It suggests that the independent model strategy is effective to estimate and thus eliminate the constant output disturbances.

Remark 5. From above simulations and discussions, we can see that optimal solutions can be really found by selecting appropriate tolerance and the order of approximation , as already stated in Remark 4. The curves of the performance index are shown in Figure 11 in the cases of wide-range load tracking and existing constant output disturbances. We can see from Figure 11 that the performance indexes can finally converge to zero and then the unit arrived at static setpoints. This fact may suggest the stability of the closed-loop systems in an intuitive way.

Remark 6. The proposed receding Galerkin optimal control method can be in fact a general approach for a wide range of nonlinear systems rather than only for the boiler-turbine unit. For any nonlinear system, one just needs to design a state observer to estimate the unmeasurable states for this nonlinear system and then embed it into the receding Galerkin optimal controller, as that has been done for the boiler-turbine unit in this paper.

5. Conclusions

In summary, this paper proposes an adaptively receding Galerkin optimal control strategy for a nonlinear boiler-turbine unit. To deal with the problem of unmeasured state variable fluid density, a state observer is designed and embedded into the receding Galerkin strategy. Meanwhile, an independent model structure is constructed in order to estimate and thus eliminate constant disturbances in output channels. Simulation results suggest that the unit can track load reference during wide-range operations with satisfactory performance via this receding Galerkin optimal control strategy.

There are still some further interests. One most possible interest is to extend the proposed method to deal with lumped disturbances.

Data Availability

All the data are available in our article, and the readers who are interested in these data can contact the authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work is supported by Natural Science Foundation of China (51676034 and 51706093) and in part by the Key Project of Yunnan Power Grid Co. Ltd. (YNYJ2016043).