Influence of stimulated Brillouin scattering on positioning accuracy of long-range dual Mach–Zehnder interferometric vibration sensors

The properties of noise induced by stimulated Brillouin scattering (SBS) in long-range interferometers and their influences on the positioning accuracy of dual Mach-Zehnder interferometric (DMZI) vibration sensing systems are studied. The SBS noise is found to be white and incoherent between the two arms of the interferometer in a 1-MHz bandwidth range. Experiments on 25-km long fibers show that the root mean square error (RMSE) of the positioning accuracy is consistent with the additive noise model for the time delay estimation theory. A low-pass filter can be properly designed to suppress the SBS noise and further achieve a maximum RMSE reduction of 6.7 dB. © The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI. [DOI: 10.1117/1.OE.55.11.116111]


Influence of stimulated Brillouin scattering on positioning accuracy of long-range dual Mach-Zehnder interferometric vibration sensors
Xiangge He, a Shangran Xie, b Shan Cao, a Fei Liu, a Xiaoping Zheng, a Min Zhang, c, * Han Yan, d and Guocai Chen d 1 Introduction Dual Mach-Zehnder interferometric (DMZI) vibration sensors can detect acoustic vibrations along the fiber, having the advantages of a long sensing distance, high sensitivity, and the ability of precision positioning. They have been widely applied in many areas, including oil pipeline monitoring, communication optical cable monitoring, and perimeter securities. [1][2][3][4] As a phase-modulated sensor, DMZI is very sensitive to external vibration events, while on the other hand it is also fairly easily affected by noise. To realize event positioning, the key function of DMZI, a cross-correlation algorithm is normally used to estimate the time delay between the bidirectional signals traveling along two interferometers caused by the same event. Previously, Quazi etc. developed the estimation model only considering additive noises. [5][6][7] By considering both additive and multiplicative noises, Xie et al. 8 have proposed the prediction theory on the positioning RMSE for DMZI vibration sensors.
In recent years, ultralong distance vibration sensing has attracted much interest to further extend the applications of DMZI. Tu et al. 9 reported a DMZI system with 160 m positioning precision on a 112-km-long sensing fiber. In this system, four single mode fibers (SMFs) in an optical cable were used. Two of them form the interference arms and others are lead fibers.
In long distance sensing applications, due to the presence of fiber loss, the power injected into the fiber needs to be high enough to get enough signal-to-noise ratio (SNR) on the photodetector, noting the minimum power arriving at the detector as P min . Given the fiber attenuation coefficient α and length L, the minimum power injected into the interference arm can be written as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 3 2 6 ; 4 5 1 Here the lengths of both sensing and lead fibers in real systems using one optical cable are considered. 9,10 On the other hand, a nonlinear effect (e.g., stimulated Brillouin scattering, SBS) starts to appear when its threshold P th is reached 11 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 3 2 6 ; 3 8 5 where K is the polarization factor, A eff is the effective mode area of the fiber, g B is the Brillouin peak gain coefficient, L eff ¼ ½1 − expð−αLÞ∕α is the effective length, and Δν p and Δν B are the linewidth of the pump light and Stokes light, respectively. For standard SMF and typical photodetectors, A eff ¼ 80 μm 2 , g B ¼ 5 × 10 −11 m∕W, α ¼ 0.21 dB∕km, and P min is typically 5 μW. Under such conditions, when P th is reached, the sensing fiber length L is ∼70 km. This means that to achieve longer distance fiber sensing, the noises contributed from a nonlinear effect should be considered. Since the interference signals in the DMZI system are bidirectional, the backscattering of a clockwise signal will bring noise into the counterclockwise signal and vice versa. Therefore, to achieve long-range vibration sensing, SBS noise could be a critical limiting factor, the contribution of which on the positioning accuracy of a DMZI system has so far not been studied.
In this paper, the properties of SBS noise in an interference system are experimentally measured. The influence of SBS noise on the positioning accuracy of the DMZI vibration sensor is further studied. It is found that SBS noise is white and incoherent in the 1-MHz bandwidth range and can be treated as an additive noise term in the positioning accuracy estimation model. We also show that a low-pass *Address all correspondence to: Min Zhang, E-mail: 1606871001@pku.edu.cn filter with a proper design can efficiently suppress SBS noise achieving a maximum positioning root mean square error (RMSE) reduction of 6.7 dB. Figure 1 shows the experimental setup we used to investigate the properties of SBS noise in interference signals. A narrow linewidth laser beam at 1550-nm wavelength (NKT Koheras BasiK E15) was injected into the interferometer after going through an isolator to suppress backreflection. Note that even though a broadband (or phase-scrambled) laser source can be used to significantly increase the SBS threshold, a narrow linewidth laser is used here to guarantee a high enough interference visibility and thus the SNR in the detection end, especially when the exact matching of the lengths of the two arms are difficult in practice. The MZI was formed by 50:50 fiber couplers C1 and C2 and two spools of 25km G652D SMFs (A and B). Photodetectors PD1 and PD2 were used to collect the backscattered and transmitted signals, respectively. Data were collected via a high-speed data acquisition card and then were further processed in an industrial personal computer (IPC). The IPC can also adjust the output power of the laser.

