Method for detecting heat exchanger tube failures and their location when using input/loss performance monitoring of a recovery boiler

Abstract

This invention relates to a recovery boiler as used by the pulp and paper industry burning black liquor, and, more particularly, to a method for rapid detection of tube failures and the location of the affect heat exchanger within the recovery boiler, without need for direct instrumentation, thereby preventing more serious equipment damage, preventing boiler explosion, preventing injury to operators and minimizing repair time on the affected heat exchanger. This method is applicable to Input/Loss methods of monitoring recovery boilers.

Description

[0004] This invention relates to a recovery boiler as used by the pulp and paper industry burning black liquor, and, more particularly, to a method for rapid detection of tube failures and the location of the affect heat exchanger within the recovery boiler, without need for direct instrumentation, thereby preventing serious equipment damage, preventing boiler explosion, preventing injury to operators and minimizing repair time on the effected heat exchanger.

BACKGROUND OF TH INVENTION

[0005] Although especially applicable to “The Input/Loss Method” as installed at recovery boilers burning black liquor, this invention may also be applied to any other of the “Input/Loss methods” installed at any thermal system burning a fossil fuel. The teachings of this invention may be implemented for monitoring of any thermal system burning a fossil-fuel, or a thermal system burning a mix of fossil fuels and inorganic fuels. Such monitoring is assumed to be conducted in a continuous manner (i.e., on-line), processing one monitoring cycle after another, each cycle includes determining stoichiometric balances of the combustion process and, specifically, the fuel's chemistry, heating value, boiler efficiency, system efficiency and evaluation for possible tube failures. Specifically, The Input/Loss Method and its associated technologies are described in the following U.S. Patents and patent applications (cited above): U.S. Pat. No. 6,522,994 (hereinafter termed '994), U.S. Pat. No. 6,584,429 (hereinafter termed '429), U.S. Pat. No. 6,560,563 (hereinafter termed '563), U.S. application Ser. No. 10/087,879 (hereinafter termed '879), U.S. Pat. No. 6,651,035 (hereinafter termed '035) and U.S. application Ser. No. 10/131,932 (hereinafter termed '932); and in any of their related provisional patent applications and Continuation-In-Parts. One of the Input/Loss methods, a rudimentary method, is described in U.S. Pat. No. 5,367,470 issued Nov. 22, 1994 (hereinafter termed '470), and in Pat. No. 5,790,420 issued Aug. 4, 1998 (hereinafter termed '420).

[0006] Conventional steam generators and recovery boilers having gas-to-working fluid heat exchangers, may be prone to tube leaks of their working fluid (typically water as liquid or steam). These tube leaks represent a potential for serious physical damage to heat exchangers due to pipe whip (i.e., mechanical movement) and/or steam cutting of metal given high leakages flowing at critical velocities. In some recovery boilers, pressures of the working fluid may exceed 2300 psia Given failure of a heat exchanger tube, such fluid will experience many times critical pressure ratio as it expands into the combustion gases; that is, mixing with the products of combustion at essentially atmospheric pressure. When undetected, the damage from such tube failures may range from $2 to $10 million/leak forcing the system down for major repairs. If detected early, tube failures may be repaired before catastrophic damage, such repairs lasting only several days and costing a fraction of the cost associated with late detection and catastrophic damage. Repair times may be further reduced if the location of the heat exchanger which has the leak is identified before repairs are initiated.

[0007] However, an unique situation found with recovery boilers is associated with the pulp producing process to which they are integrated. This process involves first de-barking and chipping wood; then digesting the wood in an aqueous solution of NaOH and Na2SO4 (or other sodium-based compounds), forming a “white liquor”; then heating the brew; then separating the pulp from the spent liquor, the spent liquor is termed “black liquor” which consists of organics, water and inorganics (mostly sodium); and then the black liquor is burned in a recovery boiler. The essential function of a recovery boiler is the reduction in the furnace of sodium sulfate (Na2SO4, present in the black liquor) to sodium sulfide (Na2S). The efficiency of this sulfur reduction process is gauged by a “Reduction Efficiency” parameter. Heat from combustion of the fired organics, originating from the wood digesting process, generates steam. Black liquor inorganics, after reduction, are collected at the bottom of the furnace as a molten smelt, removed and recycled to recover sodium. Given a high Reduction Efficiency, smelt principally consists of Na2CO3, Na2S, inerts and free carbon.

[0008] The problem of tube failures in recovery boilers, in addition conventional problems cited above, is when water comes in contact with the molten smelt (typically at over 1400 F, with heavy concentrations of sodium); explosion is likely and may occur within minutes after tube failure. Recovery boiler explosions have dogged the pulp and paper industry since inception of the pulp producing process (i.e., called the Kraft process). Recovery boiler explosions injure and kill people every year. From 1948 through 1990 the industry recorded 140 recovery boiler explosions, three-quarters of which were smelt-water explosions. To place emphasis on the problem, the industry ranks explosions by severity: by definition just a “moderate explosion” keeps the plant off-line from 10 to 50 days; whereas a severe explosion keeps the plant off-line more than 50 days (typically lasting more than 120 days).

[0009] As common with conventional steam generators, tube failures in recovery boilers are typically caused by one the following general categories:

[0010] Weld failure of heat exchanger tubes;

[0011] Metallurgical damage caused by hydrogen absorption in the metal resulting in either embrittlement or the formation of non-protective magnetite;

[0012] Caustic gouging caused by the presence of free hydroxide in the water;

[0013] Corrosion-fatigue damage from the water-side of the tube, compounded by stress;

[0014] Corrosion damage caused by impacts from solid ash particles;

[0015] Fatigue failure caused by oxidation and/or mechanical movement, compounded by stress;

[0016] Overheating (e.g., from tube blockage) causing local creep; and

[0017] Physical damage from steam cutting and/or mechanical movement associated with another failed tube in the same locale.

[0018] Commonly, the physical leak initiates as a relatively small penetration, although initial breaks may also occur. For reference and further discussion see: Chapter 18, “Failure Analysis and In-Service Experience—Fossil Boilers and Other Heat Transfer Surfaces” of The ASME Handbook on Water Technology for Thermal Power Systems, P. Cohen, Editor, The American Society of Mechanical Engineers, New York, N.Y., 1989; and J. Gommi, “Root Causes of Recovery Boiler Leaks”, 1997 Engineering and Papermakers Conference, TAPPI Proceedings, available from TAPPI Press as product code ENG97509, Atlanta, Ga.

[0019] Present industrial art associated with conventional steam generators and recovery boilers have practiced the detection of tube failures using one or more of six general methods: 1) operator interface; 2) acoustic monitoring; 3) water balance testing; 4) monitoring of effluent moisture using instrumentation located at the system's effluent boundary (i.e., Stack); 5) monitoring the concentration of chemicals added to the working fluid whose change is sensitive to leakage; and 6) through use of artificial neural network technologies. Operator interface involves the use of his/her knowledge, experience, listening skills, and visual skills using remote cameras. However, all operators do not have the same high skill-set required. Acoustic devices detect the unique noise created by fluids at high velocities. However, acoustic devices rarely work in large steam generators, are expensive and require benchmarking with known acoustical signatures. Water balance testing may be conducted periodically on the entire system through which large water losses due to tube failures might be discovered. However, water balance testing is expensive, insensitive to small leaks, and typically may not be conducted at sufficient frequency to prevent serious damage. The use of an effluent moisture instrument has been shown to be sensitive to tube failures. Effluent moisture instrumentation may not differentiate between originating sources of water (e.g., between high humidity in the combustion air, or high fuel water, or changing fuel water, or a tube leakage). However, it might be practical to detect tube leakage by monitoring the difference in signals, or the rate of change of the difference in signals, between an effluent moisture instrument and one monitoring ambient air. Note that typical black liquor fuels contain up to 35% fuel water (approximately the same amount of water as found in some Powder River Basin (PRB) coals); thus changes in effluent moisture, even referenced to an ambient measurement, may be insensitive to small tube failures.

[0020] Monitoring the concentration of chemicals added to the working fluid operates by making a chemical mass balance on the working fluid based on a combination of flow measurements and chemical concentration measurements. Computed is a mass balance of a specific stable and non-volatile species (such as phosphate or molybdate) which has been uniquely added to the working fluid of the boiler. Basically a foreign chemical is injected into the working fluid; when a tube leak occurs the concentration of the chemical will change, thus detection. This method, developed by Burgmayer, Hong and Gunther, is described in U.S. Pat. No. 6,484,108 issued Nov. 19, 2002. A similar method is described in U.S. Pat. No. 3,522,008 issued Jul. 28, 1970. A similar method is also described in U.S. Pat. No. 5,320,967 issued Jun. 14, 1994. None of these methods involve combustion gases nor any stoichiometric balance involving the combustion process. The most serious limitation of these methods is their lack of sensitivity. Burgmayer's, Table 4 presents results of actual tests indicating that for a boiler producing 500,000 lb/hr of steam, detection of a 0.76% leakage (3,800 lb/hr) took 45 minutes, while a 2.8% leakage (14,000 lb/hr) took 15 min. to detect. Such sensitivities are not adequate to safeguard operators and equipment.

[0021] U.S. Pat. No. 6,192,352 by Alouani issued Feb. 20, 2001 (and his U.S. patent application Ser. No. 2001/0001149 of May 10, 2001) discloses a method to detect tube failures using artificial neural network and fuzzy logic technology (ANN). No where in Alouani's patent is explicit thermodynamic modeling taught. Alouani's patent teaches that ANN technology may learn to predict tube failures through recognition of patterns in raw data. Such raw data may include coal pulverizer flow (fuel flow), boiler drum pressure, reheat temperature, burner tilt positions, etc. The disadvantage to this method is that it requires a database from which it may learn. Such a database, as must be associated with an actual tube leak, is not frequent, inconvenient and could not contain a defined tube leakage flow rate (explicit system mass balances are not mentioned nor inferred by Alouani). In addition, Alouani's patent FIG. 5 indicates that a number of days is required for his system to detect a tube leak. In a survey of critical tube leaks in recovery boilers, it was found that approximately half of those leaks for which the time between leak initiation and explosion was known, the explosion occurred within 15 minutes after leak initiation; approximately 75% of the recorded explosions occurred within 30 minutes. This survey's reference is: D. G. Bauer and W. B. A. Sharp, “The Inspection of Recovery Boilers to Detect Factors That Cause Critical Leaks”, TAPPI Journal, Sep. 1991, TAPPI Press, Atlanta, Ga. Although there are more than a half-dozen vendors offering one or more of the six tube leak detection methods, in practice all known methods suffer serious short-comings and are not reliable in detecting early tube failures.

[0022] The patents '470 and '420 make no mention of heat exchanger tube failures nor their detection. Although the technologies of Patents '994, '429 and '563 generally support this invention, they make no mention of tube failure detection nor their location. Applications '879, '035 and '932 support this invention directly. Although the methods of '879, '035 and '932 are useful, the present invention further improves these methods and applies them to recovery boilers. There is no established art directly related to this invention; there is clear need for early detection of tube failures and to determine their location within the recovery boiler system.

SUMMARY OF THE INVENTION

[0023] This invention relates to a recovery boiler as used by the pulp and paper industry burning black liquor, and, more particularly, to a method for rapid detection of tube failures and the location of the affect heat exchanger within the recovery boiler, without need for direct instrumentation. Tube failures are detected through use of combustion stoichiometrics, in combination with an ability to correct effluent data through use of optimization procedures. The location of the failure within a recovery boiler is determined through energy balances, high accuracy boiler efficiency and iterative techniques. Further, this invention teaches how the stoichiometric mechanism of tube failure may be identified and reported to the system operator. This invention addresses the deficiencies found in all present detection methods.

[0024] Effluent water concentration (at the Stack) may consist of any one or all of the following sources of working fluid (assuming the working fluid is water): heat exchanger tube leaks; water added at the point of combustion (e.g., steam used to atomize fuel); pollutant control processes resulting in the net in-flow of water; and/or soot blowing processes using water to clean heat exchanger surfaces (commonly used in coal-fired systems). These sources of working fluid are in addition to: water formed from the combustion of hydrocarbon fuels; free water born by the fuel; and moisture carried by combustion air including air leakage. All such sources of effluent water are addressed by this invention through combustion stoichiometrics in combination with an ability to correct effluent data through use of optimization procedures.

[0025] This invention adds to the technology associated with Input/Loss methods. Specifically The Input/Loss Method has been applied through computer software, installable on a personal computer termed a “Calculational Engine”, and has been demonstrated as being highly useful to the operators of fossil-fired systems. The Calculational Engine receives data from the system's data acquisition devices. The Calculational Engine's software consists of the EX-FOSS, FUEL and HEATRATE programs described in '994 and '429, and in FIG. 2 herein, and the ERR-CALC program described in '879 (also employed in '035 and '932) and in FIG. 3 herein. ERR-CALC now incorporates the teachings of this invention. The Calculational Engine continuously monitors system efficiency on-line, i.e., in essentially real-time, as long as the thermal system is burning fuel. The application of this invention to The Input/Loss Method significantly enhances the system operator's ability to predict tube failures and reduce outage time required for repair.

[0026] The present invention provides a procedure for determining tube leaks in a recovery boiler using combustion stoichiometrics in combination with an ability to correct effluent data such that consistent fuel chemistry is computed.

[0027] The present invention teaches the mechanism of how a tube failure has been detected stoichiometrically, such detection being important to the system operator. Also, the present invention teaches how the location of a failed tube may be determined.

[0028] Other objects and advantages of the present invention will become apparent when its general methods are considered in conjunction with the accompanying drawings and the related inventions of '994, '429, '563, '879, '035 and '932.

[0029] The principles of this invention have been reduced to practice and installed for demonstration at two conventional power plants to determine the operability and functionality of this invention. These demonstrations have produced outstanding results demonstrating several identified tube failures. Furthermore, testing at a power plant burning PRB fuel similar to black liquor in its fuel water content, in which water soot blowing flow was removed from the Calculational Engine's input (bZ=0.0) which then computed an emulated tube leakage (bT>0.0) with outstanding agreement to the measured soot blowing flow; see FIG. 4.

[0030] ///

[0031] ///

[0032] ///

[0033] ///

[0034] ///

BRIEF DESCRIPTION OF THE DRAWINGS

[0035] FIG. 1 is a schematic representation of a recovery boiler illustrating the application of stoichiometric relationships, and also contains definitions of some of the important terms used herein.

[0036] FIGS. 2A and 2B is a block diagram of the general interactions and functions of the computer programs ERR-CALC, FUEL, EX-FOSS and HEATRATE used to implement this invention; herein collectively referred to as FIG. 2. FIG. 2 illustrates the “Fuel Iterations” involving FUEL, EX-FOSS and HEATRATE.

[0037] FIG. 3 is a block diagram of the principal functions of the error analysis computer program ERR-CALC which determines optimized Choice Operating Parameters, one of which is the tube failure flow rate incorporating the teachings herein.

[0038] FIG. 2 and FIG. 3 relate the interactions of the computer programs which implement this invention. The majority of the teachings of this invention are implemented in the ERR-CALC program. The FUEL, EX-FOSS and HEATRATE programs remain substantially unchanged as are used to implement the teachings of '994 and 3 429. However the FUEL, EX-FOSS and HEATRATE programs employ the results from ERR-CALC, including its calculated tube leakage flow rate, and thus assess the impact of such leakage on the thermal system in terms of boiler efficiency, fuel flow and system efficiency. Further, through energy balances on the steam generator's working fluid, and use of iterative procedures involving these programs, determination is then made within HEATRATE as to which heat exchanger within the steam generator contains the failed tube.

[0039] FIG. 4 is a plot of a sensitivity test of this invention. This test employed water soot blowing flow, entering the combustion space of a fossil-fired steam generator, to emulate a tube failure. The soot blowing flow, normally an input to the Calculational Engine, was left an unknown to be computed by the methods of this invention as a tube failure.

DESCRIPTION OF THE PREFERRED EMBODIMENT

[0040] To assure an appropriate teaching of this invention, its description is divided by sub-sections. The first two present nomenclature, definitions of equation terms, typical units of measure, and meaning of terms used herein (such as Choice Operating Parameters and System Effect Parameters), encompassing Paragraphs 0027 through 0043. The remaining sub-sections, representing the bulk of the teachings, are divided into four general groups:

[0041] 1) the first group presents system stoichiometrics applied to recovery boilers and the determination of fuel chemistry based on effluents, these teachings support all subsequent disclosures herein (encompassing Paragraphs 0044 through 0062, employing equations numbered less than one-hundred);

[0042] 2) the next group presents the determination of boiler efficiency for a black liquor-fired boiler as it influences both tube failure flow rate and determining tube failure location (encompassing Paragraphs 0063 through 0075, employing equations in the one-hundreds);

[0043] 3) the next group teaches how a tube failure may be detected based on an ability to correct Choice Operating Parameters using multidimensional minimization techniques, this ability being dependent, in part, on system stoichiometrics, the computed fuel chemistry and boiler efficiency (encompassing Paragraphs 0076 through 0108, employing equations numbered in the two-hundreds); and

[0044] 4) the last group teaches how both the tube leakage flow rate and its location in the steam generator are determined using, as a foundation, the preceding teachings (encompassing Paragraphs 0109 through 0119, employing equations numbered in the three-hundreds). The remaining paragraphs present a conclusion, THE DRAWINGS and related teachings. Teachings of multidimensional minimization techniques, as directly applicable to this invention are also presented in '879, '035 and '932. The present invention expands the utility of Input/Loss methods to recovery boilers, and specifically builds upon and expands the utility of The Input/Loss Method described herein and in '994, '429, '879, '035 and '932, and in '563 as it teaches the L Factor. The methods described in '563 teach the foundations of the L Factor used in multidimensional minimization techniques; the L Factor is further expanded as taught herein to encompass black liquor fuels used in recovery boilers.

[0045] Definitions of Equation Terms with Typical Units of Measure

[0046] Stoichiometric Terms:

[0047] a=Moles of combustion O2 input to the system; moles/base.

[0048] aβ=O2 entering with system air leakage (typically via the air pre-heater); mole/base.

[0049] aDRY-theor=Moles of combustion O2 input to the system required for theoretical combustion associated with Dry (water free) fuel; moles/base.

[0050] AAct=Concentration of O2 in combustion air local to (and entering) the system; molar fraction.

[0051] bA=Moisture in the entering combustion air; directly proportional to the ambient air's specific humidity, ω (lb-water/lb-dry air); moles/base.

[0052] =ωa (1.0+φAct) NDRY-AIR/NH2O

[0053] bAβ=Moisture entering with system air leakage; moles/base.

[0054] bT=Moles of tube leakage; i.e., water in-leakage entering and mixing with the combustion gases from leaks in heat exchangers; moles/base.

[0055] bZ=Moles of known water in-leakage entering and mixing with the combustion gases not related to heat exchanger leaks (bT); moles/base.

[0056] DAct=Total effluent CO2 at the system's boundary (i.e., Stack); moles/base.

[0057] =(d−σNkBL)

[0058] g=Calculational effluent O2 at the system's boundary associated with a hypothetically segregated organic component of the black liquor, without air leakage; moles/base.

[0059] GAct=Total effluent oxygen at the system's boundary burning black liquor fuel; moles/base.

[0060] =(g+2σNTNS+σNkBL/2)+aβ

[0061] j=Calculational effluent H2O at the system's boundary, without air leakage; moles/base.

[0062] JAct=Total effluent water at the system's boundary (j+bAβ); moles/base.

[0063] Jtheor=Total effluent water at the boundary based on theoretical combustion; moles/base.

[0064] kBL=Effluent SO2 at the system's boundary from the reaction of Na2SO4 with CO2 obtained at the system boundary; moles/base.

[0065] ni=Molar quantities of dry gaseous effluents of combustion at the system boundary without air leakage; specifically those products associated with the following quantities: DAct, (g+2σNTNS+σNkBL/2), h, eAct, mAct, σNkBL and (1.0−σN)xα6; note: Σni≡100 moles of dry gaseous effluent at the Stack is the assumed calculational “base” for Eq.(19BL), see FIG. 1; moles/base.

[0066] nii=Molar quantities of non-gas products of combustion at the system boundary without moisture associated with air leakage, specifically those products associated with the following quantities: j, xα10, σN(TNO+kBL), σNTKO, σNTNA, σNTNS, σNTNC, and v; see FIG. 1 and Eq.(19BL); moles/base.

[0067] Nk=Molecular weight of compound k.

[0068] RAct=Ratio of moles of dry gas from the combustion process before entering the air pre-heater to the diluted gas leaving, typically: (Moles of CO2 entering the air pre-heater)/(Moles of CO2 leaving the air pre-heater), defined as the Air Pre-Heater Leakage Factor; molar fraction.

[0069] R′Act=(Moles of O2 entering the air pre-heater)/(Moles of O2 leaving the air pre-heater).

[0070] TKO=Moles of effluent K2CO3 if burning black liquor (σN=1.0); moles/base.

[0071] TNA=Moles of effluent Na2SO4 if burning black liquor (σN=1.0); moles/base.

[0072] TNC=Moles of effluent NaCl if burning black liquor (σN=1.0); moles/base.

[0073] TNS=Moles of effluent Na2S if burning black liquor (σN=1.0); moles/base.

[0074] (TNO+kBL)=Moles of effluent Na2CO3 if burning black liquor (σN=1.0); moles/base.

[0075] v=Moles of effluent free carbon, or its equivalence, typically found in the recovery boiler's smelt; moles/base.

[0076] x=Moles of As-fired fuel required per 100 moles of dry gaseous effluent; moles/base.

[0077] xtheor=Moles of As-Fired fuel associated with theoretical combustion; moles/base.

[0078] xDRY-theor=Moles of Dry fuel associated with theoretical combustion; moles/base.

[0079] xMAF-theor=Moles of Moisture-Ash-Free fuel associated with theoretical combustion; moles/base.

[0080] xα10=Inert matter in As-Fired fuel, the terms “inert” and “ash” are used interchangeably; mole-inert/base.

[0081] αk=As-Fired (wet-base) fuel chemistry constituent k per mole of fuel; Σαk=1.0, where: k=1,2,3,4,5,6,10,14,15,16; see Eq.(19BL) therein for terms; mole-k/mole-fuel.

[0082] αMAF-k=Moisture-Ash-Free (MAF) fuel constituent k per mole of MAF fuel; ΣαMAF-k=1.0, where: k=1,3,4,5,6,14,15,16; see Eq.(19BL) therein for terms; mole-k/mole-fuel.

[0083] β=Air Pre-Heater Dilution Factor (ratio of air leakage to true combustion air); molar fraction

[0084] β=≡100(RAct−1.0)/[a RAct(1.0+φAct)]

[0085] σN=Kronecker function: unity if black liquor (bearing at least sodium) is being employed in the fuel, otherwise zero; unitless.

[0086] φAct=Ratio of non-oxygen gases (N2 and Ar) to oxygen in the combustion air; molar ratio.

[0087] φAct≡(1.0−AAct)/AAct

[0088] φRef=Reference ratio of non-oxygen gases (principally N2 and Ar) to oxygen in the combustion air, taken as 3.7737245; molar ratio.

[0089] ΦR=Reduction Efficiency; see Eq.(19BL); molar ratio.

[0090] ≡TNS/(TNA+TNS)

[0091] ΦS=Sulfur to Smelt ratio; see Eq.(19BL); molar ratio.

[0092] ≡(TNA+TNS)/[σNxα10+σN(TNO+kBL)+σN(TKO+TNA+TNS+TNC)+V]

[0093] ΦN=Sodium to Carbon ratio in the fuel; see Eq.(19BL); molar ratio.

[0094] ≡αMAF-14/αMAF-4

[0095] Multidimensional Minimization Terms: F ({right arrow over (x)})=Objective function, a functional relationship using the independent variables {right arrow over (x)}; unitless. ƒ( )=> Indicates a general functional relationship; for example, the expression: HHVk3=ƒ[fuel chemistry({right arrow over (Λ)})], means that HHVk3 is a function of fuel chemistry (which in-turn is a function of the vector {right arrow over (Λ)}).

[0096] Ci=Correction factor to be applied to an initial Choice Operating Parameter i; see Eqs.(211S) through (220) for nomenclature, e.g., C1S is the correction factor for Stack CO2 as referring to Eq.(211S), etc.; unitless.

[0097] HHVk3=Higher heating value as used by the minimization techniques as a System Effect Parameter, here subscript k3 refers to either an uncorrected As-Fired (HHVAF), Dry or MAF heating value; Btu/lbmAF, Btu/lbmDRY or Btu/lbmMAF.

[0098] HHVk3-Ref=Higher heating value used as a Reference System Effect Parameter; Btu/lbmAF, Btu/lbmDRY or Btu/lbmMAF.

[0099] J0=Bessel function of the first kind of order zero.

[0100] J1=Bessel function of the first kind of order one.

[0101] LCO2=L Factor for dry effluent CO2 used by the minimization techniques as a System Effect Parameter; (lbm effluent CO2)/million-BtuFuel.

[0102] L′Fuel=L Factor for dry gaseous effluents used by the minimization techniques as a System Effect Parameter; (lbm dry gaseous effluent)/million-BtuFuel.

[0103] Lk1=Generic L Factor, here subscript k1 refers to either L′Fuel or LCO2.

[0104] LCO2-Ref=L Factor for dry effluent CO2 used as a Reference System Effect Parameter; (lbm effluent CO2)/million-BtuFuel.

[0105] L′Fuel-Ref=L Factor for dry gaseous effluents used as a Reference System Effect Parameter; (lbm dry gaseous effluent)/million-BtuFuel.

[0106] mAF=Fuel flow rate, an As-Fired mass flow quantity (i.e., wet with water and fuel mineral matter), as may be computed by Input/Loss methods; also may be used by minimization techniques as a System Effect Parameter; lbmAF/hour.

[0107] mAF-PLT=The system's measured fuel flow, an As-Fired quantity (i.e., wet with water and fuel mineral matter), also termed the system's “indicated fuel flow”; also may be used as a Reference System Effect Parameter; lbmAF/hour.

[0108] mT=Tube leakage flow rate; i.e., mass flow rate of water in-leakage entering the combustion gas path from leaks in a heat exchanger, a Choice Operating Parameter; lbm/hour.

[0109] ML=Dilution Factor applied to System Effect Parameter Lk1; ML>0.0; unitless.

[0110] MW=Dilution Factor applied to System Effect Parameter mAF; MW>0.0; unitless.

[0111] MH=Dilution Factor applied to System Effect Parameter HHVk3; MH>0.0; unitless.

[0112] MT=Dilution Factor applied to System Effect Parameter WFH2O; MT>0.0; unitless.

[0113] Si=Scaling factor for the independent variable xi; reciprocal units of measure of Λi.

[0114] si=Pre-scaling factor used to adjust Si; unitless.

[0115] WFj=As-Fired fuel mass fraction of substance j; fraction. WFH2O=As-Fired fuel water mass fraction (also termed WF2), used as a System Effect Parameter; fraction.

[0116] WFH2O-Ref=Reference As-Fired fuel water mass fraction used by the minimization techniques as a Reference System Effect Parameter; fraction.

[0117] {right arrow over (x)}=Vector of independent variables, {right arrow over (x)}=(x1, x2, X3, . . . ), as based on scaled Choice Operating Parameters (not to be confused with the term for moles of As-fired fuel, x); unitless.

[0118] Λi=Choice Operating Parameter i, see the specific parameter for units of measure, and Eqs.(211S) through (220) for definitions.

[0119] {right arrow over (Λ)}=Vector of Choice Operating Parameters, which is user selected; for example, one selection might include: {right arrow over (Λ)}=(Λ1S, Λ2S, Λ3, Λ6, Λ7B, Λ9); see Eqs.(211S) thru (220).

[0120] Λ0−i=Initial Choice Operating Parameter i, before application of a minimization technique and based on the system's raw instrumentation signal, a previous converged solution, an estimate, and/or as otherwise obtained.

[0121] ΛF−i=Converged (final) Choice Operating Parameter i, after application of a minimization technique to Λ0−i; and, thus corrected and applicable to all system thermal analyses.

[0122] Quantities Related to System Terms:

[0123] AF=Air/Fuel ratio defined by indicated air flow and mAF-PLT; unitless mass ratio.

[0124] BBTC=Energy flow to the working fluid from combustion gases (w/o tube leakage); Btu/hr.

[0125] Cp=Heat capacity; Btu/lbm-R.

[0126] hj=Specific enthalpy of substance j; Btu/lbm.

[0127] HBC≡Firing Correction; Btu/lbmAF.

[0128] HHVP=As-Fired higher heating value, based on HHVAF and used in system evaluations as corrected for a constant pressure process; Btu/lbmAF.

[0129] HNSL≡Non-Chemistry & Sensible Heat Losses; Btu/lbmAF.

[0130] HPRAct≡Enthalpy of Products from actual combustion (HHV- or LHV-based); Btu/lbmAF.

[0131] HPRIdeal≡Enthalpy of Products from ideal combustion associated with an obtained heating value at TCal (HHV- or LHV-based); Btu/lbmAF.

[0132] HRXAct≡Enthalpy of Reactants associated with actual firing conditions (HHV- or LHV-based); Btu/lbmAF.

[0133] HRXCal≡Enthalpy of Reactants associated with an obtained heating value at TCal (HHV- or LHV-based); Btu/lbmAF.

[0134] HR=System heat rate (HHV-based, HRHHV; or LHV-based, HRLHV); Btu/kWh.

[0135] HSL≡Stack Losses (HHV- or LHV-based); Btu/lbmAF.

[0136] LHVAF=Lower heating value based on the measured or calculated higher heating value (HHVAF); Btu/lbmAF.

