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

SNR enhancement of quasi-distributed weak acoustic signal detection by elastomers and MMF integrated Φ-OTDR

Open Access Open Access

Abstract

We have proposed and demonstrated a weak acoustic signal detection technology based on phase-sensitive optical time-domain reflectometry (Φ-OTDR). Non-contact acoustic signals transmitting through air gap between the sound source and the receiver are difficult to detect due to fast attenuation. In order to improve the detection ability of non-contact weak acoustic signals, we demonstrate that multi-mode fiber (MMF) is a better solution than single-mode fiber (SMF) benefiting from its larger core and higher Rayleigh backscattering (RBS) capture coefficient. The frequency signal-to-noise ratio (SNR) has been enhanced by 9.26 dB. Then, with the help of 3D printing technology, elastomers have been designed to further enhance the detection ability due to the high-sensitive response to acoustic signals. Compared with the previous reported “I” type elastomer, the location and frequency SNR enhancement caused by our new proposed “n” type elastomer are 8.39 dB and 11.02 dB in SMF based system. The values are further improved to 10.51 dB and 13.38 dB in MMF and “n” type elastomer integrated system. And a phase-pressure sensitivity of -94.62 dB re rad/µPa has been achieved at 2.5 kHz. This non-contact weak acoustic signal detection technique has great application potential in the quasi-distributed partial discharge (PD) detection of smart grid.

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

1. Introduction

Weak acoustic detection is urgently needed in many fields, such as intrusion detection [1], oil exploration [2], pipeline security [3], marine optics [4] and partial discharge (PD) in smart grid [57]. At present, a number of points sensing methods have been reported for weak acoustic signal detection. In 2018, G. M. Ma, et al. proposed a PD detection system based on phase-shifted fiber Bragg gratings to overcome the low sensitivity of power transformer PD detection [6]. In 2019, Quan. Y, et al. proposed an auxiliary Mach-Zehnder interferometer method to compensate for the influences of laser-frequency-drift in order to raise the recognition sensitivity of weak acoustic signals [8]. However, such point sensing methods are limited in the application where the remote sensing or locating is needed [9].

On the other hand, distributed optical fiber sensing technology based on Φ-OTDR with merits of remote sensing and high sensitivity has been explored in locating and monitoring acoustic signals [10,11]. However, the SNR of acoustic signal with standard SMF is fundamentally limited due to the relatively weak RBS light [12] and inherent signal fading [13]. Recently, specialty fibers integrated Φ-OTDR have attracted much attention in distributed acoustic sensing (DAS) [14]. With the assist of the specialty fibers, either the SNR of detected signal can be enhanced or the fading problem can be mitigated. In general, these specialty fibers can be classified into two catalogs by the process of fiber fabrication: the fiber that is directly fabricated from the fiber drawing process, and the fiber that requires an additional post-processing. For the directly fabricated fiber, it mainly includes multi-core fiber (MCF) [15], few-mode fiber (FMF) [16], and MMF [17]. In 2020, Y. Mao, et al. achieved a simultaneous distributed acoustic and temperature sensing based on a MMF based distributed optical fiber sensing with a 12.98-dB SNR and temperature measurements of an ±1 °C accuracy. For the fibers with additional post-processing, it mainly consists of weak fiber Bragg grating (FBG) arrays [18], microstructure optical fibers (MOFs) [19,20], cladding softened fiber [21], UV irradiation fiber [22], etc. For example, D. Liu, et al. proposed the DAS system with MOF, in which the SNR of sensing signal is increased by enhancing the microstructure light scattering. In 2021, J. Yao et al. demonstrated a phase sensitivity improvement technology of DAS by softening sensing fiber cladding. Although the aforementioned specialty fibers can improve the performance of DAS, most of fibers might not be compatible with the standard fiber communication link. High cost and the sophisticated preparation process still have to be taken into consideration. Among the aforementioned specialty fibers, MMF has been more commonly used, especially for short-distance applications [23]. Furthermore, MMF owns its inherent merits, such as high nonlinear threshold, large core for power transmission and relative low cost. Driven by the above advantages, MMF is very attractive to practical application for improving the performance of DAS.

In addition to the specialty fibers, the integration of elastomer in Φ-OTDR for weak acoustic signal detection is of great interest. In 2016, T. Zhang et al. demonstrated a Michelson interferometer-based fiber sensor with a compact and high sensitivity cylindrical elastomer for detecting acoustic emission generated from the PD [24]. However, they are lack of distributed sensing. In 2020, Z. Chen et al. realized sensitivity enhancement of acoustic emission detection during breakdown discharge process by a “I” type elastomer based on 3D printing technology incorporating with a SMF [25]. Most of the above proposed elastomers are actually sensing when they are in contact with the sound source. However, when non-contact monitoring is required, space interval exists between the detection configuration and the sound source. The current elastomer sensing technology is still not sensitive enough for this non-contact DAS application. DAS system with higher acoustic-phase response and larger SNR is still in urgent need for weak acoustic signal detection.

