Open access

Lower Order Modeling and Control of Alstom Fluidized Bed Gasifier

Written By

L. Sivakumar and X. Anithamary

Submitted: 13 December 2011 Published: 24 October 2012

DOI: 10.5772/48674

From the Edited Volume

Gasification for Practical Applications

Edited by Yongseung Yun

Chapter metrics overview

2,560 Chapter Downloads

View Full Metrics

1. Introduction

Integrated coal gasification combined cycle (IGCC) is a technology wherein coal is converted to fuel gas also referred as syngas or synthesis gas. Powdered coal is made to be in contact with a mixture of oxygen(or air) and steam to produce fuel gas. This fuel gas is burnt in a gas turbine coupled with generator to produce power. The waste heat from the gas turbine is used to produce steam and the steam is sent to a steam turbine for additional power generation (Ramezan and Stiegel, 2006).

Though, IGCC has a number of technical advantages, but until recently, its application has been limited due to its higher capital costs plus the availability of cheap natural gas. However, with pollution limits becoming more stringent and natural gas prices increasing, the performance of IGCC will become more attractive and its technical advancement will further reduce its cost.

Gasification is a technology that had its beginnings in the late 1700s. In the 19th century, gasification was widely used for the production of “town gas” especially for urban areas (Ramezan and Stiegel, 2006). But due to the widespread availability of natural gas, it got vanished in the 20th century. Today, the IGCC technology is being widely used throughout the world. 250MW IGCC demonstration plants are being constructed at Tianjin in china. In India, Andhra Pradesh Power Generation Corporation Ltd in association with Bharat heavy Electricals Limited proposed 125 MW IGCC plant at Vijayawada. In USA, 262 MW Wabash River IGCC power plants in Indiana (later acquired by Conoco Philips) and 250MW Tampa Electric Co. Polk Power Station IGCC in Florida (later acquired by GE Energy) are the two main commercial IGCC coal based power plants. Even though a number of IGCC projects exist, the UK’s Clean Coal Power Generation Group, ALSTOM has undertaken a detailed study on the development of a small-scale prototype integrated plant (PIP), based on the air blown gasification cycle with 150 MW output (Pike et al., 1998). This type of prototype plant is useful in understanding the physics of the process, designing control systems for integrated operation.

Advertisement

2. Mathematical modelling

In general, mathematical modeling has been a useful tool for performance analysis, control system design, optimization and diagnosis of plants [Sivakumar and Ganapathiraman 2006]. The approach towards mathematical modeling depends upon the purpose for which the modeling is done. A detailed nonlinear mathematical model for a power boiler had been developed [Sivakumar and Bhattacharya 1979] using first principles approach – conservation of mass, energy and momentum to study the boiler transients for different types of disturbances. A furnace model with detailed calculations on the heat flux falling on different zones of furnace had been developed to study on the water wall tube failures [Sivakumar et.al 1980]. Low order transfer function models for power plant had been developed to study the performance of the proposed controllers and to design training simulators [Sivakumar et.al 1983]. This chapter deals with the development of low order mathematical models for ALSTOM gasifier which will be available to research community to study the efficiency of different control algorithms for specified disturbances. Further the suitability of conventional PID controllers for ALSTOM gasifier is investigated by the authors.

Advertisement

3. Air blown gasification cycle

ABGC is a hybrid combined cycle power generation technology. It was first conceived by British Coal Corporation (BCC) and developed in 1990s by Clean Coal Power Generation Group (CCPGG). Later the ABGC technology is purchased by Mitsui Babcock Energy Limited (Mitsui Babcock). Advanced design for this gasification is later done by the combined industrial collaborators - GEC Alsthom, Scottish Power plc and Mitsui Babcock with support from the European Commission’s (EC’s) THERMIE Programme and Department of Trade and Industry (DTI) (Pike et al., 1998). Figure 1 shows the block diagram of ABGC.

Coal, steam and air react within the gasifier operating at 22bar pressure and 1150k temperature conditions in order to produce fuel gas with low calorific value. Limestone is also added in order to remove sulphur. This fuel gas is burnt in a gas turbine coupled with generator to produce electricity.

Approximately 20% of carbon in the coal does not react in gasifier which is extracted through ash removal system. This unburned carbon is fed to circulating fluidized bed combustor (CFBC) operating under atmospheric pressure and 1150k temperature conditions. Here the remaining unburned carbon is combusted completely. The water/steam (two phase mixture) absorbs heat from CFBC water walls. The steam separated by drum internals goes through different stages of super heaters receiving heat from exhaust gas coming from gas turbine (Pike et al., 1998). The resulting high pressure steam is given to steam turbine coupled with generator to produce additional power generation. The total capacity of commercial ABGC is 525 MW approximately.

Figure 1.

The Air Blown Gasification Cycle

Advertisement

4. Types of gasifier

There are three types of gasifier namely fixed bed, fluidized bed and entrained flow (Phillips, 2006).

4.1. Fixed bed gasifier

Here coal enters at the top of the reactor and air or oxygen enters at the bottom. As the coal moves slowly down the reactor, it is gasified and the remaining ash drops are collected at the bottom of the reactor. Example: British Gas Lurgi(BGL), Lurgi (Dry Ash) The figure 2 shows moving bed gasifier.

Figure 2.

Moving bed gasifier

4.2. Entrained flow

Finely-ground coal is injected in co-current flow with the oxidant. The coal rapidly heats up and reacts with the oxidant. Gas is collected at the bottom. Most entrained flow gasifiers use oxygen rather than air. Example: GE entrained flow gasifier(Polk Station), E-Gas, Mitsubish Figure 3 shows entrained flow gasifier.

Figure 3.

Entrained Flow Gasifier

4.3. Fluidized bed gasifier

A fluidized bed gasifier is a well-stirred reactor in which new coal particles is mixed with older, partially gasified and fully gasified particles. The mixing gives uniform temperatures throughout the bed. The flow of gas into the reactor (oxidant, steam, recycled syngas) must be sufficient to float the coal particles within the bed. However, as the particles are gasified, they will become smaller and lighter and will be entrained out of the reactor. Example: HT Winkler, KRW (Kellogg –Rust-Westinghouse) and ALSTOM gasifier.