[0137] LHVP=As-Fired lower heating value, based on LHVAF and used in system evaluations as corrected for a constant pressure process; Btu/lbmAF.

[0138] mTΔh=Energy flow from tube leakage; Ah is the enthalpy difference between the last heat exchanger effected by the leakage (typically the Final Superheater or Reheater) and the heat exchanger producing the leakage; Btu/hr.

[0139] PAmb≡Ambient pressure local to the system, psiA.

[0140] QSAH=Energy flow delivered to steam/air heaters, Btu/hr.

[0141] T=Temperature; F.

[0142] TAmb≡Ambient temperature local to the system, F.

[0143] TCal≡Calorimetric temperature to which heating value is referenced, F.

[0144] TStack≡Boundary temperature of the system effluents, defines the “Stack”; F.

[0145] WFD=Brake power associated with in-flow stream fans (e.g., Forced Draft fans); Btu/hr.

[0146] WID=Brake power associated with out-flow stream fans (e.g., Induced Draft fans); Btu/hr.

[0147] Woutput=Gross power generated from a power plant; kWe.

[0148] ηSYS=System efficiency (HHV-based, ηSYS-HHV; or LHV-based, ηSYS-LHV); unitless

[0149] ηB=Boiler efficiency (HHV-based, ηB-HHV; or LHV-based, ηB-LHV); unitless.

[0150] ηC=Combustion efficiency (HHV-based, ηC-HHV; or LHV-based, ηC-LHV); unitless.

[0151] ηA=Boiler absorption efficiency; unitless.

[0152] Subscripts and Abbreviations:

[0153] Act=Actual value obtained from the operating thermal system.

[0154] AF=As-Fired fuel at the thermodynamic boundary (i.e., wet with water and inert matter).

[0155] BL=Black liquor fuel.

[0156] DRY=Dry chemical base (i.e., free of water).

[0157] MAF=Moisture-Ash-Free chemical base (i.e., free of water and free of inert matter).

[0158] Ref=Reference value.

[0159] T=Tube failure associated with a heat exchanger leakage.

[0160] theor=Refers to conditions associated with theoretical combustion.

[0161] Meaning of Terms

[0162] The words “Operating Parameters”, as taken within the general scope and spirit of the present invention, mean common data obtained from a thermal system applicable to the thermodynamic understanding of that system. The following quantities may be included in the definition of Operating Parameters, they are not encompassing but considered typical of a minimum set of data required for such thermodynamic understanding: effluent CO2, O2, and SO2 concentrations determined at the Stack, or before the air pre-heater (Boiler side of the air pre-heater); the mass, wet-base ratio of the indicated combustion air flow at the system's fuel combustors, to the system's indicated fuel flow, termed AFAct (note that AFAct is obtained only for the determination of inerts as taught in '994); effluent H2O concentration measurement, or assumptions made (or as otherwise may be determined); effluent temperature measurement, that is the average temperature associated with the combustion gases at the system boundary (caution must be exercised in measuring non-stratified gas flows); the inlet/outlet ratio of CO2 (producing RAct as is preferred), CO or O2 (producing R′Act) across the air pre-heater where these ratios could be obtained on-line, off-line, based on periodic testing or judgement which are used for the determination of air pre-heater leakage; determination of fuel temperature at an appropriate system boundary; air psychrometric measurements, or as otherwise determined, at the system boundary (e.g., dry and wet bulb temperatures, or dry bulb and relative humidity, or dry bulb and dew point temperatures); quantities comprising the system's Firing Correction term, HBC; the discharge temperatures of the air as it exits each air heating or cooling device but before it reacts with the fuel (for example, such devices might include the air pre-heater, forced-draft fan, steam-to-air heater, etc.); the total energy flow deposition to the working fluid from the combustion gases (in typical units of measure of Btu/hr); and similar quantities. Regards the total energy flow deposition, for a typical recovery boiler, such measurements typically include feedwater flow to the steam generator, feedwater pressure and temperature, determination of the steam flow from the steam generator if different than the feedwater flow, steam pressure, steam temperature or quality (or assumed quality), and, if applicable, reheat flows, and reheat inlet and outlet pressures and temperatures. If employing a Reheater heat exchanger, determination of accurate reheat flows generally requires understanding of steam turbine flow distributions (involving high pressure turbine shaft seals, steam flows to feedwater heaters, turbine bypass leakages, attemperation spray flows and the like). The total realizable energy flow (BBTC−mTΔh) includes the effects of tube leakage flow rate, determined by this invention, as assigned to a particular heat exchanger. Specifically, when having determined there is a tube leak, the leaking heat exchanger's outlet flow of working fluid is reduced by the tube leakage, thus reducing the energy flow BBTC by mTΔh.

[0163] The words “Choice Operating Parameters”, as taken within the general scope and spirit of the present invention, mean any sub-set of Operating Parameters which directly impact system stoichiometrics, and thus may impact the determination of fuel chemistry. This invention assumes that Choice Operating Parameters may have error, said error may adversely effect the determination of fuel chemistry, but said error may be corrected as taught herein through optimization methods. In the Preferred Embodiment Choice Operating Parameters are selected by the user of this invention from an available set. This available set of Choice Operating Parameters includes the following ten: 1) effluent CO2 concentration measured at the Stack or Boiler; 2) H2O concentration measured, or as otherwise may be determined, at the Stack or Boiler; 3) the mass, wet-base ratio of the indicated combustion air flow at the system's fuel combustors, to the system's indicated fuel flow, the Air/Fuel ratio termed AFAct; 4) the Air Pre-Heater Leakage Factor, termed RAct; 5) the concentration of O2 in the combustion air local to the system, or as otherwise determined, termed AAct (leading to the determination of φAct); 6) the system's indicated limestone mass flow rate, termed mLS; 7) effluent O2 concentration measured at the Stack or Boiler; 8) mass flow rate associated with a heat exchanger tube leakage flow rate, termed mT; 9) Sodium to Carbon ratio found in the fuel (ΦN); and 10) the Reduction Efficiency (ΦR).

[0164] The words “Reference Fuel Characteristics”, as taken within the general scope and spirit of the present invention, mean the identification and development of intrinsic chemical relationships which allow resolution of fuel chemistry. Discussions of Reference Fuel Characteristics are presented in Paragraphs 0048 through 0058. For the Preferred Embodiment, Reference Fuel Characteristics include an average or typical fuel chemistry and associated MAF heating value, preferably based on historical data collections of ultimate analyses of the fuel's elemental composition, with fuel water and fuel inerts (typically reported as weight fractions summing to unity, leading to αk molar fractions), herein termed fuel chemistry constituents. If ultimate analyses are not available they may be approximated from proximate analyses as taught in Gill's text (cited below, not preferred). Reference Fuel Characteristics typically include a fuel hydrogen versus fuel carbon relationship, typically a MAF relationship; a relationship based on historical data. Reference Fuel Characteristics may also include relationships of: MAF oxygen versus MAF carbon; MAF sulfur versus MAF carbon; MAF potassium versus MAF sodium; and MAF chlorine versus MAF sodium. Reference Fuel Characteristics may also include the computed values of L′Fuel-Ref and LCO2-Ref, computed using the reference fuel chemistry. For black liquor fuel Reference Fuel Characteristics may also include the variables ΦR, ΦS and ΦN. Reference Fuel Characteristics may also include whether the variability of fuel water and fuel inert fractions in the As-Fired condition is predictable, or not. For any given fuel: fuel water may be held constant (including zero); fuel inerts may be held constant (including zero); functionalities may be observed and applied, for example, αMAF-10=ƒ (HHVMAF), and/or fuel water and/or fuel inerts may be treated as unknowns). The Preferred Embodiment for black liquor fuels is to treat fuel water as an unknown, and fuel inerts as a constant. All of these possible variations for the treatment of fuel carbon, fuel sodium, water and inerts may be included as a portion of Reference Fuel Characteristics. Reference Fuel Characteristics may also contain fitting constants associated with all correlations relating dependent fuel quantities to independent fuel quantities, typically the principle independent fuel quantity is fuel carbon. When applying the Preferred Embodiment of this invention, Reference Fuel Characteristics also contain reasonability limits (i.e., numerical minimum and maximum limits) of the following: computed fuel elemental constituents; fuel water fraction; fuel inert fraction; the correction factors being applied to initial Choice Operating Parameters; the maximum rates of change of these parameters; and the maximum rates of change of any or all of the selected Choice Operating Parameters. Such minimum and maximum limits are preferably based on engineering judgement supported by historical data collections of ultimate analyses of the fuel, historical experience of instrumentation on which the selected Choice Operating Parameters are based, and historical records of computed correction factors to the initial Choice Operating Parameters, termed Λ0−i;.

[0165] The words “System Effect Parameters”, as taken within the general scope and spirit of the present invention, mean any parameter of the thermal system or its fuel which directly impact the determination of system efficiency. In the most general sense System Effect Parameters include any parameter used in Eq.(331A), (331B), (332A) or (332B) which compute system heat rate and thus system efficiency. For the Preferred Embodiment, System Effect Parameters include the following four types of quantities: the L Factor (either L′Fuel or LCO2); the computed As-Fired fuel flow rate (mAF); the higher heating value (either HHVAF, HHVDRY or HHVMAF); and the As-Fired fuel water faction (WFH2O) which may be used to convert HHVDRY to HHVAF. The computed L Factor effects fuel chemistry which effects heating value and boiler efficiency, thus has an immediate impact on system efficiency. “Reference System Effect Parameters” are constant and targeted (i.e., desired) System Effect Parameters to which the System Effect Parameters are numerically driven by the minimization techniques through optimizing a selection of Choice Operating Parameters.

[0166] The words “Input/Loss methods”, as taken within the general scope and spirit of the present invention, mean any method or combination of methods in which one or more of the following parameters is determined based on Operating Parameters and a selection of Choice Operating Parameters: fuel flow, effluent flow, emission rates, fuel chemistry, fuel heating value, boiler efficiency, and/or system efficiency. In addition to '994, '429, '563, '879, '035 and '932, and their related provisional patent applications and Continuation-In-Parts, Input/Loss methods include the methods of '470 and '420. The words “The Input/Loss Method” refers specifically to the collection of technologies described in '994, '429, '879, '035 and '932, and in '563 as it teaches the L Factor, and to any relevant provisional patent applications and Continuation-In-Parts, in addition to the teachings disclosed herein.

[0167] As used herein, the words “Calculational Engine” refers to a computer in which software descriptive of The Input/Loss Method is installed.

[0168] As used herein, if used, the words “obtain”, “obtained”, “obtaining”, “determine”, “determined”, “determining”, “determination”, “establish”, “established” or “establishing” are defined as measuring, calculating, computing, assuming, estimating or gathering from a database.

[0169] As used herein, the words “monitoring” or “monitored” are meant to encompass both on-line monitoring (i.e., processing system data in real time) and off-line monitoring (i.e., computations involving static data). A “monitoring cycle” is meant one execution of the processes described in FIG. 2, which encompasses FIG. 3.

[0170] As used herein, the meaning of the words “smoke Stack” or “Stack” or “system boundary” are defined as the physical boundary of the thermal system where gaseous combustion effluents exit, entering the local environment; refer to 42 in FIG. 1, further discussed within THE DRAWINGS. Solid effluents not leaving the Stack (e.g., molten smelt 37) are referenced to the generic system's boundary 44 in FIG. 1.

[0171] As used herein, the meaning of the words “Boiler” or “Boiler Effluent” are defined as the region 35 in FIG. 1, or generically between the physical exit of the system's region 34 in FIG. 1 and entrance to its air pre-heater 36 in FIG. 1; see THE DRAWINGS.

[0172] As used herein, the meaning of the words “Fuel Iterations”, are defined in conjunction with a detailed description of FIG. 2 found within THE DRAWINGS.

[0173] As used herein, the meaning of the word “indicated” when used in the context of data originating from the thermal system is defined as the system's actual and uncorrected measurements of a physical process (e.g., pressure, temperature, mass flow, volumetric flow, density, and the like) whose accuracy or inaccuracy is not assumed. As examples, a system's “indicated fuel flow” or its “indicated limestone flow” denote system measurements the accuracy of which is unknown (they are “as-is”, with no judgement applied). Such indicated measurements are said to be either correctable or not. If not correctable, it may be that the associated computed value from Input/Loss methods tracks the indicated value over time (the indicated not being corrected per se). In the case of indicated limestone flow when used as a Choice Operating Parameter (Λ6), it is directly corrected as taught by this invention. In the case of indicated fuel flow when used as a System Effect Parameter, it may be shown that the computed fuel flow, mAF, tracks the indicated fuel flow, mAF-PLT, through adjustment of the Dilution Factor Mw.

[0174] As used herein, the meaning of the words “Reduction Efficiency” is defined by the variable ΦR; the meaning of “Sulfur to Smelt ratio” is defined by the variable Φs; and the meaning of “Sodium to Carbon ratio” in the fuel by the variable ΦN.

[0175] ///

[0176] ///

[0177] ///

[0178] System Stoichiometrics

[0179] Any study of the combustion of fossil fuels necessitates the formulation and use of a combustion equation. Combustion equations used by several Input/Loss methods are described in '470, '420 and '994 by their designated Eq.(29), in '429 by its Eq.(19), in '879 by its Eq.(19-corr), and in '035 and '932 by their Eq.(19B). This invention's methods are taught through a combustion equation defined by Eq.(19BL) herein. Through Eq.(19BL) stoichiometric terms are self-defined, including the use of the bT term important to this invention and inclusion of fuel components typically found in recovery boilers (hydrocarbons, sodium, potassium and chlorine). Eq.(19BL)'s nomenclature is unique in that brackets are used for clarity: for example, the expression “xα2[H2O]” means the moles of fuel water/base, algebraically simply xα2; the expression “σNxα14[Na2]” means the moles of fuel diatomic sodium/base, algebraically simply σNxα14; the expression “(d−σNkBL)Act[CO2]” means the effluent moles of carbon dioxide/base, algebraically simply (d−σNkBL)Act which defines the term DAct; etc. The stoichiometric base of Eq.(19BL) is 100 moles of dry Stack gas (i.e., at the thermodynamic boundary). Also, note that in combustion equations such as Eq.(19BL) describing commercial recovery boilers, that the assumption is typically made of possible extraneous water leakage into, and mixing with, the products of combustion. The flow of such water is assumed known (either measured or otherwise reasonably estimated) and is denoted by the symbol bZ. It is common industrial practice in the United States to use this symbol, bZ, to denote the quantity of known water in-leakage such as steam used to atomize fuel, water used for soot blowing, etc. Such bZ in-leakage is apart from water formed from combustion of hydrocarbon fuels, the xα5 term; apart from free water born by the fuel, xα2; apart from moisture carried by the combustion air and air pre-heater leakage, (1.0+β)bA; and apart from tube leakage, the bT term. 1$\begin{array}{cc}x\left[{\alpha }_1\left[N_2\right]+{\alpha }_2\left[H_{2}O\right]+{\alpha }_3\left[O_2\right]+{\alpha }_4\left[C\right]+{\alpha }_5\left[H_2\right]+{\alpha }_6\left[S\right]+{\alpha }_{10}\left[\mathrm{Inert}\right]+{\sigma }_N{\alpha }_{14}\left[{\mathrm{Na}}_2\right]+{\sigma }_N{\alpha }_{15}\left[K_2\right]+{\sigma }_N{{\alpha }_{16}\left[\mathrm{Cl}\right]}_{\mathrm{As}-\mathrm{Fired}\mathrm{Fuel}}\right]+{b_Z\left[H_{2}O\right]}_{\mathrm{In}-\mathrm{Leakage}}+{b_T\left[H_{2}O\right]}_{\mathrm{Tube}-\mathrm{Leakage}}+{\left[\left(1.0+\beta \right)\left(a\left[O_2\right]+a{\phi }_{\mathrm{Act}}\left[N_2\right]+b_A\left[H_{2}O\right]\right)\right]}_{\mathrm{Air}}=\left(d-{\sigma }_N{k}_{\mathrm{BL}}\right)\left[{\mathrm{CO}}_2\right]+\left(g+2{\sigma }_N{T}_{\mathrm{NS}}+{\sigma }_N{k}_{\mathrm{BL}}/2\right)\left[O_2\right]+h\left[N_2\right]+j\left[H_{2}O\right]+x{\alpha }_{10}\left[\mathrm{Inert}\right]+e_{\mathrm{Act}}\left[\mathrm{CO}\right]+m_{\mathrm{Act}}\left[\mathrm{NO}\right]+{\sigma }_N\left(T_{\mathrm{NO}}+k_{\mathrm{BL}}\right)\left({\mathrm{Na}}_2{\mathrm{CO}}_3\right]+{\sigma }_N{T}_{\mathrm{KO}}\left[K_2{\mathrm{CO}}_3\right]+{\sigma }_N{T}_{\mathrm{NA}}\left[{\mathrm{Na}}_2{\mathrm{SO}}_4\right]+{\sigma }_N{T}_{\mathrm{NS}}\left[{\mathrm{Na}}_{2}S\right]+{\sigma }_N{T}_{\mathrm{NC}}\left[\mathrm{NaCl}\right]+{\sigma }_N{k}_{\mathrm{BL}}\left[{\mathrm{SO}}_2\right]+\left(1.0-{\sigma }_N\right)x{\alpha }_6\left[{\mathrm{SO}}_2\right]+v\left[C_{\mathrm{Refuse}}\right]+{\left[\beta \left(a\left[O_2\right]+a{\phi }_{\mathrm{Act}}\left[N_2\right]+b_A\left[H_{2}O\right]\right)\right]}_{\mathrm{Ait}\mathrm{Leakage}}& \left(19\mathrm{BL}\right)\end{array}$

[0180] Resolution of Eq.(19BL) is had when all ni and nii quantities have been determined. Minor component terms of Eq.(19BL) are typically resolved either through direct measurement (e.g., for CO and NO), or assume zero values, or through obtained relationships. All Minor Components typically have only low parts-per-million concentrations thus have little impact. An exception to this is the v term describing refuse carbon found in the smelt. The true importance and functionality of Eq.(19BL) to The Input/Loss Method lies in the fact that consistency of molar balances is needed for successful system understanding, for conservation of mass flows and for resolution of fuel chemistry, these needed for detection of tube failures. For clarity the following major terms are associated with system stoichiometrics:

[0181] Total effluent (boundary) water≡JAct=j+bAβ

[0182] Boiler oxygen before air leakage (g′Act)≡RAct(g+2σNTNS+σNkBL/2)

[0183] Total effluent (boundary) oxygen≡GAct=(g+2σNTNS+σNkBL/2)+aβ

[0184] Total effluent (boundary) carbon dioxide≡DAct=(d−αNkBL)

[0185] Total effluents referenced to the boundary=Σni+Σnii+β(a+aφAct+bA)

[0186] Total effluents before air leakage, referenced upstream of the air pre-heater=RActΣni+RActΣnii

[0187] Dry combustion air without air leakage referenced to the boundary=(a+aφAct)

[0188] Wet combustion air without air leakage referenced to the boundary=(a+aφAct+bA) Dry air from air leakage found at the boundary=β(a+aφAct)

[0189] Total wet combustion air and air leakage found at the boundary=(1.0+β)(a+aφAct+bA).

[0190] Eq.(19BL) is unique in describing at least three features of critical importance when determining fuel chemistry using one of the Input/Loss methods. The critical features include: 1) its ability to address air pre-heater leakage through application of the Air Pre-Heater Leakage Factor, RAct, and through the Air Pre-Heater Dilution Factor, β; 2) the ability to describe effluent concentrations on either side of the air pre-heater, again through application of RAct; and 3) the use of an explicit φAct term allowing for variable O2 concentration in the system's local combustion air. Air pre-heater leakage dilutes all combustion effluents with moist air from the local environment, thus all important effluents H2O, CO2 and O2 used for this invention are altered. Furthermore, many times, although not always, a power plant's more precise effluent measurements may be found on the air pre-heater's inlet (economizer outlet or Boiler), and not at the air heater outlet (or Stack); refer to FIG. 1. Although most environmental regulation requires effluent measurements at the system's boundary, translation between the air heater inlet and outlet measurements is many times essential. Eq.(19BL) allows for such translation through the RAct term, defined above such that 100 moles of dry gas are computed both at the upstream and downstream locations of the air pre-heater; see “Boiler” of FIG. 1. Thus effluents may be used by the present invention either upstream or downstream of the air pre-heater; refer to the GAct and JAct terms defined above, allowing conversion between measurements with and without air leakage. For example, combustion gas conditions for oxygen and water upstream of the air pre-heater and after exiting the heat exchangers and combustion region, see FIG. 1, would employ the terms: (g+2σNTNS+σNkBL/2)RAct and jRAct. That is, one would actually measure (g+2σNTNS+σNkBL/2)RAct moles of dry O2 upstream of the air pre-heater and after exiting the heat exchangers and combustion region as based on 100 moles of dry gas found at that location. As a further example, the moles of dry air found upstream of the air pre-heater is given by: RAct a(1.0+φAct). Combustion gases downstream of the air pre-heater typically exit the system to the environment (i.e., Stack), in other words the gaseous effluent boundary of the system (100 moles of dry gas at the Stack includes air leakage). If limestone is injected into the combustion process to control effluent SO2 it will create additional effluent CO2; further, it could decrease the effluent H2O if the sulfate product is matrixed with water, CaSO4.zH2O. Thus such effects must be considered, and are taught in '994 referenced herein. Although the specific use of limestone is taught in '994 it is not included herein for clarity of presentation. Of course CO2 and H2O are two important effluents to the present invention. In addition to the basic stoichiometrics afforded, Eq.(19BL) allows numerous and obvious determinations of molar and mass ratios.

[0191] Based on these teachings and those presented in '994, the following further explains the importance of the Air Pre-Heater Leakage Factor, RAct, and the Air Pre-Heater Dilution Factor, β, their definitions and developments. Consider that air in-leakage associated with a fossil-fired system, and as commonly associated with in-leakage at the system's air pre-heater, is defined by the American Society of Mechanical Engineers' Performance Test Code 4.3 (1974) as the moist air leakage divided by the wet combustion gas. Typically the wet combustion gas is taken at the gas inlet of the air pre-heater (i.e., Boiler, or economizer outlet before the air pre-heater). That is, as defined herein using Eq.(19BL) nomenclature, noting that 100 moles of dry gas is the bases at the Boiler, is given by: 2$\begin{array}{cc}\mathrm{Wet}\mathrm{APH}\mathrm{Leakage}\equiv \frac{R_{\mathrm{Act}}\beta \left(a+a{\phi }_{\mathrm{Act}}+b_A\right)N_{\mathrm{MoistAir}}}{\left(100+R_{\mathrm{Act}}j\right)N_{\mathrm{WetGas}}}& \left(20\right)\end{array}$

[0192] where, as defined above:

[0193] RAct=(Moles of CO2 entering the air pre-heater)/(Moles of CO2 leaving the air pre-heater).

[0194] The expression for RAct is equivalent to [Moles of Boiler CO2] divided by [Moles of Stack CO2]. The Air Pre-Heater Dilution Factor is then developed by performing a total dry gaseous effluent molar balance at the Stack:

100 moles dry gaseous effluent at Stack=Σni+β(a+aφAct)  (21)

[0195] then solving for β: β=(100−Σni)/(a+aφAct). The stoichiometric base of Eq.(19BL) implies that 100 moles of dry gaseous effluent upstream of the air pre-heater (Boiler) is given by RActΣni, see FIG. 1; therefore: 3$\begin{array}{cc}\beta =\left(100-100/R_{\mathrm{Act}}\right)/\left[a\left(1.0+{\phi }_{\mathrm{Act}}\right)\right]\equiv 100\left(R_{\mathrm{Act}}-1.0\right)/\left[R_{\mathrm{Act}}a\left(1.0+{\phi }_{\mathrm{Act}}\right)\right].& \left(22\right)\end{array}$

[0196] If, instead of obtaining the ratio of CO2 across the air pre-heater, the ratio of O2 is obtained and, following the teachings herein and those of '994, the following may then be developed:

[0197] R′Act=(Moles of O2 entering the air pre-heater)/(Moles of O2 leaving the air pre-heater).

[0198] where, converting from R′Act to RAct, using algebraic manipulations results in, when measuring Stack O2: 4$\begin{array}{cc}{R}_{\mathrm{Act}}\equiv \frac{100-R_{\mathrm{Act}}^{\prime }{G}_{\mathrm{Act}}\left(1.0+{\phi }_{\mathrm{Act}}\right)}{100-G_{\mathrm{Act}}\left(1.0+{\phi }_{\mathrm{Act}}\right)}& \left(23\right)\end{array}$

[0199] If measuring Boiler O2 (for Eq.(24) termed g′Act): 5$\begin{array}{cc}{R}_{\mathrm{Act}}\equiv \frac{100R_{\mathrm{Act}}^{\prime }-R_{\mathrm{Act}}^{\prime }{g}_{\mathrm{Act}}^{\prime }\left(1.0+{\phi }_{\mathrm{Act}}\right)}{100R_{\mathrm{Act}}^{\prime }-g_{\mathrm{Act}}^{\prime }\left(1.0+{\phi }_{\mathrm{Act}}\right)}& \left(24\right)\end{array}$

[0200] There are, of course, a number of variations to these formulations, such as employing 100 moles of wet effluents at the Stack, thus replacing Eq.(21) with:

100 moles wet effluent at Stack=(Σni+j)+β(a+aφAct+bA)  (25)

[0201] or using an oxygen base for the wet effluents at the Stack, thus: (Σni+JAct)/a+β(1.0+φAct); or using a combustion equation which is based on a mole of fuel carbon (xα4); etc. What is important to this invention, important to The Input/Loss Method, and important to any of the Input/Loss methods, is that the Air Pre-Heater Leakage Factor (RAct) allows gaseous measurements to be employed on either side of the system air in-leakage. Typically, but not always, O2 is measured in the combustion gas stream inlet to the air pre-heater (Boiler), while CO2 is measured at the Stack (downstream from the air pre-heater).

[0202] Operators of recovery boilers will commonly determine certain intrinsic fuel and effluent ratios which assist in understanding the performance of the combustion process, i.e., performance ratios. Such ratios are determined through periodic fuel analysis, and analysis of the combustion residuals (also termed “smelt” for recovery boilers). Smelt is defined as the molten inorganic compounds leaving the furnace after combustion; using Eq.(19BL) nomenclature it consists of: xα10, σN(TNO+kBL), σNTKO, σNTNA, σNTNS, σNTNC and v. Common performance ratios include: the Reduction Efficiency (ΦR); the Sulfur to Smelt ratio (ΦS); and the Sodium to Carbon ratio in the fuel (ΦN). Reduction Efficiency judges the ability of the combustion process to further reduce Na2SO4. For theoretical combustion, ΦR=0.0; for complete reduction of Na2SO4, ΦR=1.0. The Sulfur to Smelt ratio is employed in the Preferred Embodiment of this invention. The Sodium to Carbon ratio in the fuel is a key parameter, determined periodically by laboratory analysis, as it allows stoichiometric inter-dependency between the organic and inorganic fuel components. Other such parameters may be formed which are similar to the ΦN ratio. As taught herein, what is required is some defined or otherwise obtained dependency between the organic and inorganic fuel compounds if fuel chemistry is to be determined based on gaseous effluent and smelt measurements. Although the Preferred Embodiment of this invention employs the ΦN ratio, on-line measurements of the As-Fired fuel and/or on-line measurements of the smelt could prove as useful as the ΦN ratio, and/or otherwise aid in its determination.

[0203] After establishing system stoichiometrics, the next stage of the process involves the recognition that a given fuel has an unique chemical composition, thus when burned will yield unique concentrations in its gaseous effluent. The gaseous effluent concentrations are used to compute fuel chemistry, with this chemistry heating value and boiler efficiency are computed, in turn this information allows the detection of tube failures and their location. However, key to this invention is that its methods recognize that gaseous effluent concentrations are not accurate and require correction to achieve stoichiometric consistency. The gaseous effluents from any fossil combustion process, including black liquor combustion, are N2, CO2, H2O, O2 and SO2. H2O, when effluent from combustion, is in its superheated phase, thus acting as a gas. The source of N2 is principally the air used to burn the fuel and has little chemical reactiveness, thus its sensitivity to the fuel's chemical composition is not significant. SO2 effluent concentrations are generally in the parts per million thus its impact has minor importance. SO2 effluent is formed from either Na2SO4 reduction or, if burning a pure hydrocarbon, directly from fuel sulfur (1.0−σN)xα6.

