FIBRE OPTIC FILTER REMOTE GAS CORRELATION SENSOR

Abstract

A method for detecting and identifying a remote gas, the method comprising the steps of: receiving a light wave associated with the remote gas; coupling the light wave into a single mode fibre; transmitting the light wave via the single mode fibre into a filter comprising a fibre-based tunable cavity; modulating the cavity length of the filter transmission window to cause a detected modulated signal that is proportional to the spectral feature; and processing the signal using a lock-in amplifier capable of low-pass filtering and out-of-frequency noise rejection.

Claims

1. A method for detecting and identifying a remote gas, the method comprising the steps of: receiving a light wave associated with the remote gas; coupling the light wave into a single mode fibre; transmitting the light wave via the single mode fibre into a filter comprising a fibre-based tunable cavity; modulating the cavity length of the filter transmission window to cause a detected modulated signal that is proportional to the spectral feature; and processing the signal using a lock-in amplifier capable of low-pass filtering and out-of-frequency noise rejection.

2. The method of claim 1, wherein the filter is an all-fibre Fabry-Pérot filter.

3. The method of claim 1, wherein the remote gas is detected based on its absorption or emission features in at least one of a visible portion, mid-infrared portion and near-infrared portion of an electromagnetic wave spectrum.

4. The method of claim 1, wherein the remote gas is detected based on its absorption features in a pure carbon dioxide spectrum.

5. The method of claim 1, wherein the remote gas is detected based on its absorption features in a pure carbon monoxide spectrum.

6. The method of claim 1, wherein the remote gas is detected based on its absorption features in a mixture of carbon dioxide and carbon monoxide spectrum.

7. The method of claim 1, wherein the remote gas is detected based on its absorption features in a methane spectrum.

8. The method of claim 1, wherein the remote gas is detected based on its absorption features in a measured near-infrared spectrum of Venus.

9. The method of claim 1, wherein a sensor comprising the lock-in amplifier and a single channel detector is tailored to detect a gas comprising quasi-periodic absorption lines.

10. The method of claim 1, wherein a sensor comprising the lock-in amplifier and a single channel detector is tailored to detect a gas comprising complex absorption spectra.

11. The method of claim 10, wherein the sensor comprising the lock-in amplifier and a single channel detector is tailored for monitoring emissions from remote sources.

12. The method of claim 10, wherein the sensor comprising the lock-in amplifier and a multi-channel detector is tailored for monitoring remote gases.

13. The method of claim 11, wherein the sensor is tailored for exoplanet detection and characterization.

14. The method of claim 11, wherein the sensor is tailored for remote sensing.

15. The method of claim 11, wherein the sensor extracts radial velocities of the remote gas from a phase of the lock-in correlation amplitude for identifying the remote gas.

16. A gas sensing system comprising: a fibre collimator for receiving a light wave associated with a remote gas; a single mode fibre for receiving the light wave coupled thereinto by the fibre collimator; a filter comprising a fibre-based tunable cavity, wherein the cavity receives the light wave via the single mode fibre; a modulator for modulating the cavity length; a single-channel photodetector for detecting a modulated signal that is proportional to the spectral feature; a single channel lock-in amplifier capable of low-pass filtering and out-of-frequency noise rejection.

17. The gas sensing system of claim 16, wherein the remote gas is detected based on its absorption features in at least one of a visible portion, mid-infrared portion and near-infrared portion of an electromagnetic wave spectrum.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] FIG. 1 shows a schematic of a fibre optic correlation spectroscopy system;

[0023] FIG. 2 shows a transmission spectra of a tunable fibre Fabry-Perot (FP) filter as a function of electrical bias;

[0024] FIG. 3a shows a simulated transmission spectrum of the tunable FP filter as it is tuned to 1578.66 nm without any gas present;

[0025] FIG. 3b shows a simulated transmission spectrum of the tunable FP filter as it is tuned to 1578.66 nm with carbon dioxide (CO.sub.2) present;

[0026] FIG. 4a shows a near-infrared spectrum of Venus acquired using the 3.58-meter Galileo National Telescope on Santa Cruz de La Palma with strong CO.sub.2 absorption bands centred at 1580 nm and 1610 nm;

[0027] FIG. 4b shows a detailed view of the absorption line profiles, where the absorption peaks of both Earth and Venus are overlapping, separated by approximately 60 pm;

[0028] FIG. 5 shows a modulated signal as a function of measuring steps over time for various CO.sub.2 absorption depths. The signal amplitude increases as a function of CO.sub.2 absorption depth;

