Slope-assisted fast distributed sensing in optical fibers with arbitrary Brillouin profile

Similar documents
Fast and Distributed Brillouin Sensing for Dynamic SHM

Differential measurement scheme for Brillouin Optical Correlation Domain Analysis

High-resolution long-reach distributed Brillouin sensing based on combined time-domain and correlation-domain analysis

Dynamic Distributed Brillouin Optical Fiber Sensing Based on Dual-Modulation by Combining Single Frequency Modulation and Frequency-Agility Modulation

A Hybrid Φ/B-OTDR for Simultaneous Vibration and Strain Measurement

Analysis of Stimulated Brillouin Scattering Characteristics in Frequency Domain

Analysis of pulse modulation format in coded BOTDA sensors

BOTDA using OFDM channel estimation

Brillouin optical time-domain analysis sensor with pump pulse amplification

R. J. Jones College of Optical Sciences OPTI 511L Fall 2017

LASER & PHOTONICS REVIEWS. Random-access distributed fiber sensing ARTICLE LETTER

Optical Society of America: Open Access Journals

Slow light fiber systems in microwave photonics

Timing Noise Measurement of High-Repetition-Rate Optical Pulses

Self-advanced fast light propagation in an optical fiber based on Brillouin scattering

S-band gain-clamped grating-based erbiumdoped fiber amplifier by forward optical feedback technique

Extended delay of broadband signals in stimulated Brillouin scattering slow light using synthesized pump chirp

Performance Limitations of WDM Optical Transmission System Due to Cross-Phase Modulation in Presence of Chromatic Dispersion

Elimination of Self-Pulsations in Dual-Clad, Ytterbium-Doped Fiber Lasers

Photonic time-stretching of 102 GHz millimeter waves using 1.55 µm nonlinear optic polymer EO modulators

Directly Chirped Laser Source for Chirped Pulse Amplification

SUPPLEMENTARY INFORMATION

Suppression of Stimulated Brillouin Scattering

ModBox - Spectral Broadening Unit

Sharp tunable optical filters based on the polarization attributes of stimulated Brillouin scattering

SUPPLEMENTARY INFORMATION DOI: /NPHOTON

Optical Phase Lock Loop (OPLL) with Tunable Frequency Offset for Distributed Optical Sensing Applications

PHASE TO AMPLITUDE MODULATION CONVERSION USING BRILLOUIN SELECTIVE SIDEBAND AMPLIFICATION. Steve Yao

Optical fiber-fault surveillance for passive optical networks in S-band operation window

UTILITY APPLICATIONS OF FIBER-OPTIC DISTRIBUTED STRAIN AND TEMPERATURE SENSORS

Testing with Femtosecond Pulses

All-Optical Signal Processing and Optical Regeneration

Phase-Sensitive Optical Time-Domain Reflectometry Amplified by Gated Raman Pump

ModBox-SB-NIR Near Infra Red Spectral Broadening Unit

Supplementary Figures

A broadband fiber ring laser technique with stable and tunable signal-frequency operation

Investigations on the performance of lidar measurements with different pulse shapes using a multi-channel Doppler lidar system

Simultaneous Measurements for Tunable Laser Source Linewidth with Homodyne Detection

RADIO-OVER-FIBER TRANSPORT SYSTEMS BASED ON DFB LD WITH MAIN AND 1 SIDE MODES INJECTION-LOCKED TECHNIQUE

Demonstration of multi-cavity optoelectronic oscillators based on multicore fibers

High-resolution long-range distributed Brillouin analysis using dual-layer phase and amplitude coding

Improved slow-light performance of 10 Gb/s NRZ, PSBT and DPSK signals in fiber broadband SBS

Temporal phase mask encrypted optical steganography carried by amplified spontaneous emission noise

Setup of the four-wavelength Doppler lidar system with feedback controlled pulse shaping

Stable dual-wavelength oscillation of an erbium-doped fiber ring laser at room temperature

WDM Transmitter Based on Spectral Slicing of Similariton Spectrum

DIRECT MODULATION WITH SIDE-MODE INJECTION IN OPTICAL CATV TRANSPORT SYSTEMS

Photonic Signal Processing(PSP) of Microwave Signals

MICROWAVE photonics is an interdisciplinary area

Multiwatts narrow linewidth fiber Raman amplifiers

A new picosecond Laser pulse generation method.

Time-Multiplexed Pulse Shaping

Multi-wavelength laser generation with Bismuthbased Erbium-doped fiber

