CA2015090C - Method and apparatus for detecting flame - Google Patents

Abstract

METHOD AND APPARATUS FOR DETECTING FLAME

ABSTRACT OF THE DISCLOSURE

A method of detecting flame within a region where flame is expected. Radiation emissions from the region are measured within selected portions of the visible and infra-red frequency bands. Spectral characteristics of the two measurements, including their auto spectra, coherency and transfer function, are derived. The derived spectral characteristics are compared with prestored spectral signatures representative of the spectral characteristics of radiation emitted from the region within the selected portions of the visible and infra-red frequency bands while known flame conditions prevail within the region - thereby estimating the deviation of the derived spectral characteristics from the prestored spectral signatures. The deviations aforesaid are compared with predetermined threshold alarm values to assess the presence or absence of flame.

Description

- 20~`90 FIELD OF THE INVENTION

This application pertains to a method and apparatus for detecting flame and is particularly adapted to flame detec-tion in large boilers.

BACKGROUND OF THE INVENTION

Large boilers, for example, those used in conjunction with steam turbines for electric power generation, are fired by fuels such as coal, oil, gas or liquor. Supporting igniter burners are typically associated with each of the main burners.
Because the igniter burners are typically fired with relatively expensive fuels, they are operated only intermittently. More particularly, the igniter burners are preferably fired only upon initial start up of the boiler and thereafter they are only selectably fired for short intervals to light off or support flame at the particular main burner(s) associated with the igniter burner(s).

The prior art has evolved a variety of flame detec-tion techniques for monitoring boiler fires to detect the presence or absence of flame in the boiler regions supported by the various igniter burners. If flames are extinguished in a particular region of the boiler, then the "no flame" condition must be quickly identified or else the main burners continue to supply fuel which may potentially explode if it is not evenly and continuously ignited. Accordingly, highly reliable flame monitoring techniques are required for continuously detecting the presence of flame at regions within the boiler adjacent to each of the burners which fire the boiler.

~OlS090 The apparatus to be described in this application is suitable for use with two types of boiler/burner configurations;
namely, "wall" (or "opposed") fired hoilers, and "corner" (or "vortex") fired boilers. "Wall" or "opposed" fired boilers incorporate a series of burners mounted on two opposing walls of the four vertical walls of the boiler. Sighting tubes (pipes about 5 cm. in diameter) are positioned across the boiler walls (which are typically about 1.5 meters thick) beside and nearly parallel to each burner head. The sighting tubes are pointed approximately toward the expected location of burner flame.
Flame detection apparatus is positioned to "sight" through each tube into the boiler region in which flame is expected.

"Corner" or "vortex" fired boilers incorporate vertically separated stacks of burners in each of the four corners of the boiler. The flames produced by the burners merge in a central vortex within the boiler. The burners may be individually tilted in the vertical plane in order to bPtter control the combustion characteristics and location of the fireball within the boiler. Si~hting tubes for corner fired boilers must be flexible so that the flame detection apparatus can continuously track the flame as the burners tilt.

Several prior art flame detectors examine the light emitted by the flame and, from the time variation characteris-tics of these emissions, determine whether a flame is located near to the burner ("near flame"); or, a fireball is present in the background ("far flame"); or, there is no detectable flame.
By monitoring flame flicker (i.e. time variations in the light signal emitted in the freguency band(s) under consideration) such prior art detectors attempt to derive a binary signal rep-resentative of "flame" and "no flame" conditions. Pre-determined factors such as the geometry of the detector, the wavelength band it is capable of examining, and the frequency band being monitored affect the characteristics of flame flicker and correspondingly determine the ability of such detectors to 20~09~) accurately detect the presence or absence of ~lame under varying conditions.

The best prior art flame detectors for use on opposed fired boilers appear to be those which utilize two separate linear arrays of detectors aligned horizontally and vertically to facilitate "X-Y" scanning of selected sub-regions within a region where flame is expected, through electronic selection of an appropriate detector pair. Typically, a zero-crossing waveform shaping analysis is performed on the electronic signals produced by each of the two selected detectors, to generate two bi-level output signals. The output signals are then correlated with one another ~prior art detectors of this sGrt do not however perform true signal correlation, because they work only with binary ~i.e. two level) approximations of the detector output signals, rather than with the direct analog outputs of the detec-tors). If the two signals are highly similar to one another then the correlation result approaches unity. Normally, a result which exceeds some predetermined threshold is accepted as indicating the presence of flame. If the two signals are highly dis-similar to one another then the correlation result approaches zero. A result which does not exceed the aforementioned threshold is normally taken to indicate a "no flame" condition.
In some cases, automatic tracking techni~ues are employed to locate points of maximum correlation in an effort to minimize generation of false "no flame" alarms. It will thus be under-stood that the prior art is susceptible to error, in that the cumulative approximations inherent in the operation of prior art detectors may result in a "no flame" alarm when flame is in fact present; or, may indicate that flame is present when no flame is in fact present. The prior art tends to overcompensate for these contingencies by allocating flame determination thresholds which minimize generation of false "flame present" signals. How-ever, this necessarily decreases the ability of such prior art devices to respond to flame conditions having light emission characteristics which encompass a large dynamic range.

2(~1~0~

The inventors believe that superior results may be obtained by concentrating on factors other than flame flicker.
More particularly, the inventors believe that superior results may be obtained by analyzing the time ~ frequency spectral characteristics of light emitted in different visual and infra-red wavebands from the region in which flame is expected, and comparing those characteristics with prestored spectral signa-tures representative of flame. ~he present invention accordingly compares short term estimates of the visible and infra-red au~o-spectra, the infra-red to visible transfer function, and the in-fra-red to visible coherence (all of which are hereinafter defined and explained in greater detail), with prestored signatures characteristic of "flame" and "no flame" conditions.
The auto-spectra, transfer function and coherence function are used to characterize the relationship between two signals in selected frequency bands. It is this relationship or pattern which is used to identify the flame.

SUMMARY OF THE INVENTION
In accordance with the preferred embodiment, the invention provides a method of detecting flame within a region where flame is expected. The method comprises the steps of measuring radiation emitted from the region within a selected portion of a visible frequency band, concurrently measuring radiation emitted from the region within a selected portion of an infra-red frequency band, deriving the coherency between the two measurements, comparing the derived coherency with a prestored coherency signature representative of the coherency between measurements of radiation emitted from the region within the selected portions of the visible and infra-red frequency bands while known flame conditions prevail within the region -thereby estimating the deviation of the derived coherency fromthe prestored coherency signature, and comparing the deviation with a first predetermined threshold alarm value.

201~90 The auto spectrum of the visible frequency band measurements is also derived. The visible auto spectrum measurement is then compared with prestored auto spectrum signatures representative of the auto spectrum between measure-ments of radiation emitted from the region within the selectedportion of the visible frequency band while known flame condi-tions prevail within the region - thereby estimating the deviation of the derived visible measurement auto spectrum from prestored visible auto spectrum signatures. The deviation of the derived visible measurement auto spectrum from prestored visible auto spectrum signatures is then compared with a second predeter-mined threshold alarm value.

The auto spectrum of the infra-red freguency band measurements is similarly derived. The infra-red auto spectrum measurement is then compared with prestored auto spectrum signatures representative of the auto spectrum between measure-ments of radiation emitted from the region within the selected portion of the infra-red frequency band while known flame condi-tions prevail within the region - thereby estimating the deviation of the derived infra-red measurement auto spectrum from prestored infra-red auto spectrum signatures. The deviation of the derived infra-red measurement auto spectrum from prestored infra-red auto spectrum signatures is then compared with a third predetermined threshold alarm value.

The transfer function between the visible and infra-red frequency band measurements is also derived. The derived transfer function is compared with prestored transfer function -~-signatures representative of the transfer function between measurements of radiation emitted from the region within the selected portions of the visible and infra-red frequency bands while known flame conditions prevail within the region - thereby estimating the deviation of the derived transfer function from the prestored transfer function signatures. The transfer function deviation is then compared with a fourth predetermined threshold alarm value.
' "' ~b ~ , J ~ , y ~, ~ x t i* i ~. ~X,: ~

20150~

The measurements are repeated for other separateselected portions of said visi~le and infra-red frequency bands and the various spectral signature deviations aforesaid deter-mined for each frequency band portion. A weighted least squaresfit; or, a stochastic fit; or, a bounded limits fit; or, a Gaussian fit is applied to the derived and prestored spectral signatures. The weighted spectral signatures derived from separate frequency bands are normalized, averaged and summed, then compared with a plurality of prestored corresponding spectral signatures, the prestored signatures being representa-tive of a selected flam~ conditions.

BRIEF DESCRIPTION OF THE DRAWIN~S

Figure 1 is a block diagram which illustrates the basic components of a flame detection system constructed in accord-ance with the preferred embodiment of the invention.

Figure 2 is a partially fragmented longitudinal cross-sectional illustration of a direct sighting scanner head assembly constructed in accordance with the preferred embodiment.
Figure 3 is a longitudinal cross-sectional illustra-tion of an extended direct sighting scanner head assembly constructed in accordance with the preferred embodiment.

Figure 4 is a partially fragmented longitudinal cross-sectional illustration of a fiber optic flexible scanner head assembly constructed in accordance with the preferred embodiment.

Figure 5 illustrates diagrammatically how discrete viewing windows are established by the preferred embodiment of the invention.

201~0~
Figure 6 is a cross-sectional illustration depicting the placement of an extended direct sighting scanner head assembly within a boiler wall and the range of viewing windows thereby obtained within a region of expected flame.

Figure 7 is a schematic illustration depicting the viewing window trigonometry applicable to the case in which the photocell or fiber optic termination point "P" lies on the focal plane.
Figure 8 is a schematic illustration depicting the trigonometry applicable to the situation in which the point "P"
lies in front of the focal plane.

Figure 9 is a schematic illustration depicting the trigonometry applicable to the situation in which the point "P"
lies behind the focal plane.

Figure 10 is a schematic illustration depicting the determination of windows for non-point source sensors; Figure lO(a) depicting the situation in which the sensor lies on the focal plane; and, Figure lO(b) depicting the situation in which the sensor lies behind the focal plane.

Figure 11 is a block diagram of the construction of the flame scanner head electronics of the preferred embodiment.

Figure 12 is an electronic circuit schematic diagram of the flame scanner head electronics of the preferred embodi-ment.

Figure 13 is a flowchart of the flame detection algorithm which controls the operation of the preferred em-bodiment of the invention.

:

..-...

2 a ~ o DETAILED ~ESCRIPTION OF THE PREFERRED EMBODIMENT

Principle of OPieration The primary combustion zone of a boiler flame can reach temperatures of 1800K. At this temperature the blackbody or greybody radiation emitted by the flames peaks in the near infra-red range of the spectrum. As the temperature increases, the peak energy wavelength shifts towards the visible or shorter wavelength region of the spectrum. Similarly, as the temp-erature decreases, the peak energy shifts towards the infra-red portion or longer wavelength region of the spectrum.
"Dual-colour" sensors of the type marketed by Hamamatsu Photonics k.k. of 1126 Ichino-Cho, Hamamatsu City 435, Japan under the part numbers K1713-01 (U.V. enhanced Si/PbS), K1713-02 (ULV~ enhanced Si/PbSe) and 1713-03 (U.V. enhanced Si/Ge) can simultaneously monitor both the visible and infra-red spectra emitted by individual burner flames. Suitable dual colour sensors may also be obtained from Infrared Industries Inc., of Orlando, Florida.
Although a dual-colour sensor is employed in the preferred embodiment, the invention is not limited to two colour detec-tion (i.e. sensors capable of sensing radiation in a multi-plicity of wavebands may be employed).

Combustion is a non-stationary process which can be characterized by the flicker or A.C. content observed in the infra-red and visible emissions of the primary flamefront. In the preferred embodiment, this A.C. flicker content is separ-ately monitored by the visible and infra-red sensors of a dual colour sensor over a frequency range of about 5 Hz to about 500 Hz. The resultant time dependant output signals tend to be correlated with each other. It has been found that there is a high coherency between the visible and infra-red sensor outputs in selected frequency bands when flame is present at a burner, but that the coherency is reduced when flame is not present.
This has been found to be true even in the presence of background ~ignals from other burners; and, to a lesser extent, in the presence of coal shrouding, which tends to pass the infra-red but ~01~09~

not the visible spectra. This variation in coherency can be partially explained by the fact that the infra-red and visible elements of a dual colour sensor each perceive slightly different angular windows. The divergence (i.e. difference in cross-sectional area of each window) between the infra-red sensor viewing window and the visible sensor viewing window increases with increasing distance from the scanner. This results in lower coherency between the two sensor output signals when the flame location is far from the burner being monitored (i.e. background flame or fireball). When the flame is located directly in front of the burner (i.e. near the sensor) the windows are nearly coincident and the emission spectra, as seen by the dual-cblour sensor, tends to be highly coherent (i.e. correlated). In addition, far flames have lower frequency characteristics than near flames, due to the integration effect over a larger cross-sectional window. Thus, the coherency also varies differently in different frequency bands.

Coherency between two time varying signals X(t) and Y(t) is defined as:
I ~xy( Cxy (w) ~xx ( j ~) ~YY ( i ~) where:
Cxy = squared coherency ~xy = cross-spectrum between X(t) and Y(t) ~xx = auto-spectrum of X(t) ~yy = auto-spectrum of Y(t) ~ = frequency (radians per second) j = complex root of (-l).

The coherency function varies with frequency and is limited by: -0 < Cxy(w) < 1.0 When the two signals X(t), Y(t) are linearly related the coherency tends toward unity, otherwise the coherency tends to zero.

_ g - ' 20~5090 In accordance with the preferred embodiment of the invention, short term estimates of the coherency between the visible and infra-red emissions from the flame, as detected by the dual colour sensor, are compared with prestored characteris-tic coherency signatures for the particular burner over a time domain frequency range of 5-500 Hz. The deviation of the short term coherency estimate from the prestored "ideal" signature value is integrated over the frequency range of interest using a weighted difference cost function. This integrated "cost"
estimate is then compared with a threshold alarm value, to determine the presence or absence of flame. The signature comparison approach, described above for the coherency function, is also used to compare the difference in short term estimates of the visible and infra-red auto-spectra and infra-red ~ visible transfer function gain with corresponding "ideal" prestored "flame" and "no flame" signature spectra. These short term spectral estimates may be compared with several prestored characteristic signatures to determine the most likely flame condition. The results of the comparison tests on coherency, visible auto-spectrum, infra-red auto-spectrum, and infra-red -visible transfer function may be individually weighted, by frequency and by function, and summed to form an overall measure of flame condition.
Mechanical and Optical Desiqn Criteria Flame detectors constructed in accordance with the invention preferably satisfy the following design criteria:
(1) The flame detector is compact, rugged and easily retro-fitted to existing boiler sighting tubes. The maximum front lens diameter (typically < 50 mm) is limited by the size of sensor head that can be installed in the boiler sighting tube. Practical constraints of cost and standard manufacturing sizes limit the front lens diameter to < 25 mm in most cases.

2~15090 (2) The flame detector is able to withstand moderately high temperatures (< 300C).

(3) The flame detector is able to operate in an abrasive and dirty environment without scouring or slagging of the lens assembly occurring. This is achieved by using an air supply to both cool and clean the optical components. If this approach is taken, then provision must be made to supply air to cool the apparatus and to purge and clean the optics.

(4) Lenses are easily replaceable in order to best match the optics to a specific burner design.

(5) The optics should ideally pass wavelengths in the range of 0.2 ~m < ~ < 5.0 ~m using zirconium fluoride fiber optics, although alternative embodiments of the invention may use quartz optics (which limit the upper passband to ~ 2.5 ~m).

(6) The optics permit monitoring of adjustable selected viewing windows in front of the burner. These windows are adjustable in both the longitudinal and lateral directions.
- '.

(7) The flame detector may be operated with a variety of different sensors.

(8) The signal conditioning electronics in the sensor head maximizes the signal to noise ratio from the sensor in the 5 Hz. to 500 Hz. frequency band and includes high -frequency roll-off filters to eliminate signal aliasing.
As illustrated in Figure 1, the preferred embodiment provides for one or more "scanner heads" consisting of a sighting -- 11 -- , ...... ... ... , .. . . . . ..... , , .-.. . . . .. - . . - - . ... . .. . ~ ..... . . . ... . . - . ..

tube which may be positioned with:in one of the burner viewing ports located across the boiler wall. The tube contains the viewing optics, dual colour sensor(s) and supporting electronics (each hereinafter described in greater detail). A communications S link couples the scanner head electronics to a computer. In the preferred embodiment, the computer is an IBM~ personal computer with a co-processor board adapted to monitor the flame signals and independently capable of detecting and signalling flame condition. Optionally, output signals may be provided to support the operation of a separate burner management system using the relay contact outputs provided by the co-processor board to control fuel and air flow to the burners.

Confiquration O~tions 1~ The preferred embodiment provides three different options for configuring the scanner head. These are:

(1) Direct Siqhting Head. Figure 2 shows the basic elements of a direct sighting flame scanner head 10, in which the lateral and longitudinal displacement between an array 12 of dual colour sensors 12a through 12e and the lens 14 can be varied to select the viewing window, as shown by the arrows "X" and "Y". The direct sighting head is used where sighting through a simple viewing port 16 is possible. When used with a sighting tube (not shown) the effective viewing angle (window) may be limited by the sighting tube. A camlock mechanism 18 is provided to lockably engage notches 20 on scanner head 10, to hold the head in position relative to mounting plate 22. The scanner head electronics are diagrammatically represented at 24.
Coupler 26 is provided for receiving a cable for conveying electrical signals to and from electronics 2~.
(2) Extended Siqhtinq Head. As depicted in Figure 3, the extended direct sighting head is similar to the direct 201sO~û
sighting head of Figure 2, the basic difference being the provision of fiber optic cable 30 between dual colour sensor 12' and lens 14'. Flame position may fluctuate and move out of range of the sighting angles S as limited by the sighting tube. To remove this restriction, the extended direct sighting head of Figure 3 collects light over wider angles at the front of the sighting tube. The device is air cooled by passing cooling air through port 31.

(3) Flexible Fiber optic Scanner Head (Figure 4). The flame scanner must be able to track flame in corner fired boilers at all burner tilt angles. Due to the wide range of possible flame locations, a flexible fiber optic head assembly is required to track the flame. Both the outer guide tube 32 and the inner -scanner head 34 are constructed so that they are able to flex. In all other respects the flexible fiber optic scanner head is identical to the extended 2~ sighting head.
'`,: -The three scanner head designs vary significantly in the way in which the flame emissions (visible and infra-red) are directed to the sensors. This is hereinafter explained in greater detail.
--Wideband Optics Sapphire lenses and windows are preferably used throughout. However, alternative materials, such as silicon quartz, may be used with some degradation in performance. Thus, although reference is made to the properties of sapphire lenses and windows, similar properties exist for a range of optical glasses. Similarly, zirconium fluoride fiber optics are preferred, although these too can be replaced by quartz glass equivalents with some degradation in performance.

