Background
[0001] The present disclosure pertains to industrial systems and particularly to parallel
working units in such systems. More particularly, the disclosure pertains to obtaining
models of the units for advanced process control.
Summary
[0002] The disclosure reveals an approach for modeling parallel working units of a system
for advanced process control. The approach may be a systematic solution based on structured
model order (model complexity) reduction. Two phases of it may incorporate model identification
and model combination. The first phase is where a model of each parallel unit and
a model of the remaining system without any unit may be obtained. The second phase
is where the models may be combined to obtain a model of the whole system for any
configuration needed by the advanced process control. The model of the whole system
may be subjected to a structured model reduction to obtain a reduced order model for
the advanced process control.
Brief Description of the Drawing
[0003]
Figure 1 is a diagram of parallel working units as a part of a larger system;
Figure 2 is a diagram of a workflow with two phases;
Figure 2a is a diagram of a controller for processing the workflow of the items in
the diagram of Figure 2 for an example system of Figure 1;
Figure 3 is a diagram of a boiler block scheme;
Figure 4 is a diagram of an N number of boilers feeding a single header;
Figures 5a, 5b, 5c and 5d are data graphs of header pressure versus time for fuel
flow, of header pressure versus time for steam demand, of stream flow to header versus
time for fuel flow, and of steam flow to header versus time for steam demand, respectively.
Figures 6a, 6b, 6c and 6d are data graphs of header pressure versus time for fuel
flow, of header pressure versus time for steam demand, of flow to header versus time
for fuel flow, and of flow to header versus time for steam demand, respectively;
Figure 6e is a graph of scale for steam demand in full model, structured and balanced
versions;
Figures 7a, 7b, 7c and 7d are data graphs of header pressure error versus time for
fuel flow, of header pressure versus time for steam demand, of flow to header versus
time for fuel flow, and of flow to header versus time for steam demand, respectively;
Figures 8a, 8b, 8c and 8d are data graphs of header pressure versus time for fuel
flow, of header pressure versus time for steam demand, of flow to header versus time
for fuel flow, and of flow to header versus time for steam demand, respectively;
Figures 9a, 9b, 9c and 9d are data graphs of header pressure error versus time for
fuel flow, of header pressure error versus time for steam demand, of flow to header
error versus time for fuel flow, and of flow to header error versus time for steam
demand;
Figures 10a-10d and 11a-11d are data graphs of step responses and frequency responses,
respectively, of structure vs. unstructured reduction, for the three boiler and single
header system; and
Figures 12a-12d and 13a-13d are data graphs of step responses and frequency responses,
respectively, for an example system having five boilers.
Description
[0004] Control issues in industrial settings may often involve parallel working units, such
as: 1) parallel connection of multiple boilers feeding a single header; and 2) parallel
working pumps / turbines / chemical reactors / and so forth.
[0005] These parallel units may usually be operated in multiple different on/off configurations,
where individual units are turned on/off according to process needs and optimal allocation
schemes. A design of advanced process control (APC) may require a low order model
of the full plant for each on/off configuration. This means that system identification
needs to be performed for virtually all of these configurations, which may be rather
expensive and time consuming, or even somewhat impossible, since the number of configurations
might be in the order of hundreds.
[0006] Solutions used in the related art may be either do virtually all identification experiments
or identify models of individual parallel units and remaining technology, and combine
them to get a low order plant model by using certain heuristics. The last approach
may give fairly good results, but often with inconsistent quality and no guarantees.
[0007] The present approach may avoid doing identification experiments for virtually all
configurations. This approach may replace heuristics and deliver improved models with
consistent quality. Models of individual parallel units may be identified separately
and then be combined into a global model with a reduced order and an arbitrary on/off
configuration by using an algorithm of structure preserving model order reduction.
[0008] The present approach may involve an application of a structured order reduction to
a related art issue, thus leading to models with high and consistent quality while
keeping required experimental time and resources at a minimum. The approach concerns
modeling of parallel working units for APC. An algorithm may be noted herein.
[0009] Advanced process control (APC) of systems including parallel working units (boilers,
turbines, pumps, and so forth) may require a model of system dynamics for each on/off
configuration of the parallel units, which can be used during system operation. These
models may be obtained by step testing the whole system for every on/off configuration.
However, this cannot practically be used for systems having a large number of parallel
unit configurations (tens and more).
[0010] An issue may be that (due to multiple reasons) the system models for different configurations
cannot be obtained by a simple combination of individual unit models M
i and a model of the remaining system P (Figure 1). The main reason may be that parallel
units have to be modeled with a reduced order (controllability issue) under closed
loop operation.
[0011] Certain heuristics may be used with quite good results but with no guaranties. Also,
they need to be tailored for specific systems.
[0012] The present approach may be a systematic solution based on structured model order
reduction that gives very good and consistent results. Structured model order reduction
may be known. Figure 1 is a diagram of parallel working units as a part of a larger
system. Overall a modeled system 21 is shown. Major components may incorporate parallel
units 22 and remaining system 23 (P). Parallel units may incorporate a total of an
N number of units. Unit 24 (M
1) may be the first unit and unit 25 (M
N) may be the last unit of an N number of units. There may be N-2 units (M
i) (not shown) between the first and last units. Each unit may be connected in series
with switch (S). A switch may have symbolical meaning in the Figures. The switch may
represent that an appropriate unit can be switched on/off. Physically, the switch
may be a valve, boiler coal feeder, and so forth, or just the possibility to shut
down or turn on a unit.
[0013] A switch 26 (S
1) may connected in series with unit 24 and switch 27 (S
N) may be connected in series with unit 25. The may be a switch (S
i) connected in series with each of the units (M
i) situated between the first and last units. The switches may be situated at the inputs
of the units.
[0014] The ends not connected to the inputs of the units may be connected together on line
28 (u) which is connected to an output of remaining system 23. The outputs of the
units may added at a summer 29. An output line 31 (y) of the summed inputs may go
to an input of remaining system 23. System 21 inputs (w) for model system 21 may be
of an input on a line 32 to remaining system 23. System 21 outputs (z) may be of an
output on a line 33 from remaining system 23.
[0015] Figure 2 is a diagram of a workflow with two phases. The first phase may be a model
identification phase 41, where models of each parallel unit and model of the remaining
system are obtained. These results may be used in the second phase 42 (i.e., a model
combination phase), where they may be used to get a model of the whole system for
any configuration, which is needed by an advanced process control.
[0016] A model identification phase may have the following items. A first item 43 may be
doing an identification experiment (step testing) with unit 1 enabled (switch S
1 on) and other units disabled (switches S
2,...,S
N off) . At a second item 44, model parameters may be computed from experimental data
and then be used to extract a model of unit M
1. Third and fourth items 45 and 46 may repeat the identification experiment for all
of the remaining N units (generally always with only one unit enabled). In a fifth
item 47, arbitrary previous experimental data may be used to compute parameters of
a P model revealing the dynamics of remaining system P. There may be two sources of
data for P model parameters. Either it may be possible to identify P model parameters
from experiment data for one of individual units (in many cases), or it may be necessary
to perform an individual experiment on the "remaining system".
[0017] The model combination phase 42 may incorporate an input of an arbitrary configuration
of switches S
1, ...,S
N (as required by an advanced process control controller) at a sixth item 48 where
a parallel combination of enabled models may be computed (according to a configuration
of the switches). At a seventh item 49, a structured model order reduction (selected
Mi models in a closed loop with a P model) may be applied to the combination of enabled
parallel models and a model of the remaining system. The order of parallel models
only may be reduced. At an eighth item 51, a result may be a reduced order model for
use in an advanced process control controller.
[0018] An algorithm may be used in the model combination phase. Assumptions may incorporate
known models of individual parallel units M
1, ..., M
N, and a known model of the remaining system (e.g., the remaining technology) P. Inputs
may incorporate a configuration of switches S
1, ... ,S
N.
[0020] The following items may involve a structured model reduction.
[0021] Third, controllability Gramian P and observability Gramian Q may be computed for
the whole model and separate them according to
xP and
xM dimensions,