In recent years, high sensitized optical cable has been designed and used for underwater acoustic wave detection [26]. Different from their work, we used the commercial telecommunication MMF for quasi-distributed fiber sensing. As shown in Fig. 1, such MMF with elastomer can enhance the sensitivity and SNR of weak acoustic signal detection, and has great potential in building up the integrated communication and sensing networks. In our experiment, we induce different-intensity acoustic signals between 0.04 Pa to 0.16 Pa. Among the testing experiments, MMF based system enhances 9.26 dB frequency SNR than SMF based system. The “n” type elastomer enhances the location and frequency SNR by 8.39 dB and 11.02 dB in a SMF based system. The values have been further improved to 10.51 dB and 13.38 dB in a MMF based system. Finally, the combine of the “n” type elastomer and the MMF makes the system achieve 19.85 dB location SNR and 32.25 dB frequency SNR. The new proposed sensing link greatly expands the weak acoustic signal detection ability of the Φ-OTDR system, which will benefit the application of quasi-distributed PD detection of smart grid.

 figure: Fig. 1.

Fig. 1. Illustration of sensing head for weak acoustic signal detection by the combined scheme of a MMF and an “n” type elastomer.

Download Full Size | PDF

2. Principles and elastomer design

2.1 Rayleigh backscattering in MMF

In MMF-based DAS system, the RBS signal received into a certain mode is a superposition result of the signals scattered back to this mode from all transmission modes, which is called signal aliasing phenomenon [27]. The RBS received at the input end can be derived by the following Eq. (1), when an optical pulse with a width of Δτ, a power of ${P_0}$ is launched to the MMF. Without launching limitations, the RBS light will couple to all the allowed modes. Considering light power ${P_0}$ launched into forward-propagating mode m at the input end, the temporal power of the RBS power in mode n can be expressed as [28]

$$\begin{array}{c} {{P_{mn}}(t )= {P_0}{\alpha _s}\overline {{v_g}} \Delta \tau {e^{ - 2\alpha \overline {{v_g}} t}}{F_{mn}}} \end{array}$$
where $\overline {{v_g}} = \frac{{{v_{gm}}{v_{gn}}}}{{{v_{gm + }}{v_{gn}}}}$, ${v_{gm}}$ and ${v_{gn}}$ represent the group velocities of different modes which are correlated to intermodal dispersion and lead to deteriorate the space resolution of DAS system [29]; ${\alpha _s}$ represents the attenuation constant; ${F_{mn}}$ represents the overall capture fraction. Compared with SMF, MMF has been proved to own much higher ${F_{mn}}$ [30,31]. If we only consider MMF and the multimode characteristics, ${F_{mn}}$ quantifies the ratio of the scattered power into mode n over the total scattered power, is given by [28]
$$\begin{array}{c} {{F_{mn}} \propto \frac{{\mathop \smallint \nolimits_0^\infty \mathop \smallint \nolimits_0^{2\pi } {E_m}^2({r,\phi } ){E_n}^2({r,\phi } )d\phi rdr}}{{\mathop \smallint \nolimits_0^\infty \mathop \smallint \nolimits_0^{2\pi } {E_m}^2({r,\phi } )d\phi rdr\mathop \smallint \nolimits_0^\infty \mathop \smallint \nolimits_0^{2\pi } {E_n}^2({r,\phi } )d\phi rdr}}.} \end{array}$$

It can be deducted that the ${F_{mn}}$ is maximized when $m = n$ . Based on (1) and (2), the maximum RBS power in MMF can be received when the RBS mode is the same as mode injected to the input end. Since the fundamental mode (FM) is non-degenerate, it is barely coupled to the other modes and kept FM to the end of MMF. Single mode circulator has been proved to be an effect method for exciting the FM into the MMF and only pick the FM of the RBS signal at the receiving end. Therefore, the MMF-based DAS system can act similarly as the SMF-based system, which can be called quasi-single-mode state (QSM) operation [17].

If only take the scattering and attenuation effect in to account, the maximum sensing length ${L^{max}}$ of MMF can be expressed as (3), where ${\alpha _m}$ and ${\alpha _n}$ represent the attenuation constant of mode m and n, respectively [28]. As a result, the maximum sensing length of MMF is 17 km with an attenuation of ${\alpha _m}$ = 0.058 /km (0.25 dB/km) when $m = n$. In addition, for a graded-index MMF based DAS system in which dispersion induced influence on space resolution can be alleviated, the maximum sensing length is mainly determined by the signal frequency response bandwidth [32].

