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

Recognition of centimeter-level length changes using the intensity of probe light in BOTDA

Open Access Open Access

Abstract

We proposed a method based on the Brillouin optical time domain analysis (BOTDA) system to demodulate the length changes of the heating region using the intensity of probe light, under the premise that the temperature in the heating region does not change and the Brillouin frequency shift (BFS) of the heating region is known. This method can realize the recognition of small length under the condition of wide pulse width. The theoretical analysis shows that the different lengths of the heating region will lead to different gains when the spatial resolution corresponding to the pulse width of the pulsed light is larger than the length of the heating region. And through theoretical derivation and simulation, it is concluded that the relationship between the intensity of probe light and the length of the heating region satisfies the exponential growth. Further experiments verify that the exponential growth is approximately linear in the range of small length changes. The length of the heating region can be inversely deduced by using the linear relationship by detecting the intensity of probe light. When the length of sensing fiber is 100 m and the heating region changes from 2 cm to 8 cm, 1 cm heating region changes can be identified.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Brillouin optical time domain analysis (BOTDA) [1,2] has advantages of high received signal strength, high spatial resolution, long sensing-distance and high-precision measurement. It is widely used to monitor the health of large-scale buildings, such as bridges, oil pipelines, tunnels, etc. in real time. BOTDA is based on the stimulated Brillouin scattering (SBS) effect in optical fiber, in which the pump light and probe light propagate in opposite directions, and the linear relationship between BFS and temperature and strain is used to achieve distributed sensing [35]. In the BOTDA system [69], spatial resolution, measurement accuracy, and sensing distance have always been the research hotspots of many researchers. These indicators interact with each other, and improving one indicator always comes at the expense of one or more other indicators [1025]. The spatial resolution of the traditional BOTDA system depends on the pulse width of the pulsed light, and a higher spatial resolution can be achieved as the pulse width of the pulsed light decreases. Due to the influence of the phonon lifetime, the pulse width of the pulsed light cannot be less than 10 ns, which makes it difficult for the spatial resolution of the system to reach below 1 m. In the case of narrow pulses, the Brillouin gain spectrum (BGS) is broadened, which greatly affects the accuracy of the BFS. The above situation shows that high spatial resolution and high precision measurement cannot be achieved simultaneously in the BOTDA system. Researchers have proposed many methods to improve the spatial resolution. In 2004, Seok-Beom Cho et al. [26] proposed to use the double-pulse technique to double the resolution of two adjacent events at a sensing distance of 2.1 km without reducing the dynamic range. In 2008, Che-Hien Li et al. [27] achieved a spatial resolution of 2 cm using pulse pre-pumping technology, and it is expected to achieve a sensing distance of 100 m. Anthon W. Brown et al. [28] proposed a Brillouin loss-type BOTDA based on dark pulse pumping, and achieved a spatial resolution of 50 mm at a sensing distance of 100 m. The differential pulse technology proposed by Xiaoyi Bao et al. [29,30] used the time trace obtained by two pulsed lights with different pulse widths to perform subtraction, which can significantly improve the spatial resolution of the system, and achieve a spatial resolution of 18 cm on a 1 km sensing fiber. In 2021, the rising edge demodulation algorithm [31] proposed by Yongkang Dong et al. combined with the differential pulse technology experimentally achieved a spatial resolution of 8 mm on a polarization-maintaining fiber with a sensing distance of 11 m. Most of the above technologies shorten the pulse width of pulsed light by various methods, break through the limitation of phonon lifetime, and achieve high spatial resolution. In this paper, a method is proposed to achieve high spatial resolution without reducing the pulse width of pulsed light. The length of heating region is characterized by the intensity of probe light.

We proposed a method for demodulating the length of the heating region using the intensity of probe light. The method realizes the demodulation of the length of the heating region under the premise of the known BFS of the heating region, and focuses on improving the spatial resolution of the BOTDA system. The traditional BOTDA system determines the length of the heating region by the number of sampling points of the heating region on the time trace. The method proposed in this paper realizes the identification of the length of the heating region by using the amplitude of the time trace of the heating region obtained by the wide pulse whose spatial resolution is larger than the length of the heating region. Using this method, it is easy to monitor the length change of the heating region on the order of cm. This paper theoretically analyzes the time trace obtained in two cases when the spatial resolution corresponding to the pulse width of the pulsed light is smaller than the length of the heating region and larger than the length of the heating region. The reason why the latter does not reflect the heating region in the three-dimensional BGS is explained in detail. The theoretical analysis process of this study can provide a new idea for improving the spatial resolution of the BOTDA system.

2. Principle

2.1 Analysis