A WDM passive optical network enabling multicasting with color-free ONUs

CONTROLLABLE WAVELENGTH CHANNELS FOR MULTIWAVELENGTH BRILLOUIN BISMUTH/ERBIUM BAS-ED FIBER LASER

Dependence of stimulated Brillouin scattering in pulsed fiber amplifier on signal linewidth, pulse duration, and repetition rate

Novel OBI noise reduction technique by using similar-obi estimation in optical multiple access uplink

Chapter 1. Overview. 1.1 Introduction

All-Optical Clock Division Using Period-one Oscillation of Optically Injected Semiconductor Laser

Optical phase-locked loop for coherent transmission over 500 km using heterodyne detection with fiber lasers

Channel wavelength selectable singleõdualwavelength erbium-doped fiber ring laser

Advanced Optical Communications Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay

Novel High-Q Spectrum Sliced Photonic Microwave Transversal Filter Using Cascaded Fabry-Pérot Filters

All-optical clock division at 40 GHz using a semiconductor amplifier. nonlinear interferometer

Theory and Applications of Frequency Domain Laser Ultrasonics

Noise-based Brillouin optical correlation domain analysis with mm resolution

Fiberoptic Communication Systems By Dr. M H Zaidi. Optical Amplifiers

High-Power Highly Linear Photodiodes for High Dynamic Range LADARs

Modifying Bragg Grating Interrogation System and Studying Corresponding Problems

OPTICAL CORRELATION REFLECTOMETRY WITH SYNTHESIZED COHERENCE FUNCTION AND ITS EXTRACTION BY COHERENCE SYNCHRONIZATION

Single-longitudinal mode laser structure based on a very narrow filtering technique

Slow, Fast, and Backwards Light: Fundamental Aspects

Practical Aspects of Raman Amplifier

Interferometric Distributed Sensing System With Phase Optical Time-Domain Reflectometry

o Conclusion and future work. 2

Characteristics of point-focus Simultaneous Spatial and temporal Focusing (SSTF) as a two-photon excited fluorescence microscopy

Multiwavelength Single-Longitudinal-Mode Ytterbium-Doped Fiber Laser. Citation IEEE Photon. Technol. Lett., 2013, v. 25, p.

FIBER OPTICS. Prof. R.K. Shevgaonkar. Department of Electrical Engineering. Indian Institute of Technology, Bombay. Lecture: 37

PERFORMANCE OF PHOTODIGM S DBR SEMICONDUCTOR LASERS FOR PICOSECOND AND NANOSECOND PULSING APPLICATIONS

Novel RF Interrogation of a Fiber Bragg Grating Sensor Using Bidirectional Modulation of a Mach-Zehnder Electro-Optical Modulator

Extending the Offset Frequency Range of the D2-135 Offset Phase Lock Servo by Indirect Locking

Fast Raman Spectral Imaging Using Chirped Femtosecond Lasers

Photonic Generation of Millimeter-Wave Signals With Tunable Phase Shift

Yb-doped Mode-locked fiber laser based on NLPR Yan YOU

Experimental demonstration of both inverted and non-inverted wavelength conversion based on transient cross phase modulation of SOA

Faculty of Science, Art and Heritage, Universiti Tun Hussein Onn Malaysia, Batu Pahat, Johor, Malaysia.

1550 nm Programmable Picosecond Laser, PM

A Novel Method of Evaluating the Frequency Response of a Photoacoustic Cell

PERFORMANCE ANALYSIS OF WDM AND EDFA IN C-BAND FOR OPTICAL COMMUNICATION SYSTEM

CSO/CTB PERFORMANCE IMPROVEMENT BY USING FABRY-PEROT ETALON AT THE RECEIVING SITE

Coherent temporal imaging with analog timebandwidth

Photonics (OPTI 510R 2017) - Final exam. (May 8, 10:30am-12:30pm, R307)

Communication using Synchronization of Chaos in Semiconductor Lasers with optoelectronic feedback

FI..,. HEWLETT. High-Frequency Photodiode Characterization using a Filtered Intensity Noise Technique

Loop Mirror Multi-wavelength Brillouin Fiber Laser Utilizing Semiconductor Optical Amplifier and Fiber Bragg Grating

Coherent power combination of two Masteroscillator-power-amplifier. semiconductor lasers using optical phase lock loops

SBS based Slow-Light Performance Comparison of 10-Gb/s NRZ, PSBT and DPSK Signals

