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

Coherent Φ-OTDR based on I/Q demodulation and homodyne detection

Open Access Open Access

Abstract

We demonstrate a novel distributed acoustic sensing (DAS) system based on phase-sensitive optical time-domain reflectometry (Φ-OTDR). Both the phase and the amplitude of the Rayleigh scattering (RS) light can be demodulated in real-time. The technique is based on I/Q demodulation and homodyne detection using a 90° optical hybrid. The theoretical analysis is given, and as a proof of the concept, the dynamic strain sensing is experimentally demonstrated, with a sensing range of 12.566 km and a spatial resolution of 10 m.

© 2016 Optical Society of America

1. Introduction

The distributed optical fiber sensing (DOFS) systems have been widely used in industries involved with security surveillance and structural health monitoring, owing to its large-scale monitoring range, low cost per monitored point, simple installation and geometric versatility [1]. It has been extensively studied in terms of the monitoring range, sensing parameter, spatial resolution and so on [2–7 ].

As one of the typical DOFS systems, phase sensitive optical time-domain reflectometer (Φ-OTDR) has been demonstrated as a promising technique, because of its ability to monitoring dynamic process and its high sensitivity [8]. Recently, distributed acoustic sensing (DAS) technology based on Φ-OTDR has attracted much attention due to its ability to measure the phase variation of the Rayleigh scattering (RS) signal and its promising industrial applications [9–12 ]. The reported demodulation schemes for DAS include 3 × 3 coupler demodulation [9,10 ], phase generated carrier algorithm [11], and digital coherent detection [12]. In [9–11 ], an interferometer structure is needed to retrieve the phase variation along the fiber. The interferometer structure is sensitive to the environment hence the performances of the systems may be influenced. In [12], the researchers use digital coherent detection. There is no interferometer structure in this scheme, and the phase of the signal is demodulated using software calculation, however, it requires high-speed analog to digital converter (ADC) and the large amount of data makes the data analysis quite challenging.

In the paper, an efficient phase demodulation scheme is demonstrated for the Φ-OTDR system using I/Q demodulation and homodyne detection. With the 90° hybrid, both the in-phase (I) and the quadrature (Q) components of the backscattering light field can be detected, then the phase information of the light signal can be demodulated without an interferometer structure. Moreover, the electrical signal is in baseband since the frequencies of the signal and the local-oscillator (LO) are the same, therefore only low-speed ADC is sufficient. Dynamic strain sensing is experimentally demonstrated with a sensing range of ~12.6 km and a spatial resolution of 10 m, as a proof of concept.

2. Principle

Figure 1 shows the general block diagram of I/Q demodulation used in our scheme. In the DAS system based on Φ-OTDR, the detected signal is the coherent Rayleigh scattering (RS) of the probe pulse. Assuming that the RS light can be described as

Asexp[iωst+iφ(t)+iφs]
where As, ωs,φsis the scalar amplitude, the angular frequency and initial phase of the signal, respectively. φ(t) is the phase change induced by the refraction index variation along the fiber. The signal is split into two beams by the hybrid, and the accurate 90° phase shift is induced in one beam of the signal [13] and the corresponding light field changes into

 figure: Fig. 1

Fig. 1 A general block diagram of I/Q demodulation using 90° hybrid.

Download Full Size | PDF

Asexp[iωst+iφ(t)+iφs+i(π/2)].

Then the two signal beams are respectively superposed with the LO. The LO can be described as

