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

A fully digital MIMO-OFDM scheme for fading mitigation in coherent Δϕ-OTDR

Open Access Open Access

Abstract

We introduce a coherent-MIMO Δϕ-OTDR based sensing technique that uses digital frequency division multiplexing of probing codes in order to mitigate the coherent fading issue in Δϕ-OTDR. Orthogonal frequency division multiplexing (OFDM) is used to probe the optical channel on several electrical subcarriers, and their responses are then combined according to specific reliability criterion to enhance the overall sensing sensitivity. This allows for a decrease of the noise floor in sensing measurements on standard single mode fiber (SSMF) and the mitigation of false alarms. Coherent MIMO-OFDM enables sensing over telecom fibers with a measured sensitivity level below $21\mathrm {p\varepsilon /\sqrt {Hz}}$ along 1.3km fiber over a mechanical bandwidth of 760Hz. Enhancements in the dynamic range for detection of mechanical excitation are presented with a gain over 18dB as compared to single-carrier probing.

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

1. Introduction

Distributed optical fiber sensors (DOFS) are able to perform strain measurements on long distance as well as with very precise gauge length [1], therefore they are of interest in a growing number of fields, from the biomedical domain [2] to oil and gas industry [3] or smart-city applications [4]. In particular, differential phase-sensitive optical time-domain reflectometry ($\Delta \phi$-OTDR) is a powerful technique for detection, localization and identification of events of all scales, based on Rayleigh backscaterred light. $\Delta \phi$-OTDR can be performed on top of any kind of optical fiber, including deployed telecom cables. Generally, highly coherent laser sources are used to probe the sensor, together with coherent detection at the reception to retrieve the phase information.

Coherent phase-sensitive OTDR sensors are subject to three main impairments: polarization fading, coherent fading, and laser phase noise. The latter issue depends on the laser properties and specific methods were developed [5,6] that are not further addressed here. Polarization fading is due to polarization misalignment at the transmitter and receiver, and it is addressed mainly by polarization diversity methods [7,8], or can be solved by Coherent-MIMO sensing which consists in using polarization diversity both at the transmitter and at the receiver [9,10]. On the other hand, coherent fading is inherent to coherent $\phi$-OTDR and $\Delta \phi$-OTDR since the technique relies on the speckle pattern resulting from the interrogation of the fiber using a highly coherent source [11]. As a main drawback, the sensitivity randomly fluctuates along the fiber sensor and in time. This impairment is also referred to as Rayleigh fading, interference fading, signal fading, fading noise or speckle noise, as it exists in various fields of application. In the following, we keep the appellation "coherent fading". The impact of coherent fading can be limited either thanks to trace averaging [12] at a cost of reduced mechanical sensing bandwidth, or by the use of wavelength- and frequency-diversity methods. Wavelength-diversity methods rely on the independence that exists between the Rayleigh patterns of the same optical fiber that is interrogated at different wavelengths [13] with a $50$GHz wavelength spacing, thus requiring multiplexing and demultiplexing several optical carriers at the transmitter and receiver. Frequency-diversity methods require smaller spacing, around hundreds of MHz [1416]. There are two main types of frequency-diversity interrogation: multiple pulses with different carriers sent to probe the fiber sensor [17]; single pulses with frequency chirping (linear frequency modulation, LFM) [18], from which the channel response of the fiber sensor is recovered by correlation at the receiver side. In $\phi$-OTDR, the introduction of a chirped pulse method can also allow to perform quantitative strain measurement [19] based on the amplitude of $\phi$-OTDR traces, thus avoiding some of the phase-related impairments. Recent works push forward the sensor optimization as they use a series of frequency-shifted interrogation pulses by combining frequency and polarization multiplexing in $\phi$-OTDR [20,21], thus injecting more individual pulses that would be allowed by the round-trip time of light in the fiber sensor and increasing the spectral occupancy without any impact on the probing duration. The associated optical setups may be complex which make them difficult to deploy for field tests outside the laboratories [22]. Instead, we choose here to tackle coherent fading by exploiting the Rayleigh phase pattern diversity. Our approach involves fully digital Orthogonal Frequency Division Multiplexing (OFDM) of probing codes to achieve carrier phase diversity, allowing for coherent fading mitigation while keeping a simple optical setup. First experimental results were displayed in [23].

The paper is organized as follows. In section 2., we remind the basic principles of coherent $\Delta \phi$-OTDR, the reasons for Rayleigh fading and its impact on sensing accuracy. Then, the novel carrier phase diversity channel estimation method used in this study is presented in section 3.. The method for mitigation of coherent fading is given in section 4.. Finally, section 5. provides an experimental validation of the study.

2. Coded coherent $\Delta \phi$-OTDR

2.1 Rayleigh backscattering and coherent fading

Distributed fiber sensing involves processing of the backscattered light from a probed optical fiber. The sensing signal is a summation of the contributions of the backscattered light from the probing signal $E_{in}$ all along the fiber sensor. Therefore, when $E_{in}$ consists of a rectangular pulse with an amplitude $E_0$ and a duration $T_S$, the complex optical field at the receiver side $E_{rx}$ can be expressed as:

$$E_{rx}(t) = \sum_{m=0}^{M} A_m e^{j \phi_m} e^{{-}2\alpha z_m+j\left[\omega_0(t-t_m) - 2\beta z_m\right]} E_{0}\mathrm{rect}(\frac{t-t_m}{T_S})$$
where $M$ is the number of fiber segments defined by the probe symbol duration $T_S$ and the fiber length, $A_m$ and $\phi _m$ are respectively the Rayleigh-distributed backscattered amplitude and the backscattered phase from the $m^{\mathrm {th}}$ fiber segment [24]. $\alpha$ denotes the attenuation coefficient along the fiber and $\beta$ is the propagation constant. $z_m$ is the distance of segment $m$ from the fiber start, therefore the round-trip time of the pulse traveling up to the $m^{\mathrm {th}}$ segment is $t_m =\frac {2z_m c}{n_g}$ where and $c/n_g$ is the speed of light in the fiber and $n_g$ is the refractive group index. $\omega _0$ is the optical carrier frequency (the probing signal is assumed to be quasi-monochromatic).

Each $A_m e^{j\phi _m}$ term in 1 results from the sum of elementary contributions within the $m^{\mathrm {th}}$ segment:

$$A_m e^{j\phi_m} = \sum_{i=0}^{N_{scat}} a(z_i) e^{j \varphi(z_i)}$$
where $N_{scat}$ is the number of elementary scatterers within segment $m$. Each scatterer is located at a distance $z_i$ from the fiber start, as displayed in Fig. 1, backscatters an amplitude $a(z_i)$ and rotates the phase of incident light by $\varphi (z_i)$. The model can be extended to a vector-format (for dual-polarization) using the Jones formalism, where the backscattered phase and amplitude are distributed over a $2\times 2$ Jones matrix $\mathbf {H}_m$ describing the evolution of the state of polarization (SOP) of the transmitted probe [10], as depicted in Fig. 1.

 figure: Fig. 1.

Fig. 1. Backscattering model: $M$ segments with $N_{scat}$ scatterers within each segment. $\mathbf {H}_m$ is the $2\times 2$ Jones matrix describing the SOP evolution up to segment $m$.

Download Full Size | PDF

$\phi$-OTDR and $\Delta \phi$-OTDR techniques rely on the change of interference figure $E_{rx}$ when a vibration occurs. This effect allows to precisely retrieve the phase variation induced by the strain or vibration applied to the fiber. However, coherent fading happens when the summations in 2 interfere destructively, resulting in $A_m e^{j \phi _m} \rightarrow 0$ locally. Also, coherent fading occurrences are independent of the probing method used for interrogating the fiber sensor.

2.2 Coded interrogation

The choice of the probe signal for a fiber sensor depends on the system requirements such as the maximum length of the probed fiber, the spatial resolution and the desired mechanical bandwidth.

Single light pulses or pulse trains are the most common probe techniques and they allow to directly capture the backscattered response of the fiber sensor [25] ; spatial resolution is tuned by adjusting the pulse width $T_S$. We denote the spatial resolution as $L_S=cT_S/2n_g$. The covered mechanical bandwidth $B_{meca}$ is determined by the pulse repetition rate. However, there is a trade-off between signal to noise ratio in the longer distances and spatial resolution (when increasing the pulse temporal width) or non-linearity (when increasing the peak power). Frequency sweep [26] or multiple-frequency interrogation [20], linear frequency modulation [27], as well as coded interrogation [8,9,28] address this issue by spreading the probe power in time; a narrow temporal pulse is recovered by correlation at the receiver. Coded interrogation is especially interesting as it is digitally generated thus robust to frequency distortion effects [26], and with low Peak-to-Average Power Ratio (PAPR), thus very much suited for coexistence with telecommunication data, when applied over lit fibers. With a probing code of length $T_{code}$, the achieved mechanical bandwidth is $B_{meca}=1/2T_{code}$.

Coherent-MIMO [9,10] was proposed in that respect, with a multiple-input multiple-output (MIMO) approach through polarization multiplexing at the transmitter and polarization-diversity detection at the receiver side. It allows to retrieve not only the phase and amplitude along the fiber sensor but also the evolution of the backscattered state of polarization in the form of $2\times 2$ Jones matrices $\textbf {H}_{t,d}$ (where $t$ is the slow time index and $d$ distance index (fast time)). This probing technique was proven to fully solve the polarization fading issue [10]. The probing codes consist of two mutually orthogonal complementary pairs $\{G_{a1}, G_{b1}\}$ and $\{G_{a2}, G_{b2}\}$ derived from Golay codes, mapped to binary-phase-shift-keying (BPSK) symbols that jointly and continuously modulate two orthogonal polarization states of a narrow-linewidth laser source [9]. Figure 2 depicts the general scheme of Coherent-MIMO interrogation.

 figure: Fig. 2.

Fig. 2. Pol.: polarization, MZM: Mach Zehnder Modulator, LO: Local Oscillator, DSP: Digital Signal Processing. $T_S$: symbol duration, $T_{code}$: probing code duration.

Download Full Size | PDF

Yet, any $\Delta \phi$-OTDR probing technique is still subject to coherent fading issues as described in subsection 2.1, whatever fine resolution, probing bandwidth are chosen. Hence, coherent fading mitigation without damage to the phase information is needed. The most interesting method yet is to use multiple uncorrelated channels to probe the line: the channels being either frequency channels [14] or different modes of an optical fiber [29], which in both cases will make the most of uncorrelated backscattered amplitudes and phases per channel. Therefore, we will adapt the Coherent-MIMO probing technique to solve the coherent fading issue.

3. Channel estimation scheme to mitigate coherent fading

