# Patent

## Title
System and method for calibrating remote emissions sensing instruments

## Abstract

A system and method is provided for obtaining calibration curves for CO and CO2 during laboratory calibration of one or more remote emissions sensing (RES) instruments, and for self-calibration of on-road RES instruments to compensate for changes in background CO and CO2 column densities.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This Application claims priority from U.S. Provisional Patent Application Ser. No. 60/604,011, filed Aug. 25, 2004 (now expired), which is incorporated herein by reference in its entirety.


FIELD OF THE INVENTION
The invention relates generally to remote emissions sensing (RES) systems and methods, and more particularly to a system and method for obtaining calibration curves for CO, CO2, and other gases during laboratory calibration of one or more RES instruments, and for self-calibration of on-road RES instruments to compensate for changes in background CO and CO2 concentrations.
BACKGROUND OF THE INVENTION
Remote emissions sensing (RES) systems are known. Such systems are generally utilized for remotely monitoring the exhaust gas composition of motor vehicles. Examples of RES systems and methods are described in, for example, U.S. Pat. Nos. 5,210,702, 5,319,199, 5,401,967, 5,591,975, 5,726,450, 5,797,682, and 5,831,267, each of which is hereby incorporated herein by reference in its entirety.
Despite steady advances in their sophistication and robustness, the manufacture and/or calibration of some instruments (e.g., detectors or detector arrays) utilized in RES systems may suffer from various drawbacks. For instance, one drawback associated with instrument manufacture is the need to separately calibrate each detector channel with known column densities of relevant pollutant gases.
Additionally, when remotely monitoring the exhaust gas composition of motor vehicles on a roadway or driving surface under various operating conditions, it may be expensive and/or inconvenient to have to repeatedly recalibrate on-road instruments using expensive, certified gas cylinders of surrogate exhaust pollutants. Such recalibration may be necessary for conventional instruments to compensate for changes in background (or ambient) CO and CO2 concentrations over time. CO2, for example, is reasonably abundant in the atmospheric optical path used for measurement, and may vary with wind speed, wind direction, and traffic volume, among other factors. These and other drawbacks exist with known calibration methods.
SUMMARY OF THE INVENTION
The invention solving these and other problems relates to a system and method for obtaining calibration curves for CO, CO2, and other gases during laboratory calibration of one or more RES instruments, and for self-calibration of on-road RES instruments to compensate for changes in background CO and CO2 concentrations.
As used herein, a RES system may comprise any components present when remotely monitoring the exhaust gas composition of motor vehicles on a roadway or driving surface at a given test site (or in the “field”), under various operating conditions. Examples of such components may include, but are not limited to, a source, detector (or detector array), transfer optics, certified gas cylinder, imaging unit, speed and acceleration detection unit, thermal detection unit, processor, communicator, or other components. Additionally, one or more RES components may be referred to herein interchangeably as RES instruments. Moreover, as used herein, an optical absorption-based multi-gas analyzer (OABMGA) may comprise one or more RES components (or instruments) including, for example, a source, detector (or detector array), transfer optics, processor, or other components.
One aspect of the invention relates to obtaining calibration curves for CO, CO2, and other gases during laboratory calibration of one or more RES instruments. In one embodiment, a linear absorption curve may be generated for at least one gas in a laboratory setting. In one implementation, for example, gas cell calibration for an HC detector channel may be performed by inserting a gas cell having a known concentration of propane into an optical measurement path. Propane may be used as the surrogate HC. In a laboratory setting, it may reasonably be assumed that the ambient concentration of propane is zero. The HC gas cell calibration results in a linear correlation between absorption and HC column density.
Once the linear correlation between absorption and HC column density has been obtained via HC gas cell calibration, one or more RES components (e.g., source and detector) may be isolated from outside (or ambient) air in the laboratory setting to accomplish a single gas pulse (or puff) calibration. In particular, a known gas mixture of CO, CO2, HC, N2, or other gases may be injected into the optical measurement path from a gas cylinder and raw voltage data for each gas may be acquired from a corresponding detector channel.
The linear absorption curve correlating absorption and HC column density may be used in conjunction with the time-variable puff of gas from the cylinder containing a plurality of gases at known concentration ratios to derive calibration curves for CO, CO2, and other gases, as described in detail herein. In this manner, the need to separately calibrate each detector channel with known column densities of relevant pollutant gases may be avoided.
An additional aspect of the invention relates to self-calibration of on-road RES instruments to compensate for changes in background CO and CO2 concentrations. According to an embodiment of the invention, to be “self-calibrating,” RES instruments may be adapted to continuously monitor the changes in background (or ambient) concentrations of those pollutants which have non-linear, background-dependent calibration curves, so as to remain correctly calibrated.
In particular, the system and method may utilize a combination of signals from detector channels which respond to the pollutants in question and detector channels which provide reference signals. The detector channels may be provided by either a plurality of detectors with different optical filters, or a single detector observing through a plurality of optical filters. In one implementation, the optical filters may include gas cells as in the method known as gas filter correlation spectroscopy. In one embodiment, a first detector channel may monitor CO2 at an IR wavelength of approximately 4.3 micron absorption, and a reference channel may monitor an IR wavelength of approximately 3.9 microns at which no gas absorption is expected. With this capability, a RES instrument may be calibrated once in the laboratory (as described above), and then used to obtain correct readings at virtually any location in the field without recalibration.
According to one implementation, local background pollutant concentrations at a vehicle emissions test site may be determined using a gas puff from a calibration cylinder, or the insertion of a gas cell with known pollutant concentrations. In this regard, any individual channel gain changes which might have taken place since the laboratory calibration may be accounted (or corrected) for. Background levels may be continuously updated based upon voltage signals collected in front of passing motor vehicles whose emissions are to be quantitatively detected. Pollutant background changes as a function of time may be evaluated to identify both slow and fast changes, and the software accepting vehicle emissions readings may be adjusted depending upon how stable the background concentrations are measured to be. In this manner, the expensive and/or inconvenience associated with having to repeatedly recalibrate on-road instruments using expensive, certified gas cylinders of surrogate exhaust pollutants may be avoided.
The various objects, features, and advantages of the invention will be apparent through the detailed description of the preferred embodiments and the drawings attached hereto. It is also to be understood that both the foregoing general description and the following detailed description are exemplary and not restrictive of the scope of the invention.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1A is an exemplary illustration of an optical absorption-based multi-gas analyzer (OABMGA) in a laboratory environment, according to an embodiment of the invention.
FIG. 1B is an exemplary illustration of an optical absorption-based multi-gas analyzer (OABMGA) placed under (or within) an enclosure in a laboratory environment, according to an embodiment of the invention.
FIG. 2 depicts a calibration curve illustrating the linear calibration of an HC channel.
FIG. 3 illustrates a non-linear calibration curve for CO.
FIG. 4 illustrates raw voltages obtained using a calibration method, according to an embodiment of the invention.
FIG. 5 illustrates calibration curves for CO and CO2, according to an embodiment of the invention.
FIG. 6 is an exemplary illustration of a remote emissions sensing (RES) system, according to an embodiment of the invention.
FIG. 7 illustrates a calibration curve, according to an embodiment of the invention.
FIG. 8 illustrates a correlation graph, according to an embodiment of the invention.



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
An aspect of the invention relating to obtaining calibration curves for CO, CO2, and other gases during laboratory calibration of one or more RES instruments will first be described. A description of an exemplary RES system for remotely monitoring the exhaust gas composition of motor vehicles on a roadway or driving surface at a given test site (or in the “field”), under various operating conditions will then be provided, followed by a detailed description of the self-calibration of on-road RES instruments to compensate for changes in background CO and CO2 concentrations.
I. Laboratory Calibration
According to an embodiment of the invention illustrated in FIG. 1A, an optical absorption-based multi-gas analyzer (OABMGA) 100 is provided in a laboratory environment comprising a source 50, one or more detectors or a detector array 70, and processor 110. In various implementations, as described in detail below, additional components may be provided.
Source 50 may comprise one or more sources of electromagnetic radiation (ER) which may be used in the absorption spectroscopy measurement of various gases in an optical path 40 (including, for example, components of vehicle exhaust emissions) in a known manner. Source 50 may comprise an infrared (IR) radiation source. In alternative embodiments, other types of radiation sources may be used including, for example, an ultraviolet (UV) source, a visible light source, or other suitable sources as known and understood by those having skill in the art. In some embodiments, a combination of radiation sources may be used.
Detector array 70 is preferably chosen to permit detection of electromagnetic radiation emitted by source 50. For example, detector array 70 may comprise a photodetector (e.g., a photodiode), a photomultiplier tube (PMT), a spectrometer, or any other suitable radiation detector. A mercury cadmium telluride (Hg—Cd—Te) or lead selenide photodetector may also be used to detect IR radiation. Other suitable detectors or detector arrays or combinations thereof may also be used. For instance, in some implementations, a single detector with multiple filters may be utilized instead of an array employing multiple detectors. The multiple filters may be moveable, such as spinning filters, to allow multiple components of a gas sample (e.g., vehicle exhaust components) in optical path 40 to be detected. In this regard, a single detector can be employed to detect a plurality of different components because each of the moveable filters is designed to allow only the wavelength band of interest by a particular component to pass to the detector.
According to an embodiment of the invention, processor 110 may comprise a suitable processing device such as, for example, a computer or other microprocessor. Processor 110 may optionally employ software to accomplish desired analysis of collected and/or stored data in a known manner (subject to the improvements described herein). For example, with regard to a gas sample present in optical path 40, software may be used to calculate the relative amounts of various constituents of the gas sample, concentrations of various constituents of the gas sample (e.g., HC, CO2, NOx, CO, etc.), and the decay rate (e.g., dissipation in time) of the constituents, among other things. Processor 110 may further comprise or interface to a graphical user interface (GUI) for enabling a user to access, navigate, and otherwise utilize the various software processing functions described herein.
According to an embodiment of the invention, source 50 may be configured to pass a beam of EM radiation through a gas sample in optical path 40. Detector array 70 may be configured to receive the beam after it passes through the gas sample. One or more filters (not illustrated) may be associated with detector array 70 to enable detector array 70 to determine the intensity of EM radiation having a particular wavelength or range of wavelengths. The wavelengths may be selected to correspond to wavelengths absorbed by molecular species of interest in the gas sample (e.g., hydrocarbons (HC), carbon monoxide (CO), carbon dioxide (CO2) and nitrogen oxides (NOX) such as NO and NO2). One or more detector output voltages represent the intensity of the EM radiation measured by that detector. Detector array 70 may be configured such that each molecular species of interest in the gas sample has a unique detector channel.
The detector output voltages may then be input to processor 110. Processor 110 may calculate the difference between the known intensity of source 50 and the intensity detected by the detectors to determine the amount of absorption by the particular molecular species (based on predetermined wavelengths associated with that species). Based on the measured absorption(s), the column density of one or more molecular species in the gas sample may be determined in a known manner.
As recited above, one of the less satisfactory aspects of RES instrument manufacture is the need to separately calibrate each detector channel with known column densities of relevant pollutant gases. According to an embodiment of the invention, a system and method for overcoming this drawback is provided, wherein at least one hydrocarbon detector channel is carefully calibrated with propane (or other suitable gas), and the resulting linear calibration curve is used to establish calibration equations for a predetermined number of other detector channels via a single pulse (or puff) of gas with a certified mixture of known concentrations. Alternate configurations (e.g., one or more puffs for varying predetermined time periods) may be utilized.
According to an embodiment of the invention, for calibration, a predetermined linear absorption curve may be generated for at least one gas. For example, in one implementation, gas cell calibration for an HC detector channel may be performed by inserting a gas cell (not illustrated) having a known concentration of propane (or other suitable gas), into optical path 40. Propane may be used as the surrogate HC. Gas-cell calibration for an HC detector channel may be performed with greater reliability as it can be reasonably assumed that, in a laboratory setting, the ambient concentration of propane is zero.
The HC gas cell calibration results in a linear correlation between absorption and HC concentration. FIG. 2, for example, depicts a calibration curve illustrating the linear calibration of an HC detector channel. The illustrated linear response provides the key to identifying the instantaneous CO and CO2 column densities during a single gas pulse (or puff) calibration in a laboratory setting, as described in detail below. The use of single gas pulse calibration may assist in avoiding the effort previously expended to obtain the non-linear calibration curves for CO and CO2. FIG. 3, for example, is an illustration of a non-linear calibration curve for CO. The illustrated curve was fit with a second order polynomial.
In one embodiment, once the linear correlation between absorption and HC concentration has been obtained via HC gas cell calibration, source 50, detector 70, and optical path 40 may be isolated from outside (or ambient) air in the laboratory setting to accomplish single gas pulse (or puff) calibration. As illustrated in FIG. 1B, for instance, source 50 and detector 70 may be placed under (or within) enclosure 150. Enclosure 150 may comprise a plastic enclosure or other suitable covering device or mechanism for isolating source 50 and detector 70 from outside (or ambient air). In various implementations, processor 110 may be located either within or outside of enclosure 150.
According to an embodiment of the invention, prior to single pulse calibration, enclosure 150 may be purged with a gas from purge gas source 160 over a pre-determined period of time to purge out all of the CO, CO2, HC, or other gases within optical path 40. In one implementation, the purge gas may comprise N2, and enclosure 150 may be purged with N2 during an overnight time period at a rate of approximately 1.0 L/min. Other time periods and purge rates may be utilized.
In one implementation, gas lines from gas source 170 (which may comprise a cylinder containing, for instance, a known gas mixture of CO, CO2, HC, N2, or other gases) may also be carefully purged with fresh gas not entering enclosure 150.
At the conclusion of the predetermined “purge” time period, source 50 and detector 70 may be powered on, and processor 110 may be configured to acquire raw voltage data from detector 70 for a predetermined acquisition duration. In one embodiment, the predetermined acquisition duration may comprise ten seconds. Once the predetermined acquisition duration has commenced, gas source (or cylinder) 170 may be opened (e.g., through the use of a solenoid or other valve) and a known gas mixture of CO, CO2, HC, and N2, or other gases may be passed into optical path 40 within enclosure 150. Any of the valves illustrated in FIG. 1B may be controlled manually, or in an automated manner via processor 110.
According to one embodiment of the invention, the gas mixture may be introduced into optical path 40 up to optical depths (the product of concentration times path length also termed column density) comparable to those seen on a roadway or other test surface in a real study (of vehicle emissions). This may be done such that the on-road data can be interpolated within the laboratory calibration rather than extrapolated outside the bounds of the calibration.
FIG. 4 is an exemplary illustration of raw voltages obtained from a reference detector channel, in addition to CO, HC, and CO2 detector channels over an acquisition duration of ten seconds. In this embodiment, voltage outputs were sampled every ten milliseconds for the ten second time period. Other sampling rates may be utilized. As shown, the gas mixture appears to have been introduced in a single puff into optical path 40 from approximately 1.7 seconds to 5. seconds during the total ten second time period.
Once data have been collected, processor 110 may process the data to produce a graph of absorbance vs. gas concentration in the following manner. In one implementation, all of the raw voltages may be ratioed to the reference voltage to negate any effect of electrical noise or source fluctuation using the following equation:




  
    
      
        
          R
          n
        
        =
        
          
            V
            Signal
          
          
            V
            Reference
          
        
      
    
    
      
        (
        1
        )
      
    
  




    
    
        where R, is the reference ratioed signal voltage;
        VSignal is the observed voltage on the detector channel of interest (CO, CO2, or HC); and
        VReference is the voltage on the reference channel.
    
    