SBS Noise in Interferometric System
In the experiment, the power launched into fibers A and B was varied from 3.6 to 8.4 mW with a 0.4 mW step. First, just fiber A or B was connected and the transmitted and backscattered signals of it were measured (labeled as single-arm). Then both fibers A and B were connected and their interference signals were measured (labeled as dualarm). Since the frequency of the intrusion event is often lower than tens of kHz in practical applications, 9 the bandwidth we measured is 1 MHz. Figure 2 compares the typical time-domain and frequency-domain transmitted and backscattered signals for single-arm and dual-arm cases, at an input power of 6.8 mW. For single-arm, the energy of the transmitted signal The statistical properties of the temporal behavior of the Brillouin Stokes signal have been investigated by Gaeta and Boyd,13 in which large stochastic fluctuations were resolved for SMF. In our experiment, we show that the stochastic fluctuation not only exists in a single piece of fiber, but also for the dual-arm interferometers, as plotted in Fig. 2(g). In Gaeta's theory, the SBS medium behaves as a linear amplifier for the thermal noise that initiates the SBS process. Thermal noise can be assumed as a Gaussian random process, and the spectrum of the amplified Stokes field is also predicted to be Gaussian. As a result, the Stokes signals  of the two arms of the interferometer will not interfere with each other, explaining the observed stochastic fluctuation and the Gaussian distribution of the time domain interference signal [see Fig. 2(i)]. The SBS process, therefore, introduces white noise into the system without interference noise.
To further validate that the SBS noise does not introduce interference, the average intensities of the transmitted and backscattered signals were measured at different input powers, as plotted in Fig. 3. The red-triangles, blue-triangles, and black-asterisks represent the cases of single fibers A, B, and interference signals, respectively. The pink circles represent the sum of single fibers A and B (labeled as sum AB). As seen in the figures, SBS starts to occur when the input power reaches 4.8 mW. In the transmitted case, the average intensities of sum AB are different from that of the dual-arm interference signals, indicating the coherent property of the transmission. For backscattered signals, due to the lack of interference, sum AB are the exactly the same as that of the dual-arm configuration. Note that in this measurement, an average time of 0.5 s was applied to suppress the random noise while at the same time partially retaining the coherence information of the signal.