$$\begin{array}{c} {{L^{max}} = \left\{ {\begin{array}{c} {\frac{{\ln ({\alpha_m}/{\alpha_n})}}{{{\alpha_m} - {\alpha_n}}}\; \; \; \; \; \; \; {\alpha_m} \ne {\alpha_n};}\\ {\frac{1}{{{\alpha_m}}}\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; {\alpha_m} = {\alpha_n}\; ;} \end{array}} \right.} \end{array}$$

However, high order modes still exist due to bending or other disturbance. If we take the FM as the main mode of the injection light by a SMF circulator launching to MMF, $m = 1$ and ${E_m}$ can be approximately fixed, but ${E_n}$ may suffer from multimode interference. Though FM can be picked by spatial filtering, the optical field indeed composes of a summation of many modes [29]

$$\begin{array}{c} {{E_n}({z,t} )= \mathop \sum \limits_{i = 1}^N {E_i}\textrm{cos}[{2\pi \upsilon t - {\beta_i}z + {\varphi_i}(t )} ]+ \mathop \sum \limits_{j = 1}^N {E_j}\textrm{cos}[{2\pi \upsilon t - {\beta_j}z + {\varphi_j}(t )} ]} \end{array}$$
$$\begin{array}{c} {I = {E_n}{{({z,t} )}^2} \propto \mathop \sum \limits_{i = 1}^N E_i^2 + \mathop \sum \limits_{j = 1}^N E_j^2 + 2\mathop \sum \limits_{i = 1}^N \mathop \sum \limits_{j = 1}^N {E_i}{E_j}\cos [{z\Delta {\beta_{ij}} + \Delta {\varphi_{ij}}(t )} ]} \end{array}$$
where N (i, j = 1,2,…,N) represents the amount of the modes in MMF; z represents a longitude location; ${E_{i,j}}$ represents the amplitude; $\upsilon $ represents the optical frequency; $\Delta {\beta _{ij}} = ({{\beta_i} - {\beta_j}} )$, represents the propagation constant difference which is only 1% among modes [29]; $\mathrm{\Delta }{\varphi _{ij}}(t )= {\varphi _i}(t )- {\varphi _j}(t )$ represents the phase variation between modes. Based on (4) and (5), because the $\mathrm{\Delta }{\varphi _{ij}}$ only changes with $\Delta {\beta _{ij}}$ and hardly reaches constructive interference condition [33], the multimode interference will contribute to the signal of reception signal ${E_n}$ only for the same order mode ($i = j$).

2.2 Elastomer design

In order to further enhance the SNR of weak acoustic signal detection, the acoustic sensitive elastomer is explored in this work. According to the Hooke’s law, the relationship of the stress and strain of our proposed “n” type elastomer can be expressed as (4) and (5) [20]:

$$\begin{array}{c} {{\sigma _x} = \frac{Y}{{1 - {\mu ^2}}}({{\varepsilon_x} + \mu {\varepsilon_y}} )} \end{array}$$
$$\begin{array}{c} {{\sigma _y} = \frac{Y}{{1 - {\mu ^2}}}({{\varepsilon_y} + \mu {\varepsilon_x}} )} \end{array}$$
where Y is the Young’s modulus; $\mu $ is the Poisson’s ratio; ${\sigma _x}$, ${\sigma _y}$ are the circumferential and axial stress of the elastomer; ${\varepsilon _x}$, ${\varepsilon _y}$ are the circumferential and axial strain of the elastomer. Ignoring the displacement error of the fiber wound on the elastomer, the axial strain of the fiber can be determined by ${\varepsilon _x}$, which is relative to Y and $\mu $. Thus, the sensitivity of the elastomer S can be expressed as follows [20]:
$$\begin{array}{c} {S({T,R,\; H} )= \frac{{R{\varepsilon _x}({Y,\; \mu } )}}{P}} \end{array}$$
where P is the axial pressure. Except for Y and $\mu $, S is also related to wall thickness T, middle-layer external radius R, and height H of the elastomer.

Here we propose an “n” type elastomer as shown in Fig. 2(a), which is different from previous reported “I” type elastomer (as Refs. [25]) as shown in Fig. 2(b). In the simulation, an acoustic-structure coupling model has been built up. Specifically, the sound source is arranged at 10 mm above the top face of the elastomer. The meshing size is 1/6 of the acoustic wavelength. A reference section is selected in the center of the elastomer side wall to quantitatively calibrate the sensitivity. Given by (6), (7) and (8), we firstly investigate the relationship between material characteristics ($Y$, $\mu $) and S with fixed structural parameters. Figure 2(c) shows that the sensitivity is inversely proportional to the Young’s Modulus, the lower the Young’s Modulus, the higher the sensitivity. But it does not exhibit a simple proportional or inversely proportional relationship between the sensitivity and Poisson’s ratio. And the sensitivity is slightly influenced by the Poisson’s ratio sensitivity. Limited by characteristics of the exist material, a flexible material which owns 500 MPa Young’s Modulus and 0.48 Poisson’s ratio has been chosen.

 figure: Fig. 2.

Fig. 2. The structure of the (a) “n” type elastomer and (b) “I” type elastomer. The sensitivity optimizing procedure in (c) material and (d) structural parameters of “n” type elastomer in simulation.

Download Full Size | PDF

Then, we have explored neural network (NN) algorithms for fast multi-parameters (T, R, H) structural optimization. Firstly, the outliers that have existed physically but do not exist logically should be removed. Here, based on the density method, the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) clustering algorithm is used to define objects in low density regions as outliers [34]. Here, the density around an object is defined as the number of objects within a specified distance d. When d = 5 pm/Pa, we set the lower limit of the value density to 8, that is, within the range interval of 5 pm/Pa, data points less than 8 will be regarded as outlier, and they will be deleted directly. After sorting, data with sensitivity results greater than 135 pm/Pa will be deleted. Secondly, in the multi-layer perceptron (MLP) algorithm, 3543 sets of data (70%, stratified sampling) are used for training the model, and 1519 sets of data (30%, stratified sampling) are used for evaluating the results. The prediction model uses the perceptron neural network, in which each layer consists of the perceptron layer and Rectified Linear Unit (ReLU) layer. There are four layers in total. The numbers of middle layers are 256, 512 and 32, respectively. The model takes the T, R, H as inputs, S as labels. The loss function uses Mean Absolute Error (MAE), and the optimizer uses Adaptive Moment Estimation (Adam) for training. At the same time, cross validation and hyperparameter regulation training are carried out. After training is completed, the reconstruction parameters are stored in the network, and fast elastomer-sensitivity prediction can be achieved by simply inputting the structural parameters. It is worth mentioning that the MLP algorithm saves 99.76% calculation time from 4247 minutes to 10 minutes (including training time). Figure 2(d) shows the optimization result. In order to balance the manufacturing process and sensitivity performance, the truely optimized “n” type elastomer is shown in Fig. 2(a), with T of 3.0 mm, H of 3.0 mm and R of 15.0 mm.

