FILTER SELECTION AND GAS CONCENTRATION COMPUTATION FOR SHUTTERLESS OPTICALLY STABILIZED CAPNOGRAPH
CROSS-REFERENCE TO RELATED APPLICATIONS This application is a continuation-in-part of our U.S. Patent Application Serial No. 07/522,208, filed on May 11, 1990, still pending, which is, in turn, a continuation-in- part of our U.S. Patent Application Serial No. 07/401,952, filed on September 1, 1989, now abandoned.
BACKGROUND OF THE INVENTION Field of the Invention
The present invention relates to methods and apparatus for selecting filters and for computing the concentration of the constituents of the respiratory gases of a patient measured using an infrared (IR) gas analyzer. Most preferably, the apparatus of the present invention is employed to select optical filters which allow some degree of channel overlap with minimum cross-talk and to calculate the end tidal C02 and N20 concentrations of a patient.
Description of the Prior Art It is frequently of critical importance to monitor the concentration of carbon dioxide (C02) in the gases inspired and expired from a patient under anesthesia, for expired C02 concentration is a reliable indicator of the carbon dioxide concentration in the arterial blood. In a clinical setting, monitoring expired C02 prevents malfunctions in anesthesia rebreathing apparatus from going undetected and delivering excessive amounts of C02 to the patient.
Rebreathing of anesthetic gases is very cost effective and environmentally desirable, but accurate C02 concentrations are difficult to maintain in the patient circuit without a concentration monitor. It is known that by directing infrared radiation through a sample of a gaseous mixture and measuring the incident radiation illuminating a detecting device, a measure of the infrared absorption of the gas can be obtained. Electrical signals produced by a detecting device in such a manner are thus indicative of the infrared absorption of the gas and can be processed to produce an output indicating the concentration of one or more of the constituents of the gas being analyzed. This type of gas analyzer operates on the principle that various gases exhibit substantially increased absorption characteristics at specific wavelengths in the infrared spectrum and that higher gas concentrations exhibit proportionally greater absorption.
Prior art infrared gas analyzers such as that described in U.S. Patent No. 4,648,396 to Raemer utilize thermopile detectors to analyze gas concentrations. A thermopile detector is comprised of a number of thermocouples. Thermal detectors such as thermopiles are used primarily for infrared sensing and respond to the total incident energy illuminating them. The number of thermocouples must be sufficient to develop enough voltage to develop a suitable signal to noise ratio. A thermistor is used to calculate the Seebeck coefficient, which relates the voltage developed by the thermopile to the temperature differential between the "hot" and "cold" junctions of the thermopile. This coefficient is used to scale the output signal of the detector to indicate absolute gas concentration values.
However, thermopile detectors as a class are known to suffer from thermal drift. Thermal drift causes a slow variation in the D.C. voltage output of the detectors and leads to measurement inaccuracies. Raemer overcomes this problem by using A.C. coupling between the thermopiles. This solution to the thermal drift problem, however, renders the
resulting analyzer incapable of measuring absolute, steady state, gas concentrations. Instead, such devices are merely able to respond to changes in the concentrations of the gases being measured. One method which attempts to obtain absolute gas concentrations involves modulating, or "chopping", the incident energy beam. This technique is taught in U. S. Patent No. 4,423,739 to Passaro. In order to accomplish the required modulation, mechanical means such as a chopper are employed by Passaro. The resulting analyzer is large, and is thus subject to failure as the moving parts of the beam modulation apparatus wear. Such failure may have catastrophic consequences rendering the analyzer inoperable.
An alternative method of modulating an energy beam is to construct a modulated source. Devices of this type are disclosed in U.S. Patent No. 3,745,349 to Liston, for example. However, to date, analyzers employing modulated sources have suffered from slow response times due to the relatively long period required for the infrared emissions to decay after the excitation power is removed. Another shortcoming of modulated sources, familiar to those skilled in the art, is the relative instability of their energy output.
In addition, the imbalance between the excitation signals of two detectors has been utilized as a detection mechanism in the field of infrared motion detectors. In these devices, two detectors are arranged and their signals combined so that a null signal results when the detectors are identically excited. Any change in the environment that causes the output signal of one detector to change by an amount different than that of the other detector results in a non-zero differential signal. Thus, by subtracting the two output signals from two identical detectors an imbalance in the input excitation of the two detectors can be measured. An application of the method of subtracting two detector output signals, perhaps in an attempt to stabilize thermal drift, involves totally blocking one of the thermopile detectors with metal foil. This approach, however, is
undesirable because of results in uneven heating of the substrate, thus causing an error signal to be introduced into the output.
Thus, at the present time, the devices available provide absolute concentrations with the trade-off of either decreased reliability in the case of mechanically shuttered devices, or slow response times in the case of modulated source devices. Eliminating either of these drawbacks using conventional thermopile detectors results in a device prone to inaccuracies due to thermal drift, or one incapable of determining absolute concentrations if the thermal drift is stabilized using A.C. techniques. Therefore, a reliable gas analyzer using thermopiles that is both immune to thermal drift and capable of providing absolute gas concentrations would be highly desirable.
SUMMARY OF THE INVENTION It is an object of the present invention to provide methods and apparatus for detecting the absolute concentration of the constituents of a gas stream using thermopile detectors. Another object of the present invention is to eliminate the thermal drift from which thermopile detectors are known to suffer, while retaining the capability of measuring absolute gas concentrations. Yet another object of the present invention is to develop a method of selecting overlapping optical filters which takes into account that the filters are affected by temperature. Preferably, the invention further provides a computation technique which allows such filters to be used to measure the concentration of a plurality of constituents in the measured gas. It is a further object of the present invention to accomplish the above objects without modulating or "chopping" the incident energy beam.
Accordingly, the present invention provides a shutterless capnograph comprising an apparatus which has no moving parts, which does not require a modulated source of infrared radiation, and which does not suffer from thermal
drift, thereby providing a device which overcomes the limitations of the prior art. The improvements provided by this invention are arrived at through a novel optical stabilization technique applied to thermopile detectors. A preferred embodiment of the invention provides a highly accurate, stable, reliable, light-weight and economical infrared capnograph which will permit increased clinical and diagnostic monitoring of the absolute concentrations of carbon dioxide, nitrous oxide or other gases present in the expired airstream of a patient.
The present invention provides apparatus comprising one or more channels of a shutterless infrared gas analyzer for detecting the concentration(s) of at least one gaseous component of a substantially gaseous flow stream. The apparatus comprises at least two infrared detectors which produce electrical signals when illuminated by optical energy and means for illuminating the detectors with infrared radiation. Means for directing at least a portion of the gas being analyzed between the detectors and the infrared source is provided along with means for attenuating the infrared energy illuminating at least one of the detectors and circuitry for combining the signals produced by the detectors to produce an optically stabilized signal. If desired, the apparatus may further comprise additional detectors, attenuators, and circuitry, sufficient to comprise separate channels of a multi-channel device. The apparatus may also function in conjunction with a host processor for storing, processing, and displaying the measured information.
As mentioned above, the drift usually associated with thermopile detectors is produced by a variation in the reference junction temperature. In normal construction, a thermopile detector is comprised of a series of thermocouple junction pairs. One element of each pair is attached to a substrate (typically ceramic) and called the reference junction while the other element is unattached, suspended over an opening, and is called the sensing element. Incident radiation illuminates the suspended sensing element and raises
its temperature. A voltage develops that is related to the difference in temperature between the suspended and substrate attached elements. Drift arises, however, because the incident radiation and room temperature variations change the temperature of the reference junctions. The temperature differential between the sensing (hot) and reference (cold) junctions is thus altered; this condition in turn causes the output voltage to change, or drift.
A preferred embodiment of the present invention incorporates a substantially identical pair of thermopile detectors mounted on the same ceramic substrate and connected in series opposition. In this configuration, balanced and equal incident radiation illuminating the pair will produce no signal. Because the reference junctions of both detectors are on the same ceramic substrate and at substantially the same temperature, a drift in substrate temperature will produce no discernible change in output signal, thereby overcoming the above-mentioned problems. In order to make the system respond to incident radiation, an optical filter, or attenuator, with a transmission coefficient of approximately 0.50 (i.e., 50% transmission) is placed over one of the thermopile elements in the pair. With the filter in place the system responds to incident radiation but is substantially insensitive to other thermal changes. Also, in the present invention the output signals of the two thermopile detectors are subtracted. This is done not for the purpose of measuring any imbalance in the signals exciting the detectors, but rather to eliminate the effect of a variation in the background signals. For example, variations due to thermal drift are cancelled since they are common to both detectors. The problem of uneven heating described above in connection with the method of blocking one of the detectors is eliminated in the present invention by utilizing attenuation filters that do not cause uneven heating of the substrate. These filters have known but different transmission coefficients so that the incident excitation
signal can be accurately determined from the differential signal, as set forth in greater detail below.
In order to permit the analysis of specific gases, optical bandpass filters which permit the transmission of wavelengths within a specified bandwidth (i.e., that absorbed by the gas of interest) are disposed between the detectors and the gas directing means, thereby filtering the radiation prior to its striking the detectors. Since the center wavelength and bandwidth of optical filters vary with temperature, thermal drift should be considered when designing the optical bandpass filters. A preferred method for selecting the bandpass filters is described herein.
In accordance with a preferred embodiment of the invention, the method of selecting an optical filter for use in a gas analyzer which measures the concentration of a predetermined gas constituent in accordance with the invention comprises the steps of:
(a) selecting a nominal center wavelength for the optical filter which coincides with a center wavelength readily absorbed by the constituent;
(b) selecting a nominal half-power bandwidth just wide enough to pass substantially all of the wavelengths absorbed by the constituent;
(c) varying the nominal center wavelength so as to minimize modulation change as the optical filter is shifted over a predefined range of wavelengths;
(d) varying the nominal center wavelength and its bandwidth so as to minimize cross-talk between the constituent and other constituents; and (e) varying the nominal center wavelength and its bandwidth so as to minimize changes in cross-talk as the optical filter is shifted over the predefined range of wavelengths.
In accordance with the method of the invention, the half-power bandwidth is preferably selected to have a bandwidth greater than approximately 0.lμm for manufacturability, while the modulation change is minimized
over a 0.3μm filter shift, a typical shift due to temperature change during normal operation. Also, the predetermined respiratory gas constituent is preferably C02 where the cross¬ talk modulation for C02 by N20 is less than approximately 0.7%. On the other hand, the predetermined respiratory gas constituent may be N20 where the cross-talk modulation for N20 by C02 is less than approximately 1.0%. Also, steps (c) - (e) are preferably performed by computer simulation.
Circuitry for processing the signals produced by the detectors converts these signals into a ratio which is related to the absolute concentration of the gas constituent being measured. As a result, the capnograph of the present invention is substantially unaffected by changes in the ambient temperature of the environment in which it operates. The technique described above for combining and processing the detector signals is referred to as "optical stabilization" herein.
A capnograph made in accordance with the present invention may further comprise additional pairs of thermopile detectors, attenuators, and optical bandpass filters, and such additional circuit means so as to form a multichannel device for determining the concentrations of more than one gas constituent. In such a multichannel embodiment, each of the channels has an optical filter permitting the transmission of radiation within the absorption bandwidth of the constituent being monitored.
In a preferred embodiment of the invention, the constituent concentrations of multiple gas constituents are computed in accordance with the "optical stabilization" techniques of the invention by selecting the optical filter in accordance with steps (a) - (e) for measuring the concentration of a first constituent such as C02 and repeating steps (a) - (e) for selecting an additional optical filter for measuring the concentration of a second constituent such as N20. Then, the real-time concentration of C02 and N20 may be calculated in accordance with the invention by correcting the concentration of C02 and the concentration of N20 with respect
to a reference. Preferably, such a correcting step comprises the step of separating the measured concentration of N20 and
C02 using a mathematical polynomial fitting technique or a least squares fitting technique. Thus, the gas concentration computing step in accordance with the invention may calculate
C02 and N20 concentrations in accordance with the following equations:
C02 = Ac + (Be * LC) + (Cc * LC2) + (Dc * LC3) +
(Ec * LN) + (Fc * LN ) + (Gc * LIT) + (He * CP) , and N20 = An + (Bn * LC) + (Cn * LC2) + (Dn * LC3) +
(En * LN) + (Fn * LN ) + (Gn * T) + (Hn * CP) ; where Ac, Be, Cc, Dc, Ec, Fc, Gc and He and An, Bn, Cn, Dn,
En, Fn, Gn and Hn are constants and where LC is the natural log of said generated C02 signal divided by the reference signal, LN is the natural log of said generated N20 signal divided by the reference signal and CP is the product of LC and LN.
In particular, a method of determining the real-time gas concentrations of C02 and N20 using a respiratory gas analyzer in accordance with the invention preferably comprises the steps of: detecting the concentration of C02 and N20 in the expired air of a patient by passing infrared energy through expired air and through respective optical filters which respectively pass frequencies of infrared energy readily absorbed by C02 and N20; generating signals representative of the detected concentration of C02 and N20 in the expired air of the patient; simultaneously ratioing the detected C02 and N20 signals with respect to a reference signal; and computing the real-time gas concentrations of C02 and N20 in the expired air of the patient using varied weightings and C02 and N20 product terms so as to compensate for filter drift, gas cross-talk and collision broadening of the respective optical filters due to temperature variations over time.
A capnograph system made in accordance with the present invention may also have a host system for collecting and processing data collected by the detectors. Such a capnograph systemmay also include a device for indicating the absolute concentration of the constituent or constituents beingmonitored, therebyprovidingthe necessary concentration information in a usable form.
A preferred infrared energy detector for use in a capnograph made in accordance with the present invention is also disclosed. The detector uses paired thermopiles preceded by an analytical or reference bandpass filter, with a neutral density filter also placed in the optical path of one detector in the pair. Thus, two channels — one analytical and one reference — are preferably provided to permit measurements to be taken while the background effects are cancelled.
Mostpreferably a capnograph made in accordance with the present invention will be used to determine the absolute concentrations of carbon dioxide and/or nitrous oxide in the inspired and expired airstream of a patient. However, the optically stabilized detectors disclosed need not be so limited, and may be utilized to determine the concentrations of any of a number of constituents, as will be apparent to one skilled in the art from the following description.
BRIEF DESCRIPTION OF THE DRAWINGS FIGURE 1 is a foreshortened isometric view of the thermopile detector assembly of the present invention;
FIGURE 2 is a partially schematic plan view of a gas analyzer and related circuitry, constructed in accordance with the present invention; FIGURE 3 is a schematic of further circuitry, illustrating the use of the present invention in conjunction with a host computer;
FIGURE 4 is a schematic illustration of a detector and its optical path, constructed in accordance with the principles of the present invention;
FIGURE 5 is a schematic representation of the processing performed on the signal detected by the apparatus of FIG. 4;
FIGURE 6 is a perspective view of a detector constructed in accordance with the schematic illustration of FIG. 4;
FIGURE 7 is an exploded view of the detector of FIG. 6; and
FIGURE 8 is a schematic representation of an alternate embodiment of the detector of the present invention.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS A. Optically Stabilized Thermopile Detector
FIGURE 1 depicts a thermopile detector construction in accordance with the invention. It can be seen that such a thermopile detector is actually a group of thermocouples connected in series. As shown, the present invention utilizes a first set of thermocouples, comprising a first thermopile 12, and a second set of thermocouples, comprising a second thermopile 14, both preferably mounted on a common ceramic substrate 10, as shown. The hot, or sensor, junctions of the thermocouples are typically denoted H, while the cold, or reference, junctions are denoted C. Most preferably, the thermocouples 12 and 14 are the type comprising a metallic circuit deposited upon a polyester film, such as Mylar, or other suitable substrate. A preferred embodiment utilizes thermocouples having a substrate thickness of 1 mil.
Also shown in FIGURE 1 is the placement of neutral density (attenuation) filters 16 and 18, each of which has a different transmission coefficient. In a preferred embodiment, filter 16 has a transmission coefficient of 1.00 and filter 18 a transmission coefficient of 0.50.
In a preferred embodiment of the invention, detector 12 is electrically connected to detector 14 in series opposition. The resulting difference signal will only be indicative of the incident radiation. Signals caused by substrate temperature changes will cancel, provided the
temperature is uniform across the substrate. This aspect of the invention will be discussed more fully with respect to the mathematical derivation provided below.
Thus, in a preferred embodiment, the thermopile detectors 12 and 14 are "optically stabilized". This design overcomes the requirement for source modulation described above by rendering the detector substantially immune to thermal drift. Since each of the thermopile detectors 12 and 14 is exposed to the same incident energy beam but has different attenuation filters 16 and 18 in line with that energy, the differential signal from the two detectors is therefore related only to the incident energy and the transmission coefficients of the neutral density attenuation filters. The stabilization technique disclosed by the present invention is useful with a variety of detectors, such as the Model DR34 manufactured by Dexter Research Inc., Ann Arbor, Michigan. Moreover, since the optically stabilized detectors of the present invention are sensitive to a wide range of wavelengths, and are relatively inexpensive and rugged, they can be beneficially used in a shutterless infrared capnograph for monitoring respiratory gas.
In FIGURE 4, the optical path of an embodiment of a detector 200 constructed in accordance with the principles of the present invention is shown. Infrared energy is emitted from an IR source 300 and passed through the gas being analyzed, which in this example is C02. Some of the infrared energy is absorbed by the gas in front of the IR source 300, while the remainder of the infrared energy impinges upon a group of filters, including a neutral density filter 210, an analytical C02 filter 212, and a reference filter 214. The neutral density filter 210 uniformly attenuates all wavelengths of energy which are incident upon the filter. The analytical C02 filter 212 is a bandpass filter which passes a narrow band of energy, including those wavelengths absorbed by C02. The reference filter 214 is also a bandpass filter which passes a band of energy excluding those wavelengths absorbed by C02.
As shown, the neutral density filter 210 overlaps or shadows the subsequent analytical and reference filters 212 and 214, as well as two of the thermopile detectors 216 and 218, C02 detector B 220 and reference detector B 222. Thus, each detector A 216 and 218 will receive energy which has been transmitted by its respective analytical filter 212 or reference filter 214, but not the neutral density filter 210. Each detector B 220 and 222, on the other hand, receives energy passed by its respective analytical filter 212 or reference filter 214 as well as the neutral density filter 210. Following the filters is an aperture or window 230 which transmits energy to only the detector areas of the substrate on which the detectors are formed. The outputs of the detectors are passed to a processor 240 for the production of an output signal 242 representative of the concentration of C02 in the gas being analyzed.
A main purpose of the neutral density filter 210 is to provide signals which allow the desired signal to be distinguished from signal interference such as background (local thermal) effects. A number of attempts have been made in the prior art to deal with these effects. For instance, U.S. Patent No. 4,772,790 to Aldridge attempts to isolate these effects by using ten thermopiles for each detector channel, with only a central portion of the area of the thermopiles being used to develop a usable signal. In addition, U.S. Pat. 3,539,804 to Billetdeaux et al. attempts to deal with the problem through heating and shielding of the detector, and also includes a rotatable neutral attenuator as a zeroing mechanism in assembly balancing of a C02 detector. By contrast, in an embodiment of the present invention, local thermal effects are eliminated by developing signals resulting in the cancellation of these effects when the signals are combined.
For instance, an examination of the C02 channel reveals that the C02 detectors A 216 and B 220 are in close proximity to each other on the same ceramic substrate 10 and hence experience similar effects from background and other
local radiation sources. Each detector produces an output signal containing a component I, the desired incident infrared signal, and a component b due to background effects. Absent the presence of the neutral density filter 210, the C02 detectors 216 and 220 would produce signals of the form: Det. A = R(Ia + b) , and Det. B = R(Ib + b) , where R is the response of the detectors.
It is desired to maximize the desired signal component I, which may be done by selecting a neutral density filter 210 which is about 100% transmissive at the absorption wavelength for C02. In this case A would equal B. The signal component I would be maximized but the difference of the two detector output signals would yield zero. Correspondingly, if the neutral density filter were chosen to be a 100% blocking filter (i.e., 0% transmissive, I=0) at the absorption wavelength, their difference would be A - B = RIa. This would have the undesirable effect of shadowing detector B 220 from all radiation from the IR source 300, thereby blocking detector B 220 from reception of incident energy, and would further set up an undesired thermal differential between the two detectors. Moreover, neither of these two choices leads to the ability to adequately distinguish and eliminate background effects. If, however, the neutral density filter 230 is chosen to have a transmissivity of about 50%, for instance, the desired incident energy component I would be differentiated and the background effects would substantially cancel. A neutral density filter of this character would produce signals of the form:
Det. A = R(I + b) , and
Det. B = R(0.5I + b) .
Thus, when the difference of A minus B is taken, the result is 0.5IR with the background effects b cancelling. Other transmissivity characteristics may be chosen in consideration of incident energy levels and the thermal effects of shadowing
by the neutral density filter 210. For example, a 0% transmissive filter may be used.
The details of the processor 240 of FIGURE 4 are shown in FIGURE 5. As shown, the output signals of the A and B detectors 216 and 220 for the C02 and for the reference channels 218 and 222 are respectively subtracted at subtractors 242 to eliminate background radiation and substrate temperature effects, as explained immediately above. In a preferred embodiment of the present invention, this subtraction step 242 is accomplished by connecting the thermopile detectors in series opposition, so that the subtraction of their outputs is inherent in their interconnection. A ratio 244 is then taken of the incident signals after this cancellation of thermal effects. The ratio 244 of the analytical and reference signals eliminates proportionate effects of spurious absorption, such as those resulting from the accumulation of undesired particulate matter on the components of the detectors. For example, for an incident radiation signal having a normalized value of 0.8 C02 and filters which are contaminated with substances which absorb some of the radiation being measured, the transmissivity of the optical paths is reduced to 0.9 of full transmissivity. The output signal would thus be the product of these losses, or 0.72 C02. When the reference channel is subject to the same contamination, which may be substantially achieved by close proximity of the detectors and filters, the same factor of 0.9 would be present in signals of the reference channel. Hence, the ratio of the C02 and reference signals would contain the 0.9 factor in both the numerator and the denominator for a ratio of one, eliminating the effects of contamination.
Computations are made on the calculated ratio to further refine the measurement of C02. The measurement of absolute substrate temperature, provided by the thermistor 246 in the detector, and the appropriate Seebeck coefficient 248 are factors taken into consideration in this computation. The temperature measurement is incorporated into the present
invention as shown in the functional block diagram of FIGURE 5. As will be understood by those of ordinary skill in the art, the resultant measurement of C02 is obtained by solving the well-known Bier's Law equation which is of the form I = I0e"xαc, where I0 is a constant coefficient, x is path length, is the absorption coefficient for C02, and c is concentration. The reference channel yields a value for I, since the absorption coefficient and concentration of C02 in the reference channel are zero by reason of the elimination of C02 wavelengths from the reference filter passband. The detector output signal supplies the value for I, and the equation is solved for c, the concentration of C02 in the gas being measured. The above computations are preferably performed by a host processing system of the type described below with respect to FIGURE 3.
In a preferred embodiment, a second analytical channel is employed for N20. This channel is identical to the
C02 channel in FIGURE 4, except that the analytical filter 212 is chosen to pass wavelengths representative of the absorption of infrared energy by N20 instead of C02. The outputs of the
A and B detectors for the N20 channel are combined to eliminate background and substrate thermal effects and the reference channel outputs are combined to yield a measurement of N20 concentration, as in the C02 computation set forth immediately above. In addition, such N20 concentrations may be used in accordance with the techniques taught by
Severinghaus, M.D., Larson, M.D. and Eger, M.D. in an article entitled "Correction Factors for Infrared Carbon Dioxide
Pressure Broadening by Nitrogen, Nitrous oxide and Cyclopropane" in Anesthesioloαy, May-June, 1961, pp. 429-432, to correct for collision broadening errors induced in the C02 gas concentration computation by nitrous oxide.
FIGURES 6 and 7 respectively illustrate in perspective assembled and exploded views a detector constructed in accordance with the present invention. The detector 200 employs the same principles and optical paths shown in FIGURE 4. As shown in the perspective of the
assembled detector 200 shown in FIGURE 7, the neutral density filter 210 is disposed above the reference filter 214 and the analytical filters 212 and 213. In the embodiment of the present invention shown, analytical filters 212 and 213 selected to analyze C02 and N20 are provided, as well as the reference filter 214. Each of these three filters is "shadowed" or blocked by the neutral density filter 210 as explained above with reference to FIGURE 4.
Further details of the construction of the detector are shown in the exploded view of the infrared detector 200 in FIGURE 7. Beneath the filters 210, 212, 213, 214 is the filter aperture 230. The thin film thermopile detectors 216, 218, 220 and 222 are sandwiched between two ceramic spacers 215 and 225. In FIGURE 7, the N20 detectors 224 and 226 are not visible. A foil background 245 is provided beneath the ceramic spacer/thermopile sandwich to block the further passage of the incident radiation, but it allows heat to be conducted through the device. Finally, beneath the foil background 245, a thermistor 246 is disposed which, as explained above, provides an absolute temperature indication. The above-mentioned components are mounted on a TO-8 12-pin header 250, as shown.
While the illustrated embodiments are shown to employ thermopile detectors, the principles of the present invention may also be applied through the use of other infrared detecting devices, such as thermistors, thermocouples, pyroelectric detectors, Golay cells, and PbSe photodetectors, for instance. However, the arrangement shown has the benefit that it is easy to manufacture, is small in size and light in weight.
B. Mathematical Derivation
The signal or voltage output, S, of a thermopile detector can be expressed as:
S = (T * I * R) + B, where T is the filter transmission coefficient, assuming the filter is between the detector and the energy source, I represents the incident energy of the source, R is the responsivity of the detector, and B is the component of the signal attributed to background "noise", including room temperature and detector package temperature variations. In a preferred embodiment of the present invention, two detectors 12 and 14 are mounted on a common substrate 10 as shown in FIGURE 1, creating a condition wherein the value of B for each detector is equal. The detectors 12 and 14 are used to simultaneously monitor an incident energy beam. Filters 16 and 18 having known but different transmission coefficients are placed between each detector 12 and 14 and the source of incident energy.
The two equations describing the outputs SI and S2 of the two detectors 12 and 14 are combined and solved for I, yielding the following equation, which is independent of the value of B:
Si -
I = (T, * R,) - (T2 * R2)
In a most preferred embodiment, the numerator of the above equation, S1 - S2, which represents the difference signal, is computed by connecting the detectors 12 and 14 in series opposition. It should be apparent to one skilled in the art, however, that the difference signal could be calculated using a digital computer or other electronic circuit. The equation demonstrates that the stabilization technique of the present invention eliminates the background effects term, B, from the final output of the circuit, thereby rendering a device substantially insensitive to room temperature effects.
C. Source Drift
As will be understood by those having ordinary skill in the art, variations, or drift, in the infrared source may also cause drift in the output signals of the thermopile detectors 12, 14, 22, 24, 32 and 34 illustrated in FIGURE 2. This drift may be eliminated by using two analytical filters, a first having its transmission at the absorption wavelength of the gas being analyzed, and a second designed to transmit some other unabsorbed wavelength. Energy in the unabsorbed band serves as a reference. This technique of referencing is well-known and is described in textbooks, see, e.g., Mullard Limited, "Applications of Infrared Detectors," ISBN # 901232 22X, Chapter 10, the pertinent portions of which are incorporated by reference as if fully set forth herein.
D. Optically Stabilized Gas Analyzer A typical construction of a preferred embodiment of the gas analyzer of the present invention is illustrated in FIGURE 2. As shown, the circuitry and other components which comprise the present invention will generally be assembled into a housing 50, appropriately configured to enable the various components to function as described below. Those of ordinary skill will be aware of the design parameters associated with the packaging of devices made in accordance with the present invention. Some of the considerations relevant to the design of the housing 50 are: external interface design, cooling requirements, shielding and isolation characteristics, size and weight, integration with a breathing tube, sample cell or other means of delivering a sample of the patient's expired air, and aesthetic and utility considerations. An input signal 80 from a host processor or other source is transmitted through an analog hybrid 42, the signal from which excites a source of infrared radiation 40. A reflector 52 is provided to efficiently direct an incident beam of infrared radiation. The incident radiation passes through an optical window 54 and into the region of the gas stream to be analyzed 58. The gas stream itself may be
contained in apparatus such as a clear tube or tube-like structure (not shown) , which directs a portion of the gaseous flow through the measurement region 58. After the incident infrared radiation has passed through the gas stream, it will pass through another optical window 56, bandpass filters 60, 61 and 62 and neutral density attenuation filters 16, 18, 26, 28, 36 and 38 before finally reaching the thermopile detectors 12, 14, 22, 24, 32 and 34.
The embodiment pictured in FIGURE 2 is designed to measure the concentration of two gases and may be denominated a three-channel detector, with one channel serving as a reference channel. As will be appreciated by those of ordinary skill in the art, detectors may be constructed having anywhere from one to many channels, depending upon the number of constituent gases to be detected and analyzed. Thus, for example, in FIGURE 2 in addition to detectors 12 and 14, second and third detector assemblies 22 and 24 and 32 and 34 are used to create three channels. Positioned above each detector assembly is a bandpass filter 60, 61 and 62, each channel employing a different bandpass filter. Thus, it is understood that a first channel comprises detectors 12 and 14, attenuation filters 16 and 18, and bandpass filter 60. A second channel comprises detectors 22 and 24, attenuation filters 26 and 28 and bandpass filter 61, and so forth. The output signal for each channel, i.e., the difference signal described above, is passed through its respective amplifier 20 and transmitted to a multiplexer 44, which also receives a signal from a temperature sensor 41. The multiplexer 44 then outputs a detector signal 45, which is also combined with a signal from an electronically erasable programmable read only memory (EEPROM) device 46 to fashion an output signal 90. The Seebeck coefficient, described above, is stored in the EEPROM 46, and is used for scaling the detector signal 45 to indicate absolute gas concentration values.
Attenuation filters 16, 18, 26, 28, 36 and 38 are mounted directly on the ceramic substrate. Therefore, the
output of the temperature sensor 41 will reflect the temperature of the attenuation filters as well as the reference junction temperature of the thermodetectors. A heater 30 is also shown which may further regulate the temperature of the apparatus, if desired.
E. Bandpass Filters
The bandpass filters 60, 61 and 62 are placed an equal distance in front of their associated detectors, i.e., they are placed in a common plane, so that all three thermopile detectors have an equal field of view. Each filter is characterized by its center wavelength and half-power bandwidth. As described by Conlon et al. in U.S. Patent No. 4,914,719, for example, the conventional approach to specifying the bandpass filters would be to center the filters in the absorption band of interest and to make the bandwidth as wide as the absorption band. The present invention takes a different approach than that of Conlon et al., however, and utilizes bandpass filters having center wavelengths shifted away from the center of the absorption band and bandwidths wider than the absorption band. The method of specifying the bandpass filters in accordance with the present invention is described in more detail in the following section. The shifting and broadening of the bandpass filter in accordance with the invention minimizes cross-talk between channels and minimizes the effects of temperature variations and manufacturing tolerances on the filters.
For example, in a preferred embodiment of the invention one channel is used to detect C02, another channel is used to detect N20, and a third channel is used as a reference channel. The absorption band of C02 is centered at 4.256 microns, and the half-power wavelength is 0.180 microns. In a preferred embodiment of the present invention, however, the bandpass filter used for the C02 channel is centered at 4.2624 microns with a half-power bandwidth of 0.2088 microns. Due to manufacturing tolerances, the half-power wavelengths may vary by plus or minus 0.03 microns at the short wavelength
edge, and plus 0.03 microns to minus 0.05 microns at the long wavelength edge.
The absorption band for N20 is centered at 4.5400 microns, with, a bandwidth of 0.3 microns. For the reasons stated above, in a preferred embodiment the bandpass filter for the N20 channel is chosen to have a center wavelength of 4.5397 microns and a bandwidth of 0.3190 microns. The short wavelength half-power point may by plus 0.01 microns to minus 0.03 microns. The long wavelength half-power point varies by plus or minus 0.03 microns.
The reference filter has a center wavelength of 4.8265 microns and a half-power bandwidth of 0.2088 microns. The short and long bandwidth edges may vary by plus or minus 0.03 microns. The bandwidth is designed to be as close as possible to the bandwidth of the C02 filter. This minimizes the thermal gradients across the substrate since the total energy passing through this filter will approximate the total energy passed by the C02 filter. The center wavelength is selected to be as close as possible to the center of the C02 and N20 bands, while at the same time maintaining isolation between those channels and the reference channel. The technique of the invention thus takes a filter responsive primarily to C02 but some N20 and a filter responsive primarily to N20 but some C02 and convolves their outputs. As will be described in more detail in Section F below, this is accomplished using varied weightings which allow the C02, N20 and reference signals to be processed simultaneously using, for example, a polynomial fitting algorithm or least squares fitting technique to generate the weighing factors. A calibration regression may then be run to print the weighting factors, which are fit to the polynomial to get the desired coefficients. Thus, this novel use of band shifting and widening, to minimize cross-talk and the effects of thermal shift and manufacturing tolerances, in the design of the bandpass filters 60, 61 and 62 is integral to carrying out the preferred embodiment.
1. Analytical Bandpass Filter Specification
This section describes a method for specifying bandpass filters for a preferred embodiment of an optically stabilized capnograph in accordance with the present invention. The following example describes the specification of bandpass filters employed by the C02, N02 and reference channels of a three channel analyzer. The method described does not deal with the construction of filters, but rather with a novel method of specifying bandpass filters for use in infrared gas analyzers, or capnographs. a. Conventional Technique
From the prior art, as represented by U.S. Patent No. 4,914,719 to Conlon et al., for example, a technique for specification of infrared bandpass filters can be summarized. The guidelines, or steps, can be stated as: 1. Specify a filter that has its center wavelength coinciding with the center wavelength of the gas intended to be analyzed. This is done to yield maximum specificity and minimum cross-talk.
2. Specify the half power bandwidth just wide enough to include all of the wavelengths absorbed by the gas of interest. This optimizes modulation or the amount of energy passed by the filter and absorbed by the gas.
In practice it is understood that filter center wavelengths change with temperature. This is shown at page 41 of the OCLI Infrared Handbook (1970) Optical Coating Laboratory Inc. In the construction of gas analyzers incorporating filter designs, engineers have overcome this filter shift by operating the filters in a temperature controlled environment. Such a design is taught by U.S. patent No. 4,423,739 to Passaro. This approach, however, is expensive and complicated. It is also recognized that manufacturing tolerances on infrared bandpass filters can result in a range of center wavelengths from a single production run. This is expressed as a +/- tolerance of typically 1/2% to 1% on the center wavelength. Accommodating this variation is often accomplished by calibrating each analyzer individually. This practice is also expensive.
In specifying C02 and N02 filters in accordance with the present invention, the first two steps consist of steps
1 and 2 of the conventional technique; i.e., conventional filters are specified. In the present example, gas spectra from SADTLER INDEX were used. The resulting conventional filters are characterized by the following filter specifications, titled VER 1.2. They follow in tables which show the filters as specified and as shifted by +.03μm, which is approximately 0.7%. It corresponds to the shift that would be produced by a rise in temperature of about 50° C. It also corresponds to the manufacturing tolerance specified. For reference, the spectra of C02 and N02 are also included on the graphs. b. Analysis and Computer Performance Simulation
Following the drafting of the VER 1.2 specifications, a computer simulation technique was developed by the present inventors to test the validity of the theoretical filter specification by mathematically predicting the modulation of light as if it were passed through the gas and subsequently through the specified filters. A Perkin Elmer 1640 FTIR spectrophotometer (PE 1640) was used to measure the absorption spectra of the gases of interest, C02 and NzO. The gases were at atmospheric pressure and in a 0.5 inch path length cell. 5.0% C02 and 50% N20 balance 02 were used. The spectra data were transferred into a personal computer, with data reported every O.lμm. These spectra were combined in a spreadsheet program with mathematical models of the filters to allow easy multiplication and division of the recorded values. Modulation percentages were calculated by multiplying the gas spectra by the filter models and dividing the result by the product of the spectra and a 100% transmission model. Modulation, cross-talk, and effects of filter shifts were studied. The results of this analysis will be discussed below in the section on performance comparison. The analysis showed that in order to meet the design goals of the project either individual calibration or temperature control of the filters,
or both, would have to be specified. It was decided to optimize the filter specifications to avoid these complications by mathematically altering the filter under investigation while preserving a realistic shape of the bandpass.
c. Optimal Filter Specification In accordance with the invention, the conventional filter specifications were then optimized by the following three steps: 3. The filter center wavelength was varied so as to minimize the modulation change as the filter shifts over 0.3μm.
4. The filter center wavelength and bandwidth were varied to minimize cross-talk. 5. The filter center wavelength and bandwidth were varied to minimize the change in cross-talk that occurred due to the filter shifting 0.3μm.
To provide a realistic shape to the simulated bandpass filter, the transmission of a real infrared bandpass filter was measured with the PE 1640. A data file was created which preserved the cut-on and cutoff edges of the real filter but replaced the middle section with a flat transmission of 80%. A program changed the width of the simulated filter by altering the width of the flat top to allow testing of various filter bandwidths. The simulated filter was restricted to half-power bandwidths of greater than 0.lμm to insure that the resulting optimized filter was manufacturable. With the gas spectra in memory and the adjustable ideal filter file prepared, a personal computer was programmed to sweep through all possible filter combinations for C02 and N20 filters and to select the best combination as indicated by:
1. Maximum modulation;
2. Minimum modulation change over 0.3μm filter shift,
3. Cross-talk modulation less than 0.7% for C02 (equivalent to the effect of about 1mm Hg C02 in the analyzer sample cell) ; and 4. Cross-talk modulation less than 1.0% for N20 (equivalent to the effect of about 10mm Hg N20 in the analyzer sample cell) .
d. Performance and comparison of VER 1.2 and VER 2.1 Filters
(i) Modulation
The optimized filter VER 2.1 has slightly less modulation than the conventional version.
This is due to the fact that both the VER 2.1 C02 and N20 filters are slightly narrower than the conventional
VER 1.2 filters and not centered on the center of absorption for C02 and N20, respectively. This reduction is of little or no significance in the overall design of the analyzer.
(ii) Filter Shift Error The sacrifice in modulation was offset by a much improved performance with respect to errors when the filters shift.
The optimized filter specification C02 error due to a shift of 0.3μm is reduced from 1.7mm Hg to 0.4 and the N20 error went from 112 to 26mm Hg.
(iii) Cross-talk
Cross-talk increased slightly fromVER 1.2
mmHg
0.6
e. Summary The method yielded a filter set which is slightly different from that selected by conventional methods and, according to the simulations, better suited to yield a manufacturable gas analyzer with little or no individual
adjustment. Preferred filter specifications are set forth in Table I below:
TABLE I
INFRARED FILTER SPECIFICATIONS ALL FILTERS:
Dimensions: 0.187" X 0.156" ,+1- .003"
0.020", +1- .001" thick Edge Irregularities: <0.005" Temperature Shift: 0.01% per Deg. C. Angle Shift: 0.001% per Degree of angle. Transmission at Center: 80%
Cut-on - Cutoff Slopes: 3% Humidity Withstand per MIL-C-675A Abrasion Withstand per MIL-C-675A Coating Adherence per MIL-M-13508B
Spectral Characteristics below apply at 0 Deg. incidence and 30 Deg. C.
FILTER 1: Center 4.255 micron +/- 0.03 micron
Half Power Bandwidth 0.180 micron +/- 1% FILTER 2: Center 4.540 micron +/- 0.05 micron
Half Power bandwidth 0.300 micron +/- 1%
FILTER 3: Center 3.600 micron +/- 0.03 micron
Half Power Bandwidth 0.180 micron +/- 1%
F. Optically Stabilized Gas Analyzer Used With A Host Processor
FIGURE 3 illustrates a preferred embodiment of certain signal processing circuitry 100 used to carry out the present invention in conjunction with a host system. The input and output signals to the detector apparatus 80 and 90 are passed through an isolation and filtering subsystem 140 to insure patient safety. Output signal 90 of the detector flows through the isolation subsystem 140 to the AD converter detector 130, where it is digitized, before input into the microprocessor 120. Finally, the microprocessor 120 transmits the processed detector output — as a concentration value for the respective constituent gas — to the host system interface 110, and on to the host system (not shown) . The final output signal may be monitored, recorded, stored, or further manipulated. Alternatively, a signal from the host system to
the detector apparatus shown in FIGURE 2 may be transmitted to the host system interface 110, then to the microprocessor 120, through the isolation and filtering subsystem 140, and finally on to the input 80 of the detector apparatus. The signal from the host system is not converted from digital to analog form before being input to the detector apparatus.
In accordance with a preferred embodiment of the invention, the actual gas concentration is computed from the measured values using suitable software operating on the host system. Sample software modules are described for the preferred embodiment where the real time concentration of C02 and N20 in a sample cell is computed. As will be apparent to those skilled in the art, the following modules apply all the corrections required to provide an accurate measure of gas concentration for the case where the measured concentration would otherwise be affected by temperature drift of the components. Excluding Module 0 (INITIALIZATION) , the modules are intended to be run in sequence on the host processor every 10 msec (100 HZ) .
MODULE O: INITIALIZATION
At initialization time (power up) this module reads the calibration variables from the EEPROM 46 (FIG. 2) . Variables in the EEPROM 46 are preferably specifically computed for each sensor at factory calibration. This module also initializes the operating variables of the capnometer and includes the self calibration check of the hardware gain ratios for C02/REF and N20/REF. These ratios are preferably stored as values CERAT and NERAT, respectively.
MODULE l: SIGNAL ACQUISITION This module reads the analog voltages, then updates the averages and stores the values as follows:
Preferably, the sample rates are averaged such that channels 1 - 3 are sampled at 100 Hz and averaged for an averaged rate of 50 Hz (2 averages) . On the other hand, the remaining channels may be sampled at 20 Hz and averaged for an averaged rate of 2 Hz (10 averages) .
MODULE 2: COMPUTE TEMPERATURES
This module computes the source and detector temperatures.
Step One: Compute sensor resistances in accordance with the following equations:
Rth = (VTR/VDT - 1) * RTHSERIES and
Rpt = (VTR/VST - 1) * RPTSERIES, where RTHSERIES and RPTSERIES are calibration coefficients from the EEPROM 46 and are specific to each sensor. In a preferred embodiment, an approximation of these coefficients is RTHSERIES = Detector thermistor series resistor, or, e.g., 200,000 Ohms and RPTSERIES = Source RTD series resistor, or, e.g. , 2700 Ohms.
Step Two: Compute sensor temperatures in accordance with the following equations:
Dtmp = DtO + (Dtl * Rth) + (Dt2 * Rth2) + (Dt3 + Rth3) and Stmp = StO + (Stl * RPT) + (St2 * Rpt2) + (St3 * Rpt3) , where DtO, Dtl, Dt2, Dt3 and StO, Stl, St2 and St3 are coefficients from the EEPROM 46 and are specific to each sensor. In a preferred embodiment, an approximation of these coefficients is:
DtO = 81.65593
Dtl = -.0474753
Dt2 = 1.21559 E-5
Dt3 = 1.228346 E-9 StO = -205.4541
Stl = 1.921956
St2 = 2.432033 E-3
St3 = -1.183539 E-6
Step 3: Compute board temperature in accordance with the following equation: Btmp = Vbt * BTSF, where BTSF is the board temperature scaling factor and is the same for each processor board.
MODULE 3: COMPUTE BAROMETRIC PRESSURE This module computes the barometric pressure using a pressure transducer mounted on the signal processor board in accordance with the following equation:
Pb = Vpb * BAROFCTR, where Pb has units of mmHg and BAROFCTR is a constant and is the same for each transducer.
MODULE 4: COMPUTE DETECTOR TEMPERATURE GRADIENT
This module computes the rate of temperature change of the detector dDtmp/dT. This value is used in computation of the SPAN factors in module 5 and is calculated in accordance with the following equation:
DTgrad = Dtmp(-l) - Dtmp(0) / dT where DT grad is in units of Deg. C/min. ; Dtmp(-l) is the last detector temperature; Dtmp(0) is the current detector temperature and dT is the time between readings.
MODULE 5: COMPUTE TEMPERATURE SPAN FACTORS
The span factors are used to correct the ratios of
VC02/VREF and VN20/VREF over detector temperature variations.
This corrects for IR bandpass filter drift and source emission drifts. This correction stabilizes the zero reading of the instrument. These factors are used in the following module.
The C02 span factor is computed in accordance with the following equation:
C02TSF = DTSFO + (DTSF1 * Dtmp) + (DTSF2 * Dt p2) + (DTSF3 * Dtmp ) + (DGSF1 * DTgrad) + (DGSF2 * DTgrad ) + (DGSF3 * DTgrad ) where DTSFO, DTSF1, DTSF2, DTSF3 and DGSF1, DGSF2, DGSF3 are coefficients from the EEPROM 46 and are specific to each sensor. In a preferred embodiment, an approximation of these coefficients is: DTSFO = .930592
DTSF1 = .000904
DTSF2 = 2.6 E-8
DTSF3 = -4.8 E-12
DGSF1 = -.01603 DGSF2 = -.001258
DGSF3 = -4.7 E-5
The N20 span factor is computed in accordance with the following equation:
N2OTSF = NDTSF0 + (NDTSE1 * Dtmp) + (NDTSF2 * Dtmp2) + (NDTSF3 * Dtmp ) + (NDGSF1 * DTgrad) +
(NDGSF2 * DTgrad ) + (NDGSF3 * DTgrad ) where NDTSF0, NDTSF1, NDTSF2, NDTSF3 and NDGSF1, NDGSF2,
NDGSF3 are coefficients from the EEPROM 46 and are specific to each sensor. In a preferred embodiment, an approximation of these coefficients is:
NDTSF0 = 0.40924
NDTSF1 = 2.0482 E-2
NDTSF2 = -4.7731 E-4
NDTSF3 = 3.7300 E-6 NDGSF1 = -2.7813 E-4
NDGSF2 = 2.0171 E-3
NDGSF3 = -8.0728 E-6
MODULE 6: COMPUTE GAS CONCENTRATION
This module computes the C02 and N20 gas concentrations in mmHg in accordance with the following steps:
Step One: Adjust for sensor offsets in accordance with the following equations:
Cst = Vc - Coff,
Nst = Vn - Noff, Rst = Vr - Roff;
where Coff, Noff, Roff are coefficients from the EEPROM 46 and are specific to each sensor.
Step Two: Compute ratios and adjust using span factors and electronic hardware ratios in accordance with the following equations:
Rcr - (Cst/Rst) * C02TSF * CERAT, and Rnr = (Nst/Rst) * N20TSF * NERAT;
Step Three: Compute the log of ratios and cross product terms in accordance with the following equations: LC = Ln(Rcr) ,
LN = Ln(Rnr) , and CP = LC * LR; where LN is the natural log.
Step Four: Compute the gas concentration using both the C02 ratio, LC, and the N20 ratio, LN, as well as the product term, CP. The use of all the terms rather than just the C02 in computing concentration in accordance with the invention eliminates the effect of spectral overlap and collision broadening. The gas concentration in accordance with the invention is thus calculated as follows:
Cmmhg = Ac + (Be * LC) + (Cc * LC2) + (Dc * LC3) + (Ec * LN) + (Fc * LN ) + (Gc * LN ) + (He * CP) , and
Nmmhg = An + (Bn * LC) + (Cn * LC2) + (Dn * LC3) + (En * LN) + (Fn * LN ) + (Gn * LIT) +
(Hn * CP) ; where Ac ... He and An ... Hn are probe specific coefficients stored in the sensor EEPROM 46. In a preferred embodiment, an approximate value of these coefficients is:
Hn = -116.256
MODULE 7: CHECK FOR ERRORS
This module checks for some known error conditions and sets flags leading to actions. Such error conditions may include:
1. IR signals below 1.0 volt;
2. Source temperature below 350°C (Typical is 360°C) ;
3. Detector temperature above 60°C; 4. Gas Concentration Cmmhg less than -2;
5. Gas Concentration NmmHg less than -20;
6. Power supply out of range; and
7. Board temperature out of 0 - 40°C range. Further modules may be added to compute the derived variables from the real time gas concentrations described above. In addition, inspired and end tidal C02 and N20 concentrations and respiration rate may be computed. Also, the percentages from the mm Hg values and barometric pressure may also be computed in accordance with techniques well within the level of skill of those skilled in the art. Of course, a communication system may also be provided for allowing the capnograph to communicate with the host system.
Although certain preferred embodiments of the present invention have been described herein, it should be understood that the invention is not so limited. For example, FIG. 8 illustrates a schematic representation of an alternate portion of a detector in accordance with the invention. As explained above with reference to FIG. 4, the emission from a source of infrared radiation 300 falls upon the detector 202. However, in the embodiment of FIG. 8, the functions of the neutral density and analytical filters described above are combined into the same element. Specifically, a first filter 260 having a transmissivity of 100% and a second filter 262 having a transmissivity of 50% are provided. The remainder of the detector 202 is the same as that described above with reference to FIG. 4, i.e., an aperture or window 230, detectors 216 and 220 and processor 240 are again provided.
The reference portion of the circuit may either be constructed as shown in FIG. 4, or using a pair of 50%/100% filters, as shown in FIG. 8. Thus, the analytical filters may be functionally combined with the neutral density filter as one filter element.
Thus, many such variations of the specific embodiments heretofore described are within the scope of the present invention as defined by the following claims.