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

Analysis on real-time phase delay in an interferometric FBG sensor array using polarization switching and the PGC hybrid processing method

Open Access Open Access

Abstract

The polarization switching and phase generated carrier (PS-PGC) hybrid method is typically adopted to control signal fading induced by fiber birefringence and to precisely demodulate signals in interferometric fiber Bragg grating (FBG) sensor arrays. Unlike simple PGC demodulation, the real-time phase delay between the detected interference and the PGC carrier in the hybrid method has more adverse effects as both demodulation accuracy and background noise can be deteriorated, which may invalidate the polarization switching (PS) method. Aiming at this issue, the real-time phase delay and its compensation method were analyzed in detail in this paper. The features of the real-time phase delay in the PS-PGC hybrid method were summarized and the differences among polarization channels were investigated. Theoretical analysis indicated that the real-time phase delay mainly affected the interference complex synthesis (ICS) procedure, ultimately bringing about errors in the PS-PGC algorithm. The method for demodulating the real-time phase delay from sampled interference was presented, which provided the key compensation parameter. Experimental results showed that the compensation method could greatly improve the stability of the demodulated signal and suppress the sensor background phase noise. The amplitude of the demodulated signal was stabilized with a fluctuation less than ± 0.75dB and a noise suppression of 5dB. The acceptable compensation error was also analyzed.

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

1. Introduction

The inline interferometric fiber Bragg grating (FBG) sensor array has been widely used in ultra-weak signal sensing applications, such as acoustic detection [1], seismic sensing [24]. The inherently integrated sensor array structure has only FBGs written inline, providing benefits of simple and compact structure, potential low cost and high reliability [56], making it especially attractive for applications in harsh environments. The interferometric interrogation method is of high sensitivity in its physical nature, providing an excellent foundation for ultra-weak signal sensing [7].

Compared to conventional fiber interferometers, it is more difficult for the inline FBG interferometer to obtain stable demodulation and low background noise, as no optical fiber devices can be added in the array to eliminate disturbances from the interrogation system, the leading fiber and even the sensor itself. The most prominent obstacle for the inline FBG interferometer is the polarization-induced signal fading [8], which is mainly caused by random fiber birefringence fluctuations and is usually solved using a pair of Faraday rotator mirrors in the Michelson interferometer [910]. The PS method is one solution that mainly modulates interrogation pulses [1112]. Besides the unique problem for the FBG structure, common problems in fiber interferometers should also be addressed, such as random sensor sensitivity fluctuations [13] mainly caused by environment-induced sensor phase drifting. Fortunately, those common problems have been well studied with lots of solutions proposed, such as the PGC method [14]. Considering both unique and common problems, it can be known that the interrogation design and the demodulation algorithm are the essence of an excellent inline interferometric FBG sensor system. Therefore, ingenious and precise demodulation parameters are of great significant, even determining whether the sensor can be used in practical applications. The polarization switching and phase generated carrier (PS-PGC) hybrid method was proposed under this background and has been proved useful by a series of experiments [712].

This paper focused on performance promotion for the interferometric FBG sensor array using the PS-PGC hybrid method, considering the real-time phase delay between the detected interference and the PGC carrier as a key demodulation parameter. In a simple PGC system, the real-time phase delay has been noticed and its compensation method has attracted much attention recently [1522]. Early methods tried to control the initial phase of PGC carrier from the perspective of hardware [15]. However, the real-time phase delay is related to every part of the demodulation system, difficult to be controlled accurately. Later, some researches tried to pre-estimate the real-time phase delay and compensate it in the demodulation algorithm [1622], whose focus was still the pre-estimated result of the real-time phase delay. Compared to simple PGC demodulation, the real-time phase delay in the PS-PGC hybrid method has more sophisticated effects. We have fully demonstrated the PS-PGC method in Ref. [7] and Ref. [12], showing that the PGC carrier existed in all of four polarization channels. Therefore, both the suppression of polarization-induced signal fading and the picking-up of the sensor phase were highly related to the real-time phase delay. In this paper, the real-time phase delay and its compensation method in the PS-FGC hybrid method were analyzed in detail. The principles of how the real-time phase delay affects the PS-PGC hybrid processing procedure and corresponding results were investigated. Real-time phase delays of all polarization channels were analyzed. A compensation method was proposed, which was experimentally demonstrated to be effective in promoting sensing performance.

This paper is arranged as follows. Section 2 is about theoretical demonstrations of the real-time phase delay in the PS-PGC system and its compensation method. Section 3 is on experimental results with the measured real-time phase delay in the PS-PGC system and the evaluated compensation effects. In Section 4 a brief summarization is presented.

2. Theory

2.1 Basic principle of interferometric FBG sensor array

A typical structure of an interferometric FBG sensor array is shown in Fig. 1. Three FBGs are written inline in a single fiber with absolutely uniform reflection centers. If interrogated using a pair of laser pulses, interference occurs when the second pulse reflected by the ith FBG overlaps with the first pulse reflected by the (i + 1)th FBG(i = 0, 1), as shown in Fig. 1. Two-time division multiplexed sensor elements named S1 and S2 can thus be formed, each containing two adjacent FBGs and a length of fiber between the two FBGs.

 figure: Fig. 1.

Fig. 1. Typical structure of interferometric FBG sensor array.

Download Full Size | PDF

The two reflected pulses for S1 and S2 can be written as

