Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Temperature-compensated multi-point refractive index sensing based on a cascaded Fabry-Perot cavity and FMCW interferometry

Open Access Open Access

Abstract

We proposed a novel temperature-compensated multi-point refractive index (RI) sensing system by the combination of the cascaded Fabry-Perot (FP) sensors and the frequency modulated continuous wave (FMCW) interferometry. The former is used for simultaneous sensing of RI and temperature, and the latter is used for multiplexing a series of the cascaded FP sensors to realize multi-point sensing. By means of Fourier transform-based algorithms, the interference spectra of each sub-FP sensors can be divided and demodulated independently. Experimentally, three cascaded FP sensors are multiplexed to verify multi-point RI and temperature sensing ability. RI sensitivity up to ∼1200 nm/RIU is obtained within RI range from 1.3330 to 1.3410, and temperature sensitivity up to ∼0.17 nm/°C is obtained within temperature range from 20 °C to 80 °C. The RI precision is as high as 10−5 RIU and the temperature precision is as high as 0.05 °C. In addition, the prospective multiplexing number could reach about 4000 estimated by the minimum detectable light power. The proposed sensing system has potential advantages in the practical applications that require a large number sensing points.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Refractive index (RI) measurement plays an important role in the fields of biomedicine, chemistry and environmental monitoring, etc. [1]. In the last two decades, the optical fiber RI sensors attract an increasing attention due to their advantages such as immunity to electromagnetic interference, high sensitivity, compact structure, and online monitoring capability. Many types of optical fiber RI sensors have been proposed, including optical fiber interferometer (OFI) [2], optical fiber grating (OFG) [3], surface plasmon resonance (SPR) [4], whispering gallery mode (WGM) [5], and so on. However, the RI measurement result is always affected by temperature disturbances due to the temperature dependence of RI. With the increasing demand for detection accuracy, it is realized that temperature compensation is essential in RI measurement. So far, many optical fiber sensors for simultaneous measuring RI and temperature have been reported by fabricating various structures. Examples include cascaded fiber Bragg gratings (FBG) [68], cascaded Fabry-Perot interferometers (FPI) [911], cascaded Mach-Zehnder interferometers (MZI) [1214], dual-channel SPR sensors [1517], and arbitrary combination of various fiber structures [1823]. Generally, most of these sensors requires a wider spectrum span that allows two sub-spectra responded RI and temperature respectively without cross-talk. Exceptionally, the cascaded interferometer-based sensors have low demand of the wider spectrum span. The traditional cascaded interferometer-based sensors are demodulated by the spatial frequency multiplexing method, the sub-spectra can be extracted and separated by fast Fourier transform (FFT) filter. However, it requires a large difference of optical path difference (OPD) between sub-interferometers, which results in a limited multiplexing number [911,20,24,25].

Reviewing the development of the optical fiber sensor for simultaneous measuring RI and temperature, many researchers have made significant contributions [623]. However, all of the above optical fiber sensors can only realize single point sensing. For multi-point sensing applications, these sensors are powerless. Up to now, there are few reports on the simultaneous measurement of multi-point RI and temperature. Thus, it is urgent to propose a feasible scheme to reduce the sensing cost. Fortunately, frequency-modulated-continuous-wave (FMCW) technique originated from Radar provides us helpful inspiration. The FMCW technique has been developed rapidly since it was applied to optical fiber diagnosis in 1981 [26], also named optical frequency domain reflectometer (OFDR). Afterwards, distributed temperature and strain sensing is also realized thanks to its ability of simultaneous localization and sensing [27]. Nowadays, OFDR has become a common technique for distributed optical fiber sensing, and is widely used in the sensing of temperature, strain, vibration, and shape [2831], and is also used to multiplex FBG sensor [32]. However, all of the above-mentioned sensing applications need rely on the Rayleigh backscattering (RBS) in the optical fiber for signal demodulation. The RBS is insensitive to the external RI, so the RI sensing is difficult to be achieved by the OFDR system. Recent years, some successful cases for multi-point RI sensing have been demonstrated by combining OFDR with the etched fiber [33], the macrobending fiber [34], and the tapered fiber [35]. However, the sensing point is limited the fiber length or power loss thereby difficult to reach satisfactory number. Moreover, these multi-point RI sensing systems are low-sensitive to sample RI and have not temperature compensation ability.

In this paper, we propose a novel multi-point RI sensing system with temperature compensation function. The sensing system adopts the FMCW-based interferometer, where a tunable laser source (TLS) whose wavelength in linearly sweeping is employed to replace the traditional white light source, and a photodetector is employed as receiver to replace the traditional optical spectrum analyzer. Meanwhile, the cascaded FP sensor with an open cavity is designed and fabricated for simultaneous measurement of RI and temperature. With help of the FMCW interferometry, the OPD of the sub-FP sensors is no longer limited. In other words, the OPD of the sub-FP sensors can either be the same or different, even all of the sub-FP sensors can have the same OPD, which will not only greatly reduce the manufacturing difficulty of the cascaded FP sensors, but also greatly increase the multiplexing number. Three pairs of cascaded FP sensors are connected into the sensing system in parallel for verifying the multi-point sensing capability. By means of FFT algorithm, we can distinguish each-FP sensor in the spatial domain and separate their reflection signals. By means of inverse FFT (IFFT) algorithm, the interference spectra of each sub-FP sensor can be demodulated independently in wavelength domain. Experimentally, we test the abilities of multi-point simultaneous RI and temperature sensing under RI range from 1.3330 RIU to 1.3410 RIU, and temperature range from 20 °C to 80 °C. RI sensitivities up to ∼1200 nm/RIU and temperature sensitivities up to ∼0.17 nm/°C are obtained experimentally. In addition, the prospective multiplexing number is estimated according to the minimum detectable power, reaching about 4000. The proposed sensing system not only inherits advantage of the simultaneous location and sensing capabilities of the traditional distributed sensors, but also overcomes the difficulty of the traditional single-point sensors in multi-point sensing, which is of great significance in practical applications that require a large number sensing points.