This section firstly analyzes the limitations of wide pulsed light in distributed detection, and further theoretically proves the feasibility of using the intensity of probe light to demodulate the length of the heating region. In the BOTDA system, the pump light and the probe light traveling in opposite directions will produce SBS effect in the fiber. In the gain-type BOTDA system, the energy of the pump light is transferred to the probe light. When the frequency difference is exactly the BFS of the fiber, the SBS effect is the most obvious, the gain obtained by the probe light is the largest, and the amplitude of time trace is increased. The sensing is realized by using the linear relationship between BFS and temperature and strain. In the BOTDA system, the frequency of probe light is swept to obtain the three-dimensional BGS along the fiber. Based on the acquired gain spectrum, the position of the heating region, the length of the heating region, and the BFS of the heating region are determined. The above is the practice when the spatial resolution of the pulse width of the pulsed light is smaller than the length of the heating region, and it is also the demodulation principle of the heating region of the traditional BOTDA system. When the spatial resolution of the pulse width is larger than the length of the heating region, the three-dimensional BGS obtained by sweeping the frequency of probe light is shown in Fig. 1.

 figure: Fig. 1.

Fig. 1. Top view of the three-dimensional BGS, the pulse width is 50 ns, the length of heating region is 50 cm.

Download Full Size | PDF

Figure 1 shows a top view of the three-dimensional BGS. In the figure, the 100 m-long sensing fiber has a 50 cm heating region near the end, and the heating temperature is 70°C, room temperature is 25°C. Pulsed light with a pulse width of 50 ns is used for detection with a spatial resolution of 5 m. The BFS of the fiber obtained in the figure is 10.627 GHz. The heating region is close to the end of the fiber, and the position of the heating region can be determined from the figure, but the length of the heating region and the BFS of the heating region cannot be obtained intuitively. That is often said that the use of wide pulses cannot detect the tiny lengths, and the spatial resolution of the pulse width cannot meet the requirements of detecting the length of the heating region. There are two reasons why we cannot see the heating region above the three-dimensional BGS. First, because the pulsed light with a wide pulse width is used for detection in the experiment, any heating region is always approximately equal to the length of the spatial resolution corresponding to the pulse width of the pulsed light. Second, the length of the heating region is much shorter than the spatial resolution corresponding to the width of the pulsed light, which leads to a small gain obtained by the probe light when the frequency of the probe light sweeps to the BFS of the heating region. The above two reasons lead to the inability to determine the length of the heating region and the BFS of the heating region on the three-dimensional BGS.

The BGS of the room temperature region and the heating region are respectively taken on the 100 m long fiber, as shown in Fig. 2. The BFS of the sensing fiber is 10.627 GHz, and the amplitude at the peak of the BGS is 132.4 mV. Figure 2(b) shows the BGS at the heating region of the sensing fiber near the end. From the figure, it can be read that the BFS corresponding to the maximum amplitude is 10.627 GHz. The gain spectrum at this position does not show the gain spectrum we expect to characterize the information of heating region. The reason for this is that the pulse width of the 50 ns pulse light used in the experiment is too wide, and the length of the fiber in the room temperature region covered by the pulse light with the 50 ns pulse width in the duration is much longer than the length of the heating region. As a result, when the scanning frequency of the probe light is the BFS in the room temperature region of the optical fiber, the gain obtained is still very large. Similarly, when the scanning frequency of the probe light is the BFS of the heating region, the gain obtained is very small. Because the length of the heating region is too short, the BFS at the highest amplitude in the gain spectrum is still the BFS corresponding to the room temperature region, and the BGS in the room temperature region is still present. Of course, the difference between the two figures is still obvious. In Fig. 2(a), the amplitude at the peak of the BGS in the room temperature region is greater than that at the peak of the BGS at the position of the heating region in Fig. 2(b). At the same time, there is a small bump on the right side of the gain spectrum in Fig. 2(b), which represents the BFS in the heating region. From the Fig. 2(b), the BFS in the heating region can be read as 10.674 GHz. The above situation shows that when the spatial resolution of the pulse width of the pulsed light is larger than the length of the heating region, we cannot determine the length of the heating region and the BFS of the heating region through the three-dimensional BGS. However, compared with the BGS in the room temperature region of the fiber, there will be a change in amplitude. This provides us with an idea to identify the length of the heating region by finding the relationship between the detected intensity of probe light and the length of the heating region in the case of a wide pulse width.

 figure: Fig. 2.

Fig. 2. (a) BGS of the fiber at room temperature; (b) BGS of the fiber at the heating region.

Download Full Size | PDF

2.2 Theory

