Overcoming the shrink-and-swell effect in water level control strategy on industrial boiler-drum

Similar documents
Impact of the Sampling Period on the Design of Digital PID Controllers

Two Products Manufacturer s Production Decisions with Carbon Constraint

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

A General Thermal Equilibrium Discharge Flow Model

The pn junction: 2 Current vs Voltage (IV) characteristics

Search sequence databases 3 10/25/2016

Principles of Humidity Dalton s law

Voltage, Current, Power, Series Resistance, Parallel Resistance, and Diodes

Sliding Mode Flow Rate Observer Design

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

Finite element discretization of Laplace and Poisson equations

4 x 4, and. where x is Town Square

Numbering Systems Basic Building Blocks Scaling and Round-off Noise. Number Representation. Floating vs. Fixed point. DSP Design.

Numerical considerations regarding the simulation of an aircraft in the approaching phase for landing

Why is a E&M nature of light not sufficient to explain experiments?

A Sub-Optimal Log-Domain Decoding Algorithm for Non-Binary LDPC Codes

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Characteristics of Gliding Arc Discharge Plasma

A Propagating Wave Packet Group Velocity Dispersion

Math 34A. Final Review

Elements of Statistical Thermodynamics

2. Laser physics - basics

Pipe flow friction, small vs. big pipes

EXST Regression Techniques Page 1

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

4037 ADDITIONAL MATHEMATICS

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then

The following information relates to Questions 1 to 4:

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB

Status of LAr TPC R&D (2) 2014/Dec./23 Neutrino frontier workshop 2014 Ryosuke Sasaki (Iwate U.)

2/12/2013. Overview. 12-Power Transmission Text: Conservation of Complex Power. Introduction. Power Transmission-Short Line

2008 AP Calculus BC Multiple Choice Exam

A Prey-Predator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction

5.80 Small-Molecule Spectroscopy and Dynamics

A Control Strategy for Photovoltaic-Solid Polymer Electrolysis System Based on Surface Temperature of PV Panel

dr Bartłomiej Rokicki Chair of Macroeconomics and International Trade Theory Faculty of Economic Sciences, University of Warsaw

Learning Spherical Convolution for Fast Features from 360 Imagery

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation

Einstein Equations for Tetrad Fields

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS

INTRODUCTION TO AUTOMATIC CONTROLS INDEX LAPLACE TRANSFORMS

Exchange rates in the long run (Purchasing Power Parity: PPP)

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

CE 530 Molecular Simulation

VII. Quantum Entanglement

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

Recursive Estimation of Dynamic Time-Varying Demand Models

Statistical Thermodynamics: Sublimation of Solid Iodine

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

Extraction of Doping Density Distributions from C-V Curves

Design Guidelines for Quartz Crystal Oscillators. R 1 Motional Resistance L 1 Motional Inductance C 1 Motional Capacitance C 0 Shunt Capacitance

4. Money cannot be neutral in the short-run the neutrality of money is exclusively a medium run phenomenon.

Sec 2.3 Modeling with First Order Equations

WASTE HEAT RECOVERY USING DIRECT THERMODYNAMIC CYCLE

Dynamic Characteristics Analysis of Blade of Fan Based on Ansys

15. Stress-Strain behavior of soils

Inheritance Gains in Notional Defined Contributions Accounts (NDCs)

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule

Predictive Control of Batch Reactors

0 +1e Radionuclides - can spontaneously emit particles and radiation which can be expressed by a nuclear equation.

are given in the table below. t (hours)

KINEMATIC SOIL-STRUCTURE INTERACTION EFFECTS ON MAXIMUM INELASTIC DISPLACEMENT DEMANDS OF SDOF SYSTEMS

What are those βs anyway? Understanding Design Matrix & Odds ratios

The Open Economy in the Short Run

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot

In this lecture... Subsonic and supersonic nozzles Working of these nozzles Performance parameters for nozzles