AN EFFICIENT L-BAND ERBIUM-DOPED FIBER AMPLIFIER WITH ZIRCONIA-YTTRIA-ALUMINUM CO-DOPED SILICA FIBER

Time-stretched sampling of a fast microwave waveform based on the repetitive use of a linearly chirped fiber Bragg grating in a dispersive loop

Chapter 8. Wavelength-Division Multiplexing (WDM) Part II: Amplifiers

Transcription:

Slope-assisted fast distributed sensing in optical fibers with arbitrary rillouin profile Yair Peled, * Avi Motil, Lior Yaron and Moshe Tur The Faculty of Engineering, Tel-Aviv University, Tel-Aviv 69978, Israel * yairpeled@gmail.com Abstract: We present a novel method, based on stimulated rillouin scattering (SS), for the simultaneous distributed measurement of fast strain variations along the entire length of the sensing fiber. A specially synthesized and adaptable probe wave is used to place the rillouin interaction always on the slope of the local rillouin gain spectrum, allowing a single pump pulse to sample fast strain variations along the full length of a fiber with an arbitrary distribution of the rillouin frequency shift. In this early demonstration of the method, strain vibrations of a few hundred Hz are demonstrated, simultaneously measured on two different sections of an 85m long fiber, having different static rillouin shifts and with a spatial resolution of 1.5m. 211 Optical Society of America OCIS codes: (6.237) Fiber optics sensors; (29.583) rillouin scattering, (33.188) detection; (19.19) nonlinear optics. References and links 1. M. Nikles, L. Thevenaz, and P. A. Robert, rillouin Gain Spectrum Characterization in Single-Mode Optical Fibers, IEEE J. Light. Technol. 15(1), 1842 1851 (1997). 2. A. W. rown,. G. Colpitts, and K. rown, Dark-Pulse rillouin Optical Time-Domain Sensor with 2-mm Spatial Resolution, IEEE J. Light. Technol. 25(1), 381 386 (27). 3. W. Li, X. ao, Y. Li, and L. Chen, Differential pulse-width pair OTDA for high spatial resolution sensing, Opt. Express 16(26), 21616 21625 (28). 4. S. M. Foaleng, M. Tur, J.-C. eugnot, and L. Thevenaz, High spatial and spectral resolution long-range sensing using brillouin echoes, IEEE J. Light. Tech. 28(2), 2993 33 (21). 5. K. Y. Song, Z. He, and K. Hotate, Distributed strain measurement with millimeter-order spatial resolution based on rillouin optical correlation domain analysis, Opt. Lett. 31(17), 2526 2528 (26). 6. Y. S. Kwang and K. Hotate, Distributed Fiber Strain Sensor with 1-kHz Sampling Rate ased on rillouin Optical Correlation Domain Analysis, IEEE Photon. Technol. Lett. 19(23), 1928 193 (27). 7. K. Y. Song, M. Kishi, Z. He, and K. Hotate, High-repetition-rate distributed rillouin sensor based on optical correlation-domain analysis with differential frequency modulation, Opt. Lett. 36(11), 262 264 (211). 8. L. Thévenaz, Inelastic Scatterings and Applications to Distributed Sensing in Advanced Fiber Optics - Concepts and Technology, Thévenaz L. ed, (Lausanne, Switzerland: EPFL Press, 211). 9. Z. Zhang and X. ao, Distributed optical fiber vibration sensor based on spectrum analysis of Polarization- OTDR system, Opt. Express 16(14), 124 1247 (28). 1. A. Zadok, E. Zilka, A. Eyal, L. Thévenaz, and M. Tur, Vector analysis of stimulated rillouin scattering amplification in standard single-mode fibers, Opt. Express 16(26), 21692 2177 (28). 11. A. Voskoboinik, J. Wang,. Shamee, R. S. Nuccio, L. Zhang, M. Chitgarha, E. A. Willner, and M. Tur, SS- ased Fiber Optical Sensing Using Frequency-Domain Simultaneous Tone Interrogation, IEEE J. Light.Technol. 29, 1729 1735 (211). 12. K. Hotate and S. S. L. Ong, Distributed fiber rillouin strain sensing by correlation-based continuous-wave technique ~cm-order spatial resolution and dynamic strain measurement, Proc. SPIE 492, 299 31 (22). 13. R. ernini, A. Minardo, and L. Zeni, Dynamic strain measurement in optical fibers by stimulated rillouin scattering, Opt. Lett. 34(17), 2613 2615 (29). 14. Y. Peled, A. Motil, L. Yaron, and M. Tur, Distributed and dynamical rillouin sensing in optical fibers, Proc. SPIE 7753, 775323, 775323-4 (211). 15. K. Shimizu, T. Horiguchi, and Y. Koyamada, Measurement of distributed strain and temperature in a branched optical fiber network by use of rillouin optical time-domain reflectometry, Opt. Lett. 2(5), 57 59 (1995). 16. X. ao, A. rown, M. Demerchant, and J. Smith, Characterization of the rillouin-loss spectrum of singlemode fibers by use of very short (<1-ns) pulses, Opt. Lett. 24(8), 51 512 (1999). (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19845

1. Introduction Fiber optic strain and temperature distributed sensors based on stimulated rillouin scattering (SS) have been extensively investigated in the last two decades. One of the more common approaches for implementing such a sensor, is the rillouin optical time domain analysis (OTDA). Here, a pump pulse wave, which is launched into one end of the sensing fiber, nonlinearly interacts with a counter propagating CW probe wave. This rillouin-based interaction, which involves the excitation of an acoustic field, strongly depends on the optical frequency difference between the pump and probe waves, /2, which is around 11GHz for standard single mode fibers at 155nm, with a narrow linewidth of ~3MHz. For a given fiber the exact value of the rillouin Frequency Shift (FS),, is sensitive to both the strain and temperature at the interaction position [1]. In the OTDA technique the optical frequency of either the pump or source waves is swept across 1-2MHz to determine along the sensing fiber, where the spatial resolution is determined by the pump pulse width. OTDA systems have been commercialized and proved to be very efficient for long range distributed sensing. Classical OTDA, however, has also its limitations. Due to the finite time required for the acoustic field to be excited by the interacting pump pulse and CW probe, the spatial resolution is limited to ~1m. Recently, several time domain techniques [2 4], as well as a correlation domain method (OCDA) [5] have been developed, improving the spatial resolution down to the order of cm and mm, respectively. Measurement speed is also a concern with OTDA: To achieve high strain/temperature resolution over a wide dynamic range of these two measurands, the scanned frequency range must be wide and of high granularity, resulting in a fairly slow procedure, which together with the need for averaging, limits the OTDA method to the quasi static domain. Using a correlation domain technique, 2 Hz distributed sensing (at 1kHz sampling rate) was demonstrated at a single fiber location, with 1-cm spatial resolution and 2-m measurement range [6]. A variant of the correlation technique [7] achieved strain distribution along the entire length of a 1-m fiber with 8-cm spatial resolution and 2-Hz sampling rate. While impressive, it has been argued [8] that the number of spatially resolvable points for the correlation technique is limited to ~6, which is an order of magnitude smaller than the number achieved using pulse based techniques. Another attempt to perform dynamic rillouin sensing [9] employs the dependence of the efficiency of SS on the relative states of polarization of the pump and probe waves [1]. Thus, the application of stress to a fiber segment will change the polarization states of the two interacting waves, thereby affecting the strength of the rillouin signal, although not in a way directly traceable to the magnitude of the applied strain. Still vibration frequencies of up to 5 khz were demonstrated in a 1km fiber link [9]. Another approach [11] proposes to use multiple pumps and multiple probes to avoid the time consuming frequency sweeping required by the classical OTDA technique. Measurement speed will potentially increase but at the expense of frequency granularity. It turns out, though, that a small modification of the classical OTDA setup can achieve very fast sensing, albeit with a limited dynamic range. Using the OCDA technique and working at a fixed pump-probe frequency difference on the slope of the rillouin gain spectrum (GS), Hotate and Ong [12] have measured 5Hz vibrations at 2kHz sampling rate (single fiber location). More recently, by tuning the probe frequency to the center of the rising/falling slopes of the rillouin gain spectrum (GS), Romeo et al. [13] utilized the SS interaction between two counter-propagating pump and probe pulses, meeting at a selected spatial location, to demonstrate fast strain-induced variations of the intensity of amplified probe wave with a sampling rate of 2Hz. The interrogated location was determined by the relative delay between the counter-propagating pulses. Clearly, at each interrogated location the optical frequency of the probe must be adjusted to properly sit at (or near) the center of the slope of the local GS, whose peak is likely to vary along the sensing fiber according to the local average strain/temperature and/or fiber properties. (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19846

In this paper we propose a new method for fast and distributed rillouin sensing, where the FS of the optical fiber may vary arbitrarily along the fiber length. The method is based on the classical OTDA technique, with the following modifications: instead of using a CW probe wave and a swept frequency pump pulse, it uses a pump pulse of a fixed optical frequency and a variable optical frequency CW probe wave. The time evolution of the probe frequency is designed in such a way that when the probe wave meets the counter-propagating pump pulse at location z along the fiber, the optical frequency difference between these two waves sits as close as possible to the middle of the slope of the GS, as in [12]. Any fast change in the local strain or temperature will shift the GS, and consequently, will be translated to gain variations of the probe wave. Using this technique the whole length of the fiber can be interrogated with a single pulse, or a few if averaging is required. Sec. 2 details the proposed method, first sketched in [14], while Sec. 3 describes the experimental setup and some results. We demonstrate the measurement of fast (4Hz) strain vibrations along 85m of a longitudinally rillouin-inhomogeneous fiber, with a spatial resolution of 1.5m. 2. Method Fig. 1. An example of the rillouin gain spectrum (GS) distance-frequency distribution along a sensing fiber, having 3 different GS sections with three different static 3d frequency values, d () z, of 1.9, 1.95 and 1.85GHz, respectively, see inset. For an optical fiber in a real environment, the FS, location along the fiber, can be described by:, as a function of both time and ( t, z) ( z) ( t, z) (1) where () z and ( tz, ), respectively, denote its static and dynamic components. y static we mean a temporally slow enough change that can be easily monitored by classical OTDA. Since this static (or averaged, in case of fast strain/temperature changes) FS has a non-uniform spatial distribution along the fiber length, we first employ classical OTDA (with a fixed pump frequency and swept probe frequency) to plot the static rillouin Gain Spectrum (GS), as a function of both the fiber length coordinate z and the frequency difference between the pump and probe waves. An example appears in Fig. 1, showing a fiber comprising three uniform sections, each with a different FS. Searching for the 3d contour along the fiber (say on the falling slope of the GS), we deduce the range dependence of, see inset in Fig. 2. Here, z is the point where the pump pulse enters the fiber, d () z while z L denotes the fiber end, which is also the entrance point for the probe wave. We then temporally tailor the frequency of the probe waveform so that when the pump pulse arrives at fiber location z, it meets a probe wave, whose optical frequency is exactly (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19847

below the pump frequency, ensuring rillouin interaction around the 3d gain d () z point of the GS at location z. To satisfy the requirements of the proposed method, we now consider a probe, having the following waveform: E ( t) A exp j[2 t ( t)], (2) Probe where Pump is the optical frequency of the pump, () t is defined by d g Pump t ( t) 2 V t '/ 2 dt ', (3) and V g is the fiber group velocity. ( t is chosen to be the time when the probe wave enters the fiber at z L). The propagation of the probe waveform through the fiber at t and z L obeys: E ( t, z) A exp j[2 ( t [ z L] / V ) ( t [ z L] / V )], t [ L z] / V, (4) Probe Pump g g g Therefore, the instantaneous frequency at t, z t [ L z]/ Vg is: d Probe ( t, z) Pump Pump d Vg ( t [ z L] / Vg ) / 2 dt [ V t z L] / 2 Pump d g Thus, at t the probe enters the fiber at z L with (), and it is not until / g Probe Pump d t L V that the probe waveform arrives at z, still having the frequency Probe Pump d (). At that instant ( t L / Vg ), a pump pulse, E ( t) A( t)exp{ j[2 t]}, (6) Pump whose amplitude, A() t, defined as centered at t, is launched at the fiber entrance, z. Its propagation through the fiber is governed by: Pump E ( t, z) A t [ z L]/ V exp j[2 ( t [ z L]/ V )] (7) Pump g Pump g An arbitrary fiber section at z z, characterized by its own FS, ( z d ), is reached by the pump pulse at t [ z L]/ Vg. At these values of t and z the probe waveform has a frequency of (Eq. (5)): Probe Pump d g g Pump d (5) ( t, z ) [ V [ z L]/ V z L]/ 2 ( z ) (8) which is the probe frequency ensuring mid-slope gain for the probe. Figure 2 describes the temporal evolution of the probe waveform tailored for the fiber of Fig. 1, as well as spatial snapshots of the probe and pump waves, counter-propagating through the fiber. Clearly, the waveform comprises 3 segments, each with its designed optical frequency. As dictated by the form of Eq. (5), each probe segment has twice the length of the corresponding fiber section. (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19848

Fig. 2. An example of a pump pulse, propagating against a complex probe wave, which comprises 3 different optical frequencies, corresponding to the fiber of Fig. 1. Note that each probe segment has twice the length of the corresponding fiber section, as per Eq. (5). Dynamic strain changes at z z, spectrally shift the local GS to lower or higher frequencies around its average value, ( z d ) (e.g., by 5MHz/1µS for standard single mode fibers at 155nm [15]). Thus, while meeting the pump pulse at z z and t t, where t [ z L]/ Vg, the probe at frequency Probe Pump d ( z ) will no longer experience the 3d rillouin gain; instead the gain will be lower or higher, depending on the direction of the GS shift. Accordingly, the recorded probe intensity, arriving at the fiber edge ( z ) at t [2 z L]/ Vg will reflect this gain change, which was induced by the strain at t t, z z. The magnitude of the strain-to-gain conversion factor depends not only on the strain-to-frequency sensitivity but also on the slope of the previously measured Lorentianshaped GS. As mentioned above, this linewidth is around 3MHz for standard single mode fibers at 155nm, provided the pump pulse width significantly exceeds the acoustic life time (~1ns). Unless special coding of the pump wave is employed [2,4], a shorter pump pulse gives rise to a broadened GS, resulting in a more gradual slope. Much like the classic OTDA method, each pump pulse generates strain information from the entire fiber. Since there is no need to sweep the probe frequency, the sampling rate of the strain changes is limited only by the fiber length and the need for averaging to a value 1 bounded by 2 N L / V N is the number of averages. avg g, where avg While this proposed technique enables a very fast true distributed sensing over a fiber with an arbitrary distribution of its static rillouin shift, its dynamic range for dynamic strain measurements is limited. To ensure linearity of the strain-to-gain conversion, the dynamic (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19849

strain variations to be measured must be constrained to generate frequency shifts (of the GS) smaller than the frequency span of the linear part of the GS slope (<3MHz for long pump pulses). Using the 5MHz/1µS strain-to-frequency conversion factor, the resulting dynamic range is ~6μє [13]. Using short pump pulses to broaden the GS [16] will benefit the dynamic range at the expense of smaller strain-to-gain conversion factor. Finally, slow temporal variations of () can be followed either by evaluating the d z average of the strain fluctuations coming from distance z, or by intermittent application of classical OTDA. Using the updated value for (), the frequency composition of the d z probe waveform can be appropriately readjusted using Eqs. (2-3). Fig. 3. Experimental setup: AWG: arbitrary waveform generator, EOM: electro-optic modulator, EDFA: Erbium-doped fiber amplifier, CIR: circulator, FG: fiber ragg grating, PS: polarization scrambler, IS: isolator, ATT: attenuator, FUT: fiber under test, PD: photodiode. Fig. 4. The 85m FUT, comprising five sections of SMF fiber. The two 1m sections are mounted on manually stretching stages, making it possible to adjust their static rillouin frequency shifts. Additionally, audio speakers are physically attached to these two sections in order to induce fast strain variations of various frequencies and magnitudes. All sections, with the exception of the two patch cords are made of the same fiber. 3. Experiment and results A highly coherent 155nm DF laser diode (DF-LD), with a linewidth of 1kHz, is split into pump and probe channels, Fig. 3. A complex waveform, to be described below, feeds the probe channel Mach-Zehnder modulator (EOM1), which is biased at its zero transmission point to generate two sidebands, the lower one for the probe wave and the upper one to be discarded later by the fiber ragg grating (FG) filter. The EOM1 output is then amplified by an Erbium doped fiber amplifier (EDFA1), optionally scrambled by a polarization scrambler (PS), and launched into one side of the fiber under test (FUT), Fig. 4, through the attenuator (ATT). Modulator EOM2 forms the pump pulse, which is then amplified by EDFA2 and (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 1985

launched into the other side of the FUT through a circulator (CIR1). The rillouin-amplified probe wave is finally routed to a fast photodiode (PD) by CIR1 and CIR2. A narrow bandwidth fiber ragg grating (FG) filters out pump backscattering, as well as the upper sideband generated by EOM1. Finally, the output of the photodiode is sampled at 1GSamples/s by a real-time oscilloscope with deep memory. Here d () z is a staircase function, representing the mid-slope frequency in each of the FUT sections. Since () is on the order of 11GHz, it is not easy to directly synthesize d z the probe waveform of Eqs. (2-3). Instead, a dual channel, 5MHz arbitrary waveform generator (AWG), together with a -2GHz microwave vector signal generator with I/Q inputs, are used to generate the RF input to EOM1. Thus, the signal generator frequency, f c, is set to a fixed value around, but below min[ d ( z )] (~11GHz), and we rewrite Eq. (3) as: t t ( t) 2 V t '/ 2 dt ' 2 f t 2 V t '/ 2 f dt ' d g c d g c (9) and denote the second term on the right-hand side of (9) by () t. Using this technique, the AWG 2 ft c term in (9) takes care of a fixed but high frequency component of () t, while () t is responsible for the variable part of (). Now the two channels of the AWG, AWG VI () t and VQ () t are programmed to output: d z I AWG Q AWG V ( t) V cos ( t) ; V ( t) V sin ( t) (1) These two signals of Eq. (1) are then connected to the I/Q inputs of the vector signal generator, whose output becomes: V ( t) V ( t)cos(2 f t) V ( t)sin(2 f t) V cos(2 f t ( t)) V cos( ( t)) RF I c Q c c AWG t (11) V cos 2 d Vgt '/ 2 dt ' iased at zero transmission, the EOM1 Mach-Zehnder modulator stops the pump frequency and using V () t of Eq. (11), generates a lower sideband of the form RF t EProbe ( t) EEOM1 ( t) exp j 2 Pumpt 2 d Vg t '/ 2 dt ', (12) which is exactly the desired complex probe wave, as defined by Eqs. (2-3). In all experiments reported below, 15ns wide pump pulses were used at a repetition rate of 625kHz. Once the pump pulse enters the 85m fiber, we need to collect data for 85ns ( 85m / 2 ) in order to monitor the full length of the fiber. Therefore, each pump pulse V g generates 85 (= 1GHz 85ns) recorded samples of the intensity of the rillouin-amplified probe wave, as measured at z. The sampled data were arranged in a matrix of N rows by M columns, where M 85 is the number of oscilloscope samples per pump pulse, and N is the number of pump pulse cycles used in the measurement. To maximize the frequency resolution of the measured strain variations, data were continuously collected until the scope memory was effectively exhausted, resulting in N =31,25 and a total recording time of 5ms (=N/625kHz), and, consequently, an expected frequency resolution of 2Hz. Thus, each column of the data matrix represents the time evolution of the strain at a particular location along the FUT, where the first and last columns, respectively report the gain at the beginning ( z ) and end ( z L) of the FUT. (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19851

We first demonstrate the distributed measurement of fast vibrations along fiber sections having the same FS. For that purpose, the two 1m sections at the end of the 85m fiber, Fig. 4, were identically stretched, ensuring their rillouin gain spectra overlap but away from the rest of the loose parts of the fiber. Here, EOM1 was simply fed by a single frequency RF sine signal, chosen to generate a probe wave, whose frequency coincided with the 3d point of the common GS. Using the audio speakers, the two sections were vibrated at different frequencies of 15Hz and 4Hz. Figure 5 presents the measured vibrations as a function of time along 1m out of the 85m FUT. Time sequences of two columns of the above mentioned measured data matrix, corresponding to the centers of the two 1m sections, also appear in the figure. Here gain readings were converted to frequency values using the measured slope of the GS, as found from a classical OTDA mapping of the FUT (shown below for the next experiment, Figs. 6-7), utilizing the same pump pulse and a swept CW probe wave. The obtained frequency variations can then be translated to strain values, and using the conversion factor of 5MHz/1µε and the very little observed noise (.25MHz rms, as measured at a location with no induced vibrations), quite high strain sensitivity is to be expected. Fig. 5. Top: Strain-induced gain vibrations at 15Hz and 4Hz were measured at the two 1m fiber sections, when adjusted to have the same GS. ottom: Time sequences from two columns of the measured data matrix, corresponding to the centers of the first and second sections, after 1kHz low pass filtering. Here gain variations were converted to frequency values. Note that each segment was excited by a different speaker and no effort was made to induce vibrations of the same magnitude in the two sections. To demonstrate fast rillouin sensing over a rillouin non-uniform fiber, the two 1m sections, I and II of Fig. 4, were stretched to different strains, resulting in the rillouin distance-frequency distribution of Figs. 6-7, as determined by a classical OTDA measurement using the same 15ns pump pulse. Clearly, the rillouin frequency shifts for sections I and II and for the 4m section between them are 1.9, 11.2, and 1.84GHz, respectively. Ignoring the non-vibrating 4m section, we chose to work on the falling slope of section I, (I) =1.93GHz and on the rising slope of section II, (II) =1.99GHz d (Note that working on the negative slope results in signal inversion of the measured data). In order to measure both sections with a single pump pulse, while ignoring the rest of the fiber, a composite probe waveform was synthesized, comprising a long segment, having an optical frequency falling outside the GS of all sections of Fig. 6, followed by a 5ns segment, downshifted from the pump by 1.93GHz for section I, and another 5ns segment, downshifted by 1.99GHz for section II. The pump pulse and the probe waveform were temporally synchronized so that the center of the pump pulse meets the probe frequency d (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19852

change (from 1.93GHz to 1.99GHz) at the middle of the 4m section of Fig. 4. This synchronization ensures that in the absence of vibrations, the pump pulse, propagating in either section I or II, exclusively see the exact 3d probe frequency of that section. Note that had we strictly followed the recipe of the method (section 2 above) for the 1m-4m-1m sections, the probe waveform would have had comprised of the following segments: 1ns(1.93GHz)-4ns(1.87GHz)-1ns(1.99GHz) (assuming we work on the falling slope of the GS of the 4m section). Instead, in this first demonstration of the method, we opted to concentrate on the two 1m vibrating sections and ignore the 4m section, as well as all other loose sections Still, for a theoretical zero-width pump pulse, this could be accomplished by a 3ns(1.93GHz)-3ns(1.99GHz) probe, which provides proper coverage for the 1m-4m-1m sections. However, to maximize the interaction of our finite-width, 15ns pump pulse (~2ns at its base) with the probe, we chose the longer 5ns (1.93GHz)-5ns (1.99GHz) probe segments. Fig. 6. Classical OTDA was employed to measure the static rillouin gain spectrum (GS) distance-frequency distribution along the FUT of Fig. 4 when sections I and II were nonuniformly stretched. A 15ns pump pulse was used against a CW probe wave, whose frequency was sequentially swept to cover the relevant frequency range. Fig. 7. A blow-up of the sections of interest appears this 3D drawing, indicating that the FWHM of the GS is of ~6MHz, which is due to the use of a relatively short pump pulse of 15ns. (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19853

Applying 18Hz to section I and 32Hz to section II, we obtained the results shown in Fig. 8, clearly confirming the feasibility of this new technique. The bump at the center of the 4m section (distance coordinate: 79m) is an artifact originating from the interaction of the pump pulse with the frequency transition of the probe wave from 1.93GHz to 1.99GHz that occurs exactly at that point, as dictated by our synchronization. This fast transition contains frequency components in the vicinity of the GS of the 4m section, resulting in the observed bump, which, as expected, shows no vibrations. Luckily, in a real-world scenario, sharp transitions of the GS are not expected, except where two different types of fiber are spliced together. Fig. 8. Top: Strain-induced gain vibrations at 18Hz and 32Hz were measured at the two 1m fiber sections, having different GS. ottom: Time sequences from two columns of the measured data matrix, corresponding to the centers of the first and second sections, after 1kHz low pass filtering. Here too, gain variations were converted to frequency values. Like in Fig. 5, the difference in magnitude between the 18Hz and 32Hz vibrations is due to different excitation conditions. 4. Discussion and summary In this paper we proposed and demonstrated a novel method for SS-based high-frequency distributed strain measurement for a fiber whose rillouin frequency shift varies along its length. Since no frequency sweeping of either the pump or the probe is required, the upper limit on the measurement speed is set by the pump pulse repetition rate, as determined by the fiber length and the number of required averages. Operating on the slope of the GS [13], the dynamic range of the allowable vibrations is limited by the frequency extent of this slope (~3MHz for not-too-short pump pulses), and can be effectively enlarged by employing pump pulses shorter than the acoustic life time [4]. The probe waveform can be adaptively modified to follow slowly changing static strain/temperature conditions. The proposed method is fully compatible with modern techniques for the enhancement of spatial resolution [2 4]. Further research is underway to study the various characteristics of the technique, as well as its performance in real-world scenarios. Acknowledgements This research was supported in part by the Israel Science Foundation. (C) 211 OSA 1 October 211 / Vol. 19, No. 21 / OPTICS EXPRESS 19854