Taking the occurring low frequency band inside the partial discharge into consideration, the sensitivity between 2.0 kHz to 2.9 kHz has been discussed. In particular, our proposed “n” type elastomer obtains 6.94 dB sensitivity gain at 2.5 kHz as shown in Fig. 3. It is worth mentioning that the resonant frequency of the elastomer is out of the simulation and experiment bandwidth.

 figure: Fig. 3.

Fig. 3. The sensitivity comparison of the optimized “n” type elastomer and the corresponding “I” type elastomer. The color of the insert figures represents the radial displacement field distribution.

Download Full Size | PDF

3. Experimental setup

Figure 4 shows the experimental setup of Φ-OTDR system for weak acoustic signal detection. Highly coherent light is produced by a narrow linewidth continuous-wave laser (NLL) source at 1550 nm with an isolator for protection. An optical fiber coupler with power ratio of 99:1 is used to separate the light into two parts. The part with 1% ratio is the reference path. The path with 99% light is launched to an acousto-optic modulator (AOM) which is driven by an arbitrary waveform generator (AWG). Continuous light is modulated into optical pulse with a 200-MHz frequency shift by the AOM. The output pulse has pulse width of 60 ns and a repetition rate of 25 kHz. An erbium-doped fiber amplifier (EDFA) with an amplifier spontaneous emission (ASE) noise filter is used to power boost. An optical fiber circulator is used to collect RBS signal, followed by a pre-EDFA. Mixing with the 1% reference light through a 50:50 optical fiber coupler, RBS light is converted to electrical signal by a balanced photo-detector (BPD) with a detection bandwidth of 350 MHz. Differential frequency component is recorded by a high-speed data acquisition card (DAQ) with 1 GHz for data analysis. A buzzer is induced to produce different acoustic signals.

 figure: Fig. 4.

Fig. 4. Experimental setup of the Φ-OTDR system. AOM: acousto-optic modulator; EDFA: erbium doped fiber amplifier; SMF: single mode fiber; MMF: multi-mode fiber; BPD: balanced photo-detector; DAQ: data acquisition card; PC: personal computer. Insert figure (a) is “I” type elastomer and (b) is “n” type elastomer.

Download Full Size | PDF

In the experiment, a sound level meter is used to calibrate the acoustic signal intensity of the buzzer as the variable. The corresponding sound pressure level (SPL) changes from 0.04 Pa to 0.16 Pa. In order to prove that MMF has higher SNR for weak acoustic signal detection than SMF, the performance of Φ-OTDR system with 2 spans of 1.1-km-long SMF and 2 spans of 1.2-km-long MMF without elastomers are compared in the experiment. Furthermore, both “I” type and “n” type elastomers have been inserted to the middle of the SMF and MMF fiber link, to compare and evaluate the SNR enhancement of the Φ-OTDR system.

4. Results and analysis

Based on a heterodyne coherent Φ-OTDR system, we firstly applied a 1.1-km-long SMF as the fiber under test (FUT). About 1.98 m of the fiber end has been wrapped in 7 circles. Different intervals away from the buzzer are used to simulate the contact and non-contact conditions, as shown in the insert figure of Fig. 5 (a). The buzzer emits 0.16-Pa acoustic signal at 5 kHz. Then a contact signal with the fiber circles bonding close to the buzzer and a non-contact signal with 10-mm space away from the fiber circles have been applied. The non-contact signal is extremely weak and hard to detect. Here, demodulation of the amplitude is used for location [35,36]. In details, the location SNR of the non-contact signal is 4.89 dB lower than the contact one, as shown in Fig. 5(a). And the frequency SNR of the non-contact signal is 8.12 dB lower than the contact one, as shown in Fig. 5(b). In order to obtain more accurate phase recovery, higher detection ability is needed.

 figure: Fig. 5.

Fig. 5. (a) Location and (b) frequency demodulation of contact and non-contact acoustic signals in same condition. Insert figure in (a) shows the contact and non-contact conditions.

Download Full Size | PDF