2. Multiplexing principle and simulation

2.1 Basic configuration of the multiplexing system

The basic configuration of the proposed FMCW-based multiplexing system for multi-point simultaneous measuring RI and temperature is based on an optical fiber MZ interferometer, as shown in Fig. 1(a). A TLS whose wavelength (or optical frequency) in linearly sweeping is split into reference and measurement arms. In the measurement arm, a series of cascaded FP sensors are connected in parallel by optical couplers. Each cascaded FP sensor consists of two sub-FP sensors, one of which is used for temperature sensing, and the other is used for RI sensing. The structure of the cascaded FP sensors for simultaneous measuring RI and temperature is shown in Fig. 1(b). The sub-FP sensors for RI sensing have open cavity allowing the entry of sample under test, thus they are sensitive to RI and temperature. The sub-FP sensors for temperature sensing have close cavity, thus they are only sensitive to temperature. Due to the different time delay, the reflected light from different sub-sensors would generate beat signals of different beat frequencies. The superposition result of these beat signals is received by a balanced photodetector (BPD) and a data acquisition card (DAQ). The polarization controller (PC) in the reference arm is used for maximizing the beat signal intensity. By means of FFT and IFFT algorithms, the sub-FP sensors can be distinguished in the beat frequency domain, and their interference spectra can be demodulated in the wavelength domain. The detail data processing is described in the next section.

 figure: Fig. 1.

Fig. 1. (a) Basic configuration of the proposed FMCW-based multiplexing method for multi-point simultaneous measuring RI and temperature. A series of cascaded FP sensors are connected in parallel by optical couplers. (b) The structure of the cascaded FP sensor. (c)-(e) The signal demodulation processing of the multiplexing system.

Download Full Size | PDF

2.2 Signal demodulation process

The signal demodulation processing for the proposed multiplexing method is illustrated in Figs. 1(c)–1(e). The reflected light from different sub-FP sensors have different time delays (τ) relative to the reference light, therefore, the beat signals of different frequencies are generated, as shown in Fig. 1(c). These beat signals are mixed in time domain (or wavelength domain). Applying FFT algorithm, their power spectra can be separated in frequency domain (or spatial domain), as shown in Fig. 1(d). Then, applying IFFT algorithm, the peak of interest is extracted and transformed back to time (wavelength) domain, thereby the interference spectrum of each sub-FP sensor is demodulated, as shown in Fig. 1(e). Noticeably, the time domain and wavelength domain are equivalent according to formula: Δλ=c/(γt), where λ donates wavelength, c is the speed of light in vacuum, γ donates tuning speed of the TLS in frequency unit, and t is time. The frequency domain and spatial domain are also equivalent according to formula: 2L = cfb/γ, where 2L donates optical path difference (OPD) of sub-FP sensor in measurement arm relative to reference arm, fb is beat frequency.

2.3 Simulation analysis

The theoretical model is established according to the principle of coherence for numerical simulation. The electric field through the MZ interferometer can be expressed by

$${E_{total}}(t) = {E_R}(t) + \sum\limits_{m = 1}^M {{r_m}} {E_{Pm}}({t - {\tau_m}} )$$
where ER(t) represents the electric field through the reference arm, EPm(t-τm) represents the electric field of the mth reflected point, rm represents reflectivity of the mth reflected point and τm represents time delay, τm=2Lm/c. The electric field ER(t) can be expressed by
$${E_R} = {E_{R0}}{e^{j\varphi (t)}}$$
where ER0 represents electric field amplitude, φ(t) represents electric field phase. For the linear frequency tuning light source, it can be expressed by
$$\varphi (t) = 2\pi {f_0}t + \pi \gamma {t^2}$$
where f0 represents initial optical frequency of the light source, f0=c/λ0. Replacing t with t-τm, the electric field EPm(t-τm) can also be calculated. Finally, the beat signal intensity can be calculated by
$${I_{total}} = {E_{total}}E_{total}^\ast $$

According to Eqs. (1)–(4), the original beat signal can be simulated. As an example, three cascaded FP sensors numbered as #1, #2 and #3 are adopt for simulation. Each cascaded FP sensor consists of two sub-FP sensors, one of which is used as temperature sensing channel, and the other is used as RI sensing channel. Besides, each sub-FP sensor has two reflected end-faces, a total of twelve for the whole system. Due to transmission loss, the reflectivity of the second end-face is smaller than that of the first end-face for each sub-FP sensor. These six sub-FP sensors are distributed in different locations of the multiplexing system and their cavities are filled with samples of different RIs. The detailed simulation parameters are shown in Table 1. The wavelength tuning range is from 1525 nm to 1575 nm (f0 = 190.5 THz) with tuning speed of 80 nm/s (γ = 10 THz/s). The sampling time is 0.625 s with sampling rate of 10 MSa/s.

Tables Icon

Table 1. Simulation parameters of the multiplexing system

The simulation result of the original beat signal is shown in Fig. 2(a), the data length is 6.25 MSa. Then, the power spectrum can be calculated by FFT, as shown in Fig. 2(b). Six pairs of reflected peaks appear at the preset location clearly. As same as the traditional OFDR sensing system [27], the spatial resolution can be given by Δz = c/(2ΔF), where ΔF is frequency sweeping span in Hz, here Δz = 20 µm. The close view of the reflected peaks generated by temperature sensing channels FP1, FP3 and FP5 are shown in Fig. 2(c), and that generated by RI sensing channels FP2, FP4 and FP6 are shown in Fig. 2(e). For each sub-FP sensor, the reflected peaks generated by the first end-face and the second end-face are both clearly visible, and their intensity ratios are almost equal to the preset value. Then, the reflected peaks of interest are selected by a Hanning window with width of 40 points and transformed back to wavelength domain by IFFT, thus the interference spectrum of the FP sensor of interest can be demodulated. Figure 2(d) shows the interference spectra of the temperature sensing channel FP1, FP3 and FP5, the cavity parameters of these three FP sensors are exactly consistent, resulting their interference spectra to completely overlap. Figure 2(f) shows the interference spectra of the RI sensing channel FP2, FP4 and FP6, the RIs in cavity of these three FP sensors are different, resulting their interference spectra to shift. In addition, the simulation result of the RI sensitivity is 1205 nm/RIU, which is almost consistent with other researches of traditional FP sensors [2,911].

 figure: Fig. 2.

