# Patent

## Title
Model-free adaptive control of advanced power plants

## Abstract

A novel 3-Input-3-Output (3×3) Model-Free Adaptive (MFA) controller with a set of artificial neural networks as part of the controller is introduced. A 3×3 MFA control system using the inventive 3×3 MFA controller is described to control key process variables including Power, Steam Throttle Pressure, and Steam Temperature of boiler-turbine-generator (BTG) units in conventional and advanced power plants. Those advanced power plants may comprise Once-Through Supercritical (OTSC) Boilers, Circulating Fluidized-Bed (CFB) Boilers, and Once-Through Supercritical Circulating Fluidized-Bed (OTSC CFB) Boilers.

## Background
This application claims priority to U.S. Provisional Application No. 61/473,308 filed on Apr. 8, 2011, which is herein incorporated by reference.


This invention was made with government support under SBIR grant DE-FG02-06ER84599 awarded by the U.S. Department of Energy. The government has certain rights to the invention.


INVENTION
The subject of this patent relates to automatic control of power plants, and more particularly to a method and apparatus for intelligently controlling the boiler-turbine-generator (BTG) units of power plants.
For the U.S. to reach its future energy objectives, visions to build ultra-clean and highly efficient energy plants of the future have to be realized. In parallel with the development of sensors, more robust and flexible process control technologies must be developed to build an intelligent control system that can yield a fully automated operation and be adaptive to changing process needs and fuel availability. It must be safe, reliable, and easy to install, maintain, and operate. The intelligent control system is aimed to control conventional boilers as well as advanced boilers including Once-through Supercritical Boilers, Circulating Fluidized-bed (CFB) Boilers, and Supercritical CFB Boilers in future energy plants that can deliver maximum-energy-efficiency, near-zero-emissions, fuel-flexibility, and multi-products.
First introduced in 1997, the Model-Free Adaptive (MFA) control technology overcomes the shortcomings of traditional Proportional-Integral-Derivative (PID) controllers and is able to control various complex processes that may have one or more of the following behaviors: (1) nonlinear, (2) time-varying, (3) large time delay, (4) multi-input-multi-output, (5) frequent dynamic changes, (6) open-loop oscillating, (7) pH process, and (8) processes with large load changes and disturbances.
Since MFA is “Model-Free”, it also overcomes the shortcomings of model-based advanced control methods. MFA is an adaptive and robust control technology but it does not require (1) precise process models, (2) process identification, (3) controller design, and (4) complicated manual tuning of controller parameters. A series of U.S. patents and related international patents for Model-Free Adaptive (MFA) control and optimization technologies have been issued. Some of them are listed in Table 1.








TABLE 1





U.S. Pat. No.
Patent Name







6,055,524
Model-Free Adaptive Process Control


6,556,980
Model-Free Adaptive Control for Industrial Processes


6,360,131
Model-Free Adaptive Control for Flexible Production



Systems


6,684,115
Model-Free Adaptive Control of Quality Variables (1)


6,684,112
Robust Model-Free Adaptive Control


7,016,743
Model-Free Adaptive Control of Quality Variables (2)


7,142,626
Apparatus and Method of Controlling Multi-Input-Single-



Output Systems


7,152,052
Apparatus and Method of Controlling Single-Input-Multi-



Output Systems


7,415,446
Model-Free Adaptive Optimization









Commercial hardware and software products with Model-Free Adaptive control have been successfully installed in most industries and deployed on a large scale for process control, building control, and equipment control.
In this patent, we introduce a novel 3-Input-3-Output (3×3) Model-Free Adaptive (MFA) controller, which enhances the capabilities of the MIMO MFA controllers described in U.S. Pat. Nos. 6,055,524 and 6,556,980. We will also describe a 3×3 MFA control system using the inventive 3×3 MFA controller to control key process variables including Power, Steam Throttle Pressure, and Steam Temperature of boiler-turbine-generator (BTG) units in conventional and advanced power plants.



In the accompanying drawing:
FIG. 1 is a block diagram illustrating a single variable Model-Free Adaptive (MFA) control system.
FIG. 2 is a block diagram illustrating the architecture of a single-input-single-output (SISO) Model-Free Adaptive (MFA) controller in its basic form.
FIG. 3 is a block diagram illustrating a 3-Input-3-Output (3×3) multivariable Model-Free Adaptive (MFA) control system according to this invention.
FIG. 4 is a schematic representation of a boiler-turbine-generator (BTG) unit of a power plant comprising a conventional boiler.
FIG. 5 is a schematic representation of a boiler-turbine-generator (BTG) unit of a power plant comprising a supercritical circulating fluidized-bed boiler.
FIG. 6 is a diagram illustrating the key process variables of the boiler-turbine-generator (BTG) unit of a power plant.
FIG. 7 illustrates the block diagram of a 3×3 MFA control system for controlling the 3×3 Power-Pressure-Temperature (PPT) process of a Boiler-Turbine-Generator (BTG) unit according to this invention.
FIG. 8 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoint for Power is stepped down.
FIG. 9 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoint for Steam Throttle Pressure is stepped down.
FIG. 10 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoint for Steam Temperature is stepped down.
FIG. 11 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoints for Power, Throttle Pressure, and Steam Temperature all have step changes.



In this patent, the term “mechanism” is used to represent hardware, software, or any combination thereof The term “process” is used to represent a physical system or process with inputs and outputs that have dynamic relationships.
Without losing generality, all numerical values given in controller parameters in this patent are examples. Other values can be used without departing from the spirit or scope of our invention.
DESCRIPTION
A. Single-Variable MFA Control System
FIG. 1 illustrates a single variable Model-Free Adaptive (MFA) control system as described in the U.S. Pat. No. 6,556,980 B1 comprising a single-input-single-output (SISO) MFA controller 10, a SISO process 12, and signal adders, 14, 16. The signals shown in FIG. 1 are as follows:
r(t)—Setpoint (SP)
y(t)—Measured Process Variable (PV), y(t)=x(t)+d(t).
x(t)—Process Output
u(t)—Controller Output (OP)
d(t)—Disturbance, the disturbance caused by noise or load changes.
e(t)—Error between the Setpoint and Measured Variable, e(t)=r(t)−y(t).
The control objective for the controller is to produce an output u(t) to force the process variable y(t) to track the given trajectory of its setpoint r(t) under variations of setpoint, disturbances, and process dynamics. In other words, the task of the MFA controller is to minimize the error e(t) in an online fashion. We select the objective function for the MFA control system as




  
    
      
        
          
            
              
                
                  E
                  S
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
              ⁢
              
                
                  1
                  2
                
                ⁢
                
                  
                    e
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                  2
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    
                      1
                      2
                    
                    ⁡
                    
                      [
                      
                        
                          r
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                        -
                        
                          y
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                      
                      ]
                    
                  
                  2
                
                .
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




The minimization of e(t) is achieved by (i) the regulatory control capability of the MFA controller, and (ii) the adjustment of the MFA's weighting factors allowing the controller to deal with process dynamic changes, disturbances, and other uncertainties.

FIG. 2 illustrates the core architecture of a SISO MFA controller in its basic form as described in the U.S. Pat. No. 6,556,980 B1. Used as a key component, a multilayer perceptron neural network consists of one input layer 20, one hidden layer 22 with N neurons, and one output layer 24 with one neuron. The input signal e(t) to the input layer 20 is firstly converted to a normalized error signal E1 with a range of −1 to 1 by using the normalization unit 26, where N(.) denotes a normalization function. The output of the normalization unit 26 is then scaled by a scaling function L(.) 25:




  
    
      
        
          L
          ⁢
          
            (·)
          
        
        =
        
          
            
              K
              c
            
            
              T
              c
            
          
          .
        
      
    
    
      
        (
        2
        )
      
    
  