[0029] FIG. 6 shows a correlation signal as a function of phase showing the modulation amplitude increasing with absorption depth of the Venus spectrum;

[0030] FIG. 7 shows a simulated correlation signal output from a lock-in amplifier output as a function of phase shift for various radial velocities;

[0031] FIG. 8a shows a phase shift as a function of radial velocity of various redshifted spectra;

[0032] FIG. 8b shows a voltage of a lock-in signal plotted as a function of gas pressure; and

[0033] FIG. 9 shows the radial velocity extracted from the Venus spectrum for increasing harmonics/modulation ranges of detection.

DETAILED DESCRIPTION

[0034] The following detailed description refers to the accompanying drawings. Wherever possible, the same reference numbers are used in the drawings and the following description to refer to the same or similar elements. While embodiments of the disclosure may be described, modifications, adaptations, and other implementations are possible. For example, substitutions, additions, or modifications may be made to the elements illustrated in the drawings, and the methods described herein may be modified by substituting, reordering, or adding stages to the disclosed methods. Accordingly, the following detailed description does not limit the disclosure. Instead, the proper scope of the disclosure is defined by the appended claims.

[0035] Moreover, it should be appreciated that the particular implementations shown and described herein are illustrative of the invention and are not intended to otherwise limit the scope of the present invention in any way. Indeed, for the sake of brevity, certain sub-components of the individual operating components, and other functional aspects of the systems may not be described in detail herein. Furthermore, the connecting lines shown in the various figures contained herein are intended to represent exemplary functional relationships and/or physical couplings between the various elements. It should be noted that many alternative or additional functional relationships or physical connections may be present in a practical system.

[0036] Referring to FIG. 1, there is shown a schematic of a fibre optic correlation spectroscopy system 10 capable of remote gas detection, identification and radial velocity measurements using a single channel lock-in amplification technique. System 10 comprises fibre collimator 12 which receives light 14 from free space and couples the light 14 into single mode fibre 16 for transmission into filter 18 with fibre-based tunable cavity 20. In one implementation, filter 18 is a Fabry-Perot (FP) filter with resonant cavity 20 that generates longitudinal modes at resonant wavelengths. FFP filter 18 transmission window is modulated through the spectral features, and single-channel photodetector 22 detects the modulated signal that is proportional to the spectral feature of the remote gas. Lock-in amplifier 24 processes the signal by performing low-pass filtering and out-of-frequency noise rejection.

[0037] In a simulation, the transmission spectrum of a FP mode is Lorenztian in shape [21], with a spectral profile as a function of wavelength according to

[00001]$T_{\mathrm{FP}}\left(\lambda ,{\lambda }_0,\gamma \right)=\left(\frac{1}{\pi \gamma }\right)\frac{{\gamma }^2}{{\gamma }^2+{\left(\lambda -{\lambda }_0\right)}^2},$

[0038] where γ is the half width at half max, λ is the wavelength of light, and λ.sub.0 is the centre wavelength of the peak. Other resonances are separated by over 100 nm in bandwidth so that only a single transmission window is visible over a large bandpass. The FFP filter 18 transmission spectrum was measured as a function of electrical bias, as shown in FIG. 2. The non-uniform transmission is due to the super luminescent light emitting diode (SLED) emission spectrum. In the simulation, the fibre FP filter 18 transmission spectrum is simulated using a Lorenztian profile with a spectral width that matches that of the experimental filter (15 GHz or 220 pm).

[0039] The gas absorption lines take the shape of Voigt profiles:

A(λ, T, P)=∫.sub.−∞.sup.∞G(λ.sub.0, σ)L(λ−λ.sub.0, γ)dλ

[0040] where L is the Lorenztian profile and G is the centred Gaussian profile. For simplicity and computational Voigt profile for each absorption line of a gas are generated using a rational approximation in simulation by using a rapid MATLAB subroutine [22].

[0041] The product overlap between the absorption line spectrum and the FFP transmission spectrum is expressed as:

S.sub.FP(λT)=A(λ).Math.T.sub.FP(λT),

[0042] and represents the output spectrum from the FP filter. At single-channel photodetector 22, the signal is the integral of the product above:

C(T)=∫.sub.λ.sub.1.sup.λ.sup.2S.sub.FP(λT)dλ.