A sapphire window in front of the sensors protects the sensor material, while passing all wavelengths of interest. The 2al~0so use of a sapphire lens ensures yood transmittance characteristics over the full optical range. The advantages of sapphire are: it is chemically inert and therefore not easily corroded; it is very hard and not marred by most abrasive materials; it is very strong, allowing the use of thin lenses; it withstands high temperatures; and, it has a high thermal conductivity, which aids artificial cooling.

Although sapphire does exhibit a birefringence due to its crystalline properties, this has negligible impact on the optical performance of the scanner.

Adiustable O~tical Path (Viewing Windows) The optical path from the flame to the sensor is adjustable. Five basic adjustments are possible. These are:

(1) The choice of lens focal length. The scanner head barrel length dictates that the lens focal length should be significantly less than the maximum distance that the sensor can be positioned hehind the lens.
Conventionally available plano-convex sapphire lenses have design focal lengths of 100 mm, 50 mm or 25 mm.
Other custom design focal lengths are available.

(2) The use of an aperture plate which limits the apparent lens diameter. This feature is not ordinarily employed but can effectively determine the viewing window in conjunction with item (3) below.

(3) The distance from the sensor to the lens along the viewing axis determines the viewing window in con-junction with item (2) above.

(4) The lateral position of the sensor, off the principal axis, determines the viewing window offset angle.

201~09() (5) Sensor dimensions (i.e. cross-sectional area and shape).

In each of the three scanner head design configuration options hereinbefore mentioned, the first four parameters are independently adjustable to meet particular viewing window requirements. The fifth parameter, namely the relative dimen-sions of the preferred silicon and lead selenide/sulphide dual colour sensor, also determines the size of the visible and infra-red viewing windows, but is a parameter which can only becontrolled at the time of ordering the sensor from the manufac-turer.

