Abstract

Reliable power and temperature control in pressurized water reactor (PWR) nuclear power plant is necessary to guarantee high efficiency and plant safety. Since the nuclear plants are quite nonlinear, the paper presents nonlinear fuzzy model predictive control (MPC), by incorporating the realistic constraints, to realize the plant optimization. T-S fuzzy modeling on nuclear power plant is utilized to approximate the nonlinear plant, based on which the nonlinear MPC controller is devised via parallel distributed compensation (PDC) scheme in order to solve the nonlinear constraint optimization problem. Improved performance compared to the traditional PID controller for a TMI-type PWR is obtained in the simulation.

1. Introduction

In nuclear power plants, steam energy is produced in the nuclear reactor from continuous fission of the atoms of the fuel. The steam is then used to drive the turbine and generator to produce electricity.

Since burning fossil fuels in thermal power plant can have severe environment problem, the development of nuclear power plants should be encouraged. For example in China, although the supercritical and ultrasupercritical technique, as well as integrated gasification combined cycle (IGCC) plants, have been well developed, the coal consumption still reaches about 3.91 million tons in 2012, which comprised over 70% of national primary energy. Fine particles in the air measuring less than 2.5 micrometres reached 993 micrograms per cubic meter in Beijing on 12 January 2013, compared with the limit of 25 published by the WHO guidelines.

Since nuclear power plants are the complex and nonlinear systems, it is a great challenge to control the power and temperature of the nuclear reactor, especially when wide-range power variations occur in the load following condition. Thus, various advanced control schemes have appeared during the past two decades, for example, the observer-based optimal state feedback assisted control [1, 2], the linear quadratic Gaussian with loop transfer recovery control [3, 4], the nonlinear control approach [5, 6], and the neural network/fuzzy approaches [7, 8].

Modern nuclear power plants should respond to the load demand on the power grid, which demand high plant operation performance, subject to various kinds of constraints. Meanwhile, nuclear safety and radioactive pollution prevention have long been much concerned problem. It is therefore extremely important to reach the economical and safe operation to maximize the thermal efficiency of the nuclear power plant. Under these aims, model predictive control (MPC) shows its obvious advantage, since it is an advanced model-based control scheme. It performs an optimization procedure to calculate optimal control actions at every sampling period based on an explicit process model subject to process input, output, and state constraints. So far, MPC has been well constituted for thermal power plant control [911] and also in nuclear power plant water-level control [1214]. Paper [15] reports a MPC application in a three-dimensional nuclear reactor analysis code for nuclear reactor power.

Nuclear power plants are generally nonlinear, due to the frequent changes of the operating point right across the whole operation range. In general, the nonlinear model predictive control (NMPC) also online solves an optimization problem, by using the sequential quadratic program (SQP). The resulting nonlinear programming problems are usually nonconvex, and the online computational burden is generally large. Since the nonlinearities vary with operating power levels, as well as ageing effects [5, 16], fuzzy models are used to express the plant dynamic. MPCs are realized in a parallel distributed compensation (PDC) scheme; that is, for each local model, a MPC is designed. The overall controller, which is nonlinear in nature, is similarly constructed by combining local controllers via fuzzy inference. Simulations on the three-mile island- (TMI-) type pressurized water reactor (PWR) show the effectiveness of the proposed method.

2. The Plant Description

As shown in Figure 1, nuclear fission reaction in the reactor releases tremendous energy, which is transferred from the reactor to the steam generator by the coolant, heating the water in the steam generator to generate abundant steam. Then the turbine is driven by steam to generate electricity. The main objective is to control the power and temperature of the reactor by inserting or elevating the control rods. Conventionally, reactor control system controls the average temperature of the reactor core coolant to track the reference temperature which is proportional to the turbine load, in order to guarantee well matching between the reactor power level and the load demand at the same time. This temperature control strategy, considering not only high thermal efficiency but also technical limits and economical and safe operation, is widely used in modern PWR power plants.

However, the plant is complex and highly nonlinear whose parameters vary with the operation conditions. Meanwhile, realistic constraints in the system may lead to actuator saturation. The performance of the conventional controllers is often unsatisfactory. Besides, unpredictable disturbances in the temperature measurement system will induce fluctuations on the measuring values of coolant temperature, which may lead to the abnormal control rod action and thus must be considered in the controller design. The proposed Takagi-Sugeno fuzzy modeling method is used to approximate the nonlinear plant, based on which the nonlinear MPC controller is devised via PDC scheme in order to solve the nonlinear and constraint problems in the reactor power and temperature control system.

3. T-S Fuzzy Modeling

The reactor model considered is the point kinetics with one delayed neutron group and also coolant and fuel temperature feedback. It is typical for a TMI-type PWR at the middle of the fuel cycle rated at 2500 MW [2, 7]:

The symbols in the above equations are demonstrated in the Nomenclature section.

Define where the symbol indicates a deviation about an equilibrium point. , , , , and correspond to the values of , , , , and at an equilibrium point, respectively.

Choose the state, input, and output vector, respectively, as follows:

The model can be linearized at the equilibrium point using the perturbation theory which is valid only for . The linear model can be written as the following state space form: The corresponding matrices , , and , are where is the power level and is the control rod worth for the th operating point which is defined in Table 1. The linear model (4) is defined at each operating point and the parameters in (6) are functions of the operating point [7, 17]. The remaining constant parameters are shown in Table 2. Consider the following:

The sampling period is selected to be 1 s, and then nine discrete-time local models of (4) can be described as follows: where

According to Table 1, power level and control rod worth are selected as the premise variables. and are defined in region and region [0, 0.029], respectively. and are divided into three fuzzy sets, denoted, respectively, as “low (L),” “medium (M),” and “high (H).” The corresponding univariate membership functions of or are defined by triangle membership functions, as shown in Figure 2. Thus, the space is divided into nine fuzzy subspaces, given as follows:

where denotes the th fuzzy inference rule. are premise variables for and ; denotes the fuzzy set partitioned for the th premise variable, which is valid in the th rule. ( ) represents the consequent th linear model deduced at the th equilibrium point in (7).

By using the center of gravity defuzzification method [10, 18], the final output of the nonlinear T-S model with nine fuzzy rules is where

The denotes the univariate membership function of premise variable with respect to the fuzzy set . is the multivariate membership function obtained by the product of the outputs of the two univariate membership functions, as shown in Figure 3.

Four rules are activated each time instant. At point A in Figure 2, these four rules are as follows: : if is and ,then . : if is and is ,then . : if is and is ,then . : if is and is ,then .

The membership functions are shown in Figure 4.

4. The Nonlinear Fuzzy MPC

4.1. The Local Linear MPC

The local linear model (7) can be expressed as follows:

Subtracting (13) from (7) can result in where , , and .

Hence, .

Choosing a new state variable vector can result in

Denote (15) as the following normal form: where , , , and . and are, respectively, the number of the state variables and the output variables.

The predictive control is obtained by minimizing the following cost function:

The control is obtained by minimizing the cost function at each time instant .

Define where and are the weighting matrices for the prediction error and the control, respectively, and are, respectively, the costing horizon for the prediction error and the control, and and are the optimum th step ahead prediction of the system output on data up to time instant and the th step ahead reference trajectory, respectively. is the vector composed of the future control increments which is obtained by computing the derivative of , expressed as follows: where

The optimized steps ahead control is computed, and only the first step ahead control is implemented, using a receding horizon principle [10], giving further, can be obtained by

The input constraint is expressed as follows: can be expressed by as follows:

Define and constraint (23) can be rewritten as follows:

Hence, the constrained optimization problem can be described as which can be solved by quadratic programming (QP) method.

4.2. The Nonlinear Fuzzy MPC

Based on the T-S fuzzy model, the controller is devised via the parallel distributed control (PDC) scheme [19, 20]. In this scheme, the overall controller, which is naturally nonlinear, is a fuzzy combination of each individual linear controller designed based on the corresponding local model. Here, a nonlinear MPC controller, consisting of nine local MPC controllers, is constituted as shown in Figure 5. For each local linear model, a conventional MPC controller is designed independently. Consequently, the final control is the weighted sum of the control obtained from nine local MPC controllers. Consider the following: where is the control of the th local predictive controller and is defined previously in (11). Since fuzzy logics dominate the switching between local controllers, the resulting control can be smooth by properly selecting the weights [10].

5. The Application to the PWR Nuclear Power Plant

In using the proposed fuzzy MPC, the reference temperature function of power level is defined as follows: Figure 6 shows this relationship.

Assume that the maximum speed of control rod is 0.2 fraction of core length per second. Hence, the constraint on the control is considered as which is involved in the optimization algorithm of MPC.

Define and , where and are the two coefficients. Assume the control rod worth to be 0.0145, the parameters of the nine local controllers are listed in Table 3.

The performance of the nonlinear fuzzy MPC can be evaluated in three cases.

Case A. Full power operation, that is, step changes in power level. Comparison of the performance between the proposed fuzzy MPC and the conventional PID controller is presented in Figures 7(a) and 7(b). The temperature and the power stabilize to their respect set points in an acceptable time, whereas the nonlinear fuzzy MPC reduces the fluctuation obviously. In real-time nuclear power control, inserting the control rod can slow down the release of energy in the reactor so that the temperature and power level are decreased, while lifting the control rod speed will raise the temperature and power level.

Case B. Low power operation, that is, step changes in power level. Comparison of the performance between the proposed controller and the conventional PID controller is presented in Figures 8(a) and 8(b), showing the obvious advantage of the proposed fuzzy MPC. It is noticed that the response time is longer than that in the full power operation and the resulting fluctuation is more severe.
In both cases, the constraint on the control rod speed can be handled by the proposed fuzzy MPC effectively, while the PID controller cannot achieve it, as shown in Figures 7(c) and 8(c).

Case C. Wide-range load following, that is, ramp changes in power level with 5%/min rate, as shown in Figure 9. With the constraint handling on the control rod speed by the fuzzy MPC, satisfactory tracking performance is obtained for both temperature and power.

6. Conclusions

The paper constituted a nonlinear MPC controller by incorporating fuzzy modeling technique. Feasible optimal solutions have been acquired under plant nonlinearity and constraints. The proposed nonlinear fuzzy MPC is simulated in power and temperature control for a three-mile island- (TMI-) type PWR. Better performance is obtained while it is compared with the PID controller. Load following is accomplished under constraints on the control rod.

Nomenclature

:Neutron density relative to density at rated condition
:Precursor density relative to density at rated condition
:Reactivity
:Reactivity due to the control rod
:Fraction of delayed fission neutrons
:Effective prompt neutron lifetime
:Effective precursor density
:Average reactor fuel temperature
:Fraction of reactor power deposited in fuel
:Heat capacity of the fuel
:Rated power level
:Heat transfer coefficient between fuel and coolant
:The temperature of coolant leaving the reactor
:The temperature of coolant entering the reactor
:Heat capacity of the coolant
:Mass flow rate times heat capacity of the water
:Total reactivity worth of the rod
:The control rod speed (fraction of core length per second)
:Fuel temperature reactivity coefficient
:Coolant temperature reactivity coefficient.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by National Natural Science Foundation of China under Grants nos. 60974051 and 61273144 and by Natural Science Foundation of Beijing under Grant no. 4122071.