[0043] The output from filter 18 is calculated as the product of the FFP filter 18 transmission and the normalized gas absorption spectrum. The input spectrum is modeled as a flat broadband spectrum with normally-distributed white noise. The absorption spectra of the gas is generated by importing absorption cross-sections from the high-resolution transmission molecular absorption database (HITRAN) and generating approximate Voigt profiles based on room temperature and pressure conditions. The FFP filter 18 transmission profile is controlled by setting the Voigt parameter sigma. The transmission central wavelength of filter 18 is modulated across a spectral absorption feature of CO.sub.2. The simulated overlap of the absorption feature with the FFP filter 18 transmission window is shown in FIGS. 3a and 3b. In more detail, FIG. 3a shows a simulated transmission spectrum of tunable FP filter 18 as it is tuned to 1578.66 nm without any gas present, while FIG. 3b shows a simulated transmission spectrum of tunable FP filter 18 as it is tuned to 1578.66 nm with CO.sub.2 present.

[0044] In one experimental setup, the reflection spectrum of Venus is measured. A near-infrared spectrum of Venus was acquired using the 3.58-meter Galileo National Telescope on Santa Cruz de La Palma and is shown in FIG. 4a. The spectrum reveals strong CO.sub.2 absorption bands centred at 1580 nm and 1610 nm. A more detailed view of the absorption line profiles is also shown in FIG. 4b, where the absorption peaks of both Earth and Venus are overlapping, separated by approximately 60 pm. The telluric lines are removed from spectrum to produce the approximate spectrum of Venus alone. This spectrum is used as an input to the simulations to demonstrate the FFP filter-based detection of remote CO.sub.2 detection.

[0045] For the detection of CO.sub.2 using the absorption lines around 1580 nm, a sinusoidal modulation of the simulated FFP filter 18 transmission profile is introduced. The lock-in amplification is simulated according to the equation below as the transmission window of the FFP filter 18 is swept:

[00002]$V_{out}=\frac{1}{T}{\int }_{t-T}^t\sin \left[2\pi f_{\mathrm{ref}}{t}^{*}+\phi \right]V_{in}d{t}^{*}$

[0046] FIG. 5 shows a modulated signal as a function of measuring steps over time for various CO.sub.2 absorption depths. The signal amplitude increases as a function of CO.sub.2 absorption depth.

[0047] FIG. 6 shows a correlation signal as a function of phase showing the modulation amplitude increasing with absorption depth of the Venus spectrum. As can be seem, with increasing CO.sub.2 absorption depth, the amplitude of DC signal increases.

[0048] One of the benefits of lock-in amplifier 24 is the ability to measure the input signal phase. The phase of the lock-in signal from light passing through fiber FP filter 18 is proportional to the wavelength of the spectral feature. The two components of the lock-in signal can be used to obtain the phase of the lock-in signal using the equations below.

[00003]$$\begin{gathered} V_{\mathrm{sig}}=\frac{2V_{out}}{\cos \theta } \\ X=V_{\mathrm{sig}}\cos \theta \\ Y=V_{\mathrm{sig}}\sin \theta \\ \theta ={\tan }^{-1}\left(\frac{Y}{X}\right) \end{gathered}$$

[0049] In the case of a Doppler shift, the spectral shift is proportional to the radial velocity of the object. At 1578 nm, where carbon dioxide can be detected, the spectral shift can be over 1 nm with radial velocities up to 200 km/s. The DC signal from lock-in amplifier 24 is simulated as a function of spectral shift of a pure CO.sub.2 spectrum due to non-zero radial velocity. FIG. 7 shows a simulated correlation signal output from lock-in amplifier 24 output as a function of phase shift for various radial velocities, with a change in the phase of the lock-in signal.

[0050] In the case of an ideal, single gas absorption spectrum, the radial velocity can be accurately determined by the phase shift of the lock-in signal. However with a complex spectrum with multiple absorption lines in the scanning region, the radial velocity to phase relation is affected from the spectral lines in the modulation range, leading to deviations in the radial velocity determination unless a calibration for that same spectrum is performed. Since the spectrum in most cases is not known, this can present a challenge in accurate determination of the radial velocity under complex spectra.

[0051] As there are multiple, roughly periodically spaced absorption lines for many gases, such as CO.sub.2 and carbon monoxide (CO), these additional lines can be used to increase the selectivity of the gas sensing technique and reduce errors in the radial velocity measurement. By scanning in a larger range that is a multiple of the average line spacing, a higher order harmonic of the lock-in signal would also indicate the presence of the gas. This technique is less sensitive to other spectral features from other gases since their contribution is reduced relative to the signal from the target gas. This technique is used to extract the sixth harmonic signal from a modulation range that is six times the average spacing between CO.sub.2 lines (6×0.47 nm), centred at 1580.2 nm, as shown in FIG. 8a. Gas identification is thus accomplished through their phase offsets relative to a calibrated phase for that gas. When there are both gas species, the phase of the lock-in signal is between that of both gases, with the measured phase being a function of each gases contribution to the total absorption over the modulation range of FFP filter 18.