In one implementation, all voltages may have subtracted a “dark” voltage obtained from the readings when the IR “zero” beam from source 50 is blocked. The cleanest air (e.g., the voltage readings at the beginning in FIG. 4 before the gas pulse (or puff) started), may be assumed to be when the gas concentration is equal to zero. Thus, these are referred to as the “R0” values, or the reference ratioed voltages when the air is free of pollutant.
The absorbance of the system may be calculated using the equation:




  
    
      
        A
        =
        
          ln
          ⁢
          
            
          
          ⁢
          
            
              I
              0
            
            I
          
        
      
    
    
      
        (
        2
        )
      
    
  




    
    
        where A is absorbance;
        I0 is the light intensity absent pollutant; and
        I is the measured light intensity when to pollutant is present.
    
    


In one embodiment, the absorbance data and concentration data may be established into a mathematical relationship to form a calibration equation. The concentration ratios HC/CO2 and HC/CO may be read from the gas source 170 (or cylinder), depicted in FIG. 1B, used for calibration. These ratios may be used for the calculation of real-time concentrations of CO and CO2 gases because the previously performed HC calibration enables direct calculation of the instantaneous HC levels during the pulse (or puff).
The HC calibration is the line whose slope is used to calculate the instantaneous HC concentration, and thus the CO2 concentration and CO concentration values with the aid of the gas cylinder (or gas source 170) ratios. From the previous data, the calibration curves for CO2 and CO may be obtained directly—calibrated using the known concentration values from the gas cylinder (or gas source 170). FIG. 5 is an exemplary illustration of calibration curves generated for CO2 and CO.
According to an embodiment of the invention, with this graphical representation of the data (e.g., FIG. 5), a best-fit line may be applied to the data to mathematically relate the two axes and generate a calibration equation for CO and a calibration equation for CO2. This fit may be carried out in at least two ways; the first is a second order polynomial fit (quadratic equation), and the other is an exponential growth curve (in the form of y=a+eb+cx). Other arithmetic curve fitting equations may be used as known and understood by those having skill in the art.
When the laboratory calibration is completed, specific instrument software may be updated with the appropriate laboratory calibration equation for CO and the appropriate laboratory calibration equation for CO2, and moved into a field measurement. The foregoing system and method, based on a single pulse (or puff) calibration, is less labor intensive than calibration performed one gas at a time.
Although the foregoing description was directed primarily to obtaining calibration curves for CO and CO2, it should be recognized that the system and method described may be used to obtain calibration curves for other gases of interest in the same manner, thus avoiding the need to separately calibrate each detector channel with known column densities (optical depths) of relevant pollutant gases. As such, the foregoing description should not be viewed as limiting.
II. Exemplary Remote Emissions Sensing (RES) System
Prior to setting forth a description of the self-calibration of on-road RES instruments to compensate for changes in background CO and CO2 concentrations, a description of an exemplary RES system is provided with reference to FIG. 6. The calibration methods disclosed herein can be used with this or other system configurations. Embodiments of the invention may include some or all of the RES system components as described below, or other components.
According to one implementation, the RES system measures emissions in an exhaust plume 30 (from a motor vehicle 10) in an optical (or measurement) path 40 on a roadway 20. Roadway 20 may comprise a single or multi-lane roadway, or any other roadway or driving surface suitable for the safe passage of vehicle 10 under various operating conditions. Alternatively, roadway 20 may comprise a test lane designated for vehicle emissions testing, wherein vehicle 10 may be tested under a variety of operating conditions.
As previously recited, source 50 may comprise one or more sources of electromagnetic radiation (ER) and detector array 70 is preferably chosen to permit detection of electromagnetic radiation emitted by source 50.
According to one embodiment, the RES system may comprise transfer optics 60 mounted in a manner to allow radiation from source 50 to be reflected to detector array 70 for analysis. Transfer optics 60 may comprise a mirror, flat mirror, lateral transfer mirror (LTM), vertical transfer mirror (VIM), retroflector, or other device. In one embodiment, transfer optics 60 may comprise a lateral transfer mirror to reflect radiation from source 50 along a path displaced laterally or vertically, depending on orientation, from the incident direction. Other configurations may be used.
Processor 110, as described in detail above, may employ software to accomplish desired analysis of collected and/or stored data, and to carry-out one or more of the calculations described in the laboratory calibration methods and self-calibration methods set forth herein.
With regard to an exhaust gas sample present in optical path 40, software may be used to calculate the relative amounts of various exhaust gas constituents, concentrations of various exhaust gas constituents (e.g., HC, CO2, NOx, CO, etc.), the decay rate (e.g., dissipation in time) of the exhaust constituents, and the opacity of the exhaust plume.
According to an embodiment of the invention, processor 110 may calculate the relative amounts of various exhaust gas constituents in exhaust plume 30 of vehicle 10 (as vehicle 10 passes by (or through) the RES system) by computing the ratio of the absorption for a particular exhaust gas constituent to the CO2 absorptions, in a known manner. This eliminates the need to calculate the total amount of exhaust plume present since the calculated ratios may provide sufficient information to identify vehicles which do not meet predetermined pollution criteria.
Processor 110 may also comprise software to accomplish other data analysis functions. For example, vehicle emission data may be checked for running losses. Running losses may typically include emission readings due to fuel system leaks on a vehicle (e.g., leaky fuel tank filler cap, fuel line, etc.), blow-by emissions (e.g., crank case emissions blowing by the piston rings), emissions due to other vehicles in the vicinity, or other systematic losses.
Processor 110 may also include software to accomplish various vehicle owner notification functions. For example, the owner of a vehicle that has been recorded as being in compliance with certain predetermined emission levels may receive a notification. Coordination with local authorities may be arranged to grant vehicle owners a waiver or pass of local emission certification procedures upon receiving such a notification. Likewise, vehicles that fail to meet predetermined emission levels may receive a notification requiring the owner to remedy the non-compliance. Other data processing functions are also possible. For example, processor 110 may control the release of gas from certified gas cylinder 140 for in-field calibration as described in greater detail below. Processor 110 may also interface to, control, and/or collect and reduce data from an imaging unit 80, a speed and acceleration detection unit 90, and thermal detection unit 130.
In various embodiments, the RES system may comprise an imaging unit 80 to capture and/or record an image of vehicle 10 passing by (or through) the RES system in a known manner. Imaging unit 80 may be positioned to record an image of vehicle 10 at any predetermined number of locations. Imaging unit 80 may comprise, for example, a film camera, video camera, or digital camera. Other imaging devices may also be used.
Preferably, imaging unit 80 may record an image of the identification tag (e.g., license plate) of vehicle 10. Tag information may be processed by processor 110 to provide additional information about the vehicle. For example, Motor Vehicle Department databases may be accessed to retrieve owner information, make, model type, model year, or other information. In some embodiments, this additional information may be incorporated into the emission sensing data analysis. For example, the make and model year of the vehicle may be used to determine input information for certain processing steps, including information such as whether the vehicle includes a carburetor or fuel injector, whether the car runs on diesel fuel or gasoline, etc.
According to an embodiment of the invention, the RES system may include a speed and acceleration detection unit 90. Preferably, the speed and/or acceleration of vehicle 10 may be measured as it passes through the RES system using speed and acceleration detection unit 90 in a known manner.
In one embodiment, speed and acceleration detection unit 90 may comprise an arrangement of laser beams or other light beams associated with timing circuitry. The laser or light beams may be arranged to traverse the path of vehicle 10 at various points. As vehicle 10 passes, it will cause interruptions in the laser or light beams. The times at which the beam interrupts occur may be used to calculate the vehicle's speed and/or acceleration. Other methods of determining vehicle speed and/or acceleration may also be used or incorporated into the RES system.
Alternatively, the laser or light beams may be arranged to traverse the path of vehicle 10 at a single point in the vehicle's path. For example, radar systems may be used to determine vehicle speed and acceleration. Alternatively, transducers, piezoelectric elements, or other “drive over” detectors may be placed at locations in the roadway to monitor vehicle passage. Preferably, speed and/or acceleration data may be input into processor 110 to help characterize vehicle operating conditions (e.g., accelerating or decelerating), or to determine which vehicle is to be associated with a particular sensor measurement. Other configurations and uses of speed and acceleration data are also possible.
Some embodiments of the invention may incorporate a thermal detection unit 130. Preferably, thermal detection unit 130 may comprise a non-contact thermometer system. For example, an IR thermometer may be used to optically detect the temperature of remote objects. Other temperature detection systems may also be used. Thermal detection unit 130 may, for example, be used to detect the temperature of portions of the vehicle passing through the RES system. Some embodiments may use direct sensing of the area of interest. For example, an IR thermometer may be aimed at the underside of a passing vehicle to detect the temperature(s) of vehicle components (e.g., engine, catalytic converter, muffler, etc.). Indirect sensing may also be used. For example, an IR thermometer may be aimed at the roadway to measure the heat of the passing vehicle which is reflected from the roadway surface.
Thermal information that is detected by thermal detection unit 130 may be used to indicate that an engine has just recently been started (e.g., the engine is “cold” or has not reached normal operating temperature). Such a cold engine reading may be used, for example, to initiate an alternative data processing routine. Certain embodiments of the invention may reduce the chance of a potentially misleading reading by also detecting the temperature of other portions of the vehicle. Other uses for collected thermal data are also possible. Thermal detection of the exhaust plume of a vehicle and/or ambient temperatures may also be used in connection with various aspects of the invention.
According to one embodiment of the invention, an identification tag on vehicle 10 may be read to identify the vehicle and associate particular sensed vehicle emission information with the vehicle. An identification tag, defined as a license plate above, may also comprise a transponder located on or within vehicle 10 (e.g., hung from a rear view mirror, placed on the dashboard, etc.), or that is integral within the vehicle (e.g., part of a global positioning system (“GPS”), located within the engine of the vehicle, or placed or mounted elsewhere). The transponder may transmit information about vehicle 10, including make and model of vehicle 10, engine characteristics, fuel type, the owner of vehicle 10, or other information which may be pertinent. According to an embodiment of the invention, a transponder may be used in connection with other functions. For instance, a transponder may also be used in connection with a toll pass, whereby a driver can electronically pay tolls via the transponder without stopping the vehicle.
An identification tag may also comprise a tag or decal that requires a reader. By way of example, an identification tag may comprise a decal with identifying marks (e.g., bar codes, infrared markings, etc.) containing information about vehicle 10. The decal may be located outside vehicle 10, such as on a front or rear bumper, on the under-side of vehicle 10, or any other location on vehicle 10 where the decal may be suitably read. A reader may observe the decal and thereby obtain information about vehicle 10. One embodiment employs a bar code placed on the roof of vehicle 10, which can be read by a reader placed above vehicle 10.
A receiver may be used to obtain information from an identification tag. According to an embodiment of the invention, an antenna may receive signals transmitted from an identification tag containing a transponder. Any type of conventional receiver may be used to receive signals. According to an embodiment of the invention, one reader and/or receiver may be used in connection with multiple lanes. Based on the signal received or the decal read, a receiver or reader may determine in which lane a particular vehicle is located at a particular time.
Processor 110 may receive information about vehicle 10 from a reader and/or receiver. According to an embodiment of the invention, processor 110 may receive vehicle information. Vehicle information and information obtained by sensing vehicle emissions may be stored. Processor 110 may correlate vehicle information received from an identification tag with the results from vehicle emissions sensing. Processor 110 may update a vehicle record to account for the results obtained by processing vehicle emission data, such as information regarding whether a vehicle has passed or failed predetermined emissions criteria.
According to an embodiment of the invention, the RES system may further comprise a communicator 120. Communicator 120 may communicate information such as, for example, measured vehicle emissions and identification tag information from the RES system to various other locations (e.g., Motor Vehicle Departments, a central data repository, servers, etc.) for storage, processing, viewing, or other use in a known manner. Communicator 120 may transmit and/or receive information via a wire connection, such as cable or telephone line, or a wireless connection, such as by a radio, cellular, or satellite transmitter, or via any other type of suitable wireless communication.
In some embodiments, communicator 120 may comprise appropriate hardware and/or software to enable processor 110 to be accessed remotely over a network (not illustrated) via a communications link (not illustrated). The network may include any one or more of, for instance, the Internet, an intranet, a PAN. (Personal Area Network), a LAN (Local Area Network), a WAN (Wide Area Network), a SAN (Storage Area Network), or a MAN (Metropolitan Area Network). The communications link may include any one or more of, for instance, a copper telephone line, a Digital Subscriber Line (DSL) connection, a Digital Data Service (DDS) connection, an Ethernet connection, an Integrated Services Digital Network (ISDN) line, an analog modem connection, a cable modem connection, or a wireless connection. In this regard, a user (e.g., an emissions test administrator or other individual) at a remote computer terminal can administer emissions tests, and/or analyze or process data. Thus, the RES system may, in various embodiments, comprise either manned or unmanned systems.
As recited above, alternative RES system configurations may exist incorporating some or all of the aforementioned system components. Accordingly, the RES system illustrated in FIG. 6 should not be viewed as limiting.
As an example, in certain embodiments (not illustrated), source 50 and detector 70 may be placed on opposite sides of roadway 20. Various components of speed and acceleration detection unit 90 and thermal detection unit 130 may also be positioned on opposite sides of roadway 20.
In another embodiment (not illustrated), the RES system may comprise a compact, unmanned system that may be used for unattended monitoring of vehicle emissions data (also referred to as a “bunkered” unit). In such an embodiment, source 50, detector 70, imaging unit 80, processor 110, communicator 120, and various components of speed and acceleration detection unit 90 and thermal detection unit 130 may be housed together on a first side of roadway 20, while transfer optics 60 and various other components of speed and acceleration detection unit 90 and thermal detection unit 130 may be housed together on the opposite side of roadway 20. Other configurations are possible. Emissions and other measured data may be transmitted by communicator 120 as described in detail above.
III. Self-Calibration
In certain instances, when remotely monitoring the exhaust gas composition of motor vehicles on a roadway or driving surface under various operating conditions using a RES system similar to that described in FIG. 6, it may be expensive and/or inconvenient to have to repeatedly recalibrate on-road RES instruments using expensive certified gas cylinders of surrogate exhaust pollutants.
Such recalibration may be required for conventional instruments because the same measured voltage change (even when converted to absorbance) on the CO2 and CO channels can correspond to a different amount of gas at different times of the day because the absorbance versus column density graph is non-linear (e.g., the Lambert-Beer absorption law does not hold). Stated in another way, one problem is that, where signals fall on the absorption curves referenced above, depends upon the background concentrations of the relevant pollutant. CO2 in particular is reasonably abundant in the atmospheric optical path used for measurement, and may vary with wind speed, wind direction, and traffic volume, among other factors.
According to an embodiment of the invention, to be “self-calibrating,” RES instruments may be adapted to continuously monitor the changes in background (or ambient) concentrations of those pollutants which have non-linear, background-dependent calibration curves, so as to remain correctly calibrated.
In particular, a combination of signals from detector channels which respond to the pollutants in question, and detector channels which provide reference signals, may be utilized. The detector channels may be provided by either a plurality of detectors with different optical filters, or a single detector observing through a plurality of optical filters. In one implementation, the optical filters may include gas cells as in the method known as gas filter correlation spectroscopy. In one embodiment, a first detector channel may monitor CO2 at an IR wavelength of approximately 4.3 micron absorption, and a reference channel may monitor an IR wavelength of approximately 3.9 microns at which no gas absorption is expected. With this capability, a RES instrument may be calibrated once in the laboratory (as described above), and then used to obtain correct readings at virtually any location in the field without recalibration.
According to one implementation, local background pollutant concentrations at a vehicle emissions test site may be determined using a gas puff from cylinder 140, or the insertion of a gas cell with known pollutant concentrations. In this regard, any individual channel gain changes which might have taken place since the laboratory calibration may be accounted (or corrected) for. Background levels may be continuously updated based upon voltage signals collected in front of passing motor vehicles whose emissions are to be quantitatively detected. Pollutant background changes as a function of time may be evaluated to identify both slow and fast changes, and the software accepting vehicle emissions readings may be adjusted depending upon how stable the background concentrations are measured to be.
According to an embodiment of the invention, detector array 70 may be configured such that each gas (or constituent) of exhaust plume 30 has a unique detector channel. For a gas of interest, a band-pass filter may be configured to pass radiation having a particular wavelength or range of wavelengths to a detector. The particular wavelength or range of wavelengths may be selected to correspond to wavelengths absorbed by the gas of interest. For a given band-pass filter channel for one gas, there is one absolute calibration curve of absorption vs. gas column density (e.g., similar to that shown in FIG. 5 for both CO and CO2). In a laboratory setting, the filter may be calibrated all the way from no gas in the optical path, up to a practical limit of feasible gas column density that may be present on the road.
Currently, in the field, a single gas puff calibration may be used to determine the background concentration of CO2 and the theoretical “R0” value. As recited above, the R0 value may comprise the reference ratioed voltage when the air is free of pollutant. Although not illustrated, a gas cell with a known admixture of pollutants may also be used for this application. Although an R0 value was calculated in the laboratory calibration, the CO2 voltage to Reference voltage ratio may change irrespective of one another if gain potentiometers are changed in the circuitry of the instrument, or if other changes to the hardware are made. Accordingly, it is sometimes appropriate that the new R0 be calculated in the field. This recalibration may negate any of the aforementioned changes (or other changes) from adversely affecting the absolute calibration, and may enable freedom of hardware modifications without need for redoing the laboratory calibration.
In one implementation, to perform this single field calibration, a half-second (or other duration) gas puff from certified gas cylinder 140 may be made in optical path 40 (in the field) and, using known, absolute calibration curves, processor 110 may calculate what the R0 and CO2 background “CO2B” is at the test site (in the field). Certified gas cylinder 140 may contain, for example, a known mixture of CO, CO2, and HC (similar to the lab calibration described above). The calculation procedure, according to one embodiment, is set forth below.
According to an embodiment of the invention, a preliminary assumption is that all column densities (in the case of a half-second there may be for instance fifty readings, one for each 10 ms averaging period) calculated are assumed to be in an 8 cm. cell, although the actual path length may be much larger.
All voltages may first be normalized to reference, providing Rns for all channels and readings, where Rns is (Vsignal-Vdark)/(Vref-Vdark) for each time (n) and substance/pollutant(s) [see equation 1]. Vdark readings may be measured either when the optical beam is blocked by a vehicle, or at the beginning of the day when blocked by an operator.
Next, the “clean air” or air with the smallest absorptions apparent in the puff may be located in time. This time may comprise the time of the highest CO2 voltage of the fifty readings (least absorption). These “smallest” ratios may be termed RHClow, RCOlow, and RCO2low. The time of these may be defined as Tmin.
The lowest HC voltage of the fifty readings may then be identified, which corresponds to the highest HC column density, and is termed RHChi. RHChi, together with RCOhi and RCO2hi may all be set at that same point in time, Tmax.
ΔHC may then be calculated using the following equation:




  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          HC
        
        =
        
          
            
              H
              1
            
            ⁡
            
              (
              
                
                  R
                  HClow
                
                -
                
                  R
                  HChi
                
              
              )
            
          
          
            R
            HChi
          
        
      
    
    
      
        (
        3
        )
      
    
  




    
    
        wherein ΔHC is the apparent change in HC column density;
        H1 is the laboratory slope for the linear HC calibration (e.g., FIG. 2); and
        RHCloW and RHChi are as defined above.
    
    