An electrical bandwidth $B_{elec}$ around $100$MHz, fixed by the symbol rate $f_{Symb}$ and the pulse shaping, is usually used to sense a fiber. Over such a narrow bandwidth, the elementary scatterers have a flat amplitude spectrum response, and wavelength diversity cannot be achieved. However, to reach statistically independent intensity fading, the fiber should be probed with signals separated by $\Delta f \geq f_{Symb}$ [14,16,30] where $f_{Symb} = 1/T_S$. Therefore, we introduce an orthogonal frequency division multiplexing (OFDM) scheme to independently probe the line over several orthogonal subcarriers, verifying $\Delta f = f_{Symb}$ inside the electrical bandwidth $B_{elec}$.

The dual-polarization codes $\mathbf {E}_{in} = [\begin {smallmatrix} E_x \\ E_y \end {smallmatrix}]$ of length $N_{code}$ symbols for single-carrier Coherent-MIMO sensing are mapped onto $N_{sc}$ subcarriers. Figure 3 shows a simple example for one polarization axis and two subcarriers, $N_{sc}=2$. The subcarriers carry the same code $\{G_{a1}, G_{b1}\}$ for $N_{code}=8$. For sake of simplicity, we limit our study in this paper to the case where all $N_{sc}$ subcarriers enclose the same code as we introduce this digital frequency multiplexing technique for the first time in a $\Delta \phi$-OTDR scheme. The OFDM signal is created by serial-to-parallel, inverse fast Fourier transform (iFFT), and parallel-to-serial operations. No cyclic prefix is required as temporal spread induced by chromatic dispersion is too small compared to the symbol rate.

The OFDM code alphabet is limited here to $\{-1,0,1\}$, and the symbol rate of the OFDM flow combining $N_{sc}$ subcarriers is equal to $f_{symb, OFDM} = N_{sc}\times f_{symb}$, thus expanding the electrical bandwidth when the symbol rate is left unchanged. Note that this is a simplified case of OFDM interrogation where all probing subcarriers carry the same code. To better motivate the use of OFDM interrogation, other schemes could be proposed, where different codes could be loaded onto different subcarriers. This would still allow to probe a fiber sensor with orthogonal subcarriers, only with a more complex OFDM alphabet at the modulator side. In such various configurations, different codes with possibly different lengths may modulate the subcarriers. This would achieve other types of diversity, eg. in terms of mechanical bandwidth. In that particular configuration using the same code over the $N_{sc}$ subcarriers, the OFDM interrogation amounts to performing an interrogation of the fiber sensor with a $N_{sc}$-fold expanded version of the original code (done by inserting $N_{sc}-1$ zeros between symbols). Such a configuration maximizes the PAPR of the interrogation signal compared with other generic OFDM implementations, since the mean input power is kept constant and therefore it will end up being concentrated into the non-zero symbols. Yet, such PAPR is still negligible compared with pulsed interrogation as the power is still spread in time over multiple (OFDM) symbols, thus limiting any nonlinearity-related issues. Furthermore, the OFDM formalism is especially suitable for the linear combination of the backscattered signals at the receiver side through the demultiplexing of the carriers by means of an FFT operation.

 figure: Fig. 3.

Fig. 3. OFDM combination for two subcarriers SC1 and SC2, one code duration, one polarization

Download Full Size | PDF

Our purpose is to compare the sensing performances of the OFDM probing scheme with our single carrier reference case. For sake of fairness, the comparison is applied over the same bandwidth $B_{elec}$, which is achieved by mapping subcarriers with codes shorter by a factor of $N_{sc}$ compared with the single carrier case. In such a way, the code duration $T_{code}$ is left unchanged and so is the sensed mechanical bandwidth $B_{meca}$. Finally, the OFDM setup extension allows to perform $N_{sc}$ channel estimations at the sole cost of a spatial resolution reduction by a factor of $N_{sc}$ compared with the single-carrier coherent-MIMO case and with constant $B_{elec}$. The obtained estimations of Jones matrices are here again mathematically perfect [9], the overall spectral occupancy of the modulation is preserved and the mechanical bandwidth of the analysis also. The next step aims at combining the subcarrier estimations to minimize the Rayleigh fading effect.

4. Digital carrier phase diversity combining for coherent fading reduction

The subcarrier information is constructively combined in phase, independently at each fiber segment, by means of independent rotations at each subcarrier. This technique, known as "Rotated Vector Sum" was successfully used in [18] for combining phase estimates from independent optical channels, bringing a substantial SNR increase compared to a single channel sensing scheme.

4.1 Constructive combination

MIMO-OFDM interrogation yields Jones matrices responses for each time instant $t$, segment at distance $d$ and subcarrier index $n$: $\textbf {H}_{t,d,n}$. The backscattered SOP (b-SOP, evolution of the energy on the four terms of the Jones matrix) is common to all subcarriers at any time and distance, thus it is not involved in the combination of subcarrier information. However, the intensity of the backscattered field from some subcarriers within the sensing signal is potentially degraded in practice because of coherent fading, in which case the most accurate b-SOP is the one yielding the highest backscattered intensity among all subcarriers.

The independent phase and amplitude information is therefore enclosed in the sole Jones matrices determinants $\det ( \textbf {H}_{t,d,n})$ from which the local cumulated phase is extracted as $\phi _{t,d,n} = \frac {1}{2} \angle {\det ( \textbf {H}_{t,d,n})}$, and the reliability metric derived from the backscattered intensity as $R_{t,d,n} = |\det ( \textbf {H}_{t,d,n})|$ [31]. Consequently, we choose to constructively combine the subcarriers responses from their respective Jones matrix determinants $\det ( \textbf {H}_{t,d,n})$.