[0052] Accordingly, an acquired spectrum of Venus may be characterized by system 10 using a simulated tunable FFP filter 18 transmission spectrum, a FFP spectral modulation amplitude sufficient to cover 1 nm spectral bandwidth and processing the output of the FFP analog lock-in amplifier 24. Calibration of the sensor with a non-Doppler shifted spectral lines is performed to obtain reference phase signal for radial velocity determination. A phase shift relative to the reference phase can be related to a Doppler shifted feature. In FIG. 8a, the phase of the lock-in signal is plotted as a function of simulated radial velocity for various spectra. For various CO.sub.2-optimized modulation ranges, the phase shift as a function of radial velocity is relatively unchanged. With the addition of CO, relationship is offset towards the trend for pure CO, which indicates the possibility to use the phase offset to infer relative gas ratios. The reflection spectrum of Venus is nearly indistinguishable to that of pure CO.sub.2, due to the strong CO.sub.2 absorption lines in the spectrum, showing that the technique can operate under complex spectra such as that of Venus. In FIG. 8b, the voltage signal (V) of the lock-in signal is plotted as a function of gas pressure, derived from the experimental data.

[0053] It should be noted that since the simulation uses a spectrum of Venus acquired with a spectrometer with limited spectral resolution, the spectrum contains artificially broadened spectral lines that significantly reduces the accuracy of the radial velocity determination. The radial velocity of Venus at the point of spectral acquisition by the TNG telescope was determined to be −10.92 km/s using the NASA JPL Solar system dynamics HORIZONS web interface. The measured NIR spectrum of Venus is overlapped with the transmission window of FFP filter 18 to simulate the output lock-in amplifier signal. The phase offset is extracted and compared to the phase offsets using a pure CO.sub.2 spectrum over a range of radial velocities. The radial velocity at which the CO.sub.2 phase offset matches that of Venus spectrum indicating the approximate radial velocity of Venus. The radial velocities were determined with modulation ranges from 1-6 line spacings and shown in FIG. 9. If averaged over all harmonic multipliers, the radial velocity obtained is −10.52 km/s, which has an error of approximately 3.6% from the true value. The true radial velocity of Venus at the time of acquisition is shown as the horizontal line.

[0054] In another implementation, system 10 comprises a sensor comprising the lock-in amplifier and a multi-channel detector capable of any one of: monitoring remote gases; detecting a gas comprising quasi-periodic absorption lines; detecting a gas comprising complex absorption spectra; monitoring emissions from remote sources; and exoplanet detection and characterization.

[0055] The descriptions of the various embodiments of the present disclosure have been presented for purposes of illustration, but are not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein was chosen to best explain the principles of the embodiments, the practical application or technical improvement over technologies found in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.

[0056] Embodiments are described above with reference to block diagrams and/or operational illustrations of methods, systems. While the specification includes examples, the disclosure's scope is indicated by the following claims. Furthermore, while the specification has been described in language specific to structural features and/or methodological acts, the claims are not limited to the features or acts described above. Rather, the specific features and acts described above are disclosed as example for embodiments.

REFERENCES

[0057] 1. I. Coddington, N. Newbury, and W. Swann, “Dual-comb spectroscopy,” Optica 3, 414-426 (2016).

[0058] 2. I. Coddington, W. C. Swann, and N. R. Newbury, “Coherent dual-comb spectroscopy at high signal-to-noise ratio,” Physical Review A 82, 043817 (2010).

[0059] 3. A. Perot and C. Fabry, “On the Application of Interference Phenomena to the Solution of Various Problems of Spectroscopy and Metrology,” The Astrophysical Journal 9, 87 (1899).

[0060] 4. J. Kuhn, U. Platt, N. Bobrowski, and T. Wagner, “Towards imaging of atmospheric trace gases using Fabry—Perot interferometer correlation spectroscopy in the UV and visible spectral range,” Atmos. Meas. Tech. 12, 735-747 (2019).

[0061] 5. J. Kuhn, N. Bobrowski, P. Lübcke, L. Vogel, and U. Platt, “A Fabry—Perot interferometer-based camera for two-dimensional mapping of SO.sub.2 distributions,” Atmos. Meas. Tech. 7, 3705-3715 (2014).