Fig. 2. (a) Simulation result of the original beat signal. (b) Power spectrum of the beat signal after applying FFT. The reflected peaks appear at the preset location. (c) Close view of the reflected peaks generated by temperature sensing channels FP1, FP3 and FP5, respectively. (d) The interference spectra of the temperature sensing channels after applying IFFT. (e) Close view of the reflected peaks generated by RI sensing channels FP2, FP4 and FP6, respectively. (f) The interference spectra of the RI sensing channels after applying IFFT.

Download Full Size | PDF

2.4 Effect of multiplexing number on wavelength accuracy

To investigate the effect of the multiplexing number of the cascaded FP sensors on the wavelength accuracy of an individual interference spectrum, we perform the simulation as follows. The distance between the Temperature sensing channel and the RI sensing channel is 1 mm, the distance between the adjacent cascaded FP sensors is 10 mm. Other simulation conditions are consistent with Section 2.3. Based on the wavelength tuning speed and sampling rate, the maximum sensing distance is calculated to 50 m, thus the maximum multiplexing number is 5000. We simulate the interference spectrum of the RI sensing channel of the first cascaded FP sensor when the multiplexing number are 1, 10, 100, 1000, and 5000 respectively. The simulation results are shown in Fig. 3. The interference spectrum keeps almost consistent at different multiplexing number. When the multiplexing number is 5000, the wavelength error is only 7 pm. The simulation results indicate that the accuracy of an individual sensor element is not significantly impaired even with a large number of superimposed sensors.

 figure: Fig. 3.

Fig. 3. The effect of the multiplexing number on the wavelength error of an individual interference spectrum. Inset: the interference spectrum of an individual sub-FP sensor when multiplexing number are 1, 10, 100, 1000, and 5000 respectively.

Download Full Size | PDF

3. Experiment preparation

3.1 Manufacture of cascaded fiber FP sensor

The cascaded FP sensors for experiment are self-made by bonding single-mode fiber (SMF) and glass tube together using UV glue. Simply, the manufacturing operation includes four steps. Firstly, for the temperature sensing channel, an SMF with smooth end-face is put into a glass tube, whose inner diameter is 126 µm, slightly larger than the fiber diameter of 125 µm. A drop of UV glue is put on the contact point between the SMF and the glass tube, and it then is cured for 3 minutes. Secondly, the SMF is connected to the traditional optical fiber sensing system consisting of white light source and optical spectrum analyzer (AQ6370C, Yokogawa). Another SMF is put into the glass tube, the gap between two end-faces, i.e. the cavity length, is adjusted by a micro displacement stage of 10 µm accuracy. Meanwhile, the reflected spectrum is real-time monitored until occurring suitable interference spectrum. Thirdly, the contact point is bonded by UV glue as same as the first step. The optical fiber structure is placed under UV light for a half hour until UV glue completely cured. Similarly, another glass tube with a side hole, which allows sample to flow in and out, is used to make RI sensing channel. Finally, the temperature and RI sensing channels are spliced by optical fiber fusion splicer. This manufacture method is simpler than the femtosecond laser micromachining and the special fiber micromachining [2,911]. In addition, due to the larger thermal expansion coefficient of the UV glue, the FP sensor structure will subject more obvious thermal deformation, which can be used as sensitivity enhanced temperature sensors [36]. It is noting that our sensor structure will become unstable after long term use because the UV glue is susceptible to damage. In the practical sensing applications, using the femtosecond laser micromachining is more reliable.

3.2 Experiment set-up of the multiplexing system

Figure 4 shows experiment set-up of the proposed multiplexing system, which includes three parts. Part I: Auxiliary interferometer. Generally, it is difficult to realize the ideal linear frequency sweep with the TLS, the nonlinearity will degrade the reflection peak so that it cannot be distinguished [37]. In order to compensate the nonlinearity noise, the common way is to add an auxiliary interferometer, which can realize sampling with constant frequency interval [38]. Here, we adopt a Michelson interferometer as the auxiliary interferometer, two Faraday rotating mirrors (FRM) are connected terminally to provide end reflection and eliminate polarization-induced fading [35], and a delay fiber of 100 m provides optical path difference. A clock generator is employed to convert the sinusoidal beat signal into TTL clock signal, which is used as the sampling clock of the DAQ. Part II: Absorption gas chamber. Due to ambient temperature disturbance, the wavelength of the TLS often tends to shift, which will affect the accuracy of the sensing spectrum. The absorption gas chamber can used to calibrate wavelength shift because its absorption wavelength is insensitive to temperature change. Here, we adopt acetylene gas, which has a series of absorption peaks around 1510 nm ∼ 1540 nm, to calibrate the wavelength drift of the TLS. Part III: Main interferometer. Three pairs of cascaded FP sensors are connected into the multiplexing system by three optical coupler. The detailed description of this part can be seen in section 2.1. The coupling ratio of each optical coupler is descripted in Fig. 4(a). Figure 4(b) shows the photograph of the cascaded FP sensor, limited by our equipment and fabrication method, there is an addition fiber of about 20 cm between the RI sensing channel and the temperature sensing channel. For temperature compensation, during experiment, the RI sensing channel and the temperature sensing channel are placed closely together by bending the fiber. In addition, the equipment information is as follows: the TLS is from Agilent Technologies (Agilent 81600B), the wavelength sweeping span is from 1520 nm to 1580 nm, and the wavelength sweeping speed is 80 nm/s. The output power of the TLS is 5 mW. The BPD is from Fsphotonics (PDB1008), whose bandwidth is 80 MHz. The DAQ is from ART technology (PCI8514). The sampling length is 8 MSa and the sampling time is 787 ms, so the wavelength sampling span is about 60 nm, corresponding to spatial resolution of 20 µm. Besides, according to Nyquist theory, the maximum measurable length of the main interferometer is 50 m, which can be improved by extending the delay fiber length of the auxiliary interferometer.

 figure: Fig. 4.