Influence of SBS Noise on the Positioning
Accuracy of DMZI Figure 4 sketches the experimental setup to study the positioning accuracy of DMZI with the presence of SBS noise. Couplers 2 and 3 along with fibers A and B form the DMZI structure. Both fibers A and B are 25 km long. A polarization controller was used to align the polarization state in fibers A and B. 9 Vibration events were generated by a piezoelectric ceramic transducer (PZT) near one end of fiber B. The PZT was driven by a 1-kHz sinusoidal wave generated by a signal generator. The rest of the system is the same as that illustrated in Fig. 1.
In this experiment, the input power was varied from 3.6 to 8.0 mW with 0.4 mW as the step. For each input power, 250 experiments were conducted to obtain the positioning RMSE. Crosscorrelation algorithm was applied to estimate the time delay b D with the unit of sampling point where Δt is the time delay in seconds and f s is the sampling frequency). The RMSE of b D can be expressed as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 3 2 6 ; 6 0 9 where D is the real value of the time delay, b D i is the estimated value for trial number i, and RMSE pt is the RMSE with the unit of sampling point. Equation (4) was used to convert from RMSE pt to RMSE in dB E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 3 2 6 ; 5 1 3 RMSEðdBÞ ¼ 10 log 10 ðRMSE pt Þ: (4) Figure 5 shows the experimental results. The blue-triangles and red-dots represent the SNR and RMSE (both in dB) versus the input power, respectively. Here the SNR is defined as SNR ¼ 10 × log 10 ðP signal ∕P noise Þ where P signal and P noise is the power of the 1-kHz sinusoidal signal (calculated by integrating the power spectral density from 500 to 1500 Hz) and the background noise, respectively. It can be clearly seen that above the SBS threshold (4.8-mW input power), SNR starts to decrease and RMSE starts to increase dramatically with input power, while below the SBS threshold both SNR and RMSE barely change with input power. The model developed by Quazi and Ianniello 5,6 can be used to predict the error of time delay estimation (TDE) when considering the additive noises. 8,10 Figure 6(a) plots the simulation results of σ CRLB and σ CPE when T ¼ 0.5 s and B ¼ 3 kHz (the bandwidth B is obtained by Welch power spectrum estimation method 14 in the experiment), where σ CRLB is the standard deviation of the Cramer-Rao lower bound (CRLB) for TDE, σ CPE is the standard deviation of the correlator performance estimate (CPE) for TDE, T is the observation window, and B is the bandwidth. CPE is more accurate for the crosscorrelation algorithm in the   presence of large estimation errors or anomalous estimates. 7 As can be seen in Fig. 6(a), when SNR is lower than about −6 dB, the predicted RMSE starts to increase rapidly, while when SNR is lower than −15 dB, the RMSE remains almost unchanged. Figure 6(b) shows the measured RMSE versus SNR in our experiment. It can be seen that from 0 to −23 dB, RMSE increases rapidly with the decrease of SNR, and when the SNR is lower than −23 dB the RMSE remains unchanged, which agrees well with the theoretical trend predicted in Fig. 6(a), which confirms the additive and incoherence properties of SBS noise. The experimentally measured RMSEs are larger than theoretical prediction. We attribute this difference to the residue polarization-induced fading effect caused by environmental noises accumulated along the entire bare optical fiber used in the experiment. 2 Finally, we apply a low pass digital filter to suppress the influence of SBS noise. A Chebyshev type I digital filter, with a cut-off frequency of 1.5 kHz, was inserted right before the crosscorrelation calculation in the algorithm, to filter out the noises of the bidirection signals. The applied filter can maintain the signal below 1 kHz and at the same time remove the most high frequency SBS noise. Figure 7 compares the filtered (blue) and unfiltered (red) results of positioning RMSE versus input power. It can be seen that a maximum reduction of about 6.7 dB can be achieved by using the lowpass filter. The filtered curve in Fig. 7 also shows that the positioning RMSE barely changes before and after SBS occurs (which is not the case for the unfiltered case), confirming the efficient suppression of SBS noise by the filter. A ∼6 dB of RMSE in the filtered case corresponds to ∼200-m positioning error, which is mainly limited by the residual environmental noises mentioned in Sec. 2.

Conclusion
The influences of SBS noise in DMZI vibration sensing systems are reported. Through detailed analysis of time-domain and frequency-domain properties, we show that SBS noise is a white noise and does not introduce interference between the two arms of the interferometer in a 1-MHz bandwidth range. Vibration positioning experiments on 25-km long bare SMF are conducted with the presence of SBS noise. The additive property of the SBS noise is supported by the agreement in trend of positioning RMSE versus SNR with the theoretical model. We also show that a maximum value of 6.7-dB RMSE reduction can be achieved by applying a low-pass filter with proper parameters. The results presented here can be further applied to any other fiber interferometric systems.
Xiangge He received his BE degree in optoelectronic science and engineering from Huazhong University of Science and Technology, Wuhan, China, in 2012. He is working toward his DE degree at Tsinghua University. His research interests include distributed optical fiber sensors.