[0204] As an intrinsic chemical relationship, the relative concentrations of carbon, α4, and hydrogen, α5, found in any fossil fuel will have significant impact on the relative concentrations of CO2 and H2O found in the effluent. In addition, these effluents will be influenced by the following: O2 used to burn the fuel (i.e., the Air/Fuel ratio); fuel water, α2; in-leakage of water including tube leaks; and water in the combustion air. This implies that the molar fractions of CO2, H2O and O2 present in the effluent (the system's boundary, i.e., its smoke stack or translated from air pre-heater inlet data) must be unique relative to the supplied fuel and supplied combustion air.

[0205] The following elemental molar balances may be derived from the combustion equation, Eq.(19BL). The Γk expressions are simply convenient groupings of quantities, principally comprising measured effluents (known values) which have the greatest influence on the individual elements of interest. Minor fuel terms, αMAF-j, carried within Δk expressions are multiplied, initially, by an estimated fuel moles, xMAF. These minor terms are quickly resolved when converging on xMAF. Given these groupings, the Γk expressions of Eqs.(36) through (41), with solution of the moles of combustion oxygen (the term “a”) as discussed below, may be treated as known quantities. The elemental wet fuel components typically associated with black liquor fuels are considered unknowns, as are the fuel moles, these include the following: α1, α2, α3, Γ4, α5, α6, α10, α14, α15, α16 and x. Many recovery boilers use supplementary firing with gaseous fuel such as methane or fuel oil. Inclusion of such fuels is taught in 3 994 being described by the combustion equation term: xα0[CYRHZR] of that patent. Note that The Input/Loss Method accounts for all non-black liquor fuels, such as methane or fuel oil, but they are assumed to have known chemistries with known fuel flows. Mixed fuels are input to the system simulator (e.g., EX-FOSS) as a composite, combining them by using the FUEL program, although the resultant computed fuel chemistry, after proper weighting, is assigned to the black liquor by assuming the other fuels have fixed or predictable, and known, chemistries. Alternatively, such assignment may be made to a combined subset of fuels.

xα 1 =ΓN2−aφAct  (30)

x (α5+α2)/2=ΓH2O/2  (31) 6 $$ x ⁡ ( α 3 + α 2 / 2 ) = Γ O ⁢   ⁢ 2 + 3 ⁢ ( x MAF ⁢ Γ NA + σ N ⁢ k BL ) / 2 + x MAF ⁢ α MAF - 6 ⁡ ( σ N - 2 ⁢   ⁢ σ N ⁢ φ R + 1.0 ) ( 32 )  xα4=ΓCO2+xMAFΓNA+σNkBL  (33)

xα6=ΓSO2  (34)

xα 14 =σN(TNO+kBL)+σN(TNA+TNS+TNC/2)  (35)

[0206] where: 7$\begin{array}{cc}{\Gamma }_{\mathrm{N2}}=100-\&AutoLeftMatch;\left[D_{\mathrm{Act}}+G_{\mathrm{Act}}+e_{\mathrm{Act}}+m_{\mathrm{Act}}/2+{\sigma }_N{k}_{\mathrm{BL}}+\left(1.0-{\sigma }_N\right)x_{\mathrm{MAF}}{\alpha }_{\mathrm{MAF}-6}\right]-100{\phi }_{\mathrm{Act}}\left(R_{\mathrm{Act}}-1.0\right)/\left[R_{\mathrm{Act}}\left({\phi }_{\mathrm{Act}}+1.0\right)\right]& \left(36\right)\end{array}$ ΓH2O=JAct−bz−bTbA(1.0+β)  (37) 8$\begin{array}{cc}{\Gamma }_{\mathrm{O2}}=D_{\mathrm{Act}}+G_{\mathrm{Act}}+J_{\mathrm{Act}}/2-b_Z/2-b_T/2-\left(1.0+\beta \right)\left(a+b_A/2\right)+e_{\mathrm{Act}}/2+m_{\mathrm{Act}}/2+{\sigma }_N{k}_{\mathrm{BL}}\left(2{\phi }_R-1.0\right)& \left(38\right)\end{array}$ ΓCO2=DAct+eAct+V  (39)

ΓSO2=σN(TNA+TNS+kBL)+(1.0−σN)xMAFαMAF-6  (40)

ΓNA=σN(αMAF-14+αMAF-15−αMAF-16/2−αMAF-6)  (41)

[0207] In these relationships the subscript “Act” means an effluent measurement or assumption (an “actual” value). The term JAct in Eqs.(37) and (38) relating to the moles of effluent H2O could be input as a constant value or measured. All other values in Eqs.(36) through (41) are either evaluated explicitly based on input data, internal models and/or have minor import but are carried in the formulations for consistency. The amount of inert matter (xα10) associated with black liquor fuels as used in recovery boilers is typically small, less than 5% by weight and may be held constant which is the Preferred Embodiment, or as may be determined using the Air/Fuel ratio as taught in '994, or as otherwise obtained.

[0208] As a group, these relationships are of critical importance for understanding The Input/Loss Method as applied to recovery boilers. If fuel chemistry is to be resolved, thus heating value, boiler efficiency and thus accurate fuel flow, etc., then stoichiometric relationships generally representing Eqs.(30) to (41) must be resolved. These equations are not unique in their grouping of terms, further reductions and/or complexities are possible. The grouping of terms adopted here principally follows from the right-side of Eq.(19BL).

[0209] For the following discussion, assume initially that tube leakage, fuel inerts and fuel inorganics (bT, α10, α14, α15 and α16) have zero values. With this assumption, Eqs.(30) through (34) yield five equations with eight unknowns. For this situation, unknowns include α1 through α6, and the terms “a” and “x”. The term “x” is a convenience term and could be divided through changing the base of Eq.(19BL) to unity moles of fuel, thus eliminating use of xαj terms comprising two unknowns. However, if done, then the effluent's base becomes per mole of fuel, e.g., DAct/x, adding a different complexity involving the normalization of effluent measurements. Although the requirement ΣαMAF-j=1.00 is a convenience, it affords another, and viable, equation. By making a molar nitrogen balance, and assuming 100 moles of dry gaseous effluent at the boundary, the “a” quantity (moles of combustion oxygen) may be resolved independent of Eq.(30), thus reducing the unknowns; detailed below. Again, the entire combustion equation, Eq.(19BL), could be divided through by α4, or xα4, setting a carbon base. Effluent N2 could be resolved by difference assuming 100 moles of gaseous effluent (CO2, H2O, O2, SO2, the minor pollutants being measured or assumed), or N2 could be measured directly. Or, further still, by assuming constant values for fuel nitrogen and sulfur, α1 & α6 (typically minor fuel constituents), with resolution of “a”, and say: α3=1.0−ΣαMAF-j, j≠3, the system is reduced to three equations with four unknowns; these include Eqs.(31) through (33), with α2, α4, α5 and “x”. As another example, if 3, α14, α15 and α16 are assumed constant, then the combined Eqs.(31) and (32) (with cancellation of xα2) represents one equation with two unknowns, “x” and α5; however at least α2, α3 and α14 represent major and variable constituents of black liquor fuel. And, of course, further reductions and manipulations of unknowns and equations is entirely possible. However, close examination of the physical problem of combustion stoichiometrics, in which fuel chemistry is to be determined from Choice Operating Parameters, indicates that the mathematical system invokes a stiff matrix—and of course having further complexity if the black liquor's inorganic components are fully involved. In summary, these manipulations are discussed to emphasize that, as taught by this invention, algebraic manipulations must address the physical reality of the thermal system, this means the fuel's intrinsic chemical relationships. Such intrinsic chemical relationships must recognize for example that black liquor is a highly wet fuel consisting of both hydrocarbons and inorganics, and must also recognize the influences played by the fuel on the smelt measurements ΦR and ΦS by the principal gaseous effluent measurements DAct, JAct and GAct. Invoking such intrinsic chemical relationships reduces The Input/Loss Method's sensitivity to instrumentation errors and enhances the accuracy of the applied minimization techniques which correct such measurements. The solution problem is to employ intrinsic chemical relationships within the resolution of the combustion equation.

[0210] To address the solution problem, the Preferred Embodiment capitalizes on the intrinsic chemical relationships found in organic fuels between MAF hydrogen and MAF carbon, and, for black liquor fuel, relates stoichiometrically the organic to the inorganic compounds through the Sodium to Carbon ratio (ΦN). Further, the fuel mole term, x, is used in the Preferred Embodiment for this single variable appears in all stoichiometric conversions to mass flows (the consistent determination and conservation of which is an object of this invention); and represents an iterative parameter for Eq.(19BL). Further, the Preferred Embodiment does not require that the minor fuel constituents be assumed constant. Further, as will become apparent, the Preferred Embodiment allows use of multidimensional minimization techniques which addresses instrumentation errors.

[0211] Returning to the solution problem as posed by Eq.(19BL), if fuel chemistry is to be resolved then fundamental problems require solution—independent of algebraic manipulations. These problems include the following: 1) the black liquor's inorganic components; and 2) for any set of unknown fuel quantities, αj, there is always, at least, additional unknowns: in the combustion air term, “a”; in the fuel mole term, “x”; and in the unknown tube leakage (bT) a subject of this invention. The first of these problems is resolved through the Sodium to Carbon ratio ΦN, and obtaining correlations of molar MAF fractions of fuel potassium and chlorine as a function of sodium:

αMAF-15=A15+B15αMAF-14  (42)

αMAF-16=A16+B16αMAF-14  (43)

[0212] The second of these problems is solved, in part, by reducing αj quantities to a molar MAF bases, eliminating the influence of the two components not chemically involved in the hydrocarbon nor inorganic portions of black liquor's make-up, water and inerts, and then solving for the remaining αMAF-j terms. Dependent αMAF-j terms may be algebraically resolved by obtaining correlations of molar MAF fractions of fuel hydrogen, nitrogen and oxygen to molar MAF fuel carbon, Eqs.(48), (49) & (50), and molar MAF fuel potassium and chlorine to molar MAF fuel sodium, Eqs.(42) & (43); thus establishing intrinsic chemical relationships. Fuel carbon is then resolved explicitly by: converting the left-side of Eqs.(31), (32) & (33) to a MAF base; then subtracting Eq.(31) from (32) to eliminate the water term xMAFαMAF-2; then converting inorganic compounds associated with ΓNA, using Eqs.(42) & (43) and the ΦN ratio, to MAF carbon terms; substituting Eqs.(49) & (50); resulting in two functional equations, the combined Eqs.(31) & (32) and Eq.(33), with two unknowns, xMAF and αMAF-4; solving produces:

αMAF-4=[−ΓCO2(ξC7+ξC5)+ξC4ξC2]/[ΓCO2(ξC8+ξC6)+ξC4(1.0−ξC3)]  (44)

x MAF =ΓCO2/[αMAF-4−αMAF-4ξC3−ξC2]  (45)

[0213] However, before solving Eqs.(44) & (45), the “a” quantity (moles of combustion oxygen) used in Γ02 and bA, and thus in ξC4, ξC9 & ΓH2O, is resolved via Eq.(46) by recognizing: 1) the bases of Eq.(19BL) is that, at the boundary: 100 dry gaseous moles≡DAct+GAct+(h+βaφAct)+eAct+mAct+σNkBL+(1.0−σN)xMAFαMAF-6, and is used to substitute for “h” in a nitrogen molar balance simplifying with ΓN2; and 2) that: βa=100(RAct−1.0)/[RAct(1.0+φAct)].

a =(ΓN2−xMAFαMAF-1)/φAct  (46)

[0214] In these equations the following terms are developed from algebraic manipulations and simplifications. Furthermore, these equations example the process of establishing intrinsic chemical relationships as developed with the objective of establishing functionality between Choice Operating Parameters (e.g., effluent data) and fuel chemistry. 9$\begin{array}{cc}{\xi }_{\mathrm{C2}}={\sigma }_N\left[A_{15}-A_{16}/2-{\phi }_S\left({\alpha }_{\mathrm{MAF}-10}+A_{15}+A_{16}/2+v/x_{\mathrm{MAF}}\right)\right]& \left(47A\right)\\ {\xi }_{\mathrm{C3}}={\sigma }_N{\phi }_N\left[1.0+B_{15}-B_{16}/2-{\phi }_S\left(1.0+B_{15}+B_{16}/2\right)\right]& \left(47B\right)\\ {\xi }_{\mathrm{C4}}=-{\Gamma }_{\mathrm{H2O}}/2+{\Gamma }_{\mathrm{O2}}+2{\sigma }_N{k}_{\mathrm{BL}}\left(1.0-{\phi }_R\right)+\left(1.0-{\sigma }_N\right)x_{\mathrm{MAF}}{\alpha }_{\mathrm{MAF}-6}& \left(47C\right)\\ {\xi }_{\mathrm{C5}}={\sigma }_N\left[1.5\left(A_{15}-A_{16}/2\right)+\left(0.5-2{\phi }_R\right){\phi }_S\left({\alpha }_{\mathrm{MAF}-10}+A_{15}+A_{16}/2+v/x_{\mathrm{MAF}}\right)\right]& \left(47D\right)\\ {\xi }_{\mathrm{C6}}={\sigma }_N\left[1.5{\phi }_N\left(1.0+B_{15}-B_{16}/2\right)+\left(0.5-2{\phi }_R\right){\phi }_S{\phi }_N\left(1.0+B_{15}+B_{16}/2\right)\right]& \left(47E\right)\end{array}$ ξC7=−A3+A5  (47F)

ξC8=−B3+B5  (47G)

[0215] 10 $$ ξ C9 = Γ O2 + 3 ⁢ ( x MAF ⁢ Γ NA + σ N ⁢ k BL ) / 2 + x MAF ⁢ α MAF - 6 ⁡ ( σ N - 2 ⁢ σ N ⁢ φ R + 1.0 ) ( 47 ⁢ H )

[0216] Note that the v term of Eqs.(19BL) & (39), relating free carbon found in the smelt, may be assumed constant or may be expanded as a function of αMAF-4 or αMAF-14 if its magnitude effects the accuracy of the computed αMAF-4. Given αMAF-4, xMAF is then resolved via Eq.(45), based on Eq. (33). Given αMAF-4, sodium, αMAF-14, is then resolved from (ΦNαMAF-4). Note that a system carbon balance involves all major stoichiometric terms including xMAF, αMAF-4, DAct, RAct and β. In Eqs.(42) & (43) the constants Aj and Bj may be developed as part of the Reference Fuel Characteristics associated with a specific black liquor taken from the actual As-Fired, or obtained from generic specification of a black liquor fuel, or otherwise obtained. Such data may also be collected, or otherwise obtained, for fuel nitrogen, oxygen and hydrogen. Note that if ultimate analyses are not available for the fuel, for potential used in developing Aj and Bj, and especially A5 and B5 for MAF hydrogen, then determination of elemental (ultimate) analyses may be had from so-called proximate data. Such conversions from proximate to ultimate analyses is well established art, four such conversion methods are described in A. B. Gill, Power Plant Performance, Butterworths: London, 1984, Chapter 2, pages 70-77. And although such proximate to ultimate conversions are not the Preferred Embodiment, it affords another example of establishing intrinsic chemical relationships required of Reference Fuel Characteristics.

αMAF-1=A1+B1αMAF-4  (48)

αMAF-3=A3+B3αMAF-4  (49)

αMAF-5=A5+B5αMAF-4  (50)

[0217] Fuel MAF sulfur is resolved explicitly by solving Eq.(34) and applying Eqs.(42) & (43) and the ratios Φs and ΦN: 11$\begin{array}{cc}{\alpha }_{\mathrm{MAF}-6}={\Gamma }_{\mathrm{SO2}}/x_{\mathrm{MAF}}& \left(51\right)\\ ={\sigma }_N{\phi }_S\left[{\alpha }_{\mathrm{MAF}-10}+A_{15}+A_{16}/2+v/x_{\mathrm{MAF}}\right]+{\sigma }_N{\phi }_S{\phi }_N{\alpha }_{\mathrm{MAF}-4}\left[1.0+B_{15}+B_{16}/2+v/x_{\mathrm{MAF}}\right]+{\sigma }_N{k}_{\mathrm{BL}}/x_{\mathrm{MAF}}& \left(52\right)\end{array}$

[0218] where the inorganic terms comprising ΓSO2, and whose compounds describing the inorganics produced from black liquor combustion, see Eq.(19BL), are determined as follows:

σNTKO=σNxMAFαMAF-15  (53)

σNTNC=σNxMAFαMAF-16  (54)

σNTNS=σNΦR(xMAFαMAF-6−kBL)  (55)

σNTNA=σN(1.0−ΦR)(xMAFαMAF-6−kBL)  (56)

σN(TNO+kBL)=xMAFΓNA+σNkBL−σNxMAFαMAF-15  (57)

[0219] After Eqs.(44), (49) & (50) are resolved, fuel MAF water is resolved explicitly by employing Eq.(31), or by adding Eqs.(3 1) & (32) for numerical over-check, through substitution of xMAF: 12$\begin{array}{cc}{\alpha }_{\mathrm{MAF}-2}=\left({\Gamma }_{\mathrm{H20}}-x_{\mathrm{MAF}}{\alpha }_{\mathrm{MAF}-5}\right)/x_{\mathrm{MAF}}& \left(58A\right)\\ =\left(-{\alpha }_{\mathrm{MAF}-5}{\xi }_{\mathrm{C9}}+{\alpha }_{\mathrm{MAF}-3}{\Gamma }_{\mathrm{H2O}}\right)/\left({\xi }_{\mathrm{C9}}-{\Gamma }_{\mathrm{H2O}}/2\right)& \left(58B\right)\end{array}$

[0220] Fuel MAF nitrogen concentration is small, and typically may be fixed as a constant. Fuel oxygen may be described by Eq.(49). The MAF sulfur concentration may be computed assuming: αMAF-6=A6+B6αMAF-4, as opposed to Eq.(52), if the effluent SO2 measurement is questionable (or if limestone is injected and its conversion rate is questionable). As another alternative, either the oxygen or sulfur could be used to assure that ΣαMAF-j=1.0, where: j=1,3,4,5,6,14,15,16. The sole criteria in deciding the exact methodology is the reliability and availability of effluent data and its relative impact on MAF fuel terms. The Preferred Embodiment is to solve αMAF-3 using Eq.(49), αMAF-6 using Eq.(52), and αMAF-1 by balance. For example, although nitrogen is a major combustion effluent, using Eq.(30) to solve for fuel nitrogen, typically a very minor component of the fuel's makeup, would invite even slight errors made in determining effluent N2 (made either by direct measurement or determined by difference suggested in Eq.(36)) will greatly amplify the uncertainty in αMAF-1. Indeed, the ΓN2 term is used to resolve the combustion O2 term via Eq. (46), wherein the term xMAFαMAF-1 is resolved by balance and iteration on xMAF. Besides carbon, an important term is the fuel MAF hydrogen concentration which represents an intrinsic chemical relationship with carbon. Experience has shown that a valid relationship may be achieved since the hydrogen:carbon chemical bond is predominant with carbon:carbon. With an established MAF hydrogen to carbon relationship, the variability of a specific black liquor fuel then lies with its water, inorganics and inerts.

[0221] To this point in the process, resolution of fuel chemistry has been possible by assuming a Moisture-Ash-Free (MAF) basis. In carrying this out, it is not critical to the invention that carbon be computed as the independent quantity; nor is it critical to use a MAF basis as the inert fraction is typically small in black liquors, thus a dry basis may serve equally well. Hydrogen, or any other MAF fuel element could be considered as the independent variable. For example, if solving for hydrogen first, then αMAF-4=ƒ(αMAF-5). However, the Preferred Embodiment is to place independence with carbon, where of course the greatest accuracy and sensitivity may be found in the effluent measurements. CO2 does not exist in the combustion air to any appreciable concentration; it does not leak into the system; it is generated only from combustion and as a major effluent has obvious sensitivity to resolution of αMAF-4 of Eq.(44).

[0222] Thus all fuel constituents are therefore determined on a MAF bases. From these values, the wet base molar fuel fractions are then determined, as are the wet base moles of fuel (x) and the wet base (As-Fired) weight fractions (WFj) of all fuel components j:

αj=αMAF-j/(1.0+αMAF-2+αMAF-10)  (59)

x=x MAF (1.0+αMAF-2+αMAF-10)  (60)

xαj≡xMAFαMAF-j  (61)

WF j =αjNj/(ΣαjNj)  (62)

WF DRY-j =WWF j /(1.0 −WF2)  (63)

[0223] As taught in '879 the sensitivity of computed fuel chemistries to an effluent H2O measurement is substantial. When employing any of the Input/Loss methods which compute fuel chemistry, fuel chemistry will effect the computed fuel heating value, boiler efficiency and thus system efficiency. With such Input/Loss methods, computed fuel chemistries are generally sensitive to combustion gaseous effluents, including CO2, H2O and O2. The effluent H2O measurement (or its assumption), as well as effluent CO2 and O2, may be corrected through methods taught by '879; however, such measurements will be impacted by tube leaks, and tube leaks of even moderate flows. Both '879 and this invention discuss the optimization of certain parameters, such parameters include an assumed flow of water in-leakage into and mixing with the combustion gases; thus the modeling of tube leaks using a combination of combustion stoichiometrics, Eq.(19BL), etc., and optimization methods.

[0224] Having computed black liquor's fuel chemistry, the heating value of black liquor is next computed. Following the teachings of '994, black liquor heating value is determined base on a differential analysis. Note that the term NMAF is the molecular weight of the MAF-base fuel (without water and inerts). 13$\begin{array}{cc}\begin{array}{cc}\Delta {\mathrm{HHV}}_{\mathrm{MAF}-\mathrm{delta}}=& {\mathrm{HHV}}_{\mathrm{MAF}-\mathrm{Ref}}-(-178387.18{\alpha }_{\mathrm{MAF}-3}+\\ & 183591.92{\alpha }_{\mathrm{MAF}-4}+78143.68{\alpha }_{\mathrm{MAF}-5}-\\ & {\Delta H_{C-\mathrm{Cal}/\mathrm{BL}}^0)}_{\mathrm{Ref}}/N_{\mathrm{MAF}-\mathrm{Ref}}\end{array}& \left(64\right)\\ \begin{array}{cc}{\mathrm{HHV}}_{\mathrm{MAF}-\mathrm{uncorr}}=& (-178387.18{\alpha }_{\mathrm{MAF}-3}+183591.92{\alpha }_{\mathrm{MAF}-4}+\\ & 78143.68{\alpha }_{\mathrm{MAF}-5}-\\ & {\Delta H_{C-\mathrm{Cal}/\mathrm{BL}}^0)}_{\mathrm{Actual}}/N_{\mathrm{MAF}-\mathrm{Actual}}\end{array}& \left(65\right)\\ \begin{array}{cc}\Delta H_{C-\mathrm{Cal}/\mathrm{BL}}^0=& {\sigma }_N{\alpha }_{\mathrm{MAF}-6}[N_{\mathrm{Na2SO4}}\Delta H_{f-\mathrm{Cal}/\mathrm{Na2SO4}}^0+\\ & N_{\mathrm{H2S}}\Delta H_{f-\mathrm{Cal}/\mathrm{H2S}}^0-2N_{\mathrm{NaHS}}\Delta H_{f-\mathrm{Cal}/\mathrm{NaHS}}^0]\end{array}& \left(66\right)\end{array}$ HHVMAF=HHVMAF-uncorr+ΔHHVMAF-delta  (67)

HHV DRY =HHV MAF (1.0−WFDRY-10)  (68)

HHV AF =HHV DRY (1.0−WF2)  (69)

[0225] Eq.(66) illustrates, for example, that black liquor may contain only Na2CO3 and K2CO3 (which are assumed to appear as products),and NaHS which reacts to produce Na2SO4 and H2S, thus requiring a ΔHeat of Combustion correction, ΔHC-Cal/BL0. In Eq.(66) MAF sulfur is first taken as the reference value, αMAF-6-Ref, wherein ΔHC-Cal/BL0 is then used in Eq.(64); then it is taken as the actual, αMAF-6 via Eq.(52), wherein ΔHC-Cal/BL0 is then used in Eq.(65).

[0226] Boiler Efficiency for Recovery Boilers

[0227] Boiler efficiency computations are a vital portion of any of the Input/Loss methods, and especially so when resolving tube failure leakage flow. Methods presented in '429 describe how to compute a high accuracy boiler efficiency applicable for conventional steam generators burning common hydrocarbon fuels. However, as applied to recovery boilers, given black liquor fuel contains reactive inorganic compounds subjected to reduction, the methods of '429 require modification. Methods presented herein describe how to compute a high accuracy boiler efficiency applicable for black liquor fuels.

[0228] Notably The Input/Loss Method emphasizes computational consistency, consistency principally achieved through boiler efficiency and its computation of heating value as dependent on the calorimetric temperature, TCal. The calorimetric temperature must be consistently employed in all energy terms and Heats of Formation. For recovery boilers is it common industrial practice to correct the measured heating value for Heats of Formation associated with reactions particular to recovery boilers, and commonly the chemical reduction of Na2SO4. Such industrial corrections are taught in the following documents: “Performance Test Procedure Sodium Based Recovery Units”, CA Report No. 84041601, March 1996, TAPPI Press, Atlanta Ga. (see paragraphs 5.2.2.5, 7.2.2.8 and 7.2.2.9, and Appendix B); and Steam, Its Generation and Use, 40th Edition, Edited by S. C. Stultz and J. B. Kitto, published 1992 by The Babcock & Wilcox Company, Barberton, Ohio (see Chapter 26, page 26-5). Such common industrial corrections address the difference between ideal combustion products associated with a bomb calorimeter versus actual products associated with further reducing black liquor compounds. However, as taught herein such industrial corrections are thermodynamically inconsistent, and simply wrong. Industrial corrections reduce the higher heating value with a computed ΔHR term: HHVP−ΔHR; this term effecting the denominator in common boiler efficiency formulations. Given the teachings of '429 the ΔHR term is an intrinsic portion of the numerator of boiler efficiency; an affect of describing ideal versus actual products of combustion bound in the HPRAct and HRXAct terms. Further, the principles taught about the use and computation of L Factors, an important System Effect Parameter concept, are fundamentally based on the as—measured heating value combusted under theoretical conditions with dried fuel (not a HHVP−ΔHR concept); see the section on “System Effect Parameters” and its Eqs.(272) & (273) below. '429 teaches, for a HHV-based efficiency, that: 14$\begin{array}{cc}{\eta }_{B-\mathrm{HHV}}=\frac{\mathrm{BBTC}}{m_{\mathrm{AF}}\left(\mathrm{HHVP}+\mathrm{HBC}\right)}& \left(121\right)\\ =\frac{\mathrm{HHVP}+\mathrm{HBC}-\sum \mathrm{Losses}/m_{\mathrm{AF}}}{\mathrm{HHVP}+\mathrm{HBC}}& \left(122\right)\\ =\frac{\mathrm{HHVP}+\mathrm{HBC}-\left(\mathrm{HSL}+\mathrm{HNSL}\right)}{\mathrm{HHVP}+\mathrm{HBC}}& \left(123\right)\end{array}$

[0229] where the numerator of Eq.(123), if considering only HSL losses for clarity (i.e., ignoring Non-Chemistry & Sensible Heat Losses, ηA=1.0), may be expanded as:

=HHVP+HBC−HSL

=HPRAct-HHV+HRXAct-HHV

=−HPRAct-HHV+HRXCal-HHV+HBC

=−HPRAct-HHV+HPRIdea-HHV+HHVP+HBC