The rotated-vector-sum combination method is applied to the $D_{t,d,n} = \det ( \textbf {H}_{t,d,n})$ vectors, so as to maximize the combined modulus [18]. For a given segment index $d$, the time-average of phase $\overline {\phi }_{d,n}$ is derived for each subcarrier $n$ and subtracted from the phase of all $D_{t,d,n}$, such that the resulting combined determinants are:

$$D_{t,d,combined} = \sum_{n=1}^{N_{sc}} \dfrac{D_{t,d,n} }{\exp(j\overline{\phi}_{d,n})}$$
The final differential phase between adjacent segments is obtained by subtracting the phases of adjacent $D_{t,d,combined}$, such that $\Delta \phi _{t,d,combined} = \frac {1}{2} (\angle {D_{t,d,combined}}-\angle {D_{t,d-1,combined}})$. The operation in 3 consists in removing the average cumulated phase values $\overline {\phi }_{d,n}$ and keeping the sole variations of the cumulated phase in time, per segment, since the relevant information is enclosed in the differential phases eventually. On the one hand, if the fiber segment is static (no mechanical perturbation), then the random variations of phase around the $\overline {\phi }_{d,n}$ are independent and they will fade after combination. On the other hand, if the fiber segment is mechanically excited, then all subcarriers are affected the same way such that their oscillations around the phase references $\overline {\phi }_{d,n}$ align, so magnifying the overall captured excitation. Moreover, we recently demonstrated in [31] that $|D_{t,d,n}|$ is not only the most accurate intensity estimator in MIMO sensing up to our knowledge but also a reliability criterion that informs about the ability to extract a fair phase estimation from any Jones matrix. Therefore, the combination in 3 is naturally weighted by the reliability metric $R_{t,d,n}$, hence a lower influence of the unreliable estimations on the final result.

4.2 Individual carriers versus their combination: simulation

To detect mechanical events along the sensor, the differential phase standard deviation (denoted $\sigma _\varphi$, expressed in radians) is monitored along distance: theoretically zero-valued in absence of events, and non-zero positive in case of any detected disturbance. The $\sigma _\varphi$ from static measurements indicates the threshold value or sensitivity level below which an event cannot be detected. Also, the normalized reliability metric $R_{t,d,n}/\underset {t,d,n}{\max }(R_{t,d,n})$ gives a priori information on the soundness of the estimation [31]. Its minimum value in time per segment index $d$ gives an indication on the probability to get local unstable phases estimations giving rise to false alarms.

We simulate a 4-subcarrier OFDM interrogation of a static $1$km long single mode fiber (SMF), with laser noise $\Delta \nu = 75$Hz and additive gaussian white noise (AWGN) at the receiver (shot noise and thermal noise), we set $B_{meca}=760$Hz and $B_{elec}=50$MHz. The fiber is kept short such that the attenuation with distance does not participate to the fading effect, thus focusing the study on the sole coherent fading effect. Every separate subcarrier is monitored, alongside with the subcarrier combination. Figure 4 illustrates the enhancement brought by an OFDM interrogation.

 figure: Fig. 4.

Fig. 4. Simulation OFDM4, 1km SMF, $\Delta \nu = 75$Hz, AWGN at the receiver, individual subcarrier responses compared with combined response

Download Full Size | PDF

Figure 4(a) gives the minimum reliability metric per segment versus distance, and shows that the reliability increases after subcarrier combination compared with individual subcarriers where it locally dives towards $0$. In other words, the worst case points are avoided. Figure 4(b) shows $\sigma _\varphi$ along distance, giving the noise floor above which a mechanical event can be detected. The combined measurement presents a noise floor below $0.01$rad all along the $1$km fiber, which is locally up to 7 times lower than the level given by a single subcarrier (subcarrier 3, around $850$m).

Notice that the combined $\sigma _\varphi$ in Fig. 4(b) reaches values that are even lower than the lowest $\sigma _\varphi$ per subcarrier at several locations ($300$m, $500$m, $880$m). This shows that in addition to avoiding the worst-case sensitivity situations, the OFDM combination scheme performs an averaging on the system noise and thus significantly enhances the overall sensitivity.

4.3 Impact of noises and OFDM contribution

The sensitivity gain brought by the OFDM approach is first quantified by simulation: the interrogations of a static, 1km-long SMF are simulated. Single-carrier, OFDM with 2 subcarriers and OFDM with 4 subcarriers interrogations are performed. The mechanical bandwidth is fixed to $B_{meca} = 760$Hz for all of the following measurements and simulations. We introduce two metrics to help fairly compare the sensing performance as a function of the number of subcarriers. $\mathrm {Intensity}_{dB}$ and SNR$_{\varepsilon,dB}$ are derived as follows:

$$\begin{aligned}\mathrm{Intensity}_{dB} &= 10\log_{10}(|\det(\textbf{H})|)\\ \mathrm{SNR}_{\varepsilon,dB} &= 20 \log_{10}(\dfrac{4\pi n \xi G}{\lambda \sigma_\varphi} )\end{aligned}$$
with $n$ the refractive index of the fiber, $\xi$ the photo-elastic scaling factor in isotropic material, $\lambda$ the laser wavelength of the interrogator expressed in nm, $G$ the gauge length is the system, in m. The gauge length $G$ definition is based on the receiver characteristics which define the performance of the interrogation system [32], therefore $G=cT_{S}/2n$ with $1/T_S$ the symbol rate of the system, here $G=2,05$m. The former metric is the intensity estimated from $|\det (\textbf {H})|$ which, as emphasized earlier, informs about the ability to reliably derive a phase estimate from H. The latter one, SNR$_{\varepsilon,dB}$, expresses the strain sensitivity following the relation $\sigma _\varepsilon = \dfrac {\lambda \sigma _\varphi }{4 \pi n \xi G}$ [32] in n$\varepsilon$, where the standard deviation of phase in time $\sigma _\varphi$ expresses the noise floor of the system (in radians). Therefore $N_\varepsilon = \frac {\sigma _\varepsilon ^2}{B_{meca}}$ is the strain noise power density. We define SNR$_{\varepsilon,dB}= 10\log _{10}(1/\sigma _\varepsilon ^2)$ to express the strain noise in the absence of signal (assuming a reference strain signal of power $1n\varepsilon ^2$). The maximum value of $\sigma _\varepsilon$ (or equivalently minimum value of $\mathrm {SNR}_{\varepsilon }$) in a static measurement will be considered to be the strain sensitivity threshold of the interrogator. $\mathrm {SNR}_{\varepsilon }=0$dB corresponds to a $1\mathrm {n}\epsilon$ strain sensitivity.

