Next Article in Journal
Plasmon-Enhanced Controlled Drug Release from Ag-PMA Capsules
Next Article in Special Issue
Nanostructured Thin Coatings Containing Anthriscus sylvestris Extract with Dual Bioactivity
Previous Article in Journal
Predicting Drug Resistance Using Deep Mutational Scanning
Previous Article in Special Issue
Advances in Mid-Infrared Spectroscopy-Based Sensing Techniques for Exhaled Breath Diagnostics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Analysis of Photoacoustic Spectra for the Detection of Short-Chained Hydrocarbon Isotopologues

1
Heinrich Blasius Institute of Physical Technologies, Hamburg University of Applied Sciences, Berliner Tor 21, 20099 Hamburg, Germany
2
School of Computing, Engineering and Physical Sciences, University of the West of Scotland, High Street, Paisley PA1 2BE, UK
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(9), 2266; https://doi.org/10.3390/molecules25092266
Submission received: 24 February 2020 / Revised: 1 May 2020 / Accepted: 7 May 2020 / Published: 11 May 2020
(This article belongs to the Special Issue Applications of Photoacoustic Spectroscopy)

Abstract

:
We report, to our knowledge, the first optical detection scheme for short-chained hydrocarbon isotopologues. The sensor system is based on photoacoustic spectroscopy (PAS). Two continuous wave, thermoelectrically cooled, distributed feedback interband cascade lasers (DFB-ICLs) with emission wavelengths around 3.33 and 3.38 μm, respectively, served as light sources. The investigations comprised the main stable carbon isotopologues of methane (12CH4, 13CH4), ethane (12CH3-12CH3, 13CH3-12CH3, 13CH3-13CH3), and propane (12CH3-12CH2-12CH3, 13CH3-12CH2-12CH3). They were selected because of their importance for numerous applications from climate and planetary research to natural gas exploration. Multiple measurements of single components in nitrogen and synthetic mixtures were conducted at room temperature and atmospheric pressure. Depending on the investigated hydrocarbon isotopologue, detection limits ranging from 0.043 ppmv to 3.4 ppmv were achieved. For a selective concentration determination, multivariate analysis (MVA) was applied. Partial least-squares regression (PLSR) was used to calculate concentrations from the PA spectra. The implementation of MVA has shown that the PA setup in principle works reliably and that the selective concentration determination of short-chained hydrocarbon isotopologues is possible.

1. Introduction

In environmental policy discussions, people often only talk about the carbon dioxide (CO2) emissions that need to be reduced and controlled in order to stop global warming. Due to simplified news coverage, many people are not aware that atmospheric methane (CH4) plays a major role in the greenhouse effect as it has a considerably higher global warming potential compared to CO2, i.e., CH4 emission will have 28 times the impact on temperature of a CO2 emission of the same mass over the following 100 years [1]. In addition to its own climate impact, the presence of CH4 in the atmosphere can also affect the concentration of other greenhouse gases, such as tropospheric ozone (O3), water vapor, and CO2 [2]. While CH4 does not cause direct harm to human health or crop production, its role as precursor gas contributes significantly to the negative health and agricultural impacts of O3 [3,4]. In contrast to CO2, the climate system reacts rapidly to changes in CH4 emissions and reducing CH4 emissions could provide the opportunity to immediately achieve major benefits for global and regional climate as well as for human health and agriculture [5]. A better understanding of the global methane budget is crucial to take meaningful measures. This can only be achieved by the sensitive measurement of short-chained (SC) hydrocarbons including their stable carbon isotopic ratios, as it can be used as a proxy for source and sink estimation (natural vs. anthropogenic), be an indicator for atmospheric processes, and therefore improve and correct atmospheric and climate models [6].
A second application involving hydrocarbon isotopologue detection concerns natural gas (NG) exploration. Natural gas is the fastest growing fossil fuel on the global energy market and consumption of NG worldwide is projected to increase more than 40% between 2018 and 2050. Thus, total consumption will reach nearly 200 quadrillion British thermal units (1 Btu = 1055 J) by 2050, according to the International Energy Outlook 2019 [7]. Asia especially will experience the largest rise in global NG consumption because the economies of big countries such as China and India are rapidly growing. Natural gas usage expedites the most to fulfil demand from increased industrial activity, NG-related electricity generation, and fuel-related transportation. The International Energy Outlook also stated that NG remains a key fuel in the electric power sector and in the industrial sector because in the power sector, NG is still an attractive choice for new generating plants due to its fuel efficiency. Natural gas burns cleaner than coal or oil products, and as more governments begin to implement national or regional plans to reduce CO2 emissions, they may promote the use of NG to replace more carbon-intensive coal and liquid fuels. Again, the sensitive measurement of the stable carbon isotopic composition of NG plays a key role in the production of this energy source. The 13C/12C ratio of SC hydrocarbons is strongly related to the gas type or source (i.e., bacterial, thermogenic, or mixed) and to its maturity [8]. A steadily growing share of NG production originates from shale. Several studies showed that high-production wells are often found in regions where a so called “isotopic reversal” in shale gas occurs [9,10,11]. The isotopic analysis in combination with conventional logging protocols can represent a powerful tool to control horizontal drilling and can significantly improve exploration, completion, and production technologies [12]. Ultimately, this will have both positive economic and ecological effects.
The standard technique for the measurement of carbon isotopes for the past 40 years has been the isotope ratio mass spectrometry (IRMS) together with the gas chromatography (GC) [13]. This method is highly precise but also labor-intensive, quite expensive, and typically performed in a facility, making real-time in-situ measurements impossible. GC-IRMS is also an indirect method and destroys the sample because the hydrocarbon must be converted into CO2 by burning prior to analyses of the isotope ratio. Further development of alternative isotope ratio measurement techniques is thus highly desirable. An alternative measurement method can be based on spectroscopy. The mid infrared (IR) range has proven to be particularly suitable for hydrocarbon analysis, since many hydrocarbon molecules exhibit strong rotational vibration transitions that are characteristic of the respective molecule. In the mid IR region, several laser sources are available, which include CO and CO2 gas lasers, lead-salt diode lasers, coherent sources based on difference frequency generation, optical parametric oscillators, quantum cascade lasers, and interband cascade lasers (ICLs) [14]. Distributed feedback (DFB) ICLs combine a good beam quality, single-mode operation, a small linewidth (<10 MHz) and continuous wave (cw) emission and mode-hop-free spectral tunability with a compact size, operation at room temperature and easy handling [15]. Together with their recent availability with wavelengths between 3 μm and 4 μm this makes them an excellent choice for the sensitive detection of SC hydrocarbon gases.
In a previous study, IR absorption cross sections of the ethane and propane isotopologues were measured for the first time in the mid IR range using an FTIR spectrometer [16]. Subsequently, high-resolution absorption cross sections of the same isotopologues were photo-acoustically measured using ICLs near 3.33 and 3.38 µm, respectively [15].
Photoacoustic spectroscopy (PAS) is based on the conversion of light to sound [17] and among spectroscopic techniques it offers many advantageous properties such as wide dynamic range, offset freeness, and simplicity [15]. It has been successfully applied in numerous trace gas sensing applications, which include atmospheric investigations (e.g., environmental monitoring, industrial emission measurements), detection of toxic gases, chemical warfare agents and explosives, biological and agricultural applications, industrial/manufacturing processes, and medical diagnostic (e.g., non-invasive breath analysis) [18,19,20].
In this work, we report the first PAS sensor for the detection of methane (CH4), ethane (C2H6), and propane (C3H8) including their main stable 13C-isotopologues (13CH4, 13CH3-12CH3, 13CH3-13CH3 and 13CH3-12CH2-12CH3) by using two DFB-ICLs operating in the wavelength range from 3327.7 nm to 3334.7 nm (ICL 1541, methane and ethane detection) and from 3373.3 nm to 3396.7 nm (ICL 1638, propane detection), respectively. In order to obtain the concentrations of individual substances of a measured spectrum containing multiple substances, multivariate analysis (MVA) [21,22] can be applied. This was also demonstrated in a recent work, where PA measurements of volatile organic compounds from strongly overlapping spectra in the mid IR region were evaluated [23]. Well-known multivariate regression models are the multiple linear regression, the principal component regression, and the partial least-square regression (PLSR), whereby the latter has become the standard method for spectroscopy [21,24]. If the data are noisy and strongly correlated, PLSR is considered to be superior to the other techniques [23,24]. Here, a PLSR model was trained and verified with experimental data of single and multiple components in synthetic gas mixtures. In addition, the sensor’s limit of detection (LOD) was determined.