The value of E1 at time t is computed with function L(.) and N(.):





  
    
      
        
          
            E
            1
          
          =
          
            
              
                K
                c
              
              
                T
                c
              
            
            ⁢
            
              N
              ⁡
              
                (
                
                  e
                  ⁡
                  
                    (
                    t
                    )
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




where Kc>0 is defined as controller gain and Tc is the user selected process time constant. These are important parameters for the MFA controller since Kc is used to compensate for the process steady-state gain and Tc provides information for the dynamic behavior of the process. When the error signal is scaled with these parameters, the controller's behavior can be manipulated by adjusting the parameters.

The use of Tc as part of the scaling function permits a broad choice of sample intervals, Ts, because the only restriction is that Ts must conform to the formula Ts<Tc/3 based on the principles of information theory.
The E1 signal then goes iteratively through a series of delay units 28, where z−1 denotes the unit delay operator. A set of normalized and scaled error signals E2 to EN is then generated. In this way, a continuous signal e(t) is converted to a series of discrete signals, which are used as the inputs to the neural network. These delayed error signals Ei, i=1, . . . N, are then conveyed to the hidden layer through the neural network connections. This is equivalent to adding a feedback structure to the neural network. Then the regular static multilayer neural network becomes a dynamic neural network, which is a key component for the Model-Free Adaptive controller.
Each input signal is conveyed separately to each of the neurons in the hidden layer 22 via a path weighted by an individual weighting factor wij, where i=1, 2, . . . N, and j=1, 2, . . . N. The inputs to each of the neurons in the hidden layer are summed by adder 30 to produce signal pj. Then the signal pj is filtered by an activation function 32 to produce qj, where j denotes the jth neuron in the hidden layer.
A piecewise continuous linear function f(x) mapping real numbers to [0,1] is used as the activation function in the neural network as defined by




  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          0
        
        ,
        
          
            if
            ⁢
            
              
            
            ⁢
            x
          
          <
          
            -
            
              b
              a
            
          
        
      
    
    
      
        (
        
          4
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          
            ax
            +
            b
          
        
        ,
        
          
            if
            ⁢
            
              
            
            -
            
              b
              a
            
          
          ≤
          x
          ≤
          
            b
            a
          
        
      
    
    
      
        (
        
          4
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          1
        
        ,
        
          
            if
            ⁢
            
              
            
            ⁢
            x
          
          >
          
            b
            a
          
        
      
    
    
      
        (
        
          4
          ⁢
          c
        
        )
      
    
  




where a is an arbitrary constant and b=1⁄2.

Each output signal from the hidden layer is conveyed to the single neuron in the output layer 24 via a path weighted by an individual weighting factor hj, where j=1, 2, . . . N. These signals are summed in adder 34 to produce signal z(.), and then filtered by activation function 36 to produce the output o(.) of the neural network 18 with a range of 0 to 1.
A de-normalization function 38 defined by

D(x)=100x,   (5)

maps the o(.) signal back into the real space to produce the controller signal u(t).

The algorithm governing the input-output of the controller consists of the following difference equations:




  
    
      
        
          
            
              p
              j
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          
            
              q
              j
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            f
            ⁡
            
              (
              
                
                  p
                  j
                
                ⁡
                
                  (
                  n
                  )
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          
            
              
                
                  o
                  ⁡
                  
                    (
                    n
                    )
                  
                
                =
                
                ⁢
                
                  f
                  ⁡
                  
                    (
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            
                              j
                              ⁢
                              
                                
                              
                            
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                    )
                  
                
              
              ,
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
              
              ,
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where the variables of function f(.) are in the range specified in Equation (4b) and bounded by the limits specified in Equations (4a) and (4c). The controller signal u(t) becomes





  
    
      
        
          
            
              
                u
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
              ⁢
              
                
                  
                    
                      K
                      c
                    
                    ⁡
                    
                      (
                      .
                      )
                    
                  
                  ⁢
                  
                    e
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  D
                  ⁡
                  
                    (
                    
                      o
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    )
                  
                
              
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    
                      
                        K
                        c
                      
                      ⁡
                      
                        (
                        .
                        )
                      
                    
                    ⁢
                    
                      e
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                  +
                  
                    100
                    ⁡
                    
                      [
                      
                        
                          a
                          ⁢
                          
                            
                              ∑
                              
                                j
                                =
                                1
                              
                              N
                            
                            ⁢
                            
                              
                                
                                  h
                                  j
                                
                                ⁡
                                
                                  (
                                  n
                                  )
                                
                              
                              ⁢
                              
                                
                                  q
                                  j
                                
                                ⁡
                                
                                  (
                                  n
                                  )
                                
                              
                            
                          
                        
                        +
                        b
                      
                      ]
                    
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




where n denotes the nth iteration; o(t) is the continuous function of o(n); D(.) is the de-normalization function; and Kc(.)>0, the controller gain 42, is a parameter used to adjust the magnitude of the control output. This is the same Kc used in the scaling function L(.) 25. It is useful to adjust the controller performance.

An online learning algorithm as described in the U.S. Pat. No. 6,556,980 B1 is an example of one algorithm that can be used to continuously update the values of the weighting factors of the MFA controller as follows:

Δwij(n)=a2ηe(n)Ei(n)hj(n),   (10)

Δhj(n)=aηe(n)qj(n).   (11)

The equations (1) through (11) work for both process direct-acting or reverse acting types. Direct-acting means that an increase in the process input will cause its output to increase, and vice versa. Reverse-acting means that an increase in the process input will cause its output to decrease, and vice versa. To keep the above equations working for both direct and reverse acting cases, e(t) is calculated differently based on the acting type of the process as follows:

















e(t) = r(t) − y(t),
if direct acting
(12a)



e(t) = −[r(t) − y(t)].
if reverse acting
(12b)










This is a general treatment for the process acting types. It applies to all MFA controllers to be introduced below.
The algorithm for updating the weighting factors is based on the goal of minimizing the error e(t). Since this effort is the same as the control objective, the adaptation of the weighting factors can assist the controller to minimize the error while process dynamics are changing. In addition, the neural network based MFA controller “remembers” a portion of the process data providing valuable information for the process dynamics. In comparison, a digital version of the Proportional-Integral-Derivative (PID) controller remembers only the current and previous two samples. In this regard, PID has almost no memory, and MFA possesses the memory that is essential to an “intelligent” controller.
B. 3-Input-3-Output MFA Control System
FIG. 3 is a block diagram illustrating a 3-Input-3-Output (3×3) multivariable Model-Free Adaptive (MFA) control system according to this invention. In the 3×3 MFA control system, the 3×3 MFA controller 52 consists of three main controllers C11, C22, C33, and six compensators C21, C31, C12, C32, C13, C23. The 3×3 process 54 has three main processes G11, G22, G33, and six sub-processes G21, G31, G12, G32, G13, G23. The measured process variables y1, y2, and y3 are used as the feedback signals for each of the main control loops, respectively. They are compared with the setpoints r1, r2, and r3 at adders 56 to produce errors e1, e2, and e3. The output of each main controller associated with one of the inputs e1, e2, or e3 is combined with the output of the compensators associated with the other two inputs by adders 58 to produce control signals u1, u2, and u3.
In the 3×3 MFA control system, the signals shown in FIG. 3 are as follows:
r1(t), r2(t), r3(t)—Setpoint (SP) of controllers C11, C22 and C33, respectively.
y1(t), y2(t), y3(t)—Measured Process Variables (PV) of the 3×3 process.
e1(t), e2(t), e3(t)—Error between the Setpoint and its corresponding Process Variable (PV).
u1(t), u2(t), u3(t)—Outputs of the 3×3 controller used as inputs to the 3×3 process.
d1(t), d2(t), d3(t)—Disturbance to y1, y2 and y3, respectively.
v11(t), v22(t), v33(t)—Output of the main controllers C11, C22 and C33, respectively.
v21(t), v31(t) v12(t), v32(t) v13(t), v23(t)—Output of compensators C21, C31, C12, C32, C13, and C23, respectively.
x11(t), x21(t), x31(t), x12(t), x22(t), x32(t), x13(t), x23(t), x33(t)—Output of processes G11, G21, G31, G12, G22, G32, G13, G23, and G33, respectively.
The process variable (PV) for each main control loop is affected by the process output of 3 sub-processes and a disturbance signal as follows:

y1(t)=x11(t)+x12(t)+x13(t)+d1(t),   (13a)

y2(t)=x21(t)+x22(t)+x23(t)+d2(t),   (13b)

y3(t)=x31(t)+x32(t)+x33(t)+d3(t).   (13c)

For each main control loop, the error between the Setpoint and its corresponding Process Variable (PV) is calculated based on the acting type of the corresponding process as follows:

















e1(t) = r1(t) − y1(t),
if G11 is direct acting
(14a)



e1(t) = −[r1(t) − y1(t)],
if G11 is reverse acting
(14b)



e2(t) = r2(t) − y2(t),
if G22 is direct acting
(14c)



e2(t) = −[r2(t) − y2(t)],
if G22 is reverse acting
(14d)



e3(t) = r3(t) − y3(t),
if G33 is direct acting
(14e)



e3(t) = −[r3(t) − y3(t)],
if G33 is reverse acting.
(14f)










Please note that a process is said to be direct acting, if the process input increases, the process output increases, and vice versa. A process is said to be reverse acting, if the process input increases, the process output decreases, and vice versa.

For each main control loop, the controller output is a combination of the output from the controller and its corresponding compensators as follows:

u1(t)=v11(t)+v12(t)+v13(t),   (15a)

u2(t)=v21(t)+v22(t)+v23(t),   (15b)

u3(t)=v31(t)+v32(t)+v33(t).   (15c)

For each controller C11, C22, and C33, the value of E1 at time t is computed with a normalization function N(.) as follows:




  
    
      
        
          
            
              E
              1
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                11
              
              
                T
                c
                11
              
            
            ⁢
            
              N
              ⁡
              
                (
                
                  
                    e
                    1
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          16
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              E
              1
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                22
              
              
                T
                c
                22
              
            
            ⁢
            
              N
              ⁡
              
                (
                
                  
                    e
                    2
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          16
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              E
              1
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                33
              
              
                T
                c
                33
              
            
            ⁢
            
              N
              ⁡
              
                (
                
                  
                    e
                    3
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          16
          ⁢
          c
        
        )
      
    
  




Where Kc11>0, Kc22>0, and Kc33>0 are the controller gain and Tc11>0, Tc22>0, and Tc33>0 are the time constants for C11, C22, and C33, respectively.
The algorithm governing the input-output of the 3×3 MFA controller 52 as well as the online learning algorithm of the neural network consist of the following difference equations:
For Controller C11:




  
    
      
        
          
            
              p
              j
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                  11
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                  11
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        
          17
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              q
              j
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                ap
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
            +
            b
          
        
        ,
      
    
    
      
        (
        
          17
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              v
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                11
              
              ⁢
              
                
                  e
                  1
                
                ⁡
                
                  (
                  n
                  )
                
              
            
            +
            
              100
              ⁡
              
                [
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                            11
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                            11
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
                ]
              
            
          
        
        ,
      
    
    
      
        (
        
          17
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                w
                ij
                11
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            
              a
              2
            
            ⁢
            
              η
              11
            
            ⁢
            
              
                e
                1
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                E
                i
                11
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                h
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          17
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                h
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              η
              11
            
            ⁢
            
              
                e
                1
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                q
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          17
          ⁢
          e
        
        )
      
    
  




For Controller C22 





  
    
      
        
          
            
              p
              j
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                  22
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                  22
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        
          18
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              q
              j
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                ap
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
            +
            b
          
        
        ,
      
    
    
      
        (
        
          18
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              v
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                22
              
              ⁢
              
                
                  e
                  2
                
                ⁡
                
                  (
                  n
                  )
                
              
            
            +
            
              100
              ⁡
              
                [
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                            22
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                            22
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
                ]
              
            
          
        
        ,
      
    
    
      
        (
        
          18
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                w
                ij
                22
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            
              a
              2
            
            ⁢
            
              η
              22
            
            ⁢
            
              
                e
                2
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                E
                i
                22
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                h
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          18
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            
              h
              j
              22
            
            ⁡
            
              (
              n
              )
            
          
        
        =
        
          a
          ⁢
          
            
          
          ⁢
          
            η
            22
          
          ⁢
          
            
              e
              2
            
            ⁡
            
              (
              n
              )
            
          
          ⁢
          
            
              
                q
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
            .
          
        
      
    
    
      
        (
        
          18
          ⁢
          e
        
        )
      
    
  




For Controller C33 





  
    
      
        
          
            
              p
              j
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                  33
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                  33
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              q
              j
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                ap
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
            +
            b
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              v
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                33
              
              ⁢
              
                
                  e
                  2
                
                ⁡
                
                  (
                  n
                  )
                
              
            
            +
            
              100
              ⁡
              
                [
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                            33
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                            33
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
                ]
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                w
                ij
                33
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            
              a
              2
            
            ⁢
            
              η
              33
            
            ⁢
            
              
                e
                2
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                E
                i
                33
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                h
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
    
      
        (
        
          19
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            
              h
              j
              33
            
            ⁡
            
              (
              n
              )
            
          
        
        =
        
          a
          ⁢
          
            
          
          ⁢
          
            η
            33
          
          ⁢
          
            
              e
              2
            
            ⁡
            
              (
              n
              )
            
          
          ⁢
          
            
              
                q
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
            .
          
        
      
    
    
      
        (
        
          19
          ⁢
          e
        
        )
      
    
  




In these equations, n denotes the nth iteration; a is an arbitrary constant and b=1⁄2; η11>0, η22>0, and η33>0 are the learning rate for C11, C22, and C33, respectively; pj11 , pj22 , pj33 , qj11 ,qj22, and qj33 are variables of the neural network for C11, C22 , and C33 , respectively; wij11, wij 22, wij33, hj11, hj22, and hj33 are weighting factors of the neural network for C11, C22, and C33, respectively; Kc11>0, Kc22>0, and Kc33>0 are the controller gain for C11, C22, and C33, respectively; and Ei11(n), Ei22(n), and Ei33(n) are the delayed signals of E111(n), E122(n), and E133(n), respectively. In addition, the a and b constants in Equations (17b), (17c), (18b), (18c), (19b), and (19c) are in the range specified in Equation (4b), and bounded by the limits specified in Equations (4a) and (4c).
The compensators C21, C31, C12, C32, C13, and C23 are designed to include a first-order dynamic block based on the following Laplace transfer functions:




  
    
      
        
          
            C
            21
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                21
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                11
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  21
                
                ⁢
                
                  K
                  c
                  21
                
              
              
                
                  
                    T
                    c
                    21
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            C
            31
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                31
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                11
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  31
                
                ⁢
                
                  K
                  c
                  31
                
              
              
                
                  
                    T
                    c
                    31
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            C
            12
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                12
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                22
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  12
                
                ⁢
                
                  K
                  c
                  12
                
              
              
                
                  
                    T
                    c
                    12
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            C
            32
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                32
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                22
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  32
                
                ⁢
                
                  K
                  c
                  32
                
              
              
                
                  
                    T
                    c
                    32
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            C
            13
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                13
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                33
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  13
                
                ⁢
                
                  K
                  c
                  13
                
              
              
                
                  
                    T
                    c
                    13
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          e
        
        )
      
    
  
  
    
      
        
          
            C
            23
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                23
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                33
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  23
                
                ⁢
                
                  K
                  c
                  23
                
              
              
                
                  
                    T
                    c
                    23
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          20
          ⁢
          f
        
        )
      
    
  




In these equations, S is the Laplace transform operator; V11(S), V22(S), V33(S), V21(S), V31(S), V12(S), V32(S), V13(S), and V23(S) are the Laplace transform of signals v11(t), v22(t), v33(t), v21(t), v31(t), v12(t), v32(t), v13(t), and v23(t), respectively; Kc21>0, Kc31>0, Kc12>0, Kc32>0, Kc13>0, and Kc23>0 are the compensator gain; and Tc21>0, Tc31>0, Tc12>0, Tc22>0, Tc13>0, and Tc23>0 are the compensator time constants. In the applications where only static compensation is considered, the corresponding time constant can be set to 0. If the sub-process G21=0, meaning that there is no interconnection from loop 1 to loop 2, the compensator C21 should be disabled by selecting Kc21=0. Similarly, if G12=0, one should select K12=0 to disable C12, and so on.

Kslm is the compensator sign factor, which is selected based on the acting types of the sub-processes as follows:

Kslm=1, if Gll and Glm have different acting types   (21a)

Kslm=−1, if Gll and Glm, have the same acting type   (21b)

where l=1, 2, 3; m=1, 2, 3; and l≠m. For instance, Ks12=1 if G11 and G12 have different acting types, and Ks12=−1 if G11 and G12 have the same acting types.

For control loop 1, the main controller C11 and its corresponding compensators C12 and C13 are in one group. The objective is to control the main process G11, and compensate for the interactions from sub-processes G12 and G13. In this case, the compensation signal V12 and V13 are used to deal with the interaction signals X12 and X13 from sub-processes G12 and G13, respectively. Notice that the interactions are from loop 2 and 3. For control loop 2, the main controller C22 and its corresponding compensators C21 and C23 are in one group. For control loop 3, the main controller C33 and its corresponding compensators C31 and C32 are in one group.
Our design philosophy for multivariable MFA control is not attempting a perfect cancellation for the loop interactions. Therefore, a simplified yet logical design for the estimated compensator time constant is as follows:




  
    
      
        
          
            C
            21
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                21
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                11
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  21
                
                ⁢
                
                  K
                  c
                  21
                
              
              
                
                  
                    T
                    c
                    22
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          22
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            C
            31
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                31
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                11
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  31
                
                ⁢
                
                  K
                  c
                  31
                
              
              
                
                  
                    T
                    c
                    33
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          22
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            C
            12
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                12
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                22
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  12
                
                ⁢
                
                  K
                  c
                  12
                
              
              
                
                  
                    T
                    c
                    11
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          22
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            C
            32
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                32
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                22
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  32
                
                ⁢
                
                  K
                  c
                  32
                
              
              
                
                  
                    T
                    c
                    33
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          22
          ⁢
          d
        
        )
      
    
  
  
    
      
        
          
            C
            13
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                13
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                33
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  13
                
                ⁢
                
                  K
                  c
                  13
                
              
              
                
                  
                    T
                    c
                    11
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          22
          ⁢
          e
        
        )
      
    
  
  
    
      
        
          
            C
            23
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                23
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                33
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  23
                
                ⁢
                
                  K
                  c
                  23
                
              
              
                
                  
                    T
                    c
                    22
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
    
      
        (
        
          22
          ⁢
          f
        
        )
      
    
  




In these equations, Tc11>0, Tc22>0, and Tc33>0 are the time constants for the main controllers C11, C22, and C33, respectively. They will be entered by the user when configuring the main controllers C11, C22, and C33, This design simplifies controller configuration with less parameter entry requirements of 6 compensator time constants. Since MFA controllers have inherent adaptive capability and a wider robust control range, the 3×3 MFA controller with simplified compensator design can still produce excellent control performance for seriously coupled 3×3 processes.

C. Conventional Boiler-Turbine-Generator (BTG) Units

A conventional coal, oil, or natural gas fired power plant consists of three key components: (1) Boiler, (2) Turbine, and (3) Generator. The Boiler produces superheated steam to turn the Turbine to allow the Generator to generate electricity. Operating as a set, the combined Boiler, Turbine, Generator, and all auxiliaries make up a BTG unit. FIG. 4 is a schematic representation of a Boiler-Turbine-Generator (BTG) unit of a power plant comprising a conventional boiler. Conventional here means that the boiler is sub-critical with a steam drum and also has a fixed bed in comparison to a fluidized-bed boiler.
To produce steam for power generation, water must be boiled. Feedwater first enters the Economizer where initial heating to almost boiling occurs. It then passes into a Steam Drum at the top of the boiler. From there the water recirculates through the Evaporator tubes lining inside the Furnace and back to the Steam Drum. About 20% of the water boils off and leaves the Drum through steam pipes in the top of the Drum. Superheater tubes mounted in the hottest part of the boiler heat the steam one final time.
The superheated steam from the Boiler is fed to the Turbine which is coupled with the Generator. Steam is exhausted from the Turbine at a low pressure, condensed, and then pumped back to the Boiler under pressure. The Turbine has three sections on a common shaft, the high-pressure cylinder (HP), the intermediate-pressure cylinder (IP), and one or more low-pressure cylinders (LP). Steam exhausted from the HP section is taken back to the Reheater to be reheated to its original temperature. In this way, more energy is put into the cycle without requiring further latent heat, thus improving the overall efficiency. From the Reheater, the steam goes to the IP section. Steam exhausted from the IP Section goes directly to the LP section and then passes to the Condenser which is operated at a near vacuum state. Forced draft (FD) fans are used to provide combustion air to the windbox for distribution to the burners. Induced draft (ID) fans are usually equipped to draw flue gases from the Furnace and drive them up the Stack.
The fuel such as coal is burned with a controlled amount of air to insure complete combustion and minimum pollution. Oil and natural gas are readily burned in the furnace; coal must be ground to a fine powder, known as pulverized fuel, before it can be readily burned. A large percent of boilers used for power and steam generation today are conventional subcritical pulverized coal (PC) boilers. Compared with more advanced boilers, they are lower in efficiency, higher in emissions, more fuel sensitive, and have less product versatility.
D. Advanced Power Boilers
Compared with sub-critical fixed bed conventional boilers, there are 3 types of advanced boilers: (1) Once-Through Supercritical (OTSC) Boilers, (2) Circulating Fluidized-Bed (CFB) Boilers, and (3) Once-Through Supercritical Circulating Fluidized-Bed (OTSC CFB) Boilers. Generally speaking, a power plant that is equipped with any number of advanced boilers can be called an advanced power plant.
Boilers used in energy plants are either “drum” or “once-through” types, depending on how the boiler water is circulated. Heat is transferred through the furnace tubes and into the water passing through the tubes to generate steam. In drum-type boilers, the steam-flow rate is typically controlled by the fuel-firing rate. In once-through boilers, the steam-flow rate is established by the boiler feedwater and the superheated steam temperature is controlled by the fuel-firing rate. A boiler is called supercritical when the master steam pressure is over 22.129 Mpa. In general, when water goes over the critical point (Pressure=22.129 Mpa, and Temperature=234 degree C.), it becomes steam. Therefore, a steam drum cannot be used and the Once-through design is the only choice for supercritical boilers. Once-through supercritical boilers run at higher steam temperature and pressure so that better energy efficiency is achieved. But they are difficult to control as summarized in Table 2.








TABLE 2





Challenges
Description and Comments







Severely
The relationship of throttle valve, fuel feed, and water


Nonlinear and
feed to power, steam pressure, and steam temperature


Multivariable
are nonlinear and interacting.


Serious
Because of the once-through design, there exists serious


Coupling
coupling between the boiler and turbine units.


Large
Since there is no steam drum, any changes in the throttle


Disturbances
valve position will cause a direct disturbance to the boiler



pressure and temperature.


Large load
Boiler needs to run in both subcritical and supercritical


and operating
modes causing large load and operating condition


condition
changes.


changes









Circulating fluidized-bed (CFB) boilers are becoming strategically important in power and energy generation. The unique design of CFB boilers allows fuel such as coal powders to be fluidized in the air so that they have better contact with the surrounding air for better combustion. CFB boilers can burn low-grade materials such as waste coal, wood, and refuse derived fuel. Most importantly, less emissions such as COx and NOx are produced compared to conventional boilers. The critical process variables and their control challenges for a CFB boiler are listed in Table 3.








TABLE 3





Process Variable
Control Challenges of CFB Boilers







Master
Nonlinear, tight specifications, large delay time, large


Steam
disturbance caused by load changes and poor feed


Pressure
actuation, etc.


Steam
Large time delay and time-varying.


Temperature


Bed
Multi-input-single-output process, multiple constraints,


Temperature
very critical since poor bed temp control results



in serious NOx emissions.


Excess
It is related to multiple emission constraints, varying


Oxygen
heating value of flexible fuel, and condition of the



oxygen sensors.


Furnace
Multiple fans and dampers to hold proper negative


Negative
pressure for the furnace.


Pressure


Coal or
Nonlinear, poor actuation, coal or fuel feed jams, etc.


Fuel Feed


Primary
Multiple fans and dampers to hold the proper


Air and
CFB circulating condition and fuel-air-ratio. Extremely


Secondary
sensitive to bed temperature.


Air









For a CFB boiler, the control challenges are mainly related to the combustion process of its furnace. The combustion process control of a CFB boiler is not in the scope of this patent.
E. Supercritical CFB Boilers and BTG Units
The Once-through Supercritical Circulating Fluidized-bed (OTSC CFB) boiler combines the merits of once-through supercritical and circulating fluidized-bed technologies. As a strategically important clean coal technology, OTSC CFB boilers can significantly improve combustion and energy efficiency, reduce emissions, and have fuel flexibility. It is the most promising boiler for future energy plants because of all its outstanding advantages.
A Supercritical CFB boiler based electric power plant also consists of three key components: (1) Boiler, (2) Turbine, and (3) Generator. Similar to the conventional boiler, the OTSC CFB boiler produces superheated steam to turn the turbine to allow the generator to generate electricity. Operating as a set, the combined Boiler, Turbine, Generator, and all auxiliaries make up a BTG unit. FIG. 5 is a schematic representation of a Boiler-Turbine-Generator (BTG) unit of a power plant comprising a Supercritical Circulating Fluidized-bed boiler. Feedwater first enters the Economizer where initial heating to almost boiling occurs. It then passes into the Cyclone Separator at the top of the Boiler. From there the water recirculates through the Superheaters. The superheated steam is fed directly to the Turbine which is coupled with the Generator. Steam is exhausted from the Turbine at a low pressure, condensed, and then pumped back to the boiler under pressure.
For an OTSC CFB boiler, most of the control challenges in Once-Through Supercritical boilers and in CFB boilers still exist. Since the Supercritical CFB boiler combines the chaotic operating conditions of a CFB boiler and the once-through nature of a supercritical boiler, the control challenges could double. For such a boiler, maintaining a dynamic material and energy balance becomes a big challenge. In general, for an OTSC CFB boiler, its BTG process and its CFB combustion process are much more dependent on a good automatic control system in order to keep the energy and material balance. If not careful, the entire system can get into vicious cycles causing serious consequences. For instance, when a steam demand increases, it will cause the steam pressure to go down, which will quickly affect the boiler firing and then the fluidized-bed conditions. The changed combustion condition will result in more changes in steam temperature and pressure and therefore a vicious cycle will build up causing major operation and safety problems. Conventional control methods including coordinated control of steam turbine and boiler control will have major difficulties in controlling OTSC CFB boilers.
In a power generation network, a BTG unit may be base-loaded to generate at a constant rate, or may cycle up and down as required by an automatic dispatch system. In either case, the boiler control system manipulates the firing rate of the furnace to generate the steam required to satisfy the demand for power. It is also necessary to maintain an adequate supply of feedwater and the correct mixture of fuel and air for safe and economic combustion. These requirements are actually the same for a conventional BTG unit or a BTG unit that employs an advanced power boiler such as a Supercritical boiler, a CFB boiler, or a Supercritical CFB boiler.
FIG. 6 is a diagram illustrating a simplified schematic of the BTG unit along with the key process variables, which are described in Table 4.









TABLE 4





Variable
Symbol
Description







Throttle
VT
The valve used for the Turbine governor control.


Valve


Position


Firing
RF
The firing rate of the boiler is changed by


Rate

manipulating the amounts of air and fuel to the




burners. Increasing the firing rate generates




more steam.


Water
FW
The feed water flow to the boiler.


Feed


Power
JT
The power measurement is used to indicate and


Output

control the power generation of the BTG unit.


Throttle
PT
The throttle pressure is actually the steam supply


Pressure

pressure to the turbine. It indicates the state of




balance between the supply and demand for steam.




Rising throttle pressure indicates that the steam




supply exceeds demand and falling throttle pressure




indicates that the steam demand exceeds supply.




The automatic controller for this purpose is the




Turbine Governor.


Steam
Fs
The steam flow.


Flow


Steam
T1
Temperature of superheated steam in position 1.


Temp 1









F. MFA Control System Design Method and MFA Control of BTG Units

In this section, we introduce a novel multivariable control system design method which includes the following key ideas:

    1. The control system design is based on qualitative analysis of the process input and output variables. No detailed quantitative analysis or process models are required.
    2. For a multivariable process, use S (Strong), M (Medium), and W (Weak) to represent the degree of connections between the input and output of each sub-process. Use the plus or minus sign to represent whether the process is direct or reverse acting.
    3. Properly pair the process input and output variables so that the main processes are open-loop stable and have a strong direct or reverse acting relationship to assure good controllability.
    4. The remaining sub-processes should have medium, weak, or even no connections between their input and output variables. Their acting types do not matter.
    5. If a sub-process has a strong relationship between its input and output, either improve the process or carefully launch the control system.


Based on this control system design method, an MFA control system is designed to control the critical process variables of the BTG unit including Power (JT), Steam Throttle Pressure (PT), and Steam Temperature T1. The process has 3 inputs and 3 outputs and is called a Power-Pressure-Temperature (PPT) process. We pair the variables to form 9 sub-processes G11, G21, . . . , G33 as listed in Table 5.









TABLE 5









Process Outputs -



Process Variables to be Controlled










Process Inputs -

Throttle Pressure



Manipulated Variables
Power (JT)
(PT)
Steam Temp (T1)





Firing Rate (RF)
G11
G21
G31


Throttle Valve (VT)
G12
G22
G32


Water Feed (FW)
G13
G23
G33









The importance of the variable pairing is that we want to make sure the 3 main processes G11, G22, and G33 have a strong direct or reverse acting relationship so that they have good controllability. The detailed qualitative input and output relationship among all 9 sub-processes is analyzed and presented in Table 6.










TABLE 6






Input-
Acting



Process
Output
Type
Qualitative Input and Output Relationship







G11
RF-JT
S
Strong direct acting. When the firing rate





increases, the power will increase but the





relationship is quite nonlinear.


G21
RF-PT
M to S
When the firing rate increases, the throttle





pressure will increase gradually.


G31
RF-T1
M to S
When the firing rate increases, the steam





temperature will increase gradually.


G12
VT-JT
+/− S
When the throttle valve increases, the power





will increase immediately for a short period





of time. If the firing rate does not change,





power will go down and back to its original





position.


G22
VT-PT
−S
Strong reverse acting. When the throttle





valve position increases, the throttle pressure





decreases.


G32
VT-T1
W
When the throttle valve increases, it has little





effect on the temperature.


G13
FW-JT
+/−M
When the water feed increases, the power





will increase for a short period of time





and then decrease quite significantly.


G23
FW-PT
+/−M
When the water feed increases, the throttle





pressure will increase immediately and then





decrease gradually.


G33
FW-T1
−S
Strong reverse acting. When the water feed





increases, it will cause the steam temperature





to decrease.









Based on the variable pairing and control system design, the 3×3 MFA controller illustrated in FIG. 3 can be readily applied to control the PPT process of a BTG unit. The BTG unit may include a conventional power boiler or an advanced power boiler such as a Supercritical boiler, a CFB boiler, or a Supercritical CFB boiler.
FIG. 7 illustrates the block diagram of a 3×3 MFA control system for controlling the 3×3 Power-Pressure-Temperature (PPT) process of a Boiler-Turbine-Generator (BTG) unit according to this invention. The MFA control system comprises a 3×3 MFA controller 62, a 3×3 PPT process of a BTG unit 64, a Firing Rate and Combustion Sub-system 66, a Throttle Valve and Steam Flow Sub-system 68, and a Water Flow Sub-system 70. Each of the sub-systems may include various control loops. For instance, the Water Flow Sub-system typically includes a water flow control loop. In this case, the control signal u3 will be used as the setpoint for the water flow control loop, which is the inner loop of the cascade control system. Conventional control methods could be used to control these sub-systems. Therefore, they are not in the scope of this patent.
The 3×3 PPT process has nine sub-processes G11 through G33 as listed in Table 6. The process variables y1, y2, and y3 are Power (JT), Steam Throttle Pressure (PT), and Steam Temperature T1, respectively. They are the feedback signals for each of the main control loops and compared with the setpoints r1, r2 and r3 at adders 72 to produce error signals e1, e2, and e3. The outputs of the 3×3 MFA controller u1, u2, and u3 manipulate the manipulated variables Firing Rate (RF), Throttle Valve (VT), and Water Feed (FW) to control the Power (JT), Steam Throttle Pressure (PT), and Steam Temperature T1, respectively.
FIG. 8 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoint for Power is stepped down. In FIG. 8, curves 74, 75, 76 are SP, PV, OP of the MFA Power loop, and curves 77, 78, 79 are SP, PV, OP of the PID Power loop, respectively. Curves 80, 81, 82 are SP, PV, OP of the MFA Pressure loop, and curves 83, 84, 85 are SP, PV, OP of the PID Pressure loop, respectively. Curves 86, 87, 88 are SP, PV, OP of the MFA Temperature loop, and curves 89, 90, 91 are SP, PV, OP of the PID Temperature loop, respectively. The loop interactions can be clearly seen. When the Power SP (Signals 74 and 77) is changed from 60 to 35, the controller OP (Signals 76 and 79) produces the control actions trying to force the Power PV (Signals 75 and 78) to track its setpoint. Since it is a 3×3 process, this action inevitably causes the Pressure PV and Temperature PV to change as well.
FIG. 9 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoint for Steam Throttle Pressure is stepped down. In FIG. 9, curves 92, 93, 94 are SP, PV, OP of the MFA Power loop, and curves 95, 96, 97 are SP, PV, OP of the PID Power loop, respectively. Curves 98, 99, 100 are SP, PV, OP of the MFA Pressure loop, and curves 101, 102, 103 are SP, PV, OP of the PID Pressure loop, respectively. Curves 104, 105, 106 are SP, PV, OP of the MFA Temperature loop, and curves 107, 108, 109 are SP, PV, OP of the PID Temperature loop, respectively. Please notice that the throttle valve change has little effect on the temperature loop as described for sub-process G32 in Table 6.
FIG. 10 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoint for Steam Temperature is stepped down. In FIG. 10, curves 110, 111, 112 are SP, PV, OP of the MFA Power loop, and curves 113, 114, 115 are SP, PV, OP of the PID Power loop, respectively. Curves 116, 117, 118 are SP, PV, OP of the MFA Pressure loop, and curves 119, 120, 121 are SP, PV, OP of the PID Pressure loop, respectively. Curves 122, 123, 124 are SP, PV, OP of the MFA Temperature loop, and curves 125, 126, 127 are SP, PV, OP of the PID Temperature loop, respectively. As illustrated, the MFA control system can suppress the disturbances in the Power and Pressure loops caused by the big change in the water flow. In contrast, the same disturbance caused the PID loops to swing.
FIG. 11 is a time-amplitude diagram comparing the control performance of a 3×3 MFA control system versus a PID control system for controlling two identical PPT processes of a BTG unit, where the setpoints for Power, Throttle Pressure, and Steam Temperature all have step changes. In FIG. 11, curves 128, 129, 130 are SP, PV, OP of the MFA Power loop, and curves 131, 132, 133 are SP, PV, OP of the PID Power loop, respectively. Curves 134, 135, 136 are SP, PV, OP of the MFA Pressure loop, and curves 137, 138, 139 are SP, PV, OP of the PID Pressure loop, respectively. Curves 140, 141, 142 are SP, PV, OP of the MFA Temperature loop, and curves 143, 144, 145 are SP, PV, OP of the PID Temperature loop, respectively. In this case, the Power SP (Signals 128 and 131) is firstly stepped down from 58 to 34, the Pressure SP (Signals 134 and 137) is then stepped up from 40 to 62, and the Temperature SP (Signals 140 and 143) is lastly stepped up from 24 to 42. It can be seen that each setpoint change causes disturbances to all control loops. The 3×3 MFA control system is able to deal with the disturbances and keeps the Power, Pressure, and Temperature under control. In contrast, the PID control system cannot effectively control the 3×3 process resulting in oscillations in all 3 loops.
To summarize, the control trends demonstrate outstanding control performance of the 3×3 MFA control system for both tracking and regulating capabilities. The compensators inside the 3×3 MFA controller can effectively decouple and reduce the interactions from the other loops of the 3×3 PPT process. The control trends also demonstrate un-satisfactory control performance of the PID control system. Since PID controllers are single-loop controllers and can only treat the 3-Input-3-Output (3×3) multivariable process as three single-input-single-output (SISO) processes, it is very difficult for the PID control system to be effective and achieve good control performance. When there is a setpoint change or disturbance in the process, it will take a long time for the loops to settle down due to interactions among the loops. For instance, when the setpoint of Loop 1 is changed, the PID control action in Loop 1 will disturb Loop 2 and 3 causing their PID controllers to produce control actions, which will come back to disturb Loop 1. The multiple and bi-directional interactions can cause conflicting control actions and trigger a vicious cycle resulting in loop oscillations. Therefore, when applying PID for multivariable control, most PID controllers are significantly de-tuned to avoid potential oscillations or even unstable control. In the real world, a large percentage of multi-input-multi-output (MIMO) processes are treated as single-input-single-output (SISO) processes resulting in poor control performance, inconsistent quality, wasted materials and energy, and plant safety problems. The 3×3 MFA control system described in this patent can be very useful for controlling key process variables of boiler-turbine-generator (BTG) units in conventional and advanced power plants as well as other multivariable processes.

## Claims

1. A method of controlling a 3-Input-3-Output (3×3) process having three main-processes and six sub-processes, each of said main-processes having an input and an output responsive to the said input, each of said sub-processes having an input and an output, each of said main-processes having a measured process variable, which is the summation of the output of the main-process and outputs from two corresponding sub-processes; the 3-Input-3-Output (3×3) process having unknown relationships between their inputs and outputs, without approximating or modeling said relationship, comprising:
a) selecting a setpoint representing a desired value for the measured process variable for each of the main-processes;
b) obtaining an error value which is a function of the difference between said setpoint and said measured process variable for each of the main-processes;
c) providing a hardware 3-Input-3-Output (3×3) controller that has three main controllers, each corresponding to one of the main-processes, and six compensators, each corresponding to one of the sub-processes;
d) for each of the main processes, applying said error value as the sole input to its corresponding main controller whose inputs are time-delayed functions of said error value, and whose output is a control value combined with the outputs from two adjacent compensators to produce the total control output of the main controller;
e) configuring each said main controller by entering a time constant based on the time constant of the corresponding main-process being controlled;
f) for each of the sub-processes, providing a compensator which uses the control value from the corresponding main controller as its input, and uses the time constant of the main controller as its time constant;
g) applying the total control output of each main controller to vary the input of its corresponding main-process and two sub-processes; and
h) continuously iteratively varying the parameters of each of the main controllers to help minimize said error value for each of the main-processes.