Fig. 4. (a) The sketch diagram of the proposed FMCW-based FP sensors multiplexing system; (b) photograph of the cascaded FP sensors for simultaneous RI and temperature sensing. TLS: tunable laser source; FRM: Faraday rotating mirror; PC: polarization controller; CG: clock generator; BPD: balanced photodetector; DAQ: data acquisition card; SMF: single mode fiber.

Download Full Size | PDF

In order to execute RI and temperature sensing experiments, a series of liquid samples of different RIs are pre-prepared by mixing glycerol-aqueous solutions with different concentrations. The sample RI is measured and calibrated by Abbe refractometer at 20 °C. Here, the sample RI is from 1.3330 RIU to 1.3410 RIU with the interval of 0.0010 RIU. During the experiments, the cascaded FP sensors (both RI sensing channel and temperature sensing channel) are immersed into the sample under test, meanwhile, the sample is placed into water bath to control the ambient temperature accurately.

4. Experiment results and discussions

The signal demodulation process is the same as described in the previous simulation. The original beat signal by experiment is shown in Fig. 5(a), which is messy due to the mixture of countless beat signals. The power spectrum of the original beat signal is then calculated by applying FFT, as shown in Fig. 5(b). Distinctly, there appears several peaks in the power spectrum. The peaks marked green are generated by the connectors in the optical fiber link, and the peaks marked red are generated by six sub-FP sensors. The close view of the latter (taking FP1 as example) is shown in Fig. 5(c), where two reflected peaks marked R1 and R2, which is generated by two end-face of the FP cavity, is clearly visible. The intensity of the peak R2 is smaller than that of the peak R1 due to the transmission loss of the FP structure. Then, applying IFFT to the selected peaks and the original interference spectrum of the sub-FP sensor can be demodulated, as shown in Fig. 5(d). However, different from the simulation results, the original interference spectrum is irregular in intensity due to polarization and power fluctuation of the light source, hence we employ a baseline filter to normalize the spectrum. Finally, the normalized interference spectrum of the sub-FP sensor is obtained, as shown in Fig. 5(e).

 figure: Fig. 5.

Fig. 5. (a) The original experimental beat signal; (b) power spectrum of the beat signal; (c) close view of the reflection peaks generated by FP1; (d) the original interference spectrum of the FP1 obtained after applying the IFFT process; (e) the normalized interference spectrum of the FP1 obtained after removing the baseline.

Download Full Size | PDF

Following the above signal demodulation process, the interference spectra of all sub-FP sensors can be demodulated. The close view of their reflected peaks and the corresponding interference spectra are shown in Figs. 6(a)–6(b). It can be seen that the interference spectra of every sub-FP sensors are clearly visible. The dashed lines in the Fig. 6(b) represent the normalized interference spectra measured by the traditional optical spectrum analyzer (AQ6370C). Two curves are almost consistent, verifying the spectral reliability of the FMCW-based multiplexing method. Table 2 shows the location of these six sub-FP sensors and their OPDs. The OPDs are almost the same with each other, verifying that the OPD of the cascaded FP sensors is no longer limited compared to the traditional spatial frequency multiplexing method [911,24,25].

 figure: Fig. 6.

Fig. 6. (a) Close view of the reflection peaks generated by six FP sensors from FP1 to FP6, respectively; (b) the normalized interference spectra (solid lines) of six FP sensors from FP1 to FP6 respectively, which are almost consistent with the measurement results by the traditional optical spectrum analyzer (dashed lines).

Download Full Size | PDF

Tables Icon

Table 2. Experiment parameters of the multiplexed sub-FP sensors