[0062] 6. E. Vargas-Rodriguez and H. N. Rutt, “Design of CO, CO.sub.2 and CH4 gas sensors based on correlation spectroscopy using a Fabry—Perot interferometer,” Sensors and Actuators B: Chemical 137, 410-419 (2009).

[0063] 7. W. Ma, R. Wang, Q. Rong, Z. Shao, W. Zhang, T. Guo, J. Wang, and X. Qiao, “CO.sub.2 Gas Sensing Using Optical Fiber Fabry—Perot Interferometer Based on Polyethyleneimine/Poly(Vinyl Alcohol) Coating,” IEEE Photonics Journal 9, 1-8 (2017).

[0064] 8. Y. Jiang and W. Ding, “Recent developments in fiber optic spectral white-light interferometry,” Photonic Sensors 1, 62-71 (2011).

[0065] 9. L. Zhang, Y. Jiang, H. Gao, J. Jia, Y. Cui, W. Ma, S. Wang, and J. Hu, “A diaphragm-free fiber Fabry-Perot gas pressure sensor,” Review of Scientific Instruments 90, 025005 (2019).

[0066] 10. G. Costa, P. Gouvea, L. Soares, J. M. Pereira, F. Favero, A. Braga, P. Palffy-Muhoray, A. Bruno, and I. C. Carvalho, “In-fiber Fabry-Perot interferometer for strain and magnetic field sensing,” Optics Express 24, 14690 (2016).

[0067] 11. H. Li, D. Li, C. Xiong, W. Si, C. Fu, P. Yuan, and Y. Yu, “Low-Cost, High-Performance Fiber Optic Fabry—Perot Sensor for Ultrasonic Wave Detection,” Sensors (Basel) 19, 406 (2019).

[0068] 12. X. Yin, Y. Shen, D. Su, and Z. Shao, “High-spatial-resolution ultrasonic sensor using a fiber-optic Fabry—Perot interferometer,” Optics Communications 453, 124422 (2019).

[0069] 13. B. Zhou, H. Jiang, C. Lu, and S. He, “Hot Cavity Optical Fiber Fabry-Perot Interferometer as a Flow Sensor With Temperature Self-Calibrated,” Journal of Lightwave Technology 34, 5044-5048 (2016).

[0070] 14. S. Pevec and D. Donlagic, “Miniature fiber-optic Fabry-Perot refractive index sensor for gas sensing with a resolution of 5×10-9 RIU,” Opt. Express 26, 23868-23882 (2018).

[0071] 15. S. Wu, L. Wang, X. Chen, and B. Zhou, “Flexible Optical Fiber Fabry—Perot Interferometer Based Acoustic and Mechanical Vibration Sensor,” Journal of Lightwave Technology 36, 2216-2221 (2018).

[0072] 16. A. Javernik and D. Donlagic, “Miniature, micro-machined, fiber-optic Fabry-Perot voltage sensor,” Opt. Express 27, 13280-13291 (2019).

[0073] 17. M. Mozafari, Z. Sanjabi, and N. Granpayeh, “Design and fabrication of piezo-optical Fabry-Perot voltage sensor,” in 2008 IEEE Avionics, Fiber-Optics and Photonics Technology Conference (2008), pp. 39-40.

[0074] 18. B. Zhou, H. Jiang, R. Wang, and C. Lu, “Optical Fiber Fiber Fabry—Perot Filter With Tunable Cavity for High-Precision Resonance Wavelength Adjustment,” J. Lightwave Technol. 33, 2950-2954 (2015).

[0075] 19. D. Hays, A. Zribi, S. Chandrasekaran, S. Goravar, S. Maity, L. R. Douglas, K. Hsu, and A. Banerjee, “A Hybrid MEMS—Fiber Optic Tunable Fabry—Perot Filter,” Journal of Microelectromechanical Systems 19, 419-429 (2010).

[0076] 20. D. Jáuregui-Vázquez, J. M. Estudillo-Ayala, R. Rojas-Laguna, E. Vargas-Rodriguez, J. M. Sierra-Hernandez, J. C. Hernández-Garcia, and R. I. Mata-Chavez, “An all fiber intrinsic Fabry-Perot Interferometer based on an air-microcavity,” Sensors (Basel) 13, 6355-6364 (2013).

[0077] 21. N. Ismail, C. C. Kores, D. Geskus, and M. Pollnau, “Fabry-Pérot resonator: spectral line shapes, generic and related Airy distributions, linewidths, finesses, and performance at low or frequency-dependent reflectivity,” Opt. Express 24, 16366-16389 (2016).

[0078] 22. S. Abrarov, A Rational Approximation of the Voigt Function (n.d.).