COMBUSTION BOILER CONTROL METHOD, COMBUSTION BOILER, AND BOILER COMPUTATION SYSTEM

Abstract

A combustion boiler control method includes steps of (a) monitoring the current load of a combustion boiler, (b) finding a numerical value for a current computational maximum boiler momentary load for which at least one flue gas factor computed using currently monitored process data with a numerical model of the boiler fulfills an acceptance condition, and selecting the numerical value as the current computational maximum boiler momentary load, (c) indicating the current computational maximum boiler momentary load to an operator and/or, if the current load is (c1) less than the current computational maximum boiler momentary load, (c1i) indicating to the operator that the boiler load may be increased, and/or (c1ii) automatically increasing the boiler load, and/or (c2) greater than the current computational maximum boiler momentary load, (c2i) indicating to the operator that the boiler load exceeds the current computational maximum boiler momentary load, and/or (c2ii) automatically reducing the boiler load.

Claims

1.-27. (canceled)

28. A combustion boiler control method comprising the steps of: (a) monitoring a current load (Q.sub.h) of a combustion boiler; (b) finding such a numerical value (Q.sub.h, candidate) for a current computational maximum boiler momentary load (Q.sub.h, max) for which at least one flue gas factor (df.sub.i) computed using currently monitored process data with a numerical model of the boiler fulfills an acceptance condition, and selecting the numerical value (Q.sub.h, candidate) as the current computational maximum boiler momentary load (Q.sub.h,max); (c) indicating the current computational maximum boiler momentary load (Q.sub.h,max) to a boiler operator and/or, if the current load (Q.sub.h) is (c1) less than the current computational maximum boiler momentary load (Q.sub.h,max): (c1i) indicating to the boiler operator that the boiler load (Q.sub.h) may be increased; and/or (c1ii) automatically increasing the boiler load (Q.sub.h); and/or (c2) greater than the current computational maximum boiler momentary load (Q.sub.h,max): (c2i) indicating to the boiler operator that the boiler load (Q.sub.h) exceeds the current computational maximum boiler momentary load; and/or (c2ii) automatically reducing the boiler load (Q.sub.h).

29. The method according to claim 28, wherein: (i) the currently monitored process data of the boiler includes: (ia) current flue gas exit temperature (T.sub.flue gas,exit,current) in a flue gas flow channel; and (ib) heat duty (Q.sub.fluid,i) for each heat transfer surface in the flue gas flow channel, and further wherein: (ii) monitored process data from both (ia) and (ib) is used in computation of the flue gas factor and when finding the numerical value (Q.sub.h, candidate) for the current computational maximum boiler momentary load (Q.sub.h,max).

30. The method according to claim 28, wherein the finding is performed such that, if the at least one flue gas factor (df.sub.i) computed using currently monitored process data with a numerical model of the boiler that fulfills an acceptance condition for the numerical value (Q.sub.h, candidate) for the current computational maximum boiler momentary load (Q.sub.h,max) fails to fulfill an acceptance condition, a next numerical value (Q.sub.h, candidate) is automatically selected.

31. The method according to claim 30, wherein the next numerical value (Q.sub.h, candidate) is selected iteratively.

32. The method according to claim 28, wherein the finding is carried out by performing the computational steps of: (I) computing an estimate for boiler flue gas exit temperature (T.sub.boiler, exit) that results in a computational boiler model when the thermal load of the boiler corresponds to the numerical value (Q.sub.h, candidate); (II) computing flue gas mass flow (q.sub.m,fluegas); (III) computing a heat duty (Q.sub.fluid, i, candidate) for each heat transfer surface in the flue gas flow channel with its current heat duty (Q.sub.fluid, i, current) that is corrected by using a numerical boiler model (Q.sub.fluid, i, candidate=Q.sub.fluid,i,current+S a.sub.j,I(Q.sub.steam,max).sup.jS a.sub.j,i(Q.sub.steam,current).sup.j); (IV) using the computed heat duties (Q.sub.fluid, i, candidate) for each heat transfer surface in the flue gas flow channel to compute flue gas temperatures at each heat transfer surface (T.sub.fluegas,in,i, T.sub.fluegas,out,i; i=1, . . . , k) in the flue gas flow channel in the upstream direction of flue gas flow, starting from the heat transfer surface 21.sub.k that is closest to the flue gas exit using the estimate for the boiler flue gas exit temperature (T.sub.fluegas,out,k=T.sub.FG, exit); and (V) computing a flue gas factor (df.sub.i, i=1, . . . , k) for each heat transfer surface in the flue gas flow channel.