Multiple Sensor / Array Scannina Capabilit In addition to sensor 12a (Figure 2~ which lies on the principal longitudinal axis of the scanner head, up to four more sensors 12h, 12c, 12d and 12e can be placed at progressively greater lateral distances off the principal axis to provide a linear optical array which can be selectively scanned. Sensor array 12 is able to discriminate and dynamically track the movement of the burner flame over a wider viewing angle than would ~e possible with a single sensor, while maintaining a narrow viewing acceptance angle for individual sensors (and hence retaining good A.C. flicker signal characteristics).
MultiPl~ Sensor TvPes The flame detection apparatus can be configured with three types (Si/PbS, Si/PbSe or Si/Ge) of dual-colour sensors which use four basic sensor materials. These are:
(1) Silicon (Si) (photovoltaic) sensor operating in the visible wavelength range from 0.2 ~m to 1.15 ~m; cell size ~ 2.5 mm x 2.4 mm (custom dimensions are available from the sensor manufacturer for selecting particular viewing window characteristics).

201~0~0 (2) Lead Sulphide (PbS) (photoresistive) sensor operating in the infra-red wavelength range from l.1 ~m to 2.5 ~m; sensor size ~ 2.0 mm x 2.0 mm (again, custom sensor sizes are available from the sensor manufacturer).

(3) Lead Selenide (PbSe) (photoresistive) sensor operating in the infra-red wavelength range from l.1 ~m to 4.85 ~im; sensor size ~ 2.0 mm x 2.0 mm.

(4) Germanium (Ge) (photovoltaic) sensor operating in the infra-red wavelength range from l.1 ~im to 1.9 ~m;
sensor size ~ 2.0 mm diameter.

These sensors are housed in an industry standard T05 package and are available from several sensor manufactures, including the two previously mentioned. Custom sized sensors are also available.

Dual Colour Detector The sensors are constructed as two-colour detectors.
A silicon (Si) photovoltaic sensor detects incident radiation in the visible range. This is superimposed in front of the appropriate infra-red sensor substrate. Since these sensors are thin films, they are effectively coplanar. The sensor elements are constructed to be symmetrical about a central axis, but are of different dimensions. The active sensor area of each material can be varied to achieve the desired viewing window characteris-tics. This, however, is a one time choice, made at the time the sensor is ordered from the manufacturer.

Multi Colour Detector Although the preferred embodiment herein described em ploys dual-colour (i.e. visible and infra-red) sensors as described above, three colour sensors having silicon (Si), ger-manium (Ge), and one of lead sulphide (PbS) or lead selenide (PbSe) detectors are available. The principle of detection remains the same, except that the auto-spectra, coherency and transfer function can now be estimated for three pairs of -- 2~1~0~
signals, as given by: Si - Ge; iSi -~ PbSe; and, Ge ~ PbSe. The principle of flame detection is unaltered, but the variation and sensitivity to small changes in flame state are enhanced.

Dlrect Sighting Optics Optical Layout Figure 5 shows the array scanning concept, whereby a dual colour sensor array 12 comprised of five dual colour sensors numbered 1 through 5 in Figure 5 (one of which, namely sensor 3, lies on the principal axis and the others are vertically displaced above and below the principal axis, as shown) may be electronically scanned to select one of the five sensors which "sees" through lens 14 into a particular viewing window within the boiler. For example, the dashed lines in Figure 5 illustrate the viewing window of the lowermost sensor 5, as determined by the height of the sensor, its vertical displacement off the principal axis, the distance "X" from sensor 12 to lens 14, and the lens focal length. The viewing window of sensor 5 has a mean viewing angle ~5 given by tan1(Y/X), where "Y" is the vertical displacement of the sensor relative to the principal axis. The mean viewing angles of the windows "seen" by the other four sensors are indicated in Figure 5 as ~ 2~ d3, and 94 respect-ively. -As previously explained, each dual colour sensor incor-porates separate visible and infra-red sensors. These each "see"
slightly different windows within the boiler, as illustrated in Figure 6. Fuel 50 fed through burner 52 ignites to produce flame 54. Direct sighting scanner head 56 is mounted in boiler wall 58 at an angle relative to burner 52, so that the sensors within scanner 56 can "see" the region in front of burner 52 in which flame 54 is expected. The visible sensor component of the dual colour sensor within scanner 56 "sees" a "visible window" having top and bottom visibility limits as indicated in Figure 6. The infra-red sensor component "sees" a somewhat narrower "infra-red window" having top and bottom visibility limits which are also indicated in Figure 6.

201509~

The optics will now be disc~ssed in greater detail with reference to Figures 5 through 10. Xey symbols are labelled on the drawings and are defined in the list of symbols hereinafter provided.

In the direct sighting head design (Figure 10) the dual colour sensor 12 is placed perpendicular to the principal axis IlPA'' and located a distance "X" behind the secondary principal point of lens l4. The midpoint of the sensor may also be offset a perpendicular distance "Ym" from the principal axis. Both the visible and infra-red sensors have finite dimensions + YV;S and + YIR respectively as measured from the midpoint of each sensor.

The offset 'IYm'' determines the sensor midpoint viewing angle ll~mll~ As shown in Figure 10, the surface area of each sensor absorbs incident energy that has been diffracted by lens 14. Since the dimensions of both sensor 12 and lens 14 are finite, energy sources located in front of lens 14 can be observed by sensor 12 over a range of angles. These angles are determined by the location of sensor 12 relative to lens 14 and by the lens and sensor dimensions. Projecting light rays forward from sensor 12 defines the dimensions of a viewed window at any given distance "L" in front of lens 14. Referring to Figures 7, 8 and 9, the following three window configurations are possible:

(l) The sensor is located on the focal plane (i.e. at point "P" shown in Figure 7). The viewed window diverges with respect to lens 14 due to the finite sensor dimen-sions. The sensor will only be on the foca~ plane for a particular wavelenyth, ~0, of incident light. The lens focal length decreases for shorter wavelengths (< A~) and increases for longer wavelengths. As the sensor responds over a band of wavelengths the window angle is implicitly also a function of wavelength (see "Window Design" below).

201~90 (2) If the sensor is located in front of the focal plane (i.e. at point "P" shown in Figure 8) it observes all emitted energy between the widely diverging angles ~tmin and ~h~x (see Flgure 10).

(3) If the sensor is located behind the focal plane (i.e.
at point "P" shown in Figure 9) two possible windows exist:
(i) Near the lens the window defined by et and ~b converges.
(ii) At the point beyond where the ray subtended by angle 9b crosses that subtended by ~t the window diverges. In this case, the finite sensor dimensions depicted in Figure 10 result in two (near and far) convergent points and hence define a focal range.

Siqhtinq Options for the Sensor The scanner head design allows a number of parameters to be easily changed. These design options are:

1. The material that the lens is made of. This determines the maximum optical bandwidth that can be detected.
The resulting variation in the index of refraction with wavelength affects the viewing window size, as the lens focal length is a function of wavelength.
.: .
2. The type of visible or infra-red sensor used to detect -the radiant energy. This also determines the optical bandwidth that is detected.

3. The linear dimensions of each sensor. This determines the shape of the observed window. Larger sensor dimen-sions provide a larger viewing window.
4. The ratio of the linear dimensions and areas of the visible and infra-red coplanar sensors. At any given 201~09~
wavelength the visible and infra-red window sizes are proportional to this ratio.

5. The lens diameter or intermediate aperture plate dia-meter. This determines the total energy striking the sensor and also affects the dimensions of the viewing window. A larger aperture allows more energy to strike the sensor, resulting in greater sensitivity at low energy thresholds and in a larger viewing window.

6. The lens focal length. The viewing window dimensions are inversely proportional to the focal length. A
longer focal length provides a narrower viewing angle.

7. The sensor offset 'IYm" location. This parameter determines the angle of the optical axis relative to the principal axis. This allows offset viewing angles relative to the principle mounting axis of th~ scanner head.
8. The sensor "X" location. The relationship between "X"
and the dimensions of the viewing window is nonlinear and depends on the sensor's location relative to the focal plane. (This is discussed in greater detail below under the heading "Window Design").

Viewinq Window Desiqn And Selection Criteria - -Selection criteria As indicated by Figure 6, the flame scanner head 56 is typically located in a burner viewing port tube located near the burner 52 being monitored. The tube is canted slightly towards the burner so that the axis of the tube will intersect the burner flame axis at a location near where flame 54 is expected.
Assuming the tube dimensions do not limit the viewing window, the optics can be optimized to observe a specific window area located a distance "L" in front of the sensor head for any particular wavelength. The variation in window area, for the visible and 201~0gO

infra-red, should be minimal across the desired optical bandwidth at the design distance "L". This can be approximately attained by careful design and selection of the sighting options listed above. This is an iterative procedure which may be aided by the use of a computer program to calculate the viewing window as a function of all of the relevant optical parameters. The theoretical basis for the required program is developed below under the heading "Window ~esign'l.

The signals from both the visible and infra-red sensors are sent to a remote processor. It has been determined that the A.C. amplitude signals from the visible and infra-red sensors measured over a 5-500 Hz. bandwidth contain the most useful information. The auto-spectra, transfer function and coherency of these two signals are estimated over short time intervals to determine the flame condition. ~he relative dimensions of the visible and infra-red windows may have to be adjusted in order to extract the maximum useful information from the observed flame.
Window Desiqn The optical theory underlying the invention will now be developed for a sensor assumed to be a point source or sink.
This derivation will then be extended to cover the two dimen-sional case where the sensor is assumed to be of a finite length. Finally, a three dimensional derivation, assuming a sensor having finite length and width, is presented.

Throuqhout these derivations, ray tracing techniques are used to determine the imaging characteristics of a lens. The rays possess the following properties:

(1) Rays are diffracted, or bent, only by the lens and continue unimpeded in straight lines on either side of the lens.
(2) All rays passing through the principal point (Pp) exhibit no diffraction and therefore continue with no .' ' ' . . ' . ' ' ' ' ', ,'' , , " ' , ," .' . ' ' . ', ' ~ ' . ' ' , I ' ' ' ' ''; ., ' ' ' .. . . . . ' ' , ~ 3 - ; ? j i ~ . ., :'': : :, :' . '.: : i , . . . , , , , ., " ,, ,, , , , , ,, , ,, ,. ,, , ~, ,: .:

201~09~

change in direction on ~oth sides of the lens. Each individual ray intersects the principal point (Pp) at an angle (~) relative to the principal axis.
(3) Ray paths are completely reversible, yielding the same results whether the rays are traced from in front of the lens to behind the lens or in the reverse direc-tion.
t4) All rays emanating from an arbitrary point on the focal plane and passing through the lens will be diffracted so that they continue in a parallel line in front of the lens.

The design wavelength (~0) of sapphire, at which the manufacturer specifies the optical properties of lenses, is 0.5461 ~m. At wavelengths (~) other than the design wavelength the refractive index (n) of sapphire and hence the focal length of the lens can be calculated from an empirical equation provided by the manufacturer. This equation is:

A1 ~2 A2 ~2 A3 ~2 n2 ~ 2 (1) ~ 2 ~ ~3 rearranging equation (1) gives:

n = A~ ~Z A2 ~2 A3 ~2 (2) where:
Al = 1.023798 ~1 = 0.00377588 A2 = 1.058264 ~2 = 0.0122544 A3 = 5.280792 ~3 = 321.3616 20~50~

for thin lenses:

= (n-l) _ (3) f r1 r2 for a plano-convex lens r2 = ~- Hence equation (3) reduces to:

_ = (n-l) (4) f r1 Therefore, substituting the nominal design focal length "fO" into equation (2) gives the design index of refraction "nO". Rearranging equation (4) and substituting the nominal design focal length "rO" which is:

rO = fO (nO~l) (5) The focal length "f" at any arbitrary wavelength ~ can now be calculated from:
rO
f . = fO[ (nO - lL ] (6) (n-1) ( n - 1) where the index of refraction "n" is calculated from equation (2).

Based on ray tracing techniques and using the approp-riate symbols and definitions, Figure 7 schematically illustrates the paths of the rays passing through the top "rt" middle "rm"
and bottom "rb" of the lens and converging to an arbitrary point "P" on the focal plane. The following definitions should be noted:
1. The principal axis 'IPA'I iS defined to be centred on, and perpendicular to, the surface of lens 14.
2. The principal surface is an imaginary surface where all rays parallel to the principal axis in front of 201~0~

the lens are singly re~racted to come to a focus at the rear focal point ~IPf~l.
3. The principal point"Pp" is located at the intersection of the principal surface and the principal axis IlPA''.
4. All dimensions along the principal axis IlPA'' are measured from the principal point "Pp".
5. The focal length "f" is the distance from the principal point to the rear focal point IIPf~.
6. The lens has a finite centre thickness "tc" and edge lo thickness "te".
7. The lens has a finite aperture diameter "~".
8. The lens has a design radius of curvature "rO".

In order to simplify the derivation of the optical equations the following assumptions have been made:

(1) The actual plano-convex lenses being used are quite thin; consequently, it has been assumed that the lens thickness, both "tc" and "te", has been reduced to zero.

(2) This results in the secondary principal surface being a plane centred on and perpendicular to the principal axis 'IPA" and having a diameter equal to the lens diameter "~". The lens is reduced to a single diffracting plane.

(3) Since the individual sensors can have large lateral offsets, IlYm", it has also been assumed that the focal plane is in fact a hemisphere centred at the principal point "Pp" with a spherical radius equal to the focal length "f".

The flame in front of the lens is not necessarily focused as an image behind the lens. It is only necessary to calculate the angular limits of the viewing window in front of the lens to determine which radiation sources will be viewed by 201~090 the sensor. Each sensor is activated by the total optical energy incident on its surface in the sensor bandwidth, irrespective of the source of that energy.

It will now be shown that for a point "P" arbitrarily located behind the lens, the window angles measured at the top and bottom of the lens can:

(1) result in a parallel viewing window in front of the lens if point "P" is located on the focal plane (Figure 7); or, (2) result in a diverging viewing window in front of the lens if point "P" is located in front of the focal plane (Figure 8); or, (3) result in a converging then diverging viewing window in front of the lens if point "P" is located behind the focal plane (Figure 9).
The radius "R" from the principal point to an arbitrary point "P" located ~ehind the lens is given by:
.
R = ~ x2 + y2 The corresponding angle "~" subtended by the principal point to . -the point "P" relative to the principal axis is given by ~ =
arctan (Y/X).

As defined in Figures 7, 8, 9 and 10, the following conventions hold:
:
~1) Relative to the principal axis, IIPA", positive angles are measured upwards behind the lens and downward in front of the lens.