[0230] where this last expression illustrates that actual product compounds (e.g., as effected by the process' Reduction Efficiency) are net of ideal combustion products (found in a bomb calorimeter producing the measured HHVAF). Heats of Formation, and sensible heats associated with the Firing Correction term HBC, are all thermodynamically consistent. When reducing Na2SO4 in a recovery boiler via the reactions:

Na2SO4−>Na2S+2O2  (124)

CO2+Na2SO4−>SO2+Na2CO3+O2/2  (125)

[0231] the Heats of Formation of the non-ideal product compounds Na2S, SO2 and Na2CO3 are addressed through HPRAct and its Eq.(133), whereas the ideal—and the measurement base for all heating values—are addressed through HPRIdeal. Note that Na2CO3 as produced from reducing Na2SO4 is governed by the effluent measurement of SO2 (kBL), and although may be inaccurately measured, the measurement may be corrected or otherwise obtained. The quantity BBTC/(HHVP+HBC) of Eq.(121) must be constant for a given in-situ system burning a defined fuel; thus any changes in its combustion process (e.g., changes in Reduction Efficiency, ΦR) which will of course affect fuel flow, must also affect ηB (through HPRAct) in proportion, maintaining BBTC/(HHVP+HBC) constant. A (HHVP−ΔHR) procedure is inconsistent since the term: BBTC/(HHVP−ΔHR+HBC) would not remain constant since ΔHR=ƒ(ΦR); the term would produce a variable fuel flow.

[0232] Non-Chemistry and Sensible Heat Loss term (HNSL) is computed as taught in '429. With its determination, only the three major terms HPRAct, HRXAct & HBC remain to be defined to complete boiler efficiency. These are taught in the following paragraphs. To fully understand the formulations comprising HPRAct, HRXAct and HBC, take note of the subscripts associated with the individual terms. As example, when considering water product created from combustion, nComb-H2O of Eq.(13 1), its Heat of Formation (saturated liquid phase) at TCal must be corrected for boundary (Stack) conditions, thus, hstack−hf-Cal. The Enthalpies of Reactants of Eqs.(135) & (136) are determined from ideal products at Tcal, the Firing Correction then applied.

[0233] Differences in formulations required for higher or lower heating values should also be carefully reviewed, as recovery boilers are evaluated using either. Higher heating values require use of the saturated liquid enthalpy evaluated at TCal; lower heating values require the use of the saturated vapor at TCal. The quantities which are not so corrected are the last three terms in Eqs.(131) & (132): water born by air; known in-leakages (bZ); and tube leakage (bT). These terms undergo no transformations since they have non-fuel origins. Heating values and energies used in Eqs.(131) through (137) are always associated with the system boundary, specifically defined by: the entry point of the As-Fired fuel (or the “supplied” fuel in the case of fuel rejects); the ambient air conditions; the Stack location including gaseous effluent temperature and the Continuous Emission Monitoring System (typically measuring CO2).

[0234] Boiler efficiency is defined as either HHV- or LHV-based, being composed of a combustion efficiency and a boiler absorption efficiency:

ηB-HHV=ηC-HHVηA  (126)

ηB-LHV=ηC-LHVηA  (127)

[0235] The boiler absorption is determined from the Non-Chemistry & Sensible Heat Losses term (HBC) as fully discussed in '429. It is an important teaching herein that fuel flow compute identically from either efficiency base; thus, assuming no tube leakage: 15$\begin{array}{cc}{m}_{\mathrm{AF}}=\frac{\mathrm{BBTC}}{{\eta }_{B-\mathrm{HHV}}\left(\mathrm{HHVP}+\mathrm{HBC}\right)}=\frac{\mathrm{BBTC}}{{\eta }_{B-\mathrm{LHV}}\left(\mathrm{LHVP}+\mathrm{HBC}\right)}& \left(128\right)\end{array}$

[0236] Such computations of fuel flow using either efficiency, at a defined Tcal, is an important numerical overcheck of this invention.

[0237] Combustion efficiency is determined by the following, as either a HHV-based or a LHV-based efficiency: 16$\begin{array}{cc}{\eta }_{C-\mathrm{HHV}}=\frac{-{\mathrm{HPR}}_{\mathrm{Act}-\mathrm{HHV}}+{\mathrm{HRX}}_{\mathrm{Act}-\mathrm{HHV}}}{\mathrm{HHVP}+\mathrm{HBC}}& \left(129\right)\\ {\eta }_{C-\mathrm{LHV}}=\frac{-{\mathrm{HPR}}_{\mathrm{Act}-\mathrm{LHV}}+{\mathrm{HRX}}_{\mathrm{Act}-\mathrm{LHV}}}{\mathrm{LHVP}+\mathrm{HBC}}& \left(130\right)\end{array}$

[0238] The development of the combustion efficiency term, as computed based on HPRAct & HRXAct and involving systematic use of a combustion equation, such as Eq.(19BL), is believed an improved approach versus the primary use of individual “stack loss” terms. Mis-application of terms is greatly reduced, while numerical accuracy is increased. Most importantly, valid system mass and energy balances are assured.

[0239] The Enthalpy of Products (HPRAct) term is as follows. For higher heating value calculations: 17$\begin{array}{cc}{\mathrm{HPR}}_{\mathrm{Act}-\mathrm{HHV}}\equiv \sum {\mathrm{HPR}}_i+\left[n_{\mathrm{Comb}-\mathrm{H2O}}\left(\Delta H_{f-\mathrm{Cal}/\mathrm{liq}}^0+h_{\mathrm{Stack}}-h_{f-\mathrm{Cal}}\right)+n_{\mathrm{Fuel}-\mathrm{H2O}}\left(h_{\mathrm{Stack}}{h}_{f-\mathrm{Cal}}\right)+n_{\mathrm{CAir}-\mathrm{H2O}}\left(h_{\mathrm{Stack}}{h}_{g-\mathrm{Cal}}\right)+b_Z\left(h_{\mathrm{Stack}}{h}_{\mathrm{SteamZ}}\right)+b_T\left(h_{\mathrm{Stack}}{h}_{\mathrm{SteamT}}\right)\right]N_{\mathrm{H20}}/{\left({\mathrm{xN}}_{\mathrm{AF}}\right)}_{\mathrm{H20}}& \left(131\right)\end{array}$

[0240] For lower heating value calculations: 18$\begin{array}{cc}{\mathrm{HPR}}_{\mathrm{Act}-\mathrm{LHV}}\equiv \sum {\mathrm{HPR}}_i+\left[n_{\mathrm{Comb}-\mathrm{H2O}}\left(\Delta H_{f-\mathrm{Cal}/\mathrm{vap}}^0+h_{\mathrm{Stack}}-h_{g-\mathrm{Cal}}\right)+n_{\mathrm{Fuel}-\mathrm{H2O}}\left(h_{\mathrm{Stack}}{h}_{g-\mathrm{Cal}}\right)+n_{\mathrm{CAir}-\mathrm{H2O}}\left(h_{\mathrm{Stack}}{h}_{g-\mathrm{Cal}}\right)+b_Z\left(h_{\mathrm{Stack}}{h}_{\mathrm{SteamZ}}\right)+b_T\left(h_{\mathrm{Stack}}{h}_{\mathrm{SteamT}}\right)\right]N_{\mathrm{H20}}/{\left({\mathrm{xN}}_{\mathrm{AF}}\right)}_{\mathrm{H2O}}& \left(132\right)\end{array}$

[0241] where: 19$\begin{array}{cc}\begin{array}{cc}\begin{array}{c}{\mathrm{HPR}}_i=\mathrm{Enthalpy}\mathrm{of}\mathrm{non}­\mathrm{water}\mathrm{product}i\mathrm{at}\mathrm{the}\mathrm{boundary};\\ \mathrm{refer}\mathrm{to}\mathrm{the}\mathrm{right}\mathrm{side}\mathrm{of}\mathrm{Eq}.\left(19\mathrm{BL}\right);\mathrm{for}\\ \mathrm{example},\mathrm{if}\mathrm{product}i\mathrm{is}{\mathrm{Na}}_{2}S\mathrm{then}:\\ n_{\mathrm{Na2S}}={\sigma }_N{T}_{\mathrm{NS}}={\sigma }_N{\phi }_R\left(x_{\mathrm{MAF}}{\alpha }_{\mathrm{MAF}-6}-k_{\mathrm{BL}}\right).\end{array}& \end{array}& \\ {\mathrm{HPR}}_i\equiv \left[\Delta H_{f-\mathrm{Cal}/i}^0+{\int }_{T_{\mathrm{Cal}}}^{T_{\mathrm{Stack}}}d{h}_i\right]n_i{N}_i/\left({\mathrm{xN}}_{\mathrm{AF}}\right)& \left(133\right)\end{array}$

[0242] nComb-H2O=Molar water found at the boundary formed directly from combustion of a hydrocarbon (xα5). 20$\begin{array}{cc}{n}_{\mathrm{Fuel}-\mathrm{H2O}}& =\mathrm{Molar}\mathrm{water}\mathrm{found}\mathrm{at}\mathrm{the}\mathrm{boundary}\mathrm{born}\mathrm{by}\\ & \mathrm{As­Fired}\mathrm{fuel}\left(\mathrm{as}\mathrm{total}\mathrm{inherent}\mathrm{and}\mathrm{surface}\mathrm{moisture}\right).\\ & \equiv J_{\mathrm{Act}}-\left[b_A\left(1.0+\beta \right)+b_Z+b_T+x{\alpha }_5\right]=x{\alpha }_2\\ n_{\mathrm{CAir}-\mathrm{H2O}}& =\mathrm{Molar}\mathrm{water}\mathrm{found}\mathrm{at}\mathrm{the}\mathrm{boundary}\mathrm{born}\mathrm{by}\\ & \mathrm{combustion}\mathrm{air}\mathrm{and}\mathrm{air}\mathrm{in­leakage}\\ & \equiv b_A\left(1.0+\beta \right)\end{array}$

[0243] hStack-H2O=ƒ(PStack-H2O, TStack), where PStack-H2O is water's partial pressure per total wet molar and is given by:

P Stack-H2O =P Amb (JAct+βbA)/(100+JAct+βbA).

[0244] hSteamZ=Weighted specific enthalpy of bZ in-leakage.

[0245] hSteamT=Specific enthalpy of tube leakage (bT).

[0246] ΔHƒ-Cal/i0=Heat of Formation of compound i evaluated at the calorimetric temperature, TCal. 21$\Delta H_{f-\mathrm{Cal}/i}^0=\Delta H_{f-77/i}^0+{\int }_{77}^{T-\mathrm{Cal}}d{h}_{\mathrm{Compound}/i}-\sum {\int }_{77}^{T-\mathrm{Cal}}d{h}_{\mathrm{Elements}}$

[0247] The Enthalpy of Reactants (HRXAct) term is as follows. Note that although SO2 is a common gaseous product of ideal combustion of conventional fossil fuels, it is not assumed to be produced from the ideal combustion of black liquor. For higher heating value calculations: 22$\begin{array}{cc}{\mathrm{HRX}}_{\mathrm{Act}-\mathrm{HHV}}\equiv \mathrm{HHVP}+\mathrm{HBC}+{\mathrm{HPR}}_{\mathrm{CO2}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{H2O}-\mathrm{Ideal}-\mathrm{HHV}}+{\mathrm{HPR}}_{\mathrm{Na2CO3}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{Na2SO4}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{NaCl}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{K2CO3}-\mathrm{Ideal}}& \left(135\right)\end{array}$

[0248] For lower heating value calculations: 23$\begin{array}{cc}{\mathrm{HRX}}_{\mathrm{Act}-\mathrm{LHV}}\equiv \mathrm{LHVP}+\mathrm{HBC}+{\mathrm{HPR}}_{\mathrm{CO2}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{H2O}-\mathrm{Ideal}-\mathrm{LHV}}+{\mathrm{HPR}}_{\mathrm{Na2CO3}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{Na2SO4}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{NaCl}-\mathrm{Ideal}}+{\mathrm{HPR}}_{\mathrm{K2CO3}-\mathrm{Ideal}}\begin{array}{cc}\mathrm{where}:& \\ {\mathrm{HPR}}_{\mathrm{CO2}-\mathrm{Ideal}}& =\mathrm{Energy}\mathrm{of}{\mathrm{CO}}_2\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion}\mathrm{at}\mathrm{the}\mathrm{calibration}\mathrm{temperature}.\\ & \equiv \Delta H_{f-\mathrm{Cal}/\mathrm{CO2}}^0{\alpha }_4{N}_{\mathrm{CO2}}/N_{\mathrm{AF}}\\ {\mathrm{HPR}}_{\mathrm{H2O}-\mathrm{Ideal}-\mathrm{HHV}}& \mathrm{Energy}\mathrm{of}{H}_{2}O\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion},\mathrm{ending}\mathrm{with}\mathrm{condensed}\mathrm{water},\mathrm{at}\mathrm{the}\\ & \mathrm{calibration}\mathrm{temperature}.\\ & \equiv {\left(\Delta H_{f-\mathrm{Cal}/\mathrm{liq}}^{0}N\right)}_{\mathrm{H2O}}{\alpha }_5/N_{\mathrm{AF}}\\ {\mathrm{HPR}}_{\mathrm{H2O}-\mathrm{Ideal}-\mathrm{LHV}}& \mathrm{Energy}\mathrm{of}{H}_{2}O\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion},\mathrm{ending}\mathrm{with}\mathrm{water}\mathrm{vapor},\mathrm{at}\mathrm{the}\\ & \mathrm{calibration}\mathrm{temperature}.\\ & \equiv {\left(\Delta H_{f-\mathrm{Cal}/\mathrm{vap}}^{0}N\right)}_{\mathrm{H2O}}{\alpha }_5/N_{\mathrm{AF}}\\ {\mathrm{HPR}}_{\mathrm{Na2CO3}-\mathrm{Ideal}}& \mathrm{Energy}\mathrm{of}{\mathrm{Na}}_2{\mathrm{CO}}_3\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion}\mathrm{at}\mathrm{the}\mathrm{calibration}\mathrm{temperature}.\\ & \equiv \Delta H_{f-\mathrm{Cal}/\mathrm{Na2CO3}}^0{\alpha }_N{T}_{\mathrm{NO}}{N}_{\mathrm{Na2C03}}/\left({\mathrm{xN}}_{\mathrm{AF}}\right)\\ {\mathrm{HPR}}_{\mathrm{Na2SO4}-\mathrm{Ideal}}& \mathrm{Energy}\mathrm{of}{\mathrm{Na}}_2{\mathrm{SO}}_4\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion}\mathrm{at}\mathrm{the}\mathrm{calibration}\mathrm{temperature}.\\ & \equiv \Delta H_{f-\mathrm{Cal}/\mathrm{Na2SO4}}^0{\alpha }_N{T}_{\mathrm{NA}}{N}_{\mathrm{Na2S04}}/\left({\mathrm{xN}}_{\mathrm{AF}}\right)\\ {\mathrm{HPR}}_{\mathrm{NaCl}-\mathrm{Ideal}}& \mathrm{Energy}\mathrm{of}\mathrm{NaCl}\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion}\mathrm{at}\mathrm{the}\mathrm{calibration}\mathrm{temperature}.\\ & \equiv \Delta H_{f-\mathrm{Cal}/\mathrm{NaCl}}^0{\alpha }_N{T}_{\mathrm{NC}}{N}_{\mathrm{NaCl}}/\left({\mathrm{xN}}_{\mathrm{AF}}\right)\\ {\mathrm{HPR}}_{\mathrm{K2CO3}-\mathrm{Ideal}}& \mathrm{Energy}\mathrm{of}{K}_2{\mathrm{CO}}_3\mathrm{ideal}\mathrm{product}\mathrm{from}\mathrm{complete}\\ & \mathrm{combustion}\mathrm{at}\mathrm{the}\mathrm{calibration}\mathrm{temperature}.\\ & \equiv \Delta H_{f-\mathrm{Cal}/\mathrm{K2CO3}}^0{\alpha }_N{T}_{\mathrm{KO}}{N}_{\mathrm{K2CO3}}/\left({\mathrm{xN}}_{\mathrm{AF}}\right)\end{array}& \left(136\right)\end{array}$

[0249] It should be noted that in Paragraphs 0069 and 0070 when using the fuel's calorimetric temperature, established when determining the fuel's heating value, as the thermodynamic reference energy level for the Enthalpy of Products and for the Enthalpy of Reactants, that these include corrected Heats of Formation, ΔHƒ-Cal/i0 for example as taught by Eq.(134), said ΔHƒ-Cal/i0 term being used in Eqs.(131), (132), (133), (135) and (136). Further, that ΔHƒ-Cal/i0 is corrected relative to a standard Heat of Formation for substance i, taken at 77 F (25 C) and denoted by ΔHƒ-77/i0.

[0250] The Firing Correction (HBC) term is given by the following. 24$\begin{array}{cc}\mathrm{HBC}\equiv {C_P\left(T_{\mathrm{AF}}-T_{\mathrm{Cal}}\right)}_{\mathrm{Fuel}}+\left(Q_{\mathrm{SAH}}+W_{\mathrm{FD}}\right)/m_{\mathrm{AF}}+\left[{\left(h_{\mathrm{Amb}}-h_{\mathrm{Cal}}\right)}_{\mathrm{Air}}a\left(1.0+\beta \right)\left(1.0+{\phi }_{\mathrm{Act}}\right)N_{\mathrm{Air}}+{\left(h_{g-\mathrm{Amb}}-h_{g-\mathrm{Cal}}\right)}_{\mathrm{H2O}}{b}_A\left(1.0+\beta \right)N_{\mathrm{H2O}}+{\left(h_{\mathrm{SteamZ}}-h_{f-\mathrm{Cal}}\right)}_{\mathrm{H2O}}{b}_Z{N}_{\mathrm{H2O}}+{\left(h_{\mathrm{SteamT}}-h_{f-\mathrm{Cal}}\right)}_{\mathrm{H2O}}{b}_T{N}_{\mathrm{H2O}}\right]/\left({\mathrm{xN}}_{\mathrm{AF}}\right)& \left(137\right)\end{array}$

[0251] where:

[0252] hg-Amb-H2O=Saturated water vapor enthalpy at ambient dry bulb, TAmb.

[0253] CP(TAF−TCal)Fuel=Sensible heat in the As-Fired fuel relative to TCal.

[0254] (hAmb−hCal)Air=ΔEnthalpy of combustion dry air relative to TCal.

[0255] (hg-Amb−hg-Cal)H2O=ΔEnthalpy of moisture in combustion air relative to saturated water vapor at TCal.

[0256] (hSteamZ−hf-Cal)H2O=ΔEnthalpy of bZ known water in-leakages (at an average hSteamZ) relative to saturated liquid water at TCal.

[0257] (hSteamZ−hf-Cal)H2O=ΔEnthalpy of tube leakage (at hSteamT) relative to saturated liquid water at TCal.

[0258] The above equations are dependent on common system parameters. Common system parameters are defined following their respective equations, Eqs.(131) through (137). The BBTC term, comprising common system parameters, is determined from commonly measured or determined working fluid mass flow rates, pressures and temperatures (or qualities). Further, supporting terms such as thermodynamic properties, radiation & convection loss curves, guidelines for estimating Non-Chemistry & Sensible Heat Losses, sodium compound effects, miscellaneous terms, etc. are discussed in the following codes, standards and patents: the American Society of Mechanical Engineers' (ASME) Performance Test Codes (PTC) 4.1 and 4; the “Performance Test Procedure Sodium Based Recovery Units”, CA Report No. 84041601, March 1996, TAPPI Press, Atlanta Ga.; the German standard “Acceptance Testing of Steam Generators, DIN 1942, DIN DEUTSCHES Institut Fur Normung E.V., February 1994; the European standard (draft) prEN 12952-15:1999 (also: CEN/TC 269/WG 3 N 337), “Water-Tube Boilers and Auxiliary Installations—Part 15: Acceptance Tests”, November 1999, European Committee for Standardization, Central Secretariat, rue de Stassart, 36, Brussels; the British Standard “Code for Acceptance Tests on Stationary Steam Generators of the Power Station Type”, BS 2885:1974, ISBN: 0 580 08136 2; and throughout U.S. Pat. No. '420, '470, '994 and '429. If conflicts arise between any method or procedure of these codes, standards or patents, this disclosure's methods, procedures and invention shall prevail.

[0259] It is interesting that the German standard DIN 1942 states that its reference temperature (tb) for boiler efficiency is 25 C (77 F). However, in its paragraph 6.2, DIN 1942 allows that “other temperatures may be agreed upon” specifying in its Eq.(1b) how to correct heating value for a (tb25) effect. DIN 1942's Eq.(1b) not only corrects heating value using sensible heats (CPΔT), incorporating all terms described by the Firing Correction (HBC) taught herein, but also corrects heating value for flue gas sensible heat and thus the denominator of ηB. This flue gas correction is incorrect. As taught herein a flue gas sensible heat term appears in Eq.(133), integrating the (dh) energy of combustion products from TCal, to TStack effecting the term (−HPRAct-LHV+HRXAct-LHV), and thus only the numerator of ηB. Further, DIN 1942 does not teach to correct the ΔHƒ-Cal/i0 term of Eq.(133) for a reference temperature different from TCal, or tb; the term is ignored in DIN 1942 but is taught herein through Eq.(134). DIN 1942 does not correct Heats of Formation to achieve consistency with calorimetric temperature as taught herein; indeed, nowhere in DIN 1942 does any Heat of Formation appear.

[0260] The draft European standard mimics DIN 1942; it also discusses “Input/Loss Calculations” in its Section 9.4.5.2, but these calculations are not related to “Input/Loss methods” as used herein as defined in paragraph 0035. The British Standard requires that “the standard reference temperature for the calorific value of solid, liquid and gaseous fuels shall be 25 C”(Clause 2), and offers no means to alter. The TAPPI performance test standard references energy flows to a constant 77 F, but offers no means to alter. TAPPI's Appendix B develops a “Heat of Reaction Correction”, correcting heating values as (HHVP−ΔHR), however no TCal correction or alternation to ΔHƒ0 is assumed in this development nor could one be inferred as it is not mentioned. PTC 4.1 references energy flows to an arbitrary “reference air temperature”, TRA. PTC 4 references energy flows to a constant 77 F, but offers no means to alter. Neither DIN 1942, the draft European standard, the British standard, the TAPPI standard, PTC 4.1 nor PTC 4 mention how a reference temperature should be established, other than setting it constant. There is no mention of the fuel's calorimetric temperature in any of these standards, nor the correction of reference Heats of Formation to the calorimetric temperature, nor using the fuel's calorimetric temperature as the reference temperature. Note that Heats of Formation are normally cited in the literature as being referenced to 25 C and 1.0 bar pressure. The Preferred Embodiment of '429 and for this disclosure requires the reference temperature to be taken at the calorimetric temperature established when determining heating value, and that all energy terms, including those associated with Heats of Formation, HPRAct, HRXAct and HBC, be consistently referenced to this same calorimetric temperature.

[0261] Several PTCs and “coal” textbooks employ simplifying assumptions regarding the conversion of heating values. For example, a constant is sometimes used to convert from a constant volume process HHVAF (i.e., bomb calorimeter), to a constant pressure process HHVP. The following is preferred for completeness, for solid and liquid fuels:

HHVP=HHV AF +ΔH V/P   (138)

ΔH V/P ≡R Gas T Cal,Abs (α5/2−α1)/(ξREJJMechNAF)  (139)

[0262] where, in US Customary Units: TCal,Abs is absolute calorimetric temperature (deg-R); RGas≡1545.325 ft-lbf/mole-R; and JMech≡778.169 ft-lbf/Btu. For gaseous fuels, the only needed correction is the compressibility factor (Z) assuming ideally computed heating values at the specified TCal:

HHVP=HHVIdealZ  (140)

[0263] Z and HHVIdeal may be evaluated using American Gas Association procedures. To convert from a higher heating value (also termed gross or upper) to a lower heating value (also termed net) use of Eq.(142) is exact, where Δhfg-Cal/H2O is evaluated at TCal. The oxygen in the effluent water is assumed to derive from combustion air and not from fuel oxygen: thus α3 is not included as there is no molar change, fuel oxygen is taken as the diatomic in Eq.(19BL).

LHV AF =HHV AF −ΔH L/H   (141)

ΔHL/H=Δhfg-Cal/H2O(α2+α5)NH2O/(ξREJNAF)  (142)

[0264] Within Eq.(142) the ξREJ term accounting for rejected fuel is defined as: ξREJ≡(1.0−WF′Ash-AF)/(1.0−WFAsh-Sup). ξREJ also corrects Eq.(144). These same procedures are applicable for a fuel cleaning process where the fuel's mineral matter (inert) is removed.

LHVP=LHV AF +ΔH V/P −ΔH corr-LHV   (143)

ΔH corr-LHV =ΔH L/H (ξREJ−1.0)/ξREJ  (144)

[0265] Tube Failure Detection Methods

[0266] This invention teaches to add, apart from bZ, a term descriptive of tube leakage to the typical combustion equation; its symbol herein is bT whose units are moles of liquid water or steam in-leakage per 100 moles of dry gaseous effluent. A procedure is taught by which the flow of in-leakage associated with a tube failure and its location may be determined. This procedure relies on a hydrogen stoichiometric balance incorporating bT, applying limit tests on fuel constituents and tests on the rate of change of certain parameters; which tests provide an indication of possible tube failure and its stoichiometric causality (termed a “Tube Failure Mechanism”, indicated by a numerical identification). Once a Tube Failure Mechanism has been identified, a tube failure flow rate is determined, but determined using a separative analysis technique. In summary, the Preferred Embodiment relies on modification of a traditional combustion equation, employing the bT term, it relies on limits and rates testing, and it relies on successive optimizations of Choice Operating Parameters which minimize errors in System Effect Parameters. All of these topics are discussed in detail in this disclosure. Forming a hydrogen stoichiometric balance using Eq.(19BL), and solving for bT, results in Eq.(200) and thus allows the determination of tube leakage in moles. Note that the effluent moisture (at the Stack), is defined as: JAct ≡j+βbA. Also note that Eq.(19BL) represents a mathematical model of the combustion process using a molar base. The manipulation of such a combustion equation when employing traditional fossil fuels (non-black liquor), is taught in '994; whereas methods required for black liquor fuels, considered novel, are as taught herein.

b T =J Act −x (α2+α5)−bZ−bA(1.0+β)  (200)

[0267] Eq.(200) illustrates that for bT to be positive, i.e., a tube leak being detected, that unique balance must be developed between the assumed (or measured) effluent water (JAct) and the predominating negative terms: combustion water (xα2+xα5), bZ, and moisture in the combustion air and in the air leakage bA(1.0+β). Eq.(200) demonstrates that use of an effluent H2O instrument, measuring JAct, may not detect tube failures. For example, any unusual increase in JAct could be caused by off-setting effects from high fuel water, high moisture in the combustion air, high air pre-heater leakage (a high β) and/or periodic soot blowing flow and/or use of atomizing steam (bZ). Further, a tube leak could exist when the JAct term is decreasing as caused, for example, by a large decrease in fuel water (when, at the same time, bT is increasing). To resolve such difficulties, this invention teaches the use of Eq.(200) in conjunction with one of the Input/Loss methods in which fuel chemistry, the αj terms, are determined.

[0268] When Input/Loss methods compute fuel chemistry associated with black liquor, such chemistry will include at least the determination of fuel elemental carbon (α4), fuel elemental hydrogen (α5), fuel water (α2), and fuel elemental sodium (α14). Typical Input/Loss methods will determine such quantities, in part, based on Operating Parameters including principal effluent concentrations (CO2, O2 and H2O ), combustion air psychrometrics (leading to bA), and any known water and steam flows such as soot blowing and atomizing of fuel (bZ). Further, The Input/Loss Method teaches when applied to recovery boilers, as a portion of Reference Fuel Characteristics, is to obtain intrinsic chemical relationships between fuel MAF hydrogen and carbon, and fuel sodium and carbon, for example for hydrogen: αMAF-5=A5+B5αMAF-4; αMAF-5=A5+B5αMAF-4+C5(αMAF-4)2; or in general: αMAF-5=ƒ(αMAF-4); and for example for sodium: αMAF-14=ΦNαMAF-4; or in general: αMAF-14=ƒ(αMAF-4) Such functional relationships have been found to be most useful for fossil fuels including black liquor fuels. It has been found that the constants A5, B5, C5 and ΦN are typically remarkably constant for a specific fuel. Such correlations then establish inter-dependency between fuel carbon, fuel sodium, and all principal effluents, and thus, through Eq.(200) and resolution of Eq.(19BL), between the effluent concentrations CO2, O2 and H2O , and the important bT term. In this context resolution of inorganic stoichiometrics (i.e., sodium) is vitally important to resolution of the organics, and thus to resolution of the tube failure (bT) term. Given such inter-dependencies, it is most likely that when assuming bT=0.0, when in fact a tube is leaking, one or more fuel molar quantities (α1, α2, . . . α14, α15, α16) will compute outside reasonability limits or even as negative values. Experience applying this invention has taught that fuel water will commonly compute as negative, even with a moderate leak when initially assuming bT=0.0 in Eq.(200).

[0269] In like manner, and especially for small leaks (when assuming bT=0.0), the fuel carbon, hydrogen and sodium terms could exceed reasonability limits; where for example, assuming that the constant B5 is negative and C5 is zero (which is typical for black liquor fuels): αMAF-5<αMAF-5/min and/or αMAF-4>αMAF-4/max. Such behavior when using Eq.(200), when first assuming bT=0.0 and then evaluating for reasonability limits, leads directly to an indication of tube leakage. This process is termed a possible Tube Failure Mechanism, that is an indication of possible tube leakage has been found by applying stoichiometric considerations (min/max checks); however further processing is called for to determine its validity and, if a valid leak, then to determine its mass flow rate and the location of the leak.