When the frequency difference between the probe light and the pump light is the BFS in the heating region. With the increase of the length of the heating region, the time trace under the conditions that the spatial resolution of the pulse width is smaller than the length of the heating region and the spatial resolution of the pulse width is larger than the length of the heating region as shown in Fig. 3. Figure 3(a) is the measured time trace with the increase of the length of the heating region when the spatial resolution of the pulse width of the pulsed light is smaller than the length of the heating region. Without considering the pump consumption, the measured time trace of the heating region only broadening with the increase of the length of the heating region, and the intensity does not change. Figure 3(b), the spatial resolution of the pulse width of the pulsed light is much larger than the length of the heating region. As the length of the heating region increases, in addition to the broadening phenomenon in the time trace, the intensity of probe light of the heating region continues to increase. This is because in the duration of a pulse, the longer the heating region, the longer the gain time obtained by the probe light in this period, thus showing an increase in intensity in time trace.

 figure: Fig. 3.

Fig. 3. Schematic diagram of the intensity of probe light demodulating the length of the heating region. (a) The spatial resolution corresponding to the pulse width of the pulsed light is smaller than the length of the heating region; (b) the spatial resolution corresponding to the pulse width of the pulsed light is larger than the length of the heating region.

Download Full Size | PDF

The SBS process in an optical fiber is described by the following coupling equation:

$$\frac{{d{I_S}}}{{dz}} ={-} {g_B}{I_P}{I_S} + \alpha {I_S},$$
$$\frac{{d{I_P}}}{{dz}} ={-} {g_B}{I_P}{I_S} - \alpha {I_P}.$$
Here ${I_S}$ and ${I_P}$ represent the intensities of the probe light and the pump light, respectively, ${g_B}$ represents the Brillouin gain coefficients, $\alpha $ represents the loss coefficients of the probe light and the pump light in the fiber. Without considering the pump consumption, the distribution of the intensity of pump light can be expressed using Eq. (3), as follows:
$${I_P}(z) = {I_P}(0)\exp ( - \alpha z),$$
where ${I_P}(z)$ is the intensity of pump light at any position of the fiber, and ${I_P}(0)$ is the intensity of pump light at the beginning of the fiber. Substituting Eq. (3) into Eq. (1) for integration, the distribution of the intensity of probe light along the fiber can be expressed using Eq. (4), as follows:
$${I_S}(z) = {I_S}(L)\exp ( - \frac{{{g_B}{I_P}(0)\exp ( - \alpha L)}}{\alpha } - \alpha L)\exp (\frac{{{g_B}{I_P}(0)\exp ( - \alpha z)}}{\alpha } + \alpha z),$$
${I_S}(L)$ is the intensity of probe light of $z = L$,the intensity of probe light of $z = 0$ can be expressed using Eq. (5), as follows:
$${I_s}(0) = {I_s}(L)\exp ({g_B}{I_P}(0){L_{eff}} - \alpha L),$$
where ${L_{eff}} = [1 - \exp ( - \alpha L)]/\alpha $ is the effective length of the fiber, using the relationship between optical power and the intensity of light: ${P_P} = {I_P}{A_{eff}}$, ${P_S} = {I_S}{A_{eff}}$, ${A_{eff}}$ represents the effective area of the fiber. Substituting into Eq. (3), Eq. (4) and Eq. (5) to get
$${P_P}(z) = {P_P}(0)\exp ( - \alpha z),$$
$${P_S}(z) = {P_S}(L)\exp ( - \frac{{{g_B}{P_P}(0)\exp ( - \alpha L)}}{{\alpha {A_{eff}}}} - \alpha L)\exp (\frac{{{g_B}{P_P}(0)\exp ( - \alpha z)}}{{\alpha {A_{eff}}}} + \alpha z),$$
$${P_S}(0) = {P_S}(L)\exp ({g_B}{P_P}(0){L_{eff}}/{A_{eff}} - \alpha L).$$

Transforming Eq. (8) as follows:

$${P_S}(0) = \frac{{{P_S}(L)}}{{\exp (\alpha L)}}\exp (\frac{{{g_B}{p_P}(0)}}{{{A_{eff}}}}{L_{eff}}),$$
when the total length of the fiber remains unchanged, the consumption of the pump light is ignored, and the frequency shift of the probe light is constant. ${P_S}(L)$, $\exp (\alpha L)$, ${g_B}$, ${P_P}(0)$ and ${A_{eff}}$ are all constants. Equation (8) can be transformed as follows:
$${P_S}(0) = A\exp (B{L_{eff}}),$$
where $A = {P_S}(L)/\exp (\alpha L)$, $B = {g_B}{P_P}(0)/{A_{eff}}$. Therefore, it is deduced that the intensity of probe light increasing exponentially with the increase of the length of heating region when the resolution of pulse width is larger than the length of heating region. According to the Taylor expansion formula of the exponential function, when the variation of the length of the heating region is very small, it can be approximately considered that the relationship of linear growth is satisfied near the change position.

2.3 Simulation result