[0022] Fourth, the Cholesky factor
R of controlability Gramian
PP =
RRT may be computed.
[0023] Fifth, the singular value decomposition (SVD) of
RTQPR may be computed as
RTQPR=U∑
2UT.
[0024] Sixth, the transformation matrix
T may be computed as
T =
RÛ∑̂
-1/2, where are parts of
U,∑ according to a target model order.
[0025] Seventh, one may apply the transformation

[0026] From the applying the transformation, one may get a reduced state space model
AP, BP, CP, DP of parallel units as (a truncation)

or as (a singular perturbation)

[0028] An algorithm for structured model reduction may be similar to the Sandberg algorithm
of a structured balanced reduction as indicated in
Henrik Sandberg, Richard M. Murray, "Model reduction of interconnected linear systems",
Optimal Control, Applications and Methods, Special Issue on Directions, Applications,
and Methods in Robust Control, 30:3, pp. 225-245, May/June 2009 (Sandberg algorithm). The present approach may be essentially an application of structured
model reduction algorithm to an issue of modeling parallel working units under closed-loop
for advanced process control (APC). The present approach solution may eliminate the
need to do an identification experiment for all configurations of parallel units and
require a need to do the same number of identification experiments as the number of
parallel units.
[0029] Figure 2a is a diagram of a controller 72 for processing the workflow of the items
in the diagram of Figure 2 for an example system shown in Figure 1. Controller 72
may process inputs 71 as needed for the workflow in the model identification phase
41 and the model combination phase 42 to provide a reduced order model at an output
73 for an advanced process control. The reduced order model may go from output 73
to a controller 74 which may effect the advanced process control.
[0030] An industrial example may be a set of boilers feeding a common header. Typically,
fuel flow of all of the boilers may be operated by a common signal. Each boiler may,
with significant simplification, be described by internal boiler volume and hydrodynamic
pipe resistance between the boiler and a common header. Figure 3 is a diagram of a
boiler block scheme. There may also be dynamics from fuel flow 53 to generated steam
54, which can be "normalized" by a local combustion controller (to avoid oscillations
/ pushing among the parallel units). An output 55 of flow to a header may provide
negative feedback at a junction 56. Resulting drum pressure 57 may be fed to a junction
58 along with a negative addition of header pressure 59. The result may be the flow
55 to the header.
[0031] Low order models may be obtained by performing experiments for virtually all possible
on/off combinations, and then fitting a reduced order model, or by using certain heuristics
for parallel models reduction. Performing experiments for all combinations is not
necessarily practical as the number of combinations may typically exceed tens / hundreds
for larger solutions.
[0032] Model reduction in the related art may be done by the following heuristic (based
on first principles). The parallel boilers may be replaced by a single boiler with
normalized fuel flow to steam flow dynamics. Drum volume may be computed as a sum
of individual boiler drum volumes. Then, pipe resistance may be computed as a parallel
resistance of individual boiler pipe resistances.
[0033] Although this described heuristic may appear to work quite well for boilers; it may
be preferable to replace the heuristic by a systematical approach, such as a balanced
order reduction. A straightforward naive application of balanced reduction to parallel
models and integration of a reduced model to global model, may give a significantly
biased and even an unstable global model. An issue is that the "local" reduction of
parallel models should be done with respect to a global model.
[0034] An arrangement of parallel boilers may be noted as an example. The following simulations
may assume boilers feeding steam to a single header. Figure 4 is a diagram which shows
an N number of boilers connected to a header and indicates components of boiler flow.
It may be seen that a naive balanced reduction may have a large bias and tend to be
unstable. The structure-preserving algorithm may give consistent results with minimized
degradation to a full order model.
[0035] A linearized boiler model may be simulated as:

where
Vi i-th boiler volume
Ki i-th boiler pipe to header conductivity
Ks units of steam from unit fuel
Ti i-th boiler 1
st order time constant for fuel flow -> steam flow
[0036] A linearized header model may be simulated as:

where
VH header volume
[0037] The parameters may be chosen as:
Ks=10, V=(300 500 700), K=(100 130 150), VH=100.
[0038] The global model may have inputs Fuel Flow (FF) and Steam Demand (SD) and outputs
Drum pressure (p) and overall Steam Flow to header (SF).
[0039] Figures 5a-5d are diagrams of model step responses (FF - fuel flow, SD - steam demand).
Reduction with truncation (reduction to 3
rd order) may be noted. Figures 6a-6d are diagrams of reduced models step responses.
Figures 7a-7d are diagrams of step responses differences to original system. Reduction
with singular perturbations may be noted. Figures 8a-8d are diagrams of reduced models
step responses. Figures 9a-9d are diagrams of step responses differences to original
system.
[0040] The global model may be of a 7th order with step responses in the graphs of Figures
5a-5b. The target model order may be selected as 3. Figure 5a shows header pressure
(bar) versus time (sec) for fuel flow (t/hrs). Figure 5b shows header pressure (bar)
versus time (sec) for steam demand (t/hrs). Figure 5c shows stream flow to header
(t/hrs) versus time (sec) for fuel flow (t/hrs). Figure 5d shows steam flow to header
(t/hrs) versus time for steam demand (t/hrs).
[0041] The results for truncation are shown in Figures 6a-6d and differences to full order
model are shown in Figure 7a-7d. Similar comparisons for singular perturbations are
shown in Figures 8a-8d and 9a-9d.
[0042] Figure 6a shows a graph of header pressure (bar) versus time (sec) for fuel flow
(t/hrs), in full model, structured and balanced versions. Figure 6b shows a graph
of header pressure (bar) versus time (sec) for steam demand (t/hrs), in full model,
structured and balanced versions. Figure 6c show a graph of flow to header (t/hrs)
versus time (sec) for fuel flow (t/hrs), in full model, structured and balanced versions.
Figure 6d is a graph of flow to header (t/hrs) versus time (sec) for steam demand
(t/hrs), in full model, structured and balanced versions. Figure 6e is a graph of
scale for steam demand (t/hrs), in full model, structured and balanced versions.
[0043] Figure 7a is a graph of header pressure error (bar) versus time (sec) for fuel flow
(t/hrs), in global (dark dashed line), structured (solid line) and balanced (light
dashed line) versions. Figure 7b is a graph of header pressure (bar) versus time (sec)
for steam demand (t/hrs), in global (dark dashed line), structured (solid line) and
balanced (light dashed line) versions. Figure 7c is a graph of flow to header (t/hrs)
versus time (sec) for fuel flow (t/hrs), in global (dark dashed line), structured
(solid line) and balanced (light dashed line) versions. Figure 7d is a graph of flow
to header (t/hrs) versus time (sec) for steam demand (t/hrs), in global (dark dashed
line), structured (solid line) and balanced (light dashed line) versions.
[0044] Figure 8a is a graph of header pressure (bar) versus time (sec) for fuel flow (t/hrs),
in full model (thick line), structured (thin line) and balanced (dashed line) versions.
Figure 8b is a graph of header pressure (bar) versus time (sec) for steam demand (t/hrs),
in full model (thick line), structured (thin line) and balanced (dashed line) versions.
Figure 8c is a graph of flow to header (t/hrs) versus time (sec) for fuel flow (t/hrs),
in full model (thick line), structured (thin line) and balanced (dashed line) versions.
Figure 8d is a graph of flow to header (t/hrs) versus time (sec) for steam demand
(t/hrs), in full model (thick line), structured (thin line) and balanced (dashed line)
versions.
[0045] Figure 9a is a graph of header pressure error (bar) versus time (sec) for fuel flow
(t/hrs), in global (dark dashed line), structured (solid line) and balanced (light
dashed line) versions. Figure 9b is a graph of header pressure error (bar) versus
time (sec) for steam demand (t/hrs), in global (dark dashed line), structured (solid
line) and balanced (light dashed line) versions. Figure 9c is a graph of flow to header
error (t/hrs) versus time (sec) for fuel flow (t/hrs), in global (dark dashed line),
structured (solid line) and balanced (light dashed line) versions. Figure 9d is a
graph of flow to header error (t/hrs) versus time (sec) for steam demand (t/hrs),
in global (dark dashed line), structured (solid line) and balanced (light dashed line)
versions.
[0046] Honeywell, spol. s r.o, V Parku 2326/18, Praha 4, 14800, Czech Republic, may have
industrial boiler data for demonstrating algorithm efficiency.
[0047] Figures 5a-5d, 6a-6e, 7a-7d, 8a-8d, 9a-9d, 10a-10d, 11a-11d, 12a-12d and 13a-13d
may be noted. The group of Figures 5a-5d, 6a-6e, 7a-7d, 8a-8d and 9a-9d, may appear
similar to the group of Figures 10a-10d, 11a-11d, 12a-12d and 13a-13d; but are different
in the sense that the group of Figures 5a-5d, 6a-6e, 7a-7d, 8a-8d and 9a-9d appear
to show algorithm performance on artificial data and the group of Figures 10a-10d,
11a-11d, 12a-12d and 13a-13d appear to show results on data obtained from industrial
measurements.
[0048] Figures 10a-10d and 11a-11d are graphs of data of step responses and frequency responses,
respectively, of structure vs. unstructured reduction, for the example of three boilers
and a single header. A solid line, a dashed line, a dotted line and a dash-dot line,
represent original, structured, balanced and heuristic data plots, respectively. Figure
10a shows header pressure versus time (sec) for fuel flow. Figure 10b shows header
pressure versus time (sec) for steam demand. Figure 10c shows boilers to header flow
versus time (sec) for fuel flow. Figure 10d shows boilers to header flow versus time
for steam demand. Figure 11a shows magnitude (dB) versus frequency (rad/sec) for fuel
flow
〉̶ header pressure. Figure 11b shows magnitude (dB) versus frequency (rad/sec) for steam
demand
〉̶ header pressure. Figure 11c shows magnitude (dB) versus frequency (rad/sec) for fuel
flow
〉̶ boilers to header flow. Figure 11d shows magnitude (dB) versus frequency (rad/sec)
for steam demand
〉̶ boilers to header flow.
[0049] Figures 12a-12d and 13a-13d are graphs of data of step responses and frequency responses,
respectively, for an example of five boilers (no ID). A solid line, a dashed line,
a dotted line and a dash-dot line, represent original, structured, balanced and heuristic
data plots, respectively. Figure 12a shows header pressure (MPa) versus time (sec)
for fuel flow (t/hrs). Figure 12b shows header pressure versus time (sec) for steam
demand (t/hrs). Figure 12c shows flow to header (t/hrs) versus time (sec) for fuel
flow (t/hrs). Figure 12d shows flow to header (t/hrs) versus time (sec) for steam
demand (t/hrs). Figure 13a shows magnitude (dB) versus frequency (rad/sec) for fuel
flow (t/hrs)
〉̶ header pressure (MPa). Figure 13b shows magnitude (dB) versus frequency (rad/sec)
for steam demand (t/hrs)
〉̶ header pressure (MPa). Figure 13c shows magnitude (dB) versus frequency (rad/sec)
for fuel flow (t/hrs)
〉̶flow to header (t/hrs). Figure 13d shows magnitude (dB) versus frequency (rad/sec)
for steam demand (t/hrs)
〉̶ flow to header (t/hrs) .
[0050] In general, the present approach may provide a solution to modeling of plants with
parallel working units which may allow one to: 1) Significantly reduce time and resources
needed for identification of a plant model for different on/off configurations (a
number of required identification experiments (step test) is equal to the number of
parallel units); and 2) Improve results of currently used heuristics in quality and
mainly in quality consistency.
[0051] To recap, the present system for a reduced order model for advanced process control
may incorporate a plurality of parallel units having enablement devices, a remaining
system without the plurality of units but having an output and an input, and a summer
having a plurality of inputs connected to outputs of the plurality of units and having
an output connected to the input of the remaining system. An input of each unit of
the plurality of units may be connected to the output of the remaining system.
[0052] The system may incorporate a mechanism which provides an approach for applying a
structured model reduction on the system. The approach for applying the structured
model reduction may incorporate performing an identification experiment on each unit
of the plurality of units, computing parameters of a model for each unit from the
data of a corresponding identification experiment, computing parameters of a model
for the remaining system from data of an identification experiment, combining selected
models of the plurality of units in a closed loop with the model of the remaining
system to result in a whole system model, and performing a structured model reduction
on the on the whole system model. Selected models may be determined according to a
configuration of enablement devices. The structured model reduction may result in
a reduced order model for advance process control.
[0053] An instance of the system may be where each unit is a boiler, and the remaining system
has a header and interconnections. A switch or an enablement device of each unit may
be a valve. Other instances beside a boiler setup may be implemented within the present
system.
[0054] An approach for enabling parallel working units of a system for advance process control
may have a model identification phase and a model combination phase. The model identification
phase may incorporate determining models of parallel units of a system and a model
of a remaining system which has no units. The model combination phase may incorporate
a model of a whole system having the models of the parallel units and the model of
the remaining system.
[0055] Determining models of parallel units may incorporate performing an identification
experiment of each unit of the parallel units where the unit subject to the experiment
is the only one enabled in that the other units are disabled. Parameters may be computed
from data of the identification experiment for each unit, and a model of each unit
may be extracted from the parameters.
[0056] Likewise, determining a model of the remaining system may incorporate computing parameters
from data of an experiment and extracting the model of the remaining system from the
parameters. The model combination phase may incorporate a configuration that indicates
the enabled units and non-enabled units of the parallel units.
[0057] An approach may also incorporate computing a combination of models of units enabled
according to the configuration with the model of the remaining system, and applying
a structured model reduction to the combination of parallel models according to the
configuration and the model of the remaining system to obtain a reduced order model.
The reduced order model may be provided to an advanced process control controller.
[0058] An algorithm may facilitate the model combination phase. Various algorithms may be
used in the present approach. An algorithm may deal with models of each of the parallel
units, a model of the remaining system, and inputs of the configuration of the units.
An instance of the algorithm may incorporate combining models of parallel units to
a single state-space model according to the configuration, combining the state-space
model with the model of the remaining system to a single model of dynamics of the
whole system for the configuration, computing a controllability factor and an observability
factor for the single model of the whole system, separating the controllability factor
and the observability factor according to several dimensions, computing a Cholesky
factor of the controllability factor, computing a singular value decomposition, computing
a transformation matrix, applying the transformation matrix to get a reduced state-space
model of the parallel units as a truncation or a singular perturbation, and/or combining
the reduced state space model of the parallel units with the model of the remaining
system to get a final model for advanced process control.
[0059] To summarize, it may be noted that an approach for providing a reduced order model
for an advanced process control, may involve a model identification phase and a model
combination phase. The model identification phase may incorporate obtaining models
of each unit of a plurality of parallel units of a system and a model of a remaining
system generally having no unit of the plurality of parallel units. The model combination
phase may incorporate combining models of units that constitute a configuration and
the model of the remaining system to obtain a whole system model. The whole system
model may be subject to a structured model reduction to a reduced order model for
advanced process control. The providing a reduced order model for an advanced process
control may be embedded in the system.
[0060] In the present specification, some of the matter may be of a hypothetical or prophetic
nature although stated in another manner or tense.
[0061] Although the present system and/or approach has been described with respect to at
least one illustrative example, many variations and modifications will become apparent
to those skilled in the art upon reading the specification. It is therefore the intention
that the appended claims be interpreted as broadly as possible in view of the prior
art to include all such variations and modifications.
1. A system for a reduced order model for advanced process control, comprising:
a plurality of parallel units having enablement devices;
a remaining system, without the plurality of units, having an output and an input;
and
a summer having a plurality of inputs connected to outputs of the plurality of units,
and having an output connected to the input of the remaining system; and
wherein an input of each unit of the plurality of units is connected to the output
of the remaining system.
2. The system of claim 1, further comprising a mechanism containing steps for applying
a structured model reduction on the system.
3. The system of claim 2, wherein the steps for applying the structured model reduction
comprise:
performing an identification experiment on each unit of the plurality of units;
computing parameters of a model for each unit from the data of a corresponding identification
experiment;
computing parameters of a model for the remaining system from data of an identification
experiment;
combining selected models of the plurality of units in a closed loop with the model
of the remaining system to result in a whole system model; and
performing a structured model reduction on the on the whole system model; and
wherein selected models are determined according to a configuration of enablement
devices.
4. The system of claim 3, wherein the structured model reduction results in a reduced
order model for advance process control.
5. The system of claim 4, wherein:
each unit is a boiler;
the remaining system comprises a header and interconnections; and
an enablement device of each unit is a valve.
6. An approach for enabling parallel working units of a system for advance process control,
comprising:
a model identification phase; and
a model combination phase; and
wherein:
the model identification phase comprises determining models of parallel units of a
system and a model of a remaining system which has no units; and
the model combination phase comprises a model of a whole system incorporating the
models of the parallel units and the model of the remaining system.
7. The approach of claim 6, wherein determining models of parallel units comprises:
performing an identification experiment of each unit of the parallel units, when enabled
while the other units are disabled;
computing parameters from data of the identification experiment for each unit; and
extracting a model of each unit from the parameters.
8. The approach of claim 7, wherein:
determining a model of the remaining system comprises computing parameters from data
of an experiment and extracting the model of the remaining system from the parameters;
and
the model combination phase comprises a configuration of enabled units and non-enabled
units of the parallel units.
9. The approach of claim 8, further comprising:
computing a combination of models of units enabled according to the configuration
with the model of the remaining system; and
applying a structured model reduction to the combination of parallel models according
to the configuration and the model of the remaining system to obtain a reduced order
model.
10. The approach of claim 9, further comprising providing the reduced order model to an
advanced process control controller.
11. The approach of claim 9, wherein:
an algorithm facilitates the model combination phase; and
assumptions for the algorithm comprise:
models of each of the parallel units;
a model of the remaining system; and/or
inputs for the algorithm comprising the configuration of the units.
12. The approach of claim 11, wherein the algorithm comprises:
combining models of parallel units to a single state-space model according to the
configuration;
combining the state-space model with the model of the remaining system to a single
model of dynamics of the whole system for the configuration;
computing a controllability factor and an observability factor for the single model
of the whole system;
separating the controllability factor and the observability factor according to several
dimensions;
computing a Cholesky factor of the controllability factor;
computing a singular value decomposition;
computing a transformation matrix;
applying the transformation matrix to get a reduced state-space model of the parallel
units as a truncation or a singular perturbation; and/or
combining the reduced state space model of the parallel units with the model of the
remaining system to get a final model for advanced process control.
13. A method for providing a reduced order model for an advanced process control, comprising:
a model identification phase; and
a model combination phase; and
wherein:
the model identification phase comprises obtaining models of each unit of a plurality
of parallel units of a system and a model of a remaining system which has no unit
of the plurality of parallel units; and
the model combination phase comprises combining models of units that constitute a
configuration and the model of the remaining system to obtain a whole system model.
14. The method of claim 13, wherein the whole system model is subject to a structured
model reduction to a reduced order model for advanced process control.
15. The method of claim 14, wherein a providing a reduced order model for an advanced
process control is embedded in a system.