Time and Frequency Localized Pulse Shape for Resolution Enhancement in STFT-BOTDR
Short-Time Fourier Transform-Brillouin Optical Time-Domain Reflectometry (STFT-BOTDR) implements STFT over the full frequency spectrum to measure the distributed temperature and strain along the optic fiber, providing new research advances in dynamic distributed sensing. The spatial and frequency resolution of the dynamic sensing are limited by the Signal to Noise Ratio (SNR) and the Time-Frequency (T-F) localization of the input pulse shape. T-F localization is fundamentally important for the communication system, which suppresses interchannel interference (ICI) and intersymbol interference (ISI) to improve the transmission quality in multicarrier modulation (MCM). This paper demonstrates that the T-F localized input pulse shape can enhance the SNR and the spatial and frequency resolution in STFT-BOTDR. Simulation and experiments of T-F localized different pulses shapes are conducted to compare the limitation of the system resolution. The result indicates that rectangular pulse should be selected to optimize the spatial resolution and Lorentzian pulse could be chosen to optimize the frequency resolution, while Gaussian shape pulse can be used in general applications for its balanced performance in both spatial and frequency resolution. Meanwhile, T-F localization is proved to be useful in the pulse shape selection for system resolution optimization.
Dynamic strain measurement in distributed fiber optic sensing (DFOS) is essential for structural health monitoring (SHM) of the strain changes induced by aging material and seismic or man-made activities . Since the first system was launched in 1999, due to its stability, reliability, and insensitivity to external perturbations, the distributed Brillouin sensor has become one of the most popular diagnostic tools for SHM in bridges, dams, pipelines, and piles . Among the different DFOS systems, BOTDR launches pulsed light into one end of the fiber and detects the scattered spontaneous Brillouin signal at the same end [3, 4]. Single-end detection is convenient for fiber installation and on-site measurement. However, the large averaging number for a good Signal to Noise Ratio (SNR) and the slow frequency scanning method used in classical BOTDR have become obstacles to achieving dynamic measurement .
Recently, a fast-speed Short-Time Fourier Transform (STFT) based BOTDR applies a wideband detection architecture to replace the time-consuming frequency scanning method, which provides an opportunity for dynamic measurement. In the literature, the system successfully achieved a 16.7 Hz strain variation on a 12 m fiber section with 4 m spatial resolution and 45 με strain uncertainty . The spatial resolution is the minimum detectable length of strain change and has been studied over the past decades to detect small cracks , while the strain uncertainty is the strain resolution for the measurement. A longer pulse or a longer rising/falling time will increase the interaction of the light photons and acoustic phonons and lead to a worse spatial resolution along the fiber .
BOTDR measures the peak frequency shift of the Brillouin scattering spectrum. SNR is the detective power density at spectrum’s peak divided by the power of the whole spectrum including the sidebands. The SNR of the detected Brillouin scattering spectrum at each point along the fiber reduces when the sideband of the pulse’s spectrum increases. Increased SNR will reduce average time of the measurement to enhance the speed of the analyzer to achieve real-site dynamic measurement . SNR can be affected by the pulse power  and the pulse shape, which has significant impact on the spatial resolution, temperature, and strain resolution of the system [10–12].
Previous literatures compared time-domain and frequency-domain analysis independently to evaluate the analyzer performance. In 1999, Naruse and Tateda published the relation between the pulse width and the Brillouin backscattered light power and mentioned the trade-off between spatial resolution and frequency resolution . In 2000, the same group announced that different pulse shapes would have different bandwidth of the spectrum and have different frequency-measurement error according to simulations . Hao et al. in 2013 compared three pulse shapes including trapezoidal, triangular, and rectangular pulses and expressed that the Brillouin spectrum power of the triangular pulse shape was higher than the power of other pulse shapes [10, 14]. They also compared other pulse shapes, including Lorentzian, Gaussian, super Gaussian, and triangular shape, and concluded that the Lorentzian shape generated the biggest peak power of Brillouin spectrum and has the measurement accuracy when the Brillouin spectrum peak is set to the same level .
In the above previous studies, assumptions were made that the spatial resolution was identical to the pulse width (full width at half maximum (FWHM) of the pulse) and the Brillouin spectrum was affected by frequency-domain information. Spatial resolution can be enhanced by reducing the pulse edges (in time-domain) and frequency resolution can be enhanced by reducing the sideband in pulse spectrum (in frequency-domain). However, time-domain and frequency-domain information are identically important, which should be considered simultaneously to offer the design trade-off analysis between pulse width and frequency linewidth. For instance, when the frequency shifts (due to strain or temperature variation) are narrower than the FWHM, the spectrum will have double peaks and lead to complication in comparison , in which time- and frequency-domains’ analysis needs to be considered cooperatively.
The Time-Frequency (T-F) localization has been utilized in communication and signal processing [16–18], which generates sharp localization, offers small distortion, and reduces the interference such as intercarrier interference (ICI) and intersymbol interference (ISI) in OFDM communication. The resolution in frequency and time-domain in image processing and computer visions can be enhanced in real applications using T-F localization [19–22].
We introduce the T-F localization into the design consideration of the BOTDR system. In this paper, Time-Domain Effective-Pulse Width (TEW) and Frequency-Domain Effective-Pulse Linewidth (FEL) are defined to offer a comprehensive view to optimize the analyzer performance. Different pulses with the same FWHM may have different TEW and/or different SNR requirement in the detection and the FEL can be different. In particular, the variance in laser output power, noise level, or averaging number, and so forth, in the same pulse shape, can result in different spatial and frequency resolution. Similarly, the pulses with various effective linewidth can lead to different noise immune ability in the system.
This paper demonstrates the importance for T-F localized input pulses to improve the SNR of the Brillouin scattering signal, which leads to a narrower frequency bandwidth and an efficient pulse length to enhance the spatial and frequency resolution. The pulse localization model, based on the physical description of BOTDR and time-frequency analysis of the pulses with the definition of TEW and FEL, has been built to investigate the performance of various pulse shapes, such as Lorentzian, rectangular, Gaussian, and triangular pulses, which are compared theoretically and experimentally according to their performances in the time-domain (spatial resolution) and the frequency-domain (pulse modulated Brillouin spectrum linewidth). The comparison of the four pulse shapes offered matched results in simulation and experiment in the same system and improved the understanding of detection of minimum frequency-change-length and frequency accuracy in different pulse shapes with the same FWHM. Hence the Gaussian pulse can be used in general measurement because of its balanced time and frequency resolution. Practically, the combination of different pulse shapes can be used to obtain accurate measurements in time- and frequency-domain, respectively. To enhance the analyzer performance, the methodology of the T-F localization method is a useful tool to design, compare, and implement the pulse and algorithm selection in DFOS. T-F localization should be also considered in general to evaluate the power, averaging number settings in DFOS development.
2. Physical and Mathematical Description on BOTDR
In a single mode fiber, 1550 nm light will generate Brillouin scattering signals whose peak frequency shifts about 11 GHz, due to the interaction between the lightwave and acoustic waves in the fiber [23, 24]. The acoustic wave makes the Brillouin gain spectrum (BGS), a Lorentzian shape distribution whose linewidth is determined by the phonon lifetime . Distributed fiber sensing systems detect the frequency shift of the BGS by searching the frequency peaks over time, hence producing distributed sensing information (temperature and strain) along the fiber [26, 27]. The physical and mathematical model has been discussed by several research groups. Minardo et al.  constructed the solution of the BOTDA model, while Naruse and Tateda  offered a simplified BOTDR model, and Nishiguchi et al.  introduced a more general BOTDR model.
The Brillouin scattering process can be expressed by the following equations : where , , and are the amplitudes of the pump light, stokes light, and acoustic wave, respectively. and are the polarization density of the pump and stoke waves, respectively. is the electrostrictive force per unit volume.
By simplifying the equation and assuming that noise is constant along the fiber [28–30],where is the light velocity in a vacuum, is the refractive index of the fiber, is the sum of the damping of the phonon lifetime, and the detuning parameter . is the frequency shift generated by Brillouin scattering, which is proportional to the independent thermodynamic variables of the entropy, , and pressure, :Because Spontaneous Brillouin Scattering (SpBS) does not have a stimulated effect, the SBS terms in (2) and (4) are assumed to be zero to simplify the analysis.
Then the scattered light can be expressed as The spectrum can be described asThis can be written in convolution form as where is the time-nondependent coefficient and is the Fourier transform of the pulse signal .
At a particular time, according to the model given by Naruse and Tateda , the spectrum of the Brillouin scattered signal at each point along the fiber can be expressed as where is the power spectrum of a launched pulse in which each frequency component produces a Lorentzian shape Brillouin backscattered light spectrum profile with peak frequency of and full width at half maximum (FWHM) of . The term expresses the frequency changes due to local acoustic waves, where the difference comes from the changes in the properties of the fiber or the changes in strain or temperature.
Based on the equations above, the final scattered back Brillouin spectrum is the integration of all the spectra generated by each frequency component on the pump pulse spectrum. According to (8) and (11), the time-varying Brillouin spectrum is related to the entire pulse shape. In the pulse’s rising and falling edges, the spatial resolution will be reduced because the power in the edges will worsen the spatial resolution by expanding the length of wave interaction. Hence, time localized pulse can enhance the spatial resolution.
The broadened bandwidth and expanded sideband increase the complexity of accurate peak frequency detection. The noise will be increased when the sideband power is comparable with mainlobe power in the pulse modulated SpBS spectrum . In order to remain a good SNR, the sidebands’ power should be minimized. Therefore the frequency localized pulses in BOTDR require the power to be concentrated near the peak frequency to maintain a narrow linewidth.
3. Time-Frequency () Localization
The mean square deviation of the time and frequency distribution of a signal is defined in (12) and (13) . Equation (14) describes the Heisenberg uncertainty principle : where and are the variance of the signal in the time- and frequency-domains, respectively.
In the telecommunication systems, the dispersion of the signal causes intersymbol and interchannel interference when the signal has multicarriers, such as in an Orthogonal Frequency Division Multiplex (OFDM) system . When the signal disperses over time, it may also interfere with subsequent symbols to cause “intersymbol” interference (ISI), reducing the reliability of the communications system . When the signal disperses in the frequency-domains, it may interfere with symbols in other channels, causing “interchannel” interference (ICI), reducing the SNR in each channel . A delta function in the frequency-domain produces a continuous single frequency lightwave with constant power. The uncertainty principle governs that delta functions in both time- and frequency-domains cannot simultaneously exist in the same signal.
T-F localization describes the extent to which a pulse signal is restricted in both the time- and frequency-domains. Good localization describes the minimum energy spreading over a certain time span and frequency bandwidth to avoid the power leakage into channels and interfering with neighboring symbols . Besides in telecommunication, the T-F localization is also well studied in high resolution geophysical data analysis  and in wavelets transform in signal processing .
Pulse shapes have different localization characteristic. Rectangular, Lorentzian, Gaussian, and triangular shape pulses were simulated in 3D using STFT. In order to plot the time- and frequency-domain of the pulses simultaneously, the pulses were multiplied with a square STFT window to do the time-frequency analysis of the shapes. The STFT 3D plots and the sideband in both domains are shown in Figures 1 and 2. The pulses all had 50 ns FWHM and 1 μs period and were analyzed using 50 ns window STFT.
The variances obtained by (12) and (13) in time- and frequency-domain for different pulse shapes and the product of the variances are given in Table 1. The Lorentzian shape has the smallest product with the same FWHM of the normalized pulses, showing that the Lorentz shape performs the best in both time and frequency localization, using the traditional telecommunication definition.
4. Time-Frequency Localization in the BOTDR System
The product of is a measure of the performance of a pulse, where and are defined as the variance of the pulse in the time-domain and frequency-domain, respectively. In a Brillouin optic fiber sensing system, the target is no longer to reduce the ISI and ICI in channels but to enhance the efficient pulse length and narrow the linewidth of the pulse modulated Brillouin spectrum, providing a better SNR in time- and frequency-domains, so that the system can achieve long distance and high resolution measurement.
A pulse consists of two parts in time-domain, mainlobe, and tails. Conventionally, the mainlobe is the power within the FWHM width and the tails are the parts of rising and falling edges. The final pulse modulated Brillouin spectrum consists of two parts of Brillouin scattering in the fiber: one generated by the mainlobe part of the pulse and one generated by the tails of the pulse.
When a uniform strain and temperature are applied to a fiber, the center frequencies of the SpBS remain constant; thus, the tails in pulse do not affect the monitoring result. When strain or temperature changes, the frequency of SpBS changes, and Brillouin scattering generated by mainlobe and tails are different. If the changes happen in the tails’ part, the changes will affect the result of the mainlobe part. If the power in the tail is comparable to that in the mainlobe, the peaks of the two parts will be apparent from each other, reducing the frequency accuracy due to the double peak effect. Hence, the tail of pulses will reduce the spatial resolution.
In the pulses plotted in Figure 3, the power in the mainlobe, shown in red, is and the powers in the two tails, shown in blue, are and and are assumed to be the same, as. The width of mainlobe is the FWHM of the pulse. The extinction ratio is defined asThe extinction ratios of the different pulse shapes are summarized in Table 2.
In these shapes, a rectangular pulse has ideal edges in time-domain; however, in the frequency-domain, the sideband is the worst, while the Lorentzian pulse performs best but is worst in the time-domain. The Lorentzian shape has the smallest product of pulse width and frequency bandwidth in Table 1 but has the lowest extinction ratio in Table 2. The Gaussian and triangular shape have similar extinction ratio which is worse than Lorentzian shape and better than rectangular shape. The requirements for the extinction ratio and the system resolution can differ from the pulse width with the different signal processing methods [33, 34], but the T-F localization can indicate the effectiveness of the input pulses in relation to the spatial and frequency resolution in BOTDR systems.
This paper defines two parameters for BOTDR, one in time-domain—the Time-Domain Effective-Pulse Width (TEW)—and one in frequency-domain—the Frequency-Domain Effective-Pulse Linewidth (FEL). TEW determines the effective spatial resolution detected by the monitoring system. FEL is the actual spectral linewidth required to obtain a required spectrum SNR for peak detection. The mainlobe in time-domain is defined as the part of the pulse that lies within the TEW, and the tails are defined as the parts of the pulse outside TEW in time-domain. The mainlobe in frequency-domain is the frequency slots within FEL while the sideband is formed from the frequency slots outside the FEL in frequency-domain.
In the time-domain, is defined as the ratio of the power in the mainlobe to that in the tails of the pulse. If the detection system requires (unit in dB) to obtain the precise frequency shift information, the effective pulse width is that which satisfies (16); and it can be written in the integration form as in (17): where is the normalized distribution of the shape and is the TEW.
Let us assume that, for the plotted pulses with 50 ns FWHM width, the system can only detect the accurate peak when the pulse’s mainlobe power is 3 dB greater than the tail’s power. The TEW of the plotted pulses is shown in Table 3. With the same FWHM designed width, the TEW of the Lorentzian shape is more than twice as large as the rectangular pulse. The TEW of Gaussian and triangular shape is slightly wider than rectangular but still much smaller than Lorentzian shape, showing that the rectangular shapes will give the best spatial resolution.
In the frequency-domain, FEL is defined as if the system requires that the power of the spectrum peak is dB higher than the noise level; thus, the effective Brillouin linewidth is the “bandwidth” at dB lower than the spectrum peak. Because the power level of the Brillouin spectrum is proportionate to the power of the pulse spectrum, the FEL of the pulse is the bandwidth of the power level that is dB lower than the peak level, as shown in Figure 4. Hence, it can be expressed bywhere is the power included in the frequency mainlobe which has the width of FEL, is the power in the sideband of the pulse’s spectrum, and is the FEL of the pulses.
5. Simulation and Experiment on Pulse Shapes
Based on the definition of the Time-Frequency localization in the BOTDR system, the performance of Lorentzian, Gaussian, triangular, and rectangular pulse shapes will depend on and values to quantify the spatial and frequency resolution. Linear variation of the and values is simulated from −10 dB (a noise-immune system) to 10 dB (a noise-sensitive system) to evaluate TEW and FEL. Rectangular, Lorentzian, Gaussian, and triangular shape input pulses with 50 ns pulse width and the same peak voltage are experimentally chosen to demonstrate and evaluate the T-F localization in the STFT-BOTDR system.
5.1. Simulation Result
In STFT-BOTDR, and need to be optimized to reduce and , improving the solution of the detecting system. Figures 5 and 6 show the changes in TEW and FEL when system varies from −10 dB to 10 dB. The rectangular pulse has the best performance because it contains no tails but only the mainlobe signal. Other pulses show a rising trend when the system requires an increasing extinction ratio in the time-domain. Among those pulses, the Lorentzian shape pulse is the worst when is large (which means that the system will have a larger spatial resolution with Lorentzian shape).
In the frequency-domain, the changes of FEL show opposite trends compared to the time-domain. The Lorentzian shape has the smallest linewidth and shows a slow increment with increasing, whereas the rectangular pulse has the largest linewidth and expands rapidly when increases in the frequency-domain. In the two domains, the Gaussian and triangular pulse have similar performance which are always better than Lorentzian shape and worse than rectangular shape when is larger than 0 dB. The Gaussian shape is always slightly better than triangular shape though the result is very close.
Generally, and should be minimized to enhance SNR and system resolution during system development. With a low SNR system, the Gaussian pulse gives a balanced performance in both the time- and frequency-domains, while others can have better performance in either spatial resolution or frequency resolution, sacrificing the other, too.
5.2. Experimental Setup and Results
STFT-BOTDR using frequency downconversion and STFT algorithm is shown in Figure 7. An ultra-narrow-line-width laser with the wavelength of 1554.4 nm was used as the light source. In branch A, after passing a coupler, a part of the continuous light was modulated by an electrooptic modulator (EOM) with a 20 μs period and 50 ns pulse width signal generated by an Agilent 33600 arbitrary waveform generator (AWG). The pulsed light was amplified by an Erbium-doped fiber amplifier (EDFA) and then circulated into a long sensing fiber to generate the Brillouin backscattering signal. This signal was then heterodyned with the reference CW light in branch B and transferred to a RF signal using a wideband photodetector. The frequency of the signal was then downconverted to the intermediate frequency (IF) signal, which was then sampled in the time-domain and processed using the STFT algorithm to obtain the frequency-domain information—the SpBS spectrum along the fiber under test.
The fiber under test (FUT) was a 500 m single mode fiber with heating sections, immersed in a water bath at the end. The heating part is shown in Figure 8. The heating sections were of different lengths 5 m, 4 m, 3 m, 2 m, 1 m, 0.8 m, and 0.5 m. Each section was separated by a distance of 10 m. All the sections were heated to 75°C and the rest of the fiber remained at room temperature.
Rectangular, Lorentzian, Gaussian, and triangular pulses, respectively, were experimentally generated as the source into the system, shown in Figure 9, to compare the performance in terms of the spatial resolution and pulse modulated Brillouin scattering spectrum linewidth at the same position along the testing fiber. Figure 10 shows the distribution of the center frequency along the fiber which measured the temperature distribution along the FUT. Conventionally, 50 ns pulses are considered to have a 5 m spatial resolution as the pulse width is about 50 ns [4, 34]. However, in the experimental results, all pulses could detect the spatial resolution lower than 5 m because the TEW of each pulse was less than the pulse width; hence the spatial resolution is smaller than the conventional understanding by using advanced signal processing algorithms. The Lorentzian, Gaussian, and triangular pulses could only detect the change in the 3 m section; however, the rectangular shape clearly monitored the temperature changes in the 2 m section, showing that the rectangular shape has a better TEW and gives an improved spatial resolution. In the center frequency plot, the Lorentzian shape shows a drop in 5 m section and center frequency detected decreases rapidly in small temperature change sections, given that its spatial resolution is the worst among the four shapes. The Gaussian pulse is slightly better than triangular shape with less reduction in center frequency in small sections though they both have similar peak detection in 5 m and 4 m section with rectangular pulse, but their detected center frequency deviates from correct value after 4 m section as well.
In the frequency-domain, shown in Figure 11, the rectangular pulse has a wider linewidth and a larger sideband, given the widest Brillouin scattering spectrum linewidth, which worsens the SNR of the peak power. The Lorentzian pulse performed the best, with the smallest FEL in the Brillouin spectrum. The Gaussian and triangular pulse again have middle FEL but they are close to Lorentzian pulse. When the changes are smaller or close to the spatial resolution, the Brillouin spectrum linewidth depends on the property of pulses not only in frequency-domain but also in the time-domain, because the overall frequency information over the detection range of the pulse will be compacted and represented in one localized spectrum. For instance, when the TEW of a pulse is large, such as Lorentzian shape, though it has smallest linewidth, its long tail in time-domain shared the frequency information over the entire measurement distance the pulse covered. Hence the final spectrum will include separate frequency shift components, causing the double-peak effect. As a result, the Lorentzian pulse has an error in the center frequency detection over 5 m section.
Hence, with the same FWHM of the pulses, the rectangular pulse has the best spatial resolution, while the Lorentzian gives the best linewidth. The result shows that the rectangular pulse has better time localization and the Lorentzian shape has better frequency localization. The Gaussian and triangular shape have balance performance in two domains but Gaussian pulse is slightly better in the experiment.
In this paper, we discussed the importance of Time-Frequency (T-F) localization and its contribution to improving spatial and frequency resolution. T-F localization in telecommunication systems was demonstrated and its alternative definition in the BOTDR system was introduced. With the new definition of Time-Frequency localization in BOTDR, the Time-Domain Effective-Pulse Width (TEW) and Frequency-Domain Effective-Pulse Linewidth (FEL) were introduced as important parameters for resolution optimization. The conventional relationship of pulse width and spatial resolution is challenged with advanced architecture and signal processing method. Instead, TEW and FEL are more suitable to become a bridge between spatial resolution and frequency resolution and to describe the performance of a BOTDR system. Furthermore, rather merely one domain is considered in most of the papers, which is not suitable for real cases; the T-F localization method considers two domains simultaneously. Two parameters and are defined to quantify the effect of TEW and FEL in BOTDR systems, which should be optimized to enhance the SNR and spatial and frequency resolution during system development. Commonly used input pulses, such as Lorentzian, rectangular, Gaussian, and triangular shape, were simulated and experimentally demonstrated in STFT-BOTDR to compare the Signal to Noise Ratio (SNR), the T-F localization, and the system resolution. The result indicates that the rectangular pulse should be selected when spatial resolution is crucial for the application and the Lorentzian shape should be chosen for the application requiring high frequency resolution. In practical usage, the two pulses could be combined in separate measurements to obtain a strain profile with optimized spatial and frequency resolution. Gaussian shape pulse has more balanced resolution in the time- and frequency-domains, which is suitable for general applications. The T-F localization is proved to be a powerful tool in the system design and development process to compare and implement the pulse shape, algorithm, and the selection and optimization of power, averaging number and relevant settings.
Following EPSRC policy on research data, the additional data will be accessible at https://www.repository.cam.ac.uk/handle/1810/250413.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Linqing Luo and Bo Li are supported by the China Scholarship Council and Cambridge Commonwealth, European and International Trust. EPSRC EP/K000314/1 Grant is acknowledged.
A. Funnell, X. Xu, J. Yan, and K. Soga, “Simulation of noise within BOTDA and COTDR system to study the impact on dynamic sensing,” International Journal of Smart Sensing and Intelligent Systems, vol. 8, no. 3, pp. 1576–1600, 2015.View at: Google Scholar
L. Debnath, Wavelets and Signal Processing, Springer, Berlin, Germany, 2012.
R. W. Boyd, Nonlinear Optics, Academic Press, Amsterdam, The Netherlands, 3rd edition, 2008.View at: MathSciNet
G. Agrawal, Nonlinear Fibre Optics, Academic Press, 2nd edition, 1995.
S. Xie, M. Pang, X. Bao, and L. Chen, “Polarization dependence of Brillouin linewidth and peak frequency due to fiber inhomogeneity in single mode fiber and its impact on distributed fiber Brillouin sensing,” Optics Express, vol. 20, no. 6, pp. 6385–6399, 2012.View at: Publisher Site | Google Scholar
D. Alasia, Advanced Trends in Nonlinear Opitcs Applied to Distributed Optical-Fibre Sensors, EPFL, Lausanne, Switzerland, 2006.
H. Murayama, K. Kageyama, A. Shimada, and A. Nishiyama, “Improvement of spatial resolution for strain measurements by analyzing Brillouin gain spectrum,” in 17th International Conference on Optical Fibre Sensors, vol. 5855 of Proceedings of SPIE, Bruges, Belgium, August 2005.View at: Publisher Site | Google Scholar