33. The method according to claim 32, wherein the flue gas factor includes or is: $d{f}_i={k_i\left(\frac{q_{m,\mathrm{fluegas}}}{{}_{\mathrm{fluegas},i}*A_{cross,i}}\right)}^n$ where k.sub.i is a non-zero parameter that may be chosen combustion-boiler specifically, q.sub.m,fluegas is flue gas mass flow, and n is a model parameter that may be chosen combustion-boiler specifically, .sub.fluegas,i is flue gas density at an i.sup.th heat transfer surface and A is a cross section of flue gas channel at the i.sup.th heat transfer surface.

34. The method according to claim 33, wherein n is selected to include at least one of the following: (i) in the range of 0.9 to 1.1, for using computed flue gas velocity; (ii) in the range of 2.9 to 3.5, for using computed flue gas caused erosion; or (iii) in the range of 1.8 to 2.2, for using pressure loss.

35. The method according to claim 34, wherein the value for n is changed over time.

36. The method according to claim 34, wherein the value for n is determined from a group of boilers comprising at least two separate boilers using operational data monitored for each of the boilers.

37. The method according to claim 32, wherein, in the computation in step (I), the flue gas exit temperature is substantially estimated by an equation: $T_{bo\mathrm{iler},\mathrm{exit}}=a_0+S{a_j\left(Q_{h,candidate}\right)}^j$ or its first, second, third, or higher degree approximation, and wherein the respective coefficients (a.sub.0, a.sub.1, a.sub.2, . . . ) have been obtained beforehand by fitting, after measuring flue gas exit temperature (T.sub.FG, exit) values for a number of discrete boiler load (Q.sub.steam) values.

38. The method according to claim 32, wherein, in step (II), computation of flue gas mass flow utilizes mass flow (q.sub.m,fluegas,m) of flue gas components, wherein the components include CO.sub.2, H.sub.2O, N.sub.2, SO.sub.2, and O.sub.2.

39. The method according to claim 32, wherein, in step (II), the computation of flue gas mass flow includes fuel parameters.

40. The method according to claim 28, wherein the step (b) is performed remotely from the combustion boiler.

41. The method according to claim 28, wherein the step (b) is performed locally at the combustion boiler.

42. The method according to claim 28, wherein any of the currently monitored process data and/or current load is obtained from real-time measurements, treated by filtering, treated by averaging, computing trends, or any combination of these.

43. The method according to claim 28, wherein the acceptance condition includes a hysteresis condition, requiring a predefined minimum change before changing the current computational maximum boiler momentary load (Q.sub.h,max).

44. The method according to claim 28, wherein the acceptance condition includes comparing the computed at least one flue gas factor (df.sub.i) against a respective design value, and wherein, in the method, the numerical value (Q.sub.h, candidate) is rejected if the design value is exceeded.

45. The method according to claim 28, wherein the combustion boiler is a circulating fluidized bed (CFB) or a bubbling fluidized bed (BFB) boiler, and the step (b) is carried out for the combustion boiler heat transfer surfaces, between the furnace and the stack, optionally, including the furnace.

46. A combustion boiler comprising: a furnace and associated passes defining a flue gas flow path and having a number of heat transfer surfaces that are located in the flue flow path; measurement instrumentation to monitor current load (Q.sub.h) of the combustion boiler; further, measurement instrumentation, such as sensors, to currently monitor process data; and a control system configured to carry out the combustion boiler control method according to claim 28.

47. The combustion boiler according to claim 46, wherein the control system comprises an edge server that is configured to process real-time measurement results for currently monitored process data and/or current load, namely, by filtering, averaging, and/or computing trends.

48. The combustion boiler according to claim 46, wherein the control system is configured to carry out the method step (b) to determine the current computational maximum boiler momentary load (Q.sub.h,max) locally.

49. The combustion boiler according to claim 46, wherein the control system is configured to send data to a remote computing system that is configured to carry out the method step (b) and to return the current computational maximum boiler momentary load (Q.sub.h,max) to the control system.

50. The combustion boiler according to claim 49, further comprising an edge server that is configured to reduce an amount of measurement data that is passed to the remote computing system.