$$\begin{aligned} {{\textbf{E}}_{r1}} &= r{\rho _\textrm{0}}{E_{in1}}\\ {{\textbf{E}}_{s1}} &= {t_0}{\overleftarrow {\textbf{B}} _1}r{\rho _\textrm{1}}{\overrightarrow {\textbf{B}} _1}{t_0}{{\textbf{E}}_{in0}}\\ {{\textbf{E}}_{r2}} &= {t_0}{\overleftarrow {\textbf{B}} _1}r{\rho _\textrm{1}}{\overrightarrow {\textbf{B}} _1}{t_0}{{\textbf{E}}_{in1}}\\ {{\textbf{E}}_{s2}}& = {t_0}{\overleftarrow {\textbf{B}} _1}{t_1}{\overleftarrow {\textbf{B}} _2}r{\rho _\textrm{2}}\overrightarrow {{{\textbf{B}}_2}} {t_1}\overrightarrow {{{\textbf{B}}_1}} {t_0}{{\textbf{E}}_{in0}} \end{aligned}$$
Here ρ0, ρ1 and ρ2 are the amplitude reflectivity of the three FBGs, respectively. ${E_{in0}} = {[{{E_{x0}}\textrm{ }{E_{y0}}} ]^T}$ is the Jones matrix of the first injected laser pulse. A modulation term of $C\cos (2\pi {f_0}t)]$ is introduced to the second laser pulse and ${E_{in1}} = {[{{E_{x1}}\textrm{ }{E_{y1}}} ]^T}\exp [\textrm{ - }iC\cos (2\pi {f_0}t)]$ is the Jones matrix of the second injected pulse, where ${f_0}$ and C represent the PGC modulation frequency and amplitude, respectively. ${t_i}(i = 0,1)$ is the amplitude transmission of the first and second FBGs. r is the FBG reflection matrix and $r = \,\left( {\begin{array}{cc} { - 1} &0\\ 0 &1 \end{array}} \right)$ . ${B_n}(n = 1,2)$ is the transmission matrix of S1 and S2, while the superscripts $\leftarrow $ and $\to $ indicate backward and forward transmission. Given an negligible fiber transmission loss, ${B_n}$ can be written as ${B_n} = \exp ( - i{\varphi _n}){U_n}$, where ${\varphi _n}$ is the phase delay of single transmission in the sensing fiber of Sn, and ${U_n}$ is a 2×2 unitary matrix related to the fiber birefringence transmission matrix [11]. For a single mode fiber, the following equation describes properties of ${B_n}$ [23].
$$B_n^\dagger {B_n} = B_n^TB_n^\ast{=} \left( \begin{array}{cc} 1 &0\\ 0 &1 \end{array} \right),\overleftarrow B _n^T = r{\overrightarrow B _n}r$$
The superscript * and T refer to the conjugate convertor and the matrix transposition symbol, respectively, while $\dagger$ refers to the matrix transform and conjugate convertor. Equation (1) can then be expressed as
$$\begin{aligned} {{\textbf{E}}_{r1}} &= r{\rho _\textrm{0}}{E_{in1}}\\ {{\textbf{E}}_{s1}} &= t_0^2{\rho _\textrm{1}}r{\overrightarrow {\textbf{B}} _1}^T{\overrightarrow {\textbf{B}} _1}{{\textbf{E}}_{in0}}\\ {{\textbf{E}}_{r2}} &= t_0^2{\rho _\textrm{1}}r{\overrightarrow {\textbf{B}} _1}^T{\overrightarrow {\textbf{B}} _1}{{\textbf{E}}_{in1}}\\ {{\textbf{E}}_{s2}} &= t_0^2t_1^2{\rho _\textrm{2}}r{\overrightarrow {\textbf{B}} _1}^T{\overrightarrow {\textbf{B}} _2}^T\overrightarrow {{{\textbf{B}}_2}} \overrightarrow {{{\textbf{B}}_1}} {{\textbf{E}}_{in0}} \end{aligned}$$
The two interferences, corresponding to S1 and S2, can be written as
$$\begin{aligned} {I_1} &= {({{E_{r1}} + {E_{s1}}} )^\dagger }({{E_{r1}} + {E_{s1}}} )= D{C_1} + 2t_{_0}^2{\rho _0}{\rho _1}{\textrm{Re}} ({E_{_{in1}}^\dagger {{\overrightarrow {{B_1}} }^T}\overrightarrow {{B_1}} {E_{in0}}} )\\ &= D{C_1} + 2t_{_0}^2{\rho _0}{\rho _1}\exp ( - i{\varphi _{s1}}) {\textrm{Re}} ({E_{_{in1}}^\dagger {{\overrightarrow {{U_1}} }^T}\overrightarrow {{U_1}} {E_{in0}}} )\\ {I_2} &= {({{E_{r2}} + {E_{s2}}} )^\dagger }({{E_{r2}} + {E_{s2}}} )= D{C_2} + 2t_0^4t_1^2{\rho _1}{\rho _2}{\textrm{Re}} (E_{_{in1}}^\dagger {\overrightarrow {{B_1}} ^\dagger }{\overrightarrow {{B_2}} ^T}\overrightarrow {{B_2}} \overrightarrow {{B_1}} {E_{in0}})\\ = &D{C_2} + 2t_0^4t_1^2{\rho _1}{\rho _2}\exp ( - i{\varphi _{s2}}) \Re (E_{_{in1}}^\dagger {\overrightarrow {{U_1}} ^\dagger }{\overrightarrow {{U_2}} ^T}\overrightarrow {{U_2}} \overrightarrow {{U_1}} {E_{in0}}) \end{aligned}$$
Here ${\varphi _{sn}} = 2{\varphi _n}$ is the signal of concerned phase in S1 and S2. Equation (4) indicates that the interference is a mixed result of ${\varphi _{sn}}$ and ${U_n}$. Therefore, the extraction of ${\varphi _{sn}}$ directly from interference is inevitably influenced by fiber birefringence, which can be easily disturbed by the environment. As a result, the demodulated signal amplitude and the system background noise fluctuate randomly with environmental disturbances, a phenomenon called polarization-induced signal fading [10].

2.2 Brief review of PS-PGC hybrid processing method

We have demonstrated the PS- PGC hybrid processing method in previous reports [712]. In order to analyze the influence of real-time phase delay, a brief review is given here. The sensor impulse is defined as