Figure 4 shows the results of the simulation experiment. The length of the sensing fiber is 100 m, the pulse width of the pulsed light is 200 ns, and the corresponding spatial resolution is 20 m. There is a heating region near the middle part of the entire sensing fiber. The length of heating region varies from 2 m to 10 m in steps of 1 m. With the increase of the length of the heating region, the measured time trace of probe light under the length of each heating region are shown in Fig. 4(a). The horizontal axis represents the sensing distance, and the vertical axis represents the intensity of probe light. It can be clearly seen from the curve that the time trace of the heating region broadened as the length of the heating region increase, and the intensity of the probe light also increased with the increase of the length of the heating region. It is fully explained that the gain provided to the probe light is different for different lengths of heating region. There is a positive correlation between them. Each point in Fig. 4(b) is the average of the time trace under each simulated heating region in Fig. 4(a). The horizontal axis represents the length of the heating region corresponding to each time trace of probe light, and the vertical axis represents the corresponding intensity of probe light at the length of heating region. An exponential fitting is performed on the scatter plot, as shown by the red curve in the Fig. 4(b). The fitting expression is shown in Fig. 4(b), which is consistent with that described in Eq. (10). The simulation results show that in the case of wide pulsed light, the intensity of the probe light and the length of the heating region satisfy an exponential growth.

 figure: Fig. 4.

Fig. 4. Simulation results show that the spatial resolution corresponding to the pulse width of the pulsed light is larger than the length of the heating region. (a) The change of the intensity of probe light corresponding to the heating region changing from 2 m to 10 m in steps of 1 m. (b) Exponential fitting of the intensity of probe light corresponding to the length of each heating region.

Download Full Size | PDF

At the same time, the case where the spatial resolution of the pulse width of the pulsed light is smaller than the length of the heating region is also explored. Figure 5 shows the simulation results.

 figure: Fig. 5.

Fig. 5. Simulation results show that the spatial resolution corresponding to the pulse width of the pulsed light is smaller than the length of the heating region. The change of the intensity of probe light corresponding to the heating region changing from 2 m to 10 m in steps of 1 m.

Download Full Size | PDF

The length of the sensing fiber is still 100 m, the pulse width of the pulsed light is 20 ns, the corresponding spatial resolution is 2 m, and the length of the heating region changes from 2 m to 10 m in steps of 1 m. The horizontal axis represents the sensing distance, and the vertical axis represents the intensity of probe light. Ignoring the consumption of pump light, it can be seen from the Fig. 5 that with the increase of the length of the heating region, when the spatial resolution of the pulse width of the pulsed light is smaller than the length of the heating region, the gain time obtained by the probe light at each position in the heating region is the same. It shows that the intensity of probe light does not increase with the increase of the length of the heating region, only with the increase of the length of the heating region, the time trace of probe light of the heating region is broadened.

The length of the heating region was further reduced, and the heating region was changed from 10 cm to 50 cm in steps of 10 cm. As shown in Fig. 6(a), it is the time trace corresponding to each heating region, and Fig. 6(b) shows the intensity of the probe light and the length of the heating region satisfy a linear increase, the fitting expression is also in the figure.

 figure: Fig. 6.

Fig. 6. (a) The change of the intensity of probe light corresponding to the heating region changing from 10 cm to 50 cm in steps of 10 cm. (b) Linear fitting of the intensity of probe light corresponding to the length of each heating region.

Download Full Size | PDF

3. Experimental results

Figure 7 is an experimental setup of the BOTDA system. The light emitted by a narrow linewidth laser (DFB-LD) with a wavelength of 1550 nm and the linewidth of 3 kHz is divided into two beams by a polarization-maintaining coupler (PM coupler). The upper branch is the probe light and the lower branch is the pump light. The upper branch uses an electro-optical intensity modulator (EOM) and a bias control board for double-sideband carrier suppression, and uses a microwave signal generator (MWG) to sweep the frequency of probe light. The light from the EOM is divided into two paths by a 99:1 coupler, and 1% of the light is fed back to the bias control board to stabilize the working point of the EOM, 99% of the light is amplified by an erbium-doped fiber amplifier (EDFA), and then enters the sensing fiber through a polarization scrambler (PS) and an isolator. The light in the lower branch first passes through the variable optical attenuator (VOA) and enters the semiconductor optical amplifier (SOA) for pulse modulation, then enters the pulsed optical amplifier for amplification, and enters the sensing fiber through the optical circulator (OC). The Brillouin signal is filtered by a tunable fiber Bragg grating (FBG) to filter out the low-order sidebands of the probe light, and the Brillouin signal is detected by a 300 M photodetector (PD), and the signal is collected and stored by an industrial personal computer (IPC). In the experiment, a heating plate was used to heat the fiber near the end of the 100 m-long single-mode fiber, and the heating temperature was kept at 72°C. The BFS of the fiber at room temperature 25°C is 10.627 GHz, and the BFS in the heating region is 10.674 GHz. The length of the heating region is varied from 10 cm to 50 cm in steps of 5 cm. The length of the heating region increases away from the direction of the pump light. In the experiment, pulsed light with a pulse width of 50 ns is used for detection. The repetition frequency of pulse light is 200 kHz, and the input power of the pump light and the probe light are 1.2 mW and 0.7 mW, respectively. The frequency shift of the probe light is 10.674 GHz. The sampling rate of the acquisition card is 2 GS/s, the number of average times is 5000, and the number of sampling points is 3000.

 figure: Fig. 7.