51. A combustion boiler computation system comprising: a group of combustion boilers comprising at least two separate combustion boilers according to claim 46, each boiler comprising a boiler control system comprising an edge server system that is configured to process the real-time measurement results for currently monitored process data and/or current load, namely, by filtering, averaging, and/or computing trends, and to send the processed real-time measurement results to a remote computing system; a remote computing system configured to receive data processed from real-time measurement results and to compute data using a numerical boiler model for each of the combustion boilers, and to return computation results for each of the combustion boilers, wherein the control system is configured to adapt its function based on the computation results, and the computing system is configured to find such a numerical value (Q.sub.h, candidate) for a current computational maximum boiler momentary load (Q.sub.h,max) for which at least one flue gas factor (df.sub.i) computed using currently monitored process data with a numerical model of the boiler that fulfills an acceptance condition, and selecting the numerical value (Q.sub.h, candidate) as the current computational maximum boiler momentary load (Q.sub.h,max).

52. The combustion boiler computation system according to claim 51, wherein the combustion boiler computation system is configured to adapt or to calibrate a numerical model for a combustion boiler using processed measurement data for the combustion boiler.

53. The boiler computation system according to claim 51, wherein the boiler computation system is configured to adapt or to calibrate a numerical model for a combustion boiler using processed measurement data also collected from other combustion boilers.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0084] The combustion boiler and its control method are explained in more detail below in the context of the embodiments shown in the appended drawings in FIG. 1 to FIG. 8, of which:

[0085] FIG. 1 illustrates a CFB boiler;

[0086] FIG. 2 illustrates a BFB boiler;

[0087] FIG. 3 illustrates the flow of measurement data from sensors;

[0088] FIG. 4 is a flow diagram illustrating a first method for finding the current computational maximum boiler momentary load Q.sub.h, max;

[0089] FIG. 5 is a flow diagram illustrating a second method for finding the current computational maximum boiler momentary load Q.sub.h, max;

[0090] FIG. 6 illustrates how the current computational maximum boiler momentary load Q.sub.h, max can be presented to the boiler operator;

[0091] FIG. 7 shows boiler momentary load Q.sub.h and computed current computational maximum boiler momentary load Q.sub.h, max, as well as the effect of using the method according to the invention during a test period;

[0092] FIG. 8 shows a closer look at the data of FIG. 7, showing boiler momentary load Q.sub.h computed current computational maximum boiler momentary load Q.sub.h, max where the effect of using the method according to the invention during the 10 day test period is better visible.

[0093] The same reference numerals refer to same technical features in all FIG.

DETAILED DESCRIPTION

[0094] FIG. 1 shows a combustion boiler 10 that is a CFB boiler and comprises a furnace 12 that has tube walls 13 connected to water-steam circuit of the combustion boiler 10. Water is fed from water tank (not shown) to an economizer and, from the economizer, via a steam drum to evaporative heat transfer surfaces such as the tube walls 13, and then guided via the steam drum to superheaters and then to a turbine. A flue gas channel may be provided with the economizer and/or superheater/s.

[0095] Fluidization gas (such as, air and/or oxygen-containing gas) is fed from fluidization gas supply 153 to below the grate (the grate not shown in FIG. 1) via a windbox (not shown), wherefrom the primary fluidization air enters into the furnace through nozzles (not shown) (to fluidize the bed), and secondary fluidization gas feed 152 (to feed oxygen containing gas to control combustion). The effect is that the bed materials will be fluidized and also oxygen required for the combustion is provided into the furnace 12. Further, fuel is fed into the furnace 12 via the fuel feed 22. The combustion can be adjusted by controlling the fuel feed 22 (such as, by reducing or increasing fuel feed), and by controlling the fluidization gas feed (such as, by reducing or increasing the amount of oxygen supply into the furnace 12). Fuel can be fed together with additives, in particular, with such additives that act as alkali sorbents, such as CaCO.sub.3 and/or clay for example. In addition or alternatively, NOx reduction agents, such as ammonium or urea can be fed into the combustion zone of the furnace 12, or above the combustion zone of the furnace 12.

[0096] Bed material is also fed into the furnace, which bed material may comprise sand, limestone, and/or clay, that, in particular, may comprise kaolin. One effect of the bed and, generally, of the combustion, is that, in the water-steam circuit, water and steam is heated in the tube walls 13 and water is converted to steam.

[0097] Ash may fall to the bottom of the furnace 12 and be removed via an ash chute (omitted from FIG. 1 for the sake of clarity) and part of the ash, so-called fly ash, will be carried along the flue gas.