2. The method of claim 1, further comprising adding said error value to the control value for each of the main-processes.

3. The method of claim 1, for each of the main-processes, further comprising adding the error value to said control value, in which the error value can be increased or decreased by a manually adjustable gain parameter.

4. The method of claim 1, wherein the 3-Input-3-Output (3×3) process is a Power-Pressure-Temperature (PPT) process of a Boiler-Turbine-Generator (BTG) unit of a conventional coal-fired power boiler, or a Once-Through Supercritical (OTSC) Boiler, or a Circulating Fluidized-Bed (CFB) Boiler, or a Once-Through Supercritical Circulating Fluidized-Bed (OTSC CFB) Boiler.

5. A 3-input-3-output (3×3) model-free adaptive (MFA) control system comprising:
a) a 3-Input-3-Output (3×3) process to be controlled, further comprising:
i) three main-processes G11, G22, G33 having their corresponding process outputs x11(t), x22(t), x33(t), responsive to their corresponding control signals u1(t), u2(t), u3(t) respectively, and
ii) six sub-processes G21, G31, G12, G32, G13, G23 having their corresponding process outputs x21(t), x31(t), x12(t), x32(t), x13(t), x23(t), wherein the process produces three measured process variables y1(t), y2(t), y3(t) for the corresponding main-processes, each of the process variables being affected by the output of its corresponding main-process, the outputs of two corresponding sub-processes signals as follows:
y1(t)=x11(t)+x12(t)+x13(t)+d1(t),
y2(t)=x21(t)+x22(t)+x23(t)+d2(t),
y3(t)=x31(t)+x32(t)+x33(t)+d3(t);