Figure 4.

Fluidized bed gasifier

Advertisement

5. ALSTOM gasifier model

Gasifier model is the most complex one in coal gasification. It was first started by CRE Group Ltd in 1992. Later it was continued at GEC ALSTHOM mechanical Engineering Centre. The incoming coal is dried and de-volatilized to yield char, ash and volatile gases. The oxygen in fluidized air reacts with carbon in the char to form carbon monoxide and carbon dioxide. Both exothermic and endothermic reactions occur simultaneously in the gasifier. The main equations in gasifier are

C+O2CO2E1
C+12O2 COE2

Equation 1 and 2 are exothermic gasification.

The carbon-dioxide reacts more with carbon to form carbon-monoxide. Also steam reacts with carbon to form carbon-monoxide and hydrogen.

C+CO22COE3
C+H2OCO+H2E4

Equation 3 and 4 are endothermic reactions.

The un-reacted char is added to the bed which is maintained at a constant height by char extraction system.

5.1. Alstom gasifier: Input and output variables

Alstom gasifier represents a difficult process for control because of its multivariable and non-linearity in nature with significant cross coupling between the input and output variables (Dixon 2004).

The controllable input variables to the gasifier are

  • Char off-take (u1)WCHR(kg/s)

  • Air flow rate(u2)WAIR(kg/s)

  • Coal flow rate(u3)WCOL (kg/s)

  • Steam flow rate(u4)WSTM(kg/s)

  • limestone flow rate (u5)WLS(kg/s)

The Controlled output variables are:

  • Gas calorific value (y1)CVGAS(J/kg)

  • Bed mass (y2)MASS(kg)

  • Fuel gas pressure (y3)PGAS(N/m2 )

  • Fuel gas temperature (y4)TGAS(K)

One of the inputs, limestone mass (WLS) is used to absorb sulphur in the coal and its flow rate is set to a fixed ratio of 1:10 against another input coal flow rate.(WCOL).This leaves effectively 4 degrees of freedom for the control design. Fig 5 shows gasifier with input and output variables.

Figure 5.

Gasifier with input and output variables

5.2. Load demand on gasifier

The flow rate of syngas to gas turbine is controlled through a valve at the inlet of turbine (also referred as controlled input disturbance to the gasifier). The pressure at the inlet of turbine called as PSink is the controlled variable. The control problem is to study the transient behavior of gasifier process variables such as pressure, temperature of the syngas for typical variations in gas flow drawing rate to gas turbine through appropriate changes in the throttle valve. Any proposed control system should control the pressure and temperature of the syngas at the inlet of gas turbine for any variation in gas turbine load – which in turn will affect throttle valve moment-without undue overshoots and undershoots. In fact this particular aspect has been posed as a control challenge problem for gasifier by ALSTOM.

Advertisement

6. ALSTOM benchmark challenges

The demand for clean air and stringent environmental regulations are forcing us to look for an alternate technology with reduced pollution emission and higher power generation. As a result of this, IGCC power plants are being developed all over the world. ALSTOM small-scale prototype (PIP) based on air-blown gasification cycle is one such IGCC. One of the component in ABGC called gasifier, is difficult to be controlled. For this reason, ALSTOM Power technology center issued a bench mark challenge to research community

  • To come out /propose a suitable control strategy/algorithms so as to have an efficient control of pressure and temperature of syngas without having an undue overshoot and undershoot values equal or less than those specified in the constraints by ALSTOM for specified load disturbance through the throttle value for different operating loads such as 100%, 50%and no-load.

The ALSTOM gasifier is modeled in state space form given by

X ˙= Ax+BuE5
Y=Cx+DuE6

Where

x = Internal states of gasifier, a column vector with dimension 25x1

u = Input variables, a column vector with dimension 6x1

A = system matrix governing the process dynamics, a square matrix with dimension 25x25

B = Input matrix with dimension 25x6

Y = Output variables, a column vector with dimension 4x1

C = Observable matrix with dimension 25x4

D = disturbance matrix with dimension 4x6

Towards this purpose, ALSTOM has made it available the following :

  • A, B, C, D, x(0), Y for three different loads- 100%, 50% and no-load.