[0098] Combustion products, such as flue gas, unburnt fuel and bed material proceed from the furnace 12 to a particle separator 17 that may comprise a vortex finder 103. The particle separator 17 separates flue gases from solids. Especially, in larger combustion boilers 10, there may be more than one (two, three, . . . ) separators 17 preferably arranged in parallel to each other.

[0099] Solids separated by the separator 17 pass through a loop seal 160 that preferably is located at the bottom of the separator 17. Then, the solids pass to fluidized bed heat exchanger (FBHE) 100 that is also a heat transfer surface so that the FBHE 100 collects heat from the solids to further heat the steam in the water-steam circuit. The chamber in which the FBHE 100 is located may be fluidized and the FBHE 100 itself comprises heat transfer tubes or other kinds of heat transfer surfaces. FBHE 100 may be arranged as a reheater or as a superheater. From the FBHE outlet 101, steam is passed into a high-pressure turbine (if the FBHE 100 is a superheater) or medium-pressure turbine (if the FBHE 100 is a reheater). For the sake of clarity, the turbines are not illustrated in FIG. 1. The solids may be returned from the FBHE 100 via a return channel 102 into the furnace 12. Especially, in larger combustion boilers 10, there may be more than one (two, three, . . . ) loop seals 160 and FBHE 100, and return channels 102, preferably, arranged in parallel to each other, such that for each separator 17, there will be respective loop seal 160, FBHE 100 and return channel 102. In practice, some of the FBHE 100 may be arranged as superheaters while some others may be arranged as reheaters.

[0100] The flue gases are passed from the separator 17 to horizontal pass 15 and, from there, further to backpass 16 (that, preferably, may be a vertical pass) and from there via flue gas conduit 18 to stack 19.

[0101] The backpass 16 comprises a number of heat transfer surfaces 21.sub.i (where i=1, 2, 3, . . . , k, where k is the number of heat transfer surfaces). In FIG. 1, heat transfer surfaces 21.sub.1, 21.sub.2, 21.sub.3, . . . , 21.sub.k-1, 21.sub.k are illustrated. Heat transfer surface 21.sub.k depicts air preheater. Heat transfer surfaces 21.sub.k-1, 21.sub.2 depict superheaters and heat transfer surfaces 21.sub.1, 21.sub.3 depict reheaters. The actual number of different heat transfer surfaces in each of these components, for example, may be selected for each combustion boiler differently according to actual needs. And, there may be further components as well, comprising a heat transfer surface 21.

[0102] Flue gas exiting the last heat transfer surface 21.sub.k will be at the flue gas exit temperature T.sub.FG, exit. This temperature is measured with temperature sensor 20.sub.k.

[0103] According to one aspect, the temperatures before and after each heat transfer surface 21.sub.i (T.sub.FG,in,i, T.sub.FG,in,i+1, respectively) can be measured with respective temperature sensors 20.sub.i (where i=1, 2, 3, . . . , k1, k).

[0104] According to another aspect, and, preferably, these temperatures, however, do not necessarily need to be measured. It will suffice to know the flue gas exit temperature T.sub.FG, exit. The temperatures before and after each preceding heat transfer surface 21.sub.i (T.sub.FG,in,i, T.sub.FG,in,i+1) can be obtained numerically. This will be explained further below.

[0105] A combustion boiler 10 is equipped with a plurality of sensors and computer units. Actually, one middle-size (100-150 MW.sub.th) combustion boiler 10 may produce one hundred million measurement results/day, which needs 25 GB of storage space. FIGS. 1, 2, and 3 illustrate some of the sensors and computer units. Examples of sensors are combustion gas (usually combustion air) volume flow sensors 30 (for measuring primary and secondary fluidizing gas feeds), fuel feed sensors 650 and temperature sensors 20.sub.i (i=1, 2, . . . , k), temperature sensor in FBHE and pressure sensor 116 in the return channel 102 (both only in a CFB boiler), and sensors 40 in the furnace 12.

[0106] Process data may be collected from the sensors by distributed control system (DCS) 201. The data collection may most conveniently be arranged over a field bus 290, for example. DCS 201 may have a display/monitor 202 for displaying operational status information to the operator. An EDGE server 203 may process measurement data from the obtained from sensors, such as, filter and smooth it. There may be a local storage 204 for storing data.