2. Materials and Methods

2.1. Photoacoustic Spectroscopy

The PAS technique uses the fact that absorbed electromagnetic radiation is partially converted into kinetic energy by inelastic molecular collisions. This corresponds to an increase in the temperature of the irradiated volume. If the radiation is modulated, this will lead to modulated heating and, thus, generate a pressure wave with an amplitude proportional to the concentration of the absorbing molecules, which can be recorded with a microphone. For a single absorbing compound and under the condition that the absorption is not saturated, the wavelength-dependent PA signal S (in V) as function of the wavelength λ , is given as [15,25]
S ( λ ) = C P ( λ ) N tot c σ ( λ )
where C (in V/(cm−1 W)) is known as the cell constant depending mainly on the cell geometry but also on the cell wall surface, beam profile, the modulation frequency, and the microphone sensitivity. P ( λ ) (in W) is the laser power, N tot is the total number density of molecules (in molecule/cm3), and the coefficients c and σ ( λ ) are the concentration and absorption cross section (in cm2 molecule−1) of the substance of interest, respectively. The PA signal is directly proportional to the laser power, concentration of the absorbing molecules, and the cell constant. If the laser modulation frequency corresponds to an acoustical resonance frequency of the cell, a significant signal amplification can be obtained.

2.2. Experimental Setup

A schematic diagram of the experimental setup is illustrated in Figure 1.
Two cw DFB-ICLs operating at wavelengths around 3.33 μm (ICL 1541, Nanoplus, Gerbrunn, Germany) and 3.38 μm (ICL 1638, Nanoplus), respectively, were selected as radiation sources for the isotope-selective detection of methane/ethane and propane. The laser sources (ICLs) were employed one after another. Each of them was operated by the same temperature controller (TTC001, Thorlabs, Newton, NJ, United States) and current controller (TLD001, Thorlabs). The laser beam modulated by an optical chopper (model 300CD, Scitec Instruments, Trowbridge, United Kingdom) and passes the H-shaped PA cell [26] before hitting a power meter (thermal head 3A-FS-SH, Ophir Optronics, Jerusalem, Israel). The power meter records the transmission spectrum. The linear response of the PA signal as function of laser power according to Equation (1) is in the strict sense only valid for weak absorptions. Stronger ones weaken the PA spectrum’s amplitudes and require a further correction. The spectral regions that required correction were determined by examining the transmission spectra. A mirror and two irises are required for the laser alignment. The chopper frequency at 2750 Hz is set by a dual phase digital signal processing lock-in amplifier (LIA: model 7265, Signal Recovery, Oak Ridge, TN, United States). The LIA with a time constant set to 1 s is also responsible for the signal processing of the preamplified (PAS PMV 201, PAS-Tech, Hamburg, Germany) PA signal recorded by the analogue microphone (EM-158N, Primo, Tokyo, Japan). The system control and the data acquisition are performed with a computer. Further details of the ICLs can be found elsewhere [15].
An overview of the PAS gas flow system is illustrated in Figure 2.
Four gas lines, one for nitrogen (N2) and three for sample gases, lead to a connecting tube, which is attached to a pressure meter, the inlet of the mixing cell (for preprocessing) and to the PA cell. The gas lines consist partly of polytetrafluoroethylene and partly of stainless-steel tubes. The cell pressure was measured using a ceramic pressure transducer (PM 3) connected to the vacuum gauge (DVR 5, Vacuubrand, Wertheim, Germany). The transducer exhibits a resolution of 0.1 hPa in the range between 0.1 hPa and 10 hPa and a resolution of 1 hPa for pressures higher than 10 hPa. The outlets of the cells were connected to the vacuum pump, a chemistry-hybrid pump (RC 6, Vacuubrand), which enabled the evacuation of the gas lines including the cells down to 2.0 hPa. The sample mixtures were produced injecting a defined amount of one or more sample gases into the evacuated mixing cell and then adding N2 up to atmospheric pressure. Thereafter, the evacuated PA cell was filled with a small amount of the gas mixture from the mixing cell and, if required, additional amounts of different sample gases. Finally, the PA cell was filled with N2 up to atmospheric pressure to obtain the desired gas mixture. For some target concentrations, the usage of the mixing cell is not necessary. The concentrations, purities, and suppliers of the sample gases are listed in Table 1. The gases CH4, C2H6, and C3H8 were each present at natural abundance (NA).
Assuming that isotopes are randomly distributed throughout all isotopologues, the NA of the main methane isotopologues 12CH4 and 13CH4 are 98.8% and 1.11%, respectively [27]. The NA of the main ethane and propane isotopologues was estimated (based only on the 13C/12C ratio) to be 12CH3-12CH3: 97.8% 13CH3-12CH3: 2.2%, 13CH3-13CH3: 0.013%, 12CH3-12CH2-12CH3: 96.7%, and 13CH3-12CH2-12CH3: 2.2%, respectively.