A virtual gasifier mathematical model is made available with the above quantities (http://www.ieee.org/OnComms/PN/controlauto/benchmark.cfm.) and researches can attempt different control philosophies to meet the challenge posed by ALSTOM.

The input and output variables, allowable limits on output variables during load transients for three different loads (100%, 50% and no-load) as given by ALSTOM are reproduced in Tables 1 and 2 for ready reference.

6.1. Input and output constraints

The plant inputs and outputs with their limits are given in Tables 1 and 2 respectively (Seyab et al., 2006)

InputsDescriptionMaximum
Value
RateSteady state values
100%50%0%
WCHR(kg/s)Char extraction flow rate3.50.2 kg/s20.90.890.5
WAIR (kg/s)Air flow rate201.0 kg/s217.4210.894.34
WCOL(kg/s)Coal flow rate100.2kg/s28.555.342.136
WSTM(kg/s)Steam flow rate6.01.0kg/s22.701.690.676
WLS(kg/s)Limestone flow rate1.00.02kg/s20.850.530.21

Table 1.

Input Variables and Limits

OutputsDescriptionAllowed fluctuationsSteady state values
100%50%0%
CVGAS(MJ/kg)Fuel gas calorific value± 0.014.364.494.71
MASS(kg)Bedmass± 500100001000010000
PGAS(N/m2)Fuel gas pressure± 1×1042×1061.55×1061.12×106
TGAS(K)Fuel gas temperature± 1.01223.21181.11115.1

Table 2.

Output variables and limits

6.2. Researchers attempt in the first phase (1997-2001)

The first round challenge was issued in the year 1997. It included three linear models operating under 0%, 50%and 100% load conditions respectively. The model includes state space equation with A,B,C and D values. The challenge requires a controller which controls the gasifier at three load conditions with input and output constraints in the presence of step and sinusoidal disturbances. Many controllers have been suggested for the first challenge (Dixon, 1999).

  1. Dixon (1999) used multivariable P and I controllers using multi-objective optimal tuning technique and model based predictive control design to meet the constraints.

  2. Rice et al. (2000) proposed predictive control that uses linear quadratic optimal inner loop and it is supervised by an outer predictive controller loop.

  3. Proportional integral plus (PIP) by Taylor et al. (2000) from Lancaster University was based on discrete time model of the plant.

  4. Prempain et al. (2000) demonstrated the use of loop shaping H-infinity control design method.

  5. The multi-objective Genetic algorithm (MOGA) was proposed by Griffin et al. (2000) which performed a loop-shaping H-infinity design.

  6. A sliding mode, nonlinear design approach was suggested by Sarah Spurgeon. Here switching surface is designed to move the plant from one operating point to the other.

  7. Neil Munrom decomposed the original problem into a series of much simpler schemes in an effort to divide and conquer rule.

  8. Munro (2000) combined sequential loop closing with a high –frequency decoupling approach along with divide and conquer method

But none of the controller met all the objectives specified in the challenge – more so with particular reference to the transient limits imposed on output variables during load variations.

6.3. Second challenge

The second round challenge was issued in the year 2002. In the second round challenge, ALSTOM specified nonlinear simulation model in MATLAB/SIMULINK [10] and desired the controller capability during load changes and coal quality disturbance. Recently, a group of control solutions for the benchmark problem were presented at Control-2004 Conference at Bath University, UK in September 2004. Most of controllers were reported as capable of controlling the system at disturbance tests.

The author, Dixon (2002) used multi-loop PI controller to the gasifier control. He used system identification technique to obtain the linear model from the non – linear plant data. The base line controller was used by the other researchers for comparison purposes. The following controllers were suggested to meet the performance criteria (Dixon, 2004).

  1. Multi objective optimization approach suggested by Anthony Simms from Nottingham University needs further improvement by the addition of proportional control loops.

  2. H-infinity design approach given by Sarah Gatley from Leicester University used loop shaping combined with anti-windup compensator. It produced a robust design because of its simple design process and without the need for detailed knowledge of the plant.

  3. Multiple PID controller design using penalty based multi objective genetic algorithms by Adel Farag from Technical University of Hamburg gave excellent results that satisfied reasonable input output constraints.

  4. A novel controller by Tony Wilson from Nottingham University used state estimators to improve on the base line performance. Kalman filters are used to estimate the pressure disturbance and coal quality change.

  5. Proportional integral plus controller by James Taylor of Lancaster University used discrete time linear model of the gasifier.

  6. Model Predictive controller using a linear state space model of the plant was a collaborative effort from Cranfield and Loughborough.

All the papers had achieved reasonable success in terms controlling the gasifier model. But none of the controller met the overall performance criteria and still this benchmark challenge is left for the academicians for further research.

The difficulty in meeting the performance criteria appears to necessarily work with the higher order model for control system design. This motivates the authors to derive low order transfer function models for control system study.

Advertisement

7. Low order transfer function models

On analyzing the ALSTOM gasifier model, the model is found to be more complex and it contains very high cross-coupling between input and output (Dixon 2004). It necessitates low order model for further control research. The state space equation is converted to transfer function models using MATLAB command sys = ss(a,b,c,d) and [num,den]=ss2tf(a,b,c,d,1). After conversion by Matlab command, the system is described in s- domain as follows:

y1sy2sy3sy4s= G11s G12s G13s G14sG21s G22s G23s G24sG31s G32s G33s G34sG41s G42s G43s G44(s)u1su2su3su4(s)+
Gd1sGd2sGd3sGd4(s)E7
Psink

where

yi(s) = output variables ; i={1,4 }

Gij(s) = transfer characteristic between jth output due to ith input ; i= {1,4} j={1,4}

ui(s) = input variable ; i={1,4}

Gdi(s) = describing the impact of variation in Psink on output variable; i= {1,4}

Psink = sink gas pressure at gas turbine inlet.

It is to be noted that the denominator polynomial of each element Gij is of 24th order while the numerator is of order less than or equal to 23rd. A typical transfer characteristic between an output (pressure) due to all inputs shown diagrammatically as follows:

Figure 6.

Transfer characteristic between pressure due to all inputs

Here ∆Pui is the incremental change due to different inputs ui. Thus

∆Pu1 is the incremental change in pressure due to steady state change in char extraction flow rate,

∆Pu2 is the incremental change in pressure due to steady state change in Air flow rate,

∆Pu3 is the incremental change in pressure due to steady state change in Coal flow rate and

∆Pu4 is the incremental change in pressure due to steady state change in steam extraction flow rate. The output is given below

P(t) = Psteady state+∆Pu1 + ∆Pu2 + ∆Pu3+ ∆Pu4

Now the problem boils down to the reduction of higher order transfer function models obtained by MATLAB command to lower order transfer function models by the application of different methods.

It is observed that author Haryanto et al. (2009) developed an equivalent lower order transfer function models towards the development of integrated plant simulator. In this chapter, the authors have developed lower order transfer function models using algebraic and reduced order approximation methods (Sivakumar and Anithamary, 2011).

7.1. Reduced order approximation (RSYS)

The matlab command RSYS = BALRED(SYS,ORDERS) computes a reduced order approximation(RSYS) of LTI system. The desired order (number of states) is specified by ORDERS. BALRED uses implicit balancing techniques to compute the reduced-order approximation RSYS. The second order transfer function is obtained using Henkel Singularity approximation method. The transfer function for typical block G11 corresponding to 100% load is given below:

G11= -1.197e004 s2 +330.4 s+0.001125s2+0.0008608 s+2.075e-007E8

All the transfer function blocks Gij : (i = {1,4},j={1,4}) evaluated using reduced order approximation by the authors corresponding to 100%, 50% and no-load are given in Appendix A, Appendix B and Appendix C.

7.2. Algebraic method

The higher order transfer function is equated with the lower order model:

an-1sn-1 +an-2sn-2++a0bnsm +bn-1sm-1++b0=A2S2 +A1s+A0 B2S2 +B1s+BA0E9

On cross multiplying, the equation becomes

an-1sn-1 +an-2sn-2++a0B2S2 +B1s+B0= bnsm +bn-1sm-1++b0A2S2 +A1s+A0E10

The ALSTOM transfer function for G11 is given below

G11=0.03215s22+14.45s21+1289s20+1467s19+721.8s18+208.5s17+40.02s16+5.452s15+0.547s14+0.04129s13+0.002367s12+0.0001031s11+3.388e-006s10+8.235e-008s9+1.435e-009s81.695e-011s7+1.246e-013s6+4.869e-016s5+8.221e-019s4+6.532e-022s3+2.442e-025s2+3.437e-029 s+1.126e-034s24+35.38s23+78.31s22+68.51s21+32.81s20+9.998s19+2.106s18+0.3225s17+0.03703s16+0.00325s15+0.0002203s14+1.156e-005s134.687e-007s12+1.45e-008s11+3.36e-010s10+5.64e-012s96.52e-014s8+4.785e-016s7+1.95e-18s6+3.99e-021s5+4.505e-024s4+2.982e-024s3+1.148e-030s2+2.389e-034s+2.078e-038E11

The a0 can be obtained by the formula (Poongodi et al., 2009)

ao=bm-1/bm ±an-2/an-1m±nE12
ao=35.38/1 ±14.45/0.0321524±22E13

a0 = 10.5403, 242.4178, -9.0014, -207.0325

Taking the appropriate value of a0, equating the powers of s, and solving the equation, the unknown values of B0,B1,B2,A1,A2 can be obtained. Thus,

G11=-43.210273s2-32.8849432314s+10.5403-0.0083690166s2+0.067824414s+0.0019433E14

Similarly lower order models G12 to G44 corresponding to higher order models specified by ALSTOM can be obtained.

All the transfer function blocks Gij : (i = {1,4},j={1,4}) evaluated using algebraic method by the authors corresponding to 100%, 50% and no-load are given in Appendix A, Appendix B and Appendix C.

In order to evaluate the reduced order transfer function models obtained through different methods, the unit step response of ALSTOM model has been taken as reference response and the responses obtained through different methods as in figure 7 are compared and shown in figures 8-11 for typical transfer function blocks namely

G11 – the transfer characteristic between change in calorific value due to change in char extraction flow rate.

G24 – the transfer characteristic between change in temperature due to change in air flow rate.

G33 – the transfer characteristic between change in pressure due to change in coal flow rate.

G42 – the transfer characteristic between change in bedmass due to change in steam flow rate.

Figure 7.

Matlab SIMULINK model to evaluate the IAE and ISE error

Figure 8.

Variation of calorific value(y1) with char extraction flow rate (u1) keeping u2,u3,u4 constant

Figure 9.

Variation of fuel gas temperature(y4) with air flow rate (u2) keeping u1,u3,u4 constant

Figure 10.

Variation of fuel gas pressure(y3) with coal flow rate (u3) keeping u1,u2,u4 constant

Figure 11.

Variation of Bed mass(y2) with change in steam flow rate (u4) keeping u1,u2,u3 constant

The errors on the basis of IAE (Integral Absolute Error) and ISE (Integral Squared Error) are computed for each transfer function block obtained by algebraic method, reduced order approximation and RGA loop pairing over a period of time (little above the rise time) are shown in Table 3 for 100% load.

TransferfunctionINTEGRAL ABSOLUTE ERRORINTEGRAL SQUARED ERROR
Algebraic methodReduced order approximationTF using RGA loop pairingAlgebraic methodReduced orderapproximationTF using RGA
loop pairing
G1116441.062e+0051.087e+0042.16e+0061.133e+0091.455e+007
G127.09751.52.954e+0057.6061.013e+0051.12e+010
G134.828e+0044.48e+0048.039e+0047.98e+0087.98e+0087.784e+008
G145.09688.352.308e+0055.8510336.955e+009
G212.868e+0055.23e+0068.71e+0041.157e+103.598e+128.637e+009
G2211.51.19e+00420.9720.742.549e+00757.78
G2350.564.638e+0046.8e+00410182.668e+0085.145e+008
G2473.0976.29114.21412830.32555
G319.128e+0066.606e+0068.799e+0062.166e+131.009e+132.519e+013
G320.402117476.277e+0040.03623.051e+0054.58e+008
G3335.041.78e+0059250283.14.443e+0099.1e+006
G342.549141.81.086e+0050.859836221.344e+009
G411.437e+0072.407e+0071.434e+0078.005e+131.411e+0147.98e+013
G4215.181.103e+0042.69539.140.16321.213
G43462.35.714e+0041.133e+0053.812e+0045.035e+0081.46e+009
G441.683508.80.49940.33584.662 e+0040.1532

Table 3.

Integral Absolute and Squared error criteria for 3 models

It is observed that the low order models derived using algebraic methods is much superior to one proposed by Haryanto et.al., using RGA loop pairing and reduced order approximation proposed by authors.

7.3. Lower order modeling using genetic algorithm

Out of 16 transfer functions using algebraic method, four transfer functions G21, G31, G41 and G13 (shown in bold) are found to have higher ISE and IAE error criterion than the lower order models obtained using RGA loop pairing. This observation has motivated the authors to obtain further reduced order transfer function models with minimum ISE and IAE error criterion using genetic algorithm. Appendix D gives the auxiliary scheme for low order model (Sivanandam and Deepa, 2009).

The ALSTOM higher order transfer function for G13 is given below:

G13=-1.1s23-12.44s22+2893s21+2764s20+1212s19+324.2s18+59.29s17+7.864s16+0.7827s15+0.05961s14+0.003508s13+0.00016s12+5.623e-006s11+1.505e-007s10+3e-009s9+4.293e-011s84.155e-013s7+2.4446e-015s6+6.972e-018s5+4.046e-021s4+8.036e-024s3-.1.1772e-026s2-5.48e-030 s+8.511e-034s24+35.38s23+78.31s22+68.51s21+32.81s20+9.998s19+2.106s18+0.3225s17+0.03703s16+0.00325s15+0.0002203s14+1.156e-005s134.687e-007s12+1.45e-008s11+3.36e-010s10+5.64e-012s96.52e-014s8+4.785e-016s7+1.95e-18s6+3.99e-021s5+4.505e-024s4+2.982e-024s3+1.148e-030s2+2.389e-034s+2.078e-038E15

The second approximation is given as

G13=5.48e-30s -8.511e-0341.148e-30s2 + 2.389e-034s +2.076e-038E16

The transient and steady state gain for G13 is

TG/G13(s)=1.11 = 1.1E17
SSG/G13(s)= 8.511e342.076e38= 4.0997e+04E18

The auxiliary scheme given in appendix E is used to find R(s) from G(s)

R(s)=-5.48e-030 s+8.511e-0341.148e-030s2+2.389e-034s+2.078e-038E19

The above equation should be tuned to satisfy the transient and steady state gain so that R(s) reflects the characteristics of G(s)

R(s)=-1.1 s-7.4137631e-04s2 + 2.081e-04s+1.8083624e-08E20
=
B1s+B0b2s2+b1 s+b0E21

The parameters B0 = -7.4137631e-04, b1= 2.081e-04 and b0= 1.8083624e-08 are used as seed value for genetic algorithm with ISE error as the objective function. The ISE error (E) can be obtained by taking the sum of the square of the difference between the step response of higher and lower order transfer function. The ISE error is given by

E=t=0τ(Yt=yt)2E22

where, Yt is the unit step time response of the higher order system at the tth instant in the time interval 0t ≤τ, where τ is to be chosen and yt is the unit step time response of the lower order system at the tth time instant. The matlab commands

options =gaoptimset('InitialPop', [B1 B2 B3])[x fval output reasons] = ga(@objectivefun, nvars,options)E23

are used with ISE error as objective function. Here the population is set at 20 individuals and the maximum generation is 51. The crossover fraction is 0.8. Similarly the lower order models G31, G21 and G41 corresponding to higher order models specified by ALSTOM can be obtained. Table 4 shows the IAE and ISE error using genetic algorithm is further reduced than using algebraic method. Figure 12 shows the flowchart for lower order modeling using Genetic Algorithm.

Reduced transfer function using genetic algorithmB0b1b2IAE using genetic algorithmIAE using Algebraic methodISE using genetic algorithmISE using Algebraic method
G13= -1.1 s+2.8705s2 0.0922s+0.03392.8705-0.09220.03393.001e+0044.828e+0042.575e+0087.98e+008
G21= -9207 s+4.7874s2+ 350.5581-0.35514.7874350.5581-0.35518.718e+0042.868e+0058.634e+0091.157e+10
G31=6563 s+6.3912s2+ 0.0581s+1.8084e-086.39120.05811.8084e-088.57e+0069.128e+0062.442e+0132.166e+13
G41=-8868 s+2.3705s2- 0.2803s+0.09392.3705-0.28030.09391.399e+0071.437e+0077.782e+0138.005e+13

Table 4.

Reduced errors due to genetic algorithm in the evaluation of G13,G21,G31,G41

Figure 12.

Flowchart for lower order modeling using Genetic Algorithm

Transfer function obtained using Genetic Algorithm seems to be the most effective method for obtaining lower order models. Though the transfer functions for G11, G31,G41, G22 have been obtained through genetic algorithm to illustrate the superiority over other methods, all the transfer function blocks can be obtained in the same way as explained earlier.

Advertisement

8. Gasifier control and simulation

Even though many advanced control algorithms are proposed for complex process and systems, the authors are strongly of the opinion that PID control will also meet the control requirements using appropriate controller constants and feed forwards if necessary. Hence the PID controller is considered as a tool for gasifier control and simulation studies are done.

Figure 13.

PID controller for pressure and temperature output variables

Here PID controller is used to vary the steam and coal inputs for syngas pressure and coal and air is varied for syngas temperature. Table 5 gives the PID parameters for pressure and temperature of the syngas

P-Psink errorKpKiKd
PID(temperature)0.50.250.001
PID (pressure)7.543

Table 5.

PID constants for syngas temperature and pressure

Figure 14.

SIMULINK model for syngas pressure in the presence of step and sinusoidal disturbances

Figure 15.

SIMULINK model for syngas temperature in the presence of step and sinusoidal disturbances

Figure 16.

Syngas pressure maintaining at 2*106N/m2 in the presence of disturbance

Figure 17.

Syngas temperature maintaining at 1223K in the presence of disturbance

Advertisement

9. Conclusion

The development of low order transfer function model are required due to the difficulties encountered in the development of control strategies on ALSTOM benchmark challenge. In this direction, the authors have developed low order transfer function models using Algebraic method and reduced order approximation. The performance of these models has been evaluated on the basis of ISE and IAE error criteria. It is observed that the low order models derived using algebraic methods is much superior to one proposed by Haryanto et.al., and reduced order approximation. Some lower order transfer functions obtained using algebraic method are found to have higher error criterion than RGA loop pairing. Using Genetic Algorithm these errors are minimized and it is believed that the models proposed by algebraic method with Genetic Algorithm will become basis for further research on Gasifier control.

The authors have applied PID control algorithms for gasifier control around 100% load. As desired in the challenge problem, step and sinusoidal disturbances have been given in Psink. Preliminary simulation results show that the pressure and temperature of the syngas are controlled within the permissible constraint limits. However the authors intend to do extensive simulations for 100%, 50% and no-load with error due to pressure and temperature setpoints modulating different input variables.

Appendix A: Transfer function matrix of Alstom plant for 100% load

Transfer function blocksalgebraic methodreduced order approximation
G11-43.210273s2-32.8849432314s+10.5403-0.0083690166s2+0.067824414s+0.0019433 -1.197e004 s2 +330.4 s+0.001125s2+0.0008608 s+2.075e-007
G120.67268851s2+0.22784337s+1.367398.7409426609s2-6.32996277s-0.0002336 -3.468 s2- 1.063s - 0.001214s2+ 0.0008608s + 2.075e-007
G13-29.294957767s2+58.590928399s+0.99338 0.9053252009s2+0.217203375s-0.00002424 108.4 s2+ 6.901s - 0.008504s2+0.0008608s+ 2.075e-007
G1411.42165811s2-18.197458774s+2.892-298.17003810s2+56.7756422s-2.2915 -1.851s2+0.06763 s - 2.618e-007s2+0.0008608 s + 2.075e-007
G210.7699194835s2-0.4621252416s+1.4975-0.00005375534s2-0.0000999029s-3.530*10^-5-1.068e005 s2+43.57s- 0.008799s2+0.0008608 s + 2.075e-007
G220.1119962125834s2-0.335052778707s+1.5892-10.666078439387s2-3.7028097164s-1.016069*10^-3-71.72 s2-0.1481s-0.0003245s2+0.0008608 s + 2.075e-007
G23-38.1754867787s2-606.4765403s+1.03212-0.176233518s2-0.067401899s+0.00042351.142e004 s2-14.05s-0.0005055s2+0.0008608 s + 2.075e-007
G247.589742045s2+3.20126491848s+0.8050666.2192853528s2+13.8974102235s+0.391928.026s2+0.03455s+4.229e-006s2+0.0008608 s + 2.075e-007
G317.31943261016s2-83.3609061793s+0.76028-0.011722633s2-0.0005261888s+2.49106*10^-5 1.507e005s2-171.7s+ 0.01169s2+0.0008608 s + 2.075e-007
G32-8.34856920133s2+15.2823278158s+1.482526.0070991592s2+2.5943768447s+0.000366-175.2s2+0.4962s+0.0008401s2+0.0008608 s + 2.075e-007
G33-2.8969959854S2-269.839875362S+1.645-0.1417932867S2-0.05582682S+0.538723*10^-44288s2-4.413s-0.006334s2+0.0008608 s + 2.075e-007
G34-0.149422569149s2+0.605489884s+0.8755-13.277800585s2-15.075170803s-0.05380.9117s2+0.0606s- 3.372e-006s2+0.0008608 s + 2.075e-007
G410.989892606658S2-6.37153721233S+1.50060.0005803934141S2+0.000131142815S-3.06317*10^-5 1.941e005s2+48.24s-0.01016s2+0.0008608 s + 2.075e-007
G420.863152338637S2-1.69330243903S+2.3304-12.329066005S2-3.02684037S-0.00315996 1.941e005s2+48.24s-0.01016s2+0.0008608 s + 2.075e-007
G43-15.4009119304S2-2940.056236928S+2.31138-0.566750514100S2-0.1874561152S+0.00217711.709e004s2 +6.082s+0.0002203s2+0.0008608 s + 2.075e-007
G44201.4423617140688S2+275.7771961791S+0.81865-4192.317426968S2-1162.912156389S-0.017373.079s2-0.02195s-9.775e-006s2+0.0008608 s + 2.075e-007

Table 6.

Appendix B: Transfer function obtained of 50% load

Transfer function blocksalgebraic methodreduced order approximation
G1182973.885826s+1.3118521978.394989s2+1.308653s+1.568416e-05-1330s2+395.4s+0.006024s2+0.0005765s+7.203e-008
G12-4.151117e-06s2+0.182109e-03s+2.547485-2.550036s2-2.54957s-0.91584e-042.476s2-1.01s-0.0007864s2+0.0005765s+7.203e-008
G13727.803446s2-31191.578014s+1.2665199.782304s2-7.184683s+0.001618199.6s2+5.647s+5.637e-005s2+0.0005765s+7.203e-008
G14633.894001s2-26444.64477s+1.28063511207.461122s2-350.714185s+0.2167960.3227s2+0.06097s+4.256e-007s2+0.0005765s+7.203e-008
G21-7.638891e+18s+2.4953876.3235e+14s2+4.938082e+13s-7.626576e-05-1.839e005s2+56.37s-0.002357s2+0.0005765s+7.203e-008
G221.016135e-04s2+0.00445778s+2.676878-29.50425s2+3.382317s-10.2823689166s2-0.06141s-0.0001875s2+0.0005765s+7.203e-008
G23-150.691621s2-1.446301e+12s+1.682835-1.40965e+08s2+4.809870e+06s+0.0001578739s2+7.32s+0.0007719s2+0.0005765s+7.203e-008
G240.016439s2+1.190322e+08s+1.312122.249456e+07s2+4.0072857e+06s+0.0441638.171s2+0.02981s+2.14e-006s2+0.0005765s+7.203e-008
G31470334.30765s+3.27073445.137649s2-1.79416s-8.741614e-062.322e005s2-243.4s+0.001665s2+0.0005765s+7.203e-008
G323.716764e-07s2-160.82437e-07s+2.5489663.470826s2+3.7707s+3.13062e-04-375.9s2+0.4627s+0.0005865s2+0.0005765s+7.203e-008
G33-3.1067766e+05s+2.683919-189.668841s2-551.803398s+0.0058852926s2-3.597s+3.286e-005s2+0.0005765s+7.203e-008
G3459112.633568s+1.437609-2.6507907e+05s2-3511.145849s-0.053218-0.4876s2-0.0597s+3.286e-005s2+0.0005765s+7.203e-008
G41-2.720125e+07s+2.5011142381.895719s2+740.758982s-2.580143e-053.296e005s2+24.4s-0.006982s2+0.0005765s+7.203e-008
G422.01188e-03s2-0.087054s+3.70861304.883s2+1.097948e+05s-3.3664e-0472.72s2-0.4004s-7.935e-005s2+0.0005765s+7.203e-008
G435.456232e+04s+3.63205120.535311s2+20.261304s+0.0005591.258e004s2+4.893s+0.0004882s2+0.0005765s+7.203e-008
G441.2505305e+04s+1.333895-6380.258041s2-613.370677s-0.0204171.715s2-0.02503s-4.705e-006s2+0.0005765s+7.203e-008
Gd11.427343e+06s-1.570372e+06-8.105303e+14s2-2.361753e+11s-1.1699040e+070.1224s2-2.013e-005s+9.669e-009s2+0.0005765s+7.203e-008
Gd2-2.7855e-06s2+0.122333e-03s+1.2902941901.74673s2+41909.74673s+0.3277079e-039.213e-005s2+3.538e-007s+2.198e-010s2+0.0005765s+7.203e-008
Gd31.319354e+11s+2.2179032.175717e+11s2-3.741932e+11s+2.325120.9534s2+0.0005484s+6.87e-008s2+0.0005765s+7.203e-008
Gd42.3313028e+04s+1.2987352.220288e+08s2-1.65059e+09s-4.327659e+04-3.39e-005s2-3.341e-008s-2.163e-012s2+0.0005765s+7.203e-008

Table 7.

Appendix C: Transfer function obtained of 0% load

Transfer function blocksalgebraic methodreduced order approximation
G11 6412495.306104s+59.51538781552.782731s2+14.789253s+8.741323e-05 3.828e004 s2+561.7 s + 0.006739s2+ 0.0002741 s + 9.897e-009
G125.156462e-6s2+2.9510432e-4s+3.342709638-3.34349077s2-3.3438166s-8.398176626e-05 79.85 s2- 0.955 s - 0.0003939s2+ 0.0002741 s + 9.897e-009
G13203764.661731s2-10907532.587163s+56.60632-2370.734866s2-5009.316672s+0.038167 178.3 s2+1.338 s + 1.467e-005s2+ 0.0002741 s + 9.897e-009
G148.872773e+40s2-4.859238e+42s+1.656547-1.042016e+42s2-9.464071e+40s+0.032264 3.845 s2+0.05121s+5.082e-007s2+ 0.0002741 s + 9.897e-009
G21-2.966146e+10s+114.0214441762415.861223s2-382735.79834s+0.000915 -4.377e005 s2+120.2 s + 0.001232s2+ 0.0002741 s + 9.897e-009
G22-2.21638386e-03s2+0.1268436s+3.475304171-41.36659216s2+9.7254335s+3.685291183e-04 948.2 s2+0.117 s- 9.333e-005s2+ 0.0002741 s + 9.897e-009
G23-1.783467s2-4.004262e+09s+2.15279-392574.725689s2-671071.290709s+0.000361 4701 s2+1.852 s + 5.905e-005s2+ 0.0002741 s + 9.897e-009
G24-0.000362s2-573649.046541s+1.698815-81530.566687s2-7775425361s+0.028148 11.85 s2+0.02519 s +5.973e-007s2+ 0.0002741 s + 9.897e-009
G31470334.30765s+3.27073445.137649s2-1.79416s-8.741614e-06 5.378e005 s2-408.2 s- 0.003703s2+ 0.0002741 s + 9.897e-009
G322.136783444e-07s2-122.288e-07s+3.3348674.423179153s2+4.802113698s+3.116106e-05 -1552 s2+0.2945 s+ 0.0003176s2+ 0.0002741 s + 9.897e-009
G3367905.569056s+3.379449407.551672s2-1891.8753s-0.007252 1502 s2-0.8885 s-4.612e-006s2+ 0.0002741 s + 9.897e-009
G3477345.248573s+1.937428-472193.214731s2-1741.323615s-0.024306 -5.796 s2-0.05223 s-7.889e-007s2+ 0.0002741 s + 9.897e-009
G41-1.664227e+08s+3.26639110606.927662s2+1483.255839s-9.560048e-06 9.186e005 s2-25.51 s-0.003382s2+ 0.0002741 s + 9.897e-009
G421.52637013e-03s2-0.08735416s+4.8185565-2490.9933s2-247081.2502s-15.27197e-04 45.08 s2-0.4357 s-3.123e-005s2+ 0.0002741 s + 9.897e-009
G43187977.323085s+4.50594669.134727s2+119.112126s+0.001149 6900 s2+1.238 s +3.885e-005s2+ 0.0002741 s + 9.897e-009
G4458954.826124s+1.729767-26761.155753s2-1437.956917s-0.01438 -4.063 s2-0.02662 s-1.191e-006s2+ 0.0002741 s + 9.897e-009
Gd1-1.344731e+23s-1475658.5621212.994947e+31s2+3.271975e+26s+2.050176e+07 -0.1992 s2-7.057e-005 s-7.124e-010s2+ 0.0002741 s + 9.897e-009
Gd2-5.6523279e-06s2+0.323485e-03s+3.303296936.4763s2+179827.0165s+754.357 -0.0003971 s2+6.566e-008 s +4.334e-011s2+ 0.0002741 s + 9.897e-009
Gd3-7.490644e+11s+2.88116-6.847024e+11s2-2.019168e+12s+2.921795 0.9858 s2+0.0002702 s +9.76e-009s2+ 0.0002741 s + 9.897e-009
Gd4-40689.070969s+1.677564-1.634099e+08s2-9828353.140255s-264.279412 -1.935e-005 s2-6.714e-009 s-6.282e-014s2+ 0.0002741 s + 9.897e-009

Table 8.

Appendix D: Lower order Transfer function reduction

Consider an nth higher order system represented by its transfer function

G(s) = N(s)D(s) = i=0n-1Ai sii=0naisi

=
An-1 sn-1 + An-2 sn-2 + + A2s2+ A1 s+ A0 an sn + an-1 sn-1 + + a2s2+ a1 s+ a0 E24

First Order = A0a1s+ a0 (5)

Second order = A1+ A0a2s2+ a1s+ a0…(6)

n-1 order =An-2 sn-2 + An-3 sn-3 + + A2s2+ A1 s+ A0 an-1 sn-1+ an-2 sn-2 + + a2s2+ a1 s+ a0 (7)

Equations (5) through (7) gives the lower order model for higher order system G(s). For n higher order system, (n-1) lower order models can be formulated.

Advertisement

Acknowledgement

The authors would like to express their thanks to Dr. Yongseung Yun, the book editor on Gasification for his valuable suggestions. Further the authors are grateful to Karunya University and Sri Krishna College of Engineering and Technology, Coimbatore for the support and encouragement.

References

  1. 1. AsmarB. N. W. E.JonesJaWilson. A.processengineering.approachto.thealsotm.gasifierproblem.procInst.MechEng. I. J.Systemcontrol.Eng2000441452
  2. 2. DixonR.AlstomBenchmark.challengeI. I.controlof.NonlinearGasifier.model2002http://www.iee.org/omcomms/PN/controlauto/Specification_2 pdf
  3. 3. 31021232004Dixon.R Advanced gasifier control, computing and control engineering journal IEE vol 10 N0 3 pp 93-96, 1999
  4. 4. Griffin,I.A.,, P.Schroder, AJ Chipperfield and PJ Fleming multiobjective optimization approach proc, Inst.Mech Eng.I,J.System control Eng.,2000453468
  5. 5. Haryanto,A., P.siregar, D.Kurniadi and Keum-shik Hong, Development of Integrated Alstom gasification Simulator for implementation using DCS CS 3000 proceedings of the 17th world congress The international conference federation of automatic control, seoul, korea 2009
  6. 6. LiuG. P.DixonR.DaleyS.multiobjectiveoptimal.tuningproportional.integralcontroller.designfor.theAlstom.Gasifierproblem.procInst.MechEng. I. J.Systemcontrol.Eng2000395404
  7. 7. Mitchell,M., An introduction to Genetic Algorithm, Prentice-Hall of India,New Delhi, Edition: 2004
  8. 8. MunroN.EdmundsJ. M.KontogianeesE.StImpram. A.sequentialloop.closingapproach.tothe.Alstomgasifier.problemproc.InstMech.EngI. J.Systemcontrol.Eng2000427439
  9. 9. PhillipsJ.2006Different types of gasifiers and their integration with gas turbines. In: The Gas Turbine handbook. http://www.netl.doe.gov/technologies/coalpower/turbines/refshelf/handbook/1.2.1.pdf
  10. 10. Pike A.W., Donne M.S and Dixon.R, Dynamic modeling and simulation of the air blown gasification cycle prototype integrated plant in proceedings of the international conference on simulation, IEE publication 457,York university1998354361
  11. 11. PoongodiP.S.Victor Genetic algorithm based PID controller design forLTI system via reduced order model, International Conference on Instrumentation, Control & Automation ICA2009 October 20-22, 2009Bandung, Indonesia.
  12. 12. PrempainE.PostlethwaiteI.sunX. D.Robustcontrol.ofthe.gasifierusing. a.mixedH∞.approachproc.InstMech.EngI. J.Systemcontrol.Eng2000415426
  13. 13. Ramezan and Stiegel,2006Integrated coal gasification combined cycle. In: The Gas Turbine handbook. http://www.netl.doe.gov/technologies/coalpower/turbines/refshelf/handbook/1.2.pdf
  14. 14. RiceM.RossterJ.schurmansJ.Anadvanced.predictivecontrol.approachto.theAlstom.gasifierprolem.procInst.MechEng. I. J.Systemcontrol.Eng2000405413
  15. 15. Seyab,R.K., Y.Cao and S.H Yang, The second alstom benchmark challenge on gasifier control predictive control for the ALSTOM gasifier problem IEE proceedings on control theory, vol153, N03 May 2006.
  16. 16. Sivakumar,L., Anithamary.X A low order transfer function model for MIMO ALSTOM gasifier, International conference on process Automation, Control and computing, IEEE 2011.
  17. 17. SivanandamS. N. S. N.DeepaA.ComparativeStudy.UsingGenetic.AlgorithmParticleSwarm.Optimizationfor.LowerOrder.SystemModelling.InternationalJournal.ofthe.Computerthe.InternetManagementVol.3110September- December, 2009
  18. 18. TaylorC. J.MccabeA. P.youngP. C.chotaiA.Proportionalintegral.pluscontrol.ofAlstom.Gasifierproblem.procInst.MechEng. I. J.Systemcontrol.Eng2000469480
  19. 19. Wang,X., Ke Wu, Jianhong Lu, Wenguo Xiang, Non linear identification of Alstom gasifier based on wiener model International conference on sustainable power generation and supply pg 1-7 april 2009.
  20. 20. Sivakumar. L,“Performance analysis, diagnosis and and optimisation (PADO) for power plants”, Seminar on Power Plant Automation Concepts & Applications-by The Instrumentation Systems and automation Society- Bangalore, April 22, 2006.
  21. 21. Sivakumar.L and Ganpathiraman.G “Performance analysis diagnostics and optimisation in generation”, Conference on: IT Power- Improving Performance And Productivity,New Delhi ; Sep 2006.
  22. 22. SivakumarL.ReddyK. L.SundararajanN. “.Detailedcirculation.analysisto.determinethe. D. N. B.marginin.naturalcirculation.boilers”Proceedings.ofInternational.conferenceon.HeatMassTransfer.Hyderabad.Feb1..198018
  23. 23. PonnusamyP.SivakumarL.SankaranS.V.“V.“Low-order dynamic model of a complete thermal power plant loop”, Proceedings of the Power Plant Dynamics,Control and Testing Symposium,1198310
  24. 24. SivakumarL.and.BhattacharyaR. K. “.Dynamicanalysis.ofa.powerboiler.usinga.non-linearmathematical.model”Proceedings.ofsecond.symposiumon.PowerPlant.DynamicsControlsHyderabad. .Recordof.ProceedingsFeb14-1.19792129

Written By

L. Sivakumar and X. Anithamary

Submitted: 13 December 2011 Published: 24 October 2012