Fig. 7. BOTDA experimental setup.

Download Full Size | PDF

As shown in Fig. 8, Fig. 8(a) is the time trace of the probe light obtained under different lengths of heating region, the horizontal axis is the sensing distance, and the vertical axis is the amplitude of the probe light detected by the photodetector. The position of heating region is where the time trace begins to rise. Under the condition that the parameters of the entire BOTDA system are kept unchanged, such as the injected optical power, the frequency shift of the probe light, the pulse width of the pump light, the repetition frequency of the pump light, and the sampling rate of the acquisition card, the average times, and the number of sampling points are all constant value. Only changing the length of the heating region, the amplitude of time trace at the position of the heating region increased with the increase of the length of the heating region, which is consistent with the experimental theory and simulation results. A change in length of 5 cm in the heating region can be recognized. Each point in Fig. 8(b) is the value obtained by accumulating and averaging over the entire pulse duration according to the time trace of different lengths of heating region in Fig. 8(a). The horizontal axis represents each length of the heating region, and the vertical axis represents the corresponding amplitude at the length of the heating region. The blue line in the Fig. 8(b) is the result of exponentially fitting, the fitting expression is shown in Fig. 8(b), which is consistent with the form of Eq. (10), and the curve fitting is close to 0.99. It is proved that when the resolution of pulse width is larger than the length of the heating region, the gain obtained by the probe light satisfies the exponential growth with the increase of the length of the heating region. The difference between Fig. 8(c) and Fig. 8(b) is that each point in Fig. 8(c) is linearly fitted, the fitting expression is shown in Fig. 8(c), and the curve fitting is also close to 0.99. In the case of a small range of length variation, the curve conforming to exponential growth is close to linearly growth. This provides convenience for the subsequent demodulation of the length of the heating region by detecting the intensity of probe light.

 figure: Fig. 8.

Fig. 8. Experimental results. (a) The change of the amplitude of probe light corresponding to the change of the length of heating region from 10 cm to 50 cm in steps of 5 cm; (b) the exponential fitting of the amplitude of probe light corresponding to the length of each heating region; (c) the linear fitting of the amplitude of probe light corresponding to the length of each heating region.

Download Full Size | PDF

The length of the heating region was further reduced as shown in Fig. 9, and the length of the heating region was changed from 2 cm to 8 cm in steps of 1 cm. Parameter settings of the system remain unchanged. Figure 9(a) is the time trace of the probe light at each length of heating region, the horizontal axis is the sensing distance, and the vertical axis is the amplitude of the probe light detected by the photodetector. Each point in Fig. 9(b) is the value obtained by accumulating and averaging over the entire pulse duration according to the time trace of different lengths of heating region in Fig. 9(a). The horizontal axis represents each length of the heating region, and the vertical axis represents the corresponding amplitude at the length of the heating region. The red line is the result of linearly fitting, the fitting expression is shown in Fig. 9(b), the curve fitting is close to 0.99. It shows that when the length of the heating region with a small variation range, the amplitude of the time trace in the heating region and the length of the heating region satisfy a linear relationship.

 figure: Fig. 9.

Fig. 9. Experimental results. (a) The change of the amplitude of probe light corresponding to the change of the length of heating region from 2 cm to 8 cm in steps of 1 cm; (b) the linear fitting of the amplitude of probe light corresponding to the length of each heating region.

Download Full Size | PDF

As shown in Table 1, pulse light with pulse widths of 50 ns, 70 ns, and 90 ns was used to detect the heating region from 20 cm to 50 cm in steps of 5 cm. The power of probe light is the constant of 0.705 mW. The measurement was repeated 3 times under each pulse width. The experimental results show that with different pulse widths, the fitting slopes obtained are different. The pulse width remains unchanged, and the test is repeated three times, and the slope between the amplitude of probe light and the length of the heating region is basically unchanged. In Table 2, the pulse width of the pulsed light is 50 ns, and the length of the heating region changes from 10 cm to 50 cm in steps of 5 cm. Three sets of experimental data were obtained when the power of probe light was 0.738 mW, 0.977 mW, and 1.172 mW, respectively. Under the condition that the pulse width of the pulsed light remains unchanged, the slope values obtained under different power of probe light are different, but the three sets of slope values obtained under the same power of probe light basically remain unchanged. This also proves that when the pulse width of the pulsed light is constant, the power of the probe light is constant, and the length of the fiber is constant, the amplitude of the probe light and the length of the heating region always satisfy a specific growth rate. This provides a basis for inverting the length of the heating region by using the intensity of the probe light through a linear fitting relationship.

