Abstract:
An apparatus for detecting gas concentrations includes a coded filter to oscillate proximate a resonant frequency. A photo detector is positioned below the coded filter such that the coded filter selectively blocks light that is directed at the photo detector. Optics are positioned to project spectral information on to the coded filter. A processor analyzes a signal received from the photo detector. The processor is adapted to weight a harmonic attic signal.

Description:
TECHNICAL FIELD 
       [0001]    The present invention relates to detecting natural gas concentrations and in particular detecting natural gas concentration using Dynamic Eigen Spectroscopy. 
       BACKGROUND 
       [0002]    Monitoring the atmospheric gasses in and around an oil refinery may be performed to determine whether hydrogen sulfide is present, and to qualify its concentration if it is detected. Monitoring a gas cloud heading toward an Army base may be useful to determine whether that cloud contains chemical warfare agents. Existing techniques, including TDLS (Tunable Diode Laser Spectroscopy), NDIR (Non-Dispersive InfraRed analysis), Polychromatry, and FTIR (Fourier Transform InfraRed analysis) all have limitations that, depending on the application, can limit their ability to detect atmospheric gasses at the desired level. 
         [0003]    Identifying unknown chemical contaminant in the atmosphere from afar can be very difficult because normal components of the atmosphere such as H2O and CO2 (water vapor and carbon dioxide) have spectral signatures that are similar to, and overlap with the spectral signature of many of the contaminants of interest. 
         [0004]    Existing Dynamic Eigen Spectroscopy can be designed to have superb chemical specificity, i.e., they can be designed to be very sensitive to the emission or absorption of one specific chemical or group of chemicals while being very insensitive to other chemicals or groups of chemical. Unfortunately, when such systems are manufactured they might not perform as well as modeled. 
       SUMMARY 
       [0005]    An apparatus for detecting gas concentrations includes a coded filter to oscillate proximate a resonant frequency. A photo detector is positioned below the coded filter such that the coded filter selectively blocks light that is directed at the photo detector. Optics are positioned to project spectral information on to the coded filter. A processor analyzes a signal received from the photo detector. The processor is adapted to weight a harmonic of the signal. 
         [0006]    In a further embodiment, a device for detecting gas concentrations includes a movable coded filter having multiple slits in a proof mass. An optical element is positioned to receive gas filtered light and spectrally separate the gas filtered light onto the coded filter wherein spectral bands run in the same direction as the slits, the slits positioned to cancel AC signals corresponding to at least one gas not of interest. A photo detector is positioned to receive the spectrally separated light through the oscillating slits in the moveable coded filter to provide an AC signal representative of a selected gas. A controller is coupled to receive the AC signal, convert the AC signal to a digital signal, to weight a harmonic of the received AC signal, and to correlate an amplitude of the AC signal with the weighted harmonic to a concentration of the selected gas. 
         [0007]    A method fir detecting a gas includes receiving light from a light source through a plume of gas, spectrally separating the light, oscillating a coded filter to selectively pass portions of the spectrally separated light onto a single photo detector, detecting an AC signal via the single photo detector representative of a gas of interest, weighting a harmonic of the AC signal, and combining the AC signal with the weighted harmonic of the AC signal to detect the gas. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0008]      FIG. 1  is a block diagram of a device for detecting gas according to an example embodiment, 
           [0009]      FIG. 2  is a schematic diagram of a two slit coded filter designed to detect a gas according to an example embodiment. 
           [0010]      FIG. 3  is a diagram illustrating weighted harmonics to account for coded filter misalignment according to an example embodiment. 
           [0011]      FIG. 4  is a diagram illustrated weighted harmonics to account for slit width error in the coded filter according to an example embodiment. 
           [0012]      FIG. 5  is a diagram illustrating a technique for adding weighted harmonics for signal analysis to identify gases according to an example embodiment. 
           [0013]      FIG. 6  is a top view representation of a coded filter as part of a microelectromechanical oscillator according to an example embodiment. 
           [0014]      FIG. 7  is a schematic view of a multiple lane filter according to an example embodiment. 
           [0015]      FIG. 8  is a block circuit diagram illustrating use of weighted harmonics to process signals representative of multiple gases according to an example embodiment. 
           [0016]      FIG. 9  is an alternative block circuit diagram illustrating use of weighted harmonics to process signals representative of multiple gases according to an example embodiment. 
       
    
    
     DETAILED DESCRIPTION 
       [0017]    In the following description, reference is made to the accompanying drawings that form a part hereof, and in which is shown by way of illustration specific embodiments which may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the invention, and it is to be understood that other embodiments may be utilized and that structural, logical and electrical changes may be made without departing from the scope of the present invention. The following description of example embodiments is, therefore, not to be taken in a limited sense, and the scope of the present invention is defined by the appended claims. 
         [0018]    The functions or algorithms described herein may be implemented in software or a combination of software and human implemented procedures in one embodiment. The software may consist of computer executable instructions stored on computer readable media such as memory or other type of storage devices. Further, such functions correspond to modules, which are software stored on storage devices, hardware, firmware or any combination thereof. Multiple functions may be performed in one or more modules as desired, and the embodiments described are merely examples. The software may be executed on a digital signal processor, ASIC, microprocessor, or other type of processor operating on a computer system, such as a personal computer, server or other computer system. 
         [0019]    A device receives light from a source of light after the light has passed through a gas. The gas absorbs some of the light depending on the constituents of the gas. Each gas has its own unique absorption spectra. The light reaching the device has its spectra spread out by wavelength using an optical element, such as a prism or diffraction grating. The optical element spectrally separates incoming light so that shorter wavelengths are directed in one direction and longer wavelengths are directed in a second direction. An opaque filter with slits is then oscillated at a selected frequency, with the slits moving between different frequencies or wavelengths of the spread spectra. 
         [0020]    A photo detector is positioned to receive the light passed by the slits in the filter to measure power of the total amount of light passed. The slits are designed and positioned to pass offsetting AC components of the spectrally separated light passed by at least one not of interest gas, while allowing at least one AC component of the spectrally separated light from a selected gas to be received by the photo detector. In other words, the spectra is dynamically changed in such a way that the contribution to the total (dynamic) signal from the particular spectral components of interest also change dynamically, but the contribution to the (dynamic) signal from the not-of-interest species does not change. In this manner, the signal at the detector is orthogonal to interferrents. 
         [0021]    The photo detector is responsive to the AC component or components of the selected gas to indicate the presence of the selected gas. The amplitude of an AC signal provided by the photo detector will be proportional to the amount of the selected gas that the light has passed through. This process may be referred to as DES (Dynamic Eigen Spectroscopy) due to the enhancement of the AC signal of the gas of interest while AC signals of gases not of interest cancel each other. 
         [0022]    In some embodiments, multiple slits are utilized to cancel AC components from more than one gas, while allowing an AC component from a selected gas to be detected by the photo detector. The filter may include multiple lanes with different sets of slits to detect different gases. Each lane will also be associated with a different photo detector. 
         [0023]    To offset errors in the signal introduced by various deviations from nominal design of the system, one or more harmonics of the AC signals are weighted and combined with the original AC signals. In some embodiments, one of the harmonics includes the phase of the original AC signal. In further embodiments, the phase and amplitude of higher order harmonic signals may be weighted and combined to offset different deviations from nominal, such as misalignment of the slits, a slit width error, errors due to alignment of slits and slit width error based on radius of curvature of species, deviations of the response the photo detector to different wavelengths of light, and variations in the expected absorption spectra of different species. Second, third, and fourth order harmonics may be weighted in some embodiments. Higher order harmonics may be weighted in still further embodiments. 
         [0024]    When the spectra is separated by wavelength and collimated by using real diffraction gratings, prisms, and/or other optical elements, the spectra that reaches the coded filter will often be smeared relative to the spectra that arrives at the detector. This smear will often have a Gaussian shape, but may include linear and/or other components depending on the details of the optical elements, the separation of the optical elements, and the rotation of the slits relative to separation-axis of the spectrally-separated light. While the coded filter may be designed to take some spectral smear into account, signal processing of one or more weighted harmonics may be used to account for deviations of components from nominal, and other variations that may occur from a nominal design. 
         [0025]    In one embodiment, a set of many (e.g. 25) different coded filters, each of which is orthogonal to a known set of spectra (e.g. the spectra of water vapor, carbon dioxide, methane, and ozone) to cancel out signals from such gases. Each of these ˜25 coded filters may be orthogonal and different from one another in a unique way. Each uses different portions of the spectrum to provide orthogonality. When the AC signal is measured from several of the coded filters, it becomes possible to determine the composition and the quantities of several of the chemicals that are in the atmosphere. 
         [0026]    A schematic diagram of a system  100  for detecting gas is shown in  FIG. 1 . A gas  110  may contain many different gasses that absorb light from a source of light  115 . The light source may be a blackbody that emits a spectrally broad source of light. Examples include background terrain that reflects sunlight, or an actual active source, such as a light emitting diode or other artificial source of light. Each different type of gas absorbs different wavelengths of light, while allowing other wavelengths to pass. The gas may operate as a sort of filter. The passed light is received by the system  100  via a lens  120  that directs the light toward a slit  125 . 
         [0027]    Light that proceeds through the slit is collimated by a lens  130  and then spectrally separated by an optical element  135 , such as a prism or diffraction grating. The spectrally separated light is then collimated again via a lens  140  and directed toward a coded filter  145 . The coded fitter is opaque with multiple slits  147  positioned to allow different wavelengths of the spectrally separated light to pass to a photo detector  150 . The slits are oriented parallel to the spectral lines and each extends a selected width of the spectra. The coded filter is oscillated transverse  155  to the spectrally separated light such that AC components of the spectrally separated light are incident on the photo detector  150 . A controller  160  is coupled to the photo detector to receive a signal representative of the amplitude of the light incident on the photo detector  150 . The controller may include an analog amplifier, an analog to digital converter, optional digital weighting functions, and a processor to process the digital signals derived from the photo detector  150  signal. In various embodiments, the controller may integrate the AC signals over a time period that may vary from 0.1 seconds giving a signal to noise ratio of approximately 3:1, to 10 seconds, providing a signal to noise ratio of approximately 30:1. The integration times and signal to noise ratios may vary significantly from embodiment to embodiment, with neither quoted times and ratios being limits. 
         [0028]    The controller  160  also utilizes weighted harmonics of the photodetector generated signal to account for deviations of the components and expected absorption spectra. The signal itself is referred to as a first harmonic. The phase of the first harmonic may be weighted in one embodiment. In further embodiments, the amplitude and phase of higher order harmonics may be weighted to account for a broad array of factors that may affect the accuracy and sensitivity of the system  100 . 
         [0029]    In one embodiment, the coded filter is oscillated at a frequency ω 0 , and light is also sensed at ω 0 . Controller  160  is adapted via programming or circuitry or a combination of both to assign weighting terms to the harmonic components of the optical signal at the photodectector, and add one or more components to a function that utilizes the sensed signal at ω 0 . Various harmonic components may be included in different embodiments. 
         [0030]    By adding higher harmonic k n  sin(nω 0 ) terms to the weighting function it becomes possible to deal with variabilities in the absorption spectra in the different chemicals. The Absorption function A CS  for the C th  chemical and the s th  slit can be expressed as its Taylor function expansion around the mid-point of the slit: A CS =a 0cs +a 1cs x(t)+a 2cs x 2 (t)+a 3cs x 3 (t)+a 4cs x 4 (t)+a 5cs x 5 (t)+a 6cs x 6 (t)+a 7cs x 7 (t), Here, x(t) Amplitude*sin(ω 0 t). The weighting function brings in tremendous flexibility in provide orthogonality, despite processing and packaging variabilities. 
         [0031]      FIG. 2  is a schematic diagram of a two slit coded filter generally at  200 . A light source is indicated at  215  that emits light toward a plume  220  of two gases, A and B. Light that is transmitted through the plume is spectrally separated by an optical elements  225 , and collimated via a lens  230 . A coded filter is represented at  235 , which also illustrates the light absorbed from gases A and B at  240  and  245  respectively. In one embodiment, the filter reduces the total amount of light reaching a photo detector  255  by approximately 10×. Absorption is indicated on a vertical axis with numbers corresponding to absorption at 10 −4 , and the wavelength of the light is indicated on the horizontal axis in nanometers. 
         [0032]    Two slits in the coded filter are indicated at  250  and  251  respectively. The slits are not necessarily to scale. Slit  250  is positioned to move about a portion of the spectra corresponding to gas  240  indicated at  257 , and slit  251  is positioned about a portion of the spectra corresponding to gas  240  indicated at  258 . The coded filter is then oscillated transverse to the spectra a selected distance, left and right as shown, about those portions of the spectra. Typical frequencies of oscillation are between 5 and 10 kHz, but may vary significantly in further embodiments. The oscillation of the coded filter results in the total amount of light reaching the photo detector being modulated at the oscillation frequency f 0 . Note that both of the gas  240  spectra are increasing, resulting in AC signals from gas  240  that add. The signal scales linearly with the amount of gas  240 . This signal will be detected by the photo detector  255 , and passed on to the controller. 
         [0033]    The same two slits  250  and  251  with respect to gas  245  correspond to portions  260  and  261  of the gas  245  spectra. Note that while the spectra is increasing at  260 , it is decreasing at  261 . The light passed by slits  250  and  251  from gas  245  counter each other. The resulting signals from these two areas of the gas  245  spectra effectively cancel each other out, resulting in a net zero AC signal reaching the photo detector  255 . 
         [0034]    When a system is manufactured perfectly such that there are no deviations from a nominally designed structure, constant weighing of the sensed. signal at ω 0  may be used. This results in all parts of the motion of the coded filter being weighted equally. 
         [0035]    When the coded filter is not positioned perfectly within the system, such that the slits are slightly shifted in the spectrum, the resulting signal may be shifted toward red or blue. In such a case, the controller  160  may add a sin(ω 0 ) component as illustrated generally at  300  in  FIG. 3  to the weighting function to compensate for red-blue packaging alignment errors. This changing of the phase of the first harmonic provides the ability to compensate for up to 15 um a mis-alignment with ˜1 nm equivalent physical precision and resolution in one embodiment. These lengths may vary in further embodiments depending on where the relative sizes and spacing of the system components. At  310 , the red portion of the signal is weighted more, whereas at  315 , the blue portion of the signal is weighted more. 
         [0036]    A slit width error or deviation from nominal may be addressed by adding a sin(3ω 0 ) term to the weighting function to compensate for photo-and-etch linewidth errors. This weighting of the third harmonic as illustrated at  400  in  FIG. 4  provides the ability to compensate for 1 um sizing error with ˜1 nm equivalent physical precision and resolution. Again, these sizes are system size dependent and may vary in different embodiments. Weighting the third harmonic provides the ability to weight edges of the signal more as indicated at  410 , or to weight the center of the signal more as indicated at  415 . 
         [0037]    In general, different species of gas, such as species A and species B respond similarly to red-blue shifts and to slit widening. The compensation for the differences described above may not sufficiently distinguish the species. When designing the coded filter, it is the designer&#39;s responsibility to position the slits so that spectra of species A, add spectra of species B have different radii of curvature (2 nd  derivatives) at at least one slit. The controller  160  adds a k 2  sin(2ω 0 ) component to the weighting function to enable exploitation of the different radii of curvature. 
         [0038]    An additional difficulty may arise when the photodetector is not as sensitive to red light as it is to blue light, or vice versa. Further, the light source may not be a black body emitter, and thus may emit differing amounts of red and blue light. These difficulties may be addressed by adjusting the slit widths a-priori to deal with known non-flatness in detector sensitivity and source emissivity. 
         [0039]    Even after adding in harmonics as described above, orthoganalization against all interferrants may not be sufficient. In one embodiment, the controller is adapted to add higher harmonic k n  sin(nω 0 ) terms to the weighting function as needed to deal with variabilities in the absorption spectra in the different chemicals. 
         [0040]    The Absorption function A CS  for the C th  chemical and the s th  slit can be expressed as its Taylor function expansion around the mid-point of the slit: A CS =a 0cs +a 1cs x(t)+a 2cs x 2 (t)+a 3cs x 3 (t)+a 4cs x 4 (t)+a 5cs x 5 (t)+a 6cs x 6 (t)+a 7cs x 7 (t). Here, x(t)=Amplitude*sin(ω 0 t). The weighting function W can be expressed as W=k 0 +k 1  sin(ω 0 t)+k 2  sin(2ω 0 t)+k 3  sin(3ω 0 t)+k 4  sin(4ω 0 t)+k 5  sin(5ω 0 t)+k 6  sin(6ω 0 t)+k 7  sin(7ω 0 t), or as indicated at  500  in  FIG. 5 . Because ∫ sin(nω 0 t) sin(mω 0 t)=0 sinless m=n, very few terms in the product survive. 
         [0041]    Because of Gaussian smear, the contributions of the sin(nω 0 ) terms diminish rapidly as n increases. The Weighting function brings tremendous flexibility in provide orthogonality, despite processing and packaging variabilities. 
         [0042]      FIG. 6  is a top view representation of a coded filter  600  as part of a microelectromechanical oscillator. The filter  600  includes a proof mass  610  that contains multiple slits indicated at  615 . The proof mass includes interdigitated electromagnetic drive fingers  620  to provide force to cause oscillation of the proof mass  610  in a direction illustrated at  617 , transverse to the spectra to be detected. The proof mass  610  is supported by springs  625 . The mass of the proof mass in conjunction with a spring constant of the springs  625  help define a resonant frequency of the proof mass, In various embodiments, the oscillator has a resonant frequency between 50 to 10000 Hertz. In further embodiments, the resonant frequency may be less than 50 or more than 10000 Hertz. 
         [0043]    In further embodiments, the coded filter may be moved by one or more of many different types of mechanisms, including electromagnetic actuators. The speed and frequency of the movement may vary significantly, from less than one hertz to many thousands of hertz or higher, provided resulting AC signals can be detected and processed. 
         [0044]      FIG. 7  illustrates a multiple lane coded filter  700  designed for detection of Sarin gas in air. There are eight lanes illustrated at  712 ,  714 ,  716 ,  718 ,  720 ,  722 ,  724 , and  726 . Each lane has a corresponding photo detector located beneath it, and represented by the same reference number, but not visible in this view. The lanes are formed in a proof mass of a microelectromechanical oscillator as previously described. Filter  700  is designed knowing in advance the spectra of all of the gas species likely to be present. In one embodiment, the gas is Sarin, but may also be GB, GD, and other undesirable gases that are crucial to detect prior to being exposed to humans. Gases that the filter is designed to not respond to include H 2 O, CO 2 , O 3 , CH 4 , and others. 
         [0045]    The controller  160  may also be used to drive the oscillator in some embodiments. In a block circuit diagram illustrated at  800  in  FIG. 8 , a PLL (Phase Locked Loop)  810  may be used to drive a comb oscillator  815  at its resonant frequency. In some embodiments, the resonant frequency determines ω 0 . A value α 0  may be selected to set the overall physical amplitude of the comb driven coded filters. ω 0  is used as the clock frequency for a DDS  820  (Direct Digital Synthesizer), detection electronics  825 , and an analog to digital converter  830  in the controller  160 , which has replicas to equal the number of coded filters in the comb oscillator, eight, in this embodiment. 
         [0046]    In circuit  800 , calibration coefficients k 0  . . . k 5  may be selected to weight the various harmonics components of the detection circuit. This enables a high degree of orthogonalization of the signal-of-interest to each spectral interferrant. The output signal from the DES  825  system is amplified using a high-gain TIA  835  (Trans-Impedance Amp). The amplified output signal is differentially compared it to the DDS  820  waveform using an adder. The result is converted to digital at  830 , and output the signal is provided to a microprocessor  850  or other computing device for identifying one or more gases. The control circuit used to drive the coded filter at resonance is similar to the circuit used in high-precision MEMS gyros in one embodiment. A coherent source and post-detection signal processing may be used. Signals and noise levels in the picoWatt range may be handled. In one embodiment, integration is performed until the S/N&gt;&gt;1. 
         [0047]    An alternative circuit  900  is illustrated in block form in  FIG. 9 . In this embodiment, the raw data from the DES is processed digitally. Jitter may be added to compensate LSB errors. Circuit  900  has more flexibility in dealing with the amplitude of the data, but may have less precision in dealing with timing (spectral resolution) of the data. 
         [0048]    While modeling may be used to calibrate each system, in further embodiments, individual test and calibration of each system will be performed. The physical amplitude of the oscillation of the coded filter may first be adjusted to null the primary interferrent. In one embodiment, a typical amplitude =k 0 ˜10-20 um total motion. The physical amplitude may be controllable to within 1 nm=1E-4. 
         [0049]    Next, the comb resonator is driven at ω 0 , and the signal at ω 0  is sensed and appropriately weighted higher harmonics are calculated. The 2 ω 0  component may be added to compensate for red-blue packaging alignment errors. The 3 ω 0  component may be added to compensate for photo-and-etch linewidth errors. The 4 ω 0  component may be added to compensate for wide slits behaving differently with finite amplitude oscillations than narrow slits. Finally, the 5 ω 0  components may be added as a “knob to turn” to deal with assorted non-linearities and unanticipated non-Eigen components in the spectra. 
         [0050]    The above additions and adjustment may be repeated with test samples of each of the interferrents of interest until optimal orthogonality is achieved. 
         [0051]    In the foregoing description of the embodiments, various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting that the claimed embodiments have more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive subject matter lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby incorporated into the Description of the Embodiments, with each claim standing on its own as a separate example embodiment.