4.1 SNR enhancement of weak acoustic signal by MMF based Φ-OTDR system

Based on (1) and (2), the MMF is expected to have better SNR compared with SMF [30,31]. Here the comparison experiment for weak acoustic signal detection is developed. Both of the SMF and the MMF are bent in 9-cm diameters for 7 turns as the disturbance area. As shown in Fig. 6(a), when a 0.07-Pa, 2.5 kHz acoustic signal is emitted by a buzzer, the time domain acoustic signal demodulated by MMF based system shows better phase-recovery compared to SMF based system. Frequency domain acoustic signal obtained by Fourier transform achieves 9.90-dB higher SNR in MMF based system, as shown in the insert figure in Fig. 6(a). In order to evaluate different-intensity acoustic detection performance, we do the repeat experiments. The SNR of both systems under 0.04-Pa to 0.16-Pa acoustic signals detection has been investigated in Fig. 6(b). Each of the discrete point shows the mean value of 5 times experiments, and the lengths of the vertical lines shows the standard deviations. In the whole SPL range under test, MMF based system achieves higher frequency SNR. In the area SPL lower than 0.12 Pa, the SMF based system can hardly demodulate the acoustic frequency. Obviously, the MMF based system obtains wider detection range. MMF based system achieves the highest frequency SNR enhancement of 9.26 dB when the SPL is 0.07 Pa.

 figure: Fig. 6.

Fig. 6. (a) Time domain of a 0.07-Pa, 2.5 kHz signal by SMF and MMF based Φ-OTDR system without elastomer. Insert figure is frequency domain signal. (b)Frequency SNR of SMF and MMF based Φ-OTDR system with error bar.

Download Full Size | PDF

4.2 Acoustic frequency response sensitivity enhancement based on “n” type elastomer

The comparison experiment of the previous reported “I” type elastomer and our new proposed “n” type elastomer begins with the SMF-combined system. The same setup of disturbance has been induced. Besides, the middle of the 2 spans of 1.1-km-long SMF has been wound 6 turns for each of the elastomer. As shown in Fig. 7(a), maintaining the repetition frequency, the pulse width, and the input laser power, the “I” type elastomer integrated Φ-OTDR system is hard to locate the disturbance based on amplitude demodulation, while “n” type elastomer integrated Φ-OTDR system locates the position clearly. According to Fig. 7(b), we demonstrate the frequency-demodulation ability by analyzing the frequency SNR with error bar. With the help of the elastomers, both fiber links successfully demodulate the frequency of the acoustic signal. Among them, the “n” type elastomer helps the SMF based Φ-OTDR system to achieve 8.39 dB for location SNR and 11.02 dB for frequency SNR enhancement in maximum when the SPL is 0.16 Pa, specifically.

 figure: Fig. 7.

Fig. 7. (a) Location and (b) frequency SNR of SMF integrated with “I” and “n” type elastomers based Φ-OTDR system with error bar.

Download Full Size | PDF

4.3 Frequency response sensitivity enhancement scheme based on MMF combined with “n” type elastomer

We further combine the MMF and the “n” type elastomer to enhance the SNR. The experimental setup is the same as the comparison experiment in section B. MMF combined with both of the elastomers locates the disturbance area and demodulates the acoustic frequency clearly as shown in Figs. 8(a) and (b). The “n” type elastomer improves the location SNR and the frequency SNR of the MMF based Φ-OTDR system for 10.51 dB and 13.38 dB in maximum when the SPL is 0.10 Pa. After taking a polynomial fitting to the whole range of the existed data between 0.04 Pa to 0.16 Pa, the lower limit intensity of the proposed system can reach 0.01 Pa.

 figure: Fig. 8.

Fig. 8. (a) Location and (b) frequency SNR of MMF integrated with “I” and “n” type elastomers based Φ-OTDR system with error bar.

Download Full Size | PDF

Phase change can be calculated by differential phase of 40 consecutive pulses [37]. As shown in Fig. 9(a), MMF combined with “n” type elastomer system recovers the acoustic wave most precisely. Compared with the result shown in Fig. 6(a), the new proposed system proposed in this paper achieves significant improvement. It is worth mentioning that the MMF achieves phase recovery enhancement under serious bending conditions due to its higher sensitive refractive index change with bending radius than SMF. The acoustic pressure sensitivity is key point for phase enhancement. As illustrated in Fig. 9(b), a phase-pressure sensitivity of 18.57 rad/Pa (i.e., -94.62 dB re rad/µPa) has been achieved at 2.5 kHz by linear fitting. Besides, the integration of the “n” type elastomer can also be flexibly applied along the fiber link according to the demand of locally sensitivity enhancement. We confirm the feasibility of a MMF and an “n” type elastomer integrated Φ-OTDR system and illustrate a more stable SNR enhancement performance.

 figure: Fig. 9.

Fig. 9. (a) Time domain 0.16-Pa, 2.5 kHz acoustic signal detected by SMF and MMF integrated with “I” and “n” type elastomers based Φ-OTDR system. (b) Sensitivity of MMF integrated with “I” and “n” type elastomers based Φ-OTDR system.

Download Full Size | PDF

5. Conclusion