Tables Icon

Table 1. Slopes for three experiments at three pulse widths

Tables Icon

Table 2. Slopes for three experiments at three power of probe light

4. Conclusion

In this paper, we displayed when the spatial resolution of the pulse width is larger than the length of the heating region, the reason why the length of the heating region and the BFS of the heating region cannot be distinguished on the three-dimensional BGS. The first one is that when we use pulsed light with a wide pulse width for detection, and the length of the heating region obtained on the three-dimensional gain spectrum is always approximately equal to the length of the spatial resolution corresponding to the pulse width of the pulsed light. Second, the length of the heating region is too short to provide sufficiently gain for the probe light, and the intensity of the probe light is too weak, so the size of the BFS in the heating region cannot be distinguished. Then a method to demodulate the length of the heating region by the intensity of the probe light is proposed. This method utilizes a pulsed light whose spatial resolution of the pulse width is much larger than the length of the heating region under the premise of the known BFS of the heating region. The intensity of the probe light increased with the length of the heating region. The relationship between the intensity of the probe light and the length of the heating region satisfies exponential growth through simulation, which is verified in experiments. At the same time, with a small range of length variation, the exponential growth can be approximately considered as a linear growth, which greatly simplifies the calculation process of using the intensity of probe light to infer the length of the heating region. The length change of the heating region of 1 cm was successfully discriminated through experiments. The relationship between the slope of the linear fitting and the pulse width of the pulsed light and the power of the probe light was further studied. The experimental results show that the slope changes with the pulse width and the power of probe light. Under the same pulse width or the same power of probe light, the slope of repeated experiments remains basically unchanged. The reliability of demodulating the length of the heating region according to the intensity of the probe light with a satisfied linear relationship is demonstrated. The demodulation of the length of heating region based on the intensity of probe light proposed in this paper provides a new idea for researchers to improve the spatial resolution of the BOTDA system.

Funding

National Natural Science Foundation of China (61527819, 61875146, 62075151, 62075153); Patent Promotion and Exploitation Program of Shanxi Province (20200734); Transformation of Scientific and Technological Achievements Programs of Higher Education Institutions in Shanxi; Fund for Shanxi “1331 Project” Key Innovative Research Team; Shanxi-Zheda Institute of Advanced Materials and Chemical Engineering.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. T. Horiguchi and M. Tateda, “BOTDA-Nondestructive measurement of single-mode optical fiber attenuation characteristics using Brillouin interaction: Theory,” J. Lightwave Technol. 7(8), 1170–1176 (1989). [CrossRef]  

2. D. Culverhouse, F. Farahi, C. N. Pannell, and D. A. Jackson, “Potential of stimulated Brillouin scattering as sensing mechanism for distributed temperature sensors,” Electron. Lett. 25(14), 913–915 (1989). [CrossRef]  

3. D. Zhou, Y. Dong, B. Wang, C. Pang, D. Ba, H. Zhang, Z. Lu, H. Li, and X. Bao, “Single-shot BOTDA based on an optical chirp chain probe wave for distributed ultrafast measurement,” Light: Sci. Appl. 7(1), 32–43 (2018). [CrossRef]  

4. B. Wang, B. Fan, D. Zhou, C. Pang, Y. Li, D. Ba, and Y. Dong, “High-performance optical chirp chain BOTDA by using a pattern recognition algorithm and the differential pulse-width pair technique,” Photonics Res. 7(6), 652–658 (2019). [CrossRef]  

5. Y. Chang, H. Wu, C. Zhao, L. Shen, S. Fu, and M. Tang, “Distributed Brillouin frequency shift extraction via a convolutional neural network,” Photonics Res. 8(5), 690–697 (2020). [CrossRef]  

6. P. Lu, N. Lalam, M. Badar, B. Liu, B. T. Chorpening, M. P. Buric, and P. R. Ohodnicki, “Distributed optical fiber sensing: Review and perspective,” Appl. Phys. Rev. 6(4), 041302 (2019). [CrossRef]  

7. A. Denisov, M. A. Soto, and L. Thévenaz, “Going beyond 1000000 resolved points in a Brillouin distributed fiber sensor: theoretical analysis and experimental demonstration,” Light: Sci. Appl. 5(5), e16074 (2016). [CrossRef]  