A case with constant $B_{elec} = 50$MHz is considered, in which the spatial resolution decreases with the number of OFDM carriers according to $S_r = \{2.05\mathrm {m}, 4.10\mathrm {m}, 8.20\mathrm {m}\}$ for single-carrier, OFDM2 and OFDM4 respectively. The laser used for the following experimental measurements (section 5.) has an estimated linewidth of $\Delta \nu = 75$Hz over ms-range observation windows [33]. In the simulation, its phase noise is modeled as a Wiener process of variance $\sigma ^2 = 2\pi \Delta \nu T_S$ where $\Delta \nu$ is the Lorentzian laser linewidth and $T_S$ the duration of a symbol in the probing code. Fig. 5 displays the simulation results for intensity and strain sensitivity in presence of AWGN and laser phase noise. The observed probability density functions in Fig. 5 verify a tightening of the distributions as the number of subcarriers increases. Particularly, a decrease of the lower SNR values of the distribution is noticed thanks to OFDM, meaning that there are less low-energy and low-reliability fiber segments which are the most likely to trigger false alarms.

 figure: Fig. 5.

Fig. 5. SNR probability density distributions. 1km static simulation, laser phase noise $\Delta \nu = 75$Hz, AWGN at the receiver. $\mathrm {SNR}_{\varepsilon }=0$dB corresponds to $1\mathrm {n}\epsilon$ strain sensitivity.

Download Full Size | PDF

That is, the use of phase diversity through OFDM with an adapted number of subcarriers decreases the low SNR$_{\varepsilon,dB}$ occurrences, thus improving the SNR$_{\varepsilon,dB}$ by almost $20$dB with 4 subcarriers (Fig. 5(b)).

Note that although a higher number of subcarriers can be used for channel interrogation, the sensitivity gain tends to decrease after a certain number of carriers, under the fixed constraints of equivalent electrical and mechanical bandwidths. Figure 6 shows the evolution of the SNR$_{\varepsilon,dB}$ probability density function when the number of subcarriers increases in the particular case of $B_{elec}=50$MHz on a 1km simulated distance. After interrogation on 4 subcarriers, an increase of carriers does not bring a significant sensitivity gain (lower mean sensitivity in Fig. 6(b) with 8 subcarriers); this is due to the use of a decreasing code order (smaller $N_{code}$) when increasing the number of subcarriers, which results in a lower achieved SNR [9].

 figure: Fig. 6.

Fig. 6. SNR$_{\varepsilon,dB}$ on 1km, for 1 to 16 subcarriers, simulation with 5000 points per interrogation method, $\Delta \nu = 75$Hz and AWGN noise at the receiver

Download Full Size | PDF

5. Experimental measurements

The interrogating setup presented in Fig. 7 is based on our previous MIMO interrogation setup [33]: the probing codes are sent onto two orthogonal polarization axes through a dual-polarization $25$GHz I/Q Mach-Zehnder modulator to modulate the optical wavelength. Although Fig. 7 shows the general I/Q configuration, we restrict ourselves to a BPSK modulation per polarization in this work, as described previously in subsection 2.2.

OFDM sensing is introduced by inserting additional signal processing blocks before and after the optical setup, namely a MIMO-OFDM coding block at the probing codes generation stage performing iFFT (inverse fast Fourier transform) on parallel probing codes, and a MIMO-OFDM decoding block after coherent reception and analog-to-digital conversion to perform FFT on the OFDM coded subcarriers, before correlation and combination for fading mitigation. The core optical setup of the interrogator is thus left unchanged.

 figure: Fig. 7.

Fig. 7. OFDM-MIMO interrogation setup. MIMO: Multiple-Input/Multiple-Output, OFDM: Orthogonal Frequency Division Multiplexing, LO: Local Oscillator, ADC: Analog-to-Digital Converter, PZT: PiezoElectric Transducer

Download Full Size | PDF

5.1 Static experiments

The experiments are first conducted in static mode, in laboratory conditions, in an acoustically insulated environment. The block described as "telecom cable under test" in Fig. 7 consists in a $1300$m long single mode fiber (SSMF) made of the junction of a $400$m SSMF spool with an other $900$m fiber spool. Similarly to the simulations in subsection 4.3, a configuration with $B_{elec} = 50$MHz is considered. The mechanical bandwidth is also kept at $B_{meca} = 760$Hz, and the frequencies below $60$Hz are filtered out to avoid capturing residual acoustic and electric noise from the lab.

