1. Field of the Invention
The present invention relates to an exhaust gas temperature analysis apparatus, analysis method, and analysis program.
2. Description of the Related Art
A conventional apparatus passes a laser beam through an exhaust gas and determines the concentration and temperature of specific gas components from the transmittance of the specific component (refer to, for example, Japanese Patent Application Publication No. 2004-117259 (JP-A-2004-117259)).
In particular, JP-A-2004-117259 describes processing units, such as a personal computer mounted within an automobile that analyze the concentration of gas components contained in, and the temperature of, the engine exhaust gas.
In the technology described in JP-A-2004-117259, the signal strength ratio I/Io of the intensity Io of light of a reference laser beam with the intensity I of transmitted light of a measurement laser beam is measured. When calculating the concentration of a gas component based on this signal strength ratio, it is necessary to determine the absorbance of the gas component. In calculating the gas component concentration based on the signal strength ratio, therefore, it is first necessary to analyze the temperature of the gas component.
In conventional temperature analysis of a gas component, for example as shown in
In the conventional method of fitting, from the theoretical spectra R1, R2, and so on, the absorption amount of various theoretical spectra is integrated to calculate and determine the theoretical spectrum that is closest to the measured spectrum M, and, in some cases, a determination is made using a calculation by taking into consideration the level of coincidence of the peak wavelengths.
With regard to the calculation, if it is desired to determine the temperature to the nearest degree Kelvin within the range of 1000K (Kelvin), for example, it is necessary to perform fitting among 1000 different spectra to determine the one that is closest to the measured spectrum M.
The time required to perform these calculations presents an obstacle to performing the gas component analysis quickly and accurately, and there is a need to reduce the calculation time. Also, when performing calculations with a high resolution using conventional technology, for example, in the case of determining the temperature in units of 0.1K, because the number of theoretical spectra for fitting increases even more, it becomes an even more serious problem to reduce the calculation time.
The present invention provides an approach for reducing the calculation time for calculating a temperature by fitting a measured spectrum to a theoretical spectrum
A first aspect of the present invention relates to a method for analyzing the temperature of the exhaust gas. The method fits a measured spectrum of a measured gas to theoretical spectra predefined for a series of temperatures and pressures wherein the temperature of the measured gas is defined as that of the theoretical spectrum that most closely fits the measured spectrum, and calculates an approximate temperature at least one time by fitting the measured spectrum to a portion of theoretical spectra defined in association with a first temperature range using a temperature calculated in an immediately preceding temperature analysis as a reference, and determining a calculated temperature by fitting the measured spectrum to all the theoretical spectra over a second temperature range that is narrower than the first temperature range, using the approximate temperature as a reference.
According to the first aspect of the present invention, if there is a difference between a pressure of the theoretical spectra referred to in the fitting and a pressure of the measured gas, the method may perform a discrete-calculation correction on the calculated temperature with respect to the pressure difference to calculate a corrected temperature.
In the first aspect, the temperature calculated in the immediately preceding temperature analysis may be one of an approximate temperature determined in an immediately preceding temperature analysis, a calculated temperature determined in an immediately preceding temperature analysis, or a corrected temperature determined in an immediately preceding temperature analysis.
A second aspect of the present invention is an exhaust gas temperature analysis apparatus that fits the measured spectrum of the measured gas to all theoretical spectra defined for a series of temperatures and pressures, wherein the temperature of the measured gas is defined as that of the theoretical spectrum that most closely fits the measured spectrum. The apparatus of the second aspect includes: a first fitting device that at least one time fits the measured spectrum to a portion of the theoretical spectra defined in association with a first temperature range having a temperature determined in an immediately preceding time as a reference; and a second fitting device that fits the measured spectrum to all theoretical spectra defined in association with a second temperature range that is narrower than the first temperature range to determine a calculated temperature.
The second aspect of the present invention may further include a discrete-calculation correction device that performs a discrete-calculation correction on the calculated temperature with respect to the pressure difference to determine the corrected temperature when there is a pressure difference between the pressure in a theoretical spectrum referenced in the fitting and the measured pressure.
In the second aspect of the present invention, the temperature determined in the immediately preceding temperature analysis may be one of an approximate temperature determined in the immediately preceding temperature analysis, a calculated temperature determined in the immediately preceding temperature analysis, or a corrected temperature determined in the immediately preceding temperature analysis.
A third aspect of the present invention is an exhaust gas temperature analysis program that may be executed by a computer to fit a measured spectrum of a measured gas to theoretical spectra that are pre-defined for a series of temperatures and pressures, wherein the temperature of the measured gas is defined as that of the theoretical spectrum that most closely fits the measured spectrums. The program has instructions for calculating an approximate temperature at least once by fitting the measured spectrum to a portion of the theoretical spectra defined in association with a first temperature range using a temperature determined in the immediately preceding temperature analysis as a reference, and determining a calculated temperature by fitting the measured spectrum to all the theoretical spectra defined in association with a second temperature range that is narrower than the first temperature range, using the approximate temperature as a reference.
A program of the third aspect may include instructions for performing a discrete-calculation correction on the calculated temperature with respect to a pressure difference when there is a difference between the pressure in a theoretical spectrum referenced in the fitting and the measured pressure of the measured gas to determine a corrected temperature.
In the program of the third aspect, the temperature determined in the immediately preceding temperature analysis may be one of an approximate temperature determined in an immediately preceding temperature analysis, a calculated temperature determined in an immediately preceding temperature analysis, or a corrected temperature determined in an immediately preceding temperature analysis.
According to the first, second, and third aspects of the present invention, it is possible to reduce the number of times to fit the measured spectrum to the theoretical spectra and to more quickly perform gas component analysis by reducing the calculation time, that is, by enabling real-time temperature analysis. Also, by enabling the collection of a large amount of data, it is possible to perform analysis with higher accuracy.
According to the above-described aspects of the present invention, by performing discrete-calculation correction, it is possible to eliminate the influence of pressure difference. Thus, data may be obtained with higher reliability.
According to the above-described aspects of the present invention, when it is desired to obtain information regarding the approximate temperature in a short period of time (approximate temperature) and when it is desired to obtain information regarding the temperature with relative good accuracy in a short period of time (calculated temperature), and when there is a reserve of processing capacity available and it is desired to reliably perform analysis with high accuracy (corrected temperature), the temperature that is used in the subsequent temperature analysis may be selected in accordance with the particular situation.
The foregoing and further features and advantages of the invention will become apparent from the following description of example embodiments with reference to the accompanying drawings, wherein like numerals are used to represent like elements, and wherein:
An example embodiment of the present invention is described with respect to the drawings.
As shown in
As shown in
The wavelength of the infrared laser passed through the exhaust gas in the sensor base 6 is set in accordance with the exhaust gas component to be detected. Thus, to carbon monoxide (CO), carbon dioxide (CO2), ammonia (NH3), methane (CH4), and water (H2O), infrared lasers of five different wavelengths are used. For example, 1530 nm is a suitable wavelength for detecting ammonia, 1560 nm is a suitable wavelength for detecting carbon monoxide, 1570 nm is a suitable wavelength for detecting carbon dioxide, 1680 nm is a suitable wavelength for detecting methane, and 1350 nm is a suitable wavelength for detecting water. Additionally, in order to detect other exhaust gas components, infrared lasers of additional wavelengths may need to be used. The number of laser diodes required corresponds to the number of exhaust gas components to be analyzed.
The infrared light emitted from the laser diodes LD1 to LD5 passes through the optical fibers 32 and is input to the light splitters 33, which divide the light into the measurement laser I and the reference laser Io. The measurement laser I passes through the optical fiber 34A and the light mixer 35, and then through the light-transmitting fiber 7 that guides the light to the sensor base 6. Measurement laser I is sequentially emitted from the light mixer 35 at prescribed intervals. The number of the laser diodes LD1 to LD5 provided corresponds to the number of gas components to be analyzed. For example, in the case of analyzing five gas components, five laser diodes would be provided, and in the case of analyzing ten gas components, ten laser diodes would be provided. The reference laser Io passes through the optical fiber 34B and is guided to the light mixer 36.
Laser light that is sequentially emitted from the light-transmitting optical fiber 7 and attenuated through the exhaust gas passes through the light-receiving optical fiber 8 provided in the sensor base 6 and is received by the photodiode PD1. The output of the photodiode PD1 is amplified, for example, by a pre-amplifier (not shown), passes through an A/D converter, and is input to the processing unit 20. The reference laser Io that is input to the light mixer 36 passes through an optical fiber 39 and is directly received by the photodiode PD2, the output of the photodiode PD2 being input to the processing unit 20.
The processing unit 20 synchronizes the sequentially emitted measurement laser I and the reference laser Io to correspond to the laser diodes LD1 to LD5, that is, to correspond to the gas components to be analyzed, the absorption spectrum (measured spectrum) of each gas component being measured, and the signal strength ratio (I/Io) between the intensity of the reference laser Io and the intensity of the transmitted measurement laser I also being measured. The processing unit 20 also fits the measured spectrum to theoretical spectra.
By doing the above, the signal strength ratio of each gas component is calculated, and the temperature of the exhaust gas at the time of the calculation is calculated so as to calculate the concentration of each gas component. This embodiment is directed to the calculation of the temperature of the exhaust gas, as such the description regarding the calculation of the concentration being omitted.
Next, the calculation of the temperature of the exhaust gas will be described. The exhaust gas temperature calculation may be performed separately for each gas component contained in the exhaust gas. For example, the temperature of a gas component such as H2O, which has an absorption spectrum with a prominent peak, is calculated and may be used as a representative temperature of the exhaust gas as a whole, thereby obtaining temperature values having high reliability.
In this embodiment, as shown in
If there is a pressure difference Δp between the pressure of the theoretical spectrum referenced in the above-noted fitting (matching) and the measured pressure P, a discrete-calculation correction is performed in step S5 on the calculated temperature Tβ with respect to the pressure difference Δp to determine a corrected temperature Tγ (fourth temperature Td).
The temperature determined in the immediately preceding temperature analysis may be the approximate temperature Tα determined in the immediately preceding temperature analysis, the calculated temperature Tβ determined in the immediately preceding temperature analysis, or the corrected temperature Tγ determined in the immediately preceding temperature analysis.
More specifically, the exhaust gas temperature analysis method used is that of fitting (matching) the measured spectrum to theoretical spectra for each temperature TA to TB and each pressure P1, P2, and so on, and setting the temperature defined in association with the theoretical spectrum selected by the fitting (matching) as the temperature of the exhaust gas. In step S1, a first temperature Ta is selected from the overall temperature range TA to TB with which the theoretical spectra are defined in association. In step S2, it is performed to fit (match) the measured spectrum to a plurality of theoretical spectra that are associated with the first temperature range TC to TD having the first temperature Ta as a reference and that are also defined in association with a representative pressure P2 at which the pressure difference Δp with respect to the measured pressure P is minimum, to select a theoretical spectrum Rb in which the fitting deviation (matching deviation) Δs2 with respect to the measured spectrum is minimum. In step S3, it is performed to fit (match) the measured spectrum to a plurality of theoretical spectra that are associated with the second temperature range TE to TF having the second temperature Tb as a reference and that also are defined in association with a representative pressure P2 at which the pressure difference Δp with respect to the measured pressure P is minimum, to select a theoretical spectrum Rc in which the fitting deviation Δs3 with respect to the measured spectrum is minimum. In step S4, the third temperature Tc, defined in association with the theoretical spectrum Rc selected in step S3, is determined as the temperature of the measured gas.
If there is a pressure difference Δp between the measured pressure P and the representative pressure P2, a discrete-calculation correction is performed on the third temperature Tc determined in step S4 with regard to the pressure difference Δp in step S5, and the fourth temperature Td, defined by the discrete-calculation correction performed in step S5, is determined as the temperature of the measured gas in step S6.
In step S7, one of the second temperature Tb defined in association with the theoretical spectrum Rb selected in the second step S2, the third temperature Tc defined in association with the theoretical spectrum Tc selected in the third step S3, or the fourth temperature Td determined as the temperature of the measured gas in the sixth step S6, is defined as a selection candidate for the first temperature Ta in the first step S1 of the succeeding temperature analysis.
By doing the above, it is possible to reduce the number of times to fit (match) the measured spectrum to the theoretical spectra, making it possible to perform gas component analysis at high speed by reducing the calculation time. In addition, the method also facilitates more precise analysis by allowing the collection of a large amount of data. By performing discrete-calculation correction, it is possible to eliminate the influence of the pressure difference Δp, thereby improving reliability of resultant data obtained.
Each step of the method will be described in detail below. First, in step S1, as shown in
The first temperature Ta may be one of the second temperature Tb determined in step S2, the third temperature Tc determined in step S3, or the fourth temperature Td determined in step S6 in the flow of the immediately preceding temperature analysis. In this manner, in step S1 the temperatures Tb, Tc, and Td determined in the preceding temperature analysis are referenced. Also, in the first temperature analysis, because there are no temperatures Tb, Tc, and Td to reference, instead of performing a fitting within the first temperature range TC to TD in step S2, fitting is performed over the entire temperature range TA to TB.
As shown in
Next, in step S2, as shown in
In step S2, as shown in
In step S2, as shown in
Additionally, for example, if the theoretical spectrum is defined in increments of 1K and the first temperature range TC to TD has a span of 100K (the ±50K case noted above), if an attempt is made to fit the measured spectrum to all the theoretical spectra, the calculation must be performed 100 times, this number of calculations presenting a great load. Given this, for example, ten theoretical spectra are selected in 10K units starting at the lowest temperature in association with the first temperature range TC to TD, and it is performed to fit the measured spectrum to the ten selected theoretical spectra only, thereby enabling a great reduction in the number of calculations, and a shortening of the calculation time.
That is, the theoretical spectra that are the targets of fitting in step S2 are set as a portion of the theoretical spectra defined in association with the first temperature range TC to TD, for example, taken at a prescribed temperature interval (10K intervals in this example), to select theoretical spectra and reduce the number of fittings. By using a prescribed temperature interval, it is possible for the theoretical spectrum to represent the theoretical spectra corresponding to the surrounding temperatures, thereby enabling extraction with little dispersion from the first temperature range TC to TD, which the superset of temperatures, and enabling an improvement in reliability.
Next, in step S3, as shown in
In step S3, as shown in
With regard to the fitting in step S3, for example, when the theoretical spectra are defined in increments of 1K, if the span of second temperature range TE to TF is 20K (the case of ±10K described above), fitting of the measured spectrum is performed to all the theoretical spectra, and narrowing down is performed based on the minimum temperature unit (in this example, 1K) that enables fitting by theoretical spectra.
That is, the theoretical spectra that are the targets of fitting in step S3 are all the theoretical spectra that are defined in association with the second temperature range TE to TF.
Next, as shown in
That is, as shown in
The discrete-calculation correction performed in step S5, for example, as shown in
The relationship between the “pressure difference Δp” and the “temperature offset amount” is evaluated and calculated beforehand with regard to the measured gas components for the purpose of measuring the temperature, and is stored as a function N in the processing unit 20. If the order of the function N stored in this manner in the processing unit 20 is small, it is possible to reduce the load of the calculations required to perform discrete-calculation correction of the pressure. The discrete-calculation correction of the pressure is not restricted to this embodiment. For example, the temperature may be corrected by multiplying the third temperature Tc by the pressure correction coefficient, which corresponds to the “pressure difference Δp” stored in the processing unit 20.
Next, as shown in
As shown in
In the processing unit 20, the temperatures T1, T2, and T3 (Tb, Tc, and Td) may be continuously monitored in real time, and the temperatures may be referenced effectively in response to the situation.
Additionally, a previously determined second temperature Tb may be used to detect errors in the calculation process and mechanical errors. For example, if the second temperature Tb coincides with a border value of the first temperature range TC to TD, the second temperature Tb will differ greatly from the first temperature Ta, and this fact may be used to identify a discontinuous change in the temperature and perform error detection.
Also, step S2 may be executed a plurality of times. That is, in step S3, when determining the second temperature Tb used for theoretical spectrum referencing, the narrowing down of the second temperature Tb is performed by executing, in step-wise fashion, a step such as the second step S2 a plurality of times.
With regard to the above-described method, the design can be made in accordance with the required analysis time, and by changing the design it is possible to accommodate the required analysis time.
Additionally, with regard to reducing the number of fittings and the calculation time, in the case of, for example, a different embodiment, there is a need to fit the measured spectrum to all the theoretical spectra over the overall temperature range TA to TB. All regions surrounded by TA, TB, and Tc until finally reaching the third temperature Tc may be regarded as a calculation load. In contrast, according to this embodiment, the region surrounding by TC, TE, Tc, TF, and TD may be regarded as the calculation load. Thus, it is possible to reduce the calculation load with regard to the regions A1 and A2.
With regard to the specific number of fittings in the overall temperature range TA to TB of 1000K, for example, if the theoretical spectra are defined in increments of 1K, a simple calculation shows that there is a need to perform fitting 1000 times. In contrast, in this embodiment it is possible to narrow down to the third temperature Tc by performing the fitting a total of 20 times, specifically by performing fitting 10 times in the first temperature range TC to TD in step S2, and performing fitting 10 times in the second temperature range TE to TF in step S3. If it is desired to perform temperature analysis with higher accuracy, this can be done by increasing the number of theoretical spectra that are fitted, and the number of fittings can be set in response to the required accuracy.
The above-noted analysis may be implemented using an apparatus having the following constitution. Specifically, the apparatus may be an exhaust gas temperature analysis apparatus that fits a measured spectrum of a measured gas to theoretical spectra for each temperature TA to TB and each pressure P1, P2, and so on, and sets a temperature that is defined in association with a theoretical spectrum selected by performing fitting as the temperature of the measured gas. The exhaust gas temperature analysis apparatus includes a first device that fits the measured spectrum to a portion of the theoretical spectra defined in association with the first temperature range (TC to TD) at least one time, using a temperature determined in the immediately preceding time (Ta, Tb, Tc, Td) as a reference, and that determines an approximate temperature Tα (in this embodiment, the second temperature Tb), and a second device that fits the measured spectrum to all theoretical spectra defined in association with the second temperature range (TE to TF) that is narrower than the first temperature range, and that determines a calculated temperature Tβ (in this embodiment, the third temperature Tc).
The above apparatus has a discrete-calculation correction device that performs a discrete-calculation correction on the calculated temperature Tβ with regard to the pressure difference Δp if there is a pressure difference Δp between the pressure in a theoretical spectrum referenced in the fitting and the measured pressure P of the measured gas, to determine the corrected temperature Tγ.
More specifically, as shown in
If there is a pressure difference Δp between the measured pressure P and the representative pressure P2, a fifth device performs a discrete-calculation correction on the third temperature Tc with regard to the pressure difference Δp, and a sixth device sets the fourth temperature Td as the temperature of the measured gas.
There is also a seventh device that defines any one of the second temperature Tb defined in association with the theoretical spectrum Rb selected by the second device, the third temperature Tc defined in association with the theoretical spectrum Tc selected by the third device, or the fourth temperature Td set as the temperature of the measured gas by the sixth device, as a selection candidate for the first temperature Ta in the first device in the subsequent temperature analysis.
The function of the above-noted apparatus may be implemented in the processing unit 20 of this embodiment. Specifically, the above devices may be implemented in a computer that combines a CPU, memory, and a database, and also be a dedicated microprocessing unit. It is also possible to access a dedicated database via a network to reference the table 50, and diverse variations are implementable utilizing widely known art.
The above-noted arrangement may also be implemented by a program. Specifically, this may be an exhaust gas temperature analysis program that includes instructions for fitting a measured spectrum to theoretical spectra defined in association with each temperature TA to TB and each pressure P1, P2 and so on, and setting a temperature, defined in association with the theoretical spectrum selected by the fitting, as the temperature of the measured gas. The temperature analysis program uses a computer to, at least one time, fits the measured spectrum to a portion of the theoretical spectra defined in association with the first temperature range (TC to TD), using a temperature (Ta, Tb, Tc, or Td) determined in the immediately preceding analysis as a reference, and to determine an approximate temperature Tα (in this embodiment, the second temperature Tb), and to fit the measured spectrum to all theoretical spectra defined in association with the second temperature range (TE to TF) that is narrower than the first temperature range, and then to determine a calculated temperature Tβ (in this embodiment, the third temperature Tc).
The above-noted program has an exhaust gas temperature analysis block that causes the computer to function as a discrete-calculation correction device that performs a discrete-calculation correction on the calculated temperature Tβ with regard to the pressure difference Δp if there is a pressure difference Δp between the pressure in a theoretical spectrum referenced in the fitting and the measured pressure P of the measured gas, to determine the corrected temperature Tγ (the fourth temperature Td).
More specifically, as shown in
The exhaust gas temperature analysis program instructs the computer to execute the functions of a fifth device that performs a discrete-calculation correction on the third temperature Tc with regard to the pressure difference Δp if there is a pressure difference Δp between the measured pressure P and the representative pressure P2, and a sixth device that sets the fourth temperature Td as the temperature of the measured gas.
The exhaust gas temperature analysis program uses the computer as a seventh device that may define one of the second temperature Tb defined in association with the theoretical spectrum Rb, the third temperature Tc defined in association with the theoretical spectrum Tc, or the fourth temperature Td set as the temperature of the measured gas by the sixth device, as a candidate for the first temperature Ta used in the first device in a subsequent temperature analysis.
According to the above-described temperature analysis program, it is possible to implement high-speed processing, and also to implement highly accurate exhaust gas analysis.
Number | Date | Country | Kind |
---|---|---|---|
2006-340590 | Dec 2006 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB2007/003952 | 12/17/2007 | WO | 00 | 11/3/2008 |