8. X. Bao and L. Chen, “Recent Progress in Brillouin Scattering Based Fiber Sensors,” Sensors 11(4), 4152–4187 (2011). [CrossRef]  

9. T. Horiguchi, K. Shimizu, T. Kurashima, M. Tateda, and Y. Koyamada, “Development of a Distributed Sensing Technique Using Brillouin Scattering,” J. Lightwave Technol. 13(7), 1296–1302 (1995). [CrossRef]  

10. K. Luo, B. Wang, N. Guo, K. Yu, C. Yu, and C. Lu, “Enhancing SNR by Anisotropic Diffusion for Brillouin Distributed Optical Fiber Sensors,” J. Lightwave Technol. 38(20), 5844–5852 (2020). [CrossRef]  

11. C. Feng, X. Lu, S. Preussler, and T. Schneider, “Gain Spectrum Engineering in Distributed Brillouin Fiber Sensors,” J. Lightwave Technol. 37(20), 5231–5237 (2019). [CrossRef]  

12. G. Ryu, G. T. Kim, K. Y. Song, S. B. Lee, and K. Lee, “Linearly Configured Brillouin Optical Correlation Domain Analysis System Incorporating Time-Domain Data Processing,” J. Lightwave Technol. 37(18), 4728–4733 (2019). [CrossRef]  

13. Y. Wang, M. Zhang, J. Zhang, L. Qiao, T. Wang, Q. Zhang, L. Zhao, and Y. Wang, “Millimeter-Level-Spatial-Resolution Brillouin Optical Correlation-Domain Analysis Based on Broadband Chaotic Laser,” J. Lightwave Technol. 37(15), 3706–3712 (2019). [CrossRef]  

14. J. Liu, T. Wang, Q. Zhang, H. Chen, J. Zhang, L. Qiao, S. Gao, and M. Zhang, “Rapid noise removal based dual adversarial network for the Brillouin optical time domain analyzer,” Opt. Express 29(21), 34002–34014 (2021). [CrossRef]  

15. T. Wang, J. Zhao, Q. Zhang, J. Liu, H. Chen, J. Zhang, and M. Zhang, “Extracting Brillouin frequency shift accurately based on particle swarm optimization and a cross-correlation method,” Appl. Opt. 61(2), 446–454 (2022). [CrossRef]  

16. Y. Mei, X. Xu, L. Luo, and K. Soga, “Reconstruction of Distributed Strain Profile Using a Weighted Spectrum Decomposition Algorithm for Brillouin Scattering Based Fiber Optic Sensor,” J. Lightwave Technol. 38(22), 6385–6392 (2020). [CrossRef]  

17. X. Bao, D. J. Webb, and D. A. Jackson, “32-km distributed temperature sensor based on Brillouin loss in an optical fiber,” Opt. Lett. 18(18), 1561–1563 (1993). [CrossRef]  

18. S. M. Maughan, H. H. Kee, and T. P. Newson, “57-km single-ended spontaneous Brillouin-based distributed fiber temperature sensor using microwave coherent detection,” Opt. Lett. 26(6), 331–333 (2001). [CrossRef]  

19. D. Lida and F. Ito, “Cost-effective bandwidth-reduced Brillouin optical time domain reflectometry using a reference Brillouin scattering beam,” Appl. Opt. 48(22), 4302–4309 (2009). [CrossRef]  

20. Y. T. Cho, M. Alahbabi, M. J. Gunning, and T. P. Newson, “50-km single-ended spontaneous-Brillouin-based distributed-temperature sensor exploiting pulsed Raman amplification,” Opt. Lett. 28(18), 1651–1653 (2003). [CrossRef]  

21. Y. T. Cho, M. N. Alahbabi, M. J. Gunning, and T. P. Newson, “Enhanced performance of long range Brillouin intensity based temperature sensors using remote Raman amplification,” Meas. Sci. Technol. 15(8), 1548–1552 (2004). [CrossRef]  

22. F. Rodríguez-Barrios, S. Martín-López, A. Carrasco-Sanz, P. Corredera, J. D. Ania-Castañón, L. Thévenaz, and M. González-Herráez, “Distributed Brillouin Fiber Sensor Assisted by First-Order Raman Amplification,” J. Lightwave Technol. 28(15), 2162–2172 (2010). [CrossRef]  

23. G. Ryu, G. T. Kim, K. Y. Song, S. B. Lee, and K. Lee, “50 km-Range Brillouin Optical Correlation Domain Analysis With First-Order Backward Distributed Raman Amplification,” J. Lightwave Technol. 38(18), 5199–5204 (2020). [CrossRef]  

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