201~u (2) With respect to the principal axis, IlPA'' positive, "Y"
dimensions are upward behind the lens and downward in front of the lens.

(3) Positive "X" dimensions are measured from the principal point "Pp" along the principal axis IlPA" behind the lens.

(4) Positive "L" dimensions arP measured from the principal point "Pp" along the principal axis IlPA" in front of the lens.

(5) The angular windows in front of the lens are measured from the top and bottom edges of the lens, parallel to the principal axis.

Figure 7 illustrates the case in which point "P" is arbitrarily located on the focal plane. Based on the principles of ray tracing, the middle ray "rm" traverses both point "P" and the principal point "Pp" with no change in direction. This determines the angle "~" both in front of and behind lens 14.
Both the top ray "rt" and the bottom ray "rb" converge at point "P" then continue to diverge behind point "P". In front of lens 14, all rays are parallel to the middle ray and subtend an angle "~" to the principal axis.

Figure 8 illustrates the case in which point "P" is located in front of the focal plane. As in Figure 7, the middle ray "rr" traverses both point "P" and the principal point "Pp"
with no change in direction. This determines the middle ray viewing angle "0". The top ray "rt" however intersects the focal plane at point "Pt" and the bottom ray "rb" intersects the focal plane at point IlPb~'. The angle ll~tll at which the top ray "rt"
enters the top of the lens is determined by the angle of the ray intersecting both point "Pt" and the principal point "Pp".
Similarly, the angle ll~bl' at which the bottom ray "rb" enters the bottom of the lens is determined by the angle of the ray ;' ~'.

201~090 intersecting both point IIPb" and the iecondary principal point " Pp" .

If point "P" is on the focal plane as shown in Figure 7, then:
Gt = b and the viewing windows in front of the lens are parallel to one another and therefore constant at all locations.

If point "P" is in front of the focal plane as shown in Figure 8, then:
~ t < ~ < ~b and the viewing window in front of the lens diverges.

If point "P" is behind the focal plane, as shown in Figure 9, then:
~t > > ~b and the viewing window in front of the lens converges to a focal point, then diverges.
To extend this theory to determine the viewing window for a sensor with finite "Y" dimensions (Figure 10) the above calculations are repeated for the following three "Y" locations:
:
(l) The top of the sensor at IIYm + Yvjs or Ym + YIR
(2) The mid point of the sensor at IlYm''.
(3) The bottom of the sensor at IIYm - Yvjs or Ym ~ YIR -For each of these three sensor locations the top, middle and bottom viewing window angles are calculated: t, and ~b respectively. From these calculations the following angular limits are determined:

(1) The maximum top window angle l'tm3xl' (2) The minimum top window angle ll~tminl'.
(3) The maximum bottom window angle''~bm3x''.

20~0~() (4) The minimum bottom window anyle l'O~jn".
(5) The sensor mid point viewing angle "~m" calculated at Ym If t~x < ~in~ then the viewing window angle diverges continuous-ly from the bottom of the lens at angle g~x Similarly, if ~tmin< ~lin~ then the viewing window angle diverges continuously from the top of the lens at angle tmin. Alternatively, if ~t~x > ~x~
then the viewing window angle from the bottom of the lens is determined by ~x~ until ~t~x intersects 6~x, then the window angle is determined by ~t~x- Similarly, if ~tmin > ~in~ then the viewing window angle from the top of the lens is determined by ~tmin~ until ~in intersects Gtmjn, then the viewing window angle is determined by ~in-The linear dimensions of a window (Figure 10) located on the optical axis "O" a specific distance "D" in front of the lens at a mid point angle ~mll can also be calculated. The top and bottom dimensions measured from the optical axis dt and dB, must be calculated separately then added together.
In order to extend these calculations to a three dimen-sional configuration, the two dimensional derivation is repeated for the width of a specific sensor. The resulting angular and linear lengths and widths are then multiplied together to obtain the actual observed solid angle and cross-sectional window areas.

A computer program which implements the foregoing calculations facilitates selection of the best combination of lens, sensor and sensor position for any given application.
Since any multiple lens system can be combined to yield an equivalent single lens system, this same technique is readily expandable from the lens direct sighting case to scanners having more sophisticated optics. Computer simulations have shown that the scanner head variables are interdependent. As an example, the viewing window angles vary with the wavelength of the observed radiation. This means that for a given set of input .

20~09~

variables the resulting apparent window can vary significantly over the full range of wavelengths being observed. This property is used to select different window properties for the visible and infra-red sensor elements. The windows are chosen so that they approximately coincide at the expected flame location, but diverge at other locations. Thus the sensor outputs tend to be highly coherent when flame is present, but less so otherwise.

Extension of Window Theory for Extended Direct Siahting and Fiber opticS Scanner Heads The viewing window theory developed for the direct sighting head is applicable to the extended and flexible fiber optic scanner head designs. In these cases the incoming flame radiation is focused onto a fiber bundle termination plate. The fiber optic bundle dimensions are substituted for the sensor dimensions in Figure 10 and the theory of operation is replicated exactly as long as the following conditions hold:

(l) The angle subtended by the incident radiation to the principal axis of the fiber bundle is less than the acceptance angle of the bundle (typically < + 25).

(2) All the energy transmitted by the bundle is focused on the active sensor area at the remote end of the fiber.

Given these two constraints, which are easily met in practice, the fiber optic viewing window is identical to the direct sighting window. Positioning the fiber optic termination point with respect to the plano-convex lens facilitates adjust-ment of the viewing offset angle and window.

Scanner Head Electronics overview The flame scanner head electronics (Figure 11) provide signal conditioning and channel selection for up to four dual colour sensors located in the scanner head. The printed circuit board on which the electronic components are mounted in turn mounts in the scanner head barrel, and is shielded using a mu-metal cylindrical tube which attaches to the mounting rods.

The outputs of the dual colour (visible and infra-red) sensors are routed to the inputs of a dual, one-of-four analog multiplexer 60 whose channel select address is determined by two address lines A0, Al. Two input control signals (~isible and infra-red gain/channel selects) are provided for remote selec-tion of the sensor address. A frequency encode scheme is implemented to select the desired sensor address. The presence of a 10 kHz carrier on a control line is detected by dual chan-nel tone decoder 62, which translates this carrier frequency into a TTL logic level for selecting the multiplexer address.
The outputs of the selected sensor are fed to pre-amplifier and decoupling stages. Pre-amplifiers 64, 66 provide high initial signal gain. An NE570 based compander stage 68, 70 provides further gain amplification with the overall A.C.
gains controlled by voltage controlled gain (VCG) inputs. The VCG section gains are determined remotely via two control inputs. A 60 dB gain/attenuation range is achieved, ensuring no signal saturation over extremes in flame brightness and flicker content.
-~
The outputs of the VCG stage are bandpass filtered to provide a fre~uency sensitive gain characteristic whose gain is proportional to frequency in the range 10 Hz S freq S 500 Hz.
Above 500 Hz the signals are attenuated at -30dB/octave to remove high frequency noise components. The D.C. components of the sensor outputs are fed forward to the second low-pass stage of the filter section to provide flame brightness information. The filter outputs are then buffered and routed to a remote processor ~i.e. computer) over shielded twisted pair cable.
-The scanner electronics can be configured to meet par-ticular gain characteristics by choosing intermediate stage gains 201~0~0 as required. Lead sulphide/silicon, lead selenide/silicon and germanium/silicon dual colour sensors can be accommodated, al-though a single combination is preferred in any one scanner head.

Detailed Circuit Description The design of the dual colour sensor circuit electr-onics is essentially identical for the infra-red and visible channel signal conditioning. The only significant difference is that the visible (silicon sensor) channel incorporates a dual gain mode to accommodate the wide dynamic range experienced when monitoring both coal and oil flames. Both the visible and infra-red circuit are A.C. coupled, with provision made for feeding the D.C. component forward to an output summing stage for monitoring flame intensity.
Pre-amplifier Staae As depicted in Figure 12, outputs DRA, DRB of analog multiplexer U1 are A.C. coupled via capacitors CR0 and CI0 to non-inverting amplifiers U2, U3. The sensor outputs are biased to +V by resistors RRl, RI1, with an optional dual gain mode achieved by zener diode/ resistor pairs ZR0, RR3 and ZI0, RI3, -This secondary gain mode is only operational under very bright conditions, when the zener diodes conduct. Under these condi-tions the sensor outputs are essentially attenuated by the ratios (RR3/RRl), (RI3/RI1).

The pre-amplifier stage gains are determined by feedback resistors RI4, RI2 and RR4, RR2. The pre- amplifier bandwidth is limited to about 1 kHz by feedback capacitors CRl, CI1.

Voltaae Controlled Gain Staqe fVCGl A dual channel NE570 compander integrated circuit U4 provides voltage controlled gain characteristic. Resistor, capacitor pairs RR5, CR2 (infra-red) and RI5, CI2 (visible) together with the variable impedances of the input voltage controlled stages determine the channel gains and low frequency : . . , .:: . . . :, ~ : . :

2 01 .rj ~) t3 ~

A.C. coupled response of compander U~. The inv~rting inputs o~
Compander U4 are configured as summing junctions with overall gain and high frequency roll-off determined by feedback via RR6, CR5 (infra-red) and RI6, CI5 (visible). The bias resistors RR7, RR8 and RI7, RI8 are chosen to minimize D.C. output offsets over the complete controlled gain range.

The gain control voltages are set to VDD + 1. 8 volts for minimum gain, with maximum gain at o volts. Typically, VDD
is in the range of -15V < VDD S -12V. The low-pass filtering provided by RR20, CR15 and RI20, CI15 blocks the 10 kHz carrier signal which may be present on the channel select/gain control inputs. Capacitors CR3, CI3 limit the speed of response in channel gain to changes in the D.C. level of the gain control inputs.

Si~nal Conditionin~ (Pre-emphasis Filter) The outputs of the VCG stages are bandpass filtered.
The filter characteristics are chosen suc~ that gain is approx-imately proportional to frequency in the range 5 Hz < freq < 500 Hz.

The VCG stage outputs are first high-pass filtered by U5 with the high-pass (derivative) mode time constant determined by RR10, CR6 (infra-red) and RI10, CI6 (Visible). The high-pass stage gains are limited by resistors RR9, RI9 and capacitors CR7, CI7. Provision is made for D.C. coupling the sensor outputs directly via resistors RR17, RI17. These resistor values are chosen such that + full scale D.C. output on the sensor results in +2 volt offsets on the outputs of filter U5. The second ~-stage pre-emphasis filter U6 is designed as an under-damped low-pass stage which limits the high frequency response while at the same time providing signal enhancement in the frequency range 250 Hz ~ freq ~ 500 Hz. The damping ratio is determined by capacitor pairs CR8,CR9 (infra-red) and CI8, CI9 (visible).
Overall unity D.C. gain is maintained through the VC5 stages.

201509~

Output Bufferinq The output buffer stages associated with amplifier U7 are configured as inverting buffers with 1 kHz, first order low-pass roll-off. Resistor/capacitor pairs RR17, CR9 (infra-red) and RI17, CI9 (Visible) determine the low-pass time con-stants. Resistors RR1~, RIl5 determine the stage gains.

