Next Article in Journal
Piezoelectric Energy Harvesting Solutions
Next Article in Special Issue
Recent Advances in Gas and Chemical Detection by Vernier Effect-Based Photonic Sensors
Previous Article in Journal
Knitted Strain Sensors: Impact of Design Parameters on Sensing Properties
Previous Article in Special Issue
Post-Impact Fatigue Damage Monitoring Using Fiber Bragg Grating Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synthetic Spectrum Approach for Brillouin Optical Time-Domain Reflectometry

1
Office for University-Industry Collaboration, Osaka University, Yamadaoka 2-8, Osaka 565-0871, Japan
2
Neubrex Co., Ltd., Sakaemachi-dori 1-1-24, Kobe 650-0023, Japan
*
Author to whom correspondence should be addressed.
Sensors 2014, 14(3), 4731-4754; https://doi.org/10.3390/s140304731
Submission received: 31 December 2013 / Revised: 13 February 2014 / Accepted: 26 February 2014 / Published: 7 March 2014
(This article belongs to the Special Issue Photonic Sensors for Industrial, Environmental and Health Monitoring)

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 [814]. 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 [811] 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 [1214] 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 [1618]:

( 1 υ g t + z + α 2 ) E p = i κ ρ E s
( 1 υ g t z + α 2 ) E s = i κ ρ * E p
ρ t + ( Γ B 2 + 2 π i ν B ( z ) ) ρ = i E p E s * + R ( z , t )
where Ep, Es and ρ are complex amplitudes of a pump light, a backscattered light and an acoustic wave, respectively, κ and Λ are coupling coefficients, υg is the light speed in the fiber, α is the fiber loss factor and superscript * denotes a complex conjugate. ΓB is the linewidth of the acoustic wave spectrum and νB(z) denotes the BFS at position z. Here, the spectral linewidth is defined by the full width at half maximum.

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′)] = (zz′)δ(tt′), 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:

( 1 υ g t + z + α 2 ) E p = 0
( 1 υ g t z + α 2 ) E s = i κ ρ * E p
ρ t + ( Γ + 2 π i ν B ( z ) ) ρ = R ( z , t )
The boundary conditions of Ep(z, t) and Es(z, t) are given by:
E p ( 0 , t ) = P p A eff f ( t )
E s ( z , z υ g ) = 0
where Pp and f(t) denote the power and shape function of the pump pulse injected into an optical fiber, respectively, Aeff is the effective core area of a fiber and Γ = ΓB/2 is set.

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:

E p ( z , t ) = P p A eff f ( t z υ g )
Next, the stationary solution to Equation (6) is represented as:
ρ ( z , t ) = t e ( Γ + 2 π i ν B ( z ) ) ( t s ) R ( z , s ) d s
whose autocorrelation function is given by:
E [ ρ ( z , t ) ρ * ( z , t ) ] = Q 2 Γ δ ( z z ) e 2 π i ν B ( z ) ( t t ) e Γ | t t |
Then, substituting (9) and (10) into (5) and solving it under Condition (8), we obtain:
E s ( z , t ) = i κ 1 z L f f ( t 2 z z υ g ) ρ * ( z , t z z υ g ) d z
where κ 1 = P p / A eff κ and Lf is the length of the fiber.

The backscattered light returned to the input end of an optical fiber in BOTDR is represented as:

X ( t ) = def E s ( 0 , t ) = i κ 1 0 L f f ( t 2 z υ g ) ζ * ( z , t ) d z
where we set ζ(z, t) ≡ ρ(z, tz/υg), which has the same statistical property as ρ(z, t); i.e.,
E [ ζ ( z , t ) ζ * ( z , t ) ] = Q 2 Γ δ ( z z ) e 2 π i ν B ( z ) ( t t ) e Γ | t t |
holds. We note that X(t) becomes a circular complex Gaussian (ccG) process with mean zero. The component of this signal with frequency ν is obtained by:
Y ( t , ν ) = c Y h ( t ) * [ X ( t ) e 2 π i ν t ] = i κ 2 h ( t τ ) e 2 π i ν τ 0 L f f ( τ 2 z υ g ) ζ * ( z , τ ) dzd τ
where cY is a constant, h(t) is an impulse response of a low-pass filter and * denotes a convolution, and we set κ2 = cY κ1. Y(t, ν) also becomes a ccG process with a mean of zero as X(t).