This calculation may be performed for every point along the puff data where Rn≠R0, and wherein Rn is as defined above (Rns).
Next, using the ratios of CO/HC and CO2/HC from certified gas cylinder 140 used during the field calibration, the “true” column densities of CO and CO2 in optical path 40 at the time of field calibration may be calculated. This can be achieved because from the ΔHC data and the cylinder ratios, the ΔCO and ΔCO2 may be calculated at each of the fifty time points during the calibration puff.
Overall, the ΔCO2 is equal to the CO2 column density at Tmax minus the CO2 column density at Tmin:

ΔCO2═CO2Max−CO2Min   (4)

    
    
        wherein CO2Max is the CO2 column density at Tmax; and
        CO2Min is the CO2 column density at Tmin.
    
    


Knowing that there was background CO2 or “CO2B” even at what is assumed to be the “clean air” column density at Tmin, it is assumed that CO2B is equal to the CO2 column density at Tmin:

CO2B ═CO2Min   (5)

The field puff has given a known ΔCO2 and an apparent absorbance for each of the 50 points during the half-second puff. Processor 110 may then determine at which location along the absolute laboratory calibration curve the instrument is currently operating. In the case where the absolute calibration equation is quadratic in form, the following algebra may be used:

CO2Max═C1abshi+C2abshi2   (6)