wherein d1(t), d2(t) and d3(t) are disturbance signals; and
b) a hardware 3-Input-3-Output (3×3) model-free adaptive (MFA) controller further comprising:
i) three main-controllers C11, C22, C33 responsive to error signals e1(t), e2(t) and e3(t) respectively and having their corresponding control outputs y11(t), v22(t) v33(t); and
ii) six compensators C21, C31, C12, C32, C13, C23 having their corresponding outputs v21(t), v31(t), v12(t), v32(t), v13(t), v23(t) wherein:
iii) the error signals e1(t), e2(t), e3(t) are iteratively calculated as follows:
e1(t)=r1(t)−y1(t), if G11 is direct acting
e1(t)=[r1(t)−y1(t)], if G11 is reverse acting
e2(t)=r2(t)−y2(t), if G22 is direct acting
e2(t)=−[r2(t)]−y2(t)], if G22 is reverse acting
e3(t)=r3(t)−y3(t), if G33 is direct acting
e3(t)=−[r3(t)]−y3(t)], if G33 is reverse acting

where a process is direct acting if the process output increases as the process input increases, and a process is reverse acting, where a process is direct acting if the process output increases as the process input increases, and a process is reverse acting if the process output decreases as the process input increases, and r1(t), r2(t) and r3(t) are set points;
iv) each of the compensators uses a control output from a corresponding main-controller as its input, and uses the time constant of the corresponding main controller as its time constant; and
v) each of the main-controllers has a combined control output being iteratively calculated based on its corresponding control output and the outputs from two corresponding compensators as follows:
u1(t)=v11(t)+v12(t)+v13(t),
u2(t)=v21(t)+v22(t)+v23(t),
u3(t)=v31(t)+v32(t)+v33(t),

