Research on Pressurizer Pressure Control Based on Adaptive Prediction Algorithm

According to the high control quality requirements of nuclear power plants and the features of the pressurizer pressure with large inertia, time-varying, nonlinear, multi-interference, difficulty in obtaining accurate mathematical model, and open-loop unstable dynamic characteristic, the advanced control strategy is needed for pressurizer pressure control performance optimization. To tackle the problem, an adaptive predictive control method for pressurizer pressure is devised in this paper. Firstly, the nonself-regulating system is stabilized and the adaptive dynamic matrix controller is designed by identifying the controlled object online. In order to realize the engineering application for this controller, then the control signal output is obtained. Finally, the control system simulation platform is built. Simulation results reveal a superior control performance, disturbance rejection, and adaptability. Furthermore, it provides a solution for the application of dynamic matrix control algorithm in non-self-regulating system.


Introduction
Load changes or core reactivity disturbances may cause pressure changes in the primary circuit of a nuclear power plant.If the primary circuit pressure is too high, it may lead to equipment fatigue and pipeline rupture.If the pressure is too low, the risk of melting the fuel element may increase [1].A pressurizer is a major component to control a system pressure in primary circuit for a nuclear power plant.Most of the existing pressurizer pressure control systems adopt the PID control strategy.Because of the high control quality requirements of nuclear power plants and the nonlinear characteristics of pressurizer pressure under different operating conditions, PID sometimes does not guarantee good control effect.Therefore, many studies have been done by scholars to optimize the pressurizer pressure control effect.
An internal model PID control system is applied to the pressurizer pressure control, which shows superior control performance compared with PID [2].A fuzzy PID joint control system using logic judgment and switch shifting [3] is designed, which shows good control results.
However, the selection of switching thresholds and the establishment of fuzzy controller came from experience, which is not conducive to learning and promotion.On the other hand, the pressurizer pressure has the features of large inertia, time-varying, nonlinear, multi-interference, difficulty in obtaining accurate mathematical model, and open-loop unstable dynamic characteristic.All these features lead to problem of choosing an appropriate method for control performance optimization.
Predictive control algorithm has lower requirements on model and has characteristics of rolling optimization and feedback correction.Dynamic Matrix Control (DMC) has been widely used in industrial process control, but it will produce truncation error when modeling non-self-regulating objects.The DMC algorithm is improved based on the objects' step response characteristics of an approximate straight line in the final stage [4].A stable generalized predictive controller (SGPC) is designed by decomposing the non-self-regulating model [5], which theoretically solves the non-self-regulating system control problem.Rossiter et al. and Rossiter and Kouvaritakis [6,7] make use of the particularity of SGPC structure to carry on massive research to the constrained predictive control, which obtains better stability results.Constrained predictive control algorithm can improve the dynamic performance and the disturbance rejection ability to large-inertia, large-delay controlled objects, but it has poor adaptability to time-varying and nonlinear controlled objects.
Adaptive control can adapt to changes in the dynamic characteristics of objects and disturbances.Adaptive control and fuzzy control are combined and applied to the pressurizer pressure regulation of nuclear power plant with pressurized water reactor (PWR) [8] and vessel nuclear power plant [9], which can improve the robustness, rejection of disturbances, and adaptability of model changes.Model-free adaptive control is used to optimize the steam generator water level in nuclear power plant [10].An adaptive output feedback structure is designed based on uncertain nonlinear system with time-varying delay [11].Comparing the nonlinear model predictive control and the generalized predictive control [12], it shows that the former has a better control effect, but it requires more calculation time.In order to optimize operation speed, the offline calculation is converted to online optimization and the model predictive control (MPC) is improved [13].However, few of the above methods have been applied to a pressurizer pressure control system, and few of them mentioned in the collection of engineering control signal output.
Therefore, an advanced control strategy for non-selfregulating system and its practical value of engineering are two main problems of pressurizer pressure control performance optimization.In this paper, the pressurizer pressure adaptive predictive control method based on DMC algorithm is designed.Firstly, a feedback structure is used to selfstabilize the non-self-regulating system and the adaptive predictive controller is designed by identifying the controlled object online.Then, the control signal output is solved to make this controller easy to be implemented in engineering.Finally, the feasibility of the controller design is verified by simulation comparison and disturbance test.
In Section 2, the pressurizer of a nuclear power plant and the equivalent nonparametric model of pressurizer pressure are illustrated.Section 3 presents the control method and the adaptive predictive controller.Section 4 analyzes the performance of the proposed control method.Finally, conclusions are made in Section 5. Figure 1 shows an article structure layout of the main content of this paper.