CO2B═C1abslo+C2abslo2   (7)





  
    
      
        
          
            
              CO
              2
            
            ⁢
            B
          
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            
              CO
              2
            
          
        
        =
        
          
            
              C
              1
            
            ⁡
            
              (
              
                
                  
                    R
                    0
                  
                  
                    R
                    hi
                  
                
                -
                1
              
              )
            
          
          +
          
            
              
                C
                2
              
              ⁡
              
                (
                
                  
                    
                      R
                      0
                    
                    
                      R
                      hi
                    
                  
                  -
                  1
                
                )
              
            
            2
          
        
      
    
    
      
        (
        8
        )
      
    
  




The terms C1 and C2 are coefficients derived from the laboratory calibration. It should also be recognized that A* has been replaced above with “abs” (standing for absorbance). These absorbance numbers are a function of the known Rn values in the calibration and the unknown, but calculable R0.
In equations (6)-(8):

abshi is (R0/Rhi−1);

abslo is (R0/Rlo−1); and

Rhi is as defined above.
The above derivation results in two simultaneous equations with two unknowns (R0 and CO2B). They are solved for both variables. Algebraic derivation yields a quadratic solution:




  
    
      
        
          R
          0
        
        =
        
          
            -
            b
          
          ±
          
            
              
                
                  b
                  2
                
                -
                
                  4
                  ⁢
                  a
                  ⁢
                  
                    
                  
                  ⁢
                  c
                
              
              
                2
                ⁢
                a
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        Where
        ⁢
        
          :
        
      
    
    
      
        
      
    
  
  
    
      
        a
        =
        
          
            C
            2
          
          ⁡
          
            (
            
              
                
                  R
                  lo
                
                
                  R
                  hi
                
              
              -
              
                
                  R
                  hi
                
                
                  R
                  lo
                
              
            
            )
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        b
        =
        
          
            (
            
              
                C
                1
              
              -
              
                2
                ⁢
                
                  C
                  2
                
              
            
            )
          
          ⁢
          
            
              R
              lo
            
            
              R
              hi
            
          
        
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        c
        =
        
          
            -
            
              R
              lo
            
          
          ⁢
          
            R
            hi
          
          ⁢
          Δ
          ⁢
          
            
          
          ⁢
          
            CO
            2
          
        
      
    
    
      
        (
        12
        )
      
    
  




