Synthetic Spectrum Approach for Brillouin Optical Time-Domain Reflectometry
Abstract
: We propose a novel method to improve the spatial resolution of Brillouin optical time-domain reflectometry (BOTDR), referred to as synthetic BOTDR (S-BOTDR), and experimentally verify the resolution improvements. Due to the uncertainty relation between position and frequency, the spatial resolution of a conventional BOTDR system has been limited to about one meter. In S-BOTDR, a synthetic spectrum is obtained by combining four Brillouin spectrums measured with different composite pump lights and different composite low-pass filters. We mathematically show that the resolution limit, in principle, for conventional BOTDR can be surpassed by S-BOTDR and experimentally prove that S-BOTDR attained a 10-cm spatial resolution. To the best of our knowledge, this has never been achieved or reported.1. Introduction
The spectrum of Brillouin scattering in optical fiber shifts in proportion to changes in the strain and temperature of the fiber. Brillouin optical time-domain analysis (BOTDA) [1,2] and Brillouin optical time-domain reflectometry (BOTDR) [3,4] are distributed measurement techniques utilizing this property by injecting a pulsed pump light and observing the scattered light. The amount of strain/temperature change is estimated by measuring the spectral shift, while its position is determined by the round-trip time of the light. Since BOTDR uses only one end of a fiber, it is suitable for long distance measurement, whereas BOTDA uses pump and probe lights that are injected from both ends of a fiber.
To improve the spatial resolution of both techniques, it is necessary to narrow the pulse; however, the linewidth of the observed spectrum then becomes wider and makes measurement of the spectral shift difficult. For this reason, the spatial resolution has been limited to about one meter in both techniques [5,6].
For BOTDA, Bao et al. [7] found experimentally a phenomenon in which the spectral linewidth becomes shorter when a very short pulse of about 1 ns is used. This phenomenon arises only when there is light leakage. Inspired by this discovery, various resolution improvement methods have been proposed [8–14]. The pump light of BOTDA plays two roles: phonon excitation and scattering by the phonons. The width of the pump light must be longer than the phonon lifetime (about 9 ns) for phonon excitation, whereas it must be much shorter than 10 ns, which corresponds to a one-meter resolution, for high spatial resolution. To satisfy these incompatible requirements, an idea to construct a pump light using a long and a short element was born. The long element, which is a long pulse or a CWwave, takes the role of phonon excitation, and the short element takes the role of being scattered by the phonons. A part of the spectrum is by the combination of these elements. By emphasizing this desired part, high-resolution measurement could be attained.
For the construction of pump light from two elements, the authors of [8–11] used amplitude modulation and [12,14] used phase modulation. In any construction method, the measured spectrum includes both desired and undesired parts. Parameters must be optimized to emphasize only the desired part. Instead of that, the authors of [12–14] proposed methods to suppress or cancel the undesired part of the spectrum by combining two different measurements. By these methods, high-resolution measurement of centimeter-order for BOTDA has been attained.
On the other hand, the pump light of BOTDR plays a single role, since the BOTDR utilizes spontaneous scattering, so that it is difficult to attain high resolution by the same idea as for BOTDA. For resolution improvement of BOTDR, Koyamada et al. [15] proposed a double pulse method that yields an oscillatory Brillouin spectrum and verified 20-cm spatial resolution by experiment.
In this paper, we propose a novel method, referred to as synthetic BOTDR (S-BOTDR), to improve the spatial resolution of BOTDR. In this method, a synthetic Brillouin spectrum is constructed by combining several spectrums obtained by BOTDR measurements with different pump lights and low-pass filters. The pump lights and low-pass filters are composed of short and long elements with phase differences. We mathematically show that the resolution limit in principle for conventional BOTDR could be overcome by S-BOTDR and experimentally prove that 10-cm spatial resolution, which is much smaller than the conventional BOTDR limit, is attained by S-BOTDR.
The remainder of this paper is organized as follows: Section 2 presents the mathematical formulation of BOTDR and its performance limit in principle. In Section 3, the proposed method, called S-BOTDR, is described. Sections 4 and 5 are devoted to the evaluation of the proposed method through simulations and experiments, respectively. Finally, we conclude our paper in Section 6.
2. Mathematical Model of BOTDR
2.1. Sensing Mechanism of BOTDR
A BOTDR system is shown in Figure 1. The light from the light source is divided into a pump pulse and a reference continuous wave. In an optical fiber, the acoustic waves are always excited by thermal fluctuation and the pump pulse is scattered by the acoustic waves in all parts of the fiber. Since the acoustic waves move with the sound velocity of the medium, the backscattered light has a Doppler shift, which is called a Brillouin frequency shift (BFS). As the sound velocity changes in proportion to the change in strain or temperature in each position of the fiber, the BFS also changes in proportion to them. In BOTDR, the backscattered light is heterodyned with the reference wave and the spectrum of the interfering light; i.e., the Brillouin spectrum is obtained. The BFS is estimated from the Brillouin spectrum as its center frequency. Since the round-trip time of the light corresponds to the position in the fiber, the BFS at each position in the fiber is obtained. Thus, distributed sensing of strain and temperature is possible using the Brillouin scattered light.
2.2. Basic Equations of BOTDR
Brillouin scattering in optical fiber is described by the following equations [16–18]:
Random force R(z, t) is a circular symmetric complex white noise that is white both in space and time; i.e., it is characterized by E [R(z, t)R*(z′, t′)] = Qδ(z−z′)δ(t−t′), where E[·] stands for expectation and Q is a constant. The two terms in the right-hand side (RHS) of Equation (3) correspond to stimulated and spontaneous Brillouin scattering, respectively Under usual BOTDR conditions, as analyzed in a previous paper [19], the largest part of the measured spectrum arises from spontaneous scattering, and the stimulated scattering term can be neglected. Moreover, then, since it becomes |Es| ≪ |Ep|, the RHS of Equation (1) can be neglected. Therefore, BOTDR is described by the following equations:
2.3. Analytical Solution to the BOTDR Equations
The solution to last section's BOTDR equations can be represented analytically. For simplicity, assuming that the fiber loss is small, we set α = 0.
First, the solution to Equation (1) under boundary Condition Equation (7) is represented as:
The backscattered light returned to the input end of an optical fiber in BOTDR is represented as:
2.4. Brillouin Spectrum and Point Spread Function
The Brillouin spectrum obtained by one measurement of BOTDR is represented as:
Since Y(t, ν) is a ccG process, V(t, ν) becomes a random variable with an exponential distribution for each t and ν, so its variance equals the square of its expectation:
We can calculate the expectation of V(t, ν) by using the statistical property of Equation (11) as:
We note that L(t, ν) is a time-varying Lorentzian spectrum and is determined only by the characteristics of an optical fiber, whereas ψ(t, ν) depends on the sensing mechanism of BOTDR. Since Equation (18) implies that function, ψ(t, ν), obscures the details of L(t, ν), we refer to Ψ(t, ν) as a point spread function (PSF).
2.5. Performance Limit of BOTDR
Though the ideal Brillouin spectrum is the Lorentzian spectrum, the observed Brillouin spectrum is spread by the PSF. Therefore, it is desirable for the PSF to have a shape close to a two-dimensional δ-function; that is, narrow in both time and frequency directions. However, by the following uncertainty relation (see Appendix B):
3. BOTDR by the Synthetic Approach
To overcome the resolution limitation of conventional BOTDR described in the previous section, a synthetic spectrum approach, which we refer to as synthetic BOTDR (S-BOTDR), is proposed.
3.1. Elements of the Point Spread Function
We consider a short pulse, f1(t), and a long pulse, f2(t), as elements of composite pump light:
Similar to the PPP-BOTDA method [8,9], D1 is set as less than the spatial resolution to be attained and D2 is taken as sufficiently longer than the acoustic wave lifetime, 2/ΓB, say about 9 ns, to improve frequency resolution. A low-pass filter is also composed of two elements, each of which is a matched filter of short or long pulses:
By using the above elements, a pump light and a low-pass filter are composed as:
Therefore, among the terms of the RHS of Equation (28), only the term, , has the ideal property that it is localized both in time and frequency, and the other terms are undesired elements for the BOTDR measurement.
3.2. PSF by the Synthetic Approach
To extract only element from the PSF, we prepare p phase pairs ((θj, ϕj), j = 1, 2, ⋯, p) and p pairs of pump lights and low-pass filters:
This problem is equivalent to finding θj, ϕj, cj, j = 1, 2, ⋯, p that satisfy:
The PSF by the synthetic approach is represented as:
For two pulses of Equations (22) and (23), the synthetic PSF becomes:
The PSFs for a conventional BOTDR and S-BOTDR are plotted in Figure 4. While the PSF of the conventional BOTDR spreads in either the time or frequency direction, that of S-BOTDR is concentrated at the origin and approximates a δ-function. To clarify the mechanism of generating such a δ-like function, the composition of the synthetic PSF is shown in Figure 5. Although each of the PSF, Ψ(j), j = 1, 2, 3, 4, concentrates in the time direction, it is spread in the frequency direction. This spreading is partially suppressed with two-by-two combination, Ψ(2) − Ψ(1) and Ψ(4) − Ψ(3), and is completely suppressed with the use of all PFS elements.
3.3. S-BOTDR System
As shown in the analysis of the previous section, the synthetic Brillouin spectrum has the ideal property of being close to a Lorentzian spectrum. An S-BOTDR system for obtaining the synthetic Brillouin spectrum is constructed as shown in Figure 6. The composition of the synthetic Brillouin spectrum is followed by Equation (32) and is illustrated in Figure 7. Although the spectrum of each measurement exhibits high spatial resolution, it is considerably different from the Lorentzian spectrum. Only the combination of all four spectrums produces the Lorentzian spectrum.
The ideal property described in the previous section is for the expectation of the synthetic Brillouin spectrum. Achieving this property requires averaging or summing of the spectrums obtained by a number of measurements.
3.4. Parameter Selection for S-BOTDR
The pump lights of S-BOTDR are composed of short and long pulses and include parameters t0 and r: t0 is the difference of pulse start times and r is the amplitude ratio of the two pulses. These parameters affect the performance of S-BOTDR.
Selection of t0 If the short pulse is included in the long pulse, i.e., if 0 ≤ t0 ≤ D2 − D1, the integral of Equation (40) takes the maximum value. In addition, the value of t0 affects the concentration of the PSF near the origin. The degree of concentration is estimated by the integral:
The values of C(ν; t0) are plotted for different t0's in Figure 8. Here, for t0 = (D2 − D1)/2, the short pulse is located at the center of the long pulse; for t0 = 0, the short pulse is at the edge of the long pulse, and the medium case is t0 = (D2 − D1)/4. From Figure 8, we find that when the short pulse is located at the edge of the long pulse, the degree of concentration is reduced. Thus, we adopt t0 = (D2 − D1)/2.
Selection of r It is desirable for S-BOTDR to suppress the synthetic spectrum's fluctuation as much as possible while keeping the peak level high. To meet these requirements, we define the signal-to-fluctuation ratio (SFR) as:
The SFR of S-BOTDR is obtained by substituting (46) into (76) in Appendix C. The SFR values for various D1 and D2 are plotted in Figure 9. From this figure, we find that SFR is nearly proportional to D1 for fixed D2. On the other hand, BFS estimation error variance is inversely proportional to SFR [20]. Therefore, the BFS estimation error variance is nearly proportional to 1/D1.
4. Evaluation by Simulation
Numerical simulations were performed to verify the S-BOTDR and compare it with the conventional BOTDR. The condition of the optical fiber is shown in Figure 10. In the figure, an assumed BFS is shown in (a) and the corresponding Lorentzian spectrum is shown in (b). Simulation results are shown in Figures 11, 12, and 13. The parameter values for processing are as follows:
Conventional BOTDR (Figures 11 and 12)
Pulse width is D = 10 and 1 ns.
S-BOTDR (Figure 13)
Pulse widths are D1 = 1 ns for a short pulse and D2 = 50 ns for a long pulse. The amplitude ratio is r = 0.066. The short pulse is located at the center of the long pulse.
All results were obtained by summing 212 trials. For conventional BOTDR with a 10-ns pulse in Figure 11, the linewidth of the spectrum was small and smooth; however, spatial resolution was inferior and could not detect any interval less than 50 cm. For the conventional BOTDR with a 1-ns pulse in Figure 12, because the linewidth of the spectrum was quite large (about 1 GHz), the estimation accuracy of the BFS was inferior. In the case of the S-BOTDR in Figure 13, the Brillouin spectrum was close to the Lorentzian spectrum; hence, accurate BFS estimates were obtained. Thus, S-BOTDR has a high-resolution capability exceeding that of conventional BOTDR.
The simulation was performed using the exact values of the input parameters. In an actual situation, these parameters are, however, subject to errors. As in S-BOTDR four measurements are combined, the input light amplitude fluctuation and phase-setting errors become issues of concern. Among them, the light amplitude fluctuation becomes negligible by averaging many replicate measurements. The errors in the phase differences between short and long elements are, on the other hand, anticipated. We performed a series of additional simulations in which the phase differences were subject to random errors. The results clearly demonstrated that they did not affect the spatial resolution, while the BFS estimation error increased only slightly, even if phase errors as high as 40 degrees were applied.
5. Experimental Verification
We performed BOTDR and S-BOTDR experiments using a test fiber. We constructed a test fiber by connecting two kinds of fibers: A and B (Figure 14a). The BFS between fibers A and B differs, and the difference corresponds to a strain of 1,200 μ∊. In the experiments, fibers A and B are considered strained fibers and strain-free fibers, respectively. The replicate measurement numbers were 216 in all cases. The experimental result for BOTDR is shown in Figure 15. The pulse width of the pump light for BOTDR was 5 ns. Strain estimates and the Brillouin spectrum at the 5-m strain interval are plotted in Figure 15a,b, respectively. The linewidth of the Brillouin spectrum exceeded 360 MHz (Figure 15b). Although intervals larger than or equal to one meter were accurately detected, the strain estimates of the other intervals differed from the true values (Figure 15a).
We performed S-BOTDR experiments using the following combination of short and long elements for the pump light and the low-pass filter:
Case 1 D1 = 7.8 ns, D2 = 60 ns, r = 0.26
Case 2 D1 = 2.7 ns, D2 = 50 ns, r = 0.14
Case 3 D1 = 1.6 ns, D2 = 50 ns, r = 0.093
6. Conclusions
We showed that conventional BOTDR has, in principle, a resolution limit and proposed a synthetic spectrum approach, referred to as S-BOTDR, to overcome this limit. S-BOTDR uses composite pump lights and low-pass filters that are composed of short and long elements with different phases. Four BOTDR measurements are performed using four specific pairs of phases. By combining four Brillouin spectrums with specific weights, an ideal spectrum is obtained. The spectrum is close to the Lorentzian spectrum, and the resolution limit, in principle, disappears. We experimentally evaluated S-BOTDR and proved that we attained a 10-cm resolution that exceeds the resolution limit of conventional BOTDR.
Acknowledgments
A portion of the results presented in this paper was obtained during a joint project with Mitsubishi Corporation, with financial support of Japan Oil, Gas and Metals National Corporation (JOGMEC). This support is gratefully acknowledged.
Appendix A. Derivation of (18)
In this appendix, we derive Equation (18), which is the convolutional representation of the expectation of the Brillouin spectrum. Substituting Equation (15) into Equation (16), we have:
Appendix B. Uncertainty Relation of PSF
The frequency integral of the PSF Equation (20) gives a temporal PSF as follows:
Appendix C. Signal-to-Fluctuation Ratio of S-BOTDR
The synthetic spectrum of S-BOTDR Vs(t, ν) is the weighted sum of elements V(j)(t,ν), j = 1, 2, 3, 4, as shown in Equation (32). Since all elements are statistically independent and |cj| = 1, the variance of Vs(t, ν) becomes the sum of the variances of V(j)(t, ν), j = 1, 2, 3, 4. Moreover, since each V(j)(t, ν) follows an exponential distribution, we have:
On the other hand, μpeak, the mean of the spectrum at the peak, is approximated by Equation (43) as:
Appendix D. Optimal Amplitude Ratio between Short and Long Pulses
The maximization of the SFR of Equation (76) is equivalent to the minimization of:
Conflicts of Interest
The authors declare no conflicts of interest.
References
- Horiguchi, T.; Kurashima, T.; Tateda, M. Tensile strain dependence of Brillouin frequency shift in silica optical fibers. IEEE Photon. Technol. Lett. 1989, 1, 107–108. [Google Scholar]
- Horiguchi, T.; Shimizu, K.; Kurashima, T.; Tateda, M.; Koyamada, Y. Development of a distributed sensing technique using Brillouin scattering. J. Lightw. Technol. 1995, 13, 1296–1302. [Google Scholar]
- Horiguchi, T.; Tateda, M. BOTDA–Nondestructive measurement of single-mode optical fiber attenuation characteristics using Brillouin interaction: Theory. J. Lightw. Technol. 1989, 7, 1170–1176. [Google Scholar]
- Kurashima, T.; Tsuneo Horiguchi Hisashi Izumita, S.I.F.; Koyamada, Y. Brillouin optical-fiber time domain reflectometry. IEICE Trans. Commun. 1993, E76-B, 382–390. [Google Scholar]
- Horiguchi, T.; Shimizu, K.; Kurashima, T.; Koyamada, Y. Advances in distributed sensing techniques using Brillouin scattering. Proceedings of the SPIE, Munich, Germany, 19 June 1995; Volume 2507, pp. 126–137.
- Kurashima, T.; Tateda, M.; Horiguchi, T.; Koyamada, Y. Performance improvement of a combined OTDR for distributed strain and loss measurement by randomizing the reference light polarization state. IEEE Photon. Technol. Lett. 1997, 9, 360–362. [Google Scholar]
- Bao, X.; Brown, A.; DeMerchant, M.; Smith, J. Characterization of the Brillouin-loss spectrum of single-mode fibers by use of very short (10-ns) pulses. Opt. Lett. 1999, 24, 510–512. [Google Scholar]
- Kishida, K.; Nishiguchi, K. Pulsed pre-pump method to achieve cm-order spatial resolution in Brillouin distributed measuring technique Technical Report of IEICE, OFT2004-47 (2004-10). 2004, 15–20.
- Nishiguchi, K.; Kishida, K. Perturbation analysis of stimulated Brillouin scattering in an optical fiber and a resolution improvement method. Proceedings of the 37th ISCIE International Symposium on Stochastic Systems Theory and Its Applications, (SSS2005), Osaka, Japan, 28–29 October 2005; pp. 232–239.
- Zou, L.; Bao, X.; Wan, Y.; Chen, L. Coherent probe-pump-based Brillouin sensor for centimeter-crack detection. Optics Letters 2005, 30, 370–372. [Google Scholar]
- Brown, A.W.; Colpitts, B.G.; Brown, K. Dark-pulse Brillouin optical time-domain sensor with 20-mm spatial resolution. J. Lightw. Technol. 2007, 25, 381–386. [Google Scholar]
- Thévenaz, L.; Foaleng, S.M. Distributed fiber sensing using Brillouin echoes. Proceedings of the 19th Internetional Conference Optical Fiber Sensors, SPIE, Perth, WA, Australia, 14–18 April 2008; Volume 7004, pp. N1–N4.
- Li, W.; Bao, X.; Li, Y.; Chen, L. Differential pulse-width pair BOTDA for high spatial resolution sensing. Opt. Express 2008, 16, 21616–21625. [Google Scholar]
- Horiguchi, T.; Muroi, R.; Iwasaka, A.; Wakao, K.; Miyamoto, Y. BOTDA utilizing phase-shift pulse. IEICE Trans. Commun. 2008, J91-B, 207–216. (in Japanese). [Google Scholar]
- Koyamada, Y.; Sakairi, Y.; Takeuchi, N.; Adachi, S. Novel technique to improve spatial resolution in Brillouin optical time-domain reflectometry. IEEE Photon. Technol. Lett. 2007, 19, 1919–1912. [Google Scholar]
- Boyd, R.W.; Rzazewski, K.; Narum, P. Noise initiation of stimulated Brillouin scattering. Phys. Rev. A 1990, 42, 5514–5521. [Google Scholar]
- Gaeta, A.L.; Boyd, R.W. Stochastic dynamics of stimulated Brillouin scattering in an optical fiber. Phys. Rev. A 1991, 44, 3205–3209. [Google Scholar]
- Jenkins, R.B.; Sova, R.M.; Joseph, R.I. Steady-state noise analysis of spontaneous and stimulated Brillouin scattering in optical fibers. J. Lightw. Technol. 2007, 25, 763–770. [Google Scholar]
- Nishiguchi, K. Analytical solution to equations of Brillouin optical time domain reflectometry and its properties. Trans. ISCIE 2012, 25, 375–381. [Google Scholar]
- Nishiguchi, K. Analysis of parameter estimation error for Brillouin distributed sensing. Proceedings of the 45th ISCIE International Symposium on Stochastic Systems Theory and Its Applications, 2013 (SSS'13), Okinawa, Japan, 1–2 Novenber 2013. submitted for publication.
D1 [ns] | D2 (ns) | |||
---|---|---|---|---|
30 | 40 | 50 | 60 | |
1 | 0.0896 | 0.0750 | 0.0656 | 0.0589 |
2 | 0.1499 | 0.1254 | 0.1096 | 0.0985 |
3 | 0.2020 | 0.1690 | 0.1477 | 0.1327 |
4 | 0.2493 | 0.2085 | 0.1823 | 0.1638 |
5 | 0.2930 | 0.2451 | 0.2143 | 0.1925 |
6 | 0.3342 | 0.2795 | 0.2444 | 0.2195 |
7 | 0.3731 | 0.3121 | 0.2728 | 0.2451 |
8 | 0.4102 | 0.3431 | 0.3000 | 0.2695 |
9 | 0.4458 | 0.3728 | 0.3260 | 0.2929 |
10 | 0.4800 | 0.4014 | 0.3510 | 0.3153 |
© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license ( http://creativecommons.org/licenses/by/3.0/).
Share and Cite
Nishiguchi, K.; Li, C.-H.; Guzik, A.; Kishida, K. Synthetic Spectrum Approach for Brillouin Optical Time-Domain Reflectometry. Sensors 2014, 14, 4731-4754. https://doi.org/10.3390/s140304731
Nishiguchi K, Li C-H, Guzik A, Kishida K. Synthetic Spectrum Approach for Brillouin Optical Time-Domain Reflectometry. Sensors. 2014; 14(3):4731-4754. https://doi.org/10.3390/s140304731
Chicago/Turabian StyleNishiguchi, Ken'ichi, Che-Hsien Li, Artur Guzik, and Kinzo Kishida. 2014. "Synthetic Spectrum Approach for Brillouin Optical Time-Domain Reflectometry" Sensors 14, no. 3: 4731-4754. https://doi.org/10.3390/s140304731
APA StyleNishiguchi, K., Li, C. -H., Guzik, A., & Kishida, K. (2014). Synthetic Spectrum Approach for Brillouin Optical Time-Domain Reflectometry. Sensors, 14(3), 4731-4754. https://doi.org/10.3390/s140304731