Subsequently, the experiments of simultaneous measurement of RI and temperature are carried out. Firstly, the cascaded FP sensors are successively immersed into pre-prepared samples with different RIs. In order to ensure the accuracy of the experiment results, these sensors are cleaned and dried every time when the sample is changed. Meanwhile, it is especially necessary to keep the temperature constant during the experiment. Figure 7 shows the experiment results under a series of RIs ranging from 1.3330 RIU to 1.3410 RIU with the interval of 0.0010 RIU and the constant temperature of 20 °C. The spectrum evolution (taking the cascaded FP sensor #1 as example) is shown in Figs. 7(a)–7(b). Here, only one spectrum period are plotted out for clarity. The interference spectrum of the FP1 (temperature sensing channel) remains unchanged with RI of sample because the FP cavity is isolated from the external environment. The interference spectrum of FP2 (RI sensing channel) red-shifts with increasing of RI of sample, which is consistent with the research results of traditional FP-based RI sensors. The sensing curves of the total six sub-FP sensors are shown in Fig. 7(c). The sensitivities of three sub-FP sensors as RI sensing channel are 1203.9 nm/RIU, 1279.8 nm/RIU and 1230.3 nm/RIU, respectively, and the other three sub-FP sensors as temperature sensing are all insensitive to RI change. The experiment results also verify there is no cross-talk between the sub-FP sensors.

 figure: Fig. 7.

Fig. 7. The interference spectra shift of (a) the temperature sensing channel FP1 and (b) the RI sensing channel FP2 as a function of RI of the sample under test, from 1.3330 RIU to 1.3410 RIU with the interval of 0.0010 RIU; (c) the RI sensing curves of six FP sensors from FP1 to FP6, correspondingly.

Download Full Size | PDF

Next, the cascaded FP sensors are immersed into one sample that keeps unchanged, and the ambient temperature is accurately controlled by the water bath. Figure 8 shows the experiment results under a series of temperature ranging from 20 °C to 80 °C with the interval of 10 °C. The spectrum evolution (taking the cascaded FP sensor #1 as example) is shown in Fig. 8(a) and Fig. 8(b). The interference spectra of both FP1 and FP2 red-shift with increasing of temperature. The spectrum shift is caused by the combination of the thermal refractive index change of the sample and the thermal expansion of the FP cavity. Generally, the thermo-optical coefficient of the sample is negative while the thermal expansion coefficient is positive, thus the direction of spectrum shift depends on the value of the two coefficients. For our FP cavity structure, there is a large thermal expansion coefficient because of the use of UV glue, resulting to red-shift of the spectrum. In addition, the sensitivity of the FP1 (temperature sensing channel) is larger than that of the FP2 (RI sensing channel), which is resulted from the larger thermo-optical coefficient of liquid sample than that of the air cavity. The sensing curves of the total six sub-FP sensors are shown in Fig. 8(c). The sensitivities of three sub-FP sensors as RI sensing channel are 0.042 nm/°C, 0.041 nm/°C, and 0.046 nm/°C, respectively, and the other three sub-FP sensors as temperature sensing are 0.185 nm/°C, 0.158 nm/°C, and 0.176 nm/°C, respectively.

 figure: Fig. 8.

Fig. 8. The interference spectra shift of (a) the temperature sensing channel FP1 and (b) the RI sensing channel FP2 as a function of temperature, from 20 °C to 80 °C with the interval of 10 °C; (c) the temperature sensing curves of six FP sensors from FP1 to FP6, correspondingly.

Download Full Size | PDF

Finally, the RI and temperature change ΔRIi and ΔTi of the cascaded FP sensor #i (i=1,2,3……) can be demodulated by the sensing matrix, which can be described as:

$$\begin{array}{ccc} {\left[ {\begin{array}{c} {\Delta R{I_i}}\\ {\Delta {T_i}} \end{array}} \right] = \frac{1}{{{S_{RI\_RI,i}}{S_{T\_T,i}}}}\left[ {\begin{array}{cc} {{S_{T\_T,i}}}&{ - {S_{T\_RI,i}}}\\ 0&{{S_{RI\_RI,i}}} \end{array}} \right]\left[ {\begin{array}{c} {\Delta {\lambda_{RI,i}}}\\ {\Delta {\lambda_{T,i}}} \end{array}} \right]}&,&{i = 1,2,3 \ldots \ldots } \end{array}$$
where SRI_RI,i means RI sensitivity of the RI sensing channel, ST_RI,i means temperature sensitivity of the RI sensing channel, ST_T,i means temperature sensitivity of the temperature sensing channel, ΔλRI,i means wavelength shift of the RI sensing channel, ΔλT,i means wavelength shift of the temperature sensing channel, respectively. According to the sensing matrix, the simultaneous measurement of RI and temperature is realized by the cascaded FP sensor. Meanwhile, by means of the proposed multiplexing system, the multi-point sensing is also realized.

In order to evaluate the detection accuracy of the multi-point RI sensing system, we repeat the experiment under the same conditions for 20 times. The interference spectra of 20 experiments are demodulated successively, as shown in Fig. 9(a), in which every interference spectra are all almost overlapped. Figure 9(b) shows the close view of the second interference dip at 1542.25 nm, presenting a tiny wavelength shift about 50 pm between 20 experiment results. The wavelength shift of the interference spectrum is mainly produced by the electric noise of receivers and the wavelength drift of the TLS. The former is difficult to avoid while the latter can be compensated by the acetylene gas absorption spectrum. Figure 9(c) shows the evolution of the acetylene gas absorption spectrum during 20 experiments, presenting a tiny wavelength shift about 30 pm. The wavelength shift of the acetylene gas absorption spectrum indicates the wavelength drift of the TLS, thus the interference spectrum of the FP sensors can be calibrated by subtracting the wavelength shift of the acetylene gas absorption spectrum. Figure 9(d) shows the wavelength shift of the interference spectrum before and after calibration, the standard deviation (SD) of the interference spectra is reduced from 11.6 pm to 6.6 pm. Finally, the detection precision can be obtained by substituting SD and sensitivities into Eq. (5). Table 3 shows the sensitivities and wavelength SD of three cascaded FP sensors, and the calculated results show the high temperature precision about 0.05 °C and high RI precision better than 10−5 RIU.

 figure: Fig. 9.

Fig. 9. (a) The interference spectra of 20 times repetitive experiments before wavelength calibration; (b) close view of the interference dip at 1542.25 nm; (c) acetylene gas absorption spectra during 20 experiments; (d) the wavelength shift of interference spectrum before and after wavelength calibration.

Download Full Size | PDF

Tables Icon

Table 3. Detection precision of three cascaded FP sensors

Although only three cascaded FP sensors are multiplexed for demonstration in this work, we can also increase the sensors number by connecting one coupler after another. The prospective multiplexing number is determined by the minimum detectable light power. In order to estimate the prospective multiplexing number, we research the minimum detectable light power. We place a variable optical attenuator (VOA) in the measurement arm to change the input light power. For the Nth cascaded FP sensor, the reflection light power can be calculated by

$${P_{out}} = {P_{in}}{r_{FP}}r_N^2\prod\limits_{n = 1}^N {{{({1 - {r_{n - 1}}} )}^2}}$$
where Pin represents the input light power of the measurement arm, rFP represents the reflectivity of the FP cavity (rFP ≈ 0.4% for liquid cavity), and rn represents the coupling ratio of the nth optical fiber coupler (rn = 5% in the experiment). When Pin decreases from 1 mW to 100 nW, the reflection peak of the sub-FP sensor also decreases from 12 dBm to −9 dBm, as shown in Fig. 10(a). However, the normalized interference spectrum keeps unchanged at different input light power, as shown in Fig. 10(b). According to Eq. (6), when Pin = 100 nW, Pout = 250rFP pW = 1 pW. In other words, the minimum detectable light power is smaller than 1 pW. Assuming the reflection light power of every sub-FP sensors are all 1 pW, and the light source power is 4 mW, then the maximum multiplexing number can be estimated to Nmax=[(4rFP mW)/(250rFP pW)]0.5 = 4000, which is much large than the paper [34]. In addition, the multiplexing number can be further improved by enlarging the light source power.

 figure: Fig. 10.

Fig. 10. (a) The reflection peaks at different input light power of the measurement arm. (b) The normalized interference spectra correspondingly.

Download Full Size | PDF

Table 4 shows the performance comparison among different multi-point RI sensing systems. Compared to other multi-point RI sensing system, our sensing system has higher sensitivity, which results higher RI resolution. The sensing system has also the ability of temperature compensation, which is very important in RI measurement. Most important, the sensing points’ number of the sensing system is easy to be extended by connecting one coupler after another. The only limitation is that the RI detectable range is narrow (0.008 RIU), but it can also be enlarged by reducing the cavity length of the sub-FP sensors. Noting that the shortest cavity length should larger than twice the spatial resolution of the sensing system for demodulating the interference spectrum. Concretely, the shortest cavity length should larger than 40 µm, thus the maximum RI detectable range is 0.019 RIU (Δnmax = λ/(2LFP)).

Tables Icon

Table 4. Performance comparison among different multi-point RI sensing systemsa

5. Conclusion

In summary, we proposed and demonstrated a novel multi-point RI sensing system with temperature compensation, by the combination of the cascaded FP sensors and the multiplexing method based FMCW technique. The theoretical model of the sensing system is established and the signal demodulation algorithm is built up correspondingly. The blended original signal is divided in spatial domain by FFT, and then the interference spectrum is demodulated independently in wavelength domain by IFFT. Experimentally, the multi-point RI and temperature simultaneous sensing experiments is performed and the well sensing capabilities is realized. The RI sensitivities of the three multiplexed cascaded FP sensors are high up to 1203.9 nm/RIU, 1279.8 nm/RIU and 1230.3 nm/RIU respectively, under RI range from 1.3330 RIU to 1.3410 RIU. Meanwhile, the temperature sensitivities of the three multiplexed cascaded FP sensors are high up to 0.185 nm/°C, 0.158 nm/°C, and 0.176 nm/°C respectively, under temperature range from 20 °C to 80 °C. The RI precision is as high as 10−5 RIU and the temperature precision is as high as 0.05 °C. In addition, the prospective multiplexing number is estimated according to the minimum detectable power, reaching about 4000. The proposed sensing system has potential advantages in several practical applications that require a large number sensing points.

Funding

National Key Scientific Instrument and Equipment Development Projects of China (2017YFF0108700); National Natural Science Foundation of China (61975045); Natural Science Foundation of Heilongjiang Province (LH2020F014).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. A. Urrutia, I. Del Villar, P. Zubiate, and C. R. Zamarreno, “A comprehensive review of optical fiber refractometers: toward a standard comparative criterion,” Laser Photonics Rev. 13(11), 1900094 (2019). [CrossRef]  

2. T. Wei, Y. K. Han, Y. J. Li, H. L. Tsai, and H. Xiao, “Temperature-insensitive miniaturized fiber inline Fabry-Perot interferometer for highly sensitive refractive index measurement,” Opt. Express 16(8), 5764–5769 (2008). [CrossRef]  

3. F. Chiavaioli, F. Baldini, S. Tombelli, C. Trono, and A. Giannetti, “Biosensing with optical fiber gratings,” Nanophotonics 6(4), 663–679 (2017). [CrossRef]  

4. P. Singh, “SPR biosensors: historical perspectives and current challenges,” Sens. Actuators, B 229, 110–130 (2016). [CrossRef]  

5. L. Wan, H. Chandrahalim, J. Zhou, Z. H. Li, C. Chen, S. H. Cho, H. Zhang, T. Mei, H. P. Tian, Y. Oki, N. Nishimura, X. D. Fan, and L. J. Guo, “Demonstration of versatile whispering-gallery micro-lasers for remote refractive index sensing,” Opt. Express 26(5), 5800–5809 (2018). [CrossRef]  

6. Z. Chen, J. Y. Tang, R. Y. Fan, Y. C. Zhong, J. Zhang, and S. B. Li, Side polished fiber Bragg grating sensor for simultaneous measurement of refractive index and temperature, 21st International Conference on Optical Fiber Sensors7753 (2011).

7. W. B. Hu, C. Li, S. Cheng, F. Mumtaz, C. Du, and M. H. Yang, “Etched multicore fiber Bragg gratings for refractive index sensing with temperature in-line compensation,” OSA Continuum 3(4), 1058–1067 (2020). [CrossRef]  

8. H. Y. Meng, W. Shen, G. B. Zhang, C. H. Tan, and X. G. Huang, “Fiber Bragg grating-based fiber sensor for simultaneous measurement of refractive index and temperature,” Sens. Actuators, B 150(1), 226–229 (2010). [CrossRef]  

9. W. Zhang, Y. G. Liu, T. Zhang, D. Q. Yang, Y. X. Wang, and D. K. Yu, “Integrated fiber-optic Fabry-Perot interferometer sensor for simultaneous measurement of liquid refractive index and temperature,” IEEE Sens. J. 19(13), 5007–5013 (2019). [CrossRef]  

10. Y. Ouyang, X. F. Xu, Y. J. Zhao, A. Zhou, and L. B. Yuan, “Temperature compensated refractometer based on parallel fiber Fabry-Perot interferometers,” IEEE Photon. Technol. Lett. 30(13), 1262–1265 (2018). [CrossRef]  

11. X. G. Li, S. C. Warren-Smith, L. S. Xie, H. Ebendorff-Heidepriem, and L. V. Nguyen, “Temperature-compensated refractive index measurement using a dual Fabry-Perot interferometer based on C-fiber cavity,” IEEE Sens. J. 20(12), 6408–6413 (2020). [CrossRef]  

12. H. P. Luo, Q. Z. Sun, Z. L. Xu, W. H. Jia, D. M. Liu, and L. Zhang, “Microfiber-based inline Mach-Zehnder interferometer for dual-parameter measurement,” IEEE Photonics J. 7(2), 1–8 (2015). [CrossRef]  

13. X. Q. Ni, M. Wang, D. M. Guo, H. Hao, and J. L. Zhu, “A hybrid Mach-Zehnder interferometer for refractive index and temperature measurement,” IEEE Photon. Technol. Lett. 28(17), 1850–1853 (2016). [CrossRef]  

14. H. P. Luo, Q. Z. Sun, Z. L. Xu, D. M. Liu, and L. Zhang, “Simultaneous measurement of refractive index and temperature using multimode microfiber-based dual Mach-Zehnder interferometer,” Opt. Lett. 39(13), 4049–4052 (2014). [CrossRef]  

15. Y. Q. Yuan, L. N. Wang, and J. Huang, “Theoretical investigation for two cascaded SPR fiber optic sensors,” Sens. Actuators, B 161(1), 269–273 (2012). [CrossRef]  

16. F. Xiao, D. Michel, G. Y. Li, A. S. Xu, and K. Alameh, “Simultaneous measurement of refractive index and temperature based on surface plasmon resonance sensors,” J. Lightwave Technol. 32(21), 4169–4173 (2014). [CrossRef]  

17. L. Liu, Z. H. Liu, Y. Zhang, and S. T. Liu, “V-shaped micro-structure optical fiber surface plasmon resonance sensor for the simultaneous measurement of the refractive index and temperature,” Opt. Lett. 44(20), 5093–5096 (2019). [CrossRef]  

18. T. Hu, Y. Zhao, and A. N. Song, “Fiber optic SPR sensor for refractive index and temperature measurement based on MMF-FBG-MMF structure,” Sens. Actuators, B 237, 521–525 (2016). [CrossRef]  

19. L. Zhao, Y. D. Zhang, Y. H. Chen, and J. F. Wang, “Simultaneous measurement of temperature and RI based on an optical microfiber coupler assembled by a polarization maintaining fiber,” Appl. Phys. Lett. 114(15), 151903 (2019). [CrossRef]  

20. D. Yi, Y. Z. Chen, Y. F. Geng, F. Teng, Y. Li, F. Liu, X. J. Li, and X. M. Hong, “Interrogation technique analyses of a hybrid fiber optic sensor based on SPR and MMI,” Opt. Express 28(14), 20764–20772 (2020). [CrossRef]  

21. W. H. Zhang, W. L. Gao, Z. R. Tong, Y. M. Zhong, L. F. Xue, and H. W. Zhang, “Mach-Zehnder interferometer cascaded with FBG for simultaneous measurement of RI and temperature,” Opt. Commun. 466, 125624 (2020). [CrossRef]  

22. Y. Zhao, X. G. Li, L. Cai, and Y. N. Zhang, “Measurement of ri and temperature using composite interferometer with hollow-core fiber and photonic crystal fiber,” IEEE Trans. Instrum. Meas. 65(11), 2631–2636 (2016). [CrossRef]  

23. Y. Zhao, X. Liu, R. Q. Lv, and Q. Wang, “Simultaneous measurement of RI and temperature based on the combination of Sagnac loop mirror and balloon-like interferometer,” Sensor. Actuat. B: Chem. 243, 800–805 (2017). [CrossRef]  

24. Z. H. Yu and A. B. Wang, “Fast demodulation algorithm for multiplexed low-finesse Fabry-Perot interferometers,” J. Lightwave Technol. 34(3), 1015–1019 (2016). [CrossRef]  

25. D. Tosi, “Simultaneous detection of multiple fiber-optic Fabry-Perot interferometry sensors with cepstrum-division multiplexing,” J. Lightwave Technol. 34(15), 3622–3627 (2016). [CrossRef]  

26. W. Eickhoff and R. Ulrich, “Optical frequency domain reflectometry in single-mode fiber,” Appl. Phys. Lett. 39(9), 693–695 (1981). [CrossRef]  

27. M. Froggatt and J. Moore, “High-spatial-resolution distributed strain measurement in optical fiber with Rayleigh scatter,” Appl. Opt. 37(10), 1735–1740 (1998). [CrossRef]  

28. G. Yin, L. Lu, L. Zhou, C. Shao, Q. Fu, J. Zhang, and T. Zhu, “Distributed directional torsion sensing based on an optical frequency domain reflectometer and a helical multicore fiber,” Opt. Express 28(11), 16140–16150 (2020). [CrossRef]  

29. Z. P. Zhang, X. Y. Fan, and Z. Y. He, “Long-range distributed static strain sensing with <100 nano-strain resolution realized using OFDR,” J. Lightwave Technol. 37(18), 4590–4596 (2019). [CrossRef]  

30. M. Zhu and H. Murayama, “Fast demodulation of OFDR based long length FBG sensing system for noisy signals,” Opt. Express 26(16), 19804–19814 (2018). [CrossRef]  

31. X. T. Lou, Y. B. Feng, C. Chen, and Y. K. Dong, “Multi-point spectroscopic gas sensing based on coherent FMCW interferometry,” Opt. Express 28(6), 9014–9026 (2020). [CrossRef]  

32. P. K. C. Chan, W. Jin, J. M. Gong, and M. S. Demokan, “Multiplexing of fiber Bragg grating sensors using an FMCW technique,” IEEE Photon. Technol. Lett. 11(11), 1470–1472 (1999). [CrossRef]  

33. M. Sypabekova, S. Korganbayev, W. Blanc, T. Ayupova, A. Bekmurzayeva, M. Shaimerdenova, K. Dukenbayev, C. Molardi, and D. Tosi, “Fiber optic refractive index sensors through spectral detection of Rayleigh backscattering in a chemically etched MgO-based nanoparticle-doped fiber,” Opt. Lett. 43(24), 5945–5948 (2018). [CrossRef]  

34. Y. Du, S. Jothibasu, Y. Y. Zhuang, C. Zhu, and J. Huang, “Rayleigh backscattering based macrobending single mode fiber for distributed refractive index sensing,” Sens. Actuators, B 248, 346–350 (2017). [CrossRef]  

35. Z. Y. Ding, K. L. Sun, K. Liu, J. F. Jiang, D. Yang, Z. Yu, J. Li, and T. G. Liu, “Distributed refractive index sensing based on tapered fibers in optical frequency domain reflectometry,” Opt. Express 26(10), 13042–13054 (2018). [CrossRef]  

36. G. L. Zhang, M. H. Yang, and M. Wang, “Large temperature sensitivity of fiber-optic extrinsic Fabry-Perot interferometer based on polymer-filled glass capillary,” Opt. Fiber Technol. 19(6), 618–622 (2013). [CrossRef]  

37. L. J. Qiao, D. S. Sun, X. F. Zhang, and Y. Zhao, “Linearity requirement for a linear frequency modulation lidar,” Opt. Rev. 6(2), 160–162 (1999). [CrossRef]  

38. B. J. Soller, D. K. Gifford, M. S. Wolfe, and M. E. Froggatt, “High resolution optical frequency domain reflectometry for characterization of components and assemblies,” Opt. Express 13(2), 666–674 (2005). [CrossRef]  

39. C. Zhang, S. Xu, J. F. Zhao, H. Q. Li, H. Bai, and C. Y. Miao, “Multipoint refractive index and temperature fiber optic sensor based on cascaded no core fiber-fiber Bragg grating structures,” Opt. Eng. 56(2), 027102 (2017). [CrossRef]  

40. T. Mukai and H. Fukano, “Multipoint refractive index measurement using multimode interference-based fiber-optic sensors driven by an integrable tunable laser assembly,” Jpn. J. Appl. Phys. 59(SO), SOOE02 (2020). [CrossRef]  

41. J. H. Lopez, O. Esteban, M. G. Shlyagin, and R. Martinez-Manuel, “Multipoint refractometer based on combined correlation and frequency multiplexing,” IEEE Photon. Technol. Lett. 29(17), 1479–1482 (2017). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. (a) Basic configuration of the proposed FMCW-based multiplexing method for multi-point simultaneous measuring RI and temperature. A series of cascaded FP sensors are connected in parallel by optical couplers. (b) The structure of the cascaded FP sensor. (c)-(e) The signal demodulation processing of the multiplexing system.
Fig. 2.
Fig. 2. (a) Simulation result of the original beat signal. (b) Power spectrum of the beat signal after applying FFT. The reflected peaks appear at the preset location. (c) Close view of the reflected peaks generated by temperature sensing channels FP1, FP3 and FP5, respectively. (d) The interference spectra of the temperature sensing channels after applying IFFT. (e) Close view of the reflected peaks generated by RI sensing channels FP2, FP4 and FP6, respectively. (f) The interference spectra of the RI sensing channels after applying IFFT.
Fig. 3.
Fig. 3. The effect of the multiplexing number on the wavelength error of an individual interference spectrum. Inset: the interference spectrum of an individual sub-FP sensor when multiplexing number are 1, 10, 100, 1000, and 5000 respectively.
Fig. 4.
Fig. 4. (a) The sketch diagram of the proposed FMCW-based FP sensors multiplexing system; (b) photograph of the cascaded FP sensors for simultaneous RI and temperature sensing. TLS: tunable laser source; FRM: Faraday rotating mirror; PC: polarization controller; CG: clock generator; BPD: balanced photodetector; DAQ: data acquisition card; SMF: single mode fiber.
Fig. 5.
Fig. 5. (a) The original experimental beat signal; (b) power spectrum of the beat signal; (c) close view of the reflection peaks generated by FP1; (d) the original interference spectrum of the FP1 obtained after applying the IFFT process; (e) the normalized interference spectrum of the FP1 obtained after removing the baseline.
Fig. 6.
Fig. 6. (a) Close view of the reflection peaks generated by six FP sensors from FP1 to FP6, respectively; (b) the normalized interference spectra (solid lines) of six FP sensors from FP1 to FP6 respectively, which are almost consistent with the measurement results by the traditional optical spectrum analyzer (dashed lines).
Fig. 7.
Fig. 7. The interference spectra shift of (a) the temperature sensing channel FP1 and (b) the RI sensing channel FP2 as a function of RI of the sample under test, from 1.3330 RIU to 1.3410 RIU with the interval of 0.0010 RIU; (c) the RI sensing curves of six FP sensors from FP1 to FP6, correspondingly.
Fig. 8.
Fig. 8. The interference spectra shift of (a) the temperature sensing channel FP1 and (b) the RI sensing channel FP2 as a function of temperature, from 20 °C to 80 °C with the interval of 10 °C; (c) the temperature sensing curves of six FP sensors from FP1 to FP6, correspondingly.
Fig. 9.
Fig. 9. (a) The interference spectra of 20 times repetitive experiments before wavelength calibration; (b) close view of the interference dip at 1542.25 nm; (c) acetylene gas absorption spectra during 20 experiments; (d) the wavelength shift of interference spectrum before and after wavelength calibration.
Fig. 10.
Fig. 10. (a) The reflection peaks at different input light power of the measurement arm. (b) The normalized interference spectra correspondingly.

Tables (4)

Tables Icon

Table 1. Simulation parameters of the multiplexing system

Tables Icon

Table 2. Experiment parameters of the multiplexed sub-FP sensors

Tables Icon

Table 3. Detection precision of three cascaded FP sensors

Tables Icon

Table 4. Performance comparison among different multi-point RI sensing systems a

Equations (6)

Equations on this page are rendered with MathJax. Learn more.

E t o t a l ( t ) = E R ( t ) + m = 1 M r m E P m ( t τ m )
E R = E R 0 e j φ ( t )
φ ( t ) = 2 π f 0 t + π γ t 2
I t o t a l = E t o t a l E t o t a l
[ Δ R I i Δ T i ] = 1 S R I _ R I , i S T _ T , i [ S T _ T , i S T _ R I , i 0 S R I _ R I , i ] [ Δ λ R I , i Δ λ T , i ] , i = 1 , 2 , 3
P o u t = P i n r F P r N 2 n = 1 N ( 1 r n 1 ) 2
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.