[0107] The DCS 201, display/monitor 202, EDGE server 203, local storage 204 may be in combustion boiler network 280 (local storage 204 preferably directly connected to the EDGE server). The combustion boiler network 280 is preferably separate from the field bus 290 that is used to communicate measurement results from the sensors to the DCS 201 and/or the EDGE server 203. Between the DCS 201 and EDGE server 203 there may be an open platform communications server 210 (cf. FIG. 3) to make the systems better interoperable.

[0108] Combustion boiler network 280 may be in connection with the internet 200, preferably, via a gateway 290. In this situation, measurement results may be transferred from the combustion boiler network 280 to a cloud service, such as process intelligence system 205 located in a computation cloud 206. The applicant currently operates a cloud service running an analysis platform. The cloud service may be operated on a virtualized server environment, such as on Microsoft Azure, which is a virtualized, easily scalable environment for distributed computing and cloud storage for data. Other cloud computing services may be suitable for running the analysis platform too. Further, instead of a cloud computing service, or in addition thereto, a local or a remote server can be used for running the analysis platform.

[0109] FIG. 2 illustrates a combustion boiler 10 that is a BFB boiler. A BFB boiler differs from a CFB boiler in that the fluidized bed is not a circulating bed but a bubbling bed. Thus, there is no need for the separator 17, loop seal 160, FBHE 100, and return channel 102.

[0110] There is normally at least one superheater 14 located in the furnace 12, preferably, on top of the furnace 12. Superheater 14 inlet 141 is preferably the steam drum or from another superheater and the outlet 142 is to high pressure turbine.

[0111] FIG. 4 illustrates the combustion boiler control method: [0112] (a) the current load Q.sub.h of combustion boiler 10 is monitored in step K1 (in the method illustrated in FIG. 4, also flue gas exit temperature T.sub.FG, exit is monitored and heat duty Q.sub.fluid,i for each heat transfer surface 21i in the flue gas flow channel (vertical pass 16); [0113] (b) a numerical value Q.sub.h, candidate is selected (step K3), after which heat duties at heat transfer surfaces 21.sub.i are computed and flue gas temperatures in relation to Q.sub.h, candidate. The numerical value Q.sub.h, candidate is then used to compute (step K7) at least one flue gas factor df.sub.i using currently monitored process data with a numerical model of the boiler fulfills an acceptance condition (which is tested in step K9), and selecting the numerical value; and Q.sub.h, candidate as the current computational maximum boiler momentary load Q.sub.h, max (step K11); [0114] (c) the current computational maximum boiler momentary load Q.sub.h, max is indicated to the operator (such as, by displaying on the monitor/screen 202) and/or, if the current load Q.sub.h is (c1) less than the computational boiler maximum momentary load Q.sub.h,max: (c1i) indicating the boiler operator that the boiler load Q.sub.h may be increased, and/or (c1ii) automatically increasing the boiler load Q.sub.h, and/or (c2) greater than the computational boiler maximum momentary load Q.sub.h,max: c2i) indicating the boiler operator that the boiler load Q.sub.h exceeds the boiler maximum momentary load, and/or (c2ii) automatically reducing the boiler load Q.sub.h.

[0115] The step (b) is preferably carried out for the combustion boiler 10 heat transfer surfaces 21.sub.i between furnace 12 and stack 19.

[0116] In the method, the currently monitored process data of the boiler may include (a) current flue gas exit temperature T.sub.FG,exit in a flue gas flow channel and b) heat duty Q.sub.fluid,i for each heat transfer surface 21.sub.i in the flue gas flow channel (back pass 16).

[0117] Further, in the method, monitored process data from both (a) and (b) may be used in computation of the flue gas factor df.sub.i and when finding the numerical value Q.sub.h, candidate for the current computational maximum boiler momentary load Q.sub.h,max.

[0118] The finding is performed such that, if the at least one flue gas factor df.sub.i computed using currently monitored process data with a numerical model of the boiler that fails to fulfill an acceptance condition, a next numerical value Q.sub.h, candidate is automatically selected. The automatic selection is preferably done iteratively.