wherein, the combined control output of each main-controller is the input to its corresponding main-process and two sub-processes in order to control the measured process variable of its corresponding main-process.

6. w) The control system of claim 5, wherein:
a) the three main-controllers C11, C22, and C33, each have a corresponding normalized error value E111(n), E122(n), and E133(n) being computed with a normalization function N(.) as follows:





  
    
      
        E
        1
        11
      
      ⁡
      
        (
        n
        )
      
    
    =
    
      
        
          K
          c
          11
        
        
          T
          c
          11
        
      
      ⁢
      
        N
        ⁡
        
          (
          
            
              e
              1
            
            ⁡
            
              (
              t
              )
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        E
        1
        22
      
      ⁡
      
        (
        n
        )
      
    
    =
    
      
        
          K
          c
          22
        
        
          T
          c
          22
        
      
      ⁢
      
        N
        ⁡
        
          (
          
            
              e
              2
            
            ⁡
            
              (
              t
              )
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        E
        1
        33
      
      ⁡
      
        (
        n
        )
      
    
    =
    
      
        
          K
          c
          33
        
        
          T
          c
          33
        
      
      ⁢
      
        N
        ⁡
        
          (
          
            
              e
              3
            
            ⁡
            
              (
              t
              )
            
          
          )
        
      
    
  
  ,




in which, n denotes the nth iteration; e1(t), e2(t), and e3(t) are the error signals as the controller inputs, Kc1l>0, Kc22>0, and Kc33>0 are the controller gains, and Tc11>0, Tc22>0, and Tc33>0 are the time constants for C11, C22, and C33, respectively;
b) the main-controller C11, C22, and C33 each includes a neural network with an input layer including a plurality of input neurons arranged to receive normalized, scaled and delayed forms of said error signal, a hidden layer including a plurality of hidden neurons each arranged to sum the signals received by each of said input neurons weighted by an individual first weighting factor, and an output neuron which sums the individually weighted outputs of hidden neurons, the weighting factors for said hidden neuron outputs being iteratively adjusted, and an activation function f(x) substantially of the form





  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          0
        
        ,
        
          
            if
            ⁢
            
              
            
            ⁢
            x
          
          <
          
            -
            
              b
              a
            
          
        
      
    
  
  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          
            ax
            +
            b
          
        
        ,
        
          
            if
            ⁢
            
              
            
            -
            
              b
              a
            
          
          ≤
          x
          ≤
          
            b
            a
          
        
      
    
  
  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          1
        
        ,
        
          
            if
            ⁢
            
              
            
            ⁢
            x
          
          >
          
            b
            a
          
        
      
    
  




where a is an arbitrary constant and b=1⁄2;
c) the main-controllers are arranged to produce their control outputs v11(t), v22(t), v33(t) substantially based on the following difference equations and online learning algorithms of the neural network:

For Controller C11:




  
    
      
        
          
            
              p
              j
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                  11
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                  11
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            
              q
              j
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                ap
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
            +
            b
          
        
        ,
      
    
  
  
    
      
        
          
            
              v
              11
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                11
              
              ⁢
              
                
                  e
                  1
                
                ⁡
                
                  (
                  n
                  )
                
              
            
            +
            
              100
              ⁡
              
                [
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                            11
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                            11
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
                ]
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                w
                ij
                11
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            
              a
              2
            
            ⁢
            
              η
              11
            
            ⁢
            
              
                e
                1
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                E
                i
                11
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                h
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                h
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              η
              11
            
            ⁢
            
              
                e
                1
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                q
                j
                11
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ;
      
    
  




For Controller C22




  
    
      
        
          
            
              p
              j
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                  22
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                  22
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            
              q
              j
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                ap
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
            +
            b
          
        
        ,
      
    
  
  
    
      
        
          
            
              v
              22
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                22
              
              ⁢
              
                
                  e
                  2
                
                ⁡
                
                  (
                  n
                  )
                
              
            
            +
            
              100
              ⁡
              
                [
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                            22
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                            22
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
                ]
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                w
                ij
                22
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            
              a
              2
            
            ⁢
            
              η
              22
            
            ⁢
            
              
                e
                2
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                E
                i
                22
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                h
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                h
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              η
              22
            
            ⁢
            
              
                e
                2
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                q
                j
                22
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ;
      
    
  