$${\Re _1} = \exp ( - i{\varphi _{s1}}){\overrightarrow {{U_1}} ^T}\overrightarrow {{U_1}} = \exp ( - i{\varphi _{s1}})U_1^{\prime} \quad {\Re _2} = \exp ( - i{\varphi _{s2}}){\overrightarrow {{U_1}} ^\dagger }{\overrightarrow {{U_2}} ^T}\overrightarrow {{U_2}} \overrightarrow {{U_1}} = \exp ( - i{\varphi _{s2}})U_{_2}^{\prime}$$
Here $U_n^{\prime}$ is also a 2×2 unitary matrix since it is a multiplexing result of several unitary matrixes. The sensor phase thus can be derived from $\sqrt {\det {\Re _n}} = \exp ( - i{\varphi _{sn}})$ which can eliminate all factors related to fiber birefringence. It should be noticed here that ${\Re _n}$ is a complex expression which cannot be detected directly as opto-electronic detectors only respond to the real part of interference. Thus, an interference complex synthesis (ICS) procedure is critical.

In a practical experimental system, the key device to obtain ${\Re _n}$ is a PS module. The laser input to the PS module is typically linearly polarized, and this polarization state is called the X polarization. If a half-wave voltage is acted upon the PS module, the output laser can be modulated to its orthogonal polarization state, i.e., the Y polarization. Thus, we can control the polarization of two interfering laser pulses. If the polarization states of two laser pulses are modulated into XX, XY, YY and YX successively, as shown in Fig. 2, the interferences of four polarization channels for each sensing unit can be obtained, named as IXX, IXY, IYY and IYX, respectively.

 figure: Fig. 2.

Fig. 2. Interrogation and interference pulses.

Download Full Size | PDF

The Jones matrix of the laser pulse pair in the XX polarization channel can be expressed as in Eq. (6) given that the laser intensity is normalized to be 1.

