Optimisation of Marine Boilers using Model-based Multivariable ...
Optimisation of Marine Boilers using Model-based Multivariable ...
Optimisation of Marine Boilers using Model-based Multivariable ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
36 3. SUMMARY OF CONTRIBUTIONS<br />
The model was <strong>based</strong> on first principles <strong>using</strong> mass and energy balances for describing<br />
the temperature and pressure dynamics. As flow through and efficiency <strong>of</strong> the turbocharger<br />
was required in these balance equations, approximations <strong>of</strong> the turbocharger<br />
turbine and compressor maps partly <strong>based</strong> on physical insight were made. The oxygen<br />
was modelled by describing mole balances for the oxygen in the gas generator and furnace<br />
volumes. The resulting model was <strong>of</strong> sixth order and presented in descriptor form<br />
as:<br />
F(x) dx<br />
= h(x,u,d)<br />
dt<br />
(3.7a)<br />
y = g(x,u,d) (3.7b)<br />
where u = [ ˙mfu1, ˙mfu2] T is the vector <strong>of</strong> inputs being the fuel mass flow to the gas<br />
generator and the furnace. d = [Ta,Tfu,Tm] T is the vector <strong>of</strong> disturbances being<br />
temperatures <strong>of</strong> the inlet air, <strong>of</strong> the fuel and <strong>of</strong> the metal separating the hot flue gas<br />
and the water/steam in the boiler. x = [pgg,Tgg,ω,Tfn,xgg,O2 ,xfn,O2 ]T is the state<br />
vector being pressure in the gas generator, temperature in the gas generator, shaft speed,<br />
temperature in the furnace, and oxygen fraction in the gas generator and the furnace.<br />
Finally, y = [ ˙mfu, ˙ Q,xfn,O2 ]T is the output vector with ˙mfu = ˙mfu1 + ˙mfu2 and<br />
˙Q is the energy transferred to the metal walls. Expanding, (3.7a) reveals the model<br />
structure as:<br />
⎡<br />
⎢<br />
⎣<br />
f11 f12 0 0 0 0<br />
f21 f22 0 0 0 0<br />
0 0 f33 0 0 0<br />
0 0 0 f44 0 0<br />
0 0 0 0 f55 0<br />
0 0 0 0 0 f66<br />
⎤⎡<br />
⎥⎢<br />
⎥⎢<br />
⎥⎢<br />
⎥⎢<br />
⎥⎢<br />
⎥⎢<br />
⎦⎢<br />
⎣<br />
dpgg<br />
dt<br />
dTgg<br />
dt<br />
dω<br />
dt<br />
dTfn<br />
dt<br />
dxgg,O2 dt<br />
dxfn,O2 dt<br />
⎤<br />
⎡<br />
⎥ ⎢<br />
⎥ ⎢<br />
⎥ ⎢<br />
⎥ = ⎢<br />
⎥ ⎢<br />
⎥ ⎢<br />
⎦ ⎣<br />
h1<br />
h2<br />
h3<br />
h4<br />
h5<br />
h6<br />
⎤<br />
⎥<br />
⎦<br />
(3.8)<br />
where the elements fij and hi were found in the model derivation – see [Solberg et al.,<br />
2008c]. F is never singular, hence has a well defined inverse, so (3.7a) can be written<br />
as an ordinary differential equation: ˙x = f(x,u,d) = F −1 (x)h(x,u,d).<br />
Validation against measurements collected from the test setup shows good agreement<br />
between model and measuring data in terms <strong>of</strong> capturing the dynamical behaviour.<br />
However, in terms <strong>of</strong> stationary values these are not represented well by the model<br />
for other outputs than the oxygen level, see Figure 3.5. This was accepted for now<br />
as the main focus was on oxygen control, and constraints on internal states were not<br />
considered.<br />
The performance requirements for the controller <strong>of</strong> the burner were to deliver the requested<br />
fuel flow while keeping the oxygen percentage above 3% and maximise efficiency.<br />
No requirement regarding <strong>of</strong>f-set free tracking <strong>of</strong> the fuel flow reference was<br />
introduced meaning that differences between requested and actual fuel flow should be<br />
handled by e.g. including integral action in the outer controller.