The measurements are successively single-carrier, OFDM2 and OFDM4 2-second interrogations of the same fiber with $f_{Symb}=50$MHz, yielding spatial resolutions of $2.05$m, $4.10$m and $8.20$m respectively. The lengths of the probing code are $N_{code} = 32768$ symbols for single carrier interrogation, $N_{code} = N_{code,single}/N_{SC}$ for OFDM interrogations, such that $T_{code}=655\mu s$ yielding the $B_{meca} = 760$Hz mechanical bandwidth in all three cases.

The performance of these interrogations is displayed in Fig. 8. In Fig. 8(a), we show the obtained intensity distributions. Note that the mean value of each distribution depends on the normalization at the receiver after analog to digital conversion, hence a different mean value from one measurement to another. The main observation here is the shrinking of the distributions as the number of subcarriers increases. The distributions of SNR$_{\varepsilon,dB}$ values in Fig. 8(b) behave as predicted in Fig. 5(b) as the distribution shrinks when increasing the subcarriers number. The worst case SNR$_{\varepsilon,dB}$ improves by over $11$dB when comparing single-carrier interrogation to OFDM4, yielding a minimal SNR$_{\varepsilon,dB}=4.61$ which corresponds to a $\sigma _\varepsilon =0.59$n$\epsilon$ strain within $B_{meca} = 760$Hz and therefore $21$p$\varepsilon /\sqrt {Hz}$ strain density sensitivity floor.

 figure: Fig. 8.

Fig. 8. Static experimental trace, OFDM interrogation of 1.3km SSMF, constant $B_{elec}=50$MHz configuration

Download Full Size | PDF

Thanks to subcarrier combination, the measured estimations of the fiber segments from Rayleigh-backscattered probing codes are more reliable, thus increasing the overall performance of intensity and phase estimations. Such a reduction of the noise floor (minimal SNR$_{\varepsilon,dB}$ values) is expected to decrease false alarm occurrences, together with an increase of the detection capabilities of the system in dynamic configurations.

5.2 Dynamic multicarrier sensing

A piezoelectric actuator is added to the fiber under test in the previously described laboratory configuration: $3.1$m of SSMF is wrapped around the actuator and connected between the $400$m an the $900$m fiber spools. The actuator is driven by a $440$Hz, $2$Vpp sine wave signal. Each measurement lasts for $2$ seconds similarly to the previous static interrogations.

Fig. 9 displays the experimental results for a constant $B_{elec}$ configuration: Fig. 9(a) predicts the reliability of the measurement along the fiber. The lower the minimal reliability values are, the higher the false alarm probability. Indeed, the low value at $1.14$km in Fig. 9(a) deserves some attention, as it translates into a false alarm in Fig. 9(b). Figure 9(b) shows the measured strain $\varepsilon$ in nm/m (or n$\epsilon$) along the fiber, used as the detection figure of mechanical events. The three compared measurements exhibit a significant peak around $420$m where the dynamic strain is applied, as expected. Also, the detection peak has a better dynamic range (defined as the ratio between false alarm level, namely highest peak in non-excited zones, and detection peak level) when the number of subcarriers increases, as the level of phase noise (false alarms) decreases thus gaining $3.67$n$\epsilon$ strain sensitivity in Fig. 9(b) between single carrier and OFDM4 cases thanks to OFDM coding and subcarrier combination. Figure 9(c) gives the optical phase measurement around the excitation detection in the frequency domain. Similarly to the detection stage, the identification gets better when the number of OFDM subcarriers increases, as the noise level on Fig. 9(c) shrinks by $15$dB between single-band interrogation and OFDM4. Therefore, taking both the noise level reduction and the peak level value in Fig. 9(c), the power spectral density dynamic is $18.4$dB higher in OFDM4 than in single-band interrogation, for a similar excitation level from the piezoelectric actuator.

 figure: Fig. 9.

Fig. 9. Interrogation of 1.3km SSMF excited by $440$Hz piezoelectric actuator at $420$m, 2 seconds measurement. Constant $B_{elec}=50$MHz

Download Full Size | PDF

In accordance with previous static simulations and measurements (in subsection 4.3 and subsection 5.1 respectively), enhanced sensing performances are noticed with the use of OFDM interrogation and subcarrier aggregation compared to single-carrier interrogation. The detection performance improvement is mainly driven by the false alarms reduction, which is in line with the observed reliability metric increase of the measurements.

6. Conclusion

We introduced the principle of digital carrier phase diversity to coded-$\Delta \phi$-OTDR, aiming at the mitigation of coherent fading effects. The Orthogonal Frequency Division Multiplexing (OFDM) interrogation scheme was proposed and tested in specific cases where we multiplexed the same codes, together with a subcarrier combination method for coherent fading mitigation.

Combined with Coherent-MIMO sensing which already solves the polarization fading effect, the digital OFDM interrogation scheme allows to mitigate coherent fading issues, thus providing a sensing system which is immune to strong fading. Experimental measurements confirmed the great efficiency of this approach in terms of sensitivity and reliability of the measurements. We reported an improvement of the strain sensitivity over $11$ dB SNR$_{\varepsilon,dB}$ using 4 subcarriers as compared to single-carrier interrogation, reaching a sensitivity noise floor level below $21\mathrm {p\varepsilon /\sqrt {Hz}}$ strain density on $1.3$km single mode fiber. With dynamic measurements, we demonstrated a gain up to $18.4\mathrm {dB}$ dynamic range when comparing a 4-subcarriers measurement with single-carrier one in case of a $440$Hz piezoelectric sine wave excitation.