Description of the Pressurizer
2.1.Pressurizer in a Nuclear Power Plant.The pressurizer is an important equipment in primary circuit for a nuclear power plant.The basic functions of the pressurizer are pressure control, pressure protection, and compensation of primary circuit coolant volume change.According to the different structure and operating principles, the pressurizer can be divided into gas tank pressurizer and electrothermal pressurizer.The structure of the gas tank pressurizer is simple, but there are certain nuclear safety issues in the process of compressing air or high-pressure inert gas.Therefore, the electrothermal pressurizer is often used in modern nuclear power plants.
The structure of the electrothermal pressurizer is shown in Figure 2. The upper part is the steam space, the lower part is the water space, and the bottom is connected to the primary circuit heat pipe section through the surge line.The conventional PWRs adopt a saturated steam pressure regulation method to achieve pressure control by electric heating and spraying.More specifically, the pressurizer is made up of the spray system, the electric heater unit, the safety valve group, and the measuring instrument.The spray system, which is primarily composed of spray valves, is used to spray the coolant with a low temperature to reduce the pressure of primary circuit.The electric heater unit, which consists of electric heating rods, is used to heat the coolant in pressurizer to raise the pressure of primary circuit.
Given the feedback channel gain K, the structure of the controlled object with feedback is shown in Figure 3.
The model in Figure 3 can be expressed as follows: The stability of G 1 s in ( 2) is related to the value of K.When K is smaller than the reciprocal of the product of k o and τ, G 1 s is stable [14].It can be seen that the pure delay link in the denominator is too complex to be directly used as a predictive controller internal model and needs to be simplified.According to the differential deflection linearization method of a univariate nonlinear system, the pure delay link in denominator can be expanded to the Taylor series near the steady-state point, and the high-order term is omitted.The simplification is shown in e −τs ≈ 1 − τs + τ 2 s 2  3 Substituting (3) into (2) results in an equivalent model.
If G 1 s is stable, the equivalent model G 2 s in (4) can be considered as a self-regulating controlled object to design the initialization parameters of dynamic matrix controller.Taking the opening of spray valve to adjust the pressurizer pressure as an example, this paper obtains experimental data from the nuclear power plant simulator and identifies the transfer function G s in (1) under different operating conditions.The full scope simulator has a 1 : 1 fidelity with the nuclear power plant; hence, it is equivalent to the actual data source of the nuclear power plant.G s that shows the dynamic characteristics between the spray valve opening and the pressurizer pressure is shown in Table 1.
The equivalent nonparametric model of DMC is easily obtained from G 2 s .Assuming the equivalent nonparametric model is where A 1-80 is a 100-dimensional column vector and the subscript of A represents A 1 under 80% operating condition.3 Complexity object consists of feedback channel, actuator, and pressurizer controlled object.Then, the control signal output u C k of adaptive predictive controller is obtained by (8).And the pressurizer pressure adaptive predictive controller that directly acts on non-self-regulating system is established.

Adaptive Predictive Control Method of Pressurizer Pressure
where k represents the time of simulation, u k is the controller output at step k, K is the feedback channel gain, and y k is the pressure output at step k.
The adaptive predictive controller is composed of online identification algorithm module, DMC controller module, and feedback channel module.The online identification algorithm module is an online identification layer.The DMC controller module and the feedback channel module are control layer.The online identification layer and the control layer work in parallel.The online identification algorithm uses the current and the historical values of control increment Δu k and pressure output y k to calculate the dynamic matrix A of DMC controller online, which achieves the purpose of self-adaptation of the controller.A consists of the element values in A 1 , which can be obtained by (9).The DMC algorithm is a predictive control algorithm based on device step response and applicable to progressively stable linear device [15].Its predictive model is easy to obtain from engineering, and it has less computations and strong robustness.The online cyclic operation with prediction, correction, and control achieves prediction and optimization.
where ΔU 2 k consists of the element values in ΔU 1 k to form an N-dimensional matrix, ΔU 1 k is a sequence of N control incremental historical values, Y 2 k consists of the element values in Y 1 k to form an N-dimensional column vector, and Y 1 k is a sequence of N + 1 pressure output historical values.(11).
Equation ( 11) can be rewritten as The matrix form of ( 12) can be written as The nonparametric model A 1 can be deduced in (13).5 Complexity Equation ( 14) can also be expressed as (9), and the elements in A 1 are shown as follows: It is easy to get a new dynamic matrix A at this moment from the nonparametric model A 1 of the controlled object.Then, the original A of prediction algorithm is replaced to complete the online model identification.

Self-Regulating Controlled Object Prediction Model.
The prediction time domain P and the control time domain M are selected, usually M ≤ P ≤ N. Assume that the initial predictive value of pressure at time k is ŷ0 k + i | k , i = 1, 2, … , N, where the subscript of y represents the number of control action changes and k + i|k represents the prediction of time k + i at time k.When the opening of the spray valve has an increment Δu k , the pressure predictor at the next moment is where a i is the value of the element in the nonparametric model A , i 1, 2, … , N. When M continuous spray valve opening changes Δu k , … , Δu k + M − 1 occur, the predicted value of pressurizer pressure is 17) is the predictive model of the output under the effect of continuous control increment.