In order to overcome the difficulty of weak non-contact acoustic signals detection in DAS, a MMF integrated with an optimized “n” type elastomer based Φ-OTDR system has been proposed. On one hand, the power intensity of the RBS light along the whole fiber has been improved by inducing a MMF as the FUT. Firstly, the frequency SNR obtained 9.26-dB gain in MMF based system compared with SMF based system. Secondly, compared with the previous proposed “I” type elastomer, the new proposed “n” type elastomer obtains 8.39 dB and 11.02 dB on location and frequency SNR enhancement specifically by SMF based Φ-OTDR system. Then, these two values improve to be 10.51 dB and 13.38 dB by MMF and elastomers based Φ-OTDR system. Finally, a phase-pressure sensitivity of -94.62 dB re rad/µPa has been achieved at 2.5 kHz. MMF expands the range of acoustic detection intensity to a wider range. The “n” type elastomer further enhances the acoustic frequency response and phase recover ability. The adequate analysis proved that the proposed system will not perform competently in sensitivity enhancement of DAS without any of them. We indicate its great potential for weak acoustic signal detection. With a view of the present and future, highly sensitive quasi-distributed weak acoustic signal detection will show practical value in many applications.

Funding

National Natural Science Foundation of China (62220106006); Natural Science Foundation of Guangdong Province (No. 2022A1515011434); Stable Support Program for Higher Education Institutions from Shenzhen Science, Technology & Innovation Commission (20200925162216001); Basic and Applied Basic Research Foundation of Guangdong Province (2021B1515120013); Open Fund State Key Laboratory of Information Photonics and Optical Communications Beijing University of Posts and Telecommunications (No.IPOC2020A002); Open Projects Foundation of State Key Laboratory of Optical Fiber and Cable Manufacture Technology (No. SKLD2105); General Program of Shenzhen Science, Technology & Innovation Commission (JCYJ20220530113811026); Shenzhen Research Foundation (JCYJ20220818101206015, JSGG20220831103402004).

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. J. C. Juarez, E. W. Maier, K. N. Choi, and H. F. Taylor, “Distributed fiber-optic intrusion sensor system,” J. Lightwave Technol. 23(6), 2081–2087 (2005). [CrossRef]  

2. T. Yamate, G. Fujisawa, and T. Ikegami, “Optical Sensors for the Exploration of Oil and Gas,” J. Lightwave Technol. 35(16), 3538–3545 (2017). [CrossRef]  

3. H. Wu, Y. Qian, W. Zhang, and C. Tang, “Feature extraction and identification in distributed optical-fiber vibration sensing system for oil pipeline safety monitoring,” Photonic Sens. 7(4), 305–310 (2017). [CrossRef]  

4. M. Fernández-Ruiz, M. Soto, E. Williams, S. Martin-Lopez, Z. Zhan, M. Gonzalez-Herraez, and H. Martins, “Distributed acoustic sensing for seismic activity monitoring,” APL Photonics 5(3), 030901 (2020). [CrossRef]  

5. H. Ga and A. Eyal, “On the sensitivity of distributed acoustic sensing,” Opt. Lett. 41(24), 5648–5651 (2016). [CrossRef]  

6. G. Ma, H. Zhou, C. Shi, Y. Li, Q. Zhang, C. Li, and Q. Zheng, “Distributed Partial Discharge Detection in a Power Transformer Based on Phase-Shifted FBG,” IEEE Sensors J. 18(7), 2788–2795 (2018). [CrossRef]  

7. P. Zhu, H. Wen, Q. Che, and X. Li, “Disturbed partial discharge detection system based on an improved Φ-OTDR assisted by a wFBG array,” Appl. Opt. 59(14), 4367–4370 (2020). [CrossRef]  

8. Q. Yuan, F. Wang, T. Liu, Y. Zhang, and X. Zhang, “Using an Auxiliary Mach–Zehnder Interferometer to Compensate for the Influence of Laser-Frequency-Drift in Φ-OTDR,” IEEE Photonics J. 11(1), 1–9 (2019). [CrossRef]  

9. S. Liu, F. Yu, R. Hong, W. Xu, L. Shao, and F. Wang, “Advances in phase-sensitive optical time-domain reflectometry,” Opto-Electron. Adv. 5(3), 200078 (2022). [CrossRef]  

10. I. Ashry, B. Wang, Y. Mao, F. Hveding, A. Bukhamsin, T. Ng, and B. Ooi, “A Review of Distributed Fiber-Optic Sensing in the Oil and Gas Industry,” J. Lightwave Technol. 40(5), 1407–1431 (2022). [CrossRef]  

11. D. Tosi, C. Molardi, M. Sypabekova, and W. Blanc, “Enhanced Backscattering Optical Fiber Distributed Sensors: Tutorial and Review,” IEEE Sensors J. 21(11), 12667–12678 (2021). [CrossRef]  

12. P. Westbrook, K. Feder, R. Ortiz, T. Kremp, E. Monberg, H. Wu, D. Simoff, and S. Shenk, “Kilometer length, low loss enhanced back scattering fiber for distributed sensing,” in 2017 25th International Conference on Optical Fiber Sensors (OFS), South Korea, 1–5 (2017).

