Aimed at monitoring emission of organic gases such as CH4, C2H6, C3H8, iso-C4H10, n-C4H10, C2H4, C3H6, C2H2, CO, and CO2, from coal mines, petroleum refineries, and other plants, a Fourier Transform Infrared (FT-IR) spectrometer was used to develop a portable gas analyzer for patrolling and examining gas exhaust. Firstly, structure of the instrument was introduced. Then, a spectral analysis approach was presented. Finally, instrument was tested with standard gases and with actual gases emitted from a petroleum refinery. For the latter test, a gas chromatograph (GC) was used as a reference instrument. The test results showed that the detection limit of every component of analyte was less than 10 × 10−6. The maximum test error of every analyte was less than 15 × 10−6 when its practical concentration was no more than 500 × 10−6. A final comparison showed that the result curves of analytes obtained with FT-IR spectrometer almost overlapped with those obtained with GC, and their resulting noise was less than 6.4% when the practical gas concentration was above 100 × 10−6. As a result, our instrument was suitable to be used as a portable instrument for monitoring exhaust gases.
1. Introduction
In the past decades, atmospheric pollution in urban areas, especially in many developing countries around the world [1], is a major concern for environmental protection. In recent years, greenhouse gases, including CO2, N2O, NO, NO2, CH4, NH3, and even CO, are of increasing concern because of their potential impact on climate change. As a result, some standards have recently been implemented [2]. In fact, besides these component gases, there may be other components in the exhaust emitted from petroleum refineries, coal mines, and other plants and even in those volatilized from contaminated soil [3–5]. For instance, in the gas mixture emitted from coal mines, there are more than ten components. C2H6, C3H8, iso-C4H10, n-C4H10, C2H4, C3H6, and C2H2 also exist besides CH4, CO2, and CO [5]. These gases not only bring pollution to the atmosphere but also may be ignited and cause catastrophic disasters. In developing countries with power resources such as coal, petroleum, and so on, there are always many factories that emit such waste gases secretly to reduce production costs. In order to prevent such gases from being emitted continuously, governments of some countries patrol these plants. So it is necessary to develop portable instruments to analyze these areas quantitatively for a patroller.
Although there are many types of sensors [6] and instruments used for gas detection or gas mixture analysis, infrared spectroscopy, mass spectroscopy (MS), and gas chromatography (GC) are still the most common techniques used to analyze multicomponent gas mixtures with high accuracy [7–10]. While both MS and GC provide only ex situ information of analyses, infrared spectroscopy can provide in situ information [8]. Additionally, there are several potential benefits to analyzing gas mixtures with optical analyzers, including faster analysis and update rate, no required carriers gas, and maintenance-free operation. Additionally, some types of FT-IR spectrometers, such as Spectrum Two made by Perkin Elmer and Alpha made by Bruker, are portable. This makes it possible to develop portable-type gas analyzers for operators in environmental or safety administrations to patrol and ensure compliance to emissions regulations.
In fact, we have even tried to analyze the above ten components of gas with the FT-IR spectrometer in the laboratory [11]. The testing results showed that the detection limit of all the components was less than 2 × 10−6 when the light path is 10 cm and the spectrum revolution was 1 cm−1. It was not, however, suitable for patrolling gas emissions because the time to complete one analysis was excessive. For many patrollers, the expectation for a single analysis is less than 15 seconds. Therefore, some additional measures must be taken to meet this requirement.
In this paper, we aimed to develop a fast and portable-type multifunctional gases analyzer on the basis of our previous works using FT-IR spectrometer to monitor the gas mixture emitted from coal mines, oil refineries, and other plants. For the monitored gas mixture, ten components were taken as analytes: CH4, C2H6, C3H8, iso-C4H10, n-C4H10, C2H4, C3H6, C2H2, CO, and CO2. The structure of the analyzer and techniques used for the gas mixture analysis were introduced at first. Next, both a portable gas analysis instrument developed with an FT-IR spectrometer and GC were used to monitor the gas mixtures emitted from an oil refinery at the same time to illustrate the performance of our instrument. The monitoring results showed that the result curves of FT-IR spectrometer almost overlapped with those of GC. Because GC is accepted as a standard instrument for gas analysis, it can be shown that the monitor results of the FT-IR spectrometer were accurate through the overlapping result curves.
2. Experiment2.1. Instrument Structure
The schematic diagram of the portable gas analyzer based on FT-IR spectrometer is shown in Figure 1. And its physical photo is shown in Figure 2. The gross weight of the instrument was less than 25 kg. It was made with a storage cell, micro-personal computer (MPC), FT-IR spectrometer, temperature sensor, pressure sensor, and gas pump. The storage cell provided power for the FT-IR spectrometer, gas pump, and MPC. The gas pump was used to pump the analyzed gas into gas cell. The temperature sensor and pressure sensor were used to measure the gas temperature and pressure in gas cell for compensating the effects brought to the gas concentration when the instrument was used to analyze sample. In this way, the weight of thermostat and pressostat used in our previous work [11] can be reduced to keep the instrument as portable as possible. MPC read the absorption spectrum of analyte and performed spectrum analysis.
Schematic diagram of the portable gas analyzer based on an FT-IR spectrometer. The storage cell provides power for FT-IR spectrometer, gas pump, and micro-PC. The gas pump is used to pump the analyzed gas into gas cell. Sensors include temperature and pressure. They are used to monitor the temperature and pressure of analytes for compensating the effects of gas temperature and pressure on gas concentration. A micro-PC reads the absorption spectrum of analyte and performs the spectrum analysis.
Photo of the portable gas analyzer.
The type of FT-IR spectrometer was the Spectrum Two, a portable-type FT-IR spectrometer developed and produced by Perkin Elmer. Its detector type was DTGS. The optical path was 10 cm. A Norton-Beer medium function was chosen as the apodization function since it can provide the best fitting adherence to Beer’s law for the spectra measured at moderate resolution during the process of scanning the spectra of gas samples [12]. For the spectrum resolution, 4 cm−1 was chosen to make sure that a complete absorption spectrum could be gotten in less than 15 seconds because analysis speed was also required for patrolling gas emissions while the admissible error could be reduced to 10 × 10−6.
The type of the sensor used to measure pressure and temperature was BMP085 produced by Bocsh. The storage cell was also produced by Perkin Elmer. It was bought as an accessory of Spectrum Two.
2.2. Instrument Calibration
The calibration process was the same as that performed in our previous work [11]. For this, 200 sets of samples were used to calibrate the instrument. Among these samples, 30 sets were standard gases, and the others were prepared with standard gases using a gas blending system. The concentration range of every component of analyte was shown in Table 1. Differing from our previous work [11], the range of CH4 increased from 20% to 100%, that of C2H6 from 1% to 2%, and that of C3H8 from 0.5% to 1% to meet the need of monitoring natural gas emission. The temperature and pressure of samples were kept as 20°C and 760 Torr, respectively, when the calibration sample spectrum was scanned. The gas-flow rate of the gas mixture flowing through the gas cell was kept as 1 L/min. The absorption spectra of ten components of analyte with concentration of 1000 ppm were shown in Figure 3.
Concentration range of every component of analyte (percent by volume).
Item
Gas
Range
1
CH4
0~100%
2
C2H6
0~2%
3
C3H8
0~1%
4
iso-C4H10
0~0.5%
5
n-C4H10
0~0.5%
6
C2H4
0~0.5%
7
C3H6
0~0.3%
8
C2H2
0~0.3%
9
CO
0~0.3%
10
CO2
0~20%
Absorption spectra of analytes with concentration of 1000 ppm.
After the spectra of samples had been obtained, their baselines were corrected at first using the method called as The Piecewise Two-Point Autolinear Correlated Correction Method presented by us [13]. Next, the feature variable of every component of analyte was extracted with Tikhonov and Forward Selection according to its concentration and absorptivity [14]. As shown in Figure 3, for the components such as C2H2, C2H4, C3H6, CO, and CO2, whose absorption peaks have good selectivity, Forward Selection was used. For the others, Tikhonov was used because their absorption peaks overlap with others extensively. For CH4, when its concentration range was large, some absorbances of its spectrum were high and close to saturation and more than two feature variables were extracted. One was used for low concentration; the other was used for high concentration.
For the components whose feature variables have good linearity with gas concentration, a multi-input linear analysis model was used. For the components whose feature variables have good selectivity, a polynomial model was used. For the others, Artificial Neural Network (ANN) was used since it has the potential to obtain high accuracy compared with Partial Least Square (PLS) [15]. The calibration approaches were the same as used in our previous work [11].
For quantitative analysis of gas mixture, interferents must be taken into account. In oil gas emitted from oil refinery, there are water vapor and hydrogen, no BTX (benzene, toluene, and xylene), except analytes listed above in the analyzed gas mixture. For water vapor, as discussed in our previous work [11], it was easy to avoid the interference brought from it to every component of the analyte because there was at least one absorption peak that does not overlap with the absorption regions of vapor. Hydrogen is nonpolar molecule and has no absorption peak in spectrum. So it did not bring effect to the analysis of analyte.
2.3. Testing Experiment
In order to illustrate its performance, our instrument was tested with standard gas and then was tested in situ. During the process of testing, differing from the calibration experiment, temperature and pressure of the analyzed gas mixture were not kept as constant but measured and used to compensate for the effects from the change of temperature and pressure to the analysis results. At first, the absorption spectrum was scanned with the FT-IR spectrometer. Then the baseline was recognized and corrected. If both the variances of ten absorbances near 3500 cm−1 and that near 2500 cm−1 were less than 0.01, the baseline was corrected according to [13]. Otherwise, the gas cell was cleaned with nitrogen and the background spectrum was scanned again to avoid the effect of environment parameters such as temperature to the absorption spectrum. After the corrected spectrum was analyzed with the analysis models calibrated in Section 2.2, the analysis results were compensated as(1)cci=ci×760×T+273293×P,where P is the gas pressure in torr; T is the gas temperature in °C; ci denotes the concentration of the ith component of analyte obtained directly with the analysis models, while cci is the compensated result of ci. The units of ci and cci are percent by volume.
For the test in situ, GC analyzer was used as a reference instrument to analyze the same gas mixture analyzed with FT-IR spectrometer ex situ. The diameter of the silica capillary column used in the GC is 2 mm, and its length was 4 m. Nitrogen was taken as the carrier gas. The flow rate was 30 mL/min. Both the injector and detector temperature were 60°C. The FT-IR spectrometer was placed in a petroleum refinery. Next, the GC analyzer was placed outside of the petroleum refinery for the purpose of safety because it had a hydrogen flame ionization detector. When the concentration of organic gas around the GC analyzer was high, organic gas may be ignited by hydrogen flame; ignition was possible. Normally, the concentration of organic gas in petroleum refineries is also low. In order to test the performance of our instrument over larger ranges, a gas leakage was caused purposely. The leakage gas was pumped into the gas cell of the FT-IR spectrometer and then was imported into GC analyzer with a connecting pipe.
3. Results and Analysis3.1. Results Tested with Standard Gas3.1.1. Detection Limits
Detection limit is an important index for patrol instruments. In many cases of atmospheric environmental monitoring, the concentration of analyte is low, and every absorbance of the absorption spectrum of the analyzed gas mixture is close to zero. Because every absorbance at different wavenumbers has almost constant noise levels when its value is close to zero, in this work, the detection limit was determined directly with the background spectra. First, nitrogen was taken as the background gas and was scanned continuously for 30 repetitions. Then, these ground spectra were analyzed with the calibration models developed in Section 2. Finally, both the standard deviation (root-mean-square error, RMS error) and maximum error of every component were calculated and listed in Table 2. From this table, it could be found that the maximum error of every component was less than 10 × 10−6. Because the permitted concentration limit of CO for long time-work is 24 × 10−6 and that of total alkane is 80 mg/m3 (around 31 × 10−6 for C4H10) in China [16], it was acceptable if the maximum testing error list in Table 2 was taken as detection limit of the portable gas analyzer. By comparing Table 2 with our previous work [11] where the maximum analysis error in the case of background spectra was less than 2 × 10−6, it could be found that the detection limit has been raised from 2 × 10−6 to 10 × 10−6 because of the decreasing spectrum resolution of scanning times.
Mean square root of the analyte and maximum analysis error in case of background spectra.
Item
Gas
Mean square root
Maximum error
1
CH4
3.3×10-6
5.9×10-6
2
C2H6
5.4×10-6
7.7×10-6
3
C3H8
6.5×10-6
8.7×10-6
4
iso-C4H10
3.2×10-6
6.9×10-6
5
n-C4H10
6.8×10-6
9.2×10-6
6
C2H4
2.6×10-6
6.1×10-6
7
C3H6
4.0×10-6
8.1×10-6
8
C2H2
1.2×10-6
2.5×10-6
9
CO
4.6×10-6
7.9×10-6
10
CO2
2.9×10-6
6.0×10-6
3.1.2. Results Tested with Gas Mixture
The portable gas analyzer developed in this work was also tested with the gas mixture. The test samples were the same as those used in our previous work [11]. The practical concentrations of every component of the analyte covered its range. The analyzed results of the fifteen testing samples were shown in Table 3. In Table 3, Exp means the practical concentrations of every component of the gas while Anal means the analyzed results. The maximum testing error or maximum relative error was listed in Table 4 as gas concentration rise from 0 to 100%.
Testing results with standard gas mixture.
Item
Concentration of analyte/φ(⋅) × 10−6
CH4
C2H6
C3H8
iso-C4H10
n-C4H10
C2H4
C3H6
C2H2
CO
CO2
1
Exp
976
19.8
0
0
0
0
0
0
0
760
Anal
989
18.7
5.2
4.5
2.2
0
1.3
0.9
1.8
786
2
Exp
99972
998
199
0
0
0
0
0
0
0
Anal
98823
1021
212
2.1
9.3
1.8
3.4
1.1
3.6
2.1
3
Exp
99819
498
199
50.0
50.0
50.0
50.0
50.0
10.0
230
Anal
101300
475
217
44.3
56.6
44.6
55.3
48.8
9.8
221
4
Exp
19988
498
199
100
100
498
200
50.4
15.0
1010
Anal
20180
483
218
94.6
80.1
465
214
46.3
15.3
1034
5
Exp
199790
9996
2987
2993
2497
2989
998
999
600
10020
Anal
201420
10251
3129
3043
2613
2823
1145
1002
623
10132
6
Exp
10001
497
100
50.0
30.0
0
0
0
20
490
Anal
9962
473
91.6
54.6
22.6
3.9
3.1
2.1
23
504
7
Exp
10011
501
100
50.0
100
19.8
9.7
0
200
2520
Anal
10303
474
83.5
54.9
117
17.7
11.3
1.8
187
2612
8
Exp
99892
2021
1001
500
602
501
200
100
300
2300
Anal
100321
1963
963
492
591
513
213
96
315
2387
9
Exp
501020
20012
4987
4993
5017
4989
2998
2998
5002
200020
Anal
512120
19322
4845
5076
5094
4941
3082
3049
4971
203512
10
Exp
45004
2000
1000
500
400
20.1
10.1
10.2
10.2
2200
Anal
44841
2035
1072
484
431
21.6
12.3
9.67
7.1
2215
11
Exp
198972
978
255
302
200
49.8
40.0
19.8
199
700
Anal
196742
965
273
316
223
53.6
41.3
20.5
221
694
12
Exp
198
18.1
15.0
10.0
0
0
0
0
0
504
Anal
212
15.6
18.6
16.1
5.9
2.3
4.8
0.0
2.5
513
13
Exp
502
10.1
5.0
0
0
0
0
0
0
20.2
Anal
512
7.1
0.2
0.1
5.2
0.5
4.1
3.1
0
18.3
14
Exp
39995
1987
100
40.0
30.0
20.0
19.9
0
400
998
Anal
40534
1836
92.4
35.2
24.8
16.3
18.6
0.1
364
973
15
Exp
1989
98.2
19.0
10.0
0
10.0
0
5
0
0
Anal
1942
95.2
15.5
10.8
3.2
12.8
0.6
3.5
6.1
0
Maximum testing error of relative error.
Item
Concentration range
Maximum error or relative error
1
≤10 × 10−6
9.3 × 10−6
2
10~100 × 10−6
19.9 × 10−6
3
100~1000 × 10−6
9.5%
4
1000~10000 × 10−6
4.7%
5
10000~100000 × 10−6
3.5%
6
≥100000
2.2%
From Table 4, it could be shown that the relative error decreased from 9.5% to 2.2% when gas concentration rose from 100 × 10−6 to 100%. When the gas concentration was no more than 10 × 10−6, the maximum error was 9.3 × 10−6. It was the test error of n-C4H10 in the second test sample. Compared with Table 2, it was found that it was a little greater than the maximum error of n-C4H10, which is 9.2 × 10−6. This was due to the effect of methane and ethane because of extensively overlapped absorption spectra. Certainly, the error was still less than detection limit.
From Table 3, it could also be shown that the analyzed results of C2H4, C3H6, C2H2, CO2, and CO had relatively higher accuracy. For instance, among these five components when their practical concentrations were less than 10 × 10−6, C3H6 had a maximum error of 4.8 × 10−6, which appeared in the 12th sample. As discussed in our previous work [11], this was due to their absorption peaks having good selectivity.
3.2. Test Results In Situ
In order to validate the performance of our gas analyzer, a test in situ and online was also performed together with GC. Both the test results of the FT-IR spectrometer and GC were shown in Figure 4. In this figure, curves denoted as “FT-IR” mean it was calculated with spectra from the FT-IR spectrometer. Curves denoted as “GC” mean they were derived with GC. The test time was chosen from 5 A.M. to 9 A.M. because the ambient temperature changes greatly during this period of time. Because the length of the pipe linking FT-IR spectrometer and GC analyzer was greater than 30 meters, the change of gas concentration at the location of the GC analyzer lagged behind that of FT-IR spectrometer. For the simplicity of comparison, the result curves of GC analyzer had been moved to left side to make the peaks of concentrations appear at the same time. Additionally, to save space, testing results of C2H4, C3H6, C2H2, CO2, and CO were not shown here since they had relatively higher accuracy and the others were enough to demonstrate the performance of our instrument.
Analysis results obtained with FT-IR spectrometer and with gas chromatograph. (a) Methane, (b) ethane, (c) propane, (d) iso-butane, and (e) n-butane.
From Figure 4, it was easy to find that, for alkane gas monitoring, FT-IR almost overlaps with GC. This indicated that the development scheme presented in this paper was suitable. Both the calibration approach and the compensation method of temperature and pressure were included. Consequently, the portable gas analyzer could be adopted for use.
Certainly, there were also some difference between the two systems. First, the difference between FT-IR and GC was high at times. For instance, in Figure 4(e), FT-IR was 490 × 10−6 when GC was 380 × 10−6, and the maximum testing difference of n-butane reached 110 × 10−6. The relative error was 29%, much higher than 9.5%, the maximum test relative error with gas samples when the concentration range was from 100 × 10−6 to 1000 × 10−6. This might be due to interferents such as n-pentane in the testing in situ. Because of the close molecular structure between n-butane and n-pentane, the absorption spectrum of n-pentane overlaps with that of n-butane more extensively. Thus the concentration of n-pentane was added to that of n-butane. If n-pentane was taken into account as an interferent and its samples were also prepared, the difference might be reduced greatly. Certainly, in many cases, the total concentration of alkane and the first three components were accounted for, and the difference of the n-butane concentration was not a factor.
Secondly, there was noise in FT-IR, especially for n-butane. It could be evaluated with the adjacent peak-peak value and analysis result:(2a)Enoise=maxci-ci+12×meancii=1,2,…,(2b)meanci=ci-3+ci-2+⋯+ci+37⋯,where ci is the concentration of one component at ith measurement, Enoise denotes the range of noise, and max(·) means maximization operator. By searching all the results of five components of alkane, it was found that the maximum of Enoise was 6.4%. It appeared in the test results of n-butane at 6:43 A.M. At that time, GC was 280 × 10−6 while FT-IR was 290 × 10−6. The stabilization of the FT-IR spectrometer was very high. What should be noted was that the maximum errors were chosen at the times when the concentration changes slowly to avoid the effect of the fast change of gas concentration.
Third of all, when the GC were all zeroes, for example, as shown in Figure 4(d), in the interval of 7:25 A.M. to 7:32 A.M., FT-IR were not all zeroes but decreasing slowly. This might be due to the gases being diffused into the gaps between sample cell and windows. In almost all FT-IR spectrometers, there is a cavity for placing the sample cell. In order to keep all the elements of Michelson interferometer clean, at two sides of the cavity, there are two windows for keeping the Michelson interferometer isolated from atmosphere. Between the two windows, there is a supporting stage for placing the sample cell. There are two little gaps between the sample cell and the windows. As mentioned in Section 2.3, the FT-IR spectrometer was just placed in petroleum refinery while GC was placed 50 meters away. Because the analyzed gas diffused into the gaps between the gas cell and the spectrum windows, and the cavity was badly ventilated, the concentrations of the analytes in the gaps changed slowly. Even after the concentrations of all the compositions in the gas cell had decreased to zero for several minutes, the concentrations of the gases in the gaps might be still above zero or even at high levels.
4. Conclusions
In this paper, both the structure and parameters of the portable gas analyzer based on the FT-IR spectrometer were presented at first. Next, the analysis approach for multicomponent gas mixture was introduced. Finally, the instrument developed in this work was tested with standard gas and tested in situ together with GC. The test results showed that the portable instrument developed with FT-IR spectrometer had the following performance:
This instrument was suitable for patrolling exhaust gas and for monitoring multicomponent gas mixtures including CH4, C2H6, C3H8, iso-C4H10, n-C4H10, C2H4, C3H6, C2H2, CO2, and CO;
The detection limit of every component of analyte was less than 10 × 10−6;
The maximum noise level was 6.4% when the practical gas concentration was above 100 × 10−6.
Finally, it could be concluded that the development scheme for the portable gas analyzer was practical. The portable gas analyzer developed in this work could be adopted for use in patrolling and monitoring exhaust gas. Besides oil refineries, there are many other applications. In China, both coalbed methane and coal mine gas [16] are included at least.
Conflict of Interests
The authors declare no conflict of interests.
Acknowledgments
The authors are grateful for financial support by National Great Program on Science Instrument Development of China (2012YQ240127) and by National Natural Science Foundation of China (51277144).
BegumB. A.BiswasS. K.NasiruddinM.Showkot HossainA. M.HopkeP. K.Source identification of chittagong aerosol by receptor modeling200926367968910.1089/ees.2008.00552-s2.0-61949432261National Institute of Standards and TechnologyUnited States Department of Commerce2009CarrasJ. N.DayS. J.SaghafiA.WilliamsD. J.Greenhouse gas emissions from low-temperature oxidation and spontaneous combustion at open-cut coal mines in Australia200978216116810.1016/j.coal.2008.12.0012-s2.0-61849154040Dominguez-RosadoE.PichtelJ.Phytoremediation of soil contaminated with used motor oil: II. Greenhouse studies2004212169180AdamusA.ŠancerJ.GuřanováP.ZubičekV.An investigation of the factors associated with interpretation of mine atmosphere for spontaneous combustion in coal mines201192366367010.1016/j.fuproc.2010.11.0252-s2.0-78751649480PejcicB.EadingtonP.RossA.Environmental monitoring of hydrocarbons: a chemical sensor perspective200741186333634210.1021/es07045352-s2.0-34548779091FeiS.-W.ZhangX.-B.Fault diagnosis of power transformer based on support vector machine with genetic algorithm2009368113521135710.1016/j.eswa.2009.03.0222-s2.0-67349288151UsseglioS.ThorshaugK.KarlssonA.DahlI. M.NielsenC. J.JensK.-J.TangstadE.In situ infrared emission spectroscopy for quantitative gas-phase measurement under high temperature reaction conditions: an analytical method for methane by means of an innovative small-volume flowing cell201064214114810.1366/0003702107906195642-s2.0-77649223980EslerM. B.GriffithD. W. T.WilsonS. R.SteeleL. P.Precision trace gas analysis by FT-IR spectroscopy. 1. Simultaneous analysis of CO2, CH4, N2O, and CO in air200072120621510.1021/ac99056252-s2.0-0033627048ThompsonC. J.RileyR. G.AmonetteJ. E.GassmanP. L.Quantification of volatile organics in soil aging experiments using fourier transform infrared spectroscopy200660891491910.1366/0003702067780621562-s2.0-33748313015TangX.LiangY.DongH.SunY.LuoH.Analysis of index gases of coal spontaneous combustion using fourier transform infrared spectrometer20142014841439110.1155/2014/4143912-s2.0-84911091456GriffithsP. R.De HasethJ. A.20072ndHoboken, NJ, USAJohn Wiley & SonsZhaoA.TangX.LiW.ZhangZ.LiuJ.The piecewise two points autolinear correlated correction method for fourier transform infrared baseline wander201548427427910.1080/00387010.2013.8745302-s2.0-84911940799StoutF.KalivasJ. H.HébergerK.Wavelength selection for multivariate calibration using Tikhonov regularization2007611859510.1366/0003702077797014792-s2.0-33846654921YangH.GriffithsP. R.TateJ. D.Comparison of partial least squares regression and multi-layer neural networks for quantification of nonlinear systems and application to gas phase Fourier transform infrared spectra2003489212513610.1016/s0003-2670(03)00726-82-s2.0-10644272571Emission standard of coalbed methane/coal mine gasChina, 2008