Scroll Optimization and Solution of Control Signal
Output.M changes of spray valve opening from each time k are needed, which can make the pressure predictor value ŷM k + i | k for the next P moments as close as possible to the expected pressure value ω k + i , i = 1, 2, … , P. The drastic change of Δu k is usually undesirable during control.So a soft constraint can be added to the performance optimization index.

18
where q i is the error weighting coefficient and r i is the control weighting coefficient, which, respectively, represents the restraint on tracking error and control variation.
The vector form of ( 17) can be written as where A is called the dynamic matrix.
The vector form of (18) can be written as , the error weight matrix Q = diag q 1 , ⋯, q P , and R = diag r 1 , ⋯, r M is the control weight matrix.Substituting (19) into (21), The optimal control increment ΔU M k in ( 22) can be deduced by the extremum requirement dJ k /dΔU M k = 0, The "rolling optimization" strategy means that the DMC only takes the immediate control increment Δu k of ΔU M k to form a controller output u k = u k − 1 + Δu k , and the next time, the same method is used to obtain the instantaneous control increment Δu k + 1 to complete a rolling optimization.
Since u k is the output acting on the self-regulating model, which theoretically realizes the non-self-regulating system predictive control.In order to solve the engineering problem, the control signal output C k directly acting on spray device is required by (8).Equation ( 8) represents the structural relationship of self-regulating process of controlled object.Then, the structural relationship can be expressed as where Δu k = d T ω P k − ŷP0 k , the control vector d T = C T A T QA + R −1 A T Q, and the M-dimensional column vector C T = 1, 0, ⋯, 0 which means to get the first element, i = 1, 2, … , M.

Feedback Correction.
Model mismatch, environmental disturbances (such as load changes or core reactivity disturbances) and other factors, may cause the output error The weighting method can be used to correct the error by introducing the correction vector h cor = h 1 , ⋯, h N T The corrected output prediction vector is shown in where The initial prediction value at time k + 1 is obtained by shifting, which can be expressed in vector form as follows: where S is the shift matrix

Simulation
The simulation test platform is based on Visual Studio2010 and MATLAB/Simulink.The function of adaptive DMC controller of pressurizer pressure control system in nuclear power plant is realized in Visual Studio2010.It is connected with MATLAB/Simulink through OLE for Process Control (OPC) interface.MATLAB/Simulink simulation model is shown in Figure 6.The controller output calculated by adaptive DMC controller is transmitted to the OPC Read module through OPC interface, and the control signal output acts on the controlled object.The pressure output will be returned by the OPC Write module to adaptive DMC controller for controller output calculation via the OPC interface.Finally, the tracking of pressure output value to set value is realized.
Compared with DMC and PID, the improvement of control performance is verified by the proposed non-self-regulating system adaptive predictive control method.
4.1.Tracking Performance of Control System.The feedback channel gain K = 1 meets the stability of G 1 s in (2).The disturbance is set to zero, which temporarily ignores the effects of disturbances.The initialization parameters of pressurizer pressure dynamic matrix controller are designed by A 1-90 under 90% operating condition (as shown in Table 1).The set point of the pressure steps from 0 to 1 MPa at t = 0 s, and the control characteristics of pressure adaptive predictive control, pressure dynamic matrix control, and pressure PID control are shown in Figure 7. Figure 7 shows that the rise time of ordinary DMC is 2.8 s, and the adjustment time of 7 Complexity DMC is 5 s (error is ±2%).It is observed that the rise time of adaptive prediction control is 2.3 s and the adjustment time of adaptive prediction control is 3.8 s (error is ±2%), which is, respectively, 18% and 24% lower than DMC, and the integral absolute error (IAE) is reduced by 13% compared with DMC.PID control response curve can also achieve an approximate rise time at the cost of overshoot adjustment time.The rise time of PID is 2.7 s, the peak time of PID is 7.2 s, the overshoot of PID is 15.6%, and the adjustment time of PID is 18.2 s (error is ±2%).The results show that the adaptive prediction control has faster response and better static stability.The adaptive predictive control signal output curve directly acting on spray device and the controller output curve are shown in Figure 8 (from 0 to 10 s in an interval of 30 s).It is obvious that this pressurizer pressure adaptive prediction controller can be realized in engineering, and it can also be applied to solve the kind of non-self-regulating system predictive control problems.

