Abstract
Distributed optical fibre sensors deliver a map of a physical quantity along an optical fibre, providing a unique solution for health monitoring of targeted structures. Considerable developments over recent years have pushed conventional distributed sensors towards their ultimate performance, while any significant improvement demands a substantial hardware overhead. Here, a technique is proposed, encoding the interrogating light signal by a single-sequence aperiodic code and spatially resolving the fibre information through a fast post-processing. The code sequence is once forever computed by a specifically developed genetic algorithm, enabling a performance enhancement using an unmodified conventional configuration for the sensor. The proposed approach is experimentally demonstrated in Brillouin and Raman based sensors, both outperforming the state-of-the-art. This methodological breakthrough can be readily implemented in existing instruments by only modifying the software, offering a simple and cost-effective upgrade towards higher performance for distributed fibre sensing.
Similar content being viewed by others
Introduction
Distributed optical fibre sensors (DOFS)1 offer the capability to continuously inform on the spatial distribution of environmental quantities (such as temperature/strain2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, acoustic impedance17,18, refractive index19, etc.) along with an optical fibre. Most DOFS exploit natural scattering processes present in optical fibres20, such as Rayleigh2,3,4,10, spontaneous Raman5,6,11 and spontaneous/stimulated Brillouin9,12,13,14,15,16,17,18,19 scatterings. These backscattered optical signals can be interrogated in time2,3,4,5,6,7,8,16,17,18, frequency9,10,11, correlation12,13,19 or mixed14,15 domains, each with their proper advantages targeting specific applications1.
Time-domain approaches based on optical time-domain reflectometry (OTDR)21 turn out to be intrinsically best suited for long-distance sensing (e.g. tens of km), since the sensing range is limited by the laser coherence length in frequency-domain approaches1,10, and by background noises originating from uninterrogated fibre positions in correlation-domain approaches1,12. Time-domain approaches enable to spatially resolve the local environmental information by analysing the back-reflected response from an intense optical pulse. The pulse energy is positively related to the measurement signal-to-noise ratio (SNR), which in turn ultimately scales and trades-off all critical sensing specifications (i.e. spatial resolution (SR), measurand accuracy, sensing range and measurement time)1,22,23,24. For any conventional time-domain DOFS, the maximum attainable SNR turns out to be limited by fundamental physical constraints: whilst the maximum pulse peak power is limited by the onset of nonlinear effects25,26,27,28,29, the minimum noise power is determined by the signal bandwidth thus being irreducible whatever the denoising technique used30,31.
To further improve the optimised performance of conventional time-domain DOFS, the energy of the sensing response must be enhanced without impairing the SR nor activating nonlinear effects. This can be realised by using advanced techniques, such as distributed or lumped optical amplification32,33,34,35, smart signal modulation4,16,36,37,38,39, optical pulse coding40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61 and any combination of them62,63. Among these approaches, optical pulse coding is conceptually the most cost-effective, which operates by launching optical pulse trains into the sensing fibre and demodulating (decoding) the backscattered multi-pulse response through post-processing, delivering an SNR-improved equivalent single-pulse response. Extensive studies have been carried out to maximise the coding operation efficiency40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64 (i.e. maximising the SNR improvement, so-called coding gain, whilst minimising the hardware overhead and extra measurement time); however, all currently existing code types present fundamental and/or practical limitations. For instance, although codes formed by a single sequence40,56,57,58,59,60,61,62,64 are ideal for minimising measurement time and data storage, they either exhibit an imperfect peak-to-sidelobes ratio after the decoding40,56,64 or are periodic/cyclic codes57,58,59,60,61,62,64 that suffer from impairments induced by signal-dependent noise1 (resulting in a compromised coding gain that has not yet been quantified in literature). As more robust alternatives, distortion-free aperiodic codes constituted by several sequences, such as Golay41,42,43,44,45,46,47, Simplex codes48,49,50,51,52,63 and their derivatives54,55, have been mostly used for DOFS. Both codes can provide SNR enhancement with the same number of total acquisitions as in single-pulse DOFS; however, the additional operation time (e.g. the decoding time and the codeword switching time) and the larger occupied memory (for storing distinct coded traces) compromise the coding operation efficiency. More importantly, to reduce the decoding distortion, additional devices are used to alleviate the decaying envelope of the pulse sequences imposed by the gain saturation of erbium-doped fibre amplifiers (EDFAs)46,50,53. Nevertheless, the unevenness of the amplified pulse train cannot be fully suppressed53, and remaining distortions are still detrimental for certain types of codes (e.g. Golay codes53). All above-mentioned issues make the conceptual advantages of the coding technique to be hardly realised in practice to match the vast diversity of real-world applications.
Here, we propose an approach based on the concept of deconvolution, offering the possibility of a single-sequence aperiodic code and overcoming all above-mentioned drawbacks inherent to conventional codes. This enables, for the first time to the best of our knowledge, substantial performance improvement over optimised single-pulse DOFS without any hardware extension. The aperiodic code sequence is once forever optimised by a dedicatedly developed distributed genetic algorithm (DGA), here designated as genetic-optimised code (GO-code), offering a coding gain comparable to that of conventional codes such as Golay and Simplex. More importantly, the flexibility in the code length and the adaptability to the non-uniform amplitude over the pulse sequence (namely the pulse sequence envelop) secure a maximum coding efficiency for any given experimental condition, whilst the minimised operation time and the fast decoding process enables on-line real-time measurements, all being crucially advantageous over conventional codes. The proposed technique can be readily implemented in any classical time-domain DOFS by only modifying the software in relation to the pulse generation and post-processing (decoding), showing an unmatched practical significance. All these advantages are here experimentally validated using the most standard implementations of distributed sensors based on Brillouin optical time-domain analysis (BOTDA) and Raman OTDR (ROTDR). It must be emphasised that all demonstrations for both GO-code and single-pulse schemes are performed under fully optimised conditions, i.e. using maximised signal powers25,26,28 and minimised noise bandwidth (by matching the noise bandwidth to the signal bandwidth through digital filtering31), thus securing fair comparisons.
Results
Description of the proposed coding technique
The proposed technique can retrieve the single-pulse response with an improved SNR by only modifying the software of a standard single-pulse distributed sensor in relation to pulse generation (upper row in Fig. 1) and post-processing (bottom row in Fig. 1). In this case, an Nc-point discrete-time signal c(n), representing the digitally coded pulse sequence, is uploaded into a programmable pulse generator (e.g. field-programmable gate array, FPGA), which in turn converts c(n) into an electrical signal driving the optical setup. The integer n is linearly related to the time t through the sampling rate fs, so as t = n/fs. The digital coding sequence c(n) is once forever generated from a Nu-point binary unipolar sequence u(n) following the three steps described in Fig. 1. Here, u(n) contains 0’s and 1’s elements and must be dedicatedly designed, as will be elaborated in the next subsection.
Although the amplitude of c(n) is uniform in the digital domain, the coded optical pulse sequence launched into the fibre (blue dots in the middle row of Fig. 1), designated as cf(n), may exhibit a non-uniform amplitude envelop (grey dashed curve in the middle row of Fig. 1) due to the gain saturation of the EDFA used to boost the optical sequence power. Note that all sequences defined hereafter with a subscript ‘f’ are associated with such a non-uniform amplitude envelope. The acquired coded fibre response \(r_c^m(n)\) can be commonly expressed by the linear convolution65 between the optical code sequence cf(n) launched into the fibre and the fibre impulse response h(n), merged with a zero-mean additive noise ec(n)
Defining the fibre length Lh, the number of samples in the fibre impulse response h(n) is Nh = 2Lhfs/vg, which further yields the number of samples in the coded fibre response \(r_c^m(n)\) and the noise component ec(n) to be Nr = Nc + Nh − 1. In Eq. (1) the coded optical sequence cf(n) can be considered as the linear convolution between an Nd-point sequence df(n) and the single pulse signal p(n), where df(n) = f(n)d(n), f(n) being the envelope function of cf(n) determined by the EDFA gain saturation, as shown by the grey dashed curve in the middle row of Fig. 1. This way, Eq. (1) and its discrete frequency-domain representation can be expressed as
where F denotes the Nr-point discrete Fourier transform (DFT) operation65.
Here, we designate DF(k) as the decoding function, which is known since its time-domain representation, df(n), can be precisely and easily obtained by retrieving the energy of each optical pulse in the optical code sequence cf(n) measured at the fibre input (see middle row of Fig. 1). Such calibration can be made for a long measurement session, provided that the experimental stability is sufficient and there is no uncontrolled signal drift of significant importance in the instrumentation. Based on Eq. (2) and inverse DFT65, the targeted single-pulse response rs(n) (=p(n)⊗h(n)) can be retrieved by performing the following decoding process:
where \(r_s^d(n)\) is the decoded single-pulse response that contains the targeted single-pulse response rs(n) and a noise term affected by the decoding function DF(k). Equation (3) highlights that the non-distorted single-pulse fibre response rs(n) can be recovered regardless of the optical coding sequence distortion introduced by the EDFA gain saturation, unlike other conventional codes (e.g. Simplex and Golay) that strictly require a uniform optical code sequence, i.e. f(n) ≡ 1. In addition, the decoding process can be made very fast since the DFT can be computed via fast Fourier transform65.
In Eq. (3), the measure of whether the DF(k) attenuates or magnifies the original input noise is given by a noise scaling factor Q (for details see Supplementary Eqs. (8)–(12)), which represents the ratio of the noise variance after (\(\sigma _{de}^2\), see Supplementary Eq. (8)) and before decoding (\(\sigma _e^2\), see Supplementary Eq. (4))
where UF(k) is the Nr-point DFT of the sequence uf(n) (with an amplitude envelope given by the EDFA gain response), the latter being obtained by an Nx-point downsampling of df(n), as shown in the middle row of Fig. 1. This way the coding gain (usually defined as the reduction in the noise standard deviation, i.e. σe/σde, like all other codes for DOFS) resulting from our proposed code is \(G_c = \sqrt {1{\mathrm{/}}Q} \). It appears evident that to maximise the coding gain, it is required to dedicatedly design and optimise uf(n).
Design and optimisation of the code
The design and optimisation of uf(n) should firstly maximise the optical energy increment provided by the coded optical pulse sequence cf(n) with respect to the single optical pulse p(n). Such energy increment is represented by an energy enhancement factor \(F_E = \mathop {\sum}\nolimits_{n = 1}^{N_u} {u_f\left( n \right)} \), which is a generic parameter valid for any type of coding. Note that FE equals to the number of coded pulses only when the code assumes a uniform amplitude envelope (the case of other conventional codes). In general, the upper bound of FE for any code is ultimately limited by the onset of undesired physical phenomena that get more severe for a larger FE, such as high-order pump depletion53 and/or additional Brillouin-gain dependent noises (see Supplementary Note 3) in BOTDA, and amplified shot-noise in ROTDR (see Supplementary Note 4). However, for other conventional codes, FE (the number of coded pulses) is additionally constrained by their mathematical rules42,48, so that the selected FE is restricted to specific values that could be much smaller than the ideal upper bound. This limitation makes them less flexible and efficient than the proposed code.
After determining the optimum FE (the above mentioned upper bound value), the code length Nu should be adapted accordingly (i.e. the larger the FE, the longer the Nu). This essentially means optimising the ratio m = Nu/FE (m ≥ 1), which is related to the theoretical maximum coding gain Gc (derived based on Eq. (4), see Supplementary Eqs. (12)–(19)):
To evaluate Gc, Eq. (5) is compared with standard reference coding gain \(G_r = \sqrt {F_E{\mathrm{/}}2} \) offered by other conventional unipolar codes commonly used in DOFS, thus further yielding
It must be mentioned that not all standard reference coding gain Gr can be practically realised since not every value of FE is allowed by conventional codes as mentioned before, hence here Gr only represents a numerical reference to evaluate the performance of the proposed technique. Equation (6) suggests that the larger the value of m, the larger the upper limit of Gc/Gr, as illustrated by Fig. 2a in logarithmic scale. To make Gc/Gr closest possible to this upper limit, the detailed code distribution in uf(n) must be dedicatedly designed. This is equivalent to design 1’s and 0’s distribution in u(n), since uf(n) equals to u(n) being modulated by a known decaying envelope (see Fig. 1).
It turns out that the ideal design should make the off-peak spectral region of |UF(k)|2 (the power spectrum of the designed code) as flat as possible for given FE and m, as elaborated in Supplementary Note 1, which is extremely complex and cannot be analytically solved to the best of our knowledge. As an alternative, here a DGA66 is dedicatedly developed to efficiently search for optimal solutions, so that the code proposed here is designated as GO-code. The general concept and detailed implementation of the DGA are introduced in Methods, and its performance is elaborated in Supplementary Note 2 along with some intermediate results. The optimum value of m for a given energy enhancement factor FE is empirically found to be between 3 and 4, as illustrated in Fig. 2b. The reason comes from the fact that the high computational complexity in the case of long code lengths (large values of FE and m) leads to a reduction of the genetic optimisation efficiency. Following this criterion, the optimum values of Gc achieved by the DGA searching as a function of FE are shown in Fig. 2c for m = 3 (note that similar curves have also been obtained for m = 4, showing only minor differences in terms of gain). It can be found that when FE = 20, the coding gain Gc is 0.3 dB larger than the standard reference coding gain Gr; however, the larger the energy enhancement factor FE, the more demanding the computational capability, thus leading to a larger compromise of Gc with respect to Gr. Nevertheless, the maximum difference is always <0.5 dB, which can be considered negligible in most applications.
GO-coded BOTDA
The performance of the proposed GO-code for BOTDA is investigated and compared with an optimised standard single-pulse BOTDA using the same experimental configuration (see “Methods”). In the setup, the single-pulse and GO-code schemes can be readily switched by alternating the pulse and coding sequence driving the FPGA. For both schemes, the SR is set to 2 m, corresponding to a pulse duration of 20 ns and an optimised single-pulse peak Brillouin gain of 2.5% at the fibre near-end. Other detailed experimental parameters are listed in Table 1.
For the GO-code scheme, the optimum energy enhancement factor FE is determined to be ~40, to ensure that the coding gain at the fibre far-end is not degraded by additional Brillouin-gain-dependent noises (see details in Supplementary Note 3). Then, following the procedure described in Algorithm 1, and using parameters listed in Table 2, a 135-bit GO-code sequence is found by the proposed DGA. According to Eq. (4), such GO-code sequence is expected to offer a coding gain of Gc = 6.2 dB, being only 0.3 dB lower than the standard reference coding gain (Gr = 6.5 dB, \(\sqrt {40{\mathrm{/}}2} = 4.47\) times in linear scale). The GO-coded optical pulse sequence measured at the fibre input with a duty cycle of 20% is represented by the blue curve in Fig. 3a. Note that the actual peak power of the first pulse has been adjusted equal to that of optimised single pulse for a fair comparison, and is here normalised to 1 in the figure for the sake of clarity. The corresponding df(n) used for decoding is precisely retrieved based on the actual envelop of this sequence, as the red curve shown in Fig. 3a. Using such df(n) and based on Eq. (3), decoded single-pulse responses for each scanning frequency are obtained, with a total decoding time of 1.8 s that has no critical impact on the total measurement time (~2.8 min).
The decoded single-pulse response at the Brillouin resonance is shown by the red curve in Fig. 3b, together with BOTDA traces obtained from the single-pulse scheme using 1024 (blue) and 17795 (black) averages, respectively. The latter case is set as a reference, whose SNR improvement with respect to the 1024-averaged single-pulse measurement is equivalent to the one that could be ideally obtained by the designed GO-code scheme (\(10lg\sqrt {17,795{\mathrm{/}}1024} = 6.2\,{\rm{dB}}\)). The figure demonstrates that the decoded response remains undistorted compared to the reference trace, validating the adaptability of the GO-code scheme to an uneven pulse power distribution along the sequence. Due to the additional Brillouin-gain-dependent optical noises in the GO-coded system (see Supplementary Note 3), a clear difference in the noise behaviour between the decoded response and the reference trace (visualised through the thickness of all traces) can be observed, as quantified by the SNR profiles shown in Fig. 3c. The impact of these noises reduces with distance due to the fibre attenuation, and eventually becomes negligible (i.e. showing a tiny impact on the expected coding gain) near the fibre far-end thanks to the proper optimisation of FE. Although this behaviour has not yet been clearly documented, it is evident that such enhanced Brillouin-gain-dependent optical noises are inherent to any coded-BOTDA, being a non-exclusive effect in the GO-code proposed here, as experimentally verified for Simplex-coded BOTDA in Supplementary Fig. 9.
From the measured Brillouin loss spectrum, the Brillouin frequency shift (BFS) profile along the sensing fibre is estimated through a cross-correlation method67. The obtained BFS profile around a 5-m long hotspot placed at the fibre far-end is shown by the black curve in Fig. 3d, which matches well with a reference curve (blue) obtained from the single-pulse scheme with 2 m SR and 17,795 averages, and outperforms the single-pulse measurements with 1024 averages (red). To eventually verify the performance of the proposed GO-code, the BFS uncertainty profiles are calculated over five consecutive measurement for both single-pulse (blue) and GO-code (red) schemes, as shown in Fig. 3e. The frequency uncertainty at the fibre far-end obtained by the GO-code scheme (0.63 MHz) compared to the single-pulse scheme (2.65 MHz) is 4.1-fold enhanced, in good agreement with the theoretically evaluated coding gain (6.2 dB).
A similar comparison is then carried out at 1 m SR, which allows for a larger optimal energy enhancement factor FE = 200 due to the smaller optimised single-pulse Brillouin gain (0.5%) at the fibre near-end. In this case a 723-bit GO-code sequence is found by DGA, offering a theoretically evaluated coding gain of 9.3 dB, being 0.7 dB lower than the reference coding gain. Whilst behaviours similar to the case of 2 m SR are observed (see Supplementary Fig. 13a–c), an 8.59-fold BFS uncertainty reduction (from 18.9 MHz for single-pulse scheme to 2.2 MHz for GO-code scheme) is demonstrated, as shown in Fig. 3f, matching perfectly the theoretically calculated coding gain of 9.3 dB.
GO-coded ROTDR
In this section, the performance of the proposed GO-code for ROTDR is investigated and compared with an optimised single-pulse case, both based on the same experimental layout (see “Methods”).
Comparative measurements with 2 m SRs are carried out over a 39 km sensing range. A standard single-mode fibre (SMF) is used to eliminate any multimodal dispersion, thus avoiding SR impairments near the fibre far-end1. For both single-pulse and GO-code schemes, each temporal trace is acquired over 0.408 ms and averaged 2 million times, leading to a practical measurement time of 13.6 min. In the GO-code scheme, RZ format with a duty cycle of 22.2% is used to avoid the interaction between pulses in the code sequence and their forward Raman scattering components57. The energy enhancement factor FE for the GO-code ROTDR scheme is limited by the contamination of the amplified shot-noise (as elaborated in the Supplementary Note 4), and its optimal value is calculated to be FE = 44. Taking into account the decaying envelop of the pulse sequence imposed by the EDFA gain saturation, the DGA delivers a 177-bit GO-code with a theoretical coding gain Gc = 6.37 dB, being only 0.34 dB lower than the standard reference coding gain Gr = 6.71 dB. The normalised GO-coded pulse sequence measured at the EDFA output is shown by the blue curve in Fig. 4a, along with the reconstructed df(n) (red curve) that is required for decoding. The decoding time is found to be 10 ms, showing a negligible impact on the total measurement time as well.
The resolved temperature profiles obtained from the single-pulse and GO-code schemes are shown by the blue and red curves in Fig. 4b, respectively, indicating the clear performance enhancement provided by the proposed GO-code. Such a performance enhancement is then quantified by computing the corresponding SNR profiles, as shown in Fig. 4c. A coding gain of 6.33 dB at the fibre far-end is demonstrated, being in good agreement with theory. As anticipated, the SNR improvement at the fibre near-end is lower than the theoretical value due to the signal-power-dependent amplified shot-noise in the avalanche photodetector (APD), as elaborated in Supplementary Note 4. To verify the SR achieved by the GO-code scheme with ROTDR, the temperature of a 5-m-long hotspot introduced at the fibre far-end has been increased up to 70 °C, while the rest of the fibre is kept at the room temperature (25 °C). The temperature profiles obtained using the GO-code scheme is shown by the red curve in Fig. 4d, matching well with the trend of a reference temperature profile (blue curve) obtained by the single-pulse scheme with 2 m SR. Finally, Fig. 4e demonstrates a 4.27-fold temperature uncertainty enhancement at the fibre far-end (from 8.3 °C in the single-pulse case down to 1.9 °C with GO-code) thanks to the GO-code coding gain.
Measurements with 1 m SR are then performed with the same sensing fibre and measurement time. In this case, the optimised energy enhancement factor FE is determined to be 150, leading to a 450-bit code that takes into account the decaying trend of the sequence envelops. The GO-code delivered by the DGA exhibits a coding gain of Gc = 8.67 dB. Relevant results associated with 1 m SR are illustrated in Supplementary Fig. 14a–c. Benefiting from the coding gain, the GO-code scheme offers 3.9 °C temperature uncertainty at the fibre end, representing a 7.24-fold improvement with respect to the uncertainty of 28.2 °C obtained by the single-pulse scheme, as shown in Fig. 4f.
In addition to the advantage of simple implementation, the capability of fast temperature sensing enabled by the proposed GO-code is demonstrated with 2 m SR and a 10.2-km SMF, using the same 177-bit code as illustrated in Fig. 4a. Temporal traces are acquired over 0.12 ms and averaged 9000 times, leading to a practical measurement time of ~1 s and decoding time of 1.6 ms. It must be mentioned that this short measurement and decoding time is indeed an exclusive feature of the proposed code, and cannot be reached by any other conventional codes. Due to the APD-induced amplified shot noise, the SNR improvement at the fibre far-end is found to be 5.5 dB, as quantitatively shown by the SNR profiles at 10.2 km in Fig. 4c. This corresponds to a 3.5-fold temperature uncertainty reduction, as can be observed in Fig. 4e (at 10.2 km, the GO-code reduces the temperature uncertainty from 8.3 °C in the single-pulse case down to 2.4 °C).
Based on these characterisations, a fast on-line real-time measurement that monitors the temperature of water under heating (from 25 to 100 °C) is performed. Figure 5a shows the two-dimensional map of the retrieved temperature as a function of the fibre position and the acquisition time, in which the temperature evolution of the hotspot (at 10.176 km) under heating can be identified. Figure 5b shows the measured evolution of the hotspot temperature (red) as a function of time, in good agreement with that measured by the single-pulse scheme (blue) and reference curve (black, measured by a high-precision thermometer). Figure 5c shows the temperature difference between measurements and reference, verifying the 3.5-fold performance enhancement. To explicitly show the fast decoding capabilities of the proposed GO-code scheme, a video is included in the Supplementary Material, displaying a real-time measurement of fast heating and cooling processes over two fibre sections (2 and 5 m, separated by a 5 m unperturbed section) at the fibre far-end. In the video, figures are obtained by on-line real-time decoding (synchronised with the measurements).
Discussion
Whilst providing an SNR improvement comparable to sophisticated state-of-the-art techniques, the intrinsic mechanism of the proposed GO-code leads to overall advantages listed in Table 3:
-
1.
the launch of only one sequence into the sensing fibre, minimising operation time and data storage. Note that Cyclic codes share the same advantage, while Simplex and Golay codes demand additional codeword switching time and larger memory usage, as indicated in the two top rows of Table 3;
-
2.
a negligible post-processing (decoding) time (e.g. 1.6 ms decoding time with respect to the measurement time of 1 s for ROTDR), enabling real-time on-line fast sensing. The decoding time for each coding technique is summarised in the third row of Table 3, in which it can be found that the Cyclic coding has the fastest decoding process. For instance, in the case of a 100 km-long sensing range using a sampling rate of 250 MS/s and a 255-bit code (i.e. Nh = 250,000, Nu = 255), the evaluated decoding time of Cyclic is three times shorter than that of the proposed GO-code;
-
3.
an aperiodic structure avoiding detrimental effects existing in single-sequence periodic (cyclic) codes, as indicated in the fourth and fifth rows of Table 3:
-
a.
the impaired SNR all along with the sensing fibre due to the additional signal-dependent noise originating from the first kilometres of the fibre (theoretically justified in Supplementary Notes 3 and 4). This leads to a largely compromised coding gain at all fibre positions, as experimentally demonstrated using BOTDA in Supplementary Note 3;
-
b.
the absence of the continuous background observed in the measured cyclic-coded response, hindering any signal normalisation53 required before decoding (the signal normalisation is an essential process before decoding to extract the linear coded response from the measured data53);
-
4.
an adaptability to a non-uniform power over the pulse sequence (e.g. resulting from EDFA gain saturation), as indicated in the sixth row of Table 3, thus avoiding any modification of the instrumental layout compared to an optimised single-pulse implementation. This represents the ultimate asset of the technique, in full contrast with Simplex and Golay codes that require additional devices (higher cost) or pre-distortion operation (more complexity) to alleviate the unevenness of the coded sequence imposed by the EDFA gain saturation. Even with an improved uniformity (the amplitude envelope of the sequence is never perfectly uniform), decoding distortions remain still present when using Golay codes53, resulting in a degraded performance;
-
5.
enabling the use of any number of pulses in the code sequence, thus securing the highest coding efficiency in any given condition.
Benefiting from all these advantages, the proposed GO-coding technique enables, for the first time to the best of our knowledge, a pure overall performance improvement over optimised conventional DOFS without any hardware modification. The method only requires a software modification for sequence generation and post-processing, providing a fully compatible and cost-effective solution to improve the performance of existing instruments.
The performance of the proposed GO-code has been experimentally demonstrated by comparing with fully optimised single-pulse BOTDA and ROTDR (i.e. using a maximised signal power25,26,28, and minimised noise bandwidth by applying digital Gaussian filtering31). In the case of Raman sensing, further improvement in SNR could be achieved by replacing the sensing fibre by a few-mode fibre1,68,69, which enhances the allowable peak power of the pump light, with a moderate penalty in the SR at medium-long distances.
It is worth pointing out that the SNR enhancement provided by the GO-code and its features of single-sequence and fast decoding can lead to a huge reduction of the measurement time while maintaining other specifications. For instance, the 8.7-fold SNR enhancement demonstrated in the implemented GO-coded BOTDA sensor with 1 m SR can be used to reduce the number of traces averaging, leading to a measurement time reduction of factor 75 (=8.72), while keeping the same spatial and measurand resolutions as the single-pulse BOTDA implementation for the same sensing distance. This means that a single-pulse measurement of 5.7 min can be performed by a GO-coded system in only ~5 s, thus leading to a significant enhancement in the response time of the sensor with no penalty on other specifications.
It must be mentioned that the code sequences found by the DGA may not be optimal, which is a common feature shared in most searching algorithms66,70. In addition, due to the high computational complexity when searching long code sequences, the finite computational capability limits the searching efficiency. This means that the theoretical maximum coding gain indicated in Fig. 2a, which may even exceed the gain provided by conventional codes (so far only being observed in the case of small FE as shown in Fig. 2c), is still possible to be reached if the searching is carried out with more powerful tools and approaches.
Note also that, although the code studied here has been designed in unipolar format (involving 0’s and 1’s elements), the DGA can be readily adapted to search for codes in bipolar format (involving −1’s and 1’s elements) that can be applied to phase-encodable DOFS, such as phase-sensitive optical time-domain reflectometers (ϕ-OTDR)1,2. Furthermore, the presented approach can benefit not only DOFS but also diverse fields in technical science where the retrieval of the single-pulse response with high SNR is required, such as incoherent OTDRs for fault detection, laser range finders, lidar systems, among others.
Methods
BOTDA setup
The experimental configuration used for GO-coded BOTDA is shown in Fig. 6, where the hardware layout is strictly the same as that for a standard single-pulse BOTDA. In detail, a continuous-wave (CW) light from a distributed feedback laser operating at 1551 nm is split into probe and pump branches by a 50/50 coupler. In the probe branch (upper arm in the figure), the light travels through an electro-optic modulator driven by a microwave signal, producing a carrier-suppressed double sideband modulated probe wave. The frequency of the driving radiofrequency signal is scanned to reconstruct the BGS along the sensing fibre. This dual-sideband probe is launched into a 200 km-loop configuration63 with an optical power of 3 dBm/sideband. This loop configuration consists of two fibres: the 100 km-long leading fibre delivers a probe wave with −18 dBm/sideband to the far-end of the other 100 km-long sensing fibre that acts as a distributed sensing element, thus realising a real remoteness of 100 km for the sensing range63. An isolator is placed between the leading and sensing fibres to block the pump wave, thus restricting the Brillouin interaction only to the sensing fibre. In the pump branch (lower arm in the figure), the light is intensity-modulated to a single optical pulse or a coded optical pulse train, using a high extinction ratio (>50 dB) semiconductor optical amplifier (SOA) driven by a commercial FPGA. Preceded by an EDFA and a tuneable attenuator to optimise the pump wave power to ~23 dBm (determined by the onset of MI28), a polarisation scrambler is inserted to mitigate polarisation fading and polarisation pulling effects53.
In the receiver part, an EDFA operating in linear-gain regime is used between two tuneable narrowband fibre Bragg gratings (FBGs). The FBG1 selects the high-frequency probe sideband for the detection (i.e. making the system to operate in Brillouin loss configuration), thus preventing the EDFA from being saturated by unwanted frequency components, such as the Rayleigh backscattering and the second probe sideband. The FBG2, having the same filtering window as that of the FBG1, is used to filter out the amplified spontaneous emission (ASE) noise introduced by the EDFA, and thus to mitigate the ASE–ASE beating noise at the photo-detection. The selected optical signal is finally detected by a 125 MHz photodetector and acquired by a real-time data acquisition system with a sampling rate of 250 MS/s.
ROTDR setup
Similar to the BOTDA case, the experimental configuration used for GO-coded ROTDR is strictly the same as that for standard single-pulse ROTDR, as shown in Fig. 7. A CW light from a broad-linewidth (1.2 nm) light source operating at 1551 nm is injected into a high extinction ratio (>50 dB) SOA, driven by a commercial FPGA, to modulate the light intensity with a single optical pulse or a GO-coded pulse sequence, alternatively. The lightwave is then amplified by a high-gain EDFA to a peak power of 1 W and launched into the sensing fibre through a wavelength division multiplexer (WDM). The backscattered spontaneous anti-Stokes (AS) Raman signal is extracted by the same WDM and sent to an APD with a bandwidth of 125 MHz followed by a transimpedance amplifier. The photoelectrical signal is acquired by a real-time data acquisition card with a sampling rate of 250 MS/s.
Note that such a layout corresponds to a simplified ROTDR setup that makes use of a single APD and acquisition channel to measure only the AS spontaneous backward Raman scattering. This simplified setup allows us carrying out a proof-of-concept experiment in lab-controlled conditions. Wavelength-dependent losses are made negligible, so that the Raman Stokes or Rayleigh component is not required for real-time compensation. Instead, a simple pre-calibration procedure that measures the precise AS intensity at room temperature is performed to compensate for local and distributed losses along the sensing fibre.
Distributed genetic algorithm
Genetic algorithms (GAs)70 are efficient search and optimisation techniques based on the principle of natural genetics evolution, aiming at finding optimum or sub-optimum solutions to problems that are difficult to be algebraically solved. Standard GAs perform the search through a search space involving a population of individuals, each of which is randomly set following predefined initial conditions and represents a potential solution for the given problem. During the search process, the algorithm iteratively evaluates the fitness value of each individual and breeds new individuals by carrying out operations found in natural genetics. Once given termination conditions are satisfied, the searching process stops, and the individual exhibiting the best fitness value is output as the solution. The key idea of GAs is that the new generation of population bred after each iteration should involve individuals representing better solutions thanks to the use of classical operators inspired by the advancement in nature. These operators consist of: (1) selection, which selects the potential parents of offspring based on the evaluation of fitness value in such a way that better individuals are often selected for matching; (2) crossover, in which the genetic information from two parent individuals are exchanged to generate two offsprings; and (3) mutation, which creates new genetic information that does not exist in the parents. In order to expand the searching scope and to accelerate the search speed, DGAs66 that process multiple populations in parallel are often used. Furthermore, by periodically exchanging individuals among subpopulations via a migration operator, the premature convergence problem existing in conventional GAs can be greatly alleviated.
In this work a dedicatedly designed DGA is adapted to search for the optimal/sub-optimal unipolar binary code sequence u(n) that can provide the smallest possible noise scaling factor Q (the largest possible coding gain Gc) for these given input parameters:
-
1.
the energy enhancement factor FE bound by the physical constraints, as elaborated in Supplementary Notes 3 and 4;
-
2.
the value of m, which is a number between 3 and 4 according to the empirical results shown in Fig. 2b;
-
3.
the estimated envelope function f′(n) = uf(n)/u(n) imposed by the EDFA saturation, which can be estimated from the selected m and EDFA specifications. Note that f′(n) is close to the Nx-point downsampling of f(n), the latter being rigorously retrieved from the measured cf(n).
The pseudo-code for the designed DGA is sketched in Algorithm 1 and elaborated in the following, which basically consists of three sections: initialisation (line 1–2), genetic processing (line 4–10) and stopping criteria (line 3, 11–17).
Firstly, α subpopulations denoted as Si(i = 1, 2, …, α) are defined, each of which contains β Nu-bit-long all-zero sequences (individuals) that are designated as ui,j(i = 1, 2, …, α; j = 1, 2, …, β), where Nu = mFE. Then the initialisation is performed by randomly replacing FE ‘0’ bits by ‘1’ bits for each ui,j, and the corresponding noise scaling factors Qi,j are calculated according to Eq. (7). Then all subpopulations are processed in parallel by iteratively applying φ times the following genetic operators:
Selection: The selection operator updates a given subpopulation Si by selecting individuals based on the criterion that sequences with smaller Q (as defined by our target) should have a higher chance to survive. To do so, a typical roulette wheel selection method70 is implemented, where the survival probability (i.e. the probability passing to the next generation) of each sequence ui,j in a certain subpopulation Si is expressed as:
from which it can be found that the smaller the Qi,j, the larger the P(ui,j). Using such a selection algorithm, the subpopulation is updated by carrying out the following selection steps for β rounds: (i) generate a random number Pr between 0 and 1, (ii) add P(ui,j) to the partial sum \(\mathop {\sum}\nolimits_j {P(u_{i,j})} \), starting from the first sequence ui,1of the given subpopulation till the condition \(\mathop {\sum}\nolimits_j {P(u_{i,j})} \ge P_r\) is satisfied, and (iii) select the first sequence that makes \(\mathop {\sum}\nolimits_j {P(u_{i,j})} \) exceeding Pr as the survival individual for the current round.
Crossover: The crossover operator is analogous to gene recombination, which produces offsprings by exchanging genetic material (i.e. alternating segments of code sequences) of the parents. The length of alternating segments is here set as \(N_u^c\) bits and one-point crossover70 is applied, in which two sequences in a certain subpopulation are randomly chosen and their \(N_u^c\)-bit-long segments starting from a random crossover point are swapped, generating two new offsprings. This process is repeated by \(\beta P_{c_i}\) times for a certain subpopulation Si, where β is the number of individuals (sequences) and \(P_{c_i}\) is the probability of performing the crossover action.
Mutation: Mutation is an occasional tweak in the bit value of a sequence in order to maintain and introduce diversity in the subpopulation. This is usually applied with a low probability Pm, because if the probability is set very high, the GA would regress to an undesired random search. For the unipolar binary code sequence in our proposal, mutation is carried out by flipping the value of \(N_u^m\) randomly chosen bits of a sequence, which is repeated by \(\beta P_{m_i}\) times for a certain subpopulation Si.
After finishing these φ-round iterations, the migration operation is subsequently performed to exchange the information among updated subpopulations. This is realised by transferring γ sequences with smaller Q in a source subpopulation to replace the same number of sequences with higher Q in a target subpopulation. Here each subpopulation Si is in turn defined as the source of the neighbouring subpopulation Si+1, for 1 ≤ i < α, while Sα is the source subpopulation of S1. This way the information of ‘good’ sequences is spread among subpopulations; furthermore, the immigrants (i.e. γ sequences) would effectively change the fitness landscape thus avoiding the potential problem of premature convergence.
In parallel with the migration operation, the sequence with minimum Q (Qmin) among all subpopulations is selected as the best sequence ubest for the current generation (line 11 in Algorithm 1). The searching process is terminated when the converging condition is satisfied, i.e. when there has been no improvement of ubest for μ consecutive generations. This is indicated by a counter that is initialised to zero and counts the number of the generations for which there has been no improvement of ubest. As shown in lines 12–16 of the pseudo-code, the counter increases every time a better ubest is not generated, which is however reset to zero if ubest is updated and \(\mathop {\sum}\nolimits_{n = 1}^{N_u} {f_d^\prime (n)u_{\rm{best}}} \left( n \right) \approx F_E\). The algorithm terminates when the counter value reaches μ and the ubest of the last generation is output as the final solution.
Algorithm 1
The distributed genetic algorithm.
Input: energy enhancement factor FE, total bit number of a sequence Nu = mFE, estimated envelop function f′, number of subpopulations α, number of individuals in each subpopulation β, crossover probability Pc, bit number of alternating segments \({N_u^c}\), mutation probability Pm, bit number of mutation \({N_u^m}\), migration interval φ, number of migrants γ, generation counter μ.
Output: the best sequence ubest and the minimum noise scaling factor Qmin.
1 Randomly initialise α × β sequences u, where α is the number of subpopulation and each subpopulation contains β sequences (individuals);
2 uf = f′u, and evaluate the noise scaling factor Q for each uf by Eq. (7);
3 while counter < μ do
4 for each subpopulation Si(i = 1, 2, …, α) do φ times
5 Si ← Selection (Si, Qi,1, Qi,2, …, Qi,β);
6 Si ← Crossover (Si, \({P_{c_i}}\), \({N_u^c}\));
7 Si ← Mutation (Si, \({P_{m_i}}\), \({N_u^m}\));
8 Update the noise scaling factor Qi,1, Qi,2, …, Qi,β based on uf = f′u and Eq. (7);
9 end for
10 (S, Q) ← Migration (S, Q, γ);
11 (ubest, Qmin) ← (S, Q);
12 if ubest is updated and \(\mathop {\sum}\nolimits_{{\boldsymbol{n}} = 1}^{{\boldsymbol{N}}_{\boldsymbol{u}}} {{\boldsymbol{f}}\prime ({\boldsymbol{n}}){\boldsymbol{u}}_{{\mathbf{best}}}\left( {\boldsymbol{n}} \right) \approx {\boldsymbol{F}}_{\boldsymbol{E}}} \), then
13 counter ← 0;
14 else
15 counter ← counter + 1;
16 end if
17 end while
Note that the size of subpopulation, the parameters of genetic operations and the termination condition are all adjustable when using DGA. It has been verified experimentally that the proposed DGA with appropriate parameters is efficient in the search of optimal code sequence. Benefiting from the robustness of DGA, it is not necessary to tune all these parameters for a specific sequence length. All the numerical results in this paper are obtained with the single set of parameters shown in Table 2.
Code availability
The code of distributed genetic algorithm is deposited in Zenodo (https://doi.org/10.5281/zenodo.3909600).
References
Hartog, A. H. Introduction to Distributed Optical Fiber Sensors (CRC Press, 2017).
Muanenda, Y. Recent advances in distributed acoustic sensing based on phase-sensitive optical time domain reflectometry. J. Sens. 2018, 1–16 (2018).
Sladen, A. et al. Distributed sensing of earthquakes and ocean-solid Earth interactions on seafloor telecom cables. Nat. Commun. 10, 5777 (2019).
Williams, E. F. et al. Distributed sensing of microseisms and teleseisms with submarine dark fibers. Nat. Commun. 10, 5778 (2019).
Dakin, J. P., Pratt, D. J., Bibby, G. W. & Ross, J. N. Distributed optical fibre Raman temperature sensor using a semiconductor light source and detector. Electron. Lett. 21, 569–570 (1985).
Hartog, A. H. & Leach, A. P. Distributed temperature sensing in solid-core fibres. Electron. Lett. 21, 1061–1062 (1985).
Horiguchi, T. & Tateda, M. BOTDA-nondestructive measurement of single-mode optical fiber attenuation characteristics using Brillouin interaction: theory. J. Lightwave Technol. 7, 1170–1176 (1989).
Shimizu, K., Horiguchi, T., Koyamada, Y. & Kurashima, T. Coherent self-heterodyne detection of spontaneously Brillouin-scattered light waves in a single-mode fiber. Opt. Lett. 18, 185–187 (1993).
Garus, D., Gogolla, T., Krebber, K. & Schliep, F. Distributed sensing technique based on Brillouin optical-fiber frequency-domain analysis. Opt. Lett. 21, 1402–1404 (1996).
Passy, R., Gisin, N., von der Weid, J. P. & Gilgen, H. H. Experimental and theoretical investigations of coherent OFDR with semiconductor laser sources. J. Lightwave Technol. 22, 1622–1630 (1994).
Farahani, M. A. & Gogolla, T. Spontaneous Raman scattering in optical fibers with modulated probe light for distributed temperature Raman remote sensing. J. Lightwave Technol. 17, 1379–1391 (1999).
Zadok, A. et al. Random-access distributed fiber sensing. Laser Photon 6, L1–L5 (2012).
Mizuno, Y., Hayashi, N., Fukuda, H., Song, L. Y. & Nakamura, K. Ultrahigh-speed distributed Brillouin reflectometry. Light Sci. Appl. 5, e16184 (2016).
Elooz, D., Antman, Y., Levanon, N. & Zadok, A. High-resolution long-reach distributed Brillouin sensing based on combined time-domain and correlation-domain analysis. Opt. Express 22, 6453–6463 (2014).
Denisov, A., Soto, M. & Thévenaz, L. Going beyond 1000000 resolved points in a Brillouin distributed fiber sensor: theoretical analysis and experimental demonstration. Light Sci. Appl. 5, e16074 (2016).
Zhou, D. et al. Single-shot BOTDA based on an optical chirp chain probe wave for distributed ultrafast measurement. Light Sci. Appl. 7, 32 (2018).
Chow, D. M., Yang, Z., Soto, M. A. & Thévenaz, L. Distributed forward Brillouin sensor based on local light phase recovery. Nat. Commun. 9, 2990 (2018).
Bashan, G., Diamandi, H. H., London, Y., Preter, E. & Zadok, A. Optomechanical time-domain reflectometry. Nat. Commun. 9, 2991 (2018).
Bashan, G., London, Y., Diamandi, H. H. & Zadok, A. Distributed cladding mode fiber-optic sensor. Optica 7, 85–92 (2020).
Boyd, R. W. Nonlinear Optics (Elsevier, 2003).
Barnoski, M. K., Rourke, M. D., Jensen, S. M. & Melville, R. T. Optical time domain reflectometer. Appl. Opt. 16, 2375–2379 (1977).
Soto, M. A. & Thévenaz, L. Modeling and evaluating the performance of Brillouin distributed optical fiber sensors. Opt. Express 21, 31347–31366 (2013).
Gabai, H. & Eyal, A. On the sensitivity of distributed acoustic sensing. Opt. Lett. 41, 5648–5651 (2016).
Wang, S., Yang, Z., Soto, M. A. & Thévenaz, L. Study on the signal-to-noise ratio of Brillouin optical-time domain analyzers. Opt. Express 28, 19864–19876 (2020).
Wait, P. C., De Souza, K. & Newson, T. P. A theoretical comparison of spontaneous Raman and Brillouin based fibre optic distributed temperature sensors. Opt. Commun. 144, 17–23 (1997).
Thévenaz, L., Foaleng, S. M. & Lin, J. Effect of pulse depletion in a Brillouin optical time-domain analysis system. Opt. Express 21, 14017–14035 (2013).
Martins, H. F. et al. Modulation instability-induced fading in phase-sensitive optical time-domain reflectometry. Opt. Lett. 38, 872–874 (2013).
Alem, M., Soto, M. A. & Thévenaz, L. Analytical model and experimental verification of the critical power for modulation instability in optical fibers. Opt. Express 23, 29514–29532 (2015).
Dominguez-Lopez, A. et al. Novel scanning method for distortion-free botda measurements. Opt. Express 24, 10188–10204 (2016).
Soto, M. A., Ramírez, J. & Thévenaz, L. Intensifying the response of distributed optical fibre sensors using 2D and 3D image restoration. Nat. Commun. 7, 10870 (2016).
Zaslawski, S., Yang, Z. & Thévenaz, L. On the 2D post-processing of Brillouin optical time-domain analysis. J. Lightwave Technol. https://doi.org/10.1109/JLT.2020.2967091 (2020).
Rodriguez-Barrios, F. et al. Distributed Brillouin fiber sensor assisted by first-order Raman amplification. J. Lightwave Technol. 28, 2162–2172 (2010).
Jia, X. et al. Hybrid distributed Raman amplification combining random fiber laser based 2nd-order and low-noise LD based 1st-order pumping. Opt. Express 21, 24611–24619 (2013).
Gyger, F., Rochat, E., Chin, S., Niklès, M. & Thévenaz, L. Extending the sensing range of Brillouin optical time-domain analysis up to 325 km combining four optical repeaters. In Proc. SPIE 9157, 23rd International Conference on Optical Fibre Sensors, 91576Q. https://doi.org/10.1117/12.2059590 (2014).
Urricelqui, J., Sagues, M. & Loayssa, A. Brillouin optical time-domain analysis sensor assisted by Brillouin distributed amplification of pump pulses. Opt. Express 23, 30448–30458 (2015).
Soto, M. A. et al. Time and frequency pump-probe multiplexing to enhance the signal response of Brillouin optical time-domain analyzers. Opt. Express 22, 28584–28595 (2014).
Ruiz-Lombera, R. et al. Overcoming nonlocal effects and Brillouin threshold limitations in Brillouin optical time-domain sensors. IEEE Photonics J. 7, 1–9 (2015).
Chen, D., Liu, Q., Fan, X. & He, Z. Distributed fiber-optic acoustic sensor with enhanced response bandwidth and high signal-to-noise ratio. J. Lightwave Technol. 35, 2037–2043 (2017).
Yang, Z., Soto, M. A., Chow, D. M., Ray, P. & Thévenaz, L. Brillouin distributed optical fiber sensor based on a closed-loop configuration. J. Lightwave Technol. 36, 1239–1248 (2018).
Everard, J. K. A. Novel signal processing techniques for enhanced OTDR sensors. Fiber Opt. Sens. II 0798, 42–48 (1987).
Golay, M. Sieves for low autocorrelation binary sequences. IEEE Trans. Inf. Theory 23, 43–51 (1977).
Nazarathy, M. et al. Real-time long range complementary correlation optical time domain reflectometer. J. Lightwave Technol. 7, 24–38 (1989).
Soto, M. A. et al. Distributed temperature sensor system based on Raman scattering using correlation-codes. Electron. Lett. 43, 862–864 (2007).
Soto, M. A., Le Floch, S. & Thévenaz, L. Bipolar optical pulse coding for performance enhancement in BOTDA sensors. Opt. Express 21, 16390–16397 (2013).
Yang, Z., Soto, M. A. & Thévenaz, L. Increasing robustness of bipolar pulse coding in Brillouin distributed fiber sensors. Opt. Express 24, 586–597 (2016).
Wang, F., Zhu, C., Cao, C. & Zhang, X. Enhancing the performance of BOTDR based on the combination of FFT technique and complementary coding. Opt. Express 25, 3504–3513 (2017).
Wang, Z. et al. Distributed acoustic sensing based on pulse-coding phase-sensitive OTDR. IEEE Internet Things J. 6, 6117–6124 (2019).
Jones, M. D. Using simplex codes to improve OTDR sensitivity. IEEE Photon. Technol. Lett. 5, 822–824 (1993).
Lee, D., Yoon, H., Kim, P., Park, J. & Park, N. Optimization of SNR improvement in the noncoherent OTDR based on simplex codes. J. Lightwave Technol. 24, 322–328 (2006).
Soto, M. A., Bolognini, G., Pasquale, F. D. & Thévenaz, L. Long-range Brillouin optical time-domain analysis sensor employing pulse coding techniques. Meas. Sci. Technol. 21, 094024 (2010).
Soto, M. A., Bolognini, G. & Di Pasquale, F. Analysis of pulse modulation format in coded BOTDA sensors. Opt. Express 18, 14878–14892 (2010).
Soto, M. A., Bolognini, G. & Di Pasquale, F. Long-range simplex-coded BOTDA sensor over 120 km distance employing optical preamplification. Opt. Lett. 36, 232–234 (2011).
Yang, Z., Li, Z., Zaslawski, S., Thévenaz, L. & Soto, M. A. Design rules for optimizing unipolar coded Brillouin optical time-domain analyzers. Opt. Express 26, 16505–16523 (2018).
Le Floch, S., Sauser, F., Llera, M. & Rochat, E. Novel Brillouin optical time-domain analyzer for extreme sensing range using high-power flat frequency-coded pump pulses. J. Lightwave Technol. 33, 2623–2627 (2015).
Li, Z., Yang, Z., Yan, L., Soto, M. A. & Thévenaz, L. Hybrid Golay-coded Brillouin optical time-domain analysis based on differential pulses. Opt. Lett. 43, 4574–4577 (2018).
London, Y. et al. High-resolution long-range distributed Brillouin analysis using dual-layer phase and amplitude coding. Opt. Express 22, 27144–27158 (2014).
Soto, M. A. et al. Raman-based distributed temperature sensor with 1 m spatial resolution over 26 km SMF using low-repetition-rate cyclic pulse coding. Opt. Lett. 36, 2557–2559 (2011).
Taki, M. et al. Cyclic pulse coding for fast BOTDA fiber sensors. Opt. Lett. 38, 2877–2880 (2013).
Martins, H. F. et al. Real time dynamic strain monitoring of optical links using the backreflection of live PSK data. Opt. Express 24, 22303–22318 (2016).
Shiloh, L., Levanon, N. &. Eyal, A. Highly-sensitive distributed dynamic strain sensing via perfect periodic coherent codes. In Proc. 26th International Conference on Optical Fiber Sensors, paper TuE25. https://doi.org/10.1364/OFS.2018.TuE25 (2018).
Le Floch, S. et al. Unique coded sequence for fast Brillouin distributed sensors. In Proc. 26th International Conference on Optical Fiber Sensors, paper WF20. https://doi.org/10.1364/OFS.2018.WF20 (2018) .
Iribas, H., Loayssa, A., Sauser, F., Llera, M. & Le Floch, S. Cyclic coding for Brillouin optical time-domain analyzers using probe dithering. Opt. Express 25, 8787–8800 (2017).
Soto, M. A. et al. Extending the real remoteness of long-range Brillouin optical time-domain fiber analyzers. J. Lightwave Technol. 32, 152–162 (2014).
Levanon, N., Cohen, I., Arbel, N. & Zadok, A. Non-coherent pulse compression—aperiodic and periodic waveforms. IET Radar, Sonar Navigation 10, 216–224 (2016).
Proakis, J. G. & Manolakis, D. G. Digital Signal Processing: Principles, Algorithms, and Applications 4th edn (Pearson Education, Inc., 2007).
Tanese, R. Distributed genetic algorithms. In Proc. Third International Conference on Genetic Algorithms, 434–439 (1989).
Haneef, S. M., Yang, Z., Thévenaz, L., Venkitesh, D. & Srinivasan, B. Performance analysis of frequency shift estimation techniques in Brillouin distributed fiber sensors. Opt. Express 26, 14661–14677 (2018).
Wang, M. et al. Few-mode fiber based Raman distributed temperature sensing. Opt. Express 25, 4907–4916 (2017).
Liu, Y., Ma, L., Yang, C., Tong, W. & He, Z. Long-range Raman distributed temperature sensor with high spatial and temperature resolution using graded-index few-mode fiber. Opt. Express 26, 20562–20571 (2018).
Sivanandam S. & Deepa S. Genetic Algorithms (Springer, 2008).
Acknowledgements
This study is supported by National Natural Science Foundation of China under Grant (61875018) and National Key R&D Programme of China under Grant (2017YFB0405500). M.A. Soto acknowledges the support from AC3E ANID-Basal Project FB0008. The authors thank Jianqi Hu from EPFL for the fruitful discussion on the coding optimisation.
Author information
Authors and Affiliations
Contributions
X.H. conceived the concept of deconvolution. Z.Y. and S.Z. developed the concept, established the theoretical model and designed the methodology. X.S. and Z.Y. developed the distributed genetic algorithm. X.S., Z.Y., X.H. and S.W. contributed to all experiments and analysis. M.A.S contributed to Raman experiments and background discussions. X.G. contributed to Brillouin noise analysis. X.S. and Z.Y. wrote the paper. M.A.S., X.H. and L.T. revised the paper. X.H., J.W. and L.T. supervised the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Avi Zadok and the other, anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sun, X., Yang, Z., Hong, X. et al. Genetic-optimised aperiodic code for distributed optical fibre sensors. Nat Commun 11, 5774 (2020). https://doi.org/10.1038/s41467-020-19201-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-19201-1
This article is cited by
-
On-chip silicon electro-optical modulator with ultra-high extinction ratio for fiber-optic distributed acoustic sensing
Nature Communications (2023)
-
Chaos Raman distributed optical fiber sensing
Light: Science & Applications (2023)
-
Integrated sensing and communication in an optical fibre
Light: Science & Applications (2023)
-
Physics and applications of Raman distributed optical fiber sensing
Light: Science & Applications (2022)
-
High-Performance Distributed Brillouin Optical Fiber Sensing
Photonic Sensors (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.