$${E_{in0}} = {[{1\textrm{ }0} ]^T} \quad {E_{in1}} = {[{1\textrm{ }0} ]^T}\exp [\textrm{ - }iC\cos (2\pi {f_0}t)]$$
The interference IXX thus can be written as
$${I_{XXn}} = {I_{n\_DC}} + t_n^{\prime}\Re \{ \exp ( - i{\varphi _{\textrm{s}n}}) \exp [{ - iC\cos (2\pi {f_0}t)} ]U_{11n}^{\prime}\} ( n = 1,2)$$
Here $t_1^{\prime} = 2t_{_0}^2{\rho _0}{\rho _1}, t_2^{\prime} = 2t_0^4t_1^2{\rho _1}{\rho _2}$, $U_{11n}^{\prime}$ is the first matrix element of $U_n^{\prime}$. Similarly, the laser Jones matrixes of XY polarization, YX polarization and YY polarization channels can be expressed as
$$\begin{aligned}{E_{in0}} &= {{[{1\textrm{ }0} ]}^T} \quad {E_{in1}} = {{[{0\textrm{ }1} ]}^T}\exp [\textrm{ - }iC\cos (2\pi {f_0}t)]\\ {E_{in0}} &= {{[{0\textrm{ }1} ]}^T} \quad {E_{in1}} = {{[{1\textrm{ }0} ]}^T}\exp [\textrm{ - }iC\cos (2\pi {f_0}t)]\\ {E_{in0}} &= {{[{\textrm{0 1}} ]}^T} \quad {E_{in1}} = {{[{0\textrm{ }1} ]}^T}\exp [\textrm{ - }iC\cos (2\pi {f_0}t)] \end{aligned}$$
The interference terms IXY, IYX and IXY are as follows
$$\begin{aligned} {I_{XYn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ \exp ( - i{\varphi _{sn}})\exp [{ - iC\cos (2\pi {f_0}t)} ]U_{12n}^{\prime}\} \\ {I_{YXn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ \exp ( - i{\varphi _{sn}})\exp [{ - iC\cos (2\pi {f_0}t)} ]U_{21n}^{\prime}\}\\ {I_{YYn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ \exp ( - i{\varphi _{sn}})\exp [{ - iC\cos (2\pi {f_0}t)} ]U_{22n}^{\prime}\} \end{aligned}$$
Here $U_{12n}^{\prime},\, U_{21n}^{\prime}$ and $U_{22n}^{\prime}$ are the other three matrix elements of $U_n^{\prime}$. Then, the interferences of the four polarization channels can be sampled and multiplied with $\cos (2\pi {f_0}t)$ and $\cos (4\pi {f_0}t)$, respectively. Also, a low-pass filter can be used to eliminate carrier harmonics at ${f_0}$ and above, same as in the PGC processing method. The following two items can be obtained.
$$\begin{array}{l} {{X_{1MNn}} ={-} t_n^{\prime}{J_1}(C)U{r_{MNn}}\sin ({\varphi _{sn}} + U{a_{MNn}})}\\ {{X_{2MNn}} ={-} t_n^{\prime}{J_2}(C)U{r_{MNn}}\cos ({\varphi _{sn}} + U{a_{MNn}})} \end{array}\,({MN = 11,12,21,22\textrm{ }})$$
Here $U_{MNn}^{\prime} = U{r_{MNn}}\exp ( - iU{a_{MNn}})$.

Then, the ICS can be completed as

$${\Re _{MNn}} = - \frac{{{X_{2MNn}}}}{{{J_2}(C)}} - i\frac{{{X_{1MNn}}}}{{{J_1}(C)}} = t_n^{\prime}U{r_{MNn}}\exp ( - i{\varphi _{sn}} - iU{a_{MNn}}) = t_n^{\prime}\exp ( - i{\varphi _{sn}}){U_{MNn}^{\prime}}$$
The phase can be derived using the following formula.
$$\sqrt {{\Re _{XXn}} \times {\Re _{YYn}} - {\Re _{XYn}} \times {\Re _{YXn}}} = \sqrt {\det {\Re _n}} = \exp ( - i2{\varphi _{\textrm{s}n}})$$
Now the real and imaginary parts of $\sqrt {\det {\Re _n}}$ can be obtained, which are $\cos (2{\varphi _{\textrm{s}n}})$ and $\sin (2{\varphi _{\textrm{s}n}})$, respectively. Then the different-cross-multiplexing (DCM) algorithm can be used to derive $2{\varphi _n}$, a process that has been demonstrated [1314] and also shown in Fig. 3. Please note that the final signal amplitude is half of the DCM demodulated result.

2.3 Real-time phase delay in PS-PGC hybrid processing method

As demonstrated in Eq. (7) and Eq. (9), a phase modulation of $C\cos (2\pi {f_0}t)$ is expected for the four polarization channels. Although the signal generator can provide an exact signal of $C\cos (2\pi {f_0}t)$, there is a time delay when the carrier signal is acted on ${E_{in1}}$ and detected by the photo-detector. This delay can be caused by the optical path, the circuit transmission, the optical detection process, etc. [1516]. Considering a time delay of $\Delta t$, the real carrier signal is $C\cos [2\pi {f_0}(t\textrm{ + }\Delta t)]$ and thus an original phase delay of $\theta = 2\pi {f_0}\Delta t$ exist in Eq. (7) and Eq. (9), which can be written as

$$\begin{aligned} {I_{XXn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ exp ( - i{\varphi _{\textrm{s}n}}) exp [{ - iC\cos (2\pi {f_0}t\textrm{ + }\theta )} ]U_{11}^{\prime}\} \\ {I_{XYn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ exp ( - i{\varphi _{\textrm{s}n}}) \exp [{ - iC\cos (2\pi {f_0}t\textrm{ + }\theta )} ]U_{12}^{\prime}\} \\ {I_{YXn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ \exp ( - i{\varphi _{\textrm{s}n}}) \exp [{ - iC\cos (2\pi {f_0}t\textrm{ + }\theta )} ]U_{21}^{\prime}\} \\ {I_{YYn}} &= {I_{n\_DC}} + t_n^{\prime}\Re \{ \exp ( - i{\varphi _{\textrm{s}n}}) \exp [{ - iC\cos (2\pi {f_0}t\textrm{ + }\theta )} ]U_{22}^{\prime}\} \end{aligned}$$
If the interferences are multiplied with $\cos (2\pi {f_0}t)$ and $\cos (4\pi {f_0}t)$, and then filtered with a low-pass filter, two additional terms of $\cos (\theta )$ and $\cos (2\theta )$ will be found in Eq. (10), which can be written as
$$\begin{aligned} X_{_{1MNn}}^{\prime} &={-} t_n^{\prime}{J_1}(C)\cos (\theta )U{r_{MNn}}\sin ({\varphi _{sn}} + U{a_{MNn}})\\ X_{_{2MNn}}^{\prime} &={-} t_n^{\prime}{J_2}(C)\cos (2\theta )U{r_{MNn}}\cos ({\varphi _{sn}} + U{a_{MNn}}) \end{aligned}$$
If ${X_{1MNn}}$ and ${X_{2MNn}}$ are directly substituted into Eq. (11), ${\Re _{MNn}}$ cannot be synthesized as either the real or imaginary part of ${\Re _{MNn}}$ is correct. Consequently, when extracting ${\varphi _n}$ from $\sqrt {\det {\Re _n}}$, fiber-birefringence factors cannot be eliminated. If the DCM method is applied, an additional term of $U_{MNn}^{\prime}$ is still found in the demodulated result. $U_{MNn}^{\prime}$ is related to fiber birefringence, which can be easily influenced by environmental disturbances. Specifically, the real part of $U_{MNn}^{\prime}$ is a coefficient of $\exp ( - i{\varphi _{sn}})$, causing instability to the demodulated signal amplitude. The imaginary part of $U_{MNn}^{\prime}$, which is added with ${\varphi _{sn}}$, mainly transforms environmental disturbances into system background noise. Due to the real-time phase delay, both the demodulated signal amplitude stability and background noise are deteriorated, compromising the validity of the PS method.

 figure: Fig. 3.

Fig. 3. DCM procedure.

Download Full Size | PDF

The analysis above indicates that the real-time phase delay mainly affects the ICS procedure and ultimately adds error to the demodulation. In order to compensate for the real-time phase delay, the complex synthesis can be revised as

$${\Re _{MNn}} ={-} \frac{{{X_{2MNn}}}}{{{J_2}(C)\cos (2\theta )}} - i\frac{{{X_{1MNn}}}}{{{J_1}(C)\cos (\theta )}} = t_n^{\prime}U{r_{MNn}}\exp ( - i{\varphi _{sn}} - iU{a_{MNn}}) = t_n^{\prime}\exp ( - i{\varphi _{sn}})U_{MNn}^{\prime}$$
Equation (15) indicates that by introducing $\cos (2\theta )$ and $\cos (\theta )$ to the real and imaginary part, respectively, the exact complex expression of ${\Re _{MNn}}$ can be obtained. However, this compensation method requires that the real-time phase delay should be clearly known in advance. In order to obtain $\theta$, as indicated in Eq. (13), the sampled interference is multiplied with $\cos (2\pi {f_0}t)$ and $\sin (2\pi {f_0}t)$, respectively, and the following results can be obtained after low-pass filtering.
$$\begin{aligned} X_{_{1MNn}}^{\prime\prime} &={-} t_n^{\prime}{J_1}(C)\cos (\theta )U{r_{MNn}}\sin ({\varphi _{sn}} + U{a_{MNn}})\\ X_{_{2MNn}}^{\prime\prime} &={-} t_n^{\prime}{J_1}(C)\sin (\theta )U{r_{MNn}}\sin ({\varphi _{sn}} + U{a_{MNn}}) \end{aligned}$$
Then the real-time phase delay $\theta$ can be extracted by calculating the arctangent:
$$\theta \textrm{ = arc tan(}\frac{{X_{_{2MNn}}^{\prime\prime}}}{{X_{_{1MNn}}^{\prime\prime}}}\textrm{)}$$
As the real-time phase delay $\theta$ is mainly induced by the optical path, the circuit transmission, the photo-detection process, $\theta$ generally has a fixed value, which means that it can be compensated using a certain method.

3. Experimental results and discussion

3.1 System setup

The schematic diagram of the experimental system is shown in Fig. 4. Three FBGs were written in a single-mode fiber. The center wavelength was all 1539.8nm ± 0.01nm with a 3dB bandwidth of approximately 0.36nm, and their relativities were carefully controlled to be 0.6%. The spacing between any two adjacent FGBs was 39.75m.

 figure: Fig. 4.

Fig. 4. Schematic diagram of experimental system.

Download Full Size | PDF

A narrow line-width fiber ring laser was used whose relative intensity noise was −120dB/√Hz @ 1 kHz. The laser wavelength was 1539.8nm, in accordance with the FBG reflection center. The laser output was modulated into pulses by an acoustic-optical modulator (AOM) from GOOCH&HOUSEG with an interrogation frequency of 320kHz. A Mach-Zehnder (M-Z) fiber interferometer was used as the compensation interferometer (CIF), providing a pair of interrogating laser pulses with a time delay determined by optical path imbalance (OPI). The time interval of the interrogating laser pulse pair should be exactly in accordance with the laser round trip time between two adjacent FBGs, as demonstrated in Fig. 1, thus the latter was measured firstly to provide CIF designing instruction. The laser pulse outputted by the AOM was directly injected to the FBGs array and the reflected light pulse train was recorded by the detector and the sampling module shown in Fig. 4. Then an oscilloscope from Tektronix was used to display the recorded results, as shown Fig. 5, and the laser round trip time between two adjacent FBGs was measured to be 392ns. Afterwards, CIF fiber lengths of the two arms were carefully adjusted to generate a laser pulse pair with a time interval of exactly 392ns and the final measured fiber length imbalance was 79.5m. Therefore, a fiber refractive index of 1.479 could be calculated given a light speed in vacuum of 3×108m/s. Part of the short arm of the CIF was wounded on a PZT to introduce the modulation signal for PGC demodulation. The PGC modulation amplitude and frequency were set to be 2.37 and 10kHz, respectively. The pulse pairs were then injected to a PS module from Photline, where the polarization of each pulse can be modulated individually with a PS modulation frequency of 80kHz. All parameters in the experimental system are shown in Table 1 and Table 2.

 figure: Fig. 5.

Fig. 5. Measured result of reflected pulse train.

Download Full Size | PDF

Tables Icon

Table 1. Parameters of FBG sensors array

Tables Icon

Table 2. Parameters of interrogation system

We applied a sinusoidal signal of 1kHz to S1, and the sampled interferences for S1 and S2 are shown in Fig. 6. For both S1 and S2, four polarization channels were obtained. Figure 6(a) shows the interference of S1 in the frequency domain while Fig. 6(b) shows the interference of S2 in the frequency domain. It can be seen that for S1 the 1kHz signal was carried on the PGC modulation frequency at 10kHz and its harmonics above while the PGC modulation frequency was also10kHz and no obvious signal at 1kHz was carried on for S2, indicating a small crosstalk.

 figure: Fig. 6.

Fig. 6. Sampled interferences of (a)S1 and (b) S2 in the frequency domain.

Download Full Size | PDF

3.1 Measurement of real-time phase delay

The real-time phase delays of S1 and S2 were obtained as demonstrated in Section 2.3. We sampled the interferences for 87 times, and the demodulated real-time phase delays of the four polarization channels are shown in Fig. 7 for S1 and S2, respectively. If a certain polarization channel is considered, the real-time phase delays were the same for S1 and S2, as the interrogation pulse pair was the same one though reflected by different FBGs. As shown, the phase delay for the four polarization channels showed a fixed difference of about 0.2rad, mainly determined by the time delay of interrogation pluses. As shown in Fig. 2, the time delay of interrogation pulses between two adjacent polarization channels is the reciprocal of the AOM frequency, which was 3.125µs in our experiment. Therefore, the theoretical phase delay was 0.1963rad given a PGC frequency of 10kHz, in agreement with the experimental results. Furthermore, the PS modulation frequency is generally 8 times of the PGC frequency for full period sampling, and the AOM frequency is 4 times of the PS modulation frequency, as indicated in Fig. 2. As a result, the fixed phase delay difference for the four polarization channels is $2\pi {f_{PGC}}/(32{f_{PGC}})$, where ${f_{PGC}}$ refers to the PGC frequency, which is a constant number of 0.1963rad. This discussion indicates that when the PGC frequency is changed, the measure real-time phase delay remains unchanged.

 figure: Fig. 7.

Fig. 7. Real-time phase delay of (a)S1 and (b) S2.

Download Full Size | PDF

As mentioned above, the real-time phase delay is mainly caused by the optical path, the circuit transmission and the photo-detection process, which is a fixed parameter once the system starts working. Figure 7(a) indicates that for the four polarization channels, the real-time phase delays were −1.19rad, −1rad, −0.6rad and −0.8rad, respectively. Slight fluctuations were observed which may be caused by the phase delay demodulation method. As shown in Eq. (17), the real-time phase delay is obtained by the division of $X_{_{2MNn}}^{\prime\prime}$ and $X_{_{1MNn}}^{\prime\prime}$, which eliminates a common factor of $\sin ({\varphi _{sn}} + U{a_{MNn}})$. When ${\varphi _{sn}} + U{a_{MNn}}$ is close to zero, $\sin ({\varphi _{sn}} + U{a_{MNn}})$ is also a small number close to zero. Therefore, the result of Eq. (17) is not stable since that a small number close to zero is divided. The PGC-DCM method was used to demodulate the corresponding phase in each polarization channel and the result of the real-time phase delay of S1 versus the demodulated signal phase is shown in Fig. 8.

 figure: Fig. 8.

Fig. 8. Real-time phase delay of S1 versus demodulated phase.

Download Full Size | PDF

Figure 8 clearly shows that when the phase $({\varphi _{sn}} + U{a_{MNn}})$ is smaller than 0.2rad, the real-time phase delay is not stable. Therefore, when the demodulated phase is near, $\pi \,/\,2$, the corresponding real-time phase delay can be used as the compensation parameter. Figure 8 also indicates that the demodulated sensor phases are not equal for the four polarization channels as $U{a_{MNn}}$ is related to fiber birefringence. The fluctuations of the real-time delays for the four polarization channels were not identical, showing relatively smaller fluctuations for the YX and YY channels in our experiments.

3.2 Compensation of real-time phase delay

We used the demodulated real-time phase delay of −1.19rad, −1rad, −0.6rad and −0.8rad as the compensation parameters for the XX, XY, YX, YY polarization channels, respectively. Figure 9(a) shows the demodulated results without any compensation. The four-polarization channels were demodulated using the PGC-DCM algorithm directly and the results are colored in red, blue, yellow and green, respectively. The results from the PS method are colored in black. Each polarization channel contains a factor related to fiber birefringence, as indicated in Eq. (10), thus the demodulated results of the four polarization channels differed from each other with their demodulated amplitudes at 1kHz being −18.35dB, −26.6dB, −30.11dB and −41.43dB (0dB re 1rad/ Hz1/2), respectively. The PS-PGC method delivered a signal amplitude of −32.16dB. Figure 9(b) shows the demodulated results when the real-time phase delays were completely compensated. Significant improvements can be seen, showing more consistent demodulated signal amplitudes, which were −24.35dB, −24.74 dB, −24.8 dB and −24.87dB (0dB re 1rad/ Hz1/2), respectively. The PS-PGC hybrid method obtained a signal amplitude of −24.75dB in this case.

 figure: Fig. 9.

Fig. 9. Demodulated signal at 1kHz (a) with no compensation and (b)with completely compensated phase delay.

Download Full Size | PDF

As indicated in Eq. (14), the fiber birefringence factor ${U_{MNn}}$ cannot be eliminated completely when there are real-time phase delays in the four polarization channels. As a result, the demodulated signal amplitude will fluctuate with environmental disturbances, i.e., polarization induced signal fading. We sampled the interferences of S1 for 87 times, and the demodulated signal amplitudes at 1kHz with and without real-time phase delay compensation are shown in Fig. 10. Without compensation, the demodulated signal amplitude fluctuates from −32dB to −17dB with the PS-PGC hybrid method, indicating that the PS method was invalid. When the real-time phase delays were completely compensated, the demodulated signal amplitude was −24.75dB with a fluctuation less than ± 0.75dB. The fiber disturbance-induced polarization did not affect signal demodulation and the PS-PGC hybrid method can reduce the effect of the polarization-induced signal fading well.

 figure: Fig. 10.

Fig. 10. Demodulated signal amplitude fluctuations with and without real-time phase delay compensation.

Download Full Size | PDF

Here the real-time phase delays were compensated using fixed parameters demodulated from the sampled interference, and Fig. 7 indicates that there was fluctuation when tested for a long period. For one single frame data, the compensation parameter may deviate from the true value. To evaluate the case when the real-time phase delay is not compensated completely, the demodulation results were tested when the compensation parameter deviated from the true value, as shown in Fig. 11. The x-axis shows the deviation percentage from the true value, and the y-axis is the amplitude fluctuation of the corresponding demodulation signal. Figure 11 indicates that the PS-PGC hybrid method can deliver an acceptable signal amplitude fluctuation less than ± 1.5dB when the compensation value deviation is less than 20%. It should be noted that if the signal amplitude fluctuation is larger than ± 1.5dB, the sensor cannot be used in practical environments.

 figure: Fig. 11.

Fig. 11. Demodulated signal amplitude fluctuations versus phase delay compensation deviation.

Download Full Size | PDF

Figure 12 shows the results of S2 with and without real-time phase delays. As the signal at 1kHz was applied only to S1 and the FBG reflectivities were controlled to be less than 1%, there was no obvious signal and the crosstalk was lower than 40dB. Figure 12(a) shows the demodulated results of S2 for a certain frame data with and without real-time phase delay compensation, colored in black and red, respectively. It can be found that if there was no compensation, the background noise was higher than that with complete compensation. As the phase background usually fluctuates with environmental disturbances, the result from a certain single frame data may not be convincing. Therefore, we demodulated S2 for 87times and the average results are shown in Fig. 12(b). It is clear that with the real-time phase delay compensation, the background noise level can be lowered by about 5dB than that without compensation. However, the difference was not obvious at frequency bands lower than 300Hz, which can be explained by the fact that the noise induced by environmental disturbances is much greater than by fiber-birefringence-induced polarization fading.

 figure: Fig. 12.

Fig. 12. Demodulated noise background of S2 for (a) single frame data and (b) average results accumulated for 87 times.

Download Full Size | PDF

4. Conclusions

This paper focused on performance promotion of interferometric FBG sensor arrays in which the real-time phase delay was analyzed in detail. As the PS-PGC hybrid processing method was used for sensor demodulation, the real-time phase delay must be compensated as both the demodulation accuracy and background noise could be deteriorated. The method to obtain the real-time phase delay was presented and experimental results of all polarization channels were obtained. For most of the time, the demodulated real-time phase delay remained a fixed value but fluctuations existed when the corresponding sensor phases were close to zero. There was a fixed difference between two adjacent polarization channels which was mainly caused by the interrogation time delay. As indicated by theoretical derivation, the real-time phase delay mainly affects the ICS and introduces error to the PS-PGC algorithm. The compensation method was thus presented. Experimental results showed that the compensation method could greatly improve the stability of the demodulated signal and suppress the sensor background phase noise. The demodulated signal amplitude was stabilized with a fluctuation less than ± 0.75dB and a 5dB noise suppression was observed. An acceptable compensation error of 20% was also derived by evaluating the demodulated signal amplitude. Research on real-time phase delay can facilitate more stable sensor performance to interferometric FBG sensor arrays and thus greatly promote their applications in harsh environments in practice.

Funding

National Natural Science Foundation of China (61901488).

Disclosures

The authors declare no conflicts of interest.

References

1. C. Okawara and K. Saijyou, “Fiber optic interferometric hydrophone using fiber Bragg grating with time division multiplexing,” Acoust. Sci. & Tech. 28(1), 39–42 (2007). [CrossRef]  

2. H. Nakstad and J. T. Kringlebotn, “Realization of a full-scale fibre optic ocean bottom seismic system,” Proc. SPIE 7004, 700436 (2008). [CrossRef]  

3. J. T. Kringlebotn, H. Nakstad, and M. Eriksrud, “Fibre optic bottom seismic system: from innovation to commercial success,” Proc. SPIE 7503, 75037U (2009). [CrossRef]  

4. E. Rnnekleiv, O. H. Waagaard, and H. Nakstad, “Ocean bottom seismic sensing system,” U.S. patent 7366055 (2008).

5. G. A. Cranch, G. M. H. Flockhart, and C. K. Kirkendall, “Efficient fiber Bragg grating Bragg grating and fiber Fabry–Pérot sensor multiplexing scheme using a broadband pulsed mode-locked laser,” J. Lightwave Technol. 23(11), 3798–3807 (2005). [CrossRef]  

6. O. H. Waagaard, E. Rønnekleiv, and S. Ford, “Reduction of crosstalk in inline sensor arrays using inverse scattering,” Proc. SPIE 7004, 70044Z (2008). [CrossRef]  

7. L. Ma, J. Wang, Z. Hu, X. Chen, P. Gan, and Y. Hu, “An ultra-low crosstalk and polarization independent inline interferometer fiber Bragg grating sensor array,” OFS 26, ThE77 (2018).

8. M. Ni, H. Yang, and S. Xiong, “Investigation of polarization-induced fading in fiber-optic interferometers with polarizer-based polarization diversity receivers,” Appl. Opt. 45(11), 2387–2390 (2006). [CrossRef]  

9. A. D. Kersey and M. J. Marrone, “Input-polarization scanning technique for overcoming polarization-induced signal fading in interferometric fiber sensors,” Electron. Lett. 24(15), 931–933 (1988). [CrossRef]  

10. L. A. Ferreira, J. L. Santos, and F. Farahi, “Polarization-induced noise in a fiber-optic Michelson interferometer with Faraday rotator mirror elements,” Appl. Opt. 34(28), 6399–6402 (1995). [CrossRef]  

11. O. H. Waagaard and E. Rnnekleiv, “Method and apparatus for providing polarization insensitive signal processing for interferometric sensors,” U.S. patent 7081959B2 (2006).

12. P. Jiang, L. Ma, Z. Hu, and Y. Hu, “Low-Crosstalk and polarization-independent inline interferometric fiber sensor array based on fiber Bragg gratings,” J. Lightwave Technol. 34(18), 4232–4239 (2016). [CrossRef]  

13. C. McGarrity and D. A. Kackson, “Improvement on phase generated carrier technique for passive demodulation of miniature interferometric sensors,” Opt. Commun. 109(3-4), 246–248 (1994). [CrossRef]  

14. A. Dandridge, A. B. Tveten, and T. G. Giallorenzi, “Homodyne demodulation scheme for fiber optic sensors using phase generated carrier,” IEEE Trans. Microwave Theory Tech. 30(10), 1635–1641 (1982). [CrossRef]  

15. S. Zhang, L. Yan, B. Chen, Z. Xu, and J. Xie, “Real-time phase delay compensation of PGC demodulation in sinusoidal phase-modulation interferometer for nanometer displacement measurement,” Opt. Express 25(1), 472–485 (2017). [CrossRef]  

16. N. N. Aleksandr, Y. P. Mikhail, V. V. Anton, V. M. Mikhail, and Y. K. Alexksadr, “PGC-Atan demodulation scheme with the carrier phase delay compensation for fiber-optic informetric sensors,” IEEE Sens. J. 18(5), 1985–1992 (2018). [CrossRef]  

17. S. Zhang, B. Chen, L. Yan, and Z. Xu, “Real-time normalization and nonlinearity evaluation methods of PGC-Artan demodulation in an EOM-based sinusoidal phase modulating interferometer,” Opt. Express 26(2), 605–616 (2018). [CrossRef]  

18. O. Sasaki and H. Okazaki, “Sinusoidal phase modulating interferometry for surface profile measurement,” Appl. Opt. 25(18), 3137–3140 (1986). [CrossRef]  

19. W. Mikhael and S. Michael, “Composite operational amplifiers: Generation and finite-gain applications,” IEEE Trans. Circuits Syst. 34(5), 449–460 (1987). [CrossRef]  

20. T. Lan, C. Zhang, L. Li, G. Luo, and C. Li, “Carrier phase advance technique for digital PGC demodulation,” Opto-Electronic Eng. 35(7), 49–52 (2008).

21. K. Wang, M. Zhang, F. Duan, S. Xie, and Y. Liao, “Measurement of the phase shift between intensity and frequency modulations within DFB-LD and its influences on PGC demodulation in a fiber-optic sensor system,” Appl. Opt. 52(29), 7194–7199 (2013). [CrossRef]  

22. S. C. Huang and H. Lin, “Modified phase-generated carrier demodulation compensated for the propagation delay of the fiber,” Appl. Opt. 46(31), 7594–7603 (2007). [CrossRef]  

23. P. Jiang, L. Ma, Z. Hu, and Y. Hu, “An 8-TDM inline fiber Fabry-Perot sensor array based on ultra-weak fiber Bragg gratings,” Asia Communications and Photonics Conference (ACP), AM1D.4 (2015).

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 (12)

Fig. 1.
Fig. 1. Typical structure of interferometric FBG sensor array.
Fig. 2.
Fig. 2. Interrogation and interference pulses.
Fig. 3.
Fig. 3. DCM procedure.
Fig. 4.
Fig. 4. Schematic diagram of experimental system.
Fig. 5.
Fig. 5. Measured result of reflected pulse train.
Fig. 6.
Fig. 6. Sampled interferences of (a)S1 and (b) S2 in the frequency domain.
Fig. 7.
Fig. 7. Real-time phase delay of (a)S1 and (b) S2.
Fig. 8.
Fig. 8. Real-time phase delay of S1 versus demodulated phase.
Fig. 9.
Fig. 9. Demodulated signal at 1kHz (a) with no compensation and (b)with completely compensated phase delay.
Fig. 10.
Fig. 10. Demodulated signal amplitude fluctuations with and without real-time phase delay compensation.
Fig. 11.
Fig. 11. Demodulated signal amplitude fluctuations versus phase delay compensation deviation.
Fig. 12.
Fig. 12. Demodulated noise background of S2 for (a) single frame data and (b) average results accumulated for 87 times.

Tables (2)

Tables Icon

Table 1. Parameters of FBG sensors array

Tables Icon

Table 2. Parameters of interrogation system

Equations (17)

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

E r 1 = r ρ 0 E i n 1 E s 1 = t 0 B 1 r ρ 1 B 1 t 0 E i n 0 E r 2 = t 0 B 1 r ρ 1 B 1 t 0 E i n 1 E s 2 = t 0 B 1 t 1 B 2 r ρ 2 B 2 t 1 B 1 t 0 E i n 0
B n B n = B n T B n = ( 1 0 0 1 ) , B n T = r B n r
E r 1 = r ρ 0 E i n 1 E s 1 = t 0 2 ρ 1 r B 1 T B 1 E i n 0 E r 2 = t 0 2 ρ 1 r B 1 T B 1 E i n 1 E s 2 = t 0 2 t 1 2 ρ 2 r B 1 T B 2 T B 2 B 1 E i n 0
I 1 = ( E r 1 + E s 1 ) ( E r 1 + E s 1 ) = D C 1 + 2 t 0 2 ρ 0 ρ 1 Re ( E i n 1 B 1 T B 1 E i n 0 ) = D C 1 + 2 t 0 2 ρ 0 ρ 1 exp ( i φ s 1 ) Re ( E i n 1 U 1 T U 1 E i n 0 ) I 2 = ( E r 2 + E s 2 ) ( E r 2 + E s 2 ) = D C 2 + 2 t 0 4 t 1 2 ρ 1 ρ 2 Re ( E i n 1 B 1 B 2 T B 2 B 1 E i n 0 ) = D C 2 + 2 t 0 4 t 1 2 ρ 1 ρ 2 exp ( i φ s 2 ) ( E i n 1 U 1 U 2 T U 2 U 1 E i n 0 )
1 = exp ( i φ s 1 ) U 1 T U 1 = exp ( i φ s 1 ) U 1 2 = exp ( i φ s 2 ) U 1 U 2 T U 2 U 1 = exp ( i φ s 2 ) U 2
E i n 0 = [ 1   0 ] T E i n 1 = [ 1   0 ] T exp [  -  i C cos ( 2 π f 0 t ) ]
I X X n = I n _ D C + t n { exp ( i φ s n ) exp [ i C cos ( 2 π f 0 t ) ] U 11 n } ( n = 1 , 2 )
E i n 0 = [ 1   0 ] T E i n 1 = [ 0   1 ] T exp [  -  i C cos ( 2 π f 0 t ) ] E i n 0 = [ 0   1 ] T E i n 1 = [ 1   0 ] T exp [  -  i C cos ( 2 π f 0 t ) ] E i n 0 = [ 0 1 ] T E i n 1 = [ 0   1 ] T exp [  -  i C cos ( 2 π f 0 t ) ]
I X Y n = I n _ D C + t n { exp ( i φ s n ) exp [ i C cos ( 2 π f 0 t ) ] U 12 n } I Y X n = I n _ D C + t n { exp ( i φ s n ) exp [ i C cos ( 2 π f 0 t ) ] U 21 n } I Y Y n = I n _ D C + t n { exp ( i φ s n ) exp [ i C cos ( 2 π f 0 t ) ] U 22 n }
X 1 M N n = t n J 1 ( C ) U r M N n sin ( φ s n + U a M N n ) X 2 M N n = t n J 2 ( C ) U r M N n cos ( φ s n + U a M N n ) ( M N = 11 , 12 , 21 , 22   )
M N n = X 2 M N n J 2 ( C ) i X 1 M N n J 1 ( C ) = t n U r M N n exp ( i φ s n i U a M N n ) = t n exp ( i φ s n ) U M N n
X X n × Y Y n X Y n × Y X n = det n = exp ( i 2 φ s n )
I X X n = I n _ D C + t n { e x p ( i φ s n ) e x p [ i C cos ( 2 π f 0 t  +  θ ) ] U 11 } I X Y n = I n _ D C + t n { e x p ( i φ s n ) exp [ i C cos ( 2 π f 0 t  +  θ ) ] U 12 } I Y X n = I n _ D C + t n { exp ( i φ s n ) exp [ i C cos ( 2 π f 0 t  +  θ ) ] U 21 } I Y Y n = I n _ D C + t n { exp ( i φ s n ) exp [ i C cos ( 2 π f 0 t  +  θ ) ] U 22 }
X 1 M N n = t n J 1 ( C ) cos ( θ ) U r M N n sin ( φ s n + U a M N n ) X 2 M N n = t n J 2 ( C ) cos ( 2 θ ) U r M N n cos ( φ s n + U a M N n )
M N n = X 2 M N n J 2 ( C ) cos ( 2 θ ) i X 1 M N n J 1 ( C ) cos ( θ ) = t n U r M N n exp ( i φ s n i U a M N n ) = t n exp ( i φ s n ) U M N n
X 1 M N n = t n J 1 ( C ) cos ( θ ) U r M N n sin ( φ s n + U a M N n ) X 2 M N n = t n J 1 ( C ) sin ( θ ) U r M N n sin ( φ s n + U a M N n )
θ  = arc tan( X 2 M N n X 1 M N n )
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.