Analo~ Multi~lexer & Channel Select The 10 kHz carrier frequencies for multiplexer channel select are A.C. coupled via CRll, CIll. The centre frequencies for dual channel tone decoder U8 are set by resistor/capacitor pairs RR21, CR13 (infra-red) and RI21, CI13 (Visible). The bandwidth (i.e. frequency range about the centre frequency in which the tone decoder responds) is determined by capacitors CR14, CI14 and is set to approximately + 500 Hz. The tone decoder outputs provide a 2 bit address select (A0, Al) for multiplexer Ul.

Mode of Operation A remote controller selects the input channel and adjusts the output gain via two gain/channel select input lines.
The intended mode of operation assumes gain and channel select are held constant over a measurement interval which is determined by the flame detection algorithm. If channel selects are changed then time (about 40 milliseconds) must be allowed for the channel outputs to reflect the new signal source values. This time is determined ~y the multiplexer and filter transient decay times.

Similarly, a change in channel gain, initiated by varying the input D.C. control voltage on the appropriate gain input line, results in an exponential response in the overall gain of compander 4 due to smoothing capacitors CR7, CI7.
Reducing the size of these capacitors speeds up the response of the VCG gain sections. However, to retain 60 Hz rejection it is recommended that the gain time constants be > 100 msec, where the time constants are given by:
TGAIN = CGAI~ 10 seconds , , , , , . : ,. . ~ -:

2~1~0~

The dual gain mode capability provided by RR3, ZRO on the infra-red chann~l and RI3, CL0 on the visible channel should be selected such that the circuits operate in mode 1 (high gain, diodes non-conducting) when monitoring coal flames, and in mode 2 (low gain, diodes conducting) when monitoring auxiliary flames fuelled by oil or gas. The selection of zener diode voltage and resistor values is location dependent.

Flame Detection Alqorithm Types of Flame Conditions In general there are four flame conditions or classes of flame to be detected in a multi-burner boiler. These are:
(1~ Main fuel flame from the individual burner being monitored ("MAIN FLAME").
(2) Flame from the auxiliary or igniter burner associated with the main burner ("AUX. FLAME").
(3) Fireball or background flame from other burners ("FIREBALL").
(4) Flame out condition on both the main and auxiliary - -burners ("FLAME OUT").

In most situations an attempt is made to discriminate flame for the particular burners (main and auxiliary) being monitored. This is desirable but not always possible when other burners are present and contributing to the boiler firing state.

Monitorinq Usina Multiple Scanners Although in most situations only one scanner head is required for each burner, there are situations when two scanners may be used to improve flame discrimination. It is also possible for more than one sensor to be mounted in a scanner head. In general each burner flame is characterized by "M" separate data signals, all fed to the same central processor and sampled in parallel to retain their time coherent properties. These M
signals may be obtained from one or more scanner heads, each equipped with one or more multicolour sensors.
.'., ': .

2~1~09~

S~ectral ~stimation and Notation The general case assumes "m" sensor input signals.
These input signals x1(t)..xm(t) are sampled "N" times in each of "k" block periods (k = 1, 2...), each block being of duration "T"
seconds. The ith sample point (i = 0, 1, 2 ... (N-l)) on the ith time signal Xj (j = 1, 2, ... M), in time block Tk, is denoted by xjk(i) -lo The complex discrete Fourier transform of a sampled signal x; for the kth sample block is denoted by:

Xjk[L] = DFT[Xjk(i) ]
N-l ~Ox3k(m) exp~2~iLm/N) where: I'i'' is the complex root of (-1);
"L" L = 0,1,2....N~2-1 is the Lth harmonic component at frequency (L/T) Hz;
"Xjk[L]" is complex; and, "DFT[ ]" is the discrete Fourier transform operator.
Bolded notation is used to denote frequency domain variables.

The discrete auto-power spectrum density estimate for a signal Xj on time interval Tk is given by;

~iik[L] = Xjk[L] Xjk~L] L = 0,1,2....N/2-1.
where the superscript "*" denotes the complex conjugate.

The discrete cross-power spectrum density estimate between signals x; and x; on time interval Tk is given by:

~Jik[L] = XJk~L] Xik[L] L = 0,1,2....N/2-1.

201509~1 The discrete estimate of the modulus squared transfer function between signals x; and x; on time interval Tk is given by:
Hjik[L] = (S j jk[L] S~ jk[L] )/ (S jk[L] 8jk[L] ) where L = 0, l, 2 .... N/2-1.

The discrete estimate of the modulus squared coherency function between signals x; and x; on time interval Tk is given by:
Cj jk[L] = ( (S j jk[L] ) ~ Sj jktL] )/ (Bjk[L] Sjk[L] ) where L = 0, 1, 2 ... N/2~

Estimates may be averaged over ad~acent frequency bands and/or over successive time block intervals. The software employed in the preferred embodiment allows the user to choose up to 9 separate frequency bands for frequency smoothing and to obtain long term time averaged estimates in these frequency bands using an exponential first order averaging factor (digital low pass filter). Two estimates are updated in every time block interval. Firstly, a frequency smoothed estimate is obtained from the last time block interval. This estimate is given by:

Efjlter = 1/M ~ Elast(L) where L is the harmonic number; M = ~-L1+1, ~ 2 L1; and up to nine filters are specified, each with independently specified limits L1, ~. Secondly, a frequency and time averaged estimate is obtained from the last "k" intervals as determined during setup. The time averaged estimate is given by:
-...
Eave = ~SEold + ( 1 - ~5 ) E~a5t where E2Ve is the new averaged estimate; E~ld is the previous averaged estimate; El~st is the latest estimate; and, ~ is the averaging time constant.

The frequency smoothed last block estimate E~ast, and the time averaged estimate EaVe consist of estimates of the auto-spectra, cross-spectra, squared coherency an~ modulus squared transfer gain functions in each of nine frequency bands, f~, L =
o, 1, 2 .~. 8i. These estimates are, in reality, a set of measurements including:

(1) Auto-spectra estimates 8jj[L] of signals j = 1, 2, ... M for filters L = 0, 1, ... 8.
(2) Squared Modulus Transfer function estimates Hjj[L], between signals Xj and Xj; j,i = 1, 2, ...
M.
(3) Squared coherency function estimates Cjj[L]; j,i = 1, 2~ ... M.

where L = 0,1,2,3,4,5,6,7,8 - the preset frequency bands of interest.

The measurements are corrected to account for preset channel gains which are adjusted on the scanner heads prior to commencing each block of time samples. The averaged transfer function and coherency estimates are obtained by first averaging individual estimates of the cross-spectra and the auto-spectra and then dividing the resulting averaged cross-spectra products by the appropriate auto-spectra.

In addition to the long term average and last block estimates, the variance of estimates about the long term average is also calculated as:
Evar = ~iEvar-old + (1 -- ~15) [Eave ~ Elast] .. ' where EVar is the new estimated variance; EVar-o~d is the previous estimated variance; and E~Ve, El~st and ~ are as previously defined. The averaging time constant ~ is chosen such that 0.01 S ~ S 1Ø The standard deviation of estimates is then simply calculated as:
Edev = ~
The variance and/or standard deviation can then be used to detect the onset of unstable flame conditions; usually characterized by 201~090 large flu~tuations about a normal operating point.

scanner Flame Detection se~uence The flame detector is operated in one of three modes:
(1) learn flame signature;
(2) monitor flame; or, (3) self test.

In the first ("learning") mode, the flame detector identifies the statistical properties of spectral estimates and stores these characteristic measurements as being typical of one of the four flame conditions outlined above. The amplitude probability distributions of the spectral estimates, as well as the minimum, maximum, average and variance values of these functions in each of the frequency bands are calculated. These are stored as signatures characteristic of the particular flame conditions.

In the second ("monitoring") mode, the flame detector compares latest flame spectral estimates against prestored flame signature characteristics and outputs a measure of "flame on"
confidence for the main, auxiliary and fireball flame condi-tions. These three "flame on" confidence levels are compared against individual "flame" and "no flame" setpoints to determine the corresponding flame contact output status. The setpoints have a variable dead band characteristic to avoid contact output chatter.

In the third ("self test") mode, known signals are fed to the scanner heads in a loop back mode to check system integrity.

A block overview of the scanner software logic is shown in Figure 13.

201~90 _canner Head Initialization Before c~mmencin~ flame monitoring, the flame detec-tor's co-processor selects the designated sensors in the scanner heads (1 of 4 in each head) and adjusts the sensor gains to achieve good signal to noise levels at the A/D converter. The sensor gains ar~ controlled by varying the output voltages on two D/A channels. These voltages are fed to the voltage controlled gain sections on the scanner head electronics. The sensor selection in each head is achieved by the co-processor transmitt-ing two frequency modulated carrier signals (10 kHz carriers)superimposed on the D.C. gain signals. These signals are decoded by the scanner head electronics as a two bit address for the front end multiplexer. Loop integrity is also checked by transmitting a second carrier at a lower frequency (< 500 Hz) which is then amplified by the scanner head electronics and received on the incoming data channels. The channel gain calibration can be verified as well as overall signal integrity using this secondary carrier. The channel gains are ad~usted to achieve a signal strength of approximately 2.0 volts R.M.S. from the sensor. This ensures good si~nal to noise ratios over the transmission cable, while avoiding saturation problems on the A/D converter. The A/D converter's full scale range is ~10 volts.

Data Acouisition ~-The analog data from the scanner sensors usually consists of two data channels~ x1, x2, corresponding to signals representative of the flame emissions in an infra-red and a visible wavelength band. Up to 4 signals can be accommodated.
This situation arises if:
i) more than one sensor is selected simultaneously; or, ii) a multicoloured sensor as opposed to a dual colour sensor is used (eg: Si/Ge/PbSe); or, iii) a sensor array and chromatic beam splitter are used;
or, iv) more than one scanner head is installed.
:

20150~

The discussion of the flame detection algorithms will be limited, without loss of generality to the bivariate case.
As previously explained, the signals are sampled in blocks of N
sample points, where N is usually chosen to be 2M, consistent with a radix 2 based discrete Fourier transform (DFT). The sample block mean values are calculated and subtracted. These mean levels, or D.C. components, are measures of flame brightness and may be tested as indicative of flame condition in a similar manner to the spectral estimates. The sample blocks are optionally preprocessed using a Hanning time window to suppress side-band leakage inherent in short period DFT analysis (see:
Bendat J.S., Piersol A.G., "Random Data: Analysis and Measure-ment Procedures," Wiley Interscience 1971 Library of Congress #
71-160211).
Siqnal Processing Alqorithms The spectral estimates, as described above, are estimated for the nine selected frequency bands. These bands are arbitrarily chosen and may or may not be contiguous. The spectral outputs, as estimated in these frequency bands-, are termed filter outputs. The only restrictions on the choice of filter characteristics are:

i) For each filter the low-frequency/high-frequency cutoffs must lie in the range - ~ fcutoff < sample frequency/2.
ii) The cutoff frequencies are discrete harmonics of (l/T) Hz where "T", the block sample interval, is the frequency resolution of the DFT analysis.
Both last block spectral estimates and long term time averaged spectral estimates are calculated and updated after each block of time samples has been stored.

In the "learning" mode, the signature maxima, minima, variance, and average values and the individual amplitude probability distribution functions are updated for each of the 201~0g~

spectral estimators (auto-spectra, squared modulus gain and squared coherency) in each of the filter output bands. These values are later saved as signatures indicative of the flame condition being monitored.
S