[0119] As a specific example, the finding may be carried out with performing the computational steps of: [0120] I: computing an estimate for boiler flue gas exit temperature T.sub.boiler, exit that results in a computational boiler model when the thermal load of the boiler corresponds to the numerical value Q.sub.h, candidate; [0121] II: computing flue gas mass flow q.sub.m,fluegas; [0122] III: computing a heat duty Q.sub.fluid, i, candidate for each heat transfer surface 21.sub.i in the flue gas flow channel (back pass 16) with its current heat duty Q.sub.fluid, i, current that is corrected by using a numerical boiler model Q.sub.fluid, i, candidate=Q.sub.fluid,i,current+ a.sub.j,i (Q.sub.h,candidate).sup.j par.sub.j,i (Q.sub.h,current).sup.j [0123] IV: using the computed heat duties Q.sub.fluid, i, candidate for each heat transfer surface 21.sub.i in the flue gas flow channel (back pass 16) to compute flue gas temperatures at each heat transfer surface (T.sub.fluegas,in,i, T.sub.fluegas,out,i; i=1, . . . , k) in the flue gas flow channel (back pass 16) in the upstream direction of flue gas flow, starting from the heat transfer surface 21.sub.k that is closest to the flue gas exit i.e. using the estimate for the boiler flue gas exit temperature T.sub.fluegas,out,m=T.sub.FG, exit; and V: computing a flue gas factor df.sub.i, i=1, . . . , k for each heat transfer surface 21.sub.i in the flue gas flow channel (back pass 16).

[0124] The fitting of the parameters (par.sub.j,i) can be done manually by a human or automatically by a computer utilizing historical data. Automatic update of the parameters may be done, e.g., once per month. AI and neural network based algorithms can be utilized in automatic update.

[0125] Step (II) may include computing flue gas mass flow q.sub.m,fluegas,m for selected flue gas components.

[0126] The flue gas temperatures at each heat transfer surface can be computed, for instance,

[00008]$T_{\mathrm{fluegas},in,i}=T_{\mathrm{fluegas},\mathrm{out},i}+\frac{Q_{\mathrm{fluid},i}}{q_{m,\mathrm{fluegas}}*c_p}$

wherein T.sub.fluegas,in,i is the flue gas temperature at the inlet of ith heat transfer surface, c.sub.p is specific heat capacity, and T.sub.fluegas,out,i is the flue gas temperature at the outlet of ith heat transfer surface. The flue gas temperatures could be determined with artificial intelligence tools. The flue gas temperatures could be determined with neural network.

[0127] Preferably, the flue gas factor df.sub.i includes or is:

[00009]$d{f}_i={k_i\left(q_{m,\mathrm{fluegas}}/\left({}_{\mathrm{fluegas},I}{A}_{\mathrm{cross},i}\right)\right)}^n$

where k.sub.i is a predetermined non-zero parameter that may be chosen combustion-boiler specifically, preferably, positive (non-zero) number, q.sub.m,fluegas is a flue gas mass flow, n is a positive number (which may be selected as a natural number, rational number, real number, or even as complex number), .sub.fluegas,i is flue gas density obtainable from flue gas temperature T.sub.FG, in, i at i.sup.th heat transfer surface 21.sub.i and A is a cross section of flue gas channel at i.sup.th heat transfer surface 21.sub.i.

[0128] Advantageously, n may be selected to include at least one of the following: [0129] (i) in the range of 0.9 to 1.1, preferably, equivalent or about 1.0, for using computed flue gas velocity; [0130] (ii) in the range of 2.9 to 3.5, preferably, between 3.2 and 3.35, for using computed flue gas caused erosion; or [0131] (iii) in the range of 1.8 to 2.2, preferably equivalent or about 2.0, for using pressure loss.

[0132] The value for n may be changed over time. In particular, the value for n may be determined from a group of combustion boilers, the group comprising at least two separate combustion boilers 10, such that using operational data monitored for each of the combustion boilers 10 is used in the determination.

[0133] In the computation in step (I), the computational value for flue gas exit temperature T.sub.FG, exit under any chosen numerical value Q.sub.h, candidate for boiler load can be estimated by equation:

[00010]$T_{\mathrm{FG},\mathrm{exit}}={}_0+.Math.{{}_j\left(Q_{h,\mathrm{candidate}}\right)}^i$

or, preferably, its first, second, third, or higher degree approximation. The coefficients a.sub.0, a.sub.1, a.sub.2, . . . have been obtained beforehand by fitting after measuring flue gas exit temperature T.sub.FG, exit values for a number of discrete boiler load Q.sub.steam values.