2.4. Brillouin Spectrum and Point Spread Function

The Brillouin spectrum obtained by one measurement of BOTDR is represented as:

V ( t , ν ) = | Y ( t , ν ) | 2
Since BOTDR's signal source is thermally excited acoustic waves, its spectrum fluctuates even without observation noise. Although the spectrum is averaged by many replicate measurements, we must suppress the fluctuation for one measurement.

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:

E [ V ( t , ν ) E V ( t , ν ) ] 2 = ( E V ( t , ν ) ) 2
To obtain a smooth spectrum, measurements are performed many times and the spectrums are accumulated or averaged. Therefore, the Brillouin spectrum obtained by measurements is considered to be expectation EV(t, ν).

We can calculate the expectation of V(t, ν) by using the statistical property of Equation (11) as:

E V ( t , ν ) = γ R L ( t , ν ) * t , ν Ψ ( t , ν )
where γ R = υ g κ 2 2 Q / 2 Γ 2 is a constant and * t , ν denotes a two-dimensional convolution with respect to t and ν. Functions L(t, ν) and ψ(t, ν) are defined by:
L ( t , ν ) = Γ 2 Γ 2 + ( 2 π ( ν ν B ( υ g t / 2 ) ) 2
Ψ ( t , ν ) = | F τ [ f ( τ ) h ( t τ ) ] | 2
respectively, where Sensors 14 04731i1 denotes the Fourier transform with respect to τ. Equation (18) is derived in Appendix A.

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):

2 Δ T · 2 Δ F 2 π
the product of time and frequency widths of the PSF cannot be reduced below a certain value. Since the time width corresponds to the spatial resolution, there is a limit, in principle, to improving both spatial and frequency resolution simultaneously through one measurement with whatever pump pulses and low-pass filters are used.

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:

f 1 ( t ) = I [ t 0 , t 0 + D 1 ] ( t )
f 2 ( t ) = r I [ 0 , D 2 ] ( t )
where I[a,b] is the definition function of interval [a, b], D1 and D2 are the pulse widths of two pulses and r is the amplitude ratio of the two pulses. The start time of the short pulse is denoted by t0.

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:

h 1 ( t ) = f 1 ( t )
h 2 ( t ) = f 2 ( t )

By using the above elements, a pump light and a low-pass filter are composed as:

f θ ( t ) = f 1 ( t ) + e i θ f 2 ( t )
h ϕ ( t ) = h 1 ( t ) + e i ϕ h 2 ( t )
where θ and ϕ are the phase differences of the two elements. The PSF corresponding to the composite pump light and low-pass filter becomes:
Ψ θ , ϕ ( t , ν ) = | F τ { f θ ( τ ) h ϕ ( t τ ) } | 2 = | F 11 | 2 + | F 12 | 2 + | F 21 | 2 + | F 22 | 2 + 2 [ e i θ ( F 11 F 21 * + F 12 F 22 * ) + e i ϕ ( F 11 F 12 * + F 21 F 22 * ) + e i ( θ ϕ ) F 12 F 21 * + e i ( θ + ϕ ) F 11 F 22 * ]
where ℜ[·] denotes the real part, and the RHS is composed of functions:
F k l F k l ( t , ν ) = F τ [ f k ( τ ) h l ( t τ ) ] = F τ [ f k ( τ ) f l ( τ t ) ] , k , l = 1 , 2
The spreadness of function Fkl has the following features, due to the difference in pulse widths of the two pulses (see also Figure 3):
F 11 : small in time and large in frequency , F 12 and F 21 : large in both time and frequency , F 22 : large in time and small in frequency .

Therefore, among the terms of the RHS of Equation (28), only the term, F 11 F 22 *, 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 ( F 11 F 22 * ) from the PSF, we prepare p phase pairs ((θj, ϕj), j = 1, 2, ⋯, p) and p pairs of pump lights and low-pass filters:

f ( j ) ( t ) = f 1 ( t ) + e i θ j f 2 ( t )
h ( j ) ( t ) = h 1 ( t ) + e i ϕ j h 2 ( t ) j = 1 , 2 , , p
The Brillouin spectrum and the PSF obtained by using f(j)(t) and h(j)(t) are denoted by V(j)(t, ν) and Ψ(j)(t, ν), respectively The synthetic spectrum is defined as:
V S ( t , ν ) = j = 1 p c j V ( j ) ( t , ν )
where cj, j = 1, 2, ⋯, p are weighting coefficients. The corresponding synthetic PSF becomes:
Ψ S ( t , ν ) = j = 1 p c j Ψ ( j ) ( t , ν )
The j-th PSF is represented as:
Ψ ( j ) ( t , ν ) = | F τ [ f ( j ) ( τ ) h ( j ) ( t τ ) ] | 2
which is expanded as:
Ψ ( j ) ( t , ν ) = | F 11 | 2 + | F 12 | 2 + | F 21 | 2 + | F 22 | 2 + 2 [ e i θ j ( F 11 F 21 * + F 12 F 22 * ) + e i ϕ j ( F 11 F 12 * + F 21 F 22 * ) + e i ( θ j ϕ j ) F 12 F 21 * + e i ( θ j + ϕ j ) F 11 F 22 *
If the following relation holds by combining the PSFs, an ideal PSF is obtained.
Ψ S ( t , ν ) = 2 p [ F 11 F 22 * ]

This problem is equivalent to finding θj, ϕj, cj, j = 1, 2, ⋯, p that satisfy:

( 1 1 1 e i θ 1 e i θ 2 e i θ p e i ϕ 1 e i ϕ 2 e i ϕ p e i ( θ 1 ϕ 1 ) e i ( θ 2 ϕ 2 ) e i ( θ p ϕ p ) e i ( θ 1 + ϕ 1 ) e i ( θ 2 + ϕ 2 ) e i ( θ p + ϕ p ) ) ( c 1 c 2 c p ) = ( 0 0 0 0 p )
Though the general solution has not been obtained, the following solution with p = 4 is at least a solution to Equation (37):
( θ 1 θ 2 θ 3 θ 4 ) = ( 0 π / 2 π 3 π / 2 ) , ( ϕ 1 ϕ 2 ϕ 3 ϕ 4 ) = ( 0 3 π / 2 0 π / 2 ) , ( c 1 c 2 c 3 c 4 ) = ( 1 1 1 1 )
Hereafter, we use this solution.

The PSF by the synthetic approach is represented as:

Ψ S ( t , ν ) = 8 [ F 11 F 22 * ]
and is concentrated around the origin in both the time and frequency directions, as shown in the previous section. However, in addition to this, its two-dimensional integral in the time-frequency space must be a certain positive value to approximate a two-dimensional δ-function. The two-dimensional integral of the synthetic PSF of Equation (39) becomes:
Ψ S ( t , ν ) d ν d t = 8 ( f 1 ( t ) f 2 ( t ) d t ) 2
If two pulses are separated in time, this integral becomes zero, and if they overlap partially, the integral is proportional to the width of the overlapping parts. To maximize the integral, the short pulse must be included in the long pulse, and the integral becomes 8 r 2 D 1 2.

For two pulses of Equations (22) and (23), the synthetic PSF becomes:

Ψ S ( t , ν ) = 8 r 2 cos π ( D 2 D 1 2 t 0 ) ν sin π ( D 1 | t | ) ν π ν sin π ( D 2 | t | ) ν π ν I [ D 1 , D 1 ] ( t )
If D1 is sufficiently small and D2 is sufficiently large, the PSF is approximated as:
Ψ S ( t , ν ) 8 r 2 D 1 2 δ ( ν )
and the synthetic spectrum is approximated as:
E V S ( t , ν ) 8 r 2 D 1 2 γ R L ( t , ν )

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 ≤ t0D2D1, 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:

C ( ν ; t 0 ) ν ν Ψ ( t , ν ) dtd ν / ( 8 r 2 D 1 2 )
If C(ν; t0) approaches one for small ν, the PSF is considered to be concentrated near the origin of the frequency axis.

The values of C(ν; t0) are plotted for different t0's in Figure 8. Here, for t0 = (D2D1)/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 = (D2D1)/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 = (D2D1)/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:

SFR = μ peak 2 σ peak 2
where μpeak and σ peak 2 are the mean and the variance of the synthetic spectrum at the peak, respectively As shown in Appendix C, the SFR is represented explicitly as a function of r, and the optimal r, which maximizes the SFR, is given by:
r opt = ( ( Γ D 1 ) 2 2 Γ D 1 + 2 ( 1 e Γ D 1 ) ( Γ D 2 ) 2 2 Γ D 2 + 2 ( 1 e Γ D 2 ) ) 1 / 4
Numerical examples of ropt are shown in Table 1. We used the optimal r in simulations and experiments.

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

Here, the values of r were determined to be the optimal values. The strain estimates in each case are plotted in Figure 16a–c. We attained a spatial resolution of 50 cm in Case 1, 20 cm in Case 2 and 10 cm in Case 3. Figures 16a-c also indicate that the fluctuation of estimates increases as D1 becomes small. This is because the estimation error variance is nearly proportional to 1/D1, as is described in Section 3.4. Therefore, some compromise is required to select small D1. The synthetic spectrums at the 5-m strain interval are shown in Figure 16d,e. In all cases, their linewidths were less than 40 MHz. This is because the linewidth of the synthetic spectrums is determined by the reciprocal of the long pulse's length. These experimental results prove that we attained 10-cm spatial resolution that exceeds the limit of conventional BOTDR.

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:

V ( t , ν ) = κ 2 2 h ( t τ ) h * ( t τ ) e 2 π i ν ( τ τ ) × 0 L f 0 L f f ( τ 2 z υ g ) f * ( τ 2 z υ g ) ζ * ( z , τ ) ζ ( z , τ ) dzd z d τ d τ
whose expectation is calculated using Equation (14), such that:
E V ( t , ν ) = κ 2 2 Q 2 Γ h ( t τ ) h * ( t τ ) e 2 π i ν ( τ τ ) × 0 L f f ( τ 2 z υ g ) f * ( τ 2 z υ g ) e 2 π i ν B ( z ) ( τ τ ) e Γ | τ τ | dzd τ d τ
We assume that functions that contain νB(z) vanish outside the interval [0, Lf]. The integral interval can then be replaced by [−∞, ∞], and we have:
E V ( t , ν ) = υ g κ 2 2 Q 4 Γ f ( τ s ) h ( t τ ) × f * ( τ s ) h * ( t τ ) e 2 π i ( ν ν B ( υ g s / 2 ) ) ( τ τ ) e Γ | τ τ | dsd τ d τ = γ R f ( τ s ) h ( t τ ) × f * ( τ τ s ) h * ( t τ + τ ) d τ Γ 2 e 2 π i ( ν ν B ( υ g s / 2 ) ) τ e Γ | τ | dsd τ = γ R e 2 π i ν τ { [ f ( τ s ) h ( t τ ) ] * τ [ f * ( τ s ) h * ( t + τ ) ] } Γ 2 e 2 π i ν B ( υ g s / 2 ) τ e Γ | τ | dsd τ
where we set γ R = υ g κ 2 2 Q / 2 Γ 2. Since the Fourier transform of a product is the convolution of the corresponding transforms, we have:
E V ( t , ν ) = γ R { F τ [ f ( τ s ) h ( t τ ) ] F τ [ f * ( τ s ) h * ( t + τ ) ] } * ν L ( s , v ) d s = γ R | F τ [ f ( τ s ) h ( t τ ) ] | 2 * ν L ( s , v ) d s = γ R | F τ [ f ( τ ) h ( t s τ ) ] | 2 * ν L ( s , v ) d s = γ R | F τ [ f ( τ ) h ( t τ ) ] | 2 * t , ν L ( s , v )
where L(t, ν) is defined by Equation (19) and is the Fourier transform of (Γ/2) exp(2πiνB(υgt/2)τ − Γ|τ|) with respect to τ. Moreover, defining Ψ(t, ν) by Equation (20) yields Equation (18).

Appendix B. Uncertainty Relation of PSF

The frequency integral of the PSF Equation (20) gives a temporal PSF as follows:

Ψ T ( t ) Ψ ( t , ν ) d ν = | f ( t ) | 2 * | h ( t ) | 2
Center μT and radius ΔT of ΨT(t) are obtained as:
μ T = def 1 N t Ψ T ( t ) d t = μ f + μ h
Δ T 2 = def 1 N ( t μ T ) 2 Ψ T ( t ) d t = σ f 2 + σ h 2
where:
N = Ψ T ( t ) d t
is a normalizing constant and:
μ f = 1 f 2 t | f ( t ) | 2 d t
σ f 2 = 1 f 2 ( t μ f ) 2 | f ( t ) | 2 d t
Here, ‖ · ‖ denotes a L2 norm, and μh and σh are defined in the same way. Similarly, a frequency PSF is obtained by a time integral as follows:
Ψ F ( ν ) Ψ ( t , ν ) d t = | f ^ ( ν + ξ ) | 2 · | h ^ ( ξ ) | 2 d ξ
where and ĥ are the Fourier images of f and h, respectively A center and a radius of ΨF(ν) are also obtained as:
μ F = μ f ^ μ h ^ , Δ F 2 = σ f ^ 2 + σ h ^ 2
where μ, μĥ, σ and σĥ are defined the same as μf and σf. As is well known in Fourier analysis, the following uncertainty relations hold:
σ f σ f ^ 1 4 π , σ h σ h ^ 1 4 π
the use of which, we have:
Δ T 2 Δ F 2 = ( σ f 2 + σ h 2 ) ( σ f ^ 2 + σ h ^ 2 ) 2 ( 4 π ) 2 + 1 ( 4 π ) 2 ( σ f 2 σ h 2 + σ h 2 σ f 2 )
Since σ f 2 / σ h 2 + σ h 2 / σ f 2 in the last term takes the minimum value of two, when σf = σh, we obtain:
Δ T 2 Δ F 2 1 ( 2 π ) 2

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:

E [ V S ( t , ν ) E V S ( t , ν ) ] 2 = j = 1 4 ( E V ( j ) ( t , ν ) ) 2
We now arbitrarily fix time t and obtain the fluctuation at the peak of the spectrum. For simplicity, we assume that the BFS is constant and νB(υgt/2) = 0. The peak position is then given by ν = 0. The expectation of all elements of the synthetic spectrum in S-BOTDR at ν = 0 is then given by:
E V ( j ) ( t , 0 ) = γ R [ L ( t , ν ) * t Ψ ( j ) ( t , ν ) ] d ν
Using the assumption for the BFS, the Lorentzian spectrum can be written as:
L ( t , ν ) = L 0 ( ν ) = def Γ 2 Γ 2 + ( 2 π ν ) 2
which, substituting into Equation (63), yields:
E V ( j ) = γ R L 0 ( ν ) Ψ ( j ) ( t , ν ) dtd ν
where we abbreviate the arguments, because EV(j) (t, 0) does not depend on t. Noting that:
L 0 ( ν ) = F τ [ Γ 2 e Γ | τ | ]
and using Plancherel's theorem and Equation (34), we obtain:
E V ( j ) = F ν 1 [ L 0 ( ν ) ] F ν 1 [ Ψ ( j ) ( t , ν ) d t ] d τ = γ R Γ 2 e Γ | τ | [ h ( j ) ( t ) h ( j ) * ( t + τ ) ] d t [ f ( j ) ( t ) f ( j ) * ( t τ ) ] dtd τ = γ R Γ 2 e Γ | τ | ( K 11 + r e i θ j K 12 + r e i θ j K 21 + r 2 K 22 ) × ( K 11 + r e i ϕ j K 12 + r e i ϕ j K 21 + r 2 K 22 ) d τ
where:
K k l = f k ( t ) f l ( t τ ) d t = h k ( t ) h l ( t + τ ) d t , k , l = 1 , 2
Substituting the values of Equation (38) into Equation (67), we have:
E V ( j ) = γ R Γ 2 ( a 0 + a 1 , j r 2 + a 2 r 4 )
where r is the amplitude ratio of the short and long pulses appearing in Equation (23), and:
a 0 = e Γ | τ | K 11 2 d τ
a 1 , j = 2 e Γ | τ | ( K 11 K 22 + cos ( θ j + ϕ j ) K 12 2 K 12 K 21 ) d τ
a 2 = e Γ | τ | K 22 2 d τ
From the relation of Equation (62), σ peak 2, the variance of the spectrum at the peak, is represented as:
σ peak 2 = j = 1 4 ( E V ( j ) ) 2
into which substituting Equation (63) gives:
σ peak 2 = ( γ R Γ 2 ) 2 [ 4 ( a 0 + a 1 r 2 + a 2 r 4 ) 2 + j = 1 4 ( a 1 , j 2 a 1 2 ) r 4 ]
where we set a 1 = j = 1 4 a 1 , j / 4.

On the other hand, μpeak, the mean of the spectrum at the peak, is approximated by Equation (43) as:

μ peak = max E V S ( t , ν ) 8 r 2 D 1 2 γ R
Substituting Equations (75) and (74) into Equation (45), we have:
SFR = ( 8 D 1 2 / Γ ) 2 ( a 0 r 2 + a 1 + a 2 r 2 ) 2 + 1 4 j = 1 4 ( a 1 , j 2 a 1 2 )

Appendix D. Optimal Amplitude Ratio between Short and Long Pulses

The maximization of the SFR of Equation (76) is equivalent to the minimization of:

J = a 0 r 2 + a 1 + a 2 r 2
The value of r that minimizes J is given by:
r opt = ( a 0 a 2 ) 1 / 4
which maximizes SFR. From Equations (22), (23) and (68),
k 11 ( τ ) = { D 1 | τ | , | τ | D 1 0 , otherwise
k 22 ( τ ) = { D 2 | τ | , | τ | D 2 0 , otherwise
holds, and by substituting them into Equations (70) and (72), we have:
a 0 = 2 Γ 3 [ ( Γ 2 D 1 ) 2 2 Γ D 1 + 2 ( 1 e Γ D 1 ) ]
a 2 = 2 Γ 3 [ ( Γ 2 D 2 ) 2 2 Γ D 2 + 2 ( 1 e Γ D 2 ) ]
Therefore,
r opt = ( ( Γ D 1 ) 2 2 Γ D 1 + 2 ( 1 e Γ D 1 ) ( Γ D 2 ) 2 2 Γ D 2 + 2 ( 1 e Γ D 2 ) ) 1 / 4
becomes the optimal value of r.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. 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]
  2. 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]
  3. 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]
  4. 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]
  5. 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.
  6. 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]
  7. 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]
  8. 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.
  9. 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.
  10. 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]
  11. 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]
  12. 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.
  13. 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]
  14. 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]
  15. 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]
  16. Boyd, R.W.; Rzazewski, K.; Narum, P. Noise initiation of stimulated Brillouin scattering. Phys. Rev. A 1990, 42, 5514–5521. [Google Scholar]
  17. 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]
  18. 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]
  19. Nishiguchi, K. Analytical solution to equations of Brillouin optical time domain reflectometry and its properties. Trans. ISCIE 2012, 25, 375–381. [Google Scholar]
  20. 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.