In the "monitoring" mode, the latest and/or long term average spectral estimates are compared with one or more previously stored signatures. The maximum number of signatures is limited only by the available storage memory and by real time processing constraints. Each comparison yields a probability match figure in the range of 0 < match < 1Ø The best match obtained for each of the three flame types (main flame, auxiliary flame and fireball) is used as an indication of the respective flame status. Thus, several signatures indicative of main flame may be tested and the best fit used for signalling the main flame status. If the flame "match" is greater than the "FLAME-ON"
setpoint for that type of flame the flame status is signalled "ON". If the flame "match" is less than the "FLAME-OFF" setpoint .then the flame status is signalled "OFF". If the "match" is between the "FLAME-ON" and "FLAME-OFF" setpoints the flame status remains unchanged. Initially flame status is signalled "OFF".

Flame Condition Contacts Four flame contact output relays are provided. These 5 are:
- main flame status - auxiliary flame status - fireball status - online/offline status In addition, four contact inputs are provided. These are usually designated: -:
- main fuel status - auxiliary burner fuel status - master enable/disable - self test.

201~0gO

The status of the contacts is updated after every block of data samples and after every test of flame condition.

Scanner Gain Adlust and Cell Selection The flame detector channel gains are updated after each block of data is sampled. The gains are calculated based on the signal variances measured in the previous block of samples. The channel gains are maintained constant during block sampling to avoid bias errors occurring in the spectral estimates. Similar-ly, the scanner sensor selection may be updated between sample block intervals, to better locate the position of the primary combustion zone of the burner flame. The scanner tries to locate the flame using the sensor with the viewing window closest to the burner nozzle. Where multiple sensors are installed, failure to find flame close to the burner will result in the selection of the next appropriate sensor as determined by the user prior to commencing scanning. The sensor selection sequence may be determined by spatial considerations and by contact input fuel status information. The latter is appropriate where the igniter or auxiliary burner has a very different flame pattern from the main burner and requires the use of a different viewing window to improve flame discrimination. When monitoring more than two analog channels (which allows the simultaneous monitoring of more -than one sensor), the only restrictions are the number of analog channel inputs provided (four are provided in the preferred embodiment herein described) and the real time processing delay incurred by the estimation of spectral filter outputs on multiple channels.
Learnina Flame Siqnatures Each flame condition is characterized by a spectral flame signature measured in terms of the:

- maximum filter output;
- minimum filter output;
- average filter output;

20150gO

- variance and/or standard deviation of the ~ilter output about the average; and, - amplitude probability distribution.

Each spectral function in each filter band is characterized in this manner. The flame condition can be representative of a particular firing condition or a range of firing conditions such as might be encountered by varying firing air flow or fuel flow.
Particular flame conditions of interest can be singled out if necessary to provide better flame discrimination.

Siqnature Classifications Flame signatures are classified as being indicative of one of four flame conditions:
i) Main burner flame; or, ii) Auxiliary burner flame; or, iii) Fireball flame; or, iv) Flame out.
The main burner flame is the flame associated with primary fuel burner. The auxiliary burner flame is the flame associated with the igniter or secondary burner. The fireball flame is any flame whose characteristics cannot be attributed purely to the burners being monitored. Other burners may contribute to the fireball characteristics. Flame out conditions are characterized by the absence of any of the first three flame conditions. Unfortunate-ly, the one flame condition that is of interest, must be avoided (i.e. flame out with fuel still being supplied to the burner).
This condition is not available for classification in terms of a flame out signature, as operation of the boiler under these conditions constitutes a safety hazard. Several signatures of each type of flame may be required to completely characterize the normal firing situations on the burners.

Multiple Siqnature Testinq When more than one flame signature is used to test for the flame condition, the latest spectral estimates and/or time averaged estimates are matched against each signature in turn.

201~090 The best "fit" for each flame type is return~d as the flam~
condition for that flame type. ~lowever, the probability of any flame type being "ON" is constrained to be less than (1.0 -probability of flame out) as determined by matching the flame spectral estimates against all flame out signatures. This ensures contradictory flame condition indications err on the side of safety.

Monitorin~ Burner Flames There are many ways to compare the latest spectral estimates of a flame output with previously stored characteristic signatures. However, only a limited number of comparison techniques lead to robust flame detection algorithms. As noted above, several flame signatures may be compared to detect a specific flame status. Each comparison involves deriving a measure of fit between a current measure of the flame and a previously stored flame signature. By convention it is con-venient to express this measure of fit as a normalized probabili-ty in the range 0 < prob < 1Ø The various comparison algorithms will now be described in detail.

Siqnature ComParison Algorithms The spectral estimation algorithms and the methods used to obtain smoothed estimates of the auto-spectra, transfer gains and coherence in each of nine spectral bands were described above. The discussion was presented for the general case of several daka signals. The information gathered ~hen learning a flame signature has also been described above. In particular average, maximum, minimum and the variance of spectral estimates are recorded together with the individual amplitude probability distributions of each of the estimates. Strictly speaking, a true test of the measure of match of an estimate vector X (- R~) with a previously stored signature vector Z (- Rn) requires knowledge of the n dimensional joint probability distribution, p(Z), of estimates of Z. This is only equal to the product of the individual probability distributions, p(zl).p(z2)...p(zn) (where zl..zn are members of Z), if the estimates zi (i=l..n) are 2015~90 statistically independent. It is then a simple matter to retrieve the probability of an estimate X being representative of the set Z, from the ~oint distribution characteristic previously stored. The probability is normalized to the maximum probability, Pmax, and the returned measure of fit is then in the required form. Unfortunately, the data memory storage require-ments needed to approximate the true ~oint probability distribu-tion sufficiently precisely make this approach unrealistic for large dimension, n, vectors. Furthermore, the resulting flame detection algorithm tends not to be robust.

Four alternative comparison techniques have been developed. The choice of comparison technique to apply for each signature test is made by the operator and is stored as part of the information contained with each signature record. The operator is not limited to using just one technique for all the estimates in a signature test. Each spectral estimate in each filter band may be tested using any of the four methods proposed.
The overall measure of fit is then obtained as a weighted averaged of all the individual measures of fit.

Weiqhted Least Squares Fit An estimate, x, is compared with a signature value, z, as follows: 5 Given signature values for z of:
Zave = average signature value of z.
zmin = minimum signature value of z.
zmax = maximum signature value of z.
0 Case 1 == Zm~n < X < zave e = (x - zave)2.
emax = (zmin - zave) .
p(x) = 1.0 - e/emax. : Probability of estimate x being a measure of the signature z.
Case 2 == zave ~ x < zmax e = (x - z ave ) .

: .

201~09~

emax = ( zmax -- zavc) 2 .
p(x) = l.o - e/emax.: Probability of estimate x being a measure of the signature z.

Case 3 == ~x < Zmin) or (x > zmax) p(x) = o.o The returned probability of fit is just a measure of the distance squared between the estimate x and the average 10 signature value, Zave. If the estimate x is less than the minimum value of z or greater than the maximum, then a zero probability of fit is returned. The lower and upper bound limits are usually those found by experiment, but they may be replaced or forced to other values to improve the test response where this can be 15 justified. As an example, if there is no penalty required if a measure of the auto-spectrum of a flame signal for a particular f ilter exceeds the average value, then the previously measured maximum limit can be replaced 3:y a very large value so that all estimates that exceed the average return an approximate f it 20 probability of 1Ø Similarly, the lower minimum limit might be replaced if the test is to determine a flame out characteristic, where lower amplitude estimates indicate a darker boiler with less background f lame .

25 Stochastic Fit Given a discrete approximation, p(z), to the amplitude probability distribution of a signature value z, the normalized probability of obtaining an estimate, x, is obtained as follows.
Let the discrete representation of p(z~ be denoted by the set of 0 points p ~ i ], i=0, 1, 2, . ~. n , where pti] = prob ~z, zmin + i ~ < z < zmin + (i+l)-~) ( zmax-zmin) /n == discretization resolution and ~ p[i] - 1. 0 == total probability normalized to unity.

35 The normalized probability of obtaining an estimate, x, is just (ptj ]/pmax~, where the index j is given by:
= int ( (x -- zm~n)/~) 20~090 and pmax = max (p[i], i = O..n).

If j is less than zero or greater than n then the probability is assumed to be zero. The probability distributions associated with individual signature spectral estimates are calculated during the learn flame operating mode.

Bounded Limits A third method of obtaining a measure of the fit between an estimate, x, and the signature value z is obtained by using a similar test to the least squares method described above, except that the weighting function is no longer based on the squared error law. The general formulation can be presented as follows. As before, given signature values for z of:
zave = average signat~re value of z.
Zmin = minimum signature value of zO
zmax = maximum signature value of z.

Case 1 == Zmin < X < zave let e = ¦x - zave¦.
emax . = ¦ Zmin ~ Zave¦ .
p (x) = 1. o - f(e/emax) : Probability of estimate x ~ --being a measure of the signature z;
25 where the function f ( ) is defined for all values in the range 0..1 and is normalized such that:
f(O) = o f(l) = 1 :
o < f( ) < l.o ~ ~ -Case 2 == zave < x < zmax e = ¦ X ~ Zave ¦ .
emax = ¦zmax - zave¦
ptx) = l.0 - f(e/emax) : Probability of estimate x being a measure of the signature z.

:
- 47 - ;

201~i091) Case 3 == (x < z~lin) or (x > Zm~x) p(x) = o.o The functions f( ) implemented in the preferred embodiment are:
i) Uniform weighting [f(x) = l.o (o < x < 1.0)];
ii) Square law (as given above under "Weighted Least Squares Fit");
iii) Triangular weighting ~f(x) = x; (0 < x < 1.0)]; and, iv) Cubic weighting [f (x) = X3] .
As discussed above, the zmin and z~x limits can be artificially extended to give a one sided limit test if required.

Gaussian The last test assumes a Gaussian probability distribu-tion for estimates, x, of the signature z, with standard deviation Zdev and mean value Zave. The standard deviation is calculated during the "learning" mode for each of the spectral estimates. In this case the probability of obtaining x is given by:
p(x) = exp((x - zave)2 / Zdev2) The probability is normalized to the maximum probability p(0).
The assumption of a Gaussian distribution is justified for auto-spectra estimates which are averaged over adjacent frequency points or sequential time blocks where the number of points used for averaging is large ( > 20). For estimates of the modulus squared transfer gain and squared coherence it can be shown that the distribution of log(x) for these functions is more nearly Gaussian than the distribution of x itself (see: Bendat J.S., Piersol A.G., "Random Data: Analysis and Measurement Procedures"
s~pra). The operator is given the option of testing x or log(x) for these functions.

Calculation of Estimate Weiqhtinq Functions The measures of fit returned for each of the individual spectral estimates are summed and averaged to obtain an overall 2015~0 fit probability for each signature to be tested. The overall fit probability is given by:

pfit = ~ W[i] p(X[i] ) for all estimates x[i].

The weights, w[i], are normalized so that ~ w[i] = 1Ø The choice of the weighting function determines how much importance is given, in relative terms, to the auto-spectra, transfer gain and coherence estimate errors for each filter output. Where a particular signature average estimate, Z[i]ave, for a flame "ON"
condition is very different from all measures of flame "OFF" for that estimate, the assigned weight is correspondingly large.
When the flame i'oN" to flame "OFF" difference is small the weight attached is small.
Given flame "OFF" signatures Z[j], j=0,1,2.. m, and a -flame "ON" signature X, the flame "ON" to flame "OFF" distance for an estimate x[i], (x[i] - X), is given by:

d[i] = min { ¦x[i]ave - z[i]ave¦ / zLi]dev ~

where (.Ave) denotes the average signature value and (.deY) denotes the signature standard deviation. The total distance measure d~x is determined as:
d~x = ~ d[i] - for all spectral functions x[i].