Hydrogen Atom and One Electron Ions

Setting up the problem. Constructing a Model. Stability Analysis of ODE. Cheyne-Stokes. A Specific Case of the ODE 7/20/2011

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real.

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient

Supplementary Materials

Lecture #15. Bipolar Junction Transistors (BJTs)

The influence of electron trap on photoelectron decay behavior in silver halide

Addition of angular momentum

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

A Recent Approach to Repetitive Control Strategy for Induced Draft Fan

MULTIVARIABLE CONTROLLER TUNING

Brief Notes on the Fermi-Dirac and Bose-Einstein Distributions, Bose-Einstein Condensates and Degenerate Fermi Gases Last Update: 28 th December 2008

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor

Osmium doping of UAl 2. Department of Physics and Engineering Greenville, SC Department of Physics Gainesville, FL

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Abstract Interpretation: concrete and abstract semantics

Molecular Orbital Diagrams for TM Complexes with

Optimal Design Schemes of Centrifugal Pump Compound Impeller at Super-Low Specific Speed

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved.

Unfired pressure vessels- Part 3: Design

Chapter 13 Aggregate Supply

A nonequilibrium molecular dynamics simulation of evaporation

Nuclear reactions The chain reaction

Types of Transfer Functions. Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters

Coupled Pendulums. Two normal modes.

Introduction to Condensed Matter Physics

Studies of Turbulence and Transport in Alcator C-Mod Ohmic Plasmas with Phase Contrast Imaging and Comparisons with GYRO*

Evaluating Reliability Systems by Using Weibull & New Weibull Extension Distributions Mushtak A.K. Shiker

Direct Approach for Discrete Systems One-Dimensional Elements

CHAPTER 16 HW: CONJUGATED SYSTEMS

AS 5850 Finite Element Analysis

Transcription:

Standard Scintific Rsarch and Essays Vol1 (1): 019-024, Fbruary 2013 http://www.standrsjournals.org/journals/ssre Rsarch Articl Ovrcoming th shrink-and-swll ffct in watr lvl control stratgy on industrial boilr-drum makhai Lawrnc. Email: oboscos@yahoo.com Rcivd 07 January 2013; Accptd 21 Fbruary 2013 ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Abstract Powr gnration has bcom an incrasingly comptitiv arna. Th cost of oprating powr plants coms mostly from th ful bill, which runs in th billions of dollars annually. n ordr to minimis th ful bill and maximis pant fficincy, a plant s load-following capability must b optimizd, that is, to follow th dmand of powr closly. A grat dal of attntion must b givn to th controllr that rgulats watr lvl (and stam prssur) at th stam-gnrating unit th boilr-drum. Th control stratgy is complicatd by th non-linar and non-minimum phas charactristics of th boilr-drum. This papr illustrats our work in th optimization of load-following schm by applying a 3-lmnt controllr schm to rgulat th boilr-drum s watr lvl. Th controllr is thn cascadd with anothr fdback loop that rgulats stam prssur. Th rsults from th control schm hav shown considrabl improvmnt ovr a typical boilr-drum s watr lvl control stratgy. Ky words: Shrink-and-swll, control, boilr-drum s, optimization, boilr-drum. NTRODUCTON Th cost of running a MW coal-fird unit can vary btwn USD$75,000 pr day whil th annual ful bill is narly USD$5 billion (Waddington t al., 1987). Waddington t al. (1987) showd that mor than 99% of th cost in powr plant opration coms from ful-in this cas, coal. Th typical componnts of a powr plant ar furnac, boilr-drum, suprhatd and rhatr, and turbin units. Figur 4 illustrats th boilr-drum with th down comr and risr circulation tubs. Th cylindrical drums ar not hatd. Rathr, hat is supplid to th incoming watr in th risr tubs by dirct hat from th furnac gasss. Thr ar a larg numbr of risr tubs in th drum-down comr risr circulating loop in ordr to maximiz hat transfr. Th downcomrs ar largr in siz sinc no hat transfr taks plac. Flow around th circulating loop can b ithr natural du to prssur diffrnc or forcd with pumps. n th dsign of th circuit, it is vry important that sufficint circulation occurs at all tims. Th cost of plant opration is mainly on th ful consumption. n ordr to optimiz opration, powr gnration has to follow powr dmand vry closly, that is, load-following (Rs, 1997). Th complication in th load-following schm is mainly causd by th shrink and swll phnomnon that occurs whn drum prssur changs (DiDomnico t al., 1983; Rs, 1997). Th controllr action tnds to ract ngativly du to th mislading shrink and swll ffct. Th control schm that is adoptd has to countr th ffct of th shrink and swll phnomnon in ordr to optimiz th cost of oprating th powr plant.