Figure 1. Block diagram of Brillouin optical time-domain reflectometry (BOTDR).
Figure 1. Block diagram of Brillouin optical time-domain reflectometry (BOTDR).
Sensors 14 04731f1 1024
Figure 2. Pump lights and matched filters: (a) Element 1 of the pump light; (b) Element 2 of the pump light; (c) Element 1 of the low-pass filter; and (d) Element 2 of the low-pass filter.
Figure 2. Pump lights and matched filters: (a) Element 1 of the pump light; (b) Element 2 of the pump light; (c) Element 1 of the low-pass filter; and (d) Element 2 of the low-pass filter.
Sensors 14 04731f2 1024
Figure 3. Illustration of Fkl, k, l = 1, 2 spreadness. |Fkl|'s take large values in the black regions: (a) F11; (b) F12; (c) F21; and (d) F22.
Figure 3. Illustration of Fkl, k, l = 1, 2 spreadness. |Fkl|'s take large values in the black regions: (a) F11; (b) F12; (c) F21; and (d) F22.
Sensors 14 04731f3 1024
Figure 4. Point spread functions: (a) conventional BOTDR (D = 10 ns); (b) conventional BOTDR (D = 1 ns); (c) synthetic (S)-BOTDR (D1 = 1 ns; and D2 = 50 ns).
Figure 4. Point spread functions: (a) conventional BOTDR (D = 10 ns); (b) conventional BOTDR (D = 1 ns); (c) synthetic (S)-BOTDR (D1 = 1 ns; and D2 = 50 ns).
Sensors 14 04731f4 1024
Figure 5. Composition of the synthetic point spread function of S-BOTDR.
Figure 5. Composition of the synthetic point spread function of S-BOTDR.
Sensors 14 04731f5 1024
Figure 6. Block diagram of S-BOTDR.
Figure 6. Block diagram of S-BOTDR.
Sensors 14 04731f6 1024
Figure 7. Composition of the synthetic spectrum of S-BOTDR.
Figure 7. Composition of the synthetic spectrum of S-BOTDR.
Sensors 14 04731f7 1024
Figure 8. C(ν; t0): the integral of the point spread function (PSF).
Figure 8. C(ν; t0): the integral of the point spread function (PSF).
Sensors 14 04731f8 1024
Figure 9. Signal-to-fluctuation ratio for various D1 and D2.
Figure 9. Signal-to-fluctuation ratio for various D1 and D2.
Sensors 14 04731f9 1024
Figure 10. Simulation condition: (a) assumed Brillouin frequency shift (BFS) and (b) the Lorentzian spectrum.
Figure 10. Simulation condition: (a) assumed Brillouin frequency shift (BFS) and (b) the Lorentzian spectrum.
Sensors 14 04731f10 1024
Figure 11. Simulation results of conventional BOTDR with D = 10 ns: (a) the Brillouin spectrum and (b) the estimate of BFS.
Figure 11. Simulation results of conventional BOTDR with D = 10 ns: (a) the Brillouin spectrum and (b) the estimate of BFS.
Sensors 14 04731f11 1024
Figure 12. Simulation results of conventional BOTDR with D = 1 ns: (a) the Brillouin spectrum and (b) the estimate of BFS.
Figure 12. Simulation results of conventional BOTDR with D = 1 ns: (a) the Brillouin spectrum and (b) the estimate of BFS.
Sensors 14 04731f12 1024
Figure 13. Simulation results of S-BOTDR with D1 = 1 ns and D2 = 50 ns: (a) the composite Brillouin spectrum and (b) the estimate of BFS.
Figure 13. Simulation results of S-BOTDR with D1 = 1 ns and D2 = 50 ns: (a) the composite Brillouin spectrum and (b) the estimate of BFS.
Sensors 14 04731f13 1024
Figure 14. The fiber under testing used in experiment: (a) the fiber setup and (b) the equivalent strain.
Figure 14. The fiber under testing used in experiment: (a) the fiber setup and (b) the equivalent strain.
Sensors 14 04731f14 1024
Figure 15. Experimental BOTDR results: (a) strain estimates and (b) the Brillouin spectrum at the 5-m strain interval.
Figure 15. Experimental BOTDR results: (a) strain estimates and (b) the Brillouin spectrum at the 5-m strain interval.
Sensors 14 04731f15 1024
Figure 16. Experimental S-BOTDR results: (a–c) Strain estimates of Cases 1 to 3; (d–f) composite Brillouin spectrums at the 5-m strain interval of Cases 1 to 3.
Figure 16. Experimental S-BOTDR results: (a–c) Strain estimates of Cases 1 to 3; (d–f) composite Brillouin spectrums at the 5-m strain interval of Cases 1 to 3.
Sensors 14 04731f16 1024
Table 1. Numerical examples of ropt.
Table 1. Numerical examples of ropt.
D1 [ns]D2 (ns)

30405060
10.08960.07500.06560.0589
20.14990.12540.10960.0985
30.20200.16900.14770.1327
40.24930.20850.18230.1638
50.29300.24510.21430.1925
60.33420.27950.24440.2195
70.37310.31210.27280.2451
80.41020.34310.30000.2695
90.44580.37280.32600.2929
100.48000.40140.35100.3153

Share and Cite

MDPI and ACS Style

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

AMA Style

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 Style

Nishiguchi, 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 Style

Nishiguchi, 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

Article Metrics

Back to TopTop