For Controller C33




  
    
      
        
          
            
              p
              j
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              N
            
            ⁢
            
              
                
                  w
                  ij
                  33
                
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  E
                  i
                  33
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            
              q
              j
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                ap
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
            +
            b
          
        
        ,
      
    
  
  
    
      
        
          
            
              v
              33
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              
                K
                c
                33
              
              ⁢
              
                
                  e
                  2
                
                ⁡
                
                  (
                  n
                  )
                
              
            
            +
            
              100
              ⁡
              
                [
                
                  
                    a
                    ⁢
                    
                      
                        ∑
                        
                          j
                          =
                          1
                        
                        N
                      
                      ⁢
                      
                        
                          
                            h
                            j
                            33
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                        ⁢
                        
                          
                            q
                            j
                            33
                          
                          ⁡
                          
                            (
                            n
                            )
                          
                        
                      
                    
                  
                  +
                  b
                
                ]
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                w
                ij
                33
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            
              a
              2
            
            ⁢
            
              η
              33
            
            ⁢
            
              
                e
                2
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                E
                i
                33
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                h
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ,
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                h
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
          
          =
          
            a
            ⁢
            
              
            
            ⁢
            
              η
              33
            
            ⁢
            
              
                e
                2
              
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                q
                j
                33
              
              ⁡
              
                (
                n
                )
              
            
          
        
        ;
      
    
  