[0270] To fully expand the concept of Tube Failure Mechanisms, TABLE 1A and TABLE 1B teach 35 such mechanisms, identified by a number, of how a tube failure may be detected through stoichiometrics, knowledge of such detection mechanism being important to the system operator. TABLE 1A presents static trip mechanisms, while TABLE 1B presents dynamic trip mechanisms related to rates of change. Experience in demonstrating this invention at two large power plants has indicated that making assumptions as to “apparently” impossible Tube Failure Mechanisms is not advised. Thus both minimum and maximum trip mechanisms are demonstrated in TABLE 1A, as all are tested when monitoring a thermal system on-line when invoking the methods of this invention (i.e., invoking The Input/Loss Method's “Tube Failure Model”). For example, a cursory evaluation would suggest that a low fuel water concentration (αMAF-2 or WFH2O) could not indicate a tube failure given the mechanics of Eq.(200). However, if the thermal system experiences a small but steadily increasing tube leakage The Input/Loss Method of '879 could steadily correct effluent water concentration downwards, causing failure mechanism #41 or #51; or water correction factors might exceed a lower bound causing failure mechanism ID #71. But also, unplanned scenarios of how The Input/Loss Method is correcting effluent water and other Choice Operating Parameters could create unexpected Tube Failure Mechanisms through complex stoichiometric relationships. Such considerations thus call for a blanket examination of all trip mechanisms, which is the Preferred Embodiment. TABLE 1A encompasses the most likely of stoichiometric mechanisms based on the experiences gained demonstrating this invention, while TABLE 1B encompasses likely rates of change of relevant parameters. Note that several of the rate mechanisms indicated in TABLE 1B were found by testing the methods of this invention at a power plant; specifically Tube Failure Mechanism #64, #121 and #131 were found unusually sensitive. However, in extending the teachings of this invention as taught through TABLE 1A and TABLE 1B, there are other Tube Failure Mechanisms which may become apparent without departing from the scope and spirit of the present invention.

[0271] ///

[0272] ///

[0273] ///

[0274] ///

[0275] ///

TABLE 1A
Static Tube Failure Mechanisms
Mech.
ID	Trip Mechanism	Comments
11	JAct < JAct/min	Effluent H2O concentration at Stack; min. JAct is not likely.
12	JAct > JAct/max	Effluent H2O concentration at Stack.
21	αMAF-4 < αMAF-4/min	MAF molar fraction of carbon; min. αMAF-4 is not likely.
22	αMAF-4 > αMAF-4/max	MAF molar fraction of fuel carbon.
23	Negative square root.	Resolution of MAF carbon (αMAF-4) may require solving a
		second order equation, thus the possibility of tripping on a
		negative square root (or any mathematical indeterminate).
31	αMAF-5 < αMAF-5/min	MAF molar fraction of fuel hydrogen.
32	αMAF-5 > αMAF-5/max	MAF molar fraction of hydrogen; max. αMAF-5 is not likely.
41	αMAF-2 < αMAF-2/min	MAF molar fraction of fuel water.
42	αMAF-2 > αMAF-2/max	MAF molar fraction of fuel water.
51	WFH2O < WFH2O/min	As-Fired weight fraction of fuel water.
52	WFH2O > WFH2O/max	As-Fired weight fraction of fuel water.
61	αMAF-1 < αMAF-1/min	MAF molar fraction of nitrogen in fuel; nitrogen is
		computed by balance: αMAF-1 = 1.0 − ΣαMAF-k,k≠1,2,10.
62	αMAF-1 > αMAF-1/max	MAF molar fraction of nitrogen in fuel; nitrogen is
		computed by balance: αMAF-1 = 1.0 − ΣαMAF-k,k≠1,2,10.
71	C2S < C2S/min	Correction factor for effluent H2O.
72	C2S > C2S/max	Correction factor for effluent H2O.
81	C1S < C1S/min	Correction factor for effluent CO2.
82	C1S < C1S/max	Correction factor for effluent CO2.
91	αMAF-14 < αMAF-14/min	MAF molar fraction of sodium in fuel.
92	αMAF-14 < αMAF-14/max	MAF molar fraction of sodium in fuel.
99	mT > 0.0	Tube leakage found in previous monitoring cycle.

[0276] ///

[0277] ///

[0278] ///

[0279] ///

TABLE 1B
Dynamic Tube Failure Mechanisms
Mech.
ID	Trip Mechanism	Comments
14	d(JAct)/dt > Limit	Rate of change in uncorrected effluent H2O.
15	d(C2SJAct)/dt > Limit	Rate of change in corrected effluent H2O.
24	d(DAct)/dt > Limit	Rate of change in uncorrected effluent CO2.
25	d(C1SDAct)/dt > Limit	Rate of change in corrected effluent CO2.
64	d(GAct)/dt > Limit	Rate of change in uncorrected effluent O2.
66	d(C7SGAct)/dt > Limit	Rate of change in corrected effluent O2.
67	d(C7S)/dt > Limit	Rate of change in the correction to effluent O2.
74	d(C2S)/dt > Limit	Rate of change in the correction to effluent H2O.
84	d(C1S)/dt > Limit	Rate of change in the correction to effluent CO2.
101	d(L′Fuel)/dt > Limit	Rate of change in the computed L-Factor.
103	d(HHVAF)/dt > Limit	Rate of change in the computed heating value.
109	d(mAF)/dt > Limit	Rate of change in the computed fuel flow.
111	d(WFH2O)/dt > Limit	Rate of change in the computed fuel water fraction.
121	d[C2SJAct/(C7SGAct)]/dt > Limit	Rate of change in the parameter of corrected effluent
		H2O divided by the corrected effluent O2.
131	d[mAF/HHVAF]/dt > Limit	Rate of change in the parameter of computed fuel flow
		divided by computed fuel heating value.

[0280] It is an important aspect of the present invention that it may be integrally involved with any of the Input/Loss methods which compute fuel chemistry. As seen with the used of TABLE 1A, without a determination of fuel chemistry, i.e., computing αi quantities based on consistent stoichiometrics (with or without tube leakage), then the use of Eq.(200) as taught herein to detect tube leaks would become limited. More specifically, this invention is integrally involved with The Input/Loss Method of '994 and '879 as the determination of fuel chemistry is then based on a selection of Choice Operating Parameters which might well require corrections for stoichiometric consistencies. For example, if effluent water was being measured, but whose signal was not corrected for stoichiometric consistency as taught in '879, resolution of tube leaks (even with computed αi; quantities) would be hampered; especially so if fuel chemistry was assumed constant.

[0281] Minimization Techniques, Background

[0282] The following four paragraphs teach the use of multidimensional minimization techniques, including the use of System Effect Parameters and Choice Operating Parameters. Taught is how their use allows resolution of flow rates associated with tube failures. Although only the Simulated Annealing technique is the Preferred Embodiment for determining tube failure flow rate, discussion of all techniques, including neural networks, is made to present the full scope and spirit of this invention.

[0283] Any number of multidimensional minimization techniques may be used by this invention; to provide full teaching four are discussed in detail. All techniques seek to minimize the numerical value of an objective function. The techniques discussed include: Broyden-Fletcher-Goldfarb-Shanno (BFGS), generic Conjugate Gradient, Newton-Raphson and Simulated Annealing techniques; references cited below. These techniques, and, notably, their combinations, are designed to address all situations of bias in Choice Operating Parameters. All of these techniques, except Simulated Annealing, employ derivatives of the objective function with respect to the independent variable. These techniques all require input of initial values of Choice Operating Parameters (Λ0−i). The BFGS, generic Conjugate Gradient and Newton-Raphson techniques employ unconstrained searches towards optima. Simulated Annealing employs a random but constrained search through which the Choice Operating Parameters are numerically bounded by lower and upper limits. From research and study conducted to develop this invention, the objective functions described below have proven to be superior for a wide variety of thermal systems burning any fossil fuels.

[0284] A common problem facing minimization techniques is the so-called shallow valley problem in which an appreciable change in an independent variable has a small effect on the objective function, even through that change is both real and appropriate to the physical system. This is especially true when applied to the determination of tube leakage flow rate in which a single and small tube failure of, say, 5000 lbm/hr in a large steam generator might represent ≈0.1% or less, of its feedwater flow. Study conducted for the development of this invention, and considered unique to it, has found that the Bessel function of the first kind is ideally suited to diminish the impact of the shallow valley problem. The Bessel function emulates the sensitivity that important Choice Operating Parameters have on both System Effect Parameters and on the descriptive thermal system in general. The Bessel function of the first kind of order zero (J0) has a relatively flat (shallow) functionality as its argument approaches zero. Apart from this situation, the function offers non-linearity which is advantageous in converging out-lying arguments. Of great importance to those techniques employing derivatives, is that the derivative of J0 is a Bessel function of the first kind of order one (J1), having a high degree of sensitivity as its argument approaches zero. This derivative relationship addresses a significant number of shallow valley problems presented by Choice Operating Parameters associated with fossil-fired systems. Another technique addressing the shallow valley problem and involving use of the Bessel function is the formulation of its argument, termed either λL, λW, λH or λT [i.e., J0(λL), J0(λW), J0(λH) or J0(λT)]; these arguments are fully discussed below, being defined by Eqs.(202A), (202B), (202C) and (202D).

[0285] The objective function, F, is a function of independent variables {right arrow over (x)}, i.e., F({right arrow over (x)}). Of uniqueness to this invention, to address inter-dependencies of the Choice Operating Parameters, xi is defined as a scaled Choice Operating Parameter (Λi) using the scaling factor Si; where, initially: x0−i≡SiΛ0−i; further discussed in Paragraph 0102. By design, Choice Operating Parameters which directly effect system stoichiometrics are used by The Input/Loss Method to compute certain parameters which impact the determination of system efficiency. These certain parameters are termed System Effect Parameters and, for the Preferred Embodiment, include four general types and their associated reference values: the L Factor concept (Lk1); the As-Fired fuel flow (mAF); the higher heating value (HHVk3); and the As-Fired fuel water fraction (WFH2O). Two L Factors are presented, L′Fuel and LCO2 as defined by Eqs.(272) & (273). The higher heating value is chosen as either: an uncorrected As-Fired value, HHVAF; a Dry value, HHVDRY; and/or a MAF value, HHVMAF. As the Preferred Embodiment, the As-Fired fuel water fraction is selected only when determining tube leakage flow rates. However, the system operator may select from any one or more or all of these System Effect Parameters (including any one or more or all of the heating values), whose differences with respect to reference values are minimized by altering the selected Choice Operating Parameters through minimization techniques. The minimization techniques are structured to minimize differences between a System Effect Parameter and its corresponding “Reference System Effect Parameter” (termed: L′Fuel-Ref, LCO2-Ref, mAF-PLT, HHVk3-Ref, and WFH2O-Ref). System Effect Parameters are chosen such that they reflect influences on system efficiency through Choice Operating Parameters, and, at the same time, reflect inter-dependencies of the Choice Operating Parameters. For example: changes in the concentration of effluent CO2 (defined as Λ1S or Λ1B) if caused by changes in the fuel will effect computed fuel chemistry, computed heating value and computed boiler efficiency, all of which impact system efficiency. However a change in CO2 (measured as a wet concentration) may be caused by a change in the concentration of effluent H2O (defined as Λ2S or Λ2B), or a change in combustion air effecting effluent O2 (defined as Λ7B or Λ7S), or a change in the fuel's Sodium to Carbon Ratio (defined as Λ9), where any of these changes themselves may directly effect computed fuel flow and fuel chemistry; thus the importance of inter-dependency of Choice Operating Parameters.

[0286] The following summarizes the objective functionalities for the Preferred Embodiment, demonstrating the aforementioned principles:

[0287] F({right arrow over (x)})=ΣiεIƒ[Si, J0(λL), J0(λW), J0(λH), J0(λT)]

[0288] ΛL=ƒ[Lk1, Lk1-Ref, ML]

[0289] ΛW=ƒ[mAF, mAF-PLT, ΔMAF, MW]

[0290] ΛH=ƒ[HHVk3, HHVk3-Ref, MH]

[0291] ΛT=ƒ[WFH2O, WFH2O-Ref, MT]

[0292] The symbol ΣiεI is defined following Eq.(203). Note that as F({right arrow over (x)}) is minimized the quantities {right arrow over (Λ)} are updated in turn (Λi=xi/Si), thus allowing System Effect Parameters to be computed leading directly to the computation of λL, λW, λH and λT. System Effect Parameters have general dependency on Reference Fuel Characteristics, including the following important inter-relationships: computed fuel chemistry is dependent on several or all Choice Operating Parameters, {right arrow over (Λ)}; computed heating values (HHVk3, LHVP and HHVP) are dependent on fuel chemistry, thus {right arrow over (Λ)}; and boiler efficiency (ηB) determined using the methods herein and those of '429, is dependent directly on Λi effluents CO2 and O2, is also dependent on fuel chemistry, and is also dependent on heating value, thus {right arrow over (Λ)}. All of these quantities (fuel chemistry, heating values and boiler efficiency) are also dependent on in-leakage of working fluid into the combustion path, the terms bZ and bT. Working fluid energy flow (BBTC), the Firing Correction HBC and the Δenthalpy term associated with tube failure energy flow (mTΔh), are all dependent on Operating Parameters.

[0293] xi≡SiΛi

[0294] Lk1=ƒ[fuel chemistry({right arrow over (Λ)})]

[0295] mAF=ƒ[(BBTC−mTΔh), ηB ({right arrow over (Λ)}), HHVP({right arrow over (Λ)}), HBC]

[0296] HHVk3=ƒ[fuel chemistry({right arrow over (Λ)})]

[0297] WFH2O=ƒ[tube leakage flow rate (Λ8), fuel chemistry({right arrow over (Λ)})].

[0298] System Effect Parameters

[0299] As discussed, System Effect Parameters include four general types of parameters and their associated reference values: the L Factor (Lk1); the As-Fired fuel flow (mAF); the higher heating value (HHVk3); and the As-Fired fuel water fraction (WFH2O). The most important of these is the L Factor, used routinely for most situations. The higher heating value may be employed, for example, when the thermal system is operating under controlled conditions (e.g., under a testing program), in which its fuel is well characterized. Also, during initial installation of a Calculational Engine, heating value may be used for scoping the range of reasonable correction factors. Fuel flow is discussed below. The As-Fired fuel water fraction is typically used for the special case of computing the tube leakage mass flow rate, mT.

[0300] The L Factor is important in reducing the impact of the shallow valley problem found with fossil-fired systems. An important reason for this is that both L′Fuel and LCO2 have been demonstrated to have remarkably small standard deviations for fossil fuels, including black liquor fuels. For example, some Ranks of coal have L′Fuel standard deviations as low as ±0.05%. Although black liquor dry chemistries may vary widely depending on the particular process, a given process will typically produce a highly consistent L Factor (either L′Fuel or LCO2). Use of L′Fuel as a System Effect Parameter is the Preferred Embodiment, that is when L′Fuel is computed using Eq.(272) or its equivalence. For black liquor fuels the methods taught in '563 and '879 require modification to address the fuel's inorganic compounds, although the principles taught in '563 and '879 remain and are applicable. Specifically the L Factor is computed based on theoretical combustion of dry fuel, thus: φR=αMAF-2=αDRY-2=bA=bZ=bTGAct=eAct=mAct=kBL=v=0.0; and σN=RAct=1.0; nd ΦAct=ΦRef. 25$\begin{array}{cc}{L}_{\mathrm{Fuel}}^{\prime }\equiv {10}^6\left[x_{\mathrm{DRY­theor}}{N}_{\mathrm{DRY­Fuel}}+a_{\mathrm{DRY­theor}}\left(1.0+{\phi }_{\mathrm{Ref}}\right)N_{\mathrm{Air}}-\left(J_{\mathrm{theor}}{N}_{\mathrm{H20}}+x_{\mathrm{DRY­theor}}{\alpha }_{\mathrm{MAF­}10}{N}_{\mathrm{Ash}}+{\sigma }_N{T}_{\mathrm{NC}}{N}_{\mathrm{NaCl}}+{\sigma }_N\left(T_{\mathrm{NO}}+k_{\mathrm{BL}}\right)N_{\mathrm{Na2CO3}}+{\sigma }_N{T}_{\mathrm{NA}}{N}_{\mathrm{Na2SO4}}+{\sigma }_N{T}_{\mathrm{KO}}{N}_{\mathrm{K2CO3}}\right)\right]/\left(x_{\mathrm{theor}}{N}_{\mathrm{Fuel}}{\mathrm{HHV}}_{\mathrm{AF}}\right)& \left(272\right)\end{array}$

[0301] In addition to L′Fuel, the LCO2 factor has shown promise when used as a System Effect Parameter as based exclusively on CO2, given in common units of measure as (lbm dry CO2 effluent)/million-BtuFuel. Note that Dtheor is the effluent CO2 based on theoretical combustion of dried fuel (computed with the same philosophy as used to determine Jtheor).

L CO2 ≡106[DtheorNCO2]/(xDRY-theorNDRY-FuelHHVDRY)  (273)

[0302] The following identities have been found useful in determining the L Factor, and for reducing the solution problem associated with Eq.(19BL): 26$\begin{array}{cc}{x}_{\mathrm{MAF­theor}}{N}_{\mathrm{MAF­Fuel}}{\mathrm{HHV}}_{\mathrm{MAF}}=x_{\mathrm{theor}}{N}_{\mathrm{Fuel}}{\mathrm{HHV}}_{\mathrm{AF}}& \left(274\right)\\ =x_{\mathrm{DRY­theor}}{N}_{\mathrm{DRY­Fuel}}{\mathrm{HHV}}_{\mathrm{DRY}}& \left(275\right)\end{array}$

[0303] When black liquor is burned in a bomb calorimeter for determination of its heating value, ideal solid products—as assumed when using Eq.(272) or (273) based on Eq.(19BL), consist of the following inorganics: Na2CO3, Na2SO4, NaCl and K2CO3. However, other ideal products could be produced depending on the nature of the fuel, its mix of components, combustion conditions and the nature of its chemical equilibrium. Whether potassium forms K2CO3, K2SO4 and/or K2S, or forms the same compounds as sodium in the same proportions, or whether fuel potassium is simply combined with fuel sodium as is a common chemistry assumption, does not affect the general scope and spirit of the present invention. Three tools may be used to determine such ideal compounds as would then be used to modify Eq.(272), the boiler efficiency computations, Eq.(19BL), and throughout this disclosure: 1) direct laboratory analysis of the bomb calorimeter's residue; 2) operational experience; and/or 3) modeling ideal combustion of black liquor using the software “HSC Chemistry 5” commercially available from Outokumpu Research Oy, Pori, Finland (www.HSC@Outokumpu.com).

[0304] Along with the L Factor, the system operator may also choose, in any combination, the plant's indicated fuel flow, the As-Fired heating value, the dry heating value, the MAF heating value, and/or the fuel's average water fraction as System Effect Parameters. Although the system operator has complete flexibility, with this flexibility must apply common engineering judgement. For example, optimizing effluent water against HHVMAF or HHVDRY (heating values without water) would make little sense given the lack of connectivity.

[0305] When selecting the system's indicated fuel flow (mAF-PLT) as a Reference System Effect Parameter, the recovery boiler's operator should proceed with caution. The flow measurement of black liquor feeding a recovery boiler requires an assessment as to its consistency and whether the flow measurement encompasses the total feed of fuel. This invention teaches that the minimization techniques may be used to minimize the difference between a computed fuel flow (mAF) and the system's indicated fuel flow (mAF-PLT), with an off-set ΔMAF, through converged Choice Operating Parameters. Thus, the method of this invention allows use of the system's indicated fuel flow to aid in the determination of computed fuel chemistry, fuel heating value and system efficiency. Although not required, for many situations it is the Preferred Embodiment that use of the system's fuel flow be accompanied with either L′Fuel or LCO2 of Eqs.(272) or (273), to assist with stability and reasonableness of solution. To further enhance stability and reasonableness of solution the recovery boiler's operator may option to limit the range of fuel concentrations determined by the methods of this invention. The operator may also limit the numerical range of the selected Choice Operating Parameters when using the Simulated Annealing technique. Further, to aid in addressing a possible error in MAF-PLT, both an off-set (ΔMAF) in this parameter and a Dilution Factor (MW) may be applied to the relationship between mAF and mAF-PLT; see Eq.(202B) below.

[0306] In summary, the process involving the minimization of differences in System Effect Parameters, by optimizing Choice Operating Parameters, results in converged Choice Operating Parameters ΛF−i; correction factors, Ci, are then determined as based on the ratio of the converged Choice Operating Parameter to their initial values (Λ0−i). Λ0−i are based on the system's raw instrumentation signal, a previous converged solution, an estimate, a guess or as otherwise determined.

C i ≡ΛF−i/Λ0−i  (201)

[0307] Minimization Techniques, Formulations

[0308] This sub-section presents general discussions of the multidimensional minimization techniques and detailed formulations useful to the plant operator in minimizing errors in System Effect Parameters.

[0309] The BFGS technique represents a second generation of multidimensional minimization techniques. As such, it is considered one of the most robust of techniques for a well conditioned problem. The particular BFGS technique employed by The Input/Loss Method has a superior reputation for convergence. The only input parameters the user need be concerned with are the initial relative step-length and the change in the relative step-length. A well-chosen initial relative step-length will prevent long iterations (a value of 0.100 to 0.200 is recommended). The change in the relative step-length impacts resolution of the shallow valley problem, and may be varied until proper convergence patterns are established. A value between 0.010 to 0.040 for the change in the relative step-length has been found to be satisfactory when used in conjunction with the scaling techniques taught herein. The BFGS technique is the preferred method for use on a continuous bases after the problem has been properly conditioned with scaling factors, and selections of Choice Operating and System Effect Parameters have been established. These input parameters are also applicable to the generic Conjugate Gradient technique.

[0310] The generic Conjugate Gradient technique represents a first generation of multidimensional minimization techniques. For numerical processing reasons the BFGS technique has been demonstrated to be superior in to the generic Conjugate Gradient in convergence techniques and accuracy. However, there may be situations in which a generic Conjugate Gradient may be useful as an alternative once the problem has been conditioned.

[0311] The Newton-Raphson method is one of the oldest and simplest multidimensional minimization techniques. This method requires the objective function's compounded vector gradient, resulting in a Jacobian determinant. Generally it will yield an efficient means of convergence but requires reasonable initial Choice Operating Parameters (Λ0−i); however, without such reasonableness it may fail wildly. Newton-Raphson is recommended for possible use only after the BFGS technique has failed to meet its convergence criteria It has applicability given its use of the Jacobian determinant, through which forming explicit inter-dependencies between System Effect Parameters and all Choice Operating Parameters are employed. This assures computed dependencies, if such dependencies exist. This intrinsic feature has been found to be of importance when resolving certain power plant problems. The Preferred Embodiment is to automatically default from BFGS, given failure to meet its convergence (typically due to a lack of established inter-dependencies of Choice Operating Parameters) to, first, the Newton-Raphson, and then in-turn, given its failure, to Simulated Annealing. Newton-Raphson may also be used for scoping initial installations of The Input/Loss Method given difficult combinations of System Effect and Choice Operating Parameters.