25. M. A. Soto, G. Bolognini, F. Di Pasquale, and L. Thévenaz, “Simplex-coded BOTDA fiber sensor with 1 m spatial resolution over a 50 km range,” Opt. Lett. 35(2), 259–261 (2010). [CrossRef]  

26. S. B. Cho, J. J. Lee, and I. B. Kwon, “Strain event detection using a double-pulse technique of a Brillouin scattering-based distributed optical fiber sensor,” Opt. Express 12(18), 4339–4346 (2004). [CrossRef]  

27. C. H. Li, K. Nishiguchi, and M. Miyatake, “PPP-BOTDA method to achieve 2 cm spatial resolution in Brillouin distributed measuring technique,” IEICE Tech. Rep. 32(108), 55–60 (2008).

28. A. W. Brown, B. G. Colpitts, and K. Brown, “Distributed sensor based on dark-pulse Brillouin scattering,” IEEE Photon. Technol. Lett. 17(7), 1501–1503 (2005). [CrossRef]  

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

30. Y. Dong, X. Bao, and W. Li, “Differential Brillouin gain for improving the temperature accuracy and spatial resolution in a long-distance distributed fiber sensor,” Appl. Opt. 48(22), 4297–4301 (2009). [CrossRef]  

31. D. Zhou, D. Ba, B. Wang, L. Qiu, W. L. J. Hasi, and Y. Dong, “Millimeter-level recognition capability of BOTDA based on a transient pump pulse and algorithm enhancement,” Opt. Lett. 46(14), 3440–3443 (2021). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1.
Fig. 1. Top view of the three-dimensional BGS, the pulse width is 50 ns, the length of heating region is 50 cm.
Fig. 2.
Fig. 2. (a) BGS of the fiber at room temperature; (b) BGS of the fiber at the heating region.
Fig. 3.
Fig. 3. Schematic diagram of the intensity of probe light demodulating the length of the heating region. (a) The spatial resolution corresponding to the pulse width of the pulsed light is smaller than the length of the heating region; (b) the spatial resolution corresponding to the pulse width of the pulsed light is larger than the length of the heating region.
Fig. 4.
Fig. 4. Simulation results show that the spatial resolution corresponding to the pulse width of the pulsed light is larger than the length of the heating region. (a) The change of the intensity of probe light corresponding to the heating region changing from 2 m to 10 m in steps of 1 m. (b) Exponential fitting of the intensity of probe light corresponding to the length of each heating region.
Fig. 5.
Fig. 5. Simulation results show that the spatial resolution corresponding to the pulse width of the pulsed light is smaller than the length of the heating region. The change of the intensity of probe light corresponding to the heating region changing from 2 m to 10 m in steps of 1 m.
Fig. 6.
Fig. 6. (a) The change of the intensity of probe light corresponding to the heating region changing from 10 cm to 50 cm in steps of 10 cm. (b) Linear fitting of the intensity of probe light corresponding to the length of each heating region.
Fig. 7.
Fig. 7. BOTDA experimental setup.
Fig. 8.
Fig. 8. Experimental results. (a) The change of the amplitude of probe light corresponding to the change of the length of heating region from 10 cm to 50 cm in steps of 5 cm; (b) the exponential fitting of the amplitude of probe light corresponding to the length of each heating region; (c) the linear fitting of the amplitude of probe light corresponding to the length of each heating region.
Fig. 9.
Fig. 9. Experimental results. (a) The change of the amplitude of probe light corresponding to the change of the length of heating region from 2 cm to 8 cm in steps of 1 cm; (b) the linear fitting of the amplitude of probe light corresponding to the length of each heating region.

Tables (2)

Tables Icon

Table 1. Slopes for three experiments at three pulse widths

Tables Icon

Table 2. Slopes for three experiments at three power of probe light

Equations (10)

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

d I S d z = g B I P I S + α I S ,
d I P d z = g B I P I S α I P .
I P ( z ) = I P ( 0 ) exp ( α z ) ,
I S ( z ) = I S ( L ) exp ( g B I P ( 0 ) exp ( α L ) α α L ) exp ( g B I P ( 0 ) exp ( α z ) α + α z ) ,
I s ( 0 ) = I s ( L ) exp ( g B I P ( 0 ) L e f f α L ) ,
P P ( z ) = P P ( 0 ) exp ( α z ) ,
P S ( z ) = P S ( L ) exp ( g B P P ( 0 ) exp ( α L ) α A e f f α L ) exp ( g B P P ( 0 ) exp ( α z ) α A e f f + α z ) ,
P S ( 0 ) = P S ( L ) exp ( g B P P ( 0 ) L e f f / A e f f α L ) .
P S ( 0 ) = P S ( L ) exp ( α L ) exp ( g B p P ( 0 ) A e f f L e f f ) ,
P S ( 0 ) = A exp ( B L e f f ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.