2.3. Measurement and Data Processing

Measurements were performed on single and multiple component mixtures of sample gas with N2 under ambient conditions. In total, 62 (43 with ICL 1541 and 19 with ICL 1638) different SC hydrocarbon measurements/spectra were recorded. It should be noted, that a considerably larger dataset is usually required to set up a reliable prediction model. However, the primary goal of this research project was to show the feasibility. The ICL 1541 was used for all methane and ethane measurements, and the ICL 1638 for all propane measurements. The concentrations of the hydrocarbon used vary between 50 ppmv and 1 vol%. The individual concentrations for each measurement are listed in Table A1 of the Appendix A. The measured spectra are displayed in Figure A1 and Figure A2 of the Appendix A to give an impression of the difference.
Further details of the recording of the PA spectra with the ICLs and the processing of the spectra can be found elsewhere [15]. The processing involved averaging, power normalization, and wavelength correction.

2.4. Detection Limits

The signal-to-noise ratio (SNR) is commonly defined as the ratio of the signal power to the background noise power and is often expressed in dB. In this work, the SNR is calculated as the highest possible PA signal S max of a specific substance divided by the standard deviation of the PA signal during a nitrogen measurement σ N 2 :
SNR   =   S max σ N 2
where SNR is of dimension 1 and σ N 2 can be considered the noise level. PAS is a background-free technique, and, in principle, no PA signal is generated in absence of an absorber. Nitrogen does not show any absorption in the mid IR and therefore an N2 spectrum is perfectly suited to measure the background noise. The mean value of the PA signal during a N2 measurement represents an offset, e.g., due to window absorptions, and does not influence the SNR. The LOD is specified as the lowest quantity/concentration of a substance that can be distinguished from the absence of that substance. Here, the SNR shall be used for the calculation of the LOD of a substance:
LOD   =   c SNR
where c is the concentration of the substance during the SNR measurement. The obtained LOD is idealized, because it would be hard to distinguish a meaningful signal from the noise level, if both are of the same magnitude. Therefore, the value represents only a rough estimate of the experimentally achievable LOD. In order to measure the noise level of the setup, the sample cell was flushed with pure N2 and measurements were conducted under ambient conditions with both ICLs. The standard deviations of the N2 spectrum are 1.1 × 10−5 V for the ICL 1541 and 7.9 × 10−6 V for the ICL 1638, respectively.

2.5. Partial Least Squares Regression

Multicomponent measurements were processed with a PLSR model in the numerical software Matlab to evaluate the performance of the PA setup in detecting SC hydrocarbon isotopologue concentrations. The implementation of the PLSR method in Matlab is based on the SIMPLS algorithm developed by De Jong in 1993 [28]. The Matlab PLSR calculation function requires the gas spectra, the corresponding gas concentrations, and the number of selected latent variables (LVs) as input data. As output values, loading and score matrix of X and Y, explained variance (in %) by the model by each LV, and the PLSR coefficients are computed. The PLSR coefficients are required to calculate target values (concentrations) out of initial data (spectra):
Y res =   X pre B PLSR   +   G
where B PLSR is the PLSR coefficient matrix, X pre is the predictor matrix (spectra), Y res the response matrix (concentrations), and G the error matrix [23,24].
The predictive accuracy, i.e., the standard/mean error of a multivariate calibration model, is expressed as the root-mean-square error (RMSE). For that, the predicted residual error sum of squares (PRESS) is calculated following the equation [21]:
PRESS   =   i     = 1 n ( y i y ^ i ) 2
where y i denotes the initial data (true concentrations), y ^ i the predicted values (calculated concentrations) by the model and n the number of probes. The RMSE is then determined as
RMSE   =   PRESS n

3. Results