13. P. Healey, “Fading in Heterodyne OTDR,” Electron. Lett. 20(1), 30–32 (1984). [CrossRef]  

14. Y. Sun, H. Li, C. Fan, B. Yan, J. Chen, Z. Yan, and Q. Sun, “Review of a Specialty Fiber for Distributed Acoustic Sensing Technology,” Photonics 9(5), 277–299 (2022). [CrossRef]  

15. Z. Zhao, Y. Dang, M. Tang, L. Wang, L. Gan, S. Fu, C. Yang, W. Tong, and C. Lv, “Enabling Simultaneous DAS and DTS Through Space-Division Multiplexing Based on Multicore Fiber,” J. Lightwave Technol. 36(24), 5707–5713 (2018). [CrossRef]  

16. Z. Zhao, H. Wu, J. Hu, K. Zhu, Y. Dang, Y. Yan, M. Tang, and C. Lv, “Interference fading suppression in Φ-OTDR using space-division multiplexed probes,” Opt. Express 29(10), 15452–15462 (2021). [CrossRef]  

17. Y. Mao, I. Ashry, F. Hveding, A. Bukhamsin, Y. Hong, T. Ng, and B. Ooi, “Simultaneous Distributed Acoustic and Temperature Sensing Using a Multimode Fiber,” IEEE J. Sel. Top. Quantum Electron. 26(4), 1–7 (2020). [CrossRef]  

18. C. Wang, Y. Shang, W. Zhao, X. Liu, C. Wang, H. Yu, M. Hong, and G. Peng, “Distributed Acoustic Sensor Using Broadband Weak FBG Array for Large Temperature Tolerance,” IEEE Sensors J. 18(7), 2796–2800 (2018). [CrossRef]  

19. D. Liu, T. He, Z. Xu, and Q. Sun, “New Type of Microstructure-Fiber Distributed Acoustic Sensing Technology and Its Applications,” J. Appl. Sci. 38(2), 296–309 (2020).

20. H. Li, Q. Sun, T. Liu, C. Fan, T. He, Z. Yan, D. Liu, and P. Shum, “Ultra-High Sensitive Quasi-Distributed Acoustic Sensor Based on Coherent OTDR and Cylindrical Transducer,” J. Lightwave Technol. 38(4), 929–938 (2020). [CrossRef]  

21. J. Yao, B. Han, X. Jiang, S. Cao, Y. Rao, Z. Ran, W. Wang, H. Guan, and J. Long, “Cladding softened fiber for sensitivity enhancement of distributed acoustic sensing,” Opt. Express 29(6), 8216–8222 (2021). [CrossRef]  

22. S. Loranger, M. Gagné, V. Lambin-Iezzi, and R. Kashyap, “Rayleigh scatter based order of magnitude increase in distributed temperature and strain sensing by simple UV exposure of optical fibre,” Sci. Rep. 5(1), 11177–11183 (2015). [CrossRef]  

23. X. Chen, J. Himmelreich, J. Hurley, C. Zhou, Q. Jiang, Y. Qin, J. Li, C. Wu, H. Chen, D. Coleman, and M. Li, “Universal Fiber for Short-Distance Optical Communications,” J. Lightwave Technol. 37(2), 389–395 (2019). [CrossRef]  

24. T. Zhang, F. Pang, H. Liu, J. Cheng, L. Lv, X. Zhang, N. Chen, and T. Wang, “A Fiber-Optic Sensor for Acoustic Emission Detection in a High Voltage Cable System,” Sensors 16(12), 2026–2036 (2016). [CrossRef]  

25. Z. Chen, L. Zhang, H. Liu, P. Peng, Z. Liu, S. Shen, N. Chen, S. Zheng, J. Li, and F. Pang, “3D Printing Technique-Improved Phase-Sensitive OTDR for Breakdown Discharge Detection of Gas-Insulated Switchgear,” Sensors 20(4), 1045–1556 (2020). [CrossRef]  

26. B. Lu, B. Wu, J. Gu, J. Yang, K. Gao, Z. Wang, L. Ye, Q. Ye, R. Qu, X. Chen, and H. Cai, “Distributed optical fiber hydrophone based on Φ-OTDR and its field test,” Opt. Express 29(3), 3147–3162 (2021). [CrossRef]  

27. B. Lu, J. Gu, Z. Wang, L. Ye, Y. Liu, J. Yang, B. Wu, Q. Ye, R. Qu, and H. Cai, “Ultra-Low-Noise MIMO Distributed Acoustic Sensor Using Few-Mode Optical Fibers,” J. Lightwave Technol. 40(9), 3062–3071 (2022). [CrossRef]  

28. Z. Wang, H. Wu, X. Hu, N. Zhao, Q. Mo, and G. Li, “Rayleigh scattering in few-mode optical fibers,” Sci. Rep. 6(1), 35844–35851 (2016). [CrossRef]  

29. M. J. Murray, A. Davis, and B. Redding, “Multimode fiber Φ-OTDR with holographic demodulation,” Opt. Express 26(18), 23019–23030 (2018). [CrossRef]  