Finally, in this paper, we showed that a digital multiplexing technique can be used to harness phase diversity and mitigate coherent fading. We believe that the proposed OFDM approach, applied for the first time to the sensing domain up to our knowledge, is a powerful tool to enhance the performance of future DAS systems.

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. H. Hartog, An Introduction to Distributed Optical Fibre Sensors (CRC Press, 2017).

2. A. Issatayeva, A. Beisenova, D. Tosi, and C. Molardi, “Fiber-optic based smart textiles for real-time monitoring of breathing rate,” Sensors 20(12), 3408 (2020). [CrossRef]  

3. A. H. Hartog, O. I. Kotov, and L. B. Liokumovich, “The optics of distributed vibration sensing,” Second. EAGE Work. on Perm. Reserv. Monit. p. 5 (2013).

4. G. A. Wellbrock, T. J. Xia, M.-F. Huang, Y. Chen, M. Salemi, Y.-K. Huang, P. Ji, E. Ip, and T. Wang, “First field trial of sensing vehicle speed, density, and road conditions by using fiber carrying high speed data,” in 2019 Optical Fiber Communications Conference and Exhibition (OFC) (2019), pp. 1–3.

5. M. R. Fernández-Ruiz, J. Pastor-Graells, H. F. Martins, A. Garcia-Ruiz, S. Martin-Lopez, and M. Gonzalez-Herraez, “Laser phase-noise cancellation in chirped-pulse distributed acoustic sensors,” J. Lightwave Technol. 36(4), 979–985 (2018).

6. Y. Shao, H. Liu, P. Peng, F. Pang, G. Yu, Z. Chen, N. Chen, and T. Wang, “Enhancing laser phase-noise immunity of distributed vibration sensor by phase-extraction φ-OTDR,” in 2018 Asia Communications and Photonics Conference (ACP) (2018), pp. 1–3.

7. J. C. Juarez and H. F. Taylor, “Polarization discrimination in a phase-sensitive optical time-domain reflectometer intrusion-sensor system,” Opt. Lett. 30(24), 3284 (2005). [CrossRef]  

8. J. J. Mompó, L. Shiloh, N. Arbel, N. Levanon, A. Loayssa, and A. Eyal, “Distributed dynamic strain sensing via perfect periodic coherent codes and a polarization diversity receiver,” J. Lightwave Technol. 37(18), 4597–4602 (2019). [CrossRef]  

9. C. Dorize and E. Awwad, “Enhancing the performance of coherent OTDR systems with polarization diversity complementary codes,” Opt. Express 26(10), 12878–12890 (2018). [CrossRef]  

10. S. Guerrier, C. Dorize, E. Awwad, and J. Renaudier, “Introducing coherent MIMO sensing, a fading-resilient, polarization-independent approach to ϕ-OTDR,” Opt. Express 28(14), 21081 (2020). [CrossRef]  

11. H. Gabai and A. Eyal, “On the sensitivity of distributed acoustic sensing,” Opt. Lett. 41(24), 5648 (2016). [CrossRef]  

12. P. Healey, “Fading in heterodyne OTDR,” Electron. Lett. 20(1), 30–32 (1984). [CrossRef]  

13. A. H. Hartog, “Rayleigh backscatter, issues and solutions,” in An Introduction to Distributed Optical Fibre Sensors, (CRC Press, 2017), chap. 6.5.1.

14. A. H. Hartog, L. B. Liokumovich, N. A. Ushakov, O. I. Kotov, T. Dean, T. Cuny, A. Constantinou, and F. V. Englich, “The use of multi-frequency acquisition to significantly improve the quality of fibre-optic-distributed vibration sensing,” Geophys. Prospect. 66(S1), 192–202 (2018). [CrossRef]  

15. J. Zhou, Z. Pan, Q. Ye, H. Cai, R. Qu, and Z. Fang, “Characteristics and explanations of interference fading of a $\phi$ -OTDR with a multi-frequency source,” J. Lightwave Technol. 31(17), 2947–2954 (2013). [CrossRef]  

16. M. D. Mermelstein, R. Posey, G. A. Johnson, and S. T. Vohra, “Rayleigh scattering optical frequency correlation in a single-mode optical fiber,” Opt. Lett. 26(2), 58 (2001). [CrossRef]  

17. M. Zabihi, Y. Chen, T. Zhou, J. Liu, Y. Shan, Z. Meng, F. Wang, Y. Zhang, X. Zhang, and M. Chen, “Continuous fading suppression method for ϕ-OTDR systems using optimum tracking over multiple probe frequencies,” J. Lightwave Technol. 37(14), 3602–3610 (2019).

18. D. Chen, Q. Liu, and Z. He, “Phase-detection distributed fiber-optic vibration sensor without fading-noise based on time-gated digital OFDR,” Opt. Express 25(7), 8315 (2017). [CrossRef]  

19. J. Pastor-Graells, H. Martins, A. Garcia-Ruiz, S. Martin-Lopez, and M. Gonzalez-Herraez, “Single-shot distributed temperature and strain tracking using direct detection phase-sensitive OTDR with chirped pulses,” Opt. Express 24(12), 13121–13133 (2016). [CrossRef]  

20. M. J. Murray and B. Redding, “Quantitative amplitude-measuring ϕ-OTDR with $p\varepsilon /\sqrt {Hz}$ sensitivity using a multi-frequency pulse train,” Opt. Lett. 45(18), 5226 (2020). [CrossRef]  

21. B. Redding, M. J. Murray, A. Davis, and C. Kirkendall, “Quantitative amplitude measuring φ-OTDR using multiple uncorrelated rayleigh backscattering realizations,” Opt. Express 27(24), 34952 (2019). [CrossRef]  