The linear behavior of the sensor was verified by evaluating the single-component gas measurements of the analytes. Representatively, the linearity of the single-component measurements of C2H6 and C3H8 are shown in Figure 3. The data demonstrate a good linearity with derived coefficients of determination of R2 = 0.991 (C2H6) and R2 = 0.997 (C3H8) and a small y-axis intercept of −0.162 V and −0.145 V, respectively.
The SNR and the LOD of the PAS sensor are calculated for each measured SC hydrocarbon isotopologue and are summarized in Table 2. In Figure 4, the respective PA spectra of the methane and the ethane isotopologues are shown. Although all these gases absorb across the entire wavelength range of the laser, each gas shows distinct absorption peaks, which makes them distinguishable from the others. In Figure 5, the PA spectra of the propane isotopologues are displayed. In contrast to the isotopologue spectra of methane and ethane, the propane isotopologue spectra show a broadband absorption structure with less sharp absorption peaks. In order to calculate the SNR, the maximum PA signal of each sample spectrum was determined from single component measurements with at low concentrations in N2. The SNR and LOD of the PAS sensor are different for each individual species. The PA signal depends equally on the absorption of the sample and on the laser power. This explains the noticeably higher LOD of CH4, because CH4 has a lower absorption strength (approximately one magnitude) in the spectral region of the laser source compared to the other compounds. It should be noted that due to the varying power performance of the ICLs at different wavelengths, the highest absorption peaks do not necessarily result in the maximum PA signal.
The method of choice for the validation of the PLSR model was a full cross validation (CV). Usually, all measurements have to be divided into calibration measurements to train the PLSR model (as a reference) and test measurements to verify the validity of the model. The CV is one of the most frequently used and efficient validation methods because it uses all recorded data for calibration as well as for validation [21]. In the full CV, n calibration models are created, where n stands for the number of probes (spectra). Each model uses n 1 probes for calibration and the left-out probe for validation. This is performed n times so that each probe serves once as validation. As part of the calibration, the PLSR model calculated five and two LVs for the measurements with the ICL 1541 and the ICL 1638, respectively, each LV representing the contribution of one hydrocarbon isotopologue. The 43/19 models were able to explain on average (96.7 ± 0.4)% (ICL 1541) and (99.1 ± 0.1)% (ICL 1638) of the total variance in the dataset regarding the target variable (Y).
Figure 6 shows two exemplary PA measurements of multi-component mixtures. Table 3 lists the five concentrations determined by the pressure readings and the PLRS models.
The standard approach to evaluate the prediction/quality of the model would be to specify the error (RMSE) according to Equations (5) and (6) of each individual calibration model and average it over all RMSE in the course of the CV (RMSECV). The RMSEs of Mixtures 29 and 40 equal 123 ppmv and 562 ppmv, respectively. If we leave CH4 due to its comparably high concentration out, we get an RMSE of Mixture 29 of 13 ppmv and an RMSE of Mixture 40 of 98 ppmv. Hence, we consider a separation of the error into the respective components to be more informative, since, for example, larger concentrations were selected for CH4 (with possible larger errors) in order to obtain PA signals of the same magnitude as the other components. However, it is also possible for the PLSR to predict negative concentrations. Real concentrations are logically greater or equal to zero. Here, the negative predicted concentrations were set to zero. That way, the error (RMSECV0) could be slightly improved. In addition, in order to obtain a better understanding of the error in relation to the respective concentration, the CV mean relative residual (MRR) of the measurements of the individual components was determined. All the CV results are summarized in Table 4.
The concentrations retrieved by the PLRS models can be obtained from Table A2 of the Appendix A. The residuals and the relative residuals of all predicted concentrations in the respective calibration model can be seen in Figure A3, Figure A4, Figure A5 and Figure A6 of the Appendix A, respectively.

4. Discussion

In general, the quality of the MVA strongly depends on the applied calibration method as well as on the used spectra. The PLSR is applicable as long as the individual components contributing to the overall spectrum are linearly independent of each other and possess characteristic features with signals larger than the noise level. The recorded spectra depend on the experimental setup, the absorption strength of the investigated compounds, as well as their concentration. The smallest RMSECV (28 ppmv) and MRR (7.3%) were found for 13C12C2H8 measured with the ICL 1638. For the ICL 1541, the largest RMSECV (604 ppmv) and at the same time smallest MRR (9.7%) can be assigned to CH4. This proves that the RMSE is highly dependent on the concentration, which in the case of CH4 is, on average, one magnitude larger than the other hydrocarbons. The largest MMRs (15.6–17.6%) were found for the ethane isotopologues measured with the ICL 1541. These discrepancies between the expected and the regressed concentrations are too high to be attributed only to the regression algorithm. Moreover, the estimated errors for some concentration values such as the CH4 in Mixture 29 or the C2H6 in Mixture 40 appear to be too small to justify the difference with the retrieved concentration values. Additionally, the variation for retrieved concentrations where a gas component was absent is relatively large. For CH4, it varies between −756.2 ppmv in Mixture 4 and +1241 ppmv in Mixture 40. All these results point out an underlying problem with the gas line management and/or the recorded mixed sample spectra. Especially with a limited dataset, each single measurement has a significant influence on the model. The uncertainty of the results originates (in the sense of error propagation) from the uncertainty in the sample preparation as well as from all other influences of the experimental setup. The sample preparation is associated with a certain degree of inaccuracy. The method of preparing the gas mixtures results in an uncertainty of the partial pressures and, therefore, the concentrations. After inserting sample gas into the cell and closing the valves, the gas line between the cross-adaptor and the cell was still filled with this component. When the cell was subsequently filled up to ambient pressure with N2, this remainder was pushed into the cell, thus falsifying the sample gas concentration. In addition, the uncertainty in the sample preparation of multiple components is considerably higher because the determination is solely calculated out of the pressure readings of the gas system. Furthermore, it can be assumed that the different components are not particularly well mixed in the gas system and mixing cell and that unwanted remaining substances in the system cannot be completely excluded. In fact, the different components are linearly dependent of each other because all gases at NA include small proportions of their 13C isotopologues, which also increases the error. Considering the fact that in this work only a simple PA setup was realized without substantial efforts to maximize the SNR or LOD, these results are very promising, especially since many improvements are still possible, e.g., cell optimization, signal enhancement with optimized amplifier and LIA parameters, and the application of direct modulation of the ICL instead of using the mechanical chopper.

5. Conclusions