[0312] The Simulated Annealing technique, because it employs a global, constrained search methodology, is the Preferred Embodiment for initial study of a new Input/Loss installation and for determining tube leakage flow rate. It may also be used to assist in the selection of which Choice Operating Parameters are best for a particular thermal system. This procedure simulates the annealing process of metal, requiring the controlled reduction of a pseudo-temperature (herein termed “pseudo-T”) to achieve a desired result (i.e., achieving a minimum potential energy of the metal's structure when slowly cooled, thus the minimizing of an objective function). This is a brute force approach involving random search; gradients are not used. As a global optimization procedure it may move both downhill and uphill (that is, it may move both towards and away from local optima), resulting in distinction between different local optima. Conventional optimization techniques (BFGS, generic Conjugate Gradient and Newton-Raphson) only move downhill when minimizing an objective function. Conventional techniques are blind to a global solution in the sense they immediately choose the downhill direction. When addressing fossil-fired combustion problems this may lead to optimizing on the most sensitive of a given selection of Choice Operating Parameters (most likely CO2, thus Λ1S or Λ1B). Distinction between different local optima is accomplished by first starting with initial Λ0−i values, then successively evaluating randomly acquired changes, {right arrow over (Λ)}, but which fall within user-defined step-lengths. Initially this results in coarse study of the objective function, employing large step-lengths, requiring repeated evaluations with seemingly little progress. In the process of choosing {right arrow over (Λ)} values the algorithm generally attempts to move downhill, however it also moves uphill in a probabilistic manner to escape local optima Step-lengths are dynamically chosen such that half of all uphill moves are randomly accepted, again helping to ensure that the function escapes local optima. As the annealing process proceeds and the algorithm closes on the global optimum, step-lengths decrease as the pseudo-T decreases requiring even more objective function evaluations as the optimum is approached. By viewing objective functions in general terms and with its ability to move probabilistically uphill, Simulated Annealing solves functions that are otherwise difficult to resolve, including shallow valley problems associated with fossil combustion and tube leakages. However, with such flexibility comes numerous objective function evaluations necessitating long computing times. In addition, converged solutions should be re-tested periodically with different seeds (i.e., initializations of the random number generator) to assure a global optimum.

[0313] When applied to fossil-fired combustion, the more sensitive inputs to the Simulated Annealing technique include the following: starting point Λ0−i values; the number of cycle evaluations (5 is recommended); the minimum and maximum values associated with each Λi (i.e., defining the region containing the optimum); an initial pseudo-T (0.100 is recommended); and the relative change in pseudo-T (i.e., the step-length, 0.010 to 0.020 is recommended). Each of these inputs may be established by sensitivity study to assure a robust solution, or as otherwise determined. Minimum and maximum Λi values may also be established by review of historical system data or through the experience of the system operator. The smaller the range between minimum and maximum Λi values, the tighter the search becomes with the final solution becoming narrowed. This feature is especially useful when As-Fired fuel flow or As-Fired fuel water fraction are selected as a System Effect Parameters (in combination with non-unity Dilution Factors, MW and MT).

[0314] Objective Function and Choice Operating Parameters

[0315] The following paragraphs present the preferred objective functions and their solution methodologies and specify the Choice Operating Parameters which may be employed by minimization techniques. As explained, the Bessel function is used to define the objective function. The Bessel function's argument, as taught by this invention, has been chosen to aid in addressing the shallow valley problem and in convergence of the minimization techniques. The formulations presented produce quantities which may allow numerical inter-dependencies between Choice Operating Parameters ({right arrow over (Λ)}), or not, depending on which minimization technique is chosen. This is important for addressing problems in which initial Choice Operating Parameters lie far from the optimum. This is also important where more than one System Effect Parameter is chosen which may present unique numerical convergence problems.

[0316] For the BFGS, generic Conjugate Gradient, Newton-Raphson and Simulated Annealing techniques the objective function is given by the following. Note that ML, MW, MH and MT are real numbers; ML & MH are typically equal to one, while MW & MT are typically greater than one. Again, the System Effect Parameters, Lk1, mAF, HHVk3 and WFH2O, are functions of a set of Λi.

ΛL=[(Lk1−Lk1-Ref)/Lk1-Ref]ML  (202A)

ΛW=[(mAF−mAF-PLT−ΔmAF)/(mAF-PLT+ΔmAF)]MW  (202B)

ΛH=[(HHVk3−HHVk3-Ref)/HHVk3-Ref]MH  (202C)

ΛT=[(WFH2O−WFH2O-Ref)/WFH2O-Ref]MT  (202D)

[0317] 27 $$ F ⁡ ( x → ) = ∑ i ∈ I ⁢ { S i ⁡ [ 1.0 - J 0 ⁡ ( λ L ) ] + S i ⁡ [ 1.0 - J 0 ⁡ ( λ W ) ] + S i ⁡ [ 1.0 - J 0 ⁡ ( λ H ) ] + S i ⁡ [ 1.0 - J 0 ⁡ ( λ T ) ] } ( 203 )

[0318] In Eq.(202B) ΔmAF is an off-set, or bias, observed in the indicated flow, mAF-PLT. In Eq.(203) and as used elsewhere, the symbol ΣiεI indicates a summation on the index i, where i variables are contained in the set I defined as the elements of {right arrow over (Λ)}. For example, assume the user has chosen the following: Λ1S is to be optimized to minimize the error in L′Fuel and HHVMAF, Λ2S is optimized for L′Fuel and mAF (MW=1.40, ΔmAF=0.0), Λ4 is optimized for L′Fuel, and Λ7B is optimized for L′Fuel. Therefore: {right arrow over (Λ)}=(Λ1S, Λ2S, Λ4, Λ7B), I={Λ1S, Λ2S, Λ4, Λ7B}, thus {right arrow over (x)}=(x1, x2, x3, x4); x1=S1Λ1S; x2=S2Λ2S; x3=S3Λ4; x4=S4Λ7B; where Eq.(203) for this example than becomes: 28$\begin{array}{cc}F\left(\overrightarrow{x}\right)=S_1\left\{\left[1.0-J_0\left({\lambda }_L\right)\right]+\left[1.0-J_0\left({\lambda }_H\right)\right]\right\}+S_2\left\{\left[1.0-J_0\left({\lambda }_L\right)\right]+\left[1.0-J_0\left({\lambda }_W\right)\right]\right\}+S_3\left[1.0-J_0\left({\lambda }_L\right)\right]+S_4\left[1.0-J_0\left({\lambda }_L\right)\right]& \end{array}$

[0319] Derivatives ∂F/∂xi for the BFGS and generic Conjugate Gradient techniques, based on Eq.(203), are given by the following: 29$\begin{array}{cc}\partial F/\partial x_i=\partial F/\left(S_i\partial {\Lambda }_i\right)& \\ =J_1\left({\stackrel{\_}{\lambda }}_L\right)\left[\partial {\lambda }_L/\partial {\Lambda }_i\right]+J_1\left({\stackrel{\_}{\lambda }}_W\right)\left[\partial {\lambda }_W/\partial {\Lambda }_i\right]+J_1\left({\stackrel{\_}{\lambda }}_H\right)\left[\partial {\lambda }_H/\partial {\Lambda }_i\right]+J_1\left({\stackrel{\_}{\lambda }}_T\right)\left[\partial {\lambda }_T/\partial {\Lambda }_i\right]& \left(204\right)\end{array}$

[0320] where, for example:

[∂λW/∂Λi]=MW[({overscore (m)}AF−mAF-PLT−ΔmAF)/(mAF−PLT+ΔmAF)]MW−1(∂mAF/∂Λi)/(mAF-PLT+ΔmAF);

[0321] and where {overscore (λ)}L, {overscore (λ)}W, {overscore (λ)}H, {overscore (λ)}T and {overscore (m )}AF are taken as average values. Gradients, ∂Fi/∂xj, for the Newton-Raphson method, thus defining the Jacobian determinant, are given by the following: 30$\begin{array}{cc}\partial F_i/\partial x_j=\partial F_i/\left(S_j\partial {\Lambda }_j\right)& \\ =S_i{J}_1\left({\stackrel{\_}{\lambda }}_L\right)\left[\partial {\lambda }_L/\left(S_j\partial {\Lambda }_j\right)\right]+S_i{J}_1\left({\stackrel{\_}{\lambda }}_W\right)\left[\partial {\lambda }_W/\left(S_j\partial {\Lambda }_j\right)\right]+S_i{J}_1\left({\stackrel{\_}{\lambda }}_H\right)\left[\partial {\lambda }_H/\left(S_j\partial {\Lambda }_j\right)\right]+S_i{J}_1\left({\stackrel{\_}{\lambda }}_T\right)[\partial {\lambda }_T/\left(S_j\partial {\Lambda }_j\right)& \left(205\right)\end{array}$

[0322] where, for example:

[∂λW/∂Λj]=MW[({overscore (m)}AF−mAF-PLT)/mAF-PLT]MW−1(∂mAF/∂Λj)/(mAF-PLT+ΔmAF).

[0323] In the Preferred Embodiment, selection of Choice Operating Parameters may be made by the system operator from any combination or all of the following:

Λ1S=DAct; Stack CO2 (with effects from air pre-heater leakage)  (211S)

Λ1B=DAct RAct; Boiler CO2 (without effects from air pre-heater leakage)  (211B)

Λ2S=JAct≡j+bAβ; Stack H2O (with moisture from air pre-heater leakage)  (212S)

Λ2B=j RAct; Boiler H2O (without moisture from air pre-heater leakage)  (212B)

Λ3=AF; Air/Fuel mass ratio  (213)

Λ4=RAct; Air Pre-Heater Leakage Factor  (214)

Λ5=AAct; Concentration of O2 in the boundary air  (215)

Λ6=mLS; Indicated limestone flow rate  (216)

Λ7S=GAct=(g+2σNTNS+σNkBL/2)+aβ; Stack O2 (with air pre-heater leakage)  (217S)

Λ7B=RAct(g+2σNTNS+σNkBL/2); Boiler O2 (without air pre-heater leakage)  (217B)

Λ8=mT; Tube leakage flow rate  (218)

Λ9=φN; Sodium to Carbon Ratio in the fuel  (219)

Λ10=φR; Reduction Efficiency  (220)

[0324] The selection of one or more of the Choice Operating Parameters must depend on common understanding of recovery boiler stoichiometrics and associated relationships to physical equipment. Specifically, The Input/Loss Method produces, by employing one or more of the minimization techniques (within the ERR-CALC computer program), converged Choice Operating Parameters and correction factors Ci applied to the initial values Λ0−i. The converged Choice Operating Parameters are then used within the Fuel Iterations to produce a computed fuel chemistry, discussed in conjunction with FIG. 2 and FIG. 3. A monitoring cycle, processing Fuel Iterations, may be scheduled as frequently as desired; each cycle employing correction factors produced by ERR-CALC at the same or slower frequency. For example, ERR-CALC could be processed (producing updated correction factors) once per day, while Fuel Iterations could be processed once every 2 minutes using 15 minute running averages of Operating Parameters.

[0325] In the above paragraph, the phase “common understanding of recovery boiler stoichiometrics and associated relationships to physical equipment” is meant the routine knowledge base a system operator should have concerning his/her thermal system. To thoroughly teach this invention, examples of such common understanding and their associated impacts on this invention follow: if limestone (Λ6) is not used, the system operator would not select limestone flow as a Choice Operating Parameter as such a selection would result in an unity correction factor, non-convergence, warning messages, and/or a faulted condition produced from ERR-CALC; the selection of the Air Pre-Heater Leakage Factor (Λ4) would not be made if the system uses a tubular exchanger which has no air leakage (as designed), and would result in a similar faulted condition; the selection of the air/fuel ratio (Λ3) leading to determination of the fuel inert fraction, and also invoking a constant fuel inert assumption, would not be made as such a selection would result in a similar faulted condition; the selection of Boiler CO2 (Λ1B), an Air Pre-Heater Leakage Factor (Λ4), and Boiler O2 (Λ7B), given that “correcting” the air pre-heater leakage would have no effect on the Boiler-side mix of CO2 and O2, would result in a similar faulted condition; the selection of the Sodium to Carbon Ratio in the fuel (Λ9) would not be made if periodic fuel analyses indicates a constant ratio.

[0326] The use of the exponents ML, MW, MH and MT in Eqs.(202A), (202B), (202C) and (202D), termed Dilution Factors, allows a dilution or dampening of the functionality between Reference System Effect Parameters (Lk1-Ref, mAF-PLT, HHVk3-Ref and WFH2O-Ref) and selected Choice Operating Parameters ({right arrow over (Λ)}). As an important feature of this invention, Dilution Factors allow the numerical processes to recognize that Reference System Effect Parameters may themselves have bias. Examples of such bias include: Reference Fuel Characteristics having been chosen with an out-dated database, biasing the computed reference L Factor; the reference heating value having been determined incorrectly, analyzed incorrectly in the laboratory and/or having intrinsic uncertainties; and the indicated fuel flow having serious instrumentation error. However, engineering judgement and a valid database may be reasonably anticipated and applied in the cases of reference L Factors and reference heating values. Dilution Factors ML (influencing Lk1-Ref) and MH (influencing HHVk3-Ref) may be assumed to be unity for most situations as is preferred; or they may be based on monitoring experience, sensitivity studies or as otherwise determined. Applying engineering judgement and/or a valid database is difficult in the case of the plant's indicated fuel flow; which could have bias, thus the use of ΔmAF in Eq.(202B). Thus MW (influencing mAF-PLT) should be determined based on results from The Input/Loss Method and the processes of this invention, when such results are generically compared to system data. Specifically, MW may be adjusted until Input/Loss computed total effluent flow reasonably agrees and/or tracks the measured, computed combustion air flow agrees and/or tracks the measured, computed fuel flow agrees and/or tracks the indicated fuel flow, and similar system-wide comparisons. In the context of the last sentence, “tracks” is defined as the computed value trending over time with the measured, having a constant off-set. In the case of fuel water fraction as used directly as a System Effect Parameter, WFH2O is driven towards a reference value by optimizing tube leakage flow (Λ8); when so employed, System Effect Parameter WFH2O is considered a special case. Application of Dilution Factors require that the sense of the bracketed terms of Eqs.(202A), (202B), (202C) and (202D) be always positive requiring a reversal of the derivative's sign as appropriate. The Dilution Factors ML, MW, MH and MT are real numbers; ML & MH are typically assumed to be unity, while MW & MT are typically found through sensitivity studies to be non-unity (ranging between 0.90 and 1.20 for MW, and 1.0 to 1.6 for MT).

[0327] In these relationships each Choice Operating Parameter (Λi) is scaled with the parameter Si, determined to be suitable for the BFGS, generic Conjugate Gradient and Newton-Raphson techniques. Scaling for these methods is important for proper application of this invention, as minimization techniques in general are sensitive to variations in the numerical size, and units of measure, of the Λi terms (e.g., for fossil-fired applications, an un-scaled Λ1S may be 0.14 moles-CO2/mole-Dry-Stack-Gas, while an un-scaled Λ6 may be 22,000 lbm/hr). It has been found that a good initial estimate of Si may be developed as the inverse of Λi. Further, the influence of scaling may be improved by employing a pre-scaling factor, si; which may be determined as taught in '879, or as otherwise determined by the system operator through sensitivity studies. It has been found that si for Λ8 typically of 10,000 works well if Λ8 units of measure are lbm/hr. However, it has also been found that the Newton-Raphson technique converges quickly when optimizing the combination of Λ2S or Λ2B, and Λ8, thus may be used to adjust si for Λ8 until appropriate sensitivity is reached between the Choice Operating Parameters of effluent water and tube failure flow rate (that is when one term does not predominate the other).

S i =s i /Λ0−i  (206)

xi≡SiΛi  (207)

[0328] When ERR-CALC is executed using either BFGS, generic Conjugate Gradient or Newton-Raphson techniques typically 5 to 50 iterations are required for convergence. However, when ERR-CALC is executed using Simulated Annealing typically over 1000 iterations are required for convergence. To address the problem of long computing times, associated with any minimization technique, this invention teaches to duplicate within the ERR-CALC program only those calculations which effect System Effect Parameters, and to therefore compute System Effect Parameters within ERR-CALC (which are then repeated within the Fuel Iterations). This results in a considerable reduction in computing time required to evaluate repeated objective function calculations. Specifically, these duplicated calculations include: principally HEATRATE stoichiometrics (which are also used by EX-FOSS); L Factor calculations; heating value calculations; and an approximation of the effects changing stoichiometrics and changing heating value has on boiler efficiency and thus the effects on computed fuel flow using Eq.(330A) or (330B). In summary, these duplicated calculations determine affects on the System Effect Parameters (Lk1, mAF, HHVk3 and WFH2O) of a given set of Choice Operating Parameters ({right arrow over (Λ)}).

[0329] Applicable references for the preferred minimization techniques include the following sources. For the BFGS and the generic Conjugate Gradient techniques the references are: D. F. Shanno and K. H. Phua, “Algorithm 500, Minimization of Unconstrained Multivariate Functions”, ACM Transactions on Mathematical Software, Vol.2, No.1, March 1976, pages 87-94; and D. F. Shanno and K. H. Phua, “Remark on Algorithm 500, Minimization of Unconstrained Multivariate Functions”, ACM Transactions on Mathematical Software, Vol.6, No.2, December 1980, pages 618-622. For the Simulating Annealing technique the references are: W. L. Goffe, G. D. Ferrier and J. Rogers, “Global Optimization of Statistical Functions with Simulated Annealing”, Journal of Econometrics, Vol.60, No.1/2, pp.65-100, January/February 1994; for its base technology see: A. Corana, M. Marchesi, C. Martin and S. Ridella, “Minimizing Multimodal Functions of Continuous Variables with the ‘Simulated Annealing’ Algorithm”, ACM Transactions on Mathematical Software, Vol.13, No.3, pp.262-280, September 1987; for modifications to the random number generator RANMAR which is employed by Simulating Annealing see: F. James, “A Review of Pseudorandom Number Generators”, Computer Physics Communications, Vol.60, pp.329-344, 1990. For the Newton-Raphson technique the reference is: W. H. Press, S. A. Teukolsky, W. T. Vettering & B. P. Flannery, Numerical Recipes in FORTRAN 77. The Art of Scientific Computing, Cambridge University Press, Cambridge and New York (1992), Chapter 9.6 on Newton-Raphson Method for Nonlinear Systems of Equations, and Chapter 9.7 on Globally Convergent Methods for Nonlinear Systems of Equations.

[0330] Additional minimization techniques and teachings of related mathematical procedures which may be applied to this invention, are presented in the following: J. Nocedal and S. J. Wright, Numerical Optimization, Springer-Verlag, New York (1999); G. N. Vanderplaats, Numerical Optimization Techniques for Engineering Design, McGraw-Hill Book Company, New York (1984); and W. H. Press, S. A. Teukolsky, W. T. Vettering & B. P. Flannery, Numerical Recipes in FORTRAN 77The Art of Scientific Computing, Cambridge University Press, Cambridge and New York (1992). Other common minimization techniques involving constrained or unconstrained searches may also be alternatively applied. These include Sequential Linear Programming, Direction Set using Powell's method, Simplex method, Downhill Simplex method, Simplex method with product form inverse, Quasi-Newton method, and others. Commercial products are also available, such as from Lindo Systems, Inc. of Chicago, Ill.

[0331] A further technique applicable to the reduction of errors in Choice Operating Parameters lies with use of artificial neural network technology (herein termed ANN). As traditionally employed at power plants ANN technology learns through a database how to minimize defined parameters: a change in a coal mill's air registration produces an observed result in another parameter. Generically, the minimization of defined parameters means to drive such parameters in one, but constrained, direction: the lowest combustion air flow (via bias on FD Fans) for a given power; the highest boiler efficiency by minimizing (1.0−ηB); etc. However, as used for this invention ANN technology is used to correct Choice Operating Parameters such that errors in System Effect Parameters are reduced. Explicit determination of fuel chemistry, computed heating value, computed boiler efficiency, consistent mass and energy balances, and, in general, explicit thermodynamics remain as taught in '994 and '429. ANN may be applied to recognize patterns in computed System Effect Parameters influenced by causal Choice Operating Parameters. Much like the aforementioned (and preferred) techniques, ANN technology may make corrections to initial Choice Operating Parameters to achieve a desired result [for example, to minimize the λL, λW, λH and/or λT terms of Eqs.(202A), (202B), (202C) and (202D)]. These corrections, Ci, may be based on choosing the highest probability a set of Λi will produce the lowest errors in System Effect Parameters relative to Reference System Effect Parameters. An advantage to ANN is that such corrections are learned; that is, ANN improves its correlations with an ever increasing database. Typically such learning may be done without use of an objective function, but not always as in the case of object oriented ANN. Specifically, Choice Operating Parameters used to compute fuel chemistries and heating values (leading to boiler efficiency and system efficiency), may be analyzed for their influences on System Effect Parameters, patterns then recognized which would lead directly to Ci corrections being applied via Eq.(201). Given such corrections, The Input/Loss Method would then proceed as described herein, and in '994 and '429 as applicable.

[0332] Numerous commercial ANN technology software packages are available, for example from: NeuralWare of Pittsburg, Pa.; California Scientific Software of Nevada City, Calif.; The MathWorks, Inc. of Natick, Mass.; Pegasus Technology of Mentor, Ohio a subsidiary of KFx, Inc; NeuCo, Inc. of Boston, Mass.; those available from universities; and those to be found on the internet. A particularly applicable ANN technology is available from Computer Associates of Islandia, N.Y. comprising their Neugents technology. In addition, any ANN technology which allows for object oriented programming are directly applicable to this invention as such objects, as mathematical kernels, may explicitly correct a variety of Choice Operating Parameters employed by one of the Input/Loss methods.

[0333] However, ANN technology is not the Preferred Embodiment given that such technology is historically intended for large databases, databases representing processes too complex for explicit thermodynamics and/or databases those applicable objective functions are unknown or otherwise may not be readily discerned. The teachings of the Preferred Embodiment of this invention may be applied directly using ANN technologies which have application following the general scope and spirit of the present invention.

[0334] Tube Leakage Flow Rate Computations

[0335] To determine a tube leakage flow rate using the direct application of the aforementioned techniques, that is selecting the Choice Operating Parameter for tube leakage flow (Λ8), in combination with others, as optimized for any or all System Effect Parameters Lk1, mAF, HHVk3, and/or WFH2O has been demonstrated at power plants burning fuels whose hydrocarbons and fuel water are similar to that of black liquor. However the strict application of such techniques is not preferred when determining both the flow rate and the location of the leaking tube. The preferred technique for determining a tube leakage flow rate is presented in three steps (termed “Passes”) detailed below. The technique for determining the location of the failed tube within the steam generator is detailed by separate sub-section. In summary, for the Preferred Embodiment, once a possible Tube Failure Mechanism has been identified (via TABLE 1A or TABLE 1B): a tube leakage flow rate is determined by optimizing Λ8, in combination with other Choice Operating Parameters except for Λ1S, Λ1B, Λ2S and Λ2B. Nominal correction factors to effluent CO2 and effluent water are obtained from historical evidence. This achieves stoichiometric balance, an initial fuel chemistry and heating value assuming the nominally corrected effluent CO2 and effluent water are reasonably accurate. Next a final fuel chemistry and heating value are determined but this time as influenced by the determined tube leakage flow rate and all routine Choice Operating Parameters, except Λ2S or Λ2B; that is, Choice Operating Parameters as would be routinely selected whose inter-dependencies are now effected by an established tube leakage. As defining a “separative analyses process”, it was developed to address the situation where effluent water, based on either a measurement or an assumption, was being corrected without regard to how such a correction might influence other Choice Operating Parameters, especially tube leakage and the important effluent CO2. For example, if in correcting a high effluent water signal (whose Λ0−2S value reflects an actual tube failure) to a lower nominal value, the resultant dry-based effluent CO2 may become badly skewed effecting computed heating value. The preferred process first accepts the effluent water value using an historically based correction factor, C2S-hist, i.e., not optimizing on Λ1S, Λ1B, Λ2S or Λ2B, but optimizing on Λ8 and all other Choice Operating Parameters. This optimization establishes a computed tube leakage flow rate, consistent fuel chemistry and a heating value given a tube leakage. The computed tube leakage could be essentially zero if determined to be stoichiometrically consistent. The process then repeats but including Λ1S or Λ1B and other Choice Operating Parameters, again except effluent water, and using the computed tube leakage flow rate. This final process then reflects nominal values given the constrained methodology used by Simulated Annealing; for example, computed fuel carbon established in-part from Λ1S will lie between αMAF-4/min and αMAF-4/max; effluent water is nominally corrected and is consistent with the computed tube leakage.

[0336] The preferred technique, versus a straightforward application of '879, addresses several problem areas found during the demonstration of this invention: the marked insensitivity of small tube leakages on system stoichiometrics; correction factors being adversely influenced by an actual tube leakage, but the resulting effects of converged Choice Operating Parameters on stoichiometrics would mask detection of tube leakage; shallow valley problems aggravated by tube failures; and marked differences, and possible problems, associated with scaling Choice Operating Parameters especially with widely varying tube leakage flow rates (e.g., from 2,000 to 100,000 lbm/hr). Further, this separative analyses process lends itself to the determination of the tube leak's location within the steam generator.

[0337] Pass 0 of the Preferred Embodiment represents a typical monitoring cycle using The Input/Loss Method, but where the user has optioned for tube failure checking at each execution of the ERR-CALC program. A typical monitoring cycle involves minimizing errors in System Effect Parameters (Lk1, mAF, HHVk3 and/or WFH2O) by optimizing a selection of routine Choice Operating Parameters applicable to the thermal system and its fuel, but not use of Λ8. For such routine monitoring mT=bT=0.0 is assumed for the stoichiometrics of Eqs.(19BL) and (200). The optimization employs any multidimensional minimization technique, and/or ANN technology, appropriate to the thermal system and its fuel. This type of monitoring, using routine Choice Operating Parameters, is taught in '879 and herein, and typifies normal use of ERR-CALC. This process defines Pass 0 logic described in FIG. 3.

[0338] If tube failure checking (the option having been activated in Pass 0) has detected a Tube Failure Mechanism (see TABLE 1A and TABLE 1B), a Pass 1 process is then begun to determine a computed tube leakage flow rate. Within this Pass 1 routine Choice Operating Parameters used in Pass 0 are selected with Λ8, but excluding effluent CO2 and effluent water (Λ1S, Λ1B, Λ2S and Λ2B). Choice Operating Parameter Λ8 is optimized to drive the As-Fired fuel water fraction, WFH2O, or the computed fuel flow, mAF, to their respective reference values: WFH2O-Ref or (mAF-PLT+ΔmAF) Corrections to effluent CO2 and effluent water are set to historical values associated with established instrumentation experience (or an assumption) not reflective of tube failure. For example, effluent Stack water may be equated to C2S-histΛ0-2S where C2S-hist is a nominal correction factor and where Λ0-2S is based on plant data. Corrections to effluent CO2 and effluent water are thus held constant for Pass 1. Other correction factors, Ci, established in Pass 0 are used with initial Λ0−i estimates. This process results in a computed tube leakage flow rate which satisfies stoichiometric balances of Eqs.(19BL) and (200), providing a solution in which fuel chemistry and heating value fall within numerical constraints (min/max limits). Note that other System Effect Parameters may be used, optimizing Λ8, such as HHVAF, however either WFH2O or mAF have proven to be appropriately sensitive to the stoichiometrics and represent the Preferred Embodiment. This process defines Pass 1 logic described in FIG. 3.

[0339] After convergence of Pass 1, Pass 2 then re-establishes System Effect Parameters with the previous selection of routine Choice Operating Parameters, but excluding effluent water (those effects are now replaced by the computed tube leakage flow rate). The Simulating Annealing algorithm is preferred; however for this Pass 2 BFGS has been observed to be adequate if properly scaled. Convergence results in converged Choice Operating Parameters, noting that all Choice Operating Parameters which were being used in Pass 0 are now influenced by a tube leakage flow rate (mT), including a nominally corrected effluent water. Pass 1 and Pass 2 employ Simulated Annealing technique as the Preferred Embodiment given its ability to address: shallow valley problems; the possibility of considerably different scaling associated with small to large tube failure flows; its constrained search methodology in which the Choice Operating Parameters are numerically bound; and that a variety of Choice Operating Parameters which may have wide ranging numerical values. This process defines Pass 2 logic described in FIG. 3. Note that whenever mT is determined greater than zero, procedures as taught by this invention then incorporate the influence of such leakage into combustion stoichiometrics through the bT molar quantity determined as follows:

b T =m T (xNAF)/(NH2OmAF)  (321A)

m T =b T (NH2OmAF)/(xNAF)  (321B)

[0340] The bT quantity, through use in Eq.(19BL), then effects boiler efficiency, computed fuel flow and efficiency computations in the same manner as the bZ quantity as taught in '994 and '429; bZ in '994 and '429 being replaced by the quantity (bZ+bT). Thus if the quantities: NAF; the molecular weight of the working fluid (NH2O, water is assumed); mAF (or mAF-PLT if the computed is not available); the tube leakage moles (bT); and resolution of the stoichiometric model of the combustion process, then Eq.(321B) may be solved for mT without use of a minimization or ANN technique.

[0341] To further assist in teaching this invention, TABLE 2 presents a typical scenario of routine monitoring, the identification of a possible tube leak, and then the resolution of the tube leakage flow rate. In TABLE 2, the second column denotes the selection of Choice Operating and System Effect Parameters; for example, “Λ1S min L′Fuel” means that Choice Operating Parameter Λ1S, see Eq.(211S), is selected to minimize the error in System Effect Parameter L′Fuel. The use of “Λ1S min L′Fuel”, “Λ2S min L′Fuel” and “Λ4 min L′Fuel” used in Pass 0 is typical for the assumed thermal system. If measuring Stack O2 instead of Boiler O2 as used in TABLE 2, the selection would typically consist of only “Λ1S min L′Fuel” and “Λ2S min L′Fuel”.

TABLE 2
Example of Determining Tube Failure Flow Rate using the Preferred Embodiment
The Thermal System and
Computational Sequence	Optimizations	Comments
Routine monitoring of a black liquor-	Λ1S min L′Fuel	Compute As-Fired fuel chemistry with
fired system having high water, low	Λ2S min L′Fuel	constant fuel inerts every 3 minutes;
& constant inerts; multiple O2	Λ4 min L′Fuel	optimization of Choice Operating
instruments are used at the Boiler	Λ9 min L′Fuel	Parameters every 30 minutes using
having high accuracy; CO2 & H2O		BFGS technique with Tube Failure
Stack instruments; uncertain air		Model option invoked. Historically:
leakage. Defines Pass 0.		C1S = 0.96, C2S = 1.02 and C4 = 1.10.
A possible tube failure has been	Λ4 min L′Fuel	Compute As-Fired fuel chemistry with
detected given Tube Failure	Λ8 min WFH2O	constant fuel inerts; optimization using
Mechanism 51 (see TABLE 1A). Use		Simulated Annealing. Results in the
historical values for C1S & C2S;		computed tube leakage flow rate, mT,
Λ0-8 = 500 lb/hr. Defines Pass 1.		which satisfies stoichiometrics.
Return to routine monitoring but	Λ1S min L′Fuel	Compute As-Fired fuel chemistry as in
including the computed tube leakage	Λ4 min L′Fuel	Pass 0, resulting in converged Choice
flow rate, mT, but excluding effluent		Operating Parameters associated with
water (C2S = 1.02). Defines Pass 2.		the computed tube leakage flow rate.

[0342] ///

[0343] ///

[0344] ///

[0345] Tube Leak Location

[0346] A important feature of this invention is its ability to assess the impact of a tube leak on the thermal performance of the system, and where within the steam generator the leak occurs. Once a tube leakage flow rate has been determined, its impact on the total energy flow to the working fluid and on boiler efficiency may be determined; thus its effects on fuel flow and system efficiency may be understood. If a thermal system's feedwater flow is held essentially constant, then a developing tube leak will result in less total energy flow required from the combustion gases; i.e., a reduction in the working fluid's energy flow required to meet the same working fluid boundary pressures and temperatures. If the working fluid energy flow without tube leakage is termed BBTC, then the actual energy flow, assuming a tube leak, is given by: (BBTC−mTΔh), where Δh is the enthalpy difference between the outlet of the last heat exchanger effected by the leakage, hLast, and the first exchanger so affected, hStreamT (i.e., the heat exchanger in which the leak occurs); mTΔh is the energy flow lost from the working fluid due to tube leakage. The enthalpy of the leaking fluid as it enters the combustion gas path, hStreamT, is assumed, by choice, to be the same as the heat exchanger's inlet enthalpy (any location may be applied). When applied to The Input/Loss Method of computing boiler efficiency, the enthalpy of the leaking fluid entering the combustion gas path must be properly referenced as taught by Eq.(137) using (hStreamT−hf-cal), and as taught by Eqs.(131) & (132) using (hStack−hSteamT). Quantitative effects on boiler efficiency and system heat rate are not obvious and may not be off-setting, they are preferably described through use of The Input/Loss Method of '994 and '429; they may also be described through use of other Input/Loss methods, or may be described through use of any method which determines boiler efficiency and system heat rate (and therefore system efficiency) given knowledge of the fuel's heating value and other dependent quantities. Computed As-Fired fuel flow and system heat rate are then determined by the following, assuming a tube leak (versus no leak associated with Eq.(128)): 31$\begin{array}{cc}{m}_{\mathrm{AF}}=\left(\mathrm{BBTC}-m_T\Delta h\right)/\left[{\eta }_{B\text{-}\mathrm{HHV}}\left(\mathrm{HHVP}+\mathrm{HBC}\right)\right]& \left(330A\right)\\ =\left(\mathrm{BBTC}-m_T\Delta h\right)/\left[{\eta }_{B\text{-}\mathrm{LHV}}\left(\mathrm{LHVP}+\mathrm{HBC}\right)\right]& \left(330B\right)\\ {\mathrm{HR}}_{\mathrm{HHV}}=m_{\mathrm{AF}}\left(\mathrm{HHVP}+\mathrm{HBC}\right)/W_{\mathrm{output}}& \left(331A\right)\\ =\left(\mathrm{BBTC}-m_T\Delta h\right)/\left({\eta }_{B\text{-}\mathrm{HHV}}{W}_{\mathrm{output}}\right)& \left(331B\right)\\ {\mathrm{HR}}_{\mathrm{LHV}}=m_{\mathrm{AF}}\left(\mathrm{LHVP}+\mathrm{HBC}\right)/W_{\mathrm{output}}& \left(332A\right)\\ =\left(\mathrm{BBTC}-m_T\Delta h\right)/\left({\eta }_{B\text{-}\mathrm{LHV}}{W}_{\mathrm{output}}\right)& \left(332B\right)\end{array}$ ηSYS-HHV=3412.1416/HRHHV  (333A)

ηSYS-LHV=3412.1416/HRLHV  (333B)

[0347] It is important to this invention to recognize that the location of the tube failure affects the working fluid's energy flow. The typical recovery boiler used in the paper processing industry routes the working fluid first through an economizer heat exchanger, then through a series of water wall heat exchangers and the boiler drum (if a sub-critical unit). It could happen that additional heat exchangers are also employed, routing the working fluid through a primary super-heater, through a final super-heater, and, lastly, through a reheat exchanger. Although the preceding is typical, numerous configurations of different heat exchanger names and types may be found, including, for example: lower economizers, upper economizers, primary secondary super-heaters, final secondary super-heaters, division walls, front reheaters, rear reheaters, etc. Such heat exchangers need to be identified by order. For example, if a tube leak occurs in an economizer, its loss is seen throughout the steam generator (having the greatest impact on working fluid energy flow, e.g., Δh=(hReheat-outlet−hSteamT), where hSteamT may be less than 490 Btu/lbm. If a tube leak occurs in the final reheat exchanger, its loss only affects this last exchanger (having the least impact on working fluid energy flow and thus on computed fuel flow, Δh=(hReheat-outlet−hSteamT), where hSteamT may be≈1320 Btu/lbm with a hReheat-outlet of 1520 Btu/lbm. This suggests that Eq.(330A) or (330B) has an unique solution dependent on the assigned (and actual) location of the tube leakage. When using one of the appropriate Input/Loss methods, such dependency on a location of the tube leak may be intrinsically a function of computed fuel chemistry and Firing Corrections, and thus, also a function of the resultant heating value and boiler efficiency.

[0348] Determination of the location of the tube leak is accomplished by recognizing that certain System Effect Parameters are a function of the working fluid energy flow (as affected by tube failure flow rate and its location). The System Effect Parameter of computed fuel flow, mAF, is a function of (BBTC−mTΔh) through Eq.(330A) or (330B). The System Effect Parameter of As-Fired fuel water fraction, WFH2O, is a function of the (BBTC−mTΔh) term through affects on boiler efficiency (ηB-HHV), heating value (HHVAF) and Firing Corrections (HBC). Although not obvious, The Input/Loss Method, because it determines fuel chemistry, heating value, boiler efficiency and Firing Corrections independent of fuel flow, and with great consistency, must, never-the-less, affect computed boiler efficiency consistently as taught herein and in '429; illustrating using higher heating value: 32$\begin{array}{cc}{\eta }_{B\text{-}\mathrm{HHV}}=\left(\mathrm{BBTC}-m_T\Delta h\right)/\left[m_{\mathrm{AF}}\left(\mathrm{HHVP}+\mathrm{HBC}\right)\right]& \left(334\right)\\ =\left(-{\mathrm{HPR}}_{\mathrm{Act}\text{-}\mathrm{HHV}}+{\mathrm{HRX}}_{\mathrm{Act}\text{-}\mathrm{HHV}}\right)/\left(\mathrm{HHVP}+\mathrm{HBC}\right)& \left(335\right)\end{array}$

[0349] Eq.(335) must reflect a consistently computed boiler efficiency; just as Eq.(334) as composed of a term which directly reflects tube failure location (BBTC−mTΔh), and therefore reflects a consistently computed boiler efficiency. The Enthalpy of Products and the Enthalpy of Reactants terms of Eq.(335), HPRAct and HRXAct, are computed with terms influenced by both the tube leakage flow and its location via the Δh term. HPRAct includes the enthalpy of all water exiting the system (hStack), relative to the enthalpy at associated entry points into the combustion gas path, including that for tube leaks at hSteamT; thus: (hStack−hSteamT) as found in Eqs.(131) & (132). HRXAct of Eqs.(135) & (136) includes the Firing Correction term which encompasses the entering enthalpy of all in-leakages of water, including tube leaks at hSteamT, relative to a reference enthalpy taken as the saturated liquid enthalpy at the calorimetric temperature; thus: (hSteamT−hf-Cal) as found in Eq.(137).

[0350] Determination of which heat exchanger has a tube leak is accomplished by assigning the tube leak to successive heat exchangers, in repetitive computations involving ERR-CALC and Fuel Iterations, and then examining “key comparative parameters” produced from these computations for deviations from their “reference key comparative parameter”. Reference key comparative parameters are determined from a Fuel Iteration without tube leakage, or as otherwise obtained. If using the System Effect Parameter of computed fuel flow, mAF, such key comparative parameters include: the As-Fired fuel flow, the average fuel water fraction, heating value and the Fuel Consumption Index associated with that heat exchanger. The following weightings of these key comparative parameters has been found useful in determining the lowest deviation among the j heat exchangers when using System Effect Parameter of fuel flow: 33$\begin{array}{cc}\begin{array}{c}{\mathrm{Deviation}}_j=0.02{\&LeftBracketingBar;\left(m_{\mathrm{AF}}-m_{\mathrm{AF}\text{-}\mathrm{PLT}}-\Delta m_{\mathrm{AF}}\right)/\left(m_{\mathrm{AF}\text{-}\mathrm{PLT}}+\Delta m_{\mathrm{AF}}\right)\&RightBracketingBar;}_j+\\ 0.08{\&LeftBracketingBar;\left({\mathrm{WF}}_{\mathrm{H2O}\text{-}\mathrm{Ref}}-{\mathrm{WF}}_{\mathrm{H2O}}\right)/{\mathrm{WF}}_{\mathrm{H2O}\text{-}\mathrm{Ref}}\&RightBracketingBar;}_j+\\ 0.05{\&LeftBracketingBar;\left({\mathrm{HHV}}_{\mathrm{Ref}}-{\mathrm{HHV}}_{\mathrm{AF}}\right)/{\mathrm{HHV}}_{\mathrm{Ref}}\&RightBracketingBar;}_j+\\ 0.85{\&LeftBracketingBar;\left({\mathrm{FCI}}_{\mathrm{Ref}}-\mathrm{FCI}\right)/{\mathrm{FCI}}_{\mathrm{Ref}}\&RightBracketingBar;}_j\end{array}& \left(382\right)\end{array}$

[0351] If using the System Effect Parameter of As-Fired fuel water fraction, WFH2O, such key comparative parameters include: the As-Fired fuel flow, the average fuel water fraction, heating value and the Fuel Consumption Index associated with that heat exchanger. The following weightings of these key comparative parameters has been found useful in determining the lowest deviation among the j heat exchangers when using System Effect Parameter of fuel water: 34$\begin{array}{cc}\begin{array}{c}{\mathrm{Deviation}}_j=0.08{\&LeftBracketingBar;\left(m_{\mathrm{AF}}-m_{\mathrm{AF}\text{-}\mathrm{PLT}}-\Delta m_{\mathrm{AF}}\right)/\left(m_{\mathrm{AF}\text{-}\mathrm{PLT}}+\Delta m_{\mathrm{AF}}\right)\&RightBracketingBar;}_j+\\ 0.02{\&LeftBracketingBar;\left({\mathrm{WF}}_{\mathrm{H2O}\text{-}\mathrm{Ref}}-{\mathrm{WF}}_{\mathrm{H2O}}\right)/{\mathrm{WF}}_{\mathrm{H2O}\text{-}\mathrm{Ref}}\&RightBracketingBar;}_j+\\ 0.05{\&LeftBracketingBar;\left({\mathrm{HHV}}_{\mathrm{Ref}}-{\mathrm{HHV}}_{\mathrm{AF}}\right)/{\mathrm{HHV}}_{\mathrm{Ref}}\&RightBracketingBar;}_j+\\ 0.85{\&LeftBracketingBar;\left({\mathrm{FCI}}_{\mathrm{Ref}}-\mathrm{FCI}\right)/{\mathrm{FCI}}_{\mathrm{Ref}}\&RightBracketingBar;}_j\end{array}& \left(383\right)\end{array}$

[0352] Other key comparative parameters have been studied and have been found as useful as those used in Eqs.(382) & (383), depending on the thermal system and its fuel, and the convergence criteria employed by the minimization techniques. These other key comparative parameters include those in the following list; however, this list is not meant to be exhaustive but representative of the scope and spirit of the present invention:

[0353] Boiler efficiency;

[0354] Dry heating value;

[0355] Specific working fluid energy flow [(BBTC−mTΔh)/mAF];

[0356] Weight fraction of fuel carbon as computed by one of the Input/Loss methods;

[0357] Weight fraction of fuel hydrogen as computed by one of the Input/Loss methods;

[0358] Numerical convergence associated with Minimization Techniques including ANN;

[0359] Computed overall heat transfer coefficient of the assumed leaking heat exchanger;

[0360] Computed log-mean-temperature difference of the assumed leaking heat exchanger;

[0361] Computed cleanliness factor of the assumed leaking heat exchanger; and/or

[0362] Computed relative irreversible thermodynamic loss of the assumed leaking heat exchanger.

[0363] The Preferred Embodiment of this invention is to employ Eqs.(382) and (383). However, the weightings found in these equations are presented to demonstrate the general technique of determining a minimum deviation, and/or otherwise to determine the location of the failed tube. Such weightings may easily change depending on the nature of the fuel being burned, the design of the thermal system, the key comparative parameters employed, and the specific computational procedures employed from one of the various Input/Loss methods. These weightings and the development of a set of key comparative parameters may be determined through simulations of failed tubes, computationally locating the failures within different heat exchangers of the thermal system. The concept and development of Fuel Consumption Indices is described in the following reference: F D Lang, “Fuel Consumption Index for Proper Monitoring of Power Plants—Revisited”, Am. Society of Mech. Engrs., 2002 International Joint Power Generation Conference, Scottsdale, Ariz., IJPGC2002-26097.

[0364] The Preferred Embodiment used to resolve the location of the heat exchanger containing a failed tube is to first establish a set of “reference key comparative parameters” associated with no tube leakage, then to determine the net energy flow to the working fluid assuming the tube leakage is in a particular heat exchanger, repeating such determination until all heat exchangers have been so analyzed. Finally, deviations are formed, similar to Eqs.(382) or (383), from which the lowest deviation will indicate the heat exchanger with the failed tube. Note that if using the procedures taught in Paragraphs 0109 through 0114 to determine the tube failure flow rate, then these same computational methods are used at every evaluation of the net energy flow to the working fluid and associated assumption of where the tube leakage is located. In summary, this invention teaches: to identify a set of heat exchangers descriptive of the thermal system as employed to transfer net energy flow to the working fluid from the combustion gases resulting in a set of identified heat exchangers; to then obtain a set of Operating Parameters applicable to the set of identified heat exchangers; to then determine a set of net energy flows to the working fluid from the combustion gases based on the set of identified heat exchangers, the set of Operating Parameters and the tube leakage flow, each said set of net energy flows descriptive of the thermal system and wherein the tube leakage flow is assigned to a different heat exchanger in each said set; to then determine a reference key comparative parameter for the thermal system resulting in a type of key comparative parameter; to then obtain a set of key comparative parameters associated with each identified heat exchanger and based on the set of net energy flows and the type of key comparative parameter; to then determine a set of deviations between the set of key comparative parameters and the reference key comparative parameter; to then determine an identification of the leaking heat exchanger based on the set of deviations; and to then report to the operator of the thermal system the identification of the leaking heat exchanger such that corrective action may take place.

[0365] Conclusion

[0366] Although the present invention has been described in considerable detail with regard to certain Preferred Embodiments thereof, other embodiments within the scope and spirit of the present invention are possible without departing from the general industrial applicability of the invention. For example, the tube leakage flow rate may be computed by applying the techniques discussed in Paragraphs 0076 through 0108 (i.e., using a single computational pass). The tube leakage flow rate may be determined by use of any System Effect Parameter which indicates an appropriate sensitivity to system effects (with, or in addition to, Lk1, mAF, HHVk3 and WFH2O). The separative analyses process described above may acquire numerous variations; for example: Pass 1 correction factors may be set to unity; or Pass 1 may only employ the Choice Operating Parameter of tube leakage flow rate. The method described is applicable to either higher heating value or lower heating value computations as supported by the teachings in '429. Further, aforementioned descriptions of this invention assume that a steam generator's working fluid is water, however the general procedures of this invention may be applied to any type of working fluid provided that the working fluid is definable at the boundary of the system. Examples of other working fluids are: mixtures of water and organic fluids, organic fluids, liquid metals and so forth. Further, the L Factor concept as used to optimize Choice Operating Parameters may take numerous forms, although two are demonstrated herein, others are discussed in '563, still others may be formed by one skilled in the art based on the teachings herein and in '879 and in '563. Further, the Reduction Efficiency parameter, the Sulfur to Smelt ratio and the Sodium to Carbon ratio in fuel are defined herein as molar ratios whereas they could be re-defined and incorporated within the teachings herein as mass ratios, and/or determined from laboratory analyses or from common industrial experience as mass ratios and then converted to molar ratios. Accordingly, the general theme and scope of the appended claims should not be limited to the descriptions of the Preferred Embodiment disclosed herein.

[0367] Although a Preferred Embodiment of the present invention has been illustrated in THE DRAWINGS and described in considerable detail the foregoing DESCRIPTION OF THE PREFERRED EMBODIMENT, it will be understood that the invention is not limited to the embodiments disclosed, but whose methods are capable of numerous rearrangements, modifications and substitutions without departing from the scope and spirit of the present invention as set forth and defined by the claims herein.

THE DRAWINGS

[0368] FIG. 1 is a schematic representation of a thermal system, particularly a recovery boiler system illustrating use of stoichiometric relationships important in applying this invention. It should be studied in conjunction with combustion stoichiometrics of Eq.(19BL). FIG. 1 depicts a recovery boiler system denoted as 20. In this system 20, a black liquor fuel feed 22 and combustion air 24 are all provided to the upstream side region 26 of the heat exchangers & combustion region 28. Note that this region 28 does not include the air pre-heater 36. In addition, in some types of recovery boilers 20 such as those using fluidized bed combustors, other materials may be injected into region 26, such as a flow of limestone 31 to minimize effluent SO2 by chemically binding sulfur as CaSO4. Other sorbents may be injected to control sulfur, to control other pollutants, and/or to control the combustion process. The black liquor fuel feed 22 contains, in general, combustible organic material, sodium-based inorganic material, water and mineral matter (commonly called inerts, or ash). Inerts are an unburnable component that passes through the system with little physical change, but which are heated and cooled. In the heat exchangers & combustion region 28, the black liquor fuel 22 is burned with the combustion air 24 to form hot products of combustion. Heat from the products of combustion is transferred to a working fluid that enters 23 heat exchangers 32 that are depicted as integral with the heat exchangers & combustion region 28. The heated working fluid 30 is used in a manner appropriate to a working fluid to generate a useful output 33 (for example, in a conventional power plant such useful output may be supplied to a turbine cycle for the production of electrical power). Heat exchangers 32 may consist of a series of heat exchangers as explained in Paragraph 0116. One heat exchanger of the collection of heat exchangers 32 may develop a leakage of its working fluid 27, which leakage mixes with the products of combustion 28. There may also be working fluid leakage 29 into the products of combustion 28 and into region 35, not associated with water in the fuel feed 22, or heat exchanger leakage 27, or moisture in the combustion air 24. Working fluid leakage 29 consists of known flows, or flows which may be otherwise reasonably assumed or determined; and may result from, for example, soot blowing associated with coal-fired systems, or working fluid used to atomize the black liquor fuel 22 before combustion, or used in pollutant control processes located at 35 or 42. The smelt 37 from the combustion of black liquor fuel is removed from region 28. After smelt removal, the remaining products of combustion leave the heat exchangers & combustion region 28 on its downstream region 34, the cooler products of combustion then commonly flow through ducts, region 35, which may contain fly ash removal equipment, passing then to an air pre-heater 36, where a further portion of the combustion gas energy is transferred to an incoming air stream 38, which air then becomes the combustion air 24. The total air delivered to 20 is the incoming air flow 25. In many cases, an air leakage flow 40 enters the flow of the products of combustion as it passes through the air pre-heater 36. The further cooled products of combustion leave the air pre-heater 36 and pass to the Stack 42, then being exhausted to the local environment. Within recovery boiler system 20 the combustion gas path is defined as that region encompassing the flow of products of combustion, said products also termed combustion gases, generally occupying regions 28, 35, the gas side of 36, and 42.

[0369] FIG. 1, given its general system description provided above, is applicable to a wide variety of fossil-fired systems, including recovery boilers, traditional power plants, oil-burning power plants, gas-fired power plants, biomass combustors, fluidized bed combustors, conventional electric power plants, steam generators, package boilers, combustion turbines, and combustion turbines with heat recovery boilers. This list is not meant to be exhaustive, however, and is presented to illustrate some of the areas of applicability of the present invention which encompass any thermal system burning a fossil fuel and which has at least one heat exchanger whose working fluid is being heated by the products of combustion. This invention is applicable to a wide variety of Input/Loss methods, especially its ability to identify the location of the failed tube. If a thermal system is to be characterized quantitatively using one of the Input/Loss methods, then relationships between Choice Operating Parameters to energy flow inputs and outputs as they are affected by tube leakage flow rate 27, may be understood with enhanced accuracy using this invention. This understanding, in turn, permits the operation of the fossil-fired system to be optimized for efficiency, reduction of polluting effluents, prevention of black liquor explosions and reduced down-time given a tube failure and knowledge of its location.

[0370] Within fossil-fired systems, some quantities are readily measured with adequate accuracy, and others may not be measured on-line (in real time) with accuracy sufficient to quantify the operation of the system 20 to the required accuracy to optimize efficiency and/or for the detection of a tube failure and its location. For example, working fluid flow rates, pressures and temperatures may be readily measured with good accuracy by conventional sensors located at defined boundaries such as 23, 30, 25, 33, 42, 22, 29, 31 and 37. Choice Operating Parameters all may, under idea conditions, be directly measured with high accuracy either in real-time or periodically. However, if they are not measured with high accuracy, the ability of Input/Loss methods to quantitatively improve system efficiency may then be compromised. In FIG. 1 quantities which may be (or are) Choice Operating Parameters include: the combustion gas concentrations in the regions 35 and 42 (including CO2, H2O, and O2, termed Λ1B, Λ2B, Λ7B at region 35, and Λ1S, Λ2S, Λ7S at region 42); the indicated combustion air flow 24 (when combined with indicated fuel flow then allows the Air/Fuel ratio to be determined, Λ3, which allows fuel inert fraction to be computed as taught in '994); the ratio of gas concentrations across the air pre-heater, regions 35 and 42 (either the O2 or the CO2 ratio across these regions, preferably the CO2 ratio, thus allowing the Air Pre-Heater Leakage Factor RAct to be determined, Λ4); the concentration of O2 in the combustion air local to the system 25 (termed AAct, or Λ5, allowing ΦAct to be determined); the indicated limestone flow 31 (Λ6); the Sodium to Carbon Ratio in the black liquor fuel 22 determined in real-time or periodically (Λ9); and the Reduction Efficiency term determined in real-time or periodically through analysis of the smelt 37 and the black liquor fuel 22 (Λ10). In addition to these nine, another Choice Operating Parameter is tube leakage flow rate 27 (Λ8), which, in the Preferred Embodiment, is optimized using the fuel's average water content in the fuel (WFH2O-Ref) or using the computed As-Fired fuel flow (mAF); when optimized, the tube leakage flow rate becomes defined consistent with stoichiometrics of Eqs.(19BL) and (200), using in-part Eq.(321B). Refer to Eqs.(211S) through (220). This invention teaches how to correct such measurements or their assumptions if such measurements are not available; or, in the case of a recovery boiler system having a tube leakage, a measurement of such tube leakage flow rate being impossible to obtain in any reasonable manner without using the teachings of this invention.

[0371] FIG. 2 illustrates an important portion of this invention, specifically the general calculational sequences associated with optimizing Choice Operating Parameters and subsequent Fuel Iterations when monitoring a recovery boiler system on-line, i.e., in essentially real time. Box 250 represents general data initialization including establishing Reference Fuel Characteristics, data collection, data organization and routine set-ups of all programs. Box 255 depicts the use of the ERR-CALC program, detailed in FIG. 3, which produces converged Choice Operating Parameters and, given a heat exchanger leak, the tube leakage flow rate. If ERR-CALC is not to be exercised at the same frequency as the Fuel Iterations, Box 255 is bypassed; and, if bypassed, its previously computed correction factors are applied to Λ0-i then employed within the Fuel Iterations. ERR-CALC may be employed, as-needed, for updating correction factors. Box 260 depicts the FUEL program which reduces fuel data from identified multiple sources, prepares a composite fuel, and then prepares an input file for the system simulator EX-FOSS. Reduction of fuel data involves combining the primary (computed) fuel chemistry from a previous iteration, with secondary fuels which have constant and known chemistries, producing a composite fuel. Box 270 is system data supplied to the process as on-line (or real time) including at least the following Operating Parameters (refer to Paragraph 0031 for details): working fluid pressures, temperatures and flows, air psychrometrics, useful system output and other related data Box 280 depicts the system simulator EX-FOSS which, given specification of a composite fuel from FUEL, inputs from Box 270, routine set-up data and converged Choice Operating Parameters from Box 255 (including possible tube leakage flow rate), produces the following: boiler efficiency using the methods of '429 as modified by the present invention, As-Fired fuel flow (mAF) using Eq.(330A) or (330B), complete effluent concentrations, system heat rate using Eq.(331A), (331B), (332A) or (332B), effluent flow, emission rates of all effluents including the common pollutants, and other thermal performance parameters including, for example, energy flow to the working fluid given a tube leak (BBTC−mTΔh), and the Firing Correction (HBC). The determination of many of these parameters is taught herein and in '994. Box 285 depicts the HEATRATE program within which, given the converged Choice Operating Parameters (including possible tube leakage flow rate), produces fuel chemistry, L Factors and fuel heating value for both the composite fuel (as either higher or lower heating values), and, given the known compositions of secondary fuels, the composition of the primary (unknown) fuel. Box 285 also depicts the computations within HEATRATE which lead to identification of which heat exchanger has the tube leakage; refer to Paragraphs 0115 through 0119. Designation 287 tests for convergence of the process based on composite fuel moles (x), certain effluents such as CO2 and H2O, heating value and computed fuel water fraction; if the convergence criteria is not met the process continues to iterate from Box 260. In general, convergences lie within 0.5×10−4 percent of the computed As-Fired fuel moles. Note that the iterations encompassing 260, 270, 280, 285 and 287 define what is meant by the term “Fuel Iterations”. In summary, fuel Iterations are the iterative calculations between EX-FOSS, as input with known fuel chemistry and heating value from a previous iteration but with unknown effluents (to be computed by EX-FOSS, except for effluent O2 which is input), and HEATRATE as input with known effluents (i.e., the corrected Choice Operating Parameters) but with unknown fuel chemistry and heating value (to be computed by HEATRATE).

[0372] Designation 290 of FIG. 2 illustrates an important teaching of this invention as it represents a decision to either: 1) continue towards reporting results and quitting if either the Tube Failure Model (i.e., this invention) is not invoked, or is invoked and all computations and processes have been completed; or 2) if the Tube Failure Model is invoked and the location of the failure is desired, to then continue the computations, starting afresh with Box 255 but with every pass through Box 255 to assume the tube leakage occurs in a different heat exchanger. If the tube failure flow rate and its location are both desired, then summary analysis of these computations (leading to identification of the leakage's location) is performed within the HEATRATE program, Box 285, after all heat exchangers have been analyzed. Typically, monitoring cycles are scheduled for every 2 to 3 minutes using updated data based on 15 minute running averages. Once converged and all computations have been completed, Box 294 produces reportable results from the EX-FOSS and HEATRATE programs, including system heat rate, tube leakage flow and its location, and other thermal performance parameters which are influenced by the failed tube including: Second Law analysis of the thermal system (producing Fuel Consumption Indices), fuel flow, total effluent flow, emission rates, other output and reports to system operators as to what corrective actions may take place; said reports also being provided to regulatory authorities as requested or required. Box 298 of FIG. 2 is to quit.

[0373] FIG. 3 illustrates another important portion of this invention, specifically the organization of the ERR-CALC program used to determine correction factors to the initial Choice Operating Parameters, including the computation of tube leakage flow rate if warranted. In FIG. 3 Box 310 depicts the start of the program which invokes data collection and routine program set-up associated with ERR-CALC. Routine program set-up associated with ERR-CALC includes a user option as to whether the Tube Failure Model will be invoked, or not, for a given monitoring cycle (also see Box 250 of FIG. 2). Box 320 depicts initializations of data including organization of data arrays associated with selected Choice Operating and System Effect Parameters, and determination of scaling factors, Si, and pre-scaling factors, si. Box 330 depicts the selection, collection and processing of general. input data associated with the minimization techniques, principally the selection of Choice Operating Parameters, System Effects Parameters, Reference Systems Effects Parameters and routine inputs, options and convergence criteria to the minimization techniques as are known to those skilled in the art using these techniques (such inputs and criteria are presented in cited references in Paragraphs 0104, 0105, and 0107).

[0374] Within FIG. 3, if the Tube Failure Model is optioned, Boxes 340, 350 and 360 may be executed three times for the Preferred Embodiment; these are termed Passes 0, 1 and 2. The initial execution of Boxes 340, 350 and 360 as a series is termed Pass 0. Box 340 depicts application of the minimization techniques as herein discussed, including evaluation of an objective function resulting in optimizing the selected Choice Operating Parameters. Box 350 depicts the use of a simulation principally of the HEATRATE program within ERR-CALC by which the computing time required for the supporting computations required for Box 340 are greatly reduced; refer to Paragraph 0103 for details. For Simulated Annealing, ANN and other such exhaustive procedures, having importance in implementing the Preferred Embodiment of this invention given a tube failure, Box 350 is typically caused to be executed from Box 340 thousands of times for each Pass. Inputs to Box 350 are principally Choice Operating Parameters. Output from Box 350 to Box 340 being principally System Effect Parameters from which the objective function is then evaluated. After convergence of the minimization techniques of Box 340 and thus resolution of converged Choice Operating Parameters, Box 360 depicts the calculation of correction factors associated with the selected Choice Operating Parameters. Box 360 also includes the production of appropriate warning messages associated with the ERR-CALC computations; for example: non-convergence, computational failures, the automatic switching to alternative minimization techniques, and the like. If the monitoring cycle is processing Pass 0 and the Tube Failure Model is not optioned, the program quits via Box 470.

[0375] If the Tube Failure Model is optioned, logic flows as discussed above, from Box 340 through Box 360 resulting in a converged solution (but without an assigned tube leakage). However, within Box 360 logic is then invoked which initiates Pass 1 or Pass 2, either proceeding to Box 430. For Pass 1, Box 430 tests for a trip mechanism of TABLE 1A and TABLE 1B, and if found that a tube failure is possible, then: re-sets Λ1S or Λ1B, and Λ2S or Λ2B correction factors to historically based factors (or to unity); selects the System Effect Parameter for fuel water (WFH2O) or for fuel flow (mAF) whose difference with reference values is reduced by optimizing Choice Operating Parameter for tube leakage flow (Λ8); selects other Choice Operating Parameters and associated System Effect Parameters, excluding Λ1S, Λ1B, Λ2S and Λ2B; selects Simulated Annealing as the minimization technique; employs the computed tube leakage flow (mT); and passes control to Box 340. Convergence via Boxes 340 and 350 results in a stoichiometric consistency with the computed tube leakage flow rate. After Pass 1, Box 430 logic then initiates Pass 2 proceeding to Box 430. For Pass 2, Box 430 then: re-sets all correction factors to unity; selects the same System Effect Parameters and Choice Operating Parameters established in Box 330 (excluding Λ2S and Λ2B), plus System Effect Parameter for fuel water or fuel flow whose differences with their reference values is reduced through optimization; selects the Simulated Annealing technique;. employs the computed tube leakage flow rate of Pass 1; and passes control to Box 340. Convergence via Boxes 340 and 350 results in correction factors for all selected Choice Operating Parameters which, in combination with the determined tube leakage flow rate via Eq.(321B), achieves stoichiometric consistency as stated by Eqs.(19BL) and (200). Within Box 360, logic produces appropriate warning messages that a tube failure has occurred (assuming the tube leakage flow rate is greater than zero), its mass flow rate, the associated Tube Failure Mechanism, correction factors, and other routine messages. The execution of ERR-CALC is then terminated with Box 470.

[0376] FIG. 4 illustrates a plot of both actual and computed water flow entering a combustion space of a large commercial steam generator under test conditions. The actual flow is measured soot blowing flow which was declared an unknown to the Calculational Engine operating The Input/Loss Method as disclosed herein (bZ=0.0). This test emulated a tube failure. The observed sensitivity in FIG. 4 is less than 0.1% of feedwater flow thus demonstrating the practicality of this invention when applied to recovery boilers. FIG. 4 is an outstanding example of the sensitivity and accuracy of this invention, demonstrating many of this invention's objectives.

[0377] Symbols within equations may have been italicized pursuant to Patent Office publication practices. As used in FIG. 1, FIG. 2 and FIG. 3, and throughout the DESCRIPTION OF THE PREFERRED EMBODIMENT, mathematical symbols typed in italics and mathematical symbols typed in non-italics have the same meaning when taken in context. Understanding context is afforded through the use of subscripts and/or through the normal flow of mathematical development. Thus, for example, the symbols: AAct, RAct, F({right arrow over (x)}), ƒ( ), Ci, CP, hj, h, J0, J1, JMech, L′Fuel, LCO2, mAF, x, x1, {right arrow over (x)}, AF, HR defined in Paragraphs 0027 through 0030 and elsewhere herein, have the same meaning as, respectively: AAct, RAct, F({right arrow over (x)}), f( ), Ci, CP, hj, h, J0, J1, JMech, L′Fuel, LCO2, mAF, x, xi, {right arrow over (x)}, AF, HR. The use of italic symbols is used for writing style. For example, the lower-case ƒ (or f) describes a functional relationship, defined in Paragraph 0028, and as used in Paragraphs 0033, 0059, 0064, 0069, 0077 and 0085. For example, as defined in Paragraph 0027, the letter J (or J) with the subscript “Act” describes the “Total effluent water at the system's boundary (j+bAβ); moles/base”; as opposed to as the letter J (or J) with the subscript “0” or “1”, defined in Paragraph 0028, describing the Bessel function of the first kind of order zero or one, and used in Paragraphs 0083, 0085 and 0098; as opposed to JMech which is defined in Paragraph 0075 as the mechanical equivalent of heat. For example, the letter “x” (and related terms xtheor, xDRY-theor and xMAF-theor) is used to describe the moles of fuel per base, defined in Paragraph 0027; whereas the terms xi and {right arrow over (x)} are used to describe generic independent variables associated with multidimensional minimization techniques, defined in Paragraph 0028. The moles of fuel is used principally throughout Paragraphs 0044 through 0075; while the terms xi and {right arrow over (x)} are exclusively used in Paragraphs 0081 through 0085, 0098 and 0102. For example, Ci is used as a correction factor to Choice Operating Parameters (e.g., C2S), defined in paragraph 0028 and used in paragraphs 0079, 0090, 0099, 0106, 0109, 0112 and 0114; whereas CP is heat capacity defined in paragraph 0029 and used in paragraphs 0071 and 0073. Note that correlation constants denoted, for example as A3, B15, B3, C5, etc. and used for the mathematical development of Eqs.(42), (43), (48), (49) and (50) are discussed in paragraphs 0056, 0057, 0058, 0077 and 0078. For example, the term “h” is the moles of effluent nitrogen without air leakage per base as defined in Eq.(l9BL) and exclusively used in system stoichiometrics of paragraphs 0044 through 0062; whereas hj is a specific enthalpy term defined in paragraph 0029 used in energy equations following paragraph 0063.

Claims

1. A method for quantifying the operation of a recovery boiler burning black liquor fuel bearing sodium compounds through knowledge of when its heat exchanger leaks working fluid into the combustion gas path producing a tube leakage, the method for quantifying the operation comprising the steps of: monitoring the recovery boiler burning black liquor fuel bearing sodium compounds by one of the Input/Loss methods, developing a mathematical model of the combustion process incorporating terms commonly associated with black liquor fuel combustion including sodium compounds and terms associated with sources of working fluid flows into the combustion gas path including tube leakage resulting in a stoichiometric model of the combustion process, and determining a tube leakage in moles based on the stoichiometric model of the combustion process.

2. The method according to claim 1 further comprising the steps, after determining, of: obtaining a molecular weight of the fossil fuel, obtaining a molecular weight of the working fluid, obtaining a fuel flow rate of the recovery boiler, determining a tube leakage mass flow rate based on the tube leakage in moles, the molecular weight of the fossil fuel, the molecular weight of the working fluid, the fuel flow rate, and the stoichiometric model of the combustion process, and reporting the tube leakage mass flow rate such that corrective action may take place.

3. The method of claim 1, wherein the step of determining the tube leakage in moles comprises the steps of: forming a hydrogen stoichiometric balance based on the stoichiometric model of the combustion process using a molar base, and solving the hydrogen stoichiometric balance for the tube leakage in moles. /// /// /// ///

4. The method according to claim 2 further comprising the steps, after reporting, of: identifying a set of heat exchangers descriptive of the recovery boiler as employed to transfer net energy flow to the working fluid from the combustion gases resulting in a set of identified heat exchangers, obtaining a set of Operating Parameters applicable to the set of identified heat exchangers, analyzing a set of net energy flows to the working fluid from the combustion gases based on the set of identified heat exchangers, the set of Operating Parameters and the tube leakage flow rate, each analyzed set descriptive of the recovery boiler and wherein each analyzed set the tube leakage flow rate is assigned to a different heat exchanger, resulting in an analyzed set of heat exchangers, determining a reference key comparative parameter for the recovery boiler, obtaining a set of key comparative parameters associated with each identified heat exchanger, applicable with the reference key comparative parameter, and based on the analyzed set of heat exchangers, determining a set of deviations between the set of key comparative parameters and the reference key comparative parameter, determining an identification of the leaking heat exchanger based on the set of deviations, and reporting to the operator of the recovery boiler the identification of the leaking heat exchanger such that corrective action may take place.

5. A method for quantifying the operation of a recovery boiler burning black liquor fuel bearing sodium compounds when being monitored by one of the Input/Loss methods through knowledge of when its heat exchanger leaks working fluid into the combustion gas path producing a tube leakage, the method for quantifying the operation comprising the steps of: developing a mathematical model of the combustion process incorporating terms commonly associated with the combustion of black liquor fuel including sodium compounds and terms associated with sources of working fluid flows into the combustion gas path including tube leakage resulting in a stoichiometric model of the combustion process, selecting a set of minimization techniques applicable to the recovery boiler burning black liquor fuel, and a set of routine inputs and convergence criteria to the minimization techniques, selecting a Choice Operating Parameter of tube leakage flow rate, selecting a set of routine Choice Operating Parameters, determining a set of System Effect Parameters applicable to the recovery boiler burning black liquor fuel whose functionalities are sensitive to tube leakage flow rate, determining a set of Reference System Effect Parameters applicable to the set of System Effect Parameters, determining an objective function applicable to the recovery boiler's stoichiometric situation which uses the Choice Operating Parameter of tube leakage flow rate, the set of routine Choice Operating Parameters, the set of System Effect Parameters, and the set of Reference System Effect Parameters, optimizing the Choice Operating Parameter of tube leakage flow rate and the set of routine Choice Operating Parameters based on the mathematical model of the combustion process, the set of minimization techniques and the objective function such that convergence criteria is met resulting in a set of converged Choice Operating Parameters including the tube leakage flow rate, and reporting the tube leakage flow rate such that corrective action may take place.

6. The method according to claim 5 further comprising the steps, after reporting, of: determining a set of Reference Fuel Characteristics, determining the fuel chemistry of the black liquor fuel being combusted in the recovery boiler using one of the Input/Loss methods, the mathematical model of the combustion process, the set of converged Choice Operating Parameters, and the set of Reference Fuel Characteristics, determining a fuel heating value of the system based on the fuel chemistry and the set of Reference Fuel Characteristics, obtaining a set of Operating Parameters, determining a Firing Correction base on the set of Operating Parameters, and determining a high accuracy boiler efficiency of the recovery boiler independent of fuel flow based on the set of converged Choice Operating Parameters including the tube leakage flow rate, the fuel chemistry, the fuel heating value, the Firing Correction and the set of Operating Parameters.

7. The method according to claim 6 further comprising the steps, after determining the high accuracy boiler efficiency, of: determining an energy flow to the working fluid of the recovery boiler based on the set of Operating Parameters as influenced by the tube leakage flow rate, determining a fuel flow of the fossil fuel being combusted using the energy flow to the working fluid, the fuel heating value, the Firing Correction and the high accuracy boiler efficiency, and reporting the fuel flow as influenced by the tube leakage flow rate.

8. The method according to claim 7 further comprising the steps, after reporting, of: determining a useful output from the recovery boiler, determining a system efficiency using the fuel flow, the fuel heating value, the Firing Correction and the useful output from the recovery boiler, and reporting the system efficiency as influenced by the tube leakage flow rate.

9. The method according to claim 7 further comprising the steps, after reporting, of: determining a useful output from the recovery boiler, determining a system efficiency using the energy flow to the working fluid, the high accuracy boiler efficiency and the useful output from the recovery boiler, and reporting the system efficiency as influenced by the tube leakage flow rate.

10. The method of claim 5, wherein the step of selecting the set of minimization techniques applicable to the recovery boiler burning black liquor fuel comprises a step of: incorporating a BFGS technique applicable to the recovery boiler and its fuel.

11. The method of claim 5, wherein the step of selecting the set of minimization techniques applicable to the recovery boiler burning black liquor fuel comprises a step of: incorporating a Simulated Annealing technique applicable to the recovery boiler and its fuel.

12. The method of claim 5, wherein the step of selecting the set of minimization techniques applicable to the recovery boiler burning black liquor fuel comprises a step of: incorporating a neural network technique applicable to the recovery boiler burning black liquor fuel.

13. The method of claim 5, wherein the step of selecting the set of minimization techniques applicable to the recovery boiler burning black liquor fuel comprises a step of: incorporating a Neugents technology applicable to the recovery boiler burning black liquor fuel.

14. The method of claim 5, wherein the step of selecting the set of minimization techniques applicable to the recovery boiler burning black liquor fuel comprises a step of: incorporating a Pegasus Technology applicable to the recovery boiler burning black liquor fuel.

15. The method of claim 5, wherein the step of selecting the set of minimization techniques applicable to the recovery boiler burning black liquor fuel comprises a step of: incorporating a NeuCo, Inc. technology.

16. The method of claim 5, wherein the step of selecting the set of routine Choice Operating Parameters comprises a step of: determining a set of scaling factors for the set of routine Choice Operating Parameters resulting in the set of routine Choice Operating Parameters whose values are scaled.

17. A method for quantifying the operation of a recovery boiler burning black liquor fuel when being monitored on-line by The Input/Loss Method through knowledge of when its heat exchanger leaks working fluid into the combustion gas path producing a tube leakage flow rate, the method comprising the steps of: before on-line operation, developing a mathematical model of the combustion process incorporating terms commonly associated with the combustion of black liquor fuel and terms associated with sources of working fluid flows into the combustion gas path including tube leakage resulting in a stoichiometric model of the combustion process, and selecting a set of minimization techniques applicable to the recovery boiler burning black liquor fuel, and a set of routine inputs and convergence criteria to the minimization techniques; and thereafter performing a Pass 0 monitoring cycle while operating on-line which determines a trip mechanism associated with a tube leakage flow rate, comprising the steps of: assuming a tube leakage flow rate is zero, selecting a set of routine Choice Operating Parameters, determining a set of System Effect Parameters applicable to the recovery boiler burning black liquor fuel whose functionalities effect the determination of system efficiency, and determining applicable Reference System Effect Parameters, performing a Pass 0 monitoring cycle using The Input/Loss Method resulting in a set of fuel concentrations and correction factors to the set of routine Choice Operating Parameters, testing the set of fuel concentrations and correction factors against a set of corresponding limits resulting in a trip mechanism indicating the stoichiometric reason how a heat exchanger may leak a tube leakage flow rate into the combustion gas path, and assuming the trip mechanism indicates a tube leakage flow rate is possible; and thereafter performing a Pass 1 monitoring cycle while operating on-line which determines a tube leakage flow rate, comprising the steps of: multiplying a value of uncorrected effluent water concentration by a reference correction factor resulting in a constant effluent water concentration, multiplying a value of uncorrected effluent carbon dioxide concentration by a reference correction factor resulting in a constant effluent carbon dioxide concentration, selecting a set of routine Choice Operating Parameters including tube leakage flow rate but excluding effluent water and carbon dioxide, determining a set of System Effect Parameters applicable to the recovery boiler burning black liquor fuel whose functionalities effect the determination of system efficiency and are sensitive to tube leakage flow rate, and determining applicable Reference System Effect Parameters, performing a Pass 1 monitoring cycle using The Input/Loss Method resulting in a set of converged Choice Operating Parameters including the tube leakage flow rate, and thereafter; reporting the tube leakage flow rate such that corrective action may take place.

18. The method according to claim 5 further comprising the steps, after reporting, of: identifying a set of heat exchangers descriptive of the recovery boiler as employed to transfer net energy flow to the working fluid from the combustion gases resulting in a set of identified heat exchangers, obtaining a set of Operating Parameters applicable to the set of identified heat exchangers, analyzing a set of net energy flows to the working fluid from the combustion gases based on the set of identified heat exchangers, the set of Operating Parameters and the tube leakage flow rate, each analyzed set descriptive of the recovery boiler and wherein each analyzed set the tube leakage flow rate is assigned to a different heat exchanger, resulting in an analyzed set of heat exchangers, determining a reference key comparative parameter for the recovery boiler, obtaining a set of key comparative parameters associated with each identified heat exchanger, applicable with the reference key comparative parameter, and based on the analyzed set of heat exchangers, determining a set of deviations between the set of key comparative parameters and the reference key comparative parameter, determining an identification of the leaking heat exchanger based on the set of deviations, and reporting to the operator of the recovery boiler the identification of the leaking heat exchanger such that corrective action may take place.

19. The method of claim 18, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a fuel flow as the reference key comparative parameter for the recovery boiler.

20. The method of claim 18, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a fuel water fraction as the reference key comparative parameter.

21. The method of claim 18, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a heating value as the reference key comparative parameter.

22. The method of claim 18, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a Fuel Consumption Index for each heat exchanger as the reference key comparative parameter for the recovery boiler.

23. A method for quantifying the operation of a recovery boiler burning black liquor fuel when being monitored by one of the Input/Loss methods through knowledge of a stoichiometric mechanism of how a heat exchanger could be leaking a tube leakage flow rate into the combustion gas path, the method for quantifying the operation comprising the steps of: developing a mathematical model of the combustion process incorporating terms commonly associated with fossil fuel combustion and terms associated with sources of working fluid flows into the combustion gas path including tube leakage, obtaining a set of Choice Operating Parameters, obtaining a set of Reference Fuel Characteristics, obtaining a fuel chemistry of the fuel being combusted by the recovery boiler using one of the Input/Loss methods, the mathematical model of the combustion process, the set of Choice Operating Parameters, and the set of Reference Fuel Characteristics, said fuel chemistry resulting in a set of fuel concentrations, establishing a set of concentration limits for each fuel constituent based on Reference Fuel Characteristics, testing the set of fuel concentrations against the set of concentration limits resulting in a trip mechanism indicating the stoichiometric reason how a heat exchanger leaks a tube leakage flow rate into the combustion gas path, and reporting the trip mechanism to the operator of the recovery boiler.

24. A method for quantifying the operation of a recovery boiler burning black liquor fuel when being monitored by one of the Input/Loss methods through knowledge of a stoichiometric mechanism of how a heat exchanger could be leaking a tube leakage flow rate into the combustion gas path, the method for quantifying the operation comprising the steps of: developing a mathematical model of the combustion process incorporating terms commonly associated with fossil fuel combustion and terms associated with sources of working fluid flows into the combustion gas path including tube leakage, selecting a set of minimization techniques applicable to the recovery boiler burning black liquor fuel, processing a set of routine inputs and convergence criteria to the minimization techniques, assuming a tube leakage flow rate is zero, selecting a set of routine Choice Operating Parameters, determining a set of System Effect Parameters applicable to the recovery boiler burning black liquor fuel whose functionalities effect the determination of system efficiency, determining a set of Reference System Effect Parameters applicable to the set of System Effect Parameters, determining an objective function applicable to the recovery boiler, the set of routine Choice Operating Parameters, the set of System Effect Parameters and the set of Reference System Effect Parameters, optimizing the set of routine Choice Operating Parameters based on the mathematical model of the combustion process, the set of minimization techniques, and the objective function such that convergence is met resulting in a set of converged Choice Operating Parameters, determining a fuel chemistry of the fuel being combusted by the recovery boiler using one of the Input/Loss methods, the mathematical model of the combustion process, the set of converged Choice Operating Parameters, and Reference Fuel Characteristics resulting in a fuel elemental composition, a fuel ash fraction and a fuel water fraction said composition and fractions resulting in a set of fuel concentrations, establishing a set of concentration limits for the set of fuel concentrations based on Reference Fuel Characteristics, testing the set of fuel concentrations against the set of concentration limits resulting in a trip mechanism indicating the stoichiometric reason how a heat exchanger leaks a tube leakage flow rate into the combustion gas path, and reporting the trip mechanism to the operator of the recovery boiler.

25. The method of claim 24, wherein the step of establishing the set of concentration limits for the set of fuel concentrations based on Reference Fuel Characteristics and the step of testing the set of fuel concentrations against the concentration limits, comprises the steps of: determining a set of correction factors to Choice Operating Parameters using their initial and converged values, and establishing a set of correction factor limits for the selected Choice Operating Parameters, and testing the set of correction factors against the set of correction factor limits resulting in a trip mechanism indicating the stoichiometric reason how a heat exchanger leaks a tube leakage flow rate into the combustion gas path.

26. A method for quantifying the operation of a recovery boiler burning black liquor fuel in a combustion process through knowledge of when one of its heat exchangers, whose tubes contain working fluid heated by products of combustion, has a tube leak of working fluid mixing with the products of combustion, the method for quantifying the operation comprising the steps of: selecting a neural network technology applicable to the recovery boiler, selecting a set of routine inputs and database for the neural network technology, selecting a set of Choice Operating Parameters including tube leakage flow rate, developing a mathematical model of the combustion process incorporating terms commonly associated with a combustion process and terms associated with sources of working fluid mixing with the products of combustion including tube leakage, obtaining a set of Reference System Effect Parameters associated with the recovery boiler, calculating a set of System Effect Parameters based on the mathematical model of the combustion process, optimizing the set of Choice Operating Parameters including tube leakage flow rate using the neural network technology, and the set of routine inputs and database, the set of Reference System Effect Parameters and the set of System Effect Parameters, such that convergence is met resulting in a set of converged Choice Operating Parameters including a tube leakage flow rate, and reporting the tube leakage flow rate such that corrective action may take place.

27. The method of claim 26, wherein the step of selecting the neural network technology applicable to the recovery boiler burning black liquor fuel, comprises a step of: selecting a Neugents technology applicable to the recovery boiler burning black liquor fuel.

28. The method of claim 26, wherein the step of selecting the neural network technology applicable to the recovery boiler burning black liquor fuel, comprises a step of: selecting a Pegasus Technology applicable to the recovery boiler burning black liquor fuel.

29. The method of claim 26, wherein the step of selecting the neural network technology applicable to the recovery boiler burning black liquor fuel, comprises a step of: selecting a NeuCo, Inc. technology applicable to the recovery boiler burning black liquor fuel.

30. A method for quantifying the operation of a recovery boiler burning black liquor fuel when being monitored by one of the Input/Loss methods coincident with one of its heat exchangers leaking its working fluid into the combustion gas path producing a tube leakage flow, the method for quantifying the operation by identification of the leaking heat exchanger comprising the steps of: identifying a set of heat exchangers descriptive of the recovery boiler as employed to transfer net energy flow to the working fluid from the combustion gases resulting in a set of identified heat exchangers, obtaining a set of Operating Parameters applicable to the set of identified heat exchangers, analyzing a set of net energy flows to the working fluid from the combustion gases based on the set of identified heat exchangers, the set of Operating Parameters and the tube leakage flow rate, each analyzed set descriptive of the recovery boiler and wherein each analyzed set the tube leakage flow rate is assigned to a different heat exchanger, resulting in an analyzed set of heat exchangers, determining a reference key comparative parameter for the recovery boiler, obtaining a set of key comparative parameters associated with each identified heat exchanger, applicable with the reference key comparative parameter, and based on the analyzed set of heat exchangers, determining a set of deviations between the set of key comparative parameters and the reference key comparative parameter, determining an identification of the leaking heat exchanger based on the set of deviations, and reporting to the operator of the recovery boiler the identification of the leaking heat exchanger such that corrective action may take place.

31. The method of claim 30, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a fuel flow as the reference key comparative parameter for the recovery boiler.

32. The method of claim 30, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a fuel water fraction as the reference key comparative parameter for the recovery boiler.

33. The method of claim 30, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a heating value as the reference key comparative parameter for the recovery boiler.

34. The method of claim 30, wherein the step of determining the reference key comparative parameter for the recovery boiler, comprises a step of: selecting a computed cleanliness factor for each heat exchanger as the reference key comparative parameter for the recovery boiler.

35. A method for quantifying the operation of a recovery boiler burning a fossil fuel in a combustion process through knowledge of when one of its heat exchangers, whose tubes contain working fluid heated by products of combustion, has a tube leak of working fluid mixing with the products of combustion, the method for quantifying the operation comprising the steps of: monitoring the recovery boiler using one of the Input/Loss methods, developing a mathematical model of the combustion process incorporating terms commonly associated with the combustion process and terms associated with sources of working fluid mixing with the products of combustion including tube leakage, determining a tube leakage based on the mathematical model of the combustion process, and reporting the tube leakage such that corrective action may take place.

36. The method of claim 35, wherein the step of developing a mathematical model of the combustion process comprises the steps of: forming a hydrogen stoichiometric balance of the combustion process including terms associated with sources of working fluid mixing with the combustion products including tube leakage, and solving the hydrogen stoichiometric balance for the tube leakage.

37. The method of claim 35, wherein the step of monitoring the recovery boiler using one of the Input/Loss methods, comprises the step of: monitoring the recovery boiler using The Input/Loss Method.

38. The method of claim 35, wherein the step of monitoring the recovery boiler using one of the Input/Loss methods, comprises the step of: determining a fuel chemistry based on one of the Input/Loss methods.

39. The method of claim 35, wherein the step of monitoring the recovery boiler using one of the Input/Loss methods, comprises the steps of: determining a fuel heating value based on one of the Input/Loss methods.

40. The method of claim 35 further comprising the steps, after reporting, of: identifying a set of heat exchangers descriptive of the recovery boiler resulting in a set of identified heat exchangers, obtaining a set of Operating Parameters applicable to the set of identified heat exchangers resulting in a set of heat exchanger data sufficient to determine net energy flow to the working fluid from the products of combustion for each heat exchanger in the set of identified heat exchangers, calculating a net energy flow to the working fluid of the recovery boiler as many times as there are heat exchangers in the set of identified heat exchangers, wherein each calculation of net energy flow includes all heat exchangers in the set of identified heat exchangers, wherein for each calculation of net energy flow the tube leakage is assigned to a different heat exchanger, resulting in a set of analyzed heat exchangers based on the set of heat exchanger data, determining a set of reference key comparative parameters, obtaining a set of key comparative parameters associated with the set of identified heat exchangers applicable with the set of reference key comparative parameters, and based on the set of analyzed heat exchangers, determining a set of deviations between the set of key comparative parameters and the set of reference key comparative parameters, identifying a location of the heat exchanger within the recovery boiler having the tube leak based on the set of deviations, and reporting to the operator of the recovery boiler the location of the heat exchanger within the recovery boiler having the tube leak such that corrective action may take place.

41. A method for quantifying the operation of a recovery boiler burning a black liquor fuel in a combustion process through knowledge of when one of its heat exchangers, whose tubes contain working fluid heated by products of combustion, has a tube leak of working fluid mixing with the products of combustion, the method for quantifying the operation comprising: determining a location of the heat exchanger within the recovery boiler with the tube leak based on the working fluid's energy flow by assigning the tube leak to different heat exchangers.

42. The method of claim 41 further comprising: obtaining a heating value of the black liquor fuel, obtaining a Firing Correction applicable to the recovery boiler, obtaining a high accuracy boiler efficiency, and determining a calculated fuel flow based on the working fluid's energy flow effected by the tube leak of working fluid, the location of the heat exchanger within the recovery boiler with the tube leak, the high accuracy boiler efficiency, the fossil fuel heating value, and the Firing Correction.

43. The method of claim 41 further comprising: obtaining a high accuracy boiler efficiency, obtaining a useful output produced from the recovery boiler, determining a calculated system efficiency of the recovery boiler based on the working fluid's energy flow effected by the tube leak of working fluid, the location of the heat exchanger within the recovery boiler with the tube leak, the high accuracy boiler efficiency and the useful output produced from the recovery boiler.

44. The method of claim 42, wherein the step of obtaining the high accuracy boiler efficiency comprises: using the black liquor fuel's calorimetric temperature, established when determining the fuel's heating value, as the thermodynamic reference energy level for an Enthalpy of Products term, as the thermodynamic reference energy level for an Enthalpy of Reactants term, and also as the thermodynamic reference energy level for a Firing Correction term evaluated independent of a fuel flow and an effluent flow, said terms comprising the major terms of the high accuracy boiler efficiency.

45. The method of claim 43, wherein the step of obtaining the high accuracy boiler efficiency comprises: using the black liquor fuel's calorimetric temperature, established when determining the fuel's heating value, as the thermodynamic reference energy level for an Enthalpy of Products term, as the thermodynamic reference energy level for an Enthalpy of Reactants term, and also as the thermodynamic reference energy level for a Firing Correction term evaluated independent of a fuel flow and an effluent flow, said terms comprising the major terms of the high accuracy boiler efficiency.

46. The method of claim 42 wherein the step of obtaining the heating value of the black liquor fuel comprises: obtaining a higher heating value of the fuel, obtaining a lower heating value of the fuel; and wherein the step of obtaining the high accuracy boiler efficiency comprises: obtaining a higher heating value high accuracy boiler efficiency based on the higher heating value of the fuel, obtaining a lower heating value high accuracy boiler efficiency based on the lower heating value of the fuel; and wherein the step of determining the calculated fuel flow comprises: demonstrating that a computed fuel flow based on the higher heating value high accuracy boiler efficiency and a computed fuel flow based on the lower heating value high accuracy boiler efficiency are comparable.

47. A method for quantifying the operation of a recovery boiler burning black liquor fuel having a heat exchangers/combustion region producing combustion products, the method comprising the steps of: before on-line operation, the steps of obtaining reference fuel characteristics including at least intrinsic chemical relationships between fuel hydrogen and fuel carbon, and between fuel sodium and fuel carbon, and developing explicit mathematical models of the combustion process involving at least stoichiometric balances; and thereafter operating on-line, the step of operating on-line including the steps of measuring a set of measurable operating parameters, including at least effluent concentrations of O2 and CO2, these measurements being made at a location downstream of the heat exchangers/combustion region of the recovery boiler, obtaining an effluent concentration of H2O, if reference fuel characteristics indicate fuel water is not predictable, as an obtained effluent H2O, obtaining an ambient concentration of O2, obtaining an air pre-heater leakage factor, and calculating a complete As-Fired fuel chemistry, including fuel water and fuel inerts, as a function of the reference fuel characteristics, explicit mathematical models of the combustion process, the set of measurable operating parameters, the obtained effluent H2O, the ambient concentration of O2, and the air pre-heater leakage factor.

48. The method of claim 47, wherein the step of operating on-line includes the additional step after calculating the complete As-Fired fuel chemistry, of calculating consistent moisture-ash-free, dry and As-Fired fuel heating values as a function of the complete As-Fired fuel chemistry and the reference fuel characteristics.

49. The method of claim 48, including, after the step of calculating consistent moisture-ash-free, dry and As-Fired heating values, the additional steps of obtaining a System Effect Parameter associated with the recovery boiler and its fuel, completing a multidimensional minimization analysis employing the System Effect Parameter to minimize the error associated with at least one of the measured effluent CO2, the obtained effluent H2O, the ambient concentration of O2 and the air pre-heater leakage factor, calculating and applying for subsequent on-line analysis correction factors to the measured effluent CO2, obtained effluent H2O, the ambient concentration of O2, and the air pre-heater leakage factor.

50. A method for quantifying the operation of a recovery boiler burning black liquor fuel in which a fossil fuel is supplied at a flow rate to a heat exchangers/combustion region and combusted to produce hot combustion gases, which heats a working fluid then exits through an exhaust stack, the method comprising the following steps: performing an off-line operation comprising the steps of obtaining reference fuel characteristics, obtaining current measurements of the system's operating parameters, and performing an on-line operation comprising the steps of measuring the useful output of the system, obtaining fuel data and characteristics, the step of obtaining fuel data including the step of obtaining composite fuel concentrations and composite heating value, if multiple fuels are used, introducing fuel concentrations and heating values to a mathematical model of the recovery boiler, obtaining routine systems operational parameters, obtaining values of the effluents O2, CO2, H2O and SO2, obtaining the ambient concentration of O2, obtaining air pre-heater leakage and dilution factors, computing molar moisture-ash-free fractions of fuel carbon and fuel water as explicit stoichiometric solutions, dependent at least in part on the reference fuel characteristics, the effluents O2, CO2, H2O and SO2, ambient concentration of O2, and air pre-heater leakage and dilution factors, finding the molar moisture-ash-free fractions of fuel nitrogen, oxygen, hydrogen, sulfur, sodium, potassium and chloride, converting the molar moisture-ash-free fuel concentrations to a molar dry base, then to a molar As-Fired wet base, and finally to As-Fired wet weight fractions, to obtain a complete and consistent computed As-Fired fuel chemistry, computing a heating value based on a moisture-ash-free weight base, then converted to a dry base, and then to a weight-based As-Fired heating value, and executing the mathematical model of the recovery boiler using the fuel information and the concentration of effluent O2 to produce consistent stoichiometric values of effluent CO2, SO2 and H2O values, the moles of fuel per basis moles of dry gaseous effluent, and at least the following self-consistent thermal performance parameters: As-Fired fuel flow, effluent flow, emission rates, boiler efficiency, and over-all system thermal efficiency.

51. The method of claim 50, including an additional step, after the step of executing, of performing analysis of instrumentation errors to obtain correction factors, and, if excessive, applying the correction factors to instrumentation signals such that subsequent on-line operation produces minimum errors in fuel chemistry and heating value determinations.

52. The method of claim 51, including an additional step, after the step of performing analysis of instrumentation errors, of adjusting operation of the system to improve its efficiency based upon the results.

53. A method for quantifying the operation of a recovery boiler burning black liquor fuel having a heat exchangers/combustion region producing combustion products, the method comprising the steps of: before on-line operation, the steps of obtaining a set of reference fuel characteristics, and developing explicit mathematical models of the combustion process involving at least stoichiometric balances; and thereafter operating on-line, the step of operating on-line including the steps of measuring a set of measurable operating parameters, including at least effluent concentrations of O2 and CO2, these measurements being made at a location downstream of the heat exchangers/combustion region of the recovery boiler, obtaining an effluent concentration of H2O if the set of reference fuel characteristics indicates that fuel water is not predictable, as an obtained effluent H2O, obtaining a concentration of O2 in the ambient air entering the recovery boiler, obtaining an air pre-heater leakage factor, calculating a set of fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts, as a function of the set of reference fuel characteristics, explicit mathematical models of the combustion process, the set of measurable operating parameters, the obtained effluent H2O, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

54. The method of claim 53, wherein the step of calculating the set of fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts, includes the step of calculating a set of moisture-ash-free fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts as a function of the set of reference fuel characteristics, explicit mathematical models of the combustion process, the set of measurable operating parameters, the obtained effluent H2O, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

55. The method of claim 53, wherein the step of calculating the set of fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts, includes the step of calculating a set of dry-based fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts as a function of the set of reference fuel characteristics, explicit mathematical models of the combustion process, the set of measurable operating parameters, the obtained effluent H2O, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

56. The method of claim 53, wherein the step of operating on-line includes the additional step after calculating the complete As-Fired fuel chemistry, of calculating an As-Fired fuel heating value as a function of the complete As-Fired fuel chemistry and the set of reference fuel characteristics.

57. The method of claim 56, including, after the step of calculating the As-Fired fuel heating value, the additional steps of obtaining a set of System Effect Parameters associated with the recovery boiler and its fuel, completing a multidimensional minimization analysis employing the set of System Effect Parameters to minimize the collective error associated with at least one of the measured effluent CO2, the obtained effluent H2O, the obtained fuel flow, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor, obtaining and applying for subsequent on-line analysis correction factors to the measured effluent CO2, the obtained effluent H2O, the obtained fuel flow, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

58. The method of claim 56, wherein the set of measurable operating parameters includes effluent temperature, and wherein the method includes an additional step, after the step of calculating the As-Fired fuel heating value, of obtaining a Firing Correction term, calculating a high accuracy boiler efficiency as a function of the complete As-Fired fuel chemistry, effluent temperature, the effluent concentrations, the As-Fired fuel heating value and the Firing Correction term.

59. The method of claim 53, wherein the step of operating on-line includes the additional step after calculating the complete As-Fired fuel chemistry, of calculating an As-Fired fuel heating value as a function of the complete As-Fired fuel chemistry and the set of reference fuel characteristics.

60. The method of claim 59, including, after the step of calculating the As-Fired fuel heating value, the additional steps of obtaining a set of System Effect Parameters associated with the recovery boiler and its fuel, completing a multidimensional minimization analysis employing the set of System Effect Parameters to minimize the collective error associated with at least one of the measured effluent CO2, the obtained effluent H2O, the obtained fuel flow, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor, obtaining and applying for subsequent on-line analysis correction factors to the measured effluent CO2, the obtained effluent H2O, the obtained fuel flow, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

61. The method of claim 47, wherein the step of calculating the complete As-Fired fuel chemistry includes the step of calculating explicitly a complete As-Fired fuel chemistry, including fuel water and fuel inerts, as a function of the reference fuel characteristics, explicit mathematical models of the combustion process, the set of measurable operating parameters, the obtained effluent H2O, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

62. The method of claim 47, wherein the step of obtaining the ambient concentration of O2 includes the step of using a value of 20.948 percent for the ambient concentration of O2.

63. The method of claim 47, wherein the step of obtaining the ambient concentration of O2 includes the step of using an average value at sea level determined by the National Aeronautics and Space Administration for the ambient concentration of O2.

64. The method of claim 47, wherein the step of obtaining the air pre-heater leakage factor includes the step of using a value of unity for the air pre-heater leakage factor.

65. The method of claim 53, wherein the step of calculating the set of fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts, includes the step of calculating explicitly a set of fuel chemistry concentrations including elemental fuel constituents, fuel water and fuel inerts as a function of the set of reference fuel characteristics, explicit mathematical models of the combustion process, the set of measurable operating parameters, the obtained effluent H2O, the concentration of O2 in the ambient air entering the recovery boiler, and the air pre-heater leakage factor.

66. The method of claim 53, wherein the step of obtaining the concentration of O2 in the ambient air entering the recovery boiler includes the step of using a value of 20.948 percent for the concentration of O2 in the ambient air entering the recovery boiler.

67. The method of claim 53, wherein the step of obtaining the concentration of O2 in the ambient air entering the recovery boiler includes the step of using an average value at sea level determined by the National Aeronautics and Space Administration for the concentration of O2 in the ambient air entering the recovery boiler.

68. The method of claim 53, wherein the step of obtaining the air pre-heater leakage factor includes the step of using a value of unity for the air pre-heater leakage factor.