22. H. He, L. Yan, H. Qian, Y. Zhou, X. Zhang, B. Luo, W. Pan, X. Fan, and Z. He, “Suppression of the interference fading in phase-sensitive OTDR with phase-shift transform,” J. Lightwave Technol. 39(1), 295–302 (2021). [CrossRef]  

23. C. Dorize, S. Guerrier, E. Awwad, P. A. Nwakamma, H. Mardoyan, and J. Renaudier, “An OFDM-MIMO distributed acoustic sensing over deployed telecom fibers,” in 2021 Optical Fiber Communication Conference (OFC) (2021), p. 3.

24. A. Masoudi and T. P. Newson, “Analysis of distributed optical fibre acoustic sensors through numerical modelling,” Opt. Express 25(25), 32021 (2017). [CrossRef]  

25. Yuelan Lu, Zhu Tao, Chen Liang, and Bao Xiaoyi, “Distributed vibration sensor based on coherent detection of phase-OTDR,” J. Lightwave Technol. 28(22), 3243–3249 (2010). [CrossRef]  

26. C. Dorize, E. Awwad, S. Guerrier, and J. Renaudier, “Optimal probing sequences for polarization-multiplexed coherent phase OTDR,” in 2020 27th Optical Fiber Sensors Conference (OFS), (Optical Society of America, 2020 accepted, to be published).

27. J. Zhang, H. Wu, H. Zheng, J. Huang, G. Yin, T. Zhu, F. Qiu, X. Huang, D. Qu, and Y. Bai, “80 km fading free phase-sensitive reflectometry based on multi-carrier NLFM pulse without distributed amplification,” J. Lightwave Technol. 37(18), 4748–4754 (2019). [CrossRef]  

28. R. Goldman, A. Agmon, and M. Nazarathy, “Direct detection and coherent optical time-domain reflectometry with golay complementary codes,” J. Lightwave Technol. 31(13), 2207–2222 (2013).

29. M. J. Murray, A. Davis, C. Kirkendall, and B. Redding, “Speckle-based strain sensing in multimode fiber,” Opt. Express 27(20), 28494 (2019). [CrossRef]  

30. K. Shimizu, T. Horiguchi, and Y. Koyamada, “Characteristics and reduction of coherent fading noise in rayleigh backscattering measurement for optical fibers and components,” J. Lightwave Technol. 10(7), 982–987 (1992). [CrossRef]  

31. C. Dorize, S. Guerrier, E. Awwad, and J. Renaudier, “Identification of rayleigh fading induced phase artifacts in coherent differential ϕ-OTDR,” Opt. Lett. 46(11), 2754–2757 (2021). [CrossRef]  

32. SEAFOM, DAS Parameter Definitions and Tests, http://seafom.com, measuring sensor performance-02 ed. (2018).

33. E. Awwad, C. Dorize, S. Guerrier, and J. Renaudier, “Detection-localization-identification of vibrations over long distance SSMF with coherent delta-phi-OTDR,” J. Lightwave Technol. 38(12), 3089–3095 (2020). [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 (9)

Fig. 1.
Fig. 1. Backscattering model: $M$ segments with $N_{scat}$ scatterers within each segment. $\mathbf {H}_m$ is the $2\times 2$ Jones matrix describing the SOP evolution up to segment $m$.
Fig. 2.
Fig. 2. Pol.: polarization, MZM: Mach Zehnder Modulator, LO: Local Oscillator, DSP: Digital Signal Processing. $T_S$: symbol duration, $T_{code}$: probing code duration.
Fig. 3.
Fig. 3. OFDM combination for two subcarriers SC1 and SC2, one code duration, one polarization
Fig. 4.
Fig. 4. Simulation OFDM4, 1km SMF, $\Delta \nu = 75$Hz, AWGN at the receiver, individual subcarrier responses compared with combined response
Fig. 5.
Fig. 5. SNR probability density distributions. 1km static simulation, laser phase noise $\Delta \nu = 75$Hz, AWGN at the receiver. $\mathrm {SNR}_{\varepsilon }=0$dB corresponds to $1\mathrm {n}\epsilon$ strain sensitivity.
Fig. 6.
Fig. 6. SNR$_{\varepsilon,dB}$ on 1km, for 1 to 16 subcarriers, simulation with 5000 points per interrogation method, $\Delta \nu = 75$Hz and AWGN noise at the receiver
Fig. 7.
Fig. 7. OFDM-MIMO interrogation setup. MIMO: Multiple-Input/Multiple-Output, OFDM: Orthogonal Frequency Division Multiplexing, LO: Local Oscillator, ADC: Analog-to-Digital Converter, PZT: PiezoElectric Transducer
Fig. 8.
Fig. 8. Static experimental trace, OFDM interrogation of 1.3km SSMF, constant $B_{elec}=50$MHz configuration
Fig. 9.
Fig. 9. Interrogation of 1.3km SSMF excited by $440$Hz piezoelectric actuator at $420$m, 2 seconds measurement. Constant $B_{elec}=50$MHz

Equations (4)

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

E r x ( t ) = m = 0 M A m e j ϕ m e 2 α z m + j [ ω 0 ( t t m ) 2 β z m ] E 0 r e c t ( t t m T S )
A m e j ϕ m = i = 0 N s c a t a ( z i ) e j φ ( z i )
D t , d , c o m b i n e d = n = 1 N s c D t , d , n exp ( j ϕ ¯ d , n )
I n t e n s i t y d B = 10 log 10 ( | det ( H ) | ) S N R ε , d B = 20 log 10 ( 4 π n ξ G λ σ φ )
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.