Overall, the PA setup has demonstrated to work dependable in distinguishing SC hydrocarbon isotopologues and, in combination with the evaluation algorithm, enables an isotope-selective concentration determination of SC hydrocarbons. The current way of sample preparation offers room for improvement to minimize the concentration uncertainty. This can be achieved, for instance, with a suitable gas mixer, a system based on the combination of pressure controller and flow controller, and/or with the acquisition of hydrocarbon gas samples with accurate premixed concentrations. Afterwards, a proper design of experiment should be used for hydrocarbon concentration mixture combinations. This will help to plan necessary future measurements with synthetic multicomponent gas mixtures for the calibration model (required for the sensor’s concentration determination). In the future, the research will focus on several tasks concerning the optimization of the PA setup. The achievable LOD can be further enhanced by numerous experimental improvements. The used PA cell in this work has a relatively low Q-factor. An exchange of the current cell with an optimized PA cell with a higher Q-factor directly improves the PA signal and lowers the LOD. The feasibility of using a quartz tuning fork (QTF) as acoustic transducer will be further investigated, since QTFs can have Q values > 10,000 [29], which is unreachable with traditional PA cells. It should be noted that a reorientation to the quartz enhanced photoacoustic spectroscopy (QEPAS) technique (utilizing QTFs) is associated with several substantial changes of the PA setup. Commercially available QTFs work with a resonant frequency of ca. 32.8 kHz; the use of an optical chopper is therefore no longer possible. Therefore, the direct modulation of the ICL will be studied more intensively together with certain modulation methods like the wavelength modulation spectroscopy (WMS). Both changes offer the distinct advantage of efficient noise suppression [30,31]. In addition, the enhanced frequency stability of an applied direct modulation compared to an optical chopper together with the elimination of mechanical noise can result in better sensor sensitivity. The WMS approach was shown in a recent study using an ICL for the simultaneous detection of CH4, C2H6, and C3H8 (each at NA) using QEPAS [32]. In addition, the development of QEPAS is also still ongoing; therefore, further improvements will be achieved by implementing new generations of QTFs, overtone modes operation, and new micro-resonator systems [32].

Author Contributions

Conceptualization, A.L.; methodology, A.L.; software, A.L.; validation, A.L.; formal analysis, A.L.; investigation, A.L.; resources, M.W.; data curation, A.L.; writing—original draft preparation, A.L.; writing—review and editing, M.W.; visualization, A.L.; supervision, M.W.; project administration, A.L. and M.W.; funding acquisition, M.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Hydrocarbon concentrations of the gas mixtures in ppmv. The compounds CH4, C2H6, and C3H8 are at NA. The uncertainties of the values were calculated from the purity of the sample gas and the error of the pressure readings (±1 hPa) by error propagation.
Table A1. Hydrocarbon concentrations of the gas mixtures in ppmv. The compounds CH4, C2H6, and C3H8 are at NA. The uncertainties of the values were calculated from the purity of the sample gas and the error of the pressure readings (±1 hPa) by error propagation.
ICL 1541 ICL 1638
No.C2H613C12CH613C2H6CH413CH4No.C3H813C12C2H8
1500 ± 1100001200 ± 100
2250 ± 1000002400 ± 110
310000003800 ± 110
41000 ± 12000040800 ± 100
5500 ± 11500 ± 6600050400 ± 55
6750 ± 11250 ± 3800060210 ± 34
7200 ± 10100 ± 2100070130 ± 25
8300 ± 10600 ± 7700081000
9400 ± 11300 ± 440009240 ± 10300 ± 44
10150 ± 10200 ± 3200010350 ± 11500 ± 66
1101000 ± 12200011500 ± 11440 ± 59
120400 ± 5500012150 ± 10100 ± 21
13001000 ± 1220013610 ± 11250 ± 38
1400500 ± 660014310 ± 10350 ± 49
15370 ± 11150 ± 27150 ± 270015540 ± 11540 ± 70
16210 ± 10200 ± 32200 ± 320016450 ± 11700 ± 88
170350 ± 49350 ± 490017200 ± 10100 ± 21
180100 ± 21100 ± 210018900 ± 11600 ± 77
190005000 ± 180191000 ± 1250 ± 16
20000100000
210002500 ± 140
220000500 ± 66
230007000 ± 21200 ± 32
240003000 ± 15300 ± 44
250001000 ± 12400 ± 55
26210 ± 10220 ± 35320 ± 460220 ± 35
270330 ± 47480 ± 640330 ± 47
28130 ± 10300 ± 4402000 ± 13150 ± 27
29190 ± 10200 ± 3204000 ± 16100 ± 21
30200 ± 10240 ± 3703000 ± 150
31300 ± 10360 ± 5005000 ± 180
320170 ± 29340 ± 486000 ± 190
330260 ± 39510 ± 671500 ± 120
340150 ± 27140 ± 260230 ± 36
350300 ± 44280 ± 410460 ± 62
360240 ± 37260 ± 390390 ± 54
370100 ± 21100 ± 210160 ± 28
380210 ± 34210 ± 343500 ± 15200 ± 32
390140 ± 26140 ± 264600 ± 17130 ± 25
40410 ± 110100 ± 210400 ± 55
41110 ± 100150 ± 270600 ± 77
42260 ± 100230 ± 360250 ± 38
431000160 ± 280170 ± 29
Table A2. Hydrocarbon concentrations of the gas mixtures in ppmv determined by the PLRS models.
Table A2. Hydrocarbon concentrations of the gas mixtures in ppmv determined by the PLRS models.
ICL 1541 ICL 1638
No.C2H613C12CH613C2H6CH413CH4No.C3H813C12C2H8
1419.423.6−16.8−148.64.01169.621.9
2249.318.1−7.5−30.914.42365.111.0
3121.016.41.1193.019.53835.24.7
4966.164.2−57.0−756.2−52.9442.5743.2
5480.0395.1−18.2238.4−6.4519.0359.2
6636.4272.2−28.3−439.3−16.769.3218.6
7203.179.30.4−25.018.87−0.5147.0
8281.7592.1−3.4−88.5−1.3874.425.6
9432.1293.2−12.4−8.8−5.79260.3272.2
10175.0217.10.3−9.215.910389.4466.7
1184.2834.214.4−31.823.111540.4415.3
1244.0394.76.495.524.912130.1108.6
13−35.3−37.6819.5−163.7−7.313633.0222.7
143.6−7.5456.877.411.314275.7370.9
15418.8122.3193.5117.6−15.215565.1550.5
16156.3165.7286.8528.0−2.316415.5731.5
1719.4304.9377.03.512.017210.61054.1
1827.592.4132.2236.421.118777.5619.7
1920.214.337.65750.5−1.1191038.943.9
20−14.5−25.8−33.86748.30.0
2124.23.06.52838.318.2
22−23.7−20.0−32.8−49.5369.7
23−19.4−34.3−6.87337.5154.0
24−18.3−26.2−23.03238.5272.7
25−25.4−26.8−31.71117.7372.6
26174.3300.5286.127.9174.3
2715.1386.8497.5−199.5260.3
28141.0358.4−13.22022.8171.3
29205.1196.3−19.34273.1108.8
30181.7342.6−2.12795.913.0
31299.7392.2−17.25531.2−1.1
32−13.0170.6366.36355.63.6
33−4.5250.8557.51644.18.1
34−2.7201.4132.9−171.1227.1
35−43.8316.0256.6−584.8429.3
36−37.0323.1249.2−413.9400.6
376.2138.5106.1−10.9169.8
38−22.2237.3193.83630.9226.4
39−9.8123.7114.44874.5136.0
40584.3−66.4125.71241.4459.2
4151.0−70.7164.9−135.6712.0
42330.9−24.1352.2−412.2270.7
43149.1−8.4200.3−112.3176.2
Figure A1. Measured PA spectra (43 different measurements) for multicomponent analysis with ICL 1541.
Figure A1. Measured PA spectra (43 different measurements) for multicomponent analysis with ICL 1541.
Molecules 25 02266 g0a1
Figure A2. Measured PA spectra (19 different measurements) for multicomponent analysis with ICL 1638.
Figure A2. Measured PA spectra (19 different measurements) for multicomponent analysis with ICL 1638.
Molecules 25 02266 g0a2
Figure A3. Cross validation residual of the ICL 1541 measurements. The y-axis on the right side refers to CH4 (green).
Figure A3. Cross validation residual of the ICL 1541 measurements. The y-axis on the right side refers to CH4 (green).
Molecules 25 02266 g0a3
Figure A4. Cross validation relative residual of the ICL 1541 measurements.
Figure A4. Cross validation relative residual of the ICL 1541 measurements.
Molecules 25 02266 g0a4
Figure A5. Cross validation residual of the ICL 1638 measurements.
Figure A5. Cross validation residual of the ICL 1638 measurements.
Molecules 25 02266 g0a5
Figure A6. Cross validation relative residual of the ICL 1638 measurements.
Figure A6. Cross validation relative residual of the ICL 1638 measurements.
Molecules 25 02266 g0a6