Aoexp(iωot+iφo)
where Ao, ωo, φo is the scalar amplitude, the angular frequency and initial phase of LO, respectively. Each of the two outputs of the hybrid is then detected by a bandwidth-limited photodetector, and the photo-currents should be [14]:
{I(t)=12R[AsAocos((ωsωo)t+φ(t)+φsφo)+A+s2A]o2Q(t)=12R[AsAosin((ωsωo)t+φ(t)+φsφo)+A+s2A]o2
In the equation, the coefficient 1/2 comes from the splitting of the signal and R is the responsivity of the two detectors (providing the responsivities of the two detectors are the same), I/Q mean the in-phase and the quadrature components of the electrical signal, respectively. The AC components of the photo-currents are:
{IAC12RAsAocos((ωsωo)t+φ(t)+φsφo)QAC12RAsAosin((ωsωo)t+φ(t)+φsφo)
In our experiment, the signal and the LO come from the same laser, and no frequency different is present, so ωsωo=0,φsφo=φd. φd is the difference between the initial phases of the LO and the signal, and we assume that it is a constant. Then we have:
{IAC12RAsAocos(φ(t)+φd)QAC12RAsAosin(φ(t)+φd)
Then the amplitude and the phase of the RS light can be obtained:
{AsIAC2+QAC2φ(t)arctan(QACIAC)-φd
Here, the value of the arctan function is from -π/2 toπ/2. By inducing the phase unwrapping algorithm [15], the range of the demodulated phase becomes from negative infinite to positive infinite. Then the phase and amplitude of the RS along the fiber can be demodulated. By subtracting phases of two certain positions, the phase change of the corresponding part of the fiber can be demodulated as
Δφ(t)=φZ1(t)φZ2(t)
which is linearly associated with the strain induced by the acoustic field [16].

3. Experimental setup

The experimental scheme is depicted in Fig. 2 . An ultra-narrow-linewidth (100 Hz) laser operating at 1550.1 nm is used as the light source with polarization-maintaining (PM) output. The laser output is split into two branches by a 1:99 PM coupler. The 99% branch is modulated by an acoustic-optic modulator (AOM) with 80 MHz frequency shift to generate the pulsed probe wave (with 100 ns pulse width, allowing for 10 m spatial resolution). The repetition rate of the probe pulse is 5 kHz. The probe pulse is injected into the sensing fiber through a circulator and the RS signal is injected into the 90° hybrid. The 1% branch of the laser output is used as the LO. An acoustic frequency shifter (AFS) is used in order to introduce 80 MHz frequency shift into LO, to make sure the signal and the LO are with the same frequency. Before the LO is injected into the hybrid, a polarization controller (PC) is inserted in order to match the LO with the selected polarization branch of the hybrid. The two outputs of the hybrid are converted into electrical signals by a two ports AC coupled detector and then sampled by a 25 MS/s ADC. The process of phase unwrapping algorithm and signal reconstruction are completed in real-time with a personal computer. The sensing fiber is 12.56 km long. As represented in Fig. 2, at the position of ~12.4 km, 10 m bare fiber is coiled over a cylindrical PZT, used as the test point. A two-channel arbitrary waveform generator (AWG) is used in our experiment. One channel of the AWG is used to provide RF signal to the AOM driver and the trigger of the ADC, and the other one is used to drive the PZT.

 figure: Fig. 2

Fig. 2 Experimental setup (Cir: circulator; SMF: single mode fiber).

Download Full Size | PDF

It should be noted that a dual polarizations 90° hybrid is used in the system. In optical communication systems, the coherent receiver with a dual polarization 90° hybrid has been proved to be an effective method to measure the amplitude and phase of optical signals with arbitrary polarization [17]. Therefore, essentially this scheme is able to cope with the polarization fading phenomenon in OTDR systems, if all the outputs of the hybrid are used. The polarization dependent attributes of the approach are the topic of a further study. In this work, only one polarization channel is used since we only focus on the demonstration of the phase recovery.

4. Results and discussions

From Eq. (8) one can see that the interval between Z1 and Z2 can actually limit the spatial resolution. Typically, the interval is set to a length no longer than the spatial resolution limited by the pulse width of the probe light. In our experiment, the pulse width is 100 ns and it allows the highest spatial resolution of 10 m. Besides, the sampling rate of data acquisition card is 25 MS/s, which means the spatial interval between every two adjacent data points is 4 m. To match the 10 m spatial resolution we deliberately choose an interval of 8 m, using the #3100 data point and the #3102 data point to demodulate the phase difference.

Figure 3 shows the trace of the RS light along the fiber without PZT perturbation. The trace shows the speckle pattern as in conventional Φ-OTDR. To test the ability of strain change detection of the system, sinusoidal signals are applied on the PZT with different voltage magnitudes or different frequencies. By applying the phase demodulation approach proposed in Section 2, the phase difference at the relevant part of the fiber can be demodulated. The results are filtered by a low pass filter with a bandwidth of 2.5 kHz.

 figure: Fig. 3

Fig. 3 Rayleigh backscattering trace of the DAS.

Download Full Size | PDF

First, the sinusoidal signal that we apply on the PZT has a magnitude of 2 V and a frequency of 500 Hz. The results are showed in Fig. 4 . The dotted blue line in Fig. 4(a) represents the experimental data, while the red solid line represents the sinusoidal fitting curve. The demodulation phase differences present a clear sinusoidal variation over time. The amplitude of the fitting curve is about 2.51 rad. The pink dash line in Fig. 4(b) shows the FFT of the demodulated signals. Clear peak can be seen at 500 Hz and the calculated SNR is about 34.1 dB (SNR is calculated as 20*log(As/A n),where As is the amplitude of the signal, An is the root mean square of the background noise [11]).

 figure: Fig. 4

Fig. 4 (a) The time-domain trace and (b) the spectrum of the demodulated phase variation.

Download Full Size | PDF

Then we change the magnitude of the sinusoidal signal applied on the PZT, in order to test the relationship between the magnitude of the applied voltage on the PZT and the demodulated phase change (at 500 Hz). The magnitude of the applied voltage is set from 1 V to 2.4 V with a step of 0.2 V. The experimental results are shown in Fig. 5 . The blue dots are the amplitudes of the sinusoidal fitting curve of the phase difference traces and the red solid line is the linear fitting line of the amplitudes. The R-square of this fitting is 0.974, which shows a reasonable linear response of the system.

 figure: Fig. 5

Fig. 5 Amplitude of the phase variation v.s. voltage applied to the PZT.

Download Full Size | PDF

To further demonstrate the ability of the system to sense the dynamic events, we change the frequency of the sinusoidal signal applied on the PZT as 400 Hz and 600 Hz, respectively, and set the magnitude to 2 V. As it can be seen clearly in Fig. 6 , similar peaks stand at 400Hz and 600 Hz, respectively. The blue solid line represents experimental result with frequency of 400 Hz, and green dotted line represents experimental result with frequency of 600 Hz. Accurate frequency responses of the system is verified. Moreover, the SNRs in frequency domain are 42.0 dB, 34.1 dB and 34.8 dB corresponding to signals of 400 Hz, 500 Hz and 600 Hz, respectively. The differences between the SNRs are considered to be mainly induced by the non-uniform frequency response of the PZT.

 figure: Fig. 6

Fig. 6 Spectrums with different perturbation frequencies.

Download Full Size | PDF

Figure 4-6 illustrate the capability of the system of detecting the acoustic event by measuring the phase and amplitude of the RS signals in real-time. The amplitude and the frequency of the applied signal can be measured correctly and the results fit the applied signals well.

There are several advancements in this scheme. First, the data processing load is reduced. From the principle shown in Section 2 one can see that, the I/Q demodulation is mainly realized in optical domain rather than in electrical domain. Consequently, calculating amount in the scheme can be relatively small compared with the digital coherent detection in [12]. In addition, compared with the heterodyne coherent detection [11, 12 ], the signals in the system with our scheme are located in baseband, which allows the use of an ADC with relatively low sampling rate, hence the data amount is significantly reduced. Second, the mature telecom technique can assure the hybrid to induce accurate 90° phase shift and it is much less sensitive to environmental vibration and temperature drift compared with the interferometer structure. Moreover, the demodulation proposed in this paper also benefits from the amplification from the LO. The LO can improve the sensitivity of the system, as a result, no pre-amplifier is needed in the experimental system and a sensing range of more than 10 km is achieved.

6. Conclusion

In this paper, we propose a DAS system by introducing an efficient phase demodulation scheme in Φ-OTDR. Both the amplitude and the phase of the RS light are recovered by I/Q demodulation and homodyne detection utilizing a 90° hybrid. The dynamic phase difference induced by the PZT is detected. A DAS system with a sensing range of 12.56 km and a spatial resolution of 10 m is demonstrated.

Acknowledgments

This work is supported by Natural Science Foundation of China (61205048, 61290312, 41527805), Research Fund for the Doctoral Program of Higher Education of China (20120185120003), the PCSIRT project (IRT1218), and the 111 project (B14039).

References and links

1. L. Thévenaz, “Next generation of optical fibre sensors: new concepts and perspectives,” Proc. SPIE 9157, 23rd International Conference on Optical Fibre Sensors, 9157AN (2014).

2. Z. N. Wang, J. J. Zeng, J. Li, M. Q. Fan, H. Wu, F. Peng, L. Zhang, Y. Zhou, and Y. J. Rao, “Ultra-long phase-sensitive OTDR with hybrid distributed amplification,” Opt. Lett. 39(20), 5866–5869 (2014). [CrossRef]   [PubMed]  

3. S. Martin-Lopez, M. Alcon-Camas, F. Rodriguez, P. Corredera, J. D. Ania-Castañon, L. Thévenaz, and M. Gonzalez-Herraez, “Brillouin optical time-domain analysis assisted by second-order Raman amplification,” Opt. Express 18(18), 18769–18778 (2010). [CrossRef]   [PubMed]  

4. L. Zhang, Z. N. Wang, J. Li, J. J. Zeng, X. H. Jia, and Y. J. Rao, “Ultra-long dual-sideband BOTDA with balanced detection,” Opt. Laser Technol. 68, 206–210 (2015). [CrossRef]  

5. Y. J. Rao, J. Z. Li, Z. L. Ran, and K. L. Xie, “Distributed intrusion detection based on combination of φ-OTDR and POTDR,” Proc. SPIE 7004, 19th International Conference on Optical Fibre Sensors, 700461 (2008). [CrossRef]  

6. W. Li, X. Bao, Y. Li, and L. Chen, “Differential pulse-width pair BOTDA for high spatial resolution sensing,” Opt. Express 16(26), 21616–21625 (2008). [CrossRef]   [PubMed]  

7. H. Martins, S. Martin-Lopez, P. Corredera, J. D. Ania-Castanon, O. Frazao, and M. Gonzalez-Herraez, “Distributed vibration sensing over 125km with enhanced SNR using phi-OTDR over an URFL cavity,” J. Lightwave Technol. 33(12), 2628–2632 (2015). [CrossRef]  

8. H. F. Tarlor and C. E. Lee, “Apparatus and method for fiber optic intrusion sensing,” U.S. Patent 5194847 (1993).

9. A. Masoudi, M. Belaland, and T. P. Newson, “A distributed optical fibre dynamic strain sensor based on phase-OTDR,” Meas. Sci. Technol. 24(8), 085204 (2013). [CrossRef]  

10. C. Wang, C. Wang, Y. Shang, S. Ying, X. H. Liu, and G. D. Peng, “Distributed acoustic mapping based on interferometry of phase optical time-domain reflectometry,” Opt. Commun. 346, 172–177 (2015). [CrossRef]  

11. G. S. Fang, T. W. Xu, S. W. Feng, and F. Li, “Phase-sensitive optical time domain reflectometer based on phase-generated carrier algorithm,” J. Lightwave Technol. 33(13), 2811–2816 (2015). [CrossRef]  

12. Z. Pan, K. Liang, Q. Ye, H. Cai, R. Qu, and Z. Fang, “Phase-sensitive OTDR system based on digital coherent detection,” in Optical Sensors and Biophotonics, J. Popp, D. Matthews, J. Tian, and C. Yang, eds., Proc. SPIE 8311 (Optical Society of America, 2011), paper 83110S.

13. M. G. Taylor, “Phase estimation methods for optical coherent detection using digital signal processing,” J. Lightwave Technol. 27(7), 901–914 (2009). [CrossRef]  

14. J. Zhou, Z. Q. Pan, Q. Ye, H. W. Cai, R. H. Qu, and Z. J. Fang, “Phase demodulation technology using a multi- frequency source discrimination of interference-fading induced false alarm in a Φ-OTDR system,” Chin. J. Lasers 40(9), 0905003 (2013). [CrossRef]  

15. R. Hui and M. S. O’Sullivan, Fiber Optic Measurement Techniques (Academic Press, 2009).

16. G. J. Tu, X. P. Zhang, Y. X. Zhang, F. Zhu, L. Xia, and B. Nakarmi, “The Development of an Φ-OTDR System for Quantitative Vibration Measurement,” IEEE Photonics Technol. Lett. 27, 1349–1352 (2015). [CrossRef]  

17. B. Szafraniec, T. S. Marshall, and B. Nebendahl, “Performance monitoring and measurement techniques for coherent optical systems,” J. Lightwave Technol. 31(4), 648–663 (2013). [CrossRef]  

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

Fig. 1
Fig. 1 A general block diagram of I/Q demodulation using 90° hybrid.
Fig. 2
Fig. 2 Experimental setup (Cir: circulator; SMF: single mode fiber).
Fig. 3
Fig. 3 Rayleigh backscattering trace of the DAS.
Fig. 4
Fig. 4 (a) The time-domain trace and (b) the spectrum of the demodulated phase variation.
Fig. 5
Fig. 5 Amplitude of the phase variation v.s. voltage applied to the PZT.
Fig. 6
Fig. 6 Spectrums with different perturbation frequencies.

Equations (8)

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

A s exp [ i ω s t + i φ ( t ) + i φ s ]
A s exp [ i ω s t + i φ ( t ) + i φ s + i ( π / 2 ) ] .
A o exp ( i ω o t + i φ o )
{ I ( t ) = 1 2 R [ A s A o cos ( ( ω s ω o ) t + φ ( t ) + φ s φ o ) + A + s 2 A ] o 2 Q ( t ) = 1 2 R [ A s A o sin ( ( ω s ω o ) t + φ ( t ) + φ s φ o ) + A + s 2 A ] o 2
{ I A C 1 2 R A s A o cos ( ( ω s ω o ) t + φ ( t ) + φ s φ o ) Q A C 1 2 R A s A o sin ( ( ω s ω o ) t + φ ( t ) + φ s φ o )
{ I A C 1 2 R A s A o cos ( φ ( t ) + φ d ) Q A C 1 2 R A s A o sin ( φ ( t ) + φ d )
{ A s I A C 2 + Q A C 2 φ ( t ) arc tan ( Q A C I A C ) - φ d
Δ φ ( t ) = φ Z 1 ( t ) φ Z 2 ( t )
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.