[0134] In step (II), the computation of the components q.sub.m,fluegas,m preferably includes at least some, most preferably, all of the following: m=CO.sub.2, H.sub.2O, N.sub.2, SO.sub.2, O.sub.2 so as to determine flue gas mass flow. In other words, in step (IV) of the computation, as q.sub.m,fluegas,m values some or all of q.sub.m,fluegas,CO2, q.sub.m,fluegas,H2O, q.sub.m,fluegas,N2, q.sub.m,fluegas,SO2, q.sub.m,fluegas,O2 may be used. They are preferably measured in flue gas conduit 18 or in flute 19, for which reason suitable sensors are installed in the flue gas passage. In step (II), the component values may further include fuel parameters.

[0135] Flue gas mass flow may be based on computation of sums of flue gas component mass flows q.sub.m,fluegas,m that are calculated based on fuel analysis (proximate and ultimate analysis of fuel), combustion air flow and/or recirculation gas flow according to boiler mass and energy balance calculation.

[0136] Preferably, the flue gas mass flow may be computed:

q.sub.m,fluegas=q.sub.m,fluegas,i [0137] i.e., for example, the sums of the following flue gas mass flow components CO.sub.2, H.sub.2O, N.sub.2, SO.sub.2 and O.sub.2:

[00011]$q_{m,\mathrm{fluegas},CO2}=x_{C,\mathrm{fuel}}*\frac{M_{\mathrm{CO}2}}{M_C}*q_{m,\mathrm{fuel}}$$q_{m,fluegas,H2O}=0.5*x_{H,\mathrm{fuel}}*\frac{M_{H2O}}{M_H}*q_{m,\mathrm{fuel}}+x_{H2O,\mathrm{fuel}}*q_{m,\mathrm{fuel}}+x_{moist,\mathrm{air}}*q_{m,\mathrm{air}}$$q_{m,\mathrm{fluegas},N2}=0.5*x_{N,fuel}*q_{m,\mathrm{fuel}}+x_{N2,\mathrm{air}}*q_{m,\mathrm{air}}$$q_{m,\mathrm{fluegas},\mathrm{SO}2}=x_{S,fuel}*\frac{M_{\mathrm{SO}2}}{M_S}*q_{m,fuel}$$q_{m,flue\mathrm{gas},O2}=x_{O2,\mathrm{air}}*q_{m,\mathrm{air}}-q_{m,fluegas,CO2}*\frac{M_{O2}}{M_{CO2}}-0.25*x_{H,fuel}*\frac{M_{H2O}}{M_H}*q_{m,fuel}-q_{m,flue\mathrm{gas},SO2}*\frac{M_{O2}}{M_{\mathrm{SO}2}}$

where, for instance, x.sub.C,fuel represents carbon in fuel i.e. first subscript denotes component and second subscript is either fuel or combustion air referred, q.sub.m,fuel is a fuel flow, q.sub.m,air is combustion air flow and M.sub.x denotes molar mass. Advantageously, fuel properties as utilized in flue gas mass flow components and combustion air properties. Fuel moisture may be measured or calculated.

[0138] The step (b) may be performed remotely to the combustion boiler, such as, in the process intelligence system 205. Alternatively, the step (b) may be performed locally at the combustion boiler, preferably, at the EDGE server 203.

[0139] Any of the currently monitored process data and/or current load may be obtained from real-time measurements, treated by filtering, treated by averaging, computing trends, or any combination of these.

[0140] The acceptance condition may include a hysteresis condition, requiring a predefined minimum change before changing the current computational maximum boiler momentary load Q.sub.h,max.

[0141] The acceptance condition preferably includes comparing the computed at least one flue gas factor df.sub.i against a respective maximum value df.sub.max,i. The maximum value df.sub.max,i is a preset value and preferably boiler specific. The numerical value Q.sub.h, candidate is rejected if the maximum value df.sub.max,i is exceeded.

[0142] In the combustion boiler 10, the furnace 12 and associated passes (horizontal pass 15 and back pass 16) define a flue gas flow path. The furnace 12 and the passes 15, 16 have a number of heat transfer surfaces 21.sub.i in the flue gas flow path. The combustion boiler 10 also has measurement instrumentation to monitor current load Q.sub.h of the combustion boiler, and further measurement instrumentation to currently monitor process data.

[0143] The control system (DCS 201, and EDGE server 203, or DCS 201 remote process intelligence system 205, possibly, under the participation of the EDGE server 203) is configured to carry out the boiler control method.

[0144] The EDGE server 203 may be configured to process the real-time measurement results for currently monitored process data and/or current load, namely, by filtering, averaging, and/or computing trends.