30. A. R. Mickelson and M. Eriksrud, “Theory of the backscattering process in multimode optical fibers,” Appl. Opt. 21(11), 1898–1909 (1982). [CrossRef]  

31. M. Nakazawa, “Rayleigh backscattering theory for single-mode optical fibers,” J. Opt. Soc. Am. 73(9), 1175–1180 (1983). [CrossRef]  

32. Y. Zhang, S. Fu, Y. Chen, Z. Ding, Y. Shan, F. Wang, M. Chen, X. Zhang, and Z. Meng, “A visibility enhanced broadband phase-sensitive OTDR based on the UWFBG array and frequency-division-multiplexing,” Opt. Fiber Technol. 53, 101995 (2019). [CrossRef]  

33. W. Wang, “The interference in multimode optical fiber interferometer,” Appl. Laser 10(2), 52–59 (1989).

34. E. Schubert, J. Sander, M. Ester, H. Kriegel, and X. Xu, “DBSCAN revisited, revisited: Why and how you should (still) use DBSCAN,” ACM Trans. Database Syst. 42(3), 1–21 (2017). [CrossRef]  

35. X. Mei, F. Pang, H. Liu, G. Yu, Y. Shao, T. Qian, C. Mou, L. Lv, and T. Wang, “Fast coarse-fine locating method for φ-OTDR,” Opt. Express 26(3), 2659–2667 (2018). [CrossRef]  

36. F. Pang, M. He, H. Liu, X. Mei, J. Tao, T. Zhang, X. Zhang, N. Chen, and T. Wang, “A Fading-Discrimination Method for Distributed Vibration Sensor Using Coherent Detection of φ -OTDR,” IEEE Photonics Technol. Lett. 28(23), 2752–2755 (2016). [CrossRef]  

37. H. Liu, F. Pang, L. Lv, X. Mei, Y. Song, J. Chen, and T. Wang, “True Phase Measurement of Distributed Vibration Sensors Based on Heterodyne φ -OTDR,” IEEE Photonics J. 10(1), 1–9 (2018). [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. Illustration of sensing head for weak acoustic signal detection by the combined scheme of a MMF and an “n” type elastomer.
Fig. 2.
Fig. 2. The structure of the (a) “n” type elastomer and (b) “I” type elastomer. The sensitivity optimizing procedure in (c) material and (d) structural parameters of “n” type elastomer in simulation.
Fig. 3.
Fig. 3. The sensitivity comparison of the optimized “n” type elastomer and the corresponding “I” type elastomer. The color of the insert figures represents the radial displacement field distribution.
Fig. 4.
Fig. 4. Experimental setup of the Φ-OTDR system. AOM: acousto-optic modulator; EDFA: erbium doped fiber amplifier; SMF: single mode fiber; MMF: multi-mode fiber; BPD: balanced photo-detector; DAQ: data acquisition card; PC: personal computer. Insert figure (a) is “I” type elastomer and (b) is “n” type elastomer.
Fig. 5.
Fig. 5. (a) Location and (b) frequency demodulation of contact and non-contact acoustic signals in same condition. Insert figure in (a) shows the contact and non-contact conditions.
Fig. 6.
Fig. 6. (a) Time domain of a 0.07-Pa, 2.5 kHz signal by SMF and MMF based Φ-OTDR system without elastomer. Insert figure is frequency domain signal. (b)Frequency SNR of SMF and MMF based Φ-OTDR system with error bar.
Fig. 7.
Fig. 7. (a) Location and (b) frequency SNR of SMF integrated with “I” and “n” type elastomers based Φ-OTDR system with error bar.
Fig. 8.
Fig. 8. (a) Location and (b) frequency SNR of MMF integrated with “I” and “n” type elastomers based Φ-OTDR system with error bar.
Fig. 9.
Fig. 9. (a) Time domain 0.16-Pa, 2.5 kHz acoustic signal detected by SMF and MMF integrated with “I” and “n” type elastomers based Φ-OTDR system. (b) Sensitivity of MMF integrated with “I” and “n” type elastomers based Φ-OTDR system.

Equations (8)

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

P m n ( t ) = P 0 α s v g ¯ Δ τ e 2 α v g ¯ t F m n
F m n 0 0 2 π E m 2 ( r , ϕ ) E n 2 ( r , ϕ ) d ϕ r d r 0 0 2 π E m 2 ( r , ϕ ) d ϕ r d r 0 0 2 π E n 2 ( r , ϕ ) d ϕ r d r .
L m a x = { ln ( α m / α n ) α m α n α m α n ; 1 α m α m = α n ;
E n ( z , t ) = i = 1 N E i cos [ 2 π υ t β i z + φ i ( t ) ] + j = 1 N E j cos [ 2 π υ t β j z + φ j ( t ) ]
I = E n ( z , t ) 2 i = 1 N E i 2 + j = 1 N E j 2 + 2 i = 1 N j = 1 N E i E j cos [ z Δ β i j + Δ φ i j ( t ) ]
σ x = Y 1 μ 2 ( ε x + μ ε y )
σ y = Y 1 μ 2 ( ε y + μ ε x )
S ( T , R , H ) = R ε x ( Y , μ ) P
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.