The signature weighting functions w[i] are either calculated as:

w[i] = (d[i]/d~x) or as w[i~ = (d[i]/dmax)2 ':
depending on whether a modulus or square law weighting is required. The operator may override the weighting function for any or all spectral estimates if required and manually set alternative values.
. .

2~1S~9~
Measure of Fl me Match Signatures are obtained for four flame conditions as explained above. Each flame condition may be characterized by one or more signatures. The probability of an estimate, X, belonging to a particular flame type is given by:

Probability of flame type = max ~ pfit [ i] ) where pfit[i] = Prob ~ X is an estimate of signature Z[i] } for all signatures Z[i] which characterize the flame type tested.

In other words, the best fit is considered to be the probability of a particular flame type. ~owever, as noted above, notwithstanding the probabilities obtained above, the maximum probability of any type of flame "ON" condition is constrained to be less than or equal to (1.0 - max. probability of flame "OFF"). This ensures flame "OFF" takes precedence over flame "ON" and that conflicts result in a flame "OFF" condition being signalled.
Orthoqonality Test Modifier Each spectral estimate X is in fact an 'n' vector.
Similarly, the averaged signature estimates Zave~i] are also 'n' vectors. A simple measure of the cosine of the solid angle ei between the estimate X and each signature Zave[i] is obtained by taking the dot vector product as follows:
X ~ Zave[ i]
cos (ei) r-3 0 `¦ ¦ ¦ X ¦ ¦ ¦ ¦ Zave [ i ] ¦ ¦ :
where "¦¦ ¦¦" denotes the euclidean norm, a scalar product representing the magnitude squared of the vector. The above estimate has the desired property 0 S cos (ei) s 1 . o, and indicates the degree to which the estimate X is orthogonal to 35 each signature Zave[i]. The overall measure of fit as obtained and used in the "probability of flame type" equation given above, is modified by this orthogonality factor to enforce a stricter :, , .. ~ ~. . . , , . -201~0~

classification on X than is obtained purely by using the tests previously outlined.

SummarY
The method of flame detection herein described depends on characterization of the different flame conditions in terms of characteristic spectral signatures; and on the calculation of a weighted measure of fit between a latest spectral estimate and previously stored signatures. No knowledge of the burner condition is assumed except that available from the spectral estimation process.

The procedure followed to detect flame condition is as follows:
i) Select spectral filter characteristics suitable for monitoring the flame and fuel type.
ii) Select the burner operating range and conditions for which flame is to be detected.
iii) Obtain the signatures characteristic of all the flame conditions to be monitored.
iv) Select the type of comparison tests to be used with each signature.
v) Calculate the weighting function associated with each signature.
vi) Select the spectral estimate averaging mode to be used.
vii) Collect blocks of sample points of the relevant data channels and estimate the spectral function outputs from these data.
viii) Compare the latest spectral function outputs with the previously stored signatures and obtain a measure of flame fit.
ix) Output flame condition.
x) Repeat steps (vii) through (ix).

The filter characteristics of the spectral functions may be chosen arbitrarily as low-pass, high-pass, or, bandpass with overlap between different filters if desired. The only restric-tions on the choice of filter corner fre~uencies are those imposed by data sampling rates and the number of samples in each data block. The sampling rate should be chosen to be greater than twice the frequency of the highest frequency component in the data signals. The tests for flame "fit" may be conducted using the last block frequency smoothed estimates and/or the time averaged estimates.

The methods for estimating the measure of fit can be readily extended to include alternative algorithms such as:

- maximum likelihood - minimax fit criteria - true joint amplitude probability testing as outlined above.
- various function fit criteria along the lines of those discussed above under the heading "Bounded Limits".

As will be apparent to those skilled in the art in the light of the foregoing disclosure, many alterations and modi-fications are possible in the practice of this invention without departing from the spirit or scope thereof. Accordingly, the scope of the invention is to be construed in accordance with the substance defined by the following claims.

LIST OF YMBOLS
dt The perpendicular distance from the optical axis to the top of the viewing window.

db The perpendicular distance from the optical axis to the bottom of the viewing window.

D The distance along the optical axis (O) in front of the lens from the principal point (Pp) to the viewing window.

f The focal length of the lens at any given wav~length. All parallel rays entering the front of the lens will come to a focus this distance behind the lens.

15 fO The nominal design focal length of the lens at a specific wavelength (~O) specified by the manufacturer.

L The distance along the principal axis (PA) in front of the lens from the principal point (Pp) to any arbitrary loca-tion.
.
n The actual index of refraction of the lens at any given wavelength (~).

25 nO The nominal design index of refraction of the lens at a specific wavelength (~O) specified by the manufacturer.

O The optical axis is the path traced by a ray intersecting any point (P), usually the mid point of a finite sensor, and the principal point (Pp).

P A point located at an arbitrary location behind the lens. ~;

PA The principal axis is the axis of symmetry passing through the centre of a circular lens.

2015~'30 Pb The point of intersection on the focal plane n~ade by a ray intersecting the bottom edge of the lens and any arbitrary point ~P).

5 Pf The focal point is the point on the principal axis ( PA) where all rays of any given wavelength (~) entering the front of the lens parallel to this axis come to a focus behind the lens.
0 Pp The principal point is the point where the principal plane is intersected by the principal axis (PA)' Light rays passing through this point are not diffractedO

Pt The point of intersection on the focal plane made by a ray intersecting the top edge of the lens at any arbitrary point (P)-R The distance from the principal point (Pp) to any arbitrary point (P).
rO The actual radius of curvature of a lens.

r1 The radius of curvature of the front surface of a convex lens.
r2 The radius of curvature of the rear surface of a convex lens.

rb The path of the light ray passing through the bottom of the lens and through the principal point (Pp).

rt The path of the light ray passing through the top of the lens and through the principal point (Pp).
5 rm The path of the light ray passing through the middle of the lens and through the principal point (Pp).

tc The centre thickness of the lens.

te The edge thickness of the lens.

X The distance along the principal axis (PA) behind the lens from the principal point (Pp) to any arbitrary location.

Y The perpendicular distance from the principal axis (PA) to any arbitrary location behind the lens. "Y" is positive when above the principal axis ( PA) and negative when below.

Ym The perpendicular distance from the principal axis (PA) to the mid point of a sensor surface.

YVjs The distance from the visible range sensor mid point to its edge.

YIR The distance from the infra-red range sensor mid point to its edge.
The angle that any light ray makes relative to the principal axis (PA). Behind the principal point (Pp) "~" is positive above the principal axis ( PA) and negative below. The re-verse is true in front of the principal point (Pp).

b The angle at which the bottom light ray (rb) enters the front of the lens.

~m The angle at which the sensor mid point light ray (rm) passes through the principal point (Pp).

~t The angle at which the top light ray (rt) enters the front of the lens.

5bmax The maximum angle that a ray entering the bottom of the lens will impinge on the surface of a finite sensor.

201~09~

G~in The minimum angle that a ray entering the bottom of the lens will impinge on the surface of a finite sensor.

~tmax The maximum angle that a ray entering the top of the lens will impinge on the surface of a finite sensor.

~tmin The minimum angle that a ray entering the top of the lens will impinge on the surface of a finite sensor.
O ~ The nominal lens aperture diameter. This is usually assumed to be egual to the actual lens diameter.

The nominal design wavelength of a lens as specified by the manufacturer.
The actual wavelength.

Claims (52)

1. A method of detecting flame within a region, comprising the steps of:
(a) measuring radiation emitted from said region within a selected portion of a visible frequency band;
(b) concurrently measuring radiation emitted from said region within a selected portion of an infra-red frequency band;
(c) deriving the coherency between said measurements;
(d) comparing said coherency with a prestored coherency signature representative of the coherency between measurements of radiation emitted from said region within said selected portions of said visible and infra-red frequency bands while known flame conditions prevail within said region, thereby estimating the deviation of said derived coherency from said prestored coherency signature; and, (e) comparing said deviation with a first predetermined threshold alarm value.

2. A method as defined in claim 1, further comprising:
(a) deriving the auto spectrum of said visible frequency band measurements;
(b) comparing said visible measurement auto spectrum with a prestored auto spectrum signature representative of the auto spectrum between measurements of radiation emitted from said region within said selected portion of said visible frequency band while known flame conditions prevail within said region, thereby estimating the deviation of said derived visible measurement auto spectrum from said prestored visible auto spectrum signature; and, (c) comparing said deviation with a second predetermined threshold alarm value.

3. A method as defined in claim 2, further comprising:
(a) deriving the auto spectrum of said infra-red frequency band measurements;
(b) comparing said infra-red measurement auto spectrum with a prestored auto spectrum signature representative of the auto spectrum between measurements of radiation emitted from said region within said selected portion of said infra-red frequency band while known flame conditions prevail within said region, thereby estimating the deviation of said derived infra-red measurement auto spectrum from said prestored infra-red auto spectrum signature; and, (c) comparing said deviation with a third predetermined threshold alarm value.

4. A method as defined in claim 3, further comprising:
(a) deriving the transfer function of said visible and infra-red frequency band measurements;
(b) comparing said transfer function with a prestored transfer function signature representative of the transfer function between measurements of radiation emitted from said region within said selected portions of said visible and infra-red frequency bands while known flame conditions prevail within said region, thereby estimating the deviation of said derived transfer function from said prestored transfer function signature; and, (c) comparing said deviation with a fourth predetermined threshold alarm value.

5. A method as defined in claim 1, further comprising repeating said measurements for other portions of said visible and infra-red frequency bands.

6. A method as defined in claim 2, further comprising repeating said measurements for other portions of said visible and infra-red frequency bands.

7. A method as defined in claim 3, further comprising repeating said measurements for other portions of said visible and infra-red frequency bands.

8. A method as defined in claim 4, further comprising repeating said measurements for other portions of said visible and infra-red frequency bands.

9. A method as defined in claim 5, wherein said coherency comparing step comprises applying a weighted least squares fit to said derived coherency and said prestored coherency signature.

10. A method as defined in claim 6, wherein said visible auto spectrum comparing step comprises applying a weighted least squares fit to said derived visible auto spectrum and said prestored visible auto spectrum signature.

11. A method as defined in claim 7, wherein said infra-red auto spectrum comparing step comprises applying a weighted least squares fit to said derived infra-red auto spectrum and said prestored infra-red auto spectrum signature.

12. A method as defined in claim 8, wherein said transfer function comparing step comprises applying a weighted least squares fit to said derived transfer function and said prestored transfer function signature.

13. A method as defined in claim 5, wherein said coherency comparing step comprises applying a stochastic fit to said derived coherency and said prestored coherency signature.

14. A method as defined in claim 6, wherein said visible auto spectrum comparing step comprises applying a stochastic fit to said derived visible auto spectrum and said prestored visible auto spectrum signature.

15. A method as defined in claim 7, wherein said infra-red auto spectrum comparing step comprises applying a stochastic fit to said derived infra-red auto spectrum and said prestored infra-red auto spectrum signature.

16. A method as defined in claim 8, wherein said transfer function comparing step comprises applying a stochastic fit to said derived transfer function and said prestored transfer function signature.

17. A method as defined in claim 5, wherein said coherency comparing step comprises applying a bounded limits fit to said derived coherency and said prestored coherency signature.

18. A method as defined in claim 6, wherein said visible auto spectrum comparing step comprises applying a bounded limits fit to said derived visible auto spectrum and said prestored visible auto spectrum signature.

19. A method as defined in claim 7, wherein said infra-red auto spectrum comparing step comprises applying a bounded limits fit to said derived infra-red auto spectrum and said prestored infra-red auto spectrum signature.

20. A method as defined in claim 5, wherein said coherency comparing step comprises applying a Gaussian fit to said derived coherency and said prestored coherency signature.

21. A method as defined in claim 6, wherein said visible auto spectrum comparing step comprises applying a Gaussian fit to said derived visible auto spectrum and said prestored visible auto spectrum signature.

22. A method as defined in claim 7, wherein said infra-red auto spectrum comparing step comprises applying a Gaussian fit to said derived infra-red auto spectrum and said prestored infra-red auto spectrum signature.

23. A method as defined in claim 8, wherein said transfer function comparing step comprises applying a Gaussian fit to said derived transfer function and said prestored transfer function signature.

24. A method as defined in claim 5, further comprising:
(a) weighting said coherency deviation estimates;
(b) summing said weighted estimates;
(c) averaging said summed, weighted estimates; and, (d) normalizing said averaged, summed, weighted estimates.

25. A method as defined in claim 6, further comprising:
(a) weighting said visible auto spectrum deviation estimates;
(b) summing said weighted estimates;
(c) averaging said summed, weighted estimates; and, (d) normalizing said averaged, summed, weighted estimates.

26. A method as defined in claim 7, further comprising:
(a) weighting said infra-red auto spectrum deviation estimates;
(b) summing said weighted estimates;
(c) averaging said summed, weighted estimates; and, (d) normalizing said averaged, summed, weighted estimates.

27. A method as defined in claim 8, further comprising:
(a) weighting said transfer function deviation estimates;
(b) summing said weighted estimates;
(c) averaging said summed, weighted estimates; and, (d) normalizing said averaged, summed, weighted estimates.

28. A method as defined in claim 1, further comprising repeat-ing said coherency comparing step with prestored coherency signatures representative of selected flame conditions.

29. A method as defined in claim 2, further comprising repeat-ing said visible auto spectrum comparing step with prestored coherency signatures representative of selected flame conditions.

30. A method as defined in claim 3, further comprising repeat-ing said infra-red auto spectrum comparing step with prestored infra-red auto spectrum signatures representative of selected flame conditions.

31. A method as defined in claim 4, further comprising repeat-ing said transfer function comparing step with prestored transfer function signatures representative of selected flame conditions.

32. A method of detecting flame within a region, comprising the steps of:
(a) deriving "m" data signals xi(t), where i = 1, 2, ... m, each of said data signals characterizing radiation emitted from said region within a corresponding portion of the visible frequency band;
(b) sampling each of said data signals "N" times in each of "k" block periods, where "k" is an integer, each block having duration "T" seconds, to derive a plurality of "c" signal samples characterized by Xik (C);
(c) deriving the discrete auto-power spectrum density estimate Siik[L] for xi(t) as:
Siik[L] = Xik[L]* ? Xik[L]
where: Xik[L] is the complex discrete Fourier transform of sampled signal xi(t) for the kth sample block;
"L" = O, 1, ... N/2-1 is the Lth harmonic component at frequency (L/T) Hertz; and, "*" denotes the complex conjugate;

(d) comparing said auto-power spectrum density estimate with a prestored auto-power spectrum density signature representative of the auto-power spectrum density between measurements of radiation emitted from said region within said corresponding portions of said visible frequency band while known flame conditions prevail within said region, thereby estimating the deviation of said derived auto-power spectrum density from said prestored auto-power spectrum density signature; and, (e) comparing said deviation with a first predetermined threshold alarm value.

33. A method as defined in claim 32, further comprising:
(a) while deriving said data signals xi(t), concurrently deriving "m" data signals xj(t), where j = 1, 2, ...
m, each of said data signals xj(t) characterizing radiation emitted from said region within a correspond-ing portion of the infra-red frequency band;
(b) sampling each of said data signals xj(t) "N" times in each of "k" block periods, where "k" is an integer, each block having duration "T" seconds, to derive a plurality of "c" signal samples characterized by xak(c), where a = 1, 2, ... m;
(c) deriving the discrete auto-power spectrum density estimate Sjjk[L] for Xj(t) as:
Sjjk[L] = Xjk[L]* ? Xjk[L]
where Xjk[L] is the complex discrete Fourier transform of sampled signal xj(t) for the kth sample block;
(e) comparing said auto-power spectrum density estimate Sjjk[L] with a prestored auto-power spectrum density signature representative of the auto-power spectrum density between measurements of radiation emitted from said region within said corresponding portions of said infra-red frequency band while known flame conditions prevail within said region, thereby estimating the deviation of said derived auto-power spectrum density Sjjk[L] from said prestored auto-power spectrum density signature; and, (e) comparing said deviation with a second predetermined threshold alarm value.

34. A method as defined in claim 33, further comprising:
(a) deriving the discrete modulus squared transfer function estimate Hjjk[L] for xi(t) and xj(t) as:
Hjik[L] (Sjik[L]* ? Sjik[L])/(Sjk[L] ? Sijk[L]) (b) comparing said discrete modulus squared transfer function estimate with a prestored discrete modulus squared transfer function signature representative of the discrete modulus squared transfer function between measurements of radiation emitted from said region within said corresponding portions of said visible and infra-red frequency bands while known flame conditions prevail within said region, thereby estimating the deviation of said derived discrete modulus squared transfer function from said prestored discrete modulus squared transfer function signature; and, (c) comparing said deviation with a third predetermined threshold alarm value.

35. A method as defined in claim 34, further comprising-(a) deriving the discrete modulus squared coherency function estimate Cjik[L] for xi(t) and xj(t) as:
Cjik[L] = ((Sjik[L]) ? Sjik[L])/(Sjk[L] ? Sik[L]) (b) comparing said discrete modulus squared coherency function estimate with a prestored discrete modulus squared coherency function signature representative of the discrete modulus squared coherency function between measurements of radiation emitted from said region within said corresponding portions of said visible and infra-red frequency bands while known flame conditions prevail within said region, thereby estimating the deviation of said derived discrete modulus squared coherency function from said prestored discrete modulus squared coherency function signature; and, (c) comparing said deviation with a fourth predetermined threshold alarm value.

36. A method as defined in claim 35, further comprising repeating said data signal derivation steps within a plurality of portions of said visible and infra-red frequency bands and then repeating said comparison steps for each of said derived data signals.

37. A method as defined in claim 36, wherein said prestored signatures "z" comprise the average "zave", minimum "zmin", and maximum "zmax", values of spectral estimates derived for said known flame conditions, said method further comprising deriving the probability "p(x)" that each of said derived estimates "x" is a measure of the corresponding prestored signature "z", as follows:
(a) if zmin < x < zave:
deriving e = (x - zave)2;
deriving emax = (zmin - zave)2;
deriving p(x) = 1.0 - e/emax; or, (b) if zave < x < zmax:
deriving e = (x - zave) 2;
deriving emax = (zmax - zave)2;
deriving p(x) = 1.0 - e/emax; or, (c) if (x < zmin) or (x > zmax):
setting p(x) = 0Ø

38. A method as defined in claim 36, wherein said prestored signatures "z" comprise the average "zave", minimum "zmin", and maximum "zmax" values of spectral estimates derived for said known flame conditions, said method further comprising deriving the normalized probability "(p[j]/pmax)" that each of said derived estimates "x" is a measure of the cor-responding prestored signature "z", where:
p[i] = prob (z, zmin + i?.delta. < z < zmin + (i+1)?.delta.);

j = int ((x - Zmin)/.delta.) .delta. = (Zmax-Zmin)/n;
.SIGMA. p[i] = 1.0; and, pmax = max (p[i], i = 0..n).

39. A method as defined in claim 36, wherein said prestored signatures "z" comprise the average "Zave", minimum "Zmin", and maximum "Zmax", values of spectral estimates derived for said known flame conditions, said method further comprising deriving the probability "p(x)" that each of said derived estimates "x" is a measure of the corresponding prestored signature "z", as follows:
(a) if Zmin < x < zave:
deriving e = ¦x - Zave¦;
deriving emax = ¦Zmin - Zave¦;
deriving p(x) = 1.0 - ?(e/emax); or, (b) if Zave < X < Zmax:
deriving e = ¦X - Zave¦;
deriving emax = ¦Zmax - Zave¦;
deriving p(x) = 1.0 - ?(e/emax); or, (c) if (X < Zmin) or (X > Zmax):
setting p(x) = 0.0;
where the function ?(?) is defined for all values in the range 0, ..., 1 and is normalized such that:
?(O) = o;
?(1) = l; and, O ? ?(?) ? 1Ø

40. A method as defined in claim 39, wherein said function ?(?) is a uniform weighting function [?(x) = 1.0 (O < x < 1.0)].

41. A method as defined in claim 39, wherein said function ?(?) is a triangular weighting function [?(x) = x; (O < x <
1.O)].

42. A method as defined in claim 39, wherein said function ?(?) is a cubic weighting function [?(x) = X3].

43. A method as defined in claim 36, wherein said prestored signatures "z" comprise the average "zave" standard deviation "zdev" values of spectral estimates derived for said known flame conditions, said method further comprising deriving the probability "p(x)" that each of said derived estimates "x" is a measure of the corresponding prestored signature "z", as follows:
p(x) = exp((x - zave)2 / zdev2).

44. A method as defined in claim 37, further comprising weighting, summing and normalizing said probabilities to obtain an overall fit probability "pfit", where:
pfit = .SIGMA. w[i]?p(x[i]), for all estimates x[i]; and, .SIGMA. w[i] = 1Ø

45. A method as defined in claim 38, further comprising weighting, summing and normalizing said probabilities to obtain an overall fit probability "pfit", where:
pfit = .SIGMA. w[i]?p(x[i]), for all estimates x[i]; and, .SIGMA. w[i] = 1Ø

46. A method as defined in claim 39, further comprising weighting, summing and normalizing said probabilities to obtain an overall fit probability "pfit", where:
pfit = .SIGMA. w[i]?p(x[i]), for all estimates x[i]; and, .SIGMA. w[i] = 1Ø

47. A method as defined in claim 43, further comprising weighting, summing and normalizing said probabilities to obtain an overall fit probability "pfit", where:
pfit = .SIGMA. W[i]?p(x[i]), for all estimates x[i]; and, .SIGMA. w[i] = 1Ø

48. A method as defined in claim 44, wherein said deviation comparison step comprises deriving said deviation "d[i]" as:
d[i] = min ( ¦x[i]ave - z[i]ave¦ / z[i]dev ) where: Z[j], j=0,1,2...m, is a predetermined flame "off"
signature;
X is a predetermined flame "on" signature;
[i]ave denotes the average signature value; and, [i]dev denotes the signature standard deviation.

49. A method as defined in claim 48, wherein:
(a) said signature weighting function w[i] =
(d[i]/dmax); and, (b) dmax = .SIGMA. d[i] for all spectral functions x[i].

50. A method as defined in claim 48, wherein:
(a) said signature weighting function w[i] =
(d[i]/dmax)2; and, (b) dmax = .SIGMA. d[i] for all spectral functions x[i].

51. Flame detection apparatus, comprising:
(a) a first sensor for sensing radiation emissions in the visible frequency band;
(b) a second sensor for sensing radiation emissions in the infra-red frequency band; and, (c) a lens;
wherein:
(i) said first and second sensors lie coaxially on an axis passing through the principal point of said lens; and, (ii) said sensors' position relative to said lens, said sensors' shape and size, said sensors' material, said lens size and said lens material are selected such that said first sensor receives radiation emissions from a first zone and said second sensor receives radiation emissions from a second zone, said first and second zones coinciding within a selected region and diverging outside said region.

52. Flame detection apparatus as defined in claim 51, further comprising optical filter means for filtering light passing from said lens to said first and second sensors, thereby selectably limiting said sensors to sensing radiation emissions in selected portions of said visible and infra-red frequency bands respectively.