[0145] The control system may be configured to carry out the method step (b) to determine the current computational maximum boiler momentary load Q.sub.h,max locally at the combustion boiler 10, and/or to send data to a remote, preferably cloud-based (such as, computation cloud 206), computing system (such as, process intelligence system 205) which is configured to carry out the method step (b) and return the current computational maximum boiler momentary load Q.sub.h,max to the control system. The control system may then use the display/monitor to indicate the information, such as in method step (c), to the boiler operator, such as, by displaying the information.

[0146] The EDGE server 203 may be configured to reduce amount of measurement data that is passed to the remote computing system.

[0147] A combustion boiler computation system comprises a group of combustion boilers 10, each combustion boiler 10 comprising a boiler control system (CS) comprising an EDGE server (203) system that is configured to process the real-time measurement results for currently monitored process data and/or current load, namely, by filtering, averaging, and/or computing trends, and sending the processed real-time measurement results to a remote computing system. The remote computing system is preferably a cloud-based computing system, configured to receive data processed from real time measurement results and to compute data using a numerical boiler model for each of the combustion boilers 10, and to return computation results for each of the combustion boilers 10. The boiler control system may be configured to adapt its function based on the computation results.

[0148] The computing system is preferably configured to find such a numerical value Q.sub.h, candidate for a current computational maximum boiler momentary load Q.sub.h,max for which at least one flue gas factor df.sub.i computed using currently monitored process data with a numerical model of the boiler that fulfills an acceptance condition, and selecting the numerical value Q.sub.h, candidate as the current computational maximum boiler momentary load Q.sub.h,max.

[0149] The boiler computation system may be configured to adapt or to calibrate a numerical model for a boiler using processed measurement data for the combustion boiler 10.

[0150] Alternatively, or in addition, the boiler computation system may be configured to adapt or to calibrate a numerical model for a combustion boiler 10 using processed measurement data collected also from other combustion boilers 10.

[0151] FIG. 5 shows a modification of the method shown in FIG. 4. Steps L1, L3, L7, L9 are the same as steps K1, K3, K9, K11, respectively, but, in step L5, the flue gas factors df.sub.i can be directly computed for all heat transfer surfaces 201: if the temperatures T.sub.FG,in,i are measured using the respective temperature sensors 21.sub.i, the back-calculation will not be necessary and thus the step K7 can be omitted in the method illustrated in FIG. 5.

[0152] FIG. 6 shows in step N1 the use of possible inputs to the numerical boiler model. In step N3 the Q.sub.h,max is computed numerically using the boiler model, and in step N5, the estimated maximum load Q.sub.h,max is presented to boiler operator via a specific user interface (UI), preferably, via display/monitor 202.

[0153] FIG. 7 shows boiler momentary load Q.sub.h and computed current computational maximum boiler momentary load Q.sub.h, max, as well as the effect of using the method according to the invention during a test period. During the ten day test period, the 120 MW.sub.th boiler power obtained in average a 3 to 6 MW.sub.th higher load as outside the test period. FIG. 8 illustrates the ten day test period in more detail.

[0154] In other words, in the boiler control method, the current computational maximum boiler momentary load Q.sub.h,max of the combustion boiler is estimated using a numerical model using determined fluidized bed combustion boiler operating parameters. The current boiler load Q.sub.h is computed using steam circuit measurement data.

[0155] Then, if the boiler load Q.sub.h is less than the current computational maximum boiler momentary load Q.sub.h,max, it is (i) indicated to the boiler operator that the boiler load may be increased, and/or ii) the boiler load is automatically increased. Alternatively or in addition, if the boiler load Q.sub.h is larger than the boiler maximum momentary load Q.sub.h,max, it is (i) indicated to the boiler operator that the boiler load exceeds the boiler maximum momentary load, and/or (ii) the boiler load is automatically reduced.

[0156] It is obvious to the skilled person that, along with the technical progress, the basic idea of the invention can be implemented in many ways. The invention and its embodiments are thus not limited to the examples and samples described above but they may vary within the contents of patent claims and their legal equivalents.

[0157] In addition, or instead of using above mentioned specific empirical equations, it is possible to utilize artificial intelligence tools and/or neural network in the numerical model computations.

[0158] In the claims that follow and in the preceding description of the invention, except where the context requires otherwise due to express language or necessary implication, the word comprise or variations such as comprises or comprising is used in an inclusive sense, i.e., to specify the presence of the stated feature, but not to preclude the presence or addition of further features in various embodiments of the invention.