In equations (10)-(12), Rlo is as defined for Rlow above.
The numerical solution for R0 (the theoretical CO2 to reference voltage ratio that would have been observed had the CO2 column density equaled zero) may be used to solve for CO2B at the time of field calibration. The solution becomes:




  
    
      
        
          
            CO
            2
          
          ⁢
          B
        
        =
        
          
            
              C
              1
            
            ⁡
            
              (
              
                
                  
                    R
                    0
                  
                  -
                  
                    R
                    lo
                  
                
                
                  R
                  lo
                
              
              )
            
          
          +
          
            
              
                C
                2
              
              ⁡
              
                (
                
                  
                    
                      R
                      0
                    
                    -
                    
                      R
                      lo
                    
                  
                  
                    R
                    lo
                  
                
                )
              
            
            2
          
        
      
    
    
      
        (
        13
        )
      
    
  




The CO2B and R0 values may be stored in a data file or database (not illustrated) operatively connected to processor 110. R0 is used and CO2B is updated until the instrument is recalibrated, which may occur, for instance, when it is next setup at another site.
According to an embodiment of the invention, the “cleanest air” CO2to reference voltage ratio is used to determine a specific CO2B for each reading of a vehicle's exhaust emissions. Thus, CO2B can change as each car passes. This allows for a quadratic or other non-linear equation to be used to correct for CO2 for each passing vehicle. This remains an empirical method. The absorption vs. column density data may be manipulated in any mathematically equation to conveniently relate the two. A quadratic equation may not always be a good fit to the empirical data. An exponential growth curve, for instance, may turn out to be a better fit.
FIG. 7 illustrates a calibration curve wherein the exponential growth curve used as the empirical fit is almost indistinguishable from the calibration data, whereas the quadratic fit of the data grossly underestimates the column density of CO2 near x=0. The poor quadratic fit may arise partly because most of the 1000 data points are above a column density of 0.4 atm*cm of CO2.
According to an embodiment of the invention, an exponential growth curve can be fit with the change of absorbance and the concentration input from certified gas cylinder 140, and iteratively solved for a CO2B and an R0 in the field much more quickly and effectively than with closed form solutions. As such, processor 110 may be input with the laboratory-calculated curves, and may iteratively fit the field calibration data to solve for the CO2B and an R0.
The treatment of the data is quite similar in both quadratic and exponential methods, in that the zero-corrected voltage data may all be reference ratioed to acquire Rn values. From these Rn values, the “cleanest” and “dirtiest” air data points may be found and set as the respective Rlo and Rhi values. Natural log (Rlo/Rn) may then be calculated as the absorbance value to be used in calibration. The HC data may be used, and calculated back against the laboratory line. ΔHC may be calculated at each point, thus giving the amount of HC added during the calibration puff. Using the cylinder ratios of HC/CO2 and HC/CO, the ΔCO2 and ΔCO may be calculated as well.
In the first iteration to be tested, it may be assumed that Rlo═R0 for the instrument. Using the exponential growth curve from the laboratory and the known Rhi & Rlo (assumed temporarily to be R0), the first approximation of ΔCO2 may be calculated using the equation:




  
    
      
        
          [
          
            CO
            2
          
          ]
        
        =
        
          a
          +
          
            e
            
              b
              +
              
                c
                ⁡
                
                  (
                  
                    ln
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        R
                        lo
                      
                      
                        R
                        hi
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




where [CO2] is the CO2 column density. Since the assumed CO2B is zero, CO2hi−ΔCO2—CO2B<0. If that answer comes out, as expected, to be negative, the iteration is continued, and a new iteration may be started at a point CO2B═CO2Bprevious+0.0064 atm*cm. Using this new CO2B, the inverse calibration equation is solved:




  
    
      
        
          ln
          ⁡
          
            (
            
              
                R
                lo
              
              
                R
                0
              
            
            )
          
        
        =
        
          
            
              ln
              ⁡
              
                (
                
                  
                    [
                    
                      CO
                      2
                    
                    ]
                  
                  -
                  a
                
                )
              
            
            -
            b
          
          c
        
      
    
    
      
        (
        15
        )
      
    
  




A new R0 is thus determined. Using this R0, Rlo, and Rhi, the ΔCO2 can be calculated. Again, if CO2hi—CO2B—CO2added<0, the iteration is continued, and a new one begins once 0.0064 atm*cm is added to the initial CO2B.
This iteration may be continued until CO2hi—CO2B−CO2added≧0. At this point, the iteration may be stopped and the R0 used for that iteration is logged into a computer calibration file, as it is the accepted value for the R0. Once the proper R0 is calculated, an on-screen display may illustrate the full curve of [CO2]calculated vs. [CO2]input. Since the instrument calculates the column density of CO2 based on the calculated R0, the “calculated” side represents the instrument's analysis of what is in the cylinder. The “input” axis may be obtained using the HC and the known cylinder ratios to determine how much CO2 was actually put into the instrument. Therefore, the graphs (both for CO2 and CO) should be linear with a slope of one.
FIG. 8 is an exemplary illustration of a correlation graph that may appear on-screen during a successful on-road gas puff calibration. If such a graph does not appear linear, with a slope of one, an operator may be alerted that there is a problem to be fixed. The operator may accordingly fix and recalibrate the instrument before measurements are taken. Once this field calibration is accepted and data are taken, for each vehicle reading, a new CO2B is calculated from the CO2 to Reference voltage ratio. In this regard, processor 110 may properly use the absorption data to derive the ΔCO2 for the exhaust plume seen. A similar process may be used for CO.
Other embodiments, uses and advantages of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. The specification should be considered exemplary only, and the scope of the invention is accordingly intended to be limited only by the following claims.

## Claims

1. A method of generating calibration curves for implementation in an optical absorption-based multi-gas analyzer (OABMGA), comprising:
injecting, into an optical path between an electromagnetic radiation source and a detector, a time-variable puff of gas from a container, wherein the container comprises a mixture of gases at known concentration ratios, and wherein the mixture of gases includes a first gas and other gases;
determining a concentration of the first gas in the puff of gas in the optical path based on a measured absorption of electromagnetic radiation, from the electromagnetic radiation source, due to the presence of the first gas in the optical path and a predetermined calibration curve for the first gas;
determining a concentration of one or more of the other gases in the optical path based on the determined concentration of the first gas in the optical path and the known concentration ratios for the mixture of gases; and
generating calibration curves for one or more of the other gases based on the determined concentrations of the one or more other gases.

2. The method of claim 1, wherein the first gas comprises propane, which is used as a surrogate hydrocarbon.

3. The method of claim 1, wherein the predetermined calibration curve for the first gas is obtained by inserting, into the optical path, a gas cell having a known concentration of the first gas.

4. The method of claim 1, wherein the other gases in the mixture of gases include at least CO and CO2.

5. The method of claim 1, wherein the electromagnetic radiation source, detector, and the optical path there-between are isolated from ambient air via an enclosure, and wherein the time-variable puff of gas is injected into the optical path within the enclosure.

6. The method of claim 5, wherein the enclosure is first purged with a purge gas for a pre-determined time period to purge any gases present within the optical path prior to injecting the time-variable puff of gas into the optical path.

7. The method of claim 6, wherein the purge gas comprises N2.

8. A system for generating calibration curves for implementation in an optical absorption-based multi-gas analyzer (OABMGA), comprising:
an optical path extending between an electromagnetic radiation source and a detector;
a container that injects a time-variable puff of gas into the optical path, wherein the container comprises a mixture of gases at known concentration ratios, and wherein the mixture of gases includes a first gas and other gases; and
a processor that is configured to:
determine a concentration of the first gas in the puff of gas in the optical path based on a measured absorption of electromagnetic radiation, from the electromagnetic radiation source, due to the presence of the first gas in the optical path and a predetermined calibration curve for the first gas;
determine a concentration of one or more of the other gases in the optical path based on the determined concentration of the first gas in the optical path and the known concentration ratios for the mixture of gases; and
generate calibration curves for one or more of the other gases based on the determined concentrations of the one or more other gases.

9. The system of claim 8, wherein the first gas comprises propane, which is used as a surrogate hydrocarbon.

10. The system of claim 8, wherein the predetermined calibration curve for the first gas is obtained by inserting, into the optical path, a gas cell having a known concentration of the first gas.

11. The system of claim 8, wherein the other gases in the mixture of gases include at least CO and CO2.

12. The system of claim 8, wherein the electromagnetic radiation source, detector, and the optical path there-between are isolated from ambient air via an enclosure, and wherein the time-variable puff of gas is injected into the optical path within the enclosure.

13. The system of claim 12, wherein the enclosure is first purged with a purge gas for a pre-determined time period to purge any gases present within the optical path prior to injecting the time-variable puff of gas into the optical path.

14. The system of claim 13, wherein the purge gas comprises N2.