Stand. Sci. Rs. Essays Lawrnc 020 METHODOLOGY Mod of oprations n ordr to work around th non-linar bhaviour of th plant along its oprating rgion (that is, 0-100% load), it will b subdividd into thr diffrnt rgions-dtrmind by thr diffrnt load lvls. Th valus shown in Tabl 1 ar dtrmind by th powr plant nginrs in such a way that thy rprsnt th most common oprating mods of th boilr-drum. Th PD-basd control mods of th boilr-drum. Th PD-basd control schm will only work if w trat th non-linar systm as bing linar (locally) at ach rgion. Tabl 1. Mod of Oprations Paramtr Low load load Mdium load High load Prssur, p (mpa) 8.70 9.35 10.00 Powr Lvl (MW) 24.00 80.00 136.00 Opn loop rsponss To illustrat th dynamic bhaviour of th modl, w will simulat th rsponss-to-stp changs in th inputs. Th modl was simulatd using th following paramtrs, basd on Orsund Powr Station in Swdn: m t =,000 kg; m r = 20,000 kg; A d = 20 m 2 ; V d = 40 m 3 ; v r = 37 m 3 ; v dc = 11 m 3 ; V sd0 = 8 m 3 ; C p = 650; C fw = 4.18; k = 25;b =0.3. Th stam tabls usd in th simulation ar approximatd using quadratic approximation. Th quadratic approximation implmntation is givn in Fawnizu (1). As dscribd arlir, th drum variabls ar non-linar. nitial incras of drum watr lvl du to shrink and swll phnomnon. Th comparison of rsponss to 10 kg/stam flow rat chang at diffrnt oprating conditions. Th swlling ffct on th watr lvl is largst at low load, as illustratd by th gratst ovrshoot. This is mainly causd by th gratr variation of stam contribution at low load than at high load. Th non-minimum phas charactristics of th drum watr lvl prov to b non-trivial. Any control schm that w us to control th watr lvl will b affctd by th non-minimum phas charactristics. n addition, th varying snsitivity of th paramtrs at diffrnt oprating condition will impos xtra difficulty to th control dsign. Singl lmnt control Singl lmnt fdwatr control uss only th drum lvl procss variabl as fdback. Th masurd drum lvl is compard to th drum lvl st point. Whn th diffrnc dviats from zro (that is, non-zro rror), th fd watr control valv will b adjustd by th proportional-plus-intgral (P) controllr to compnsat for th rror. Th intgral componnt of th controllr rgulats th drum lvl rror to zro. This control mthod nsurs that th watr lvl stays as clos as possibl to st point valu. This schm prforms satisfactorily undr constant or small changs in load (stam flow) and saturation prssur in th drum. f th stam flow incrass fastr than th hat input, th drum prssur drops quickly and causs th saturation condition to chang rapidly from on stat to anothr (Astrom t al., 0). This phnomnon is known as swll ffct bcaus of th swlling of stam bubbls undrnath watr surfac. This would caus th watr lvl to ris initially. On th contrary, a suddn dcras in stam flow would caus th stam bubbls undrnath watr surfac. This would caus th watr lvl to ris initially. On th contrary, a suddn dcras in stam flow would caus th stam bubbls to shrink and th watr lvl to drop. f th transint is not too svr, th lvl will vntually rturn to th stpoint. n ths circumstancs, mor complicatd control is ncssary. Thr-lmnt control Thr-lmnt control schm (Figur 1) uss fd watr flow masurmnts and stam flow masurmnts as inputs to th controllr in addition to watr lvl fdback signals. This improvd control schm adds prdictability by anticipating chang in load by using stam flow as fd forward and fd watr flow rat as fdback rgulation. Figur 2 shows th comparison btwn singl lmnt control and thr lmnt control rsults on watr lvl and stam prssur. Th improvmnt shown by th improvd control mthod is vry significant. Thr is hardly any oscillation shown by th solid lin, xcpt during th shrink and swll (non-minimum phas) ffct.