in which, n denotes the nth iteration; a is an arbitrary constant and b=1⁄2 which are bounded by the activation function f(x); η11>0, η22>0, and η33>0 are the learning rates for C11, C22, and C33, respectively; pj11, pj22, pj33, qj11, qj22, and qj33 are variables of the neural network for C11, C22, and C33, respectively; wij11, wij22, wij33, hj11, hj22, and hj33 are weighting factors of the neural network for C11, C22, and C33, respectively; Kc11>0, Kc22>0, and Kc33>0 are the controller gains for C11, C22, and C33, respectively; and Ei11(n), Ei22(n), Ei33(n) are the delayed signals of E111(n), E122(n), and E133(n), respectively; and
d) the compensators have their corresponding outputs C21(S), C31(S), C12(S), C32(S), C13(S), C23(S) substantially in the form of the following Laplace transfer functions:





  
    
      
        
          
            C
            21
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                21
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                11
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  21
                
                ⁢
                
                  K
                  c
                  21
                
              
              
                
                  
                    T
                    c
                    21
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
  
  
    
      
        
          
            C
            31
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                31
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                11
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  31
                
                ⁢
                
                  K
                  c
                  31
                
              
              
                
                  
                    T
                    c
                    31
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
  
  
    
      
        
          
            C
            12
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                12
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                22
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  12
                
                ⁢
                
                  K
                  c
                  12
                
              
              
                
                  
                    T
                    c
                    12
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
  
  
    
      
        
          
            C
            32
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                32
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                22
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  32
                
                ⁢
                
                  K
                  c
                  32
                
              
              
                
                  
                    T
                    c
                    32
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
  
  
    
      
        
          
            C
            13
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                13
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                33
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  13
                
                ⁢
                
                  K
                  c
                  13
                
              
              
                
                  
                    T
                    c
                    13
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
  
  
    
      
        
          
            C
            23
          
          ⁡
          
            (
            S
            )
          
        
        =
        
          
            
              
                V
                23
              
              ⁡
              
                (
                S
                )
              
            
            
              
                V
                33
              
              ⁡
              
                (
                S
                )
              
            
          
          =
          
            
              
                
                  K
                  s
                  23
                
                ⁢
                
                  K
                  c
                  23
                
              
              
                
                  
                    T
                    c
                    23
                  
                  ⁢
                  S
                
                +
                1
              
            
            .
          
        
      
    
  




in which, S is the Laplace transform operator; V11(S), V22(S), V33(S), V21(S), V31(S), V12(S), V32(S), V13(S), and V23(S) are the Laplace transform of signals v11(t), v22(t), v33(t), v21(t), v31(t), v12(t), v32(t), v13(t), and v23(t), respectively; Kc21>0, Kc31>0, Kc12>0, Kc32>0, Kc13>0, and Kc23>0 are the compensator gains; Ks21, Ks31, Ks12, Ks32, Ks13, Ks23 are the compensator sign factors, and Tc11>0, Tc22>0, and Tc33>0 are the time constants for the main controllers C11, C22, and C33, respectively.

7. The control system of claim 6, in which the compensators have said compensator sign factors being selected based on the acting types of the main-processes Gll and sub-processes Glm as follows:
Kslm=1, if Gll and Glm have different acting types
Kslm=−1,if Gll and Glm, have the same acting type

where l=1, 2, 3; m=1, 2, 3; and l≠m.

8. The control system of claim 5, in which the 3-Input-3-Output (3×3) process is a Power-Pressure-Temperature (PPT) process of a Boiler-Turbine-Generator (BTG) unit of an electrical power or energy generation plant, in which the 3×3 process has the following main-processes and sub-processes:
for G11, the input is Firing Rate and output is Power,
for G22, the input is Throttle Valve Position and output is Steam Throttle Pressure,
for G33, the input is Water Feed and output is Master Steam Temperature,
for G21, the input is Firing Rate and output is Steam Throttle Pressure,
for G31, the input is Firing Rate and output is Master Steam Temperature,
for G12, the input is Throttle Valve Position and output is Power,
for G32, the input is Throttle Valve Position and output is Master Steam Temperature,
for G13, the input is Water Feed and output is Power,
for G23, the input is Water Feed and output is Steam Throttle Pressure.

9. The control system of claim 8, in which the 3-Input-3-Output (3×3) process is a Power-Pressure-Temperature (PPT) process of a Boiler-Turbine-Generator (BTG) unit of a conventional coal-fired power boiler, or a Once-Through Supercritical (OTSC) Boiler, or a Circulating Fluidized-Bed (CFB) Boiler, or a Once-Through Supercritical Circulating Fluidized-Bed (OTSC CFB) Boiler.