References

  1. Myhre, G.; Shindell, D.; Bréon, F.-M.; Collins, W.; Fuglestvedt, J.; Huang, J.; Koch, D.; Lamarque, J.-F.; Lee, D.; Mendoza, B.; et al. Anthropogenic and Natural Radiative Forcing. In Climate Change 2013: The Physical Science Basis; Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Ed.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar]
  2. Warneck, P. Chemistry of the Natural Atmosphere, 2nd ed.; International geophysics series; Academic: San Diego, CA, USA; London, UK, 2000; ISBN 0-12-735632-0. [Google Scholar]
  3. Chen, T.-M.; Kuschner, W.G.; Gokhale, J.; Shofer, S. Outdoor air pollution: Ozone health effects. Am. J. Med. Sci. 2007, 334, 244–248. [Google Scholar] [CrossRef]
  4. van Dingenen, R.; Dentener, F.; Raes, F.; Krol, M.; Emberson, L.; Cofala, J. The global impact of ozone on agricultural crop yields under current and future air quality legislation. Atmos. Environ. 2009, 43, 604–618. [Google Scholar] [CrossRef]
  5. Shindell, D.; Kuylenstierna, J.C.I.; Vignati, E.; Van Dingenen, R.; Amann, M.; Klimont, Z.; Anenberg, S.C.; Muller, N.; Janssens-Maenhout, G.; Raes, F.; et al. Simultaneously mitigating near-term climate change and improving human health and food security. Science 2012, 335, 183–189. [Google Scholar] [CrossRef] [Green Version]
  6. Gensch, I.; Kiendler-Scharr, A.; Rudolph, J. Isotope ratio studies of atmospheric organic compounds: Principles, methods, applications and potential. Int. J. Mass Spectrom. 2014, 365, 206–221. [Google Scholar] [CrossRef]
  7. U.S. Energy Information Administration. International Energy Outlook 2019: With Projections to 2050; U.S. Energy Information Administration: Washington, DC, USA, 2019. Available online: https://www.eia.gov/outlooks/ieo/pdf/ieo2019.pdf (accessed on 10 May 2020).
  8. Faramawy, S.; Zaki, T.; Sakr, A.-E. Natural gas origin, composition, and processing: A review. J. Nat. Gas Sci. Eng. 2016, 34, 34–54. [Google Scholar] [CrossRef]
  9. Zumberge, J.E.; Ferworn, K.; Brown, S. Isotopic reversal (‘rollover’) in shale gases produced from the Mississippian Barnett and Fayetteville formations. Mar. Pet. Geol. 2012, 31, 43–52. [Google Scholar] [CrossRef]
  10. Xia, X.; Chen, J.; Braun, R.; Tang, Y. Isotopic reversals with respect to maturity trends due to mixing of primary and secondary products in source rocks. Chem. Geol. 2013, 339, 205–212. [Google Scholar] [CrossRef]
  11. Tilley, B.; Muehlenbachs, K. Isotope reversals and universal stages and trends of gas maturation in sealed, self-contained petroleum systems. Chem. Geol. 2013, 339, 194–204. [Google Scholar] [CrossRef]
  12. Ellis, L.; Brown, A.; Schoell, M.; Uchytil, S. Mud gas isotope logging (MGIL) assists in oil and gas drilling operations. Oil Gas J. 2003, 101, 32–41. [Google Scholar]
  13. Sessions, A.L. Isotope-ratio detection for gas chromatography. J. Sep. Sci. 2006, 29, 1946–1961. [Google Scholar] [CrossRef]
  14. Vurgaftman, I.; Weih, R.; Kamp, M.; Meyer, J.R.; Canedy, C.L.; Kim, C.S.; Kim, M.; Bewley, W.W.; Merritt, C.D.; Abell, J.; et al. Interband cascade lasers. J. Phys. D Appl. Phys. 2015, 48, 123001. [Google Scholar] [CrossRef]
  15. Loh, A.; Wolff, M. High resolution spectra of 13C ethane and propane isotopologues photoacoustically measured using interband cascade lasers near 3.33 and 3.38 µm, respectively. J. Quant. Spectrosc. Radiat. Transf. 2019, 227, 111–116. [Google Scholar] [CrossRef]
  16. Loh, A.; Wolff, M. Absorption cross sections of 13C ethane and propane isotopologues in the 3 µm region. J. Quant. Spectrosc. Radiat. Transf. 2017, 203, 517–521. [Google Scholar] [CrossRef]
  17. Bell, A.G. On the production and reproduction of sound by light. Am. J. Sci. 1880, 3, 305–324. [Google Scholar] [CrossRef] [Green Version]
  18. Harren, F.J.M.; Christescu, S.M. Photoacoustic Spectroscopy in Trace Gas Monitoring: (2019 version). In Encyclopedia of Analytical Chemistry; Meyers, R.A., Ed.; Wiley: Hoboken, NJ, USA, 2006; ISBN 9780470027318. [Google Scholar]
  19. Elia, A.; Lugarà, P.M.; Di Franco, C.; Spagnolo, V. Photoacoustic techniques for trace gas sensing based on semiconductor laser sources. Sensors 2009, 9, 9616–9628. [Google Scholar] [CrossRef] [Green Version]
  20. Patimisco, P.; Scamarcio, G.; Tittel, F.K.; Spagnolo, V. Quartz-enhanced photoacoustic spectroscopy: A review. Sensors 2014, 14, 6165–6206. [Google Scholar] [CrossRef] [Green Version]
  21. Kessler, W. Multivariate Datenanalyse in der Bio- und Prozessanalytik: Mit Beispielen aus der Praxis; Wiley-VCH: Weinheim, Germany, 2006; ISBN 9783527312627. [Google Scholar]
  22. Hair, J.F.; Black, W.C.; Babin, B.J.; Anderson, R.E. Multivariate Data Analysis, 7th ed.; Pearson Custom Library Prentice Hall: Upper Saddle River, NJ, USA, 2010; ISBN 978-1-292-02190-4. [Google Scholar]
  23. Saalberg, Y.; Wolff, M. Multivariate Analysis as a Tool to Identify Concentrations from Strongly Overlapping Gas Spectra. Sensors 2018, 18, 1562. [Google Scholar] [CrossRef] [Green Version]
  24. Wold, S.; Sjöström, M.; Eriksson, L. PLS-regression: A basic tool of chemometrics. Chemom. Intell. Lab. Syst. 2001, 58, 109–130. [Google Scholar] [CrossRef]
  25. Sigrist, M.W. Laser photoacoustic spectrometry for trace gas monitoring. Analyst 1994, 119, 525–531. [Google Scholar] [CrossRef]
  26. Nodo, E. Optimization of resonant cell design for optoacoustic gas spectroscopy (H-type). Appl. Opt. 1978, 17, 1110–1119. [Google Scholar] [CrossRef]
  27. Stolper, D.; Sessions, A.L.; Ferreira, A.; Neto, E.S.; Schimmelmann, A.; Shusta, S.; Valentine, D.L.; Eiler, J. Combined 13C–D and D–D clumping in methane: Methods and preliminary results. Geochim. Cosmochim. Acta 2014, 126, 169–191. [Google Scholar] [CrossRef]
  28. de Jong, S. SIMPLS: An alternative approach to partial least squares regression. Chemom. Intell. Lab. Syst. 1993, 18, 251–263. [Google Scholar] [CrossRef]
  29. Patimisco, P.; Sampaolo, A.; Dong, L.; Tittel, F.K.; Spagnolo, V. Recent advances in quartz enhanced photoacoustic sensing. Appl. Phys. Rev. 2018, 5, 011106. [Google Scholar] [CrossRef]
  30. Schilt, S.; Thévenaz, L. Wavelength modulation photoacoustic spectroscopy: Theoretical description and experimental results. Infrared Phys. Technol. 2006, 48, 154–162. [Google Scholar] [CrossRef]
  31. Li, J.; Du, Z.; An, Y. Frequency modulation characteristics for interband cascade lasers emitting at 3 μm. Appl. Phys. B 2015, 121, 7–17. [Google Scholar] [CrossRef]
  32. Sampaolo, A.; Csutak, S.; Patimisco, P.; Giglio, M.; Menduni, G.; Passaro, V.; Tittel, F.K.; Deffenbaugh, M.; Spagnolo, V. Methane, ethane and propane detection using a compact quartz enhanced photoacoustic sensor and a single interband cascade laser. Sens. Actuators B Chem. 2019, 282, 952–960. [Google Scholar] [CrossRef]