Lawrnc 021 hat input Lvl + - P + + - P Fdwatr Valu Boilr Drum Modl Stam flow Stam prssur stpoint Watr lvl Figur 1. Thr-lmnt controllr structur. Hat control loop is not shown - 0.05 Drum Watr Lvl, 1 0 - - 900 Stam Prssur, P (MPa) - 9.4 9.35 9.3-9.25-9.2-9.15 900 1000 Tim (s) Figur 2. Singl-lmnt controllr (dottd) and 3-lmnt controllr(solid). Stp rsponss to 10 kg/s chang of stam flow rat on boilr-drum modl at mdium load 100 - nput Fdwatr, wfw (kg/s 95 - - 90 85-80 75 900 150 - nput powr, Q (MW) 145 - - 140 135-130 - 125 900 Tim (s) Figur 3. Singl-lmnt controllr (dottd) and 3-lmnt controllr (solid). Stp rsponss to 10 kg/s chang of stam flow rat on boilr-drum modl with dlayd input at mdium Load Figur 3 shows th fdwatr and powr input corrsponding to th watr lvl and stam prssur rsponss in Figur 2. Th thr-lmnt controllr input fdwatr liminats th initial rduction shown by th dottd lin. This ffct is

Lawrnc 022 du to th fd forward (prdiction) stam flow signal, which incrass controllr flow vn whn watr lvl is incrasing. Th controllr knows that watr lvl will vntually dcras. Thrfor, it did not ract to th mislading changs in watr lvl. Th oscillation in input fd watr right aftr th stp input is du to th multi-cascadd structur of th controllr fdback and fd forward loops. DSCUSSON Shrink and swll phnomnon Figur 4 shows th simplifid diagram of th boilr-drum and sown comr-risr circulation loop. Whn powr dmand incrass, stam flow rat is rapidly incrasd, causing th stam prssur to drop momntarily. This drop of prssur causs th air bubbls to incras in siz and th watr lvl to incras. Th phnomnon is trmd swll ffct. Th principl of mass balanc, howvr, dictats that th incras of stam flow rat laving th drum will caus rduction of th total mass insid th drum. Thus, by kping th fdwatr input constant, th mass of watr insid th drum will vntually b dcrasd, causing th watr lvl to drop. This is shown in th opn-loop rspons in Figur 5. On th othr hand, if stam dmand is rducd, stam bubbls shrink initially and th watr lvl will vntually incras du to th incrasing mass of watr and stam in th drum. Th combind phnomnon is trmd shrink and swll. FEEDWATER BOLER DRUM STEAM FLOW RSER HEAT DOWNCOMER Figur 4. Drum Downcomr-Risr circulation loop Boilrs-drum dynamics Drum boilr modl (Astrom t al., 0) is formd using th basic thrmodynamics mass and nrgy balanc quations. Blow ar th stat quations of th boilr-drum dynamics. Th four stats ar: drum prssur p, total watr volum v wt, stam quality at th risr outlt sd and volum of stam undr th liquid lvl in th drum v sd. Th quations hav bn drivd using mass and nrgy balancs togthr. Th rsulting physical quations hav thn bn manipulatd into th following form: 11 dv wt + 12 dp = q f q s dt dt 21 dv wt + 22 dp = Q + q f h f q s h s dt dt 32 dp + 33 d r = Q - r h r q dc dt dt 42 dp + 43 d r + 44 dv sd = p s (V 0 sd V sd ) + h f h w q s (1) dt dt dt T d h c Whr th cofficints ij ar givn by: 11 = p w - p s

Lawrnc 023 12 = V wt dp w + V st dp s dp dp 21 = p w h w - p s h s 22 V wt h w dp w + p w dp w + V st h s dp s + p s dp s - V t + m t C p dt s dp dp dp dp dp 32 p w dp w + r h c dp w ( 1- v ) V r + (1- r ) h c dp s + p s dh s v V r dp dp dp dp + (p s + (p w -p s ) r ) h c v r d v - v r + m r C p dt s dp dp 33 = ((1- r ) p s + r p w ) h c v r d v dp 42 V sd dp s + 1 p s V sd dh s + p w V wd dh w v sd + m d C p dt s dp h c dp dp dp + r (1+β) V r v dp s + (1- v ) dp w + (p s p w ) d v dp dp dp dp 43 = r (1+β) (p s -p w ) V r d v d r 44 = p s (2) Th paramtrs notations ar, xampl, pm is spcific dnsity of mtal. Paramtrs V: volum; p: spcific dnsity; u: Spcific intrnal nrgy; h: Spcific Enthalpy; t: Tmpratur; q: mass flow rat; q: hat input; r : stam quality; v : stam volum ratio; C p : spcific hat of mtal; M t : total mass for mtal tub and drum. Subscripts s: Stam; w: watr; f: fdwatr; m: mtal; d: drum; dc: Downcomr; t: total. Modlling with simulink This modl is implmntd in Matlab in th form of an S function, in ordr to b usd in Simulink to build th non-linar systm in block diagrams. Th working S function is givn in Fawnizu t al. (1). Th S-function is thn maskd to bcom a 4 input, 2 output MMO systms shown in Figur 5. Hat input Fdwatr flow Q qf Drum watr lvl Fdwatr tmpratur Stam flow tf qs Stam prssur Boilr-Drum Figur 5. Maskd MMO systm of a drum-boilr Using S function

Lawrnc 024 CONCLUSON Th simpl, pur fdback singl-lmnt PD controllr prforms satisfactorily to rgulat watr lvl rror to zro whn oprating lvl stays constant. n today s powr industry, comptition has drivn powr plant oprators to maximiz profit by optimizing powr gnratd-powr dmandd ratio. n ordr to follow th changing dmand, powr plant boilr has to chang its stam production rapidly; thrfor th load-following problm. Thr-lmnt control schm is a good altrnativ to th simplr singl lmnt control structur. Th prformanc shown by Figur 2 indicats how thr-lmnt controllr is a much bttr altrnativ. Rfrncs Astrom K J, Bll R D (0). Drum Boilr Dynamics, Automatica vol. 36 Pp 363-378. Di Domnico, Ptr N (1983). Practical Application of Fdwatr Controls for a Utility Typ Drum Boilr. Amrican Control Confrnc, San Francisco, Pp 22-24. Fawnizu A H, Rs N W (1). Drum Watr Lvl Control: A study by simulation. MEngSc thsis submittd to UNSW, Australia. Rs N W (1997).Advancd Powr Plant Control For Larg Load Changs and Disturbancs, FAC/CGRE Symposium on Control of Powr Systms and Powr Plant, Bijing, China. Pp18-21. Waddington J, Mapls G C (1987). Th Control of Larg Coal-And Oil-Fird Gnrating Units. Powr Enginring Journal. Pp 102-135