Responses with Disturbance.
Disturbance is inevitable in an industrial process.Load changes and core reactivity disturbances are ultimately reflected in changes of pressure.It is necessary to test the disturbance rejection ability of the adaptive predictive controller with a pressure disturbance.The feedback channel gain K is set to 1.The initialization parameters of pressurizer pressure dynamic matrix controller are designed by A 1-90 under 90% operating condition (as shown in Table 1).In an interval of 50 s, a step pressure disturbance of 0.4 MPa is added at t = 20 s after the system is stable.The response curves with different control systems under disturbance are shown in Figure 9.It is observed that the adjustment time of ordinary DMC is 5.5 s (error is ±2%) and the adjustment time of adaptive predictive control is 4.8 s (error is ±2%).The latter is 13% shorter than the former.And the adjustment time of PID is 15 s (error is ±2%), which is about three times than that of adaptive predictive control and DMC.The results show that the disturbance rejection ability of adaptive predictive controller is better than that of ordinary dynamic matrix controller and PID controller, and it can achieve a faster and more stable control effect.

Robustness of Control System.
In order to further verify the adaptability of adaptive predictive controller to the change of controlled object model, take 90% operating condition (as shown in Table 1) as initial condition, 100% and 80% (as shown in Table 1) as comparison conditions to compare the adaptive predictive control effect with DMC and PID.
The feedback channel gain K is set to 1, and the disturbance is set to zero.The initialization parameters of pressurizer pressure dynamic matrix controller are designed by A 1-90 under 90% operating condition (as shown in Table 1).The operating condition is changed from 90% to 100% at t = 0 s, and the characteristic curves of control systems are shown in Figure 10.It is observed that the rise time of DMC is 2.5 s, the peak time of DMC is 5.3 s, the overshoot of DMC is 9.3%, and the adjustment time of DMC is 8 s (error 8 Complexity is ±2%).The rise time of adaptive predictive control is 3.5 s, the adjustment time of adaptive predictive control is 6 s (error is ±2%) which is 25% shorter than DMC, and the IAE is reduced by 2% compared with DMC.The rise time of PID is 4.7 s, the peak time of PID is 12.2 s, the overshoot of PID is 25%, and the adjustment time of PID is 24.8 s (error is ±2%).It is clear that the adaptive predictive control characteristic curve changes more smoothly and there is no overshoot.The operating condition is changed from 90% to 80% at t = 0 s, and the characteristic curves of control systems are shown in Figure 11.It shows that the rise time of ordinary DMC is 3 s, the adjustment time of DMC is 5.3 s (error is ±2%), the rise time of adaptive prediction control is 2 s which is 33% shorter than DMC, the adjustment time of adaptive prediction control is 3.3 s (error is ±2%) which is 38% shorter than DMC, and the IAE is reduced by 21% compared with DMC.Similar to the performance in Figure 10, the adaptability of PID controller to the change of controlled object model is relatively poor.It can be seen that adopting adaptive predictive controller not only makes the control characteristic curve more stable but also has better tracking ability and adaptability.

Conclusion
Stimulated by the control performance optimization of pressurizer pressure, the adaptive predictive control method based on DMC algorithm is proposed in this paper.Some conclusions could be drawn as follows.
(1) The pressurizer pressure adaptive predictive controller has good control performance and dynamic characteristic.It is characterized by rapid response, small overshoot, strong disturbance rejection ability, and guaranteed robustness.
(2) The stabilization of non-self-regulating system overcomes the disadvantages of DMC in modeling non-self-regulating object.
(3) The acquisition of control signal output shows the adaptive predictive controller is effective and practical value of engineering application.
(4) The online identification of models helps to improve the ability of predictive controller to adapt to changes in controlled objects.
(5) The physical constraints of valve have been ignored for the time being, and further research and improvement are needed.

10 3. 2 .
Online Identification of Adaptive Predictive Controller 3.2.1.Online Identification Algorithm Structure.The online identification algorithm continuously detects the instantaneous Δu k and the current y k and updates the vectors ΔU 1 k , Y 1 k , and Y 2 k and the matrix ΔU 2 k in real time.Then, A 1 of the controlled object is calculated, and finally, the dynamic matrix A is obtained for adaptive predictive control, which completes online identification of the controlled object model.The algorithm flowchart is shown in Figure 5.At each time k, the control increments of first k − 1 times in ΔU 1 k are shifted up by one, then the online identification module detects the current Δu k to replace the original

Figure 6 :
Figure 6: Simulation model based on pressurizer pressure adaptive predictive control.

Table 1 :
Mathematical model of the controlled object. .And the new vector ΔU 1 k at that moment is obtained that can be used to construct a new matrix ΔU 2 k .In the same way, the new vectors Y 1 k and Y 2 k at this moment are obtained.Multiply ΔU 2 k and Y 2 k to get A 1 , then A is obtained.
3.2.2.Online Solution of Dynamic Matrix.In every moment, the relationships between control increment, pressure output, and nonparametric model elements are shown in