Figure 1. Experimental photoacoustic (PA) setup (block diagram).
Figure 1. Experimental photoacoustic (PA) setup (block diagram).
Molecules 25 02266 g001
Figure 2. Photoacoustic spectroscopy (PAS) gas flow system. MVs are metering valves, PMs are pressure meters, RV is a pressure reducing valve, PVs are plug valves, and NVs are needle valves. The metering valves MV 1 and MV 3–5 are part of the pressure reducers. Pressure meters PM 1–2 and the reducing valve RV 1 belong to the N2 gas cylinder’s pressure reducer. The plug valves PV 1–6 separate the individual gas lines and PV 7–10 seal the mixing and the PA cell, respectively. The MV 2 and needle valves NV 1–4 allow a gentle filling of the cells. The MV 6 allows the adjustment of the vacuum pump power and the PV 11 and PV 12 close/open the purging outlet and the vacuum pump, respectively.
Figure 2. Photoacoustic spectroscopy (PAS) gas flow system. MVs are metering valves, PMs are pressure meters, RV is a pressure reducing valve, PVs are plug valves, and NVs are needle valves. The metering valves MV 1 and MV 3–5 are part of the pressure reducers. Pressure meters PM 1–2 and the reducing valve RV 1 belong to the N2 gas cylinder’s pressure reducer. The plug valves PV 1–6 separate the individual gas lines and PV 7–10 seal the mixing and the PA cell, respectively. The MV 2 and needle valves NV 1–4 allow a gentle filling of the cells. The MV 6 allows the adjustment of the vacuum pump power and the PV 11 and PV 12 close/open the purging outlet and the vacuum pump, respectively.
Molecules 25 02266 g002
Figure 3. Linearity of the normalized PA signal (integrated over the whole spectrum) obtained of individual C2H6 and C3H8 measurements at different concentrations.
Figure 3. Linearity of the normalized PA signal (integrated over the whole spectrum) obtained of individual C2H6 and C3H8 measurements at different concentrations.
Molecules 25 02266 g003
Figure 4. PA spectra measured with interband cascade lasers (ICL) 1541 of sample gases with concentrations according to Table 2.
Figure 4. PA spectra measured with interband cascade lasers (ICL) 1541 of sample gases with concentrations according to Table 2.
Molecules 25 02266 g004
Figure 5. PA spectra measured with ICL 1638 of sample gases with concentrations according to Table 2.
Figure 5. PA spectra measured with ICL 1638 of sample gases with concentrations according to Table 2.
Molecules 25 02266 g005
Figure 6. PA spectra measured with ICL 1541 of two multi-component mixtures with concentrations according to Table 3.
Figure 6. PA spectra measured with ICL 1541 of two multi-component mixtures with concentrations according to Table 3.
Molecules 25 02266 g006
Table 1. Sample gases used in this work.
Table 1. Sample gases used in this work.
Sample GasConcentrationPuritySupplier
CH4 (at NA)1 vol%3.5Air Liquide S.A., Paris, France
13CH4pure2.0Sigma-Aldrich, Inc., St. Louis, MO, United States
C2H6 (at NA)100 ppmv6.0Linde AG., Dublin, Ireland
1 vol%3.5Air Liquide S.A.
13C12CH6pure2.0Sigma-Aldrich, Inc.
13C2H6pure2.0Sigma-Aldrich, Inc.
C3H8 (at NA)100 ppmv6.0Linde AG.
1 vol%3.5Air Liquide S.A.
13C12C2H8pure2.0Cambridge Isotope Laboratories, Inc.
Tewksbury, MA, United States
N2pure5.0Linde AG.
Table 2. Idealized detection limits of the PA setup. The signal-to-noise ratio (SNR) of each sample spectrum was obtained at these wavelengths.
Table 2. Idealized detection limits of the PA setup. The signal-to-noise ratio (SNR) of each sample spectrum was obtained at these wavelengths.
Sample GasConcentration (ppmv)Wavelength (nm)SNRLOD (ppmv)
CH4 (at NA)25003328.4728.93.4
13CH45003329.94292.50.12
C2H6 (at NA)1003333.01449.70.069
13C12CH64003329.93515.00.11
13C2H65003332.23968.50.13
C3H8 (at NA)1003381.91564.30.064
13C12C2H81303382.23039.80.043
Table 3. Concentrations of two exemplary chosen mixtures determined by the PLRS models.
Table 3. Concentrations of two exemplary chosen mixtures determined by the PLRS models.
Sample GasConcentration
of Mixture 29 (ppmv)
PLRS
Concentration (ppmv)
Concentration
of Mixture 40 (ppmv)
PLRS
Concentration (ppmv)
CH4 (at NA) 4000 ± 16427301241
13CH4100 ± 21108400 ± 55459
C2H6 (at NA)190 ± 10205410 ± 11584
13C12CH6200 ± 321960−66
13C2H60−19100 ± 21126
Table 4. Summary of the cross validation (CV) results.
Table 4. Summary of the cross validation (CV) results.
Sample GasRMSECV (ppmv)RMSECV0 (ppmv)MRR (%)
CH4 (at NA)6045729.7
13CH4353412.4
C2H6 (at NA)464417.3
13C12CH6494515.6
13C2H6434117.6
C3H8 (at NA)40409.7
13C12C2H828287.3

Share and Cite

MDPI and ACS Style

Loh, A.; Wolff, M. Multivariate Analysis of Photoacoustic Spectra for the Detection of Short-Chained Hydrocarbon Isotopologues. Molecules 2020, 25, 2266. https://doi.org/10.3390/molecules25092266

AMA Style

Loh A, Wolff M. Multivariate Analysis of Photoacoustic Spectra for the Detection of Short-Chained Hydrocarbon Isotopologues. Molecules. 2020; 25(9):2266. https://doi.org/10.3390/molecules25092266

Chicago/Turabian Style

Loh, Alain, and Marcus Wolff. 2020. "Multivariate Analysis of Photoacoustic Spectra for the Detection of Short-Chained Hydrocarbon Isotopologues" Molecules 25, no. 9: 2266. https://doi.org/10.3390/molecules25092266

APA Style

Loh, A., & Wolff, M. (2020). Multivariate Analysis of Photoacoustic Spectra for the Detection of Short-Chained Hydrocarbon Isotopologues. Molecules, 25(9), 2266. https://doi.org/10.3390/molecules25092266

Article Metrics

Back to TopTop