Vol. 25, No. 18 | 4 Sep 2017 | OPTICS EXPRESS 21094
Phase demodulation of interferometric fiber sensor based on fast Fourier analysis XIN FU,1 PING LU,1,* WENJUN NI,1 HAO LIAO,1 DEMING LIU,1 AND JIANGSHAN ZHANG2,3 1Wuhan
National Laboratory for Optoelectronics (WNLO) and National Engineering Laboratory for Next Generation Internet Access System, School of Optics and Electronic Information, Huazhong University of Science and Technology, Wuhan 430074, China 2Department of Electronics and Information Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
[email protected] *
[email protected]
Abstract: A demodulation method for interferometric fiber sensors (IFSs) is proposed in this article. The phase variation induced by the measurands can be estimated by calculating the Fourier phase at the intrinsic spatial frequencies of the fiber sensor. Theoretical analysis of the demodulation method is discussed in detail. Numerical simulations are put forward to demonstrate the consistency of the demodulation results under different wavelength sampling interval and noise level, showing a better stability compared with the conventional peak wavelength tracking technique. The proposed method is also experimentally demonstrated by an inline multimode interferometer based on a single-mode fiber (SMF) offset-splicing structure. Experimental results indicate that the phase response of different cladding modes can be analyzed simultaneously. Simultaneous measurement of strain and temperature is realized in our confirmatory experiment by analyzing the phase sensitivities of two selected cladding modes. © 2017 Optical Society of America OCIS codes: (060.2310) Fiber optics; (060.2370) Fiber optics sensors; (280.4788) Optical sensing and sensors.
References and links 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13.
D. Kreit, R. C. Youngquist, and D. E. N. Davies, “Two-mode fiber interferometer/amplitude modulator,” Appl. Opt. 25(23), 4433–4438 (1986). T. Okamoto and I. Yamaguchi, “Multimode fiber-optic Mach-Zehnder interferometer and its use in temperature measurement,” Appl. Opt. 27(15), 3085–3087 (1988). M. Imai, T. Ohashi, and Y. Ohtsuka, “Fiber-optic Michelson interferometer using an optical power divider,” Opt. Lett. 5(10), 418–420 (1980). H. Liu, D. Wang, J. Liu, and S. Liu, “Range tunable optical fiber micro-Fabry–Pérot interferometer for pressure sensing,” IEEE Photonics Technol. Lett. 28(4), 402–405 (2016). M. Hou, F. Zhu, Y. Wang, Y. Wang, C. Liao, S. Liu, and P. Lu, “Antiresonant reflecting guidance mechanism in hollow-core fiber for gas pressure sensing,” Opt. Express 24(24), 27890–27898 (2016). L. Cheng, C. Wang, Y. Huang, H. Liang, and B. O. Guan, “Silk fibroin diaphragm-based fiber-tip Fabry-Perot pressure sensor,” Opt. Express 24(17), 19600–19606 (2016). X. Zhang and W. Peng, “Bent fiber interferometer,” J. Lightwave Technol. 33(15), 3351–3356 (2015). Y. Zhao, F. Xia, and J. Li, “Sensitivity-enhanced photonic crystal fiber refractive index sensor with two waistbroadened tapers,” J. Lightwave Technol. 34(4), 1373–1379 (2016). J. Tian, Z. Lu, M. Quan, Y. Jiao, and Y. Yao, “Fast response Fabry-Perot interferometer microfluidic refractive index fiber sensor based on concave-core photonic crystal fiber,” Opt. Express 24(18), 20132–20142 (2016). X. Ni, S. Fu, and Z. Zhao, “Thin-fiber-based Fabry–Pérot cavity for monitoring microfluidic refractive index,” IEEE Photonics J. 8(3), 1–7 (2016). C. Han, H. Ding, X. Li, and S. Dong, “Temperature insensitive refractive index sensor based on single-mode micro-fiber Sagnac loop interferometer,” Appl. Phys. Lett. 104(18), 181906 (2014). J. N. Dash, R. Jha, J. Villatoro, and S. Dass, “Nano-displacement sensor based on photonic crystal fiber modal interferometer,” Opt. Lett. 40(4), 467–470 (2015). J. Su, X. Dong, and C. Lu, “Property of bent few-mode fiber and its application in displacement sensor,” IEEE Photonics Technol. Lett. 28(13), 1387–1390 (2016).
#301288 Journal © 2017
https://doi.org/10.1364/OE.25.021094 Received 29 Jun 2017; revised 15 Aug 2017; accepted 17 Aug 2017; published 21 Aug 2017
Vol. 25, No. 18 | 4 Sep 2017 | OPTICS EXPRESS 21095
14. G. K. B. Costa, P. M. P. Gouvêa, L. M. B. Soares, J. M. B. Pereira, F. Favero, A. M. B. Braga, P. PalffyMuhoray, A. C. Bruno, and I. C. S. Carvalho, “In-fiber Fabry-Perot interferometer for strain and magnetic field sensing,” Opt. Express 24(13), 14690–14696 (2016). 15. T. Geng, J. He, W. Yang, X. Chen, M. An, C. Sun, X. Jin, and L. Yuan, “Modal interferometer using three-core fiber for simultaneous measurement strain and temperature,” IEEE Photonics J. 8(4), 1–8 (2016). 16. S. Liu, K. Yang, Y. Wang, J. Qu, C. Liao, J. He, Z. Li, G. Yin, B. Sun, J. Zhou, G. Wang, J. Tang, and J. Zhao, “High-sensitivity strain sensor based on in-fiber rectangular air bubble,” Sci. Rep. 5(1), 7624 (2015). 17. X. Dong, H. Y. Tam, and P. Shum, “Temperature-insensitive strain sensor with polarization-maintaining photonic crystal fiber based Sagnac interferometer,” Appl. Phys. Lett. 90(15), 151113 (2007). 18. S. Wang, P. Lu, L. Mao, D. Liu, and S. Jiang, “Cascaded interferometers structure based on dual-pass MachZehnder interferometer and Sagnac interferometer for dual-parameter sensing,” Opt. Express 23(2), 674–680 (2015). 19. K. Tian, Y. Xin, W. Yang, T. Geng, J. Ren, Y.-X. Fan, G. Farrell, E. Lewis, and P. Wang, “A curvature sensor based on twisted single-mode–multimode–single-mode hybrid optical fiber structure,” J. Lightwave Technol. 35(9), 1725–1731 (2017). 20. Z. Wu, P. P. Shum, X. Shao, H. Zhang, N. Zhang, T. Huang, G. Humbert, J. L. Auguste, F. Gérome, J. M. Blondy, and X. Q. Dinh, “Temperature- and strain-insensitive curvature sensor based on ring-core modes in dual-concentric-core fiber,” Opt. Lett. 41(2), 380–383 (2016). 21. G. Salceda-Delgado, A. Van Newkirk, J. E. Antonio-Lopez, A. Martinez-Rios, A. Schülzgen, and R. Amezcua Correa, “Compact fiber-optic curvature sensor based on super-mode interference in a seven-core fiber,” Opt. Lett. 40(7), 1468–1471 (2015). 22. L. Gounaridis, P. Groumas, E. Schreuder, R. Heideman, H. Avramopoulos, and C. Kouloumentas, “New set of design rules for resonant refractive index sensors enabled by FFT based processing of the measurement data,” Opt. Express 24(7), 7611–7632 (2016). 23. L. Gounaridis, P. Groumas, E. Schreuder, G. Tsekenis, A. Marousis, R. Heideman, H. Avramopoulos, and C. Kouloumentas, “High performance refractive index sensor based on low Q-factor ring resonators and FFT processing of wavelength scanning data,” Opt. Express 25(7), 7483–7495 (2017). 24. B. Dong, Y. Peng, Y. Wang, and C. Yu, “Mode division multiplexing in a fiber modal interferometer for dualparameters measurement,” IEEE Photonics Technol. Lett. 28(2), 143–146 (2016). 25. Y. Sun, D. Liu, P. Lu, Q. Sun, W. Yang, S. Wang, L. Liu, and J. Zhang, “Dual-parameters optical fiber sensor with enhanced resolution using twisted MMF based on SMS structure,” IEEE Sens. J. 17(10), 3045–3051 (2017). 26. C. Gouveia, P. A. S. Jorge, J. M. Baptista, and O. Frazao, “Temperature-independent curvature sensor using FBG cladding modes based on a core misaligned splice,” IEEE Photonics Technol. Lett. 23(12), 804–806 (2011). 27. X. Yu, D. Bu, X. Chen, J. Zhang, and S. Liu, “Lateral stress sensor based on an in-fiber Mach–Zehnder interferometer and Fourier analysis,” IEEE Photonics J. 8(2), 1–10 (2016).
1. Introduction Interferometric fiber sensors (IFSs) are widely investigated due to the remarkable superiorities of high sensitivity and precision. The interferences are formed by two or more beams that have certain phase difference [1–3]. Since the optical wavelength is very short, usually several hundreds of nanometers or a few micrometers, very small optical path difference (OPD) variation can be detected by IFSs, offering relatively high sensitivity and resolution in measurement of parameters such as pressure [4–6], refractive index (RI) [7–11], displacement [12,13], strain [14–18], and curvature [19–21], etc. Wavelength tracking demodulation technique is mostly used for interferometric fiber sensors due to the simple operation and linear response. The mechanism is based on the resonant wavelength shift due to the OPD variation caused by physical parameters. The wavelength shift usually shows a linear response to the measurands. By calibrating the sensitivity, the measurand can be interrogated from the wavelength shift. However, a significant limitation of wavelength demodulation is the accuracy. Since the resonant wavelength is determined by searching the maximum (for resonant peak) or minimum (for resonant dip) intensity value in a certain wavelength window, the results may suffer from unpredictable errors that depend on wavelength sampling rate or spectral noise [22,23]. In order to improve the accuracy, Fast Fourier Transform (FFT) is applied to the optical spectrum to demodulate the sensor signal in the frequency spectrum domain [22–25]. Traditional FFT method to demodulate the fiber sensor is to monitor the peak variations (including both peak intensity variation and peak frequency variation) in the FFT spectrum. It
Vol. 25, No. 18 | 4 Sep 2017 | OPTICS EXPRESS 21096
usually requires a relatively large amount of variation of the measurands. But in some sensing applications, parameters that have very small amount of changes should be detected, such as micro pressure or strain sensing. In these applications, the OPD variation introduced by the measurands is small enough to be neglected compared with the initial OPD value. Therefore, the peak location variation in the FFT spectrum is not able to be recognized. In this article, we propose a demodulation method to calculate the phase variation of the sensor in the Fourier domain. For fiber interferometer, most of the energy is contained in the intrinsic spatial frequencies that are corresponding to the interference components. For this reason, we believe that the information of phase variation between each interference component can be obtained by the Fourier phase at these spatial frequencies. The mechanism is theoretically analyzed. Numerical simulations prove that the demodulation result is insensitive to the spectral noise and sampling interval. Experiment is conducted by testing multimode interferometers to demonstrate the proposed scheme. For multimode interference, phase response of each cladding mode can be obtained by calculating the Fourier phase at each related spatial frequency, so it can be applied to multi-parameters measurement or quasidistributed sensing system. 2. Theoretical analysis To theoretically explain the principle of the demodulation method, we firstly use a model of two-beam interference for analysis. The condition of multimode interference is generalized theoretically and demonstrated in the subsequent experiment. For a two-beam interferometer, the spectrum function can be expressed by Eq. (1). The interference function can be regarded as a trigonometric function of the wavelength within a certain range. In the equation, d represents the optical path difference (OPD). A and B are two constants determined by the optical intensity and coupling ratio between two optical paths. T (λ ) = A + B cos(
2π
(1) ⋅ d ) ≈ A + B cos( βλ ) λ When the OPD of the interferometer is modulated by the environmental parameters, the whole spectrum will show a wavelength shift, as written by Eqs. (2) and (3). d is the OPD variation, while λ is the wavelength shift value. T '(λ )= A + B cos 2π ⋅(d −∆d )= T (λ + ∆λ )
= ∆λ
λ
λ ⋅ ∆d
≈
(2)
λ0 ⋅ ∆d
(3) d − ∆d d − ∆d From Eq. (3), the shift value is dependent on the wavelength. Within a certain wavelength window that is not too wide, the shift can be regarded as a constant. The free spectrum range (FSR) variation is approximately presented by Eq. (4). Usually the OPD variation is too small compared with the initial OPD value, thus the FSR variation is small enough to be neglected. According to this, we can conclude that the spectrum experiences a constant wavelength shift without shape distortion, as demonstrated by Eqs. (5) and (6), and φ is the additional phase variation item introduced by the wavelength shift. It can be seen from the equation that the wavelength is discretely sampled, just like the spectrum data acquired for analysis in the experiment or practical applications. N and λ0 stand for the number and interval of wavelength sampling points. = ∆FSR λ0 2 (
1 1 ∆d − )≈ ⋅ FSR