Theoretical and Experimental Investigation of the Effect of Pump Laser Frequency Fluctuations on Signal-to-Noise Ratio of Brillouin Dynamic Grating Measurement with Coherent FMCW Reflectometry

Signal-dependent speckle-like noise has constituted a serious factor in Brillouin-grating based frequency-modulated continuous-wave (FMCW) reflectometry and it has been indispensable for improving the signal-to-noise ratio (S/N) of the Brillouin dynamic grating measurement to clarify the noise generation mechanism. In this paper we show theoretically and experimentally that the noise is generated by the frequency fluctuations of the pump light from a laser diode (LD). We could increase the S/N from 36 to 190 merely by driving the LD using a current source with reduced technical noise. On the basis of our experimental result, we derived the theoretical formula for S/N as a function of distance, which contained the second and fourth-order moments of the frequency fluctuations, by assuming that the pump light frequency was modulated by the technical noise. We calculated S/N along the 1.35 m long optical fiber numerically using the measured power spectral density of the frequency fluctuations, and the resulting distributions agreed with the measured values in the 10 to 190 range. Since higher performance levels are required if the pump light source is to maintain the S/N as the fiber length increases, we can use the formula to calculate the light source specifications including the spectral width and rms value of the frequency fluctuations to achieve a high S/N while testing a fiber of a given length.


Introduction
The Brillouin-enhanced four-wave mixing induced by counter-propagating pump lights and one probe light produces backward Stokes light at every location in an optical fiber under test [1], which is assumed to be the reflection of the probe light by the acoustic wave or Brillouin dynamic grating generated in the fiber by the two pump lights. Optical time-domain reflectometry (OTDR) has been used to detect the power of the Stokes light as a function of distance while changing the frequency difference between the pump lights to obtain the Brillouin spectrum distribution in the optical fiber [2][3][4][5][6]. A micrometer-scale spatial resolution is indispensable for diagnosing planar lightwave circuits (PLCs) [7] with the same four-wave mixing technique. Since the spatial resolution of the OTDR is determined by the temporal width of the employed optical pulses, a picosecond optical pulse should be launched into the optical fiber and detect return pulses without deformation by using a high-speed optical detector. Such an attempt to increase the spatial resolution often resulted in increased electrical noise due to the ultra-wide detection bandwidth over 10 GHz, and this degraded the signal-to-noise ratio (S/N). To our knowledge, therefore, there have been no reports on Brillouin dynamic grating measurement with a micrometer-scale spatial resolution using the OTDR method. To overcome the S/N degradation, we have proposed Brillouin grating-based coherent frequency-modulated continuous-wave (FMCW) reflectometry [8], which is the frequencydomain counterpart of optical low coherence reflectometry (OLCR) [9] and it detects the power of the Stokes light by utilizing its interference with local oscillator (LO) light at a detection bandwidth of less than 1 MHz. We have succeeded in incorporating a fiber loop mechanism in conventional coherent FMCW reflectometry to generate the Stokes light [10][11][12]. Since the distance to the fiber segment where the Stokes light is generated is derived from the beat frequency, we do not have to use a mechanical stage such as that used in OLCR and thus we can extend the available distance range much further than can be achieved by the translation of the stage. Hereafter in this paper, the distribution of the Stokes light power along a fiber is referred to simply as a reflectogram from the fiber. Once we acquire reflectograms at various frequency differences, the Brillouin spectrum distribution along the optical fiber is obtained by changing the horizontal grid from equal distances to equal frequency intervals in the two-dimensional power data, which provides us with the strain distribution of the fiber.
Although we succeeded in demonstrating reflection measurement from a Brillouin dynamic grating with coherent FMCW reflectometry, we found signal-dependent multiplicative noise [13] to be the dominant noise, which meant that we could not improve the S/N merely by increasing the powers of the tunable laser output and pump lights and/or by narrowing the detection bandwidth. Hereafter in this paper, for convenience such multiplicative noise is referred to as speckle-like noise. Since the reflectograms that we derived from different measurements or made at different times had unavoidable speckle-like noise, we had no choice but to repetitively sweep the tunable laser source and add them to obtain a smoothed profile at every frequency difference. That is, we could obtain the desired Brillouin spectrum distribution and thus the strain distribution only after making a vast number of the repetitive sweeps over a long period of time. Therefore, to reduce the number of sweeps and the measurement time it was indispensable that we clarify the origin of the speckle-like noise and reduce it greatly. This paper shows theoretically and experimentally that the speckle-like noise was generated by frequency fluctuations contained in the output light wave from the distributed feedback laser diode (DFB LD) used as the pump light source. Since the current from the employed current source injected into the LD had components that fluctuated with time, or technical noise, the resultant instantaneous frequency of the light output also fluctuated, and this generated the speckle-like noise in the reflectogram. First, we describe the experimental setup we employed for the coherent FMCW reflectometry system. Since silica-based PLC chips are usually between few millimeters and several tens of centimeters long, we adjusted the length of the optical fiber under test to be around 1 m. We constructed an unbalanced Mach-Zehnder interferometer (MZI) to obtain the power spectral density of the frequency fluctuations [14]. We employed two commercially available current sources from different manufacturers to drive the same LD, which had different power spectral density distributions. We used the density data to calculate the theoretical S/N and compared it with the experimental result.
Next, we theoretically derive the S/N dependence on the power spectral density of the frequency fluctuations on the condition that the absolute value of the complex amplitude of the acoustic wave or Brillouin dynamic grating was constant along the fiber. We calculated the second and fourth-order terms of the variance of the fluctuating Stokes light signal, which depended on the moments of the instantaneous frequency deviation of the pump light wave, and we derived the latter term as the correction of the former. Then, we compare the theoretical results with the experimental data obtained from a 1.35 m long optical fiber to confirm that the speckle-like noise was generated by the frequency fluctuations of the pump light waves. We showed that the second-order term provided us with the correct values for the S/N ranging from 10 to 190 when the coherence effect of the pump light wave was negligible. Then we measured the S/N distributions for 10 and 40 m long optical fibers. As the fiber length approached the coherence length of the pump light, Sensors 2021, 21, 2870 3 of 33 the S/N values calculated with both the second and fourth-order terms no longer agreed with the measured value. Finally, we introduced the coherence function of the pump light wave into the second and fourth-order terms of the variance to explain the coherent effect and compared the calculated results with the data.

Experimental Setup for the Reflectometer
The Brillouin grating-based coherent FMCW reflectometer comprised the conventional coherent FMCW reflectometry setup [10][11][12] for detecting the reflection from a device under test (DUT) and an optical fiber loop [2][3][4][5][6] for generating an acoustic wave or Brillouin dynamic grating in the DUT by counter-propagating pump lights, as shown in Figure 1. The conventional FMCW part was designed to detect the reflection whose optical frequency was down-converted by the Brillouin dynamic grating induced in the DUT. To achieve this detection, we incorporated a single-sideband modulator (SSBM) in the LO arm to down-convert the LO light. The power of the probe light launched into the DUT was 60 mW. The photocurrent output from the balanced mixer was converted to a voltage with a transimpedance amplifier (TIA) to obtain the beat signal waveform produced by interference between the Stokes light and the LO light. We drove the LiNbO 3 phase modulators of PM1 and PM3 at f 0 = 150 kHz and f 1 = 190 kHz, respectively, with sawtooth voltage waveforms supplied by a two-channel function generator so that the carrier frequency of the beat signal waveform was f 1 − f 0 = 40 kHz. sors 2021, 21, x FOR PEER REVIEW 3 of 33 pump light wave was negligible. Then we measured the S/N distributions for 10 and 40 m long optical fibers. As the fiber length approached the coherence length of the pump light, the S/N values calculated with both the second and fourth-order terms no longer agreed with the measured value. Finally, we introduced the coherence function of the pump light wave into the second and fourth-order terms of the variance to explain the coherent effect and compared the calculated results with the data.

Experimental Setup for the Reflectometer
The Brillouin grating-based coherent FMCW reflectometer comprised the conventional coherent FMCW reflectometry setup [10][11][12] for detecting the reflection from a device under test (DUT) and an optical fiber loop [2][3][4][5][6] for generating an acoustic wave or Brillouin dynamic grating in the DUT by counter-propagating pump lights, as shown in Figure 1. The conventional FMCW part was designed to detect the reflection whose optical frequency was down-converted by the Brillouin dynamic grating induced in the DUT. To achieve this detection, we incorporated a single-sideband modulator (SSBM) in the LO arm to down-convert the LO light. The power of the probe light launched into the DUT was 60 mW. The photocurrent output from the balanced mixer was converted to a voltage with a transimpedance amplifier (TIA) to obtain the beat signal waveform produced by interference between the Stokes light and the LO light. We drove the LiNbO3 phase modulators of PM1 and PM3 at f0 = 150 kHz and f1 = 190 kHz, respectively, with sawtooth voltage waveforms supplied by a two-channel function generator so that the carrier frequency of the beat signal waveform was f1 − f0 = 40 kHz.  [15], which consisted of the conventional coherent FMCW reflectometry setup for detecting the reflection from a device under test (DUT) and an optical fiber loop for generating a Brillouin dynamic grating in the DUT. DFB LD: Distributed feedback laser diode, CL1 and CL2: Optical circulators, CP1, CP2 and CP3: 3-dB optical fiber couplers, EDFA: Erbium-doped fiber amplifier, PC1 and PC2: Polarization controllers, PBS1 and PBS2: Polarization beam splitters, SSBM: Single-sideband modulator (T.SBXH1.5-20PD-ADC, Sumitomo Osaka Cement), FG: 2-channel function generator, PM1, PM2 and PM3: LiNbO3 phase modulators (LN65S-FC, Thorlabs), LO: Local oscillator, TIA: Transimpedance amplifier, ωp: Pump light frequency, Ω: Up-conversion frequency of the pump light which is equal to the down-conversion frequency of the LO light, f0 = 150 kHz, f1 = 190 kHz. Picosecond pulses were launched into the fiber loop to locate the position at zc where the optical path lengths of the counter-propagating pump light waves were equal.
Although not shown in Figure 1, some of the tunable laser output was supplied to an auxiliary unbalanced fiber-optic MZI to change the grid of the interference fringe produced during the frequency sweep from equal time to equal frequency increments. The time delay between the two arms of the MZI was adjusted to 58.449 ns. We sampled the beat signal waveform from the balanced mixer twice when the reference beat signal from the MZI traveled one cycle. The DUT, which is shown in blue in Figure 1, consisted of four Although not shown in Figure 1, some of the tunable laser output was supplied to an auxiliary unbalanced fiber-optic MZI to change the grid of the interference fringe produced during the frequency sweep from equal time to equal frequency increments. The time delay between the two arms of the MZI was adjusted to 58.449 ns. We sampled the beat signal waveform from the balanced mixer twice when the reference beat signal from the MZI traveled one cycle. The DUT, which is shown in blue in Figure 1, consisted of four pieces of single-mode fiber. Two of them were fiber pigtails attached to optical circulators CL1 and CL2 and the rest were fiber patch cords with angled physical contact (APC) connectors. We spliced each pigtail to either patch cord using fusion splicing and mated the connectors together with an SC/APC adapter. The total length of the DUT was L = 1.35 m and the bases of the fiber pigtails at the CL1 and CL2 were denoted as the input and output ends of the DUT, respectively. We used a configuration consisting of a pair of polarization controllers (PC1 and PC2) and a polarization beam splitter (PBS2) to block the pump light from entering the balanced mixer [15].
In the optical fiber loop, we used a DFB LD operating at 1550 nm as the pump light source, which was driven by two commercially available current sources that we refer to as current sources A and B. We used current source A for the reflection measurements described in our previous papers [8,15], whereas we prepared current source B specifically for use in this experiment. The LD was packaged in an aluminum case where a D-sub connector was attached to supply injection current to the LD and control current to a thermoelectric cooler installed in the LD module. We selected either current source as the LD driver by connecting the D-sub connector at the LD case with that of the current source via RS232 cable.
The LD output was divided into two with an optical fiber coupler (CP3) for use as two counter-propagating pump lights. We describe the fixed frequency of the laser output as ω p in the figure. To up-convert the frequency of the pump light by the same frequency Ω as the down-conversion frequency of the probe light, we applied phase modulation at Ω to the pump light with a LiNbO 3 phase modulator (PM2) and extracted the up-converted light component with an optical narrow-band filter. The resultant pump light at ω p + Ω was amplified with an erbium-doped fiber amplifier (EDFA), combined with the probe light at a polarization beam splitter (PBS1), and launched into the DUT after passing through the optical circulator CL1, where the launched power was 120 mW. The other laser output was launched into the DUT from the opposite direction for use as a pump light at ω p after passing through the optical circulator CL2.
We drove the pump LD with the current source A and acquired 30 reflectograms from the 1.35 m long fiber by repeatedly sweeping the tunable laser, as shown in Figure 2a, where the up-conversion frequency was 10.861 GHz. The optical fiber had two spliced points and one mated point where the Brillouin frequency shifts had been moved downward so that the Stokes signals at those points decreased rapidly. With current source A, the signal changed greatly with every sweep, and we had to obtain a smooth profile by adding individual reflectograms. We then drove the same LD with another current source, B, while maintaining the same measurement condition and found that the signal variations were greatly reduced, as shown in Figure 2b. By comparing these results, we considered that the noise contained in the current injected into the LD produced the fluctuations in the optical frequency of the pump light, which resulted in the signal variations. We placed the optical fiber under test straight on an optical bench in a normal laboratory environment, but the Stokes signal changed gradually along the fiber. This was because the optical fiber was not a polarization-maintaining fiber and the states of polarization (SOPs) of the pump lights and the probe light changed along the fiber, resulting in changes in both the amplitude and SOP of the Stokes light.

Basic Formulation Using Instantaneous Frequency Shift
In Figure 3, the light propagation of the two pump light waves in the fiber loop is shown with green dotted lines, and the optical fiber as the DUT is shown in blue. The origin of the distance z is chosen at the point where the optical path lengths of the LO light and the reflection from the point are equal. We assume that the fiber distributed from the input end at z = zi to the output end at z = ze so that the total length was ze − zi =L = 1.35 m. The acoustic field at time t and distance z excited by material density fluctuations in the DUT is described as ρ(z,t) = ρ0 + {R(z)exp[i(qz − Ωt)] + c.c.} [1], where ρ0 is mean density, q is the wavenumber of the acoustic field, Ω is the frequency difference between the counter-propagating pump light waves, i is pure imaginary and c.c. denotes the complex conjugate. Figure 3. Schematic of the propagation (highlighted in green) of counter-propagating pump light waves from the LD in the fiber loop. t5 and t6 are the propagation times of the pump light waves from the LD to optical circulators CL1 and CL2, respectively. z is a coordinate of the distance along the fiber under test highlighted in blue and the origin of the distance is positioned at the point where the produced reflection has the same optical path length as the LO light at the balanced mixer. The input and output ends of the fiber were assumed to be located at zi and ze, respectively. τ is the round-trip time from the origin to any position at z as defined by τ = 2nz/c, where n is the refractive index of the fiber and c is the velocity of light in a vacuum. Similarly, τi and τe are defined as τi = 2nzi/c and τe = 2nze/c. Ap1 and Ap2 are the complex amplitudes of the electric fields of the pump light waves propagating in the clockwise and counterclockwise directions, respectively. The center of pumping is defined by the position where the optical path lengths of the two pump light waves are equal in the optical fiber loop. CL1 and CL2: Optical circulators, CP3: Fiber coupler, LD: Laser diode.

Basic Formulation Using Instantaneous Frequency Shift
In Figure 3, the light propagation of the two pump light waves in the fiber loop is shown with green dotted lines, and the optical fiber as the DUT is shown in blue. The origin of the distance z is chosen at the point where the optical path lengths of the LO light and the reflection from the point are equal. We assume that the fiber distributed from the input end at z = z i to the output end at z = z e so that the total length was z e − z i = L = 1.35 m. The acoustic field at time t and distance z excited by material density fluctuations in the DUT is described as ρ(z,t) = ρ 0 + {R(z)exp[i(qz − Ωt)] + c.c.} [1], where ρ 0 is mean density, q is the wavenumber of the acoustic field, Ω is the frequency difference between the counterpropagating pump light waves, i is pure imaginary and c.c. denotes the complex conjugate.

Basic Formulation Using Instantaneous Frequency Shift
In Figure 3, the light propagation of the two pump light waves in the fiber loop is shown with green dotted lines, and the optical fiber as the DUT is shown in blue. The origin of the distance z is chosen at the point where the optical path lengths of the LO light and the reflection from the point are equal. We assume that the fiber distributed from the input end at z = zi to the output end at z = ze so that the total length was ze − zi =L = 1.35 m. The acoustic field at time t and distance z excited by material density fluctuations in the DUT is described as ρ(z,t) = ρ0 + {R(z)exp[i(qz − Ωt)] + c.c.} [1], where ρ0 is mean density, q is the wavenumber of the acoustic field, Ω is the frequency difference between the counter-propagating pump light waves, i is pure imaginary and c.c. denotes the complex conjugate.  Under a steady-state condition, the amplitude of the acoustic wave is denoted as: where ε 0 is the vacuum permittivity, γ e is the electrostrictive constant of the DUT material, Ω B is the center frequency of the Brillouin spectrum as a function of z, which is to be measured with the reflectometer, and Γ B is the Brillouin linewidth. When the amplitude of the electric field of the LD output at t is denoted as A(t), the amplitudes of the counterpropagating pump light waves at (t,z) are represented by where c 1 and c 2 are constants, and t 5 and t 6 are the times required for these light waves to propagate from the LD to the optical circulators CL1 and CL2, respectively. It should be noted that we had already introduced parameters t 1 t 2 , t 3 and t 4 as the propagation times in the reflectometer setup in our previous paper [8] and thus in this paper we introduce new parameters t 5 and t 6 to avoid confusion. The phases of both waves always agree with each other at a point z = z c = L/2 + c(t 6 − t 5 )/2n, which is hereafter referred to as the center of pumping. Here we introduce τ and τ c as the round-trip times from the origin to a point at z and to the center of pumping at z c , which are given by τ = 2nz/c and τ c = 2nz c /c, respectively. By shifting the origin of time t by t 5 , the amplitudes of the counter-propagating pump light waves at (t, z) are denoted as If the phase of the output light wave from the LD is modulated by the technical noise contained in the injection current, we can denote its amplitude as A(t) = A 0 exp(−iφ(t)) with a constant A 0 , and thus A p1 A p2 * in Equation (1)

is written as
That is, the absolute value of the complex amplitude of the acoustic wave is unchanged, but its temporal phase is modulated by the technical noise. When we make an approximation where , the approximation error is of the order of the product of the root-mean-square value of φ (t) and the maximum value of |(τ − τ c )τ c |/2.
Since the phase induced by the current noise is generally represented by φ(t) = ∑φ k cos2 πf k t, we have φ (t) = −2π∑δν k sin2πf k t and thus φ (t) = −(2π) 2 ∑δν k f k cos2πf k t, where δν k is the maximum frequency deviation at f k . Since φ (t)/2π is the instantaneous frequency deviation of the pump light wave from the center frequency, its mean square value, which is denoted as δν rms 2 , is calculated to be δν rms 2 = ∑δν k 2 /2. Similarly, we have <(φ (t)) 2 > = (2π) 4 ∑δν k 2 f k 2 /2, which should be less than (2π) 4 f u 2 ∑δν k 2 /2 = (2π) 4 f u 2 δν rms 2 , where f u is the upper limit of the frequency contributing to the speckle-like noise. The result shows that we can estimate the root mean square value of φ (t) as √ < (φ (t)) 2 > ≈ (2π) 2 f u δν rms . On the other hand, the factor |(τ−τ c )τ c |/2 is of the order of τ e 2 /8 because τ e ≈ 0 and thus τ c ≈ τ e /2 and 0 < τ < τ e . Therefore, we find that the approximation error is of the order of { √ < (φ (t)) 2 >}τ e 2 /8 = π 2 f u δν rms τ e 2 /2. In our experiment we adopt τ e = 13.5 ns, f u ≈ 1 kHz, and δν rms ≈ 2 MHz for the current source A, as will be described in Section 4, and thus the approximation error is estimated to be 1.8 × 10 −6 rad, meaning that the accuracy of the approximation is very high and the error is negligibly small.
Since φ (t) represents the instantaneous angular frequency deviation from ω p , we denote it as ∆(t) so that we have a good approximation of . We have already derived the expression for the power of the Stokes light as a function of τ [8]. With this approximation, the power is obtained by calculating the absolute square of the Fourier inverse transform of the following analytic signal I(t) with respect to ω(t): where the proportionality coefficient of the signal is included in r(τ) for simplicity so that r(τ) is the complex function of τ which is independent of the phase modulation due to technical noise. ω 1 and ω 1 + ω(t) are the start frequency and the instantaneous frequency at time t of the tunable laser output, respectively. The phase of the analytic signal changed with time due to ∆(t), resulted in the fluctuations of the Stokes light signal.

Experimental Setup for Measuring Power Spectral Density
We define the signal-to-noise ratio or S/N of the reflection measurement by the ratio of the mean level of the fluctuating powers of the Stokes light signal to their standard deviation, which is defined by the square root of the variance. Since the variance depends on the power spectral density of the angular frequency fluctuations ∆(t), we constructed an unbalanced fiber-optic MZI as shown in Figure 4 to measure the power spectral density [14,[16][17][18][19]. Since we did not have two more LiNbO 3 phase modulators for frequency up-conversion, we installed acousto-optic frequency modulators (AOM1 and AOM2) in the two arms of the interferometer and drove them so that the frequency difference between the up-conversion frequencies in both arms became the same at 40 kHz. By connecting a fiber patch cord as a fiber delay line to one arm of the interferometer, we adjusted the time delay between the two arms to τ MZI = 8.912 ns.
Sensors 2021, 21, x FOR PEER REVIEW 7 of 33 at time t of the tunable laser output, respectively. The phase of the analytic signal changed with time due to Δ(t), resulted in the fluctuations of the Stokes light signal.

Experimental Setup for Measuring Power Spectral Density
We define the signal-to-noise ratio or S/N of the reflection measurement by the ratio of the mean level of the fluctuating powers of the Stokes light signal to their standard deviation, which is defined by the square root of the variance. Since the variance depends on the power spectral density of the angular frequency fluctuations Δ(t), we constructed an unbalanced fiber-optic MZI as shown in Figure 4 to measure the power spectral density [14,[16][17][18][19]. Since we did not have two more LiNbO3 phase modulators for frequency upconversion, we installed acousto-optic frequency modulators (AOM1 and AOM2) in the two arms of the interferometer and drove them so that the frequency difference between the up-conversion frequencies in both arms became the same at 40 kHz. By connecting a fiber patch cord as a fiber delay line to one arm of the interferometer, we adjusted the time delay between the two arms to τMZI = 8.912 ns. Schematic of an unbalanced Mach-Zehnder interferometer (MZI) used to measure frequency fluctuations of the light that was incident from the input port of the MZI. AOM1 and AOM2: Acousto-optic frequency shifters, FG: Function generator, PC3: Polarization controller, TIA: Transimpedance amplifier. f2 and f3 are up-conversion frequencies by the AMO1 and AOM2, respectively. PC3 was adjusted for the amplitude of the beat signal to reach its maximum value. The inset shows a typical phase change waveform, which was retrieved from the acquired beat signal waveform.
We launched the output from the DFB LD into the MZI and sampled the beat signal from the balanced mixer at a rate of 256 k samples/s for 11.7 s. We calculated the Fourier transform of the acquired waveform from which we removed the carrier frequency and calculated the Fourier inverse transform of the result to obtain the analytic signal. Then we unwrapped the phase term of the analytic signal and obtained the temporal change of the phase, as shown by the inset in Figure 4. The overall profile of the retrieved phase change tended to decay gradually and then remained constant mainly due to environmental perturbations applied to the MZI during the measurement. Considering that the delay time τMZI was shorter than τe, the phase of the analytic signal is given by Δ(t)τMZI, and so we could obtain the actual change of Δ(t) by dividing the measured phase change by the delay τMZI. It is noted that the temporal change of the light frequency, δν(t), is given by Δ(t)/2π.
We calculated the power spectral density H(f) of the temporal frequency change by applying Fourier transformation to δν(t) which was given by Δ(t)/2π, as shown in Figure  5a,b, where the LD was driven with current sources A and B, respectively, and the sampling interval was 0.0853 Hz. For comparison, in each figure we plotted the spectral density, which we obtained from a low noise continuous wave diode-pumped solid state (DPSS) laser operating at 1.34 μm. Each spectral density from the LD contained a vast number of peaks that were superimposed on the gradually decaying background component. Since all the components ranging from 1 Hz to 10 kHz were larger than the background noise level of the spectral density of the DPSS laser, we concluded that they all originated from the frequency fluctuations of the LD itself, not from noise in the detection system. Transimpedance amplifier. f 2 and f 3 are up-conversion frequencies by the AMO1 and AOM2, respectively. PC3 was adjusted for the amplitude of the beat signal to reach its maximum value. The inset shows a typical phase change waveform, which was retrieved from the acquired beat signal waveform.
We launched the output from the DFB LD into the MZI and sampled the beat signal from the balanced mixer at a rate of 256 k samples/s for 11.7 s. We calculated the Fourier transform of the acquired waveform from which we removed the carrier frequency and calculated the Fourier inverse transform of the result to obtain the analytic signal. Then we unwrapped the phase term of the analytic signal and obtained the temporal change of the phase, as shown by the inset in Figure 4. The overall profile of the retrieved phase change tended to decay gradually and then remained constant mainly due to environmental perturbations applied to the MZI during the measurement. Considering that the delay time τ MZI was shorter than τ e , the phase of the analytic signal is given by ∆(t)τ MZI , and so we could obtain the actual change of ∆(t) by dividing the measured phase change by the delay τ MZI . It is noted that the temporal change of the light frequency, δν(t), is given by ∆(t)/2π.
We calculated the power spectral density H(f ) of the temporal frequency change by applying Fourier transformation to δν(t) which was given by ∆(t)/2π, as shown in Figure 5a,b, where the LD was driven with current sources A and B, respectively, and the sampling interval was 0.0853 Hz. For comparison, in each figure we plotted the spectral density, which we obtained from a low noise continuous wave diode-pumped solid state (DPSS) laser operating at 1.34 µm. Each spectral density from the LD contained a vast number of peaks that were superimposed on the gradually decaying background component. Since all the components ranging from 1 Hz to 10 kHz were larger than the background noise level of the spectral density of the DPSS laser, we concluded that they all Over a 100 Hz to 1 kHz frequency range, we observed more intense peaks in the spectral density (a) than in (b), which we considered to be the main origin of the larger fluctuations in the reflectograms in (a) than in (b) of Figure 2. As observed in the inset in Figure 4, the phase changed gradually for 11.7 s mainly due to external perturbations applied to the interferometer, and thus we concluded that the spectral components ranging from 0.1 to 1 Hz in both Figure 5a,b arose not from the actual frequency fluctuations, but from the environmental perturbations.

Calculation
In this section we derive the mean level and standard deviation of the fluctuating power of the Stokes light signal to obtain the S/N of the reflection measurement. We assume that the frequency of the tunable laser output sweeps from the start frequency ω1 sufficiently linearly with time to meet ω(t) = βt with a constant β. To simplify the calculation of the Fourier inverse transform of Equation (2) with respect to ω(t), we introduce a new variable ῶ, which is defined by ω1 + ω(t) − ωp = ῶ, so that t is a function of ῶ represented by t =(ῶ + ωp − ω1)/β. By expanding the exponential function in the integrand of I(t) in power series up to the fourth order, I(t) is approximated as: The Fourier inverse transform of I(t) with respect to ῶ is then written as: where X (0) = r(τ) and: and where t1 = (ῶ1 + ωp − ω1)/β. Hereafter we abbreviate r(τ) as r when there is no confusion.
is independent of the random process, whereas X (k) (k ≥ 1) is a random variable that contains the kth power of Δ(t) in the integrand. From Equation (4) the absolute square is given by: which provides the mean value of Z as: Over a 100 Hz to 1 kHz frequency range, we observed more intense peaks in the spectral density (a) than in (b), which we considered to be the main origin of the larger fluctuations in the reflectograms in (a) than in (b) of Figure 2. As observed in the inset in Figure 4, the phase changed gradually for 11.7 s mainly due to external perturbations applied to the interferometer, and thus we concluded that the spectral components ranging from 0.1 to 1 Hz in both Figure 5a,b arose not from the actual frequency fluctuations, but from the environmental perturbations.

Calculation
In this section we derive the mean level and standard deviation of the fluctuating power of the Stokes light signal to obtain the S/N of the reflection measurement. We assume that the frequency of the tunable laser output sweeps from the start frequency ω 1 sufficiently linearly with time to meet ω(t) = βt with a constant β. To simplify the calculation of the Fourier inverse transform of Equation (2) with respect to ω(t), we introduce a new variableω, which is defined by ω 1 + ω(t) − ω p =ω, so that t is a function ofω represented by t =(ω + ω p − ω 1 )/β. By expanding the exponential function in the integrand of I(t) in power series up to the fourth order, I(t) is approximated as: The Fourier inverse transform of I(t) with respect toω is then written as: where X (0) = r(τ) and: and where t 1 = (ω 1 + ω p − ω 1 )/β. Hereafter we abbreviate r(τ) as r when there is no confusion.
is independent of the random process, whereas X (k) (k ≥ 1) is a random variable that contains the kth power of ∆(t) in the integrand. From Equation (4) the absolute square is given by: which provides the mean value of Z as: when we decompose Z into <Z> +δZ, the variance σ 2 of Z is <(δZ) 2 >, which is represented by: The calculations for deriving Equations (7) and (8) are detailed in Appendix A, where the summands in Equation (8) are defined by Equations (A5) to (A10). The superscript k in Y (k) denotes that Y (k) contains the kth power of ∆(t) in the integral.
In accordance with the Wiener-Khinchin theorem [20,21], we introduce the power spectral density G(χ) of ∆(t), which is defined by: where G(χ) should be a real-valued and even function of the angular frequency χ. By letting ∆(t) = 2πδν(t) and χ = 2πf, Equation (9) is changed to <δν(t 1 )δν( is the instantaneous frequency deviation of the pump light and f is frequency, and thus we find the relational expression of G(χ) = 2πH(χ/2π) between G(χ) and the power spectral density H(f ) of δν(t). It is noted that each plot in Figure 5 shows H(f ) measured when we drove the LD with either current source A or current source B. The mean square value of δν(t), which we denote as δν rms 2 , is represented by the infinite integral of H(f ) with respect to f. Actually, however, the detection bandwidth for the Stokes light signal was limited to f u = βτ e /2π so that δν rms 2 must be evaluated by taking the finite integral over (−f u , f u ), which is changed to the integral over (0, f u ) as: because H(f ) is an even function of f. Equation (10) means that δν rms depends on the length of the fiber under test and is related to ∆ rms , which is the root mean square value of ∆(t) via ∆ rms = 2πδν rms . To express <Z> using G(χ), we must calculate the following double integral at k = 2 and 4: By assuming that ∆(t) is a zero-mean Gaussian random variable, we obtain <∆ 4 (t)> = 3∆ rms 4 with ∆ rms 2 = <∆ 2 (t)> according to the moment theorem [21][22][23]. In the double integral in Equation (11), first we take the iterated integral with respect toω 1 to obtain 2πδ(τ − τ 1 ), where δ(τ) is the delta function. Then we integrate the resultant integrand including the delta function with respect to τ 1 to obtain <X (2) > = − r(τ)∆ rms 2 (τ c − τ) 2 /2 and <X (4) > = r(τ)∆ rms 4 (τ c − τ) 4 /8. When we neglect the last two terms <|X (1) | 2 > + <|X (2) | 2 > in Equation (7), we have <Z> ≈ |r(τ) which is approximated by the Gaussian function as: It is noted that the actual third terms of the power series expansion of the Gaussian function should be ∆ rms ) 2 > + higher-order terms, where the second term vanishes because the integrand of it contains the third-order moment of ∆(t). Then the variance is approximated by adding the second-order term σ (2) and fourth-order term σ (4) as: where: σ (4) is considered to be a correction term when we estimate the variance from Equation (14). We calculated σ (2) and σ (4) by assuming that ∆(t) obeys a Gaussian random process as described in Appendices B and C, respectively. Because |σ (4) | < < σ (2) as long as we measure the 1.35 m long optical fiber as will be shown in Section 4, here we describe the result of σ (2) only: In general, the variance depends on the distribution r(τ), which changes along the fiber in an actual field application, and this means that it would be rather difficult to calculate the integrals in Equation (16). Since our aim in this paper is to clarify the origin of the speckle-like noise that we observed in the reflection measurement, we assume that the fiber under test is uniform and tension-free throughout the fiber and thus r(τ) is a constant function of τ, or r(τ) = r 0 in the range τ i ≤ τ≤ τ e and r(τ) = 0 elsewhere. The integrands of the first and second terms in Equation (16) contain the functions of r(τ − χ/β) and r(τ − χ/β)r(τ + χ/β), which are finite only when the variable χ satisfies the conditions of τ i ≤ τ−χ/β ≤ τ e , and τ i ≤ τ ± χ/β ≤ τ e , respectively. Then <Z> and √ σ (2) are calculated to be proportional to |r 0 | 2 , which is dropped when taking the ratio of <Z> to √ σ (2) to obtain S/N. For this reason, we redefined σ (2) as the one which is obtained by letting |r 0 | = 1 in the original expression. The explicit forms of σ (2) for τ i < τ < (τ i + τ e )/2 and (τ i + τ e )/2 < τ < τ e are given by Equations (A42) and (A43), respectively, in Appendix D. The S/N including the fourth-order term σ (4) is calculated in Appendix E.
By introducing the variable u and parameters u c and u i , which are defined as u = τ/τ e , u c = τ c /τ e and u i = τ i /τ e , the ranges of τ for τ i < τ < (τ i + τ e )/2 and (τ i + τ e )/2 < τ < τ e are changed to those of u for u i < u < (u i + 1)/2, and (u i + 1)/2 < u < 1, respectively. Then the mean signal level described by Equation (12) is rewritten as: By letting |r(τ)| = 1 and thus S/N for u i < u < (u i + 1)/2 is described with the variance as: where: and where η = χ/(βτ e ) = f /f u , f u = βτ e /2π and G(βτ e η) = 2πH(f). To obtain the expression for (u i + 1)/2 < u < 1, in Equation (19) we should change the interval of the first integral to (1 − u, u − u i ) and that of the second integral to (0, 1 − u), and change the sign of the variable η in the integrand of the first term.
Since the power spectral density H(f ) of the frequency fluctuations had a vast number of peaks due to the technical noise from the current sources, as shown in Figure 5, to evaluate ∆ rms and σ (2) we should calculate Equation (19) numerically using the measured power spectral density. In Appendix F, we calculated σ (2) analytically on the assumption that the noise was concentrated on a narrow region expressed by a Gaussian function G(χ) = G 0 exp[− (χ/δχ e ) 2 ] using δχ e for the spectral half width at 1/e maximum and that the fiber under test was sufficiently long for the fiber input end and the center of pumping to be τ i ≈ 0 and τ c ≈ τ e /2, respectively, so that we can let u i = 0 and u c = 1/2. S/N at an arbitrary value of α = βτ e /δχ e is obtained by substituting σ (2) represented by Equation (A78) into Equation (18) for 0 < u < 1/2, or by substituting that represented by Equation (A80) into Equation (18) for 1/2 < u < 1.
Although the analytical solution to σ (2) had a complicated form, we could approximate it as represented by Equation (A83) and it is valid that: (20) at α ≥ α min = 8 within a calculation error of around 20%, as will be shown numerically in Section 4. In Equation (20), f h and γ, respectively, are the spectral half width at a half maximum of the Gaussian spectrum defined by f h = δχ e ( √ ln2)/2π, and the sweep rate of the tunable laser output defined by γ = β/2π. α is denoted as ( √ ln2)γτ e /f h , which should be less than α min so that the approximation of Equation (20) holds. To meet the condition, the range allowed for f h should be: Because 0 < u < 1, the numerator in Equation (20) takes γ at u = u c = 1/2 and decreases to γexp[− (∆ rms τ e /2) 2 ] at both ends, which is written as γexp[− (πδν rms τ e ) 2 ] using the relational expression ∆ rms = 2πδν rms . Thus, we impose the second condition: to regard the numerator as the constant γ so that the mean signal level is kept constant throughout the fiber. Then, (S/N) ∞ is simplified to √ (ln2/2)γ/(2πf h δν rms ), which should be equal to or greater than the target value of (S/N) tgt , and so we found that δν rms and f h must satisfy the third condition: at given values of γ and (S/N) tgt .

Experimental Results
We fixed the up-conversion frequency of the pump light wave at 10.861 GHz, which was the center frequency of the Brillouin spectrum of the 1.35 m long optical fiber under test. The parameters u i and u c were originally defined by u i = τ i /τ e and u c = τ c /τ e and were rewritten as u i = z i /z e and u c = z c /z e because of the linearity between the propagation time and distance, or τ i = 2nz i /c, τ c = 2nz c /c and τ e = 2nz e /c. We calculated the mean reflectogram from the 30 reflectograms shown in Figure 2a and measured the distances to the points where the Stokes signal raised and fell to give z i = 4.11 cm and z e = 139.8 cm.
With these values, we obtained u i = 0.0294. To locate the center of pumping, we launched a picosecond optical pulse from another input port of fiber coupler CP3 into the optical fiber loop and measured the propagation times during which the counter-propagating optical pulses reached the two angled-polished end faces of the optical fiber under test. Since the center of pumping was defined by the position where the optical path lengths of the counter-propagating pump lights were equal in the fiber loop, the distance z c to the center of pumping could be determined by the difference between the propagation times of the pulses. Since the distance increased by connecting a longer optical fiber delay line between polarizer #2 and optical circulator CL2, we precisely increased the length of the fiber delay line to shift the center of pumping to u c = −0.042, 0.45 and 0.93 in sequence while measuring the time differences between the counter-propagating optical pulses.
We drove the DFB LD with current source A. At each u c value, we swept the tunable laser source at a rate of 0.5 nm/s 30 times and derived the magnitude of the Stokes light signal as a function of z. We calculated the mean value and standard deviation at every distance from the 30 reflectograms and plotted the distribution of the S/N in Figure 6 after changing the grids of the horizontal axis from z to u, which was defined as u = z/z e . We drove the DFB LD with current source A. At each uc value, we swept the tunable laser source at a rate of 0.5 nm/s 30 times and derived the magnitude of the Stokes light signal as a function of z. We calculated the mean value and standard deviation at every distance from the 30 reflectograms and plotted the distribution of the S/N in Figure 6 after changing the grids of the horizontal axis from z to u, which was defined as u = z/ze. Assuming that statistically independent random noise was superimposed on the Stokes light signal, we could increase the S/N by √N by calculating the mean value of N individual signals [24]. Therefore, we should acquire at least 30 beat signal waveforms by sweeping the tunable laser 30 times and calculate the mean reflectogram of the individual reflectograms to achieve an S/N of the order of 200 at every up-conversion frequency. We measured the S/N as a function of u when the DFB LD was driven with current source B. The results are shown in Figure 7a-c when we set uc at −0.042, 0.45 and 0.93, respectively. By comparing each pair in Figures 6 and 7, e.g., Figures 6a and 7a, it was clear that we could greatly increase the S/N merely by changing the current source from A to B. For example, the best S/N was increased to 190 at u = 0.18 and uc = 0.45. This meant that a single sweep was sufficient to obtain a smooth reflectogram, which could be obtained only after sweeping the laser 30 times when we used current source A. For comparison with the measured values, we calculated S/N as a function of u numerically using Equations (18) and (19). The theoretical S/N was given by the ratio of the (2)  Assuming that statistically independent random noise was superimposed on the Stokes light signal, we could increase the S/N by √ N by calculating the mean value of N individual signals [24]. Therefore, we should acquire at least 30 beat signal waveforms by sweeping the tunable laser 30 times and calculate the mean reflectogram of the individual reflectograms to achieve an S/N of the order of 200 at every up-conversion frequency. We measured the S/N as a function of u when the DFB LD was driven with current source B.
The results are shown in Figure 7a-c when we set u c at −0.042, 0.45 and 0.93, respectively. By comparing each pair in Figures 6 and 7, e.g., Figures 6a and 7a, it was clear that we could greatly increase the S/N merely by changing the current source from A to B. For example, the best S/N was increased to 190 at u = 0.18 and u c = 0.45. This meant that a single sweep was sufficient to obtain a smooth reflectogram, which could be obtained only after sweeping the laser 30 times when we used current source A. We performed the numerical integration to obtain the value of δνrms by approximating the integral of H(f) with respect to f as the sum of the areas HkΔf of the equally-spaced rectangles under the curve H(f) for all values of k satisfying fL < fk < fu, where the value of H(f) taken at fk is denoted as Hk. In addition, we introduced the lower limit of integration at fL = 1 Hz to avoid the effect of the environmental perturbations. Following the numerical calculation, we obtained δνrms = 2.17 MHz and therefore we had Δrms = 1.36×10 7 rad/s from the relational expression of Δrms = 2πδνrms. With the values of τe and Δrms, we had Δrmsτe = 0.184 and thus we found that [Δrmsτe(uc − u)] 2 < 0.0368 for all values of u and uc. The result meant that we could consider the Gaussian function as unity when we tested the 1.35 m long optical fiber.
We calculated σ (2) as a function of u ranging from 0 to 1 in steps of 0.01, where βτe 3 = 9.66 × 10 −13 s, ui = 0.0294 and the value of G(χ) at χ was obtained from that of 2πH(f) at f = χ/2π. That is, at each value of u for ui < u < (ui + 1)/2 and ui = 0.0294, we approximated the first integral in Equation (19) as the sum of 2π(uc -u − ηk) 2 HkΔη for all values of k satisfying max{fL,(u − ui)fu} < fk < (1 − u)fu, where ηk = fk/fu, Δη = Δf/fu and max{ , } returns the greater of two values. Similarly, the second integral in Equation (19) was approximated as the sum of 2πηk 2 HkΔη for all values of k satisfying fL < fk < (u − ui)fu. With the two numerical integrations we calculated the value of σ (2) and the resultant S/N according to Equation (18), where the numerator was approximated as unity. At each value of u for (ui + 1)/2 < u < 1, we performed the same numerical calculation of the integrals in Equation (A45) in Appendix D to obtain σ (2) and thus S/N as a function of u.
By repeating these series of numerical integrations at uc = −0.042, 0.45 and 0.93, we obtained S/N as a function of u as shown in Figure 6a-c, respectively. Similarly, we calculated S/N as a function of u when the DFB LD was driven with current source B. The results are shown in Figure 7a-c when we set uc at uc = −0.042, 0.45 and 0.93, respectively. The calculation showed that when we set uc at 0.45 or near the center of the optical fiber under test, the distribution of S/N had peaks on both sides of the center, whereas one or both peaks were suppressed at uc = −0.042 and 0.93, which was consistent with the measurement result obtained using either current source A or current source B. We should set uc close to 0.5 to avoid the suppression of the peaks and achieve a higher S/N throughout the fiber. It is noted that the value of uc automatically approaches 0.5 as the length of the optical fiber under test is increased to 10 m or more. This is because the distance to the For comparison with the measured values, we calculated S/N as a function of u numerically using Equations (18) and (19). The theoretical S/N was given by the ratio of the Gaussian function to the square root of the variance σ (2) , where the Gaussian function could be regarded as unity when we tested the 1.35 m long optical fiber. This was concluded by calculating the value of δν rms as follows. δν rms was given by integrating H(f ) over the interval (0, f u ) with respect to f as described in Equation (10). Here we had already measured the power spectral density H(f ) as shown in Figure 5a, where the sampling interval was ∆f = 0.0853 Hz and each sampling frequency was represented by f k = k∆f (k = 0,1,2, . . . ). The upper limit of integration was f u , which was the detection bandwidth and given by f u = τ e β/2π. We obtained τ e = 13.5 ns by substituting z e = 139.8 cm and n = 1.45 into the relational expression of τ e = 2nz e /c. Since we set the sweep speed of the tunable laser source at 0.5 nm/s at 1.55 µm, we obtained β/2π = 62.5 GHz/s. With these values we found that the upper limit of integration was f u = 844 Hz.
We performed the numerical integration to obtain the value of δν rms by approximating the integral of H(f ) with respect to f as the sum of the areas H k ∆f of the equally-spaced rectangles under the curve H(f ) for all values of k satisfying f L < f k < f u , where the value of H(f ) taken at f k is denoted as H k . In addition, we introduced the lower limit of integration at f L = 1 Hz to avoid the effect of the environmental perturbations. Following the numerical calculation, we obtained δν rms = 2.17 MHz and therefore we had ∆ rms = 1.36 × 10 7 rad/s from the relational expression of ∆ rms = 2πδν rms . With the values of τ e and ∆ rms , we had ∆ rms τ e = 0.184 and thus we found that [∆ rms τ e (u c − u)] 2 < 0.0368 for all values of u and u c . The result meant that we could consider the Gaussian function as unity when we tested the 1.35 m long optical fiber.
We calculated σ (2) as a function of u ranging from 0 to 1 in steps of 0.01, where βτ e 3 = 9.66 × 10 −13 s, u i = 0.0294 and the value of G(χ) at χ was obtained from that of 2πH(f ) at f = χ/2π. That is, at each value of u for u i < u < (u i + 1)/2 and u i = 0.0294, we approximated the first integral in Equation (19) as the sum of 2π(u c -u − η k ) 2 H k ∆η for all values of k satisfying max{f L ,(u − u i )f u } < f k < (1 − u)f u , where η k = f k /f u , ∆η = ∆f /f u and max{,} returns the greater of two values. Similarly, the second integral in Equation (19) was approximated as the sum of 2πη k 2 H k ∆η for all values of k satisfying f L < f k < (u − u i )f u . With the two numerical integrations we calculated the value of σ (2) and the resultant S/N according to Equation (18), where the numerator was approximated as unity. At each value of u for (u i + 1)/2 < u < 1, we performed the same numerical calculation of the integrals in Equation (A45) in Appendix D to obtain σ (2) and thus S/N as a function of u.
By repeating these series of numerical integrations at u c = −0.042, 0.45 and 0.93, we obtained S/N as a function of u as shown in Figure 6a-c, respectively. Similarly, we calculated S/N as a function of u when the DFB LD was driven with current source B. The results are shown in Figure 7a-c when we set u c at u c = −0.042, 0.45 and 0.93, respectively. The calculation showed that when we set u c at 0.45 or near the center of the optical fiber under test, the distribution of S/N had peaks on both sides of the center, whereas one or both peaks were suppressed at u c = −0.042 and 0.93, which was consistent with the measurement result obtained using either current source A or current source B. We should set u c close to 0.5 to avoid the suppression of the peaks and achieve a higher S/N throughout the fiber. It is noted that the value of u c automatically approaches 0.5 as the length of the optical fiber under test is increased to 10 m or more. This is because the distance to the center of the optical fiber is much longer than the distance to the original center of pumping which is located if we connect the input end directly to the output end of the fiber.
We plotted the difference between the calculated and measured S/N values in percent as a function of u in Figure 8, where (a) and (b) show the distributions calculated from the data shown in Figures 6b and 7b, respectively. The optical fiber under test consisted of four pieces of the single-mode fiber which were fusion spliced at u = 0.31 and 0.78 and mated with APC connectors at 0.53, where the Stokes light signal decreased and thus the deviation changed rapidly. Excluding these parts, however, the differences were limited to within 20% in most parts of the fiber as highlighted in yellow in each figure. We derived the theoretical expression for S/N by imposing a condition whereby the absolute value of the complex amplitude of the acoustic wave was constant along the fiber, whereas its temporal phase was modulated by the frequency fluctuations of the pump light. In spite of the fact that we could not experimentally generate the Brillouin dynamic grating uniformly throughout the fiber, the calculated values were in good agreement with those measured in the 10 to 190 range as shown in Figures 6 and 7. Therefore, we confirmed that the theoretical expression we derived could represent the actual S/N of the reflection measurement and that the main source which substantially determined the S/N in our experiment was technical noise from the current source used to drive the DFB LD. Although S/N fell rapidly to 70 at both ends of the fiber as shown in Figure 7b even when we used current source B, we will be able to increase it to the same level as the peak by using a different current source with much reduced technical noise.
Since the calculated values did not exactly agree with the measured ones obtained with either current source A or current source B, we calculated the correction term σ (4) when using current source A at u c = 0.45 by numerically calculating the constituent terms A through H, which are defined by Equations (A59) to (A66), while using the power spectral density shown in Figure 5a. The calculated values of σ (2) and σ (4) as a function of u are plotted in Figure 9 together with the change in each term. σ (2) had a minimum value of 10 −3 on both sides of the center of pumping, whereas σ (4) had a maximum value of 3 × 10 −5 at 0.1 < u < 0.9. Thus, we found that the correction term was at best 3% of σ (2) , and we could not bring the calculated values very much closer to the measured ones. The theoretical expressions for σ (2) and σ (4) were based on the expansion of the exponential function in the integrand of I(t) in a power series of ∆(t). The result showed that σ (4) was negligibly smaller than σ (2) , and the first-order expansion was sufficient to calculate S/N as far as the reflection measurement of the 1.35 m long optical fiber was concerned. measured in the 10 to 190 range as shown in Figures 6 and 7. Therefore, we confirmed that the theoretical expression we derived could represent the actual S/N of the reflection measurement and that the main source which substantially determined the S/N in our experiment was technical noise from the current source used to drive the DFB LD. Although S/N fell rapidly to 70 at both ends of the fiber as shown in Figure 7b even when we used current source B, we will be able to increase it to the same level as the peak by using a different current source with much reduced technical noise. Since the calculated values did not exactly agree with the measured ones obtained with either current source A or current source B, we calculated the correction term σ (4) when using current source A at uc = 0.45 by numerically calculating the constituent terms A through H, which are defined by Equations (A59) to (A66), while using the power spectral density shown in Figure 5a. The calculated values of σ (2) and σ (4) as a function of u are plotted in Figure 9 together with the change in each term. σ (2) had a minimum value of 10 −3 on both sides of the center of pumping, whereas σ (4) had a maximum value of 3 × 10 −5 at 0.1 < u < 0.9. Thus, we found that the correction term was at best 3% of σ (2) , and we could not bring the calculated values very much closer to the measured ones. The theoretical expressions for σ (2) and σ (4) were based on the expansion of the exponential function in the integrand of I(t) in a power series of Δ(t). The result showed that σ (4) was negligibly smaller than σ (2) , and the first-order expansion was sufficient to calculate S/N as far as the reflection measurement of the 1.35 m long optical fiber was concerned. . Numerical calculation of the constituent terms A to H of the correction term σ (4) together with σ (2) when current source A was used, where the center of pumping was located at uc = 0.45. We calculated σ (4) to improve the discrepancy between the measured and calculated data observed in Figure 6b.
Since the phases of the counter-propagating pump light waves always coincide with each other at the center of pumping, the temporal phase of the acoustic wave generated at the center is independent of the frequency fluctuations of the pump light waves. Considering that the phase of the Stokes light that was generated there was also unaffected by the fluctuations, it appears that the variance due to the frequency fluctuations become negligibly small and the resultant S/N should be at its maximum exactly at the center. Actually, however, when the center was located at uc = 0.45, the S/N value was not at its maximum at the center and had noticeable peaks on both sides. The reflectometry method we employed decomposed the beat signal waveform into components with different frequencies and distributed the powers of the decomposed components along the fiber. When the frequency distribution of the beat signal waveform overlapped that of Δ(t), as in our experiment, the Stokes lights generated on either side of the center had components with the same frequency as that produced at the center via up and down frequency conversions by the phase modulation, which were superimposed on and recognized as the Stokes light generated at the center. This meant that even the signal detected at the center had larger fluctuations than we expected.
Planar lightwave circuits (PLCs) such as an 8 × 8 optical matrix switch and an arrayed-waveguide grating [7] have fiber pigtails several meters long, which are connected to their input and output ends with adhesive for practical applications. To investigate the possibility of using our reflectometer to diagnose them by means of strain distribution, we tested a 5 m long polarization-maintaining (PM) fiber. To excite the fast and slow axes of the PM fiber with the probe light and the counter-propagating pump lights, respectively, we spliced 250 μm buffered single-mode fibers with APC connectors to both ends of the PM fiber with a mechanical fusion splicer and mated the APC connectors to those of the fiber pigtails, which were attached to optical circulators CL1 and CL2. We inserted Figure 9. Numerical calculation of the constituent terms A to H of the correction term σ (4) together with σ (2) when current source A was used, where the center of pumping was located at u c = 0.45. We calculated σ (4) to improve the discrepancy between the measured and calculated data observed in Figure 6b.
Since the phases of the counter-propagating pump light waves always coincide with each other at the center of pumping, the temporal phase of the acoustic wave generated at the center is independent of the frequency fluctuations of the pump light waves. Considering that the phase of the Stokes light that was generated there was also unaffected by the fluctuations, it appears that the variance due to the frequency fluctuations become negligibly small and the resultant S/N should be at its maximum exactly at the center. Actually, however, when the center was located at u c = 0.45, the S/N value was not at its maximum at the center and had noticeable peaks on both sides. The reflectometry method we employed decomposed the beat signal waveform into components with different frequencies and distributed the powers of the decomposed components along the fiber. When the frequency distribution of the beat signal waveform overlapped that of ∆(t), as in our experiment, the Stokes lights generated on either side of the center had components with the same frequency as that produced at the center via up and down frequency conversions by the phase modulation, which were superimposed on and recognized as the Stokes light generated at the center. This meant that even the signal detected at the center had larger fluctuations than we expected.
Planar lightwave circuits (PLCs) such as an 8 × 8 optical matrix switch and an arrayedwaveguide grating [7] have fiber pigtails several meters long, which are connected to their input and output ends with adhesive for practical applications. To investigate the possibility of using our reflectometer to diagnose them by means of strain distribution, we tested a 5 m long polarization-maintaining (PM) fiber. To excite the fast and slow axes of the PM fiber with the probe light and the counter-propagating pump lights, respectively, we spliced 250 µm buffered single-mode fibers with APC connectors to both ends of the PM fiber with a mechanical fusion splicer and mated the APC connectors to those of the fiber pigtails, which were attached to optical circulators CL1 and CL2. We inserted each 250 µm buffered single-mode fiber into an in-line polarization controller with a rotatable fiber squeezer mechanism to adjust the SOPs of the probe and pump lights propagating through the fiber.
After setting the up-conversion frequency at 10.881 GHz and changing to current source B, we measured 30 reflectograms from the 5 m long PM fiber and calculated the S/N distribution along the fiber as shown in Figure 10. As a result of the measurement range being extended from 1 to 5 m, the S/N was degraded to the 20 to 45 range. Considering that except around the center of pumping at z c = 3.8 m the measured and calculated values agreed, we concluded that we needed to use a current source with greatly reduced noise to increase the S/N. Although not shown in this paper, we observed a background component in the beat signal waveform that undulated periodically during the frequency sweep and which we did not observe when we tested the 1.35 m long single-mode fiber. We changed to another PM fiber with a different cutoff wavelength, but we still observed such a component. Therefore, we believe that some fraction of the launched probe and pump lights propagated without being attenuated through the fiber in the cladding modes such as the modes in the stress-applying parts or between it and the core. There was a probability that the periodic change in the S/N that we observed around the center was produced by the light component, which propagated in the cladding modes and entered the balanced mixer to produce the fluctuated components by interference between it and the LO light. such a component. Therefore, we believe that some fraction of the launched probe and pump lights propagated without being attenuated through the fiber in the cladding modes such as the modes in the stress-applying parts or between it and the core. There was a probability that the periodic change in the S/N that we observed around the center was produced by the light component, which propagated in the cladding modes and entered the balanced mixer to produce the fluctuated components by interference between it and the LO light. We investigated the specifications of the pump light wave required for achieving the S/N at 200 even when we tested a 5 m long PM fiber. When the power spectral density was a Gaussian function, we obtained the analytical solution to σ (2) , as detailed in Appendix F. We calculated the ratio of S/N to (S/N)∞ as functions of α and u using Equations (A81) and (A82) in Appendix F and plotted the result in Figure 11a. When α was close to 0, the ratio had the same two peaks as in Figures 6b and 7b and then approached unity for almost all values of u as α increased. To observe the precise change in the ratio around α = 10, we plotted the ratio as a function of u by setting α at 4 to 14 in steps of 2 in Figure  11b. It was clear from the figure that at α = 8 the maximum deviation of the ratio from unity was only 23% in the 0.1 < u < 0.9 range and the deviation decreased rapidly as α increased from 8. Since a difference of around 20% was unavoidable between the calculated and measured data, as observed in Figure 8, we confirmed that we could regard the Figure 10. Comparison of the measured and calculated S/N along a 5 m long polarizationmaintaining (PM) fiber, where the up-conversion frequency was 10.881 GHz and the DFB LD was driven with current source B whose power spectral density is shown in Figure 5b. Both ends of the PM fiber were spliced to 250 µm buffered single-mode fibers with APC connectors which were mated to those of the fiber pigtails of optical circulators CL1 and CL2. z c (= 3.8 m) is the distance to the center of pumping. The fast and slow axes of the PM fiber were excited by the respective probe light and counter-propagating pump lights by using in-line polarization controllers (PCs). CL1 and CL2 are optical circulators. ω p and ω p + Ω are the angular frequencies of the counter-propagating pump light waves. APC: Angled physical contact.
We investigated the specifications of the pump light wave required for achieving the S/N at 200 even when we tested a 5 m long PM fiber. When the power spectral density was a Gaussian function, we obtained the analytical solution to σ (2) , as detailed in Appendix F. We calculated the ratio of S/N to (S/N) ∞ as functions of α and u using Equations (A81) and (A82) in Appendix F and plotted the result in Figure 11a. When α was close to 0, the ratio had the same two peaks as in Figures 6b and 7b and then approached unity for almost all values of u as α increased. To observe the precise change in the ratio around α = 10, we plotted the ratio as a function of u by setting α at 4 to 14 in steps of 2 in Figure 11b. It was clear from the figure that at α = 8 the maximum deviation of the ratio from unity was only 23% in the 0.1 < u < 0.9 range and the deviation decreased rapidly as α increased from 8. Since a difference of around 20% was unavoidable between the calculated and measured data, as observed in Figure 8, we confirmed that we could regard the ratio as unity as long as α ≥ 8, where we could use the simple expression given by Equation (20). In the relational expressions of (21) and (22), τe was the round-trip time from the origin to the output end of the fiber under test and was given by τe =2nL/c, where we assumed that the fiber was long enough to allow zi = 0. With the values L = 5 m and n = 1.45, we estimated τe to be 48.4 ns. By substituting τe = 48.4 ns and γ = 62.5 GHz (or 0.5 nm/s) into the right-hand sides of expressions (21) and (22), we obtained upper limits for fh and δνrms of 0.315 kHz and 3.3 MHz, respectively. In Figure 12 we plot red and green dotted lines, which are upper boundaries defined by fh = 0.315 kHz and δνrms = 3.3 MHz, respectively, in two-dimensional space with a (fh,δνrms) coordinate system. By substituting γ = 62.5 GHz and (S/N)tgt = 200 into the right-hand side of the relational expression (23), we found that the point (fh, δνrms) should be included in the domain located below the red solid line, which was the boundary curve defined by δνrms = 0.0293/fh and was plotted in a log-log scale, where the units of fh and δνrms are kHz and MHz, respectively. Therefore, when we sweep the laser at the rate of γ = 62.5 GHz, the point (fh, δνrms) should be included in the allowable domain surrounded by the three lines highlighted in yellow in Figure 12 to achieve (S/N)tgt = 200.  In the relational expressions of (21) and (22), τ e was the round-trip time from the origin to the output end of the fiber under test and was given by τ e =2nL/c, where we assumed that the fiber was long enough to allow z i = 0. With the values L = 5 m and n = 1.45, we estimated τ e to be 48.4 ns. By substituting τ e = 48.4 ns and γ = 62.5 GHz (or 0.5 nm/s) into the right-hand sides of expressions (21) and (22), we obtained upper limits for f h and δν rms of 0.315 kHz and 3.3 MHz, respectively. In Figure 12 we plot red and green dotted lines, which are upper boundaries defined by f h = 0.315 kHz and δν rms = 3.3 MHz, respectively, in two-dimensional space with a (f h , δν rms ) coordinate system. By substituting γ = 62.5 GHz and (S/N) tgt = 200 into the right-hand side of the relational expression (23), we found that the point (f h , δν rms ) should be included in the domain located below the red solid line, which was the boundary curve defined by δν rms = 0.0293/f h and was plotted in a log-log scale, where the units of f h and δν rms are kHz and MHz, respectively. Therefore, when we sweep the laser at the rate of γ = 62.5 GHz, the point (f h , δν rms ) should be included in the allowable domain surrounded by the three lines highlighted in yellow in Figure 12 to achieve (S/N) tgt = 200.
For example, when the spectral half width is reduced to f h = 50 Hz by installing a low pass filter in the current source circuit, the rms frequency fluctuations should be equal to or lower than δν rms = 0.59 MHz. This is determined by drawing a vertical line at f h = 50 Hz and finding the coordinate of the point A, which is the intersection between the vertical line and the red solid line. When the spectral half width is much wider due to the potential difficulty involved in installing such a low pass filter, one way to achieve (S/N) tgt is to sweep the tunable laser faster in such a way that the point (f h , δν rms ) is included in the extended domain determined by the new sweep rate. For example, by letting γ = 625 GHz or 5 nm/s, the domain highlighted in blue is added as the allowable one so that the range of the spectral half width becomes ten times wider. γ = 62.5 GHz and (S/N)tgt = 200 into the right-hand side of the relational expression (23), we found that the point (fh, δνrms) should be included in the domain located below the red solid line, which was the boundary curve defined by δνrms = 0.0293/fh and was plotted in a log-log scale, where the units of fh and δνrms are kHz and MHz, respectively. Therefore, when we sweep the laser at the rate of γ = 62.5 GHz, the point (fh, δνrms) should be included in the allowable domain surrounded by the three lines highlighted in yellow in Figure 12 to achieve (S/N)tgt = 200.

Discussion
This section describes the effect of the coherence time of the pump light on the S/N of the reflection measurement. When the counter-propagating pump light waves collided at z along the 1.35 m long optical fiber, the difference between the propagation times of the two pump light waves was given by |τ c − τ| where τ c = 2nz c /c and τ = 2nz/c, whose upper limit was τ e = 13.5 ns for all possible combinations of (τ c , τ) satisfying 0 < τ c < τ e and 0 < τ < τ e . According to the data supplied by the manufacturer of the LD, the nominal spectral linewidth was of the order of 1 MHz, and so the coherence time of the pump light waves was estimated to be a few microseconds. That is, the coherence time was much longer than any possible values of the time differences between the pump light waves, and they were considered to be coherent with each other anywhere along the 1.35 m long optical fiber. In the frequency domain, the individual spectral components in the pump light waves should receive the same phase modulation from the technical noise and generate the same acoustic wave along the entire length of the fiber.
To observe the change in the Stokes light signal level as the fiber length increased, we spliced 250 µm buffered single-mode fibers with APC connectors to both ends of a 10 m long PM fiber and installed it as the DUT in the reflectometer. After testing the 10 m long PM fiber, we prepared a 40 m long PM fiber as the DUT and tested it again. Here we increased the carrier frequency for detecting the beat signal waveform to 300 kHz by driving the phase modulators PM1 and PM3 at f 0 = 500 kHz and f 1 = 800 kHz, respectively. This was because the maximum frequency of the beat signal waveform was increased to 25 kHz. In accordance with the higher detection frequency, we increased the carrier frequency of the unbalanced MZI, which we used to measure the power spectral density, to the same frequency at 300 kHz. As the fiber length increased, unfortunately the number of data needed for the numerical integration also increased, and this resulted in a longer calculation time and a fatal memory overflow. To avoid such results, we reduced the record length of the beat signal waveform to one-tenth and increased the sampling interval of the power spectral density to 0.853 Hz.
We obtained 30 reflectograms from the 10 m long PM fiber and plotted them all and their mean reflectogram in Figure 13a,b, respectively. It was clear from the figures that all these reflectograms had peaks at the center of pumping and decayed as the position deviated from the center. In the first place, the mean Stokes light signal level <Z> should decrease due to the technical noise according to the Gaussian function as given by Equation (17). We numerically integrated H(f) with respect to f using Equation (10) and obtained δνrms = 0.60 MHz and thus Δrms = 3.74×10 6 rad/s. With the value of Δrms, τe = 115 ns and uc = 0.545, we calculated <Z> as a function of u as drawn by the red solid curve in Figure 13b and found that the decay induced by the technical noise was too small to fit the measured signal change. We employed an external-cavity tunable laser source and a DFB LD source in our experiment and both light sources were capable of causing the signal decay. If the finite coherence time of the output light from the former laser was the dominant origin of the decay, each waveform in Figure 13a,b should have a peak at z = 0 and decay as the distance increases. This was because the optical path length difference between the LO light and the Stokes light generated at z = 0 was zero and increased with distance. Actually, we observed the peak at the center of pumping, which was consistent with the fact that the optical path length difference between the counter-propagating pump light waves was zero at the center and increased as the position deviated from it. Therefore, it was clear that the decay we observed was caused by the finite coherence time of the pump light waves.
Since we measured the frequency fluctuations due to the technical noise by using an unbalanced MZI with a short delay of τMZI = 8.912 ns, we considered that the individual spectral components of the pump light passed through the MZI coherently and produced the same signal change with time, and thus the resultant power spectral density did not contain the effect of the finite coherence time. In Section 2, we denoted the amplitude of the electric field of the LD output as A(t) = A0exp(− iϕ(t)), where A0 was a constant and ϕ(t) represented the phase change due to the technical noise. To introduce the effect of the finite coherence time, we express the amplitude as A(t) = A0Ã(t)exp(− iϕ(t)) and let Ã(t) be the amplitude fluctuation caused by the finite coherence time while keeping <|Ã(t)| 2 > = 1. Then the analytic signal represented by Equation (2) should be changed to Ic(t) as follows: where t = (ῶ + ωp − ω1)/β. Considering that the fluctuation Ã(t) was statistically independent of the frequency fluctuations Δ(t) due to the technical noise, a simple way to incorporate the coherence effect into the reflection measurement is to approximate the fluctuating product Ã(t − τ1/2)Ã *(t − τc + τ1/2) in the integrand of Equation (24) by the ensemble average <Ã(t − τ1/2)Ã *(t − τc + τ1/2)>, which is denoted as V(τc − τ1) using the coherence function V(τ) defined by: We numerically integrated H(f ) with respect to f using Equation (10) and obtained δν rms = 0.60 MHz and thus ∆ rms = 3.74 × 10 6 rad/s. With the value of ∆ rms , τ e = 115 ns and u c = 0.545, we calculated <Z> as a function of u as drawn by the red solid curve in Figure 13b and found that the decay induced by the technical noise was too small to fit the measured signal change. We employed an external-cavity tunable laser source and a DFB LD source in our experiment and both light sources were capable of causing the signal decay. If the finite coherence time of the output light from the former laser was the dominant origin of the decay, each waveform in Figure 13a,b should have a peak at z = 0 and decay as the distance increases. This was because the optical path length difference between the LO light and the Stokes light generated at z = 0 was zero and increased with distance. Actually, we observed the peak at the center of pumping, which was consistent with the fact that the optical path length difference between the counter-propagating pump light waves was zero at the center and increased as the position deviated from it. Therefore, it was clear that the decay we observed was caused by the finite coherence time of the pump light waves.
Since we measured the frequency fluctuations due to the technical noise by using an unbalanced MZI with a short delay of τ MZI = 8.912 ns, we considered that the individual spectral components of the pump light passed through the MZI coherently and produced the same signal change with time, and thus the resultant power spectral density did not contain the effect of the finite coherence time. In Section 2, we denoted the amplitude of the electric field of the LD output as A(t) = A 0 exp(− iφ(t)), where A 0 was a constant and φ(t) represented the phase change due to the technical noise. To introduce the effect of the finite coherence time, we express the amplitude as A(t) = A 0 Ã(t)exp(− iφ(t)) and let Ã(t) be the amplitude fluctuation caused by the finite coherence time while keeping <|Ã(t)| 2 > = 1. Then the analytic signal represented by Equation (2) should be changed to I c (t) as follows: where t = (ω + ω p − ω 1 )/β. Considering that the fluctuation Ã(t) was statistically independent of the frequency fluctuations ∆(t) due to the technical noise, a simple way to incorporate the coherence effect into the reflection measurement is to approximate the fluctuating product Ã(t − τ 1 /2)Ã *(t − τ c + τ 1 /2) in the integrand of Equation (24) by the ensemble average <Ã(t − τ 1 /2)Ã *(t − τ c + τ 1 /2)>, which is denoted as V(τ c − τ 1 ) using the coherence function V(τ) defined by: G p (ω p +ω) was the optical power spectrum of the pump light, and we assumed that the spectrum was normalized in such a way that −∞ +∞ G p (ω p +ω)dω = 1. Then the analytic signal I c (t) was approximated by: Equation (26) means that the mean signal level <Z c >, the variance σ c (2) and the correction term σ c (4) including the coherence effect are obtained by replacing r(τ) with r(τ)V(τ c − τ) in Equations (12), (16), and (A31)-(A38) in Appendix C, respectively.
Supposing that the pump light has a Lorentzian spectrum with a full width at half maximum of δν L , we obtain a real function of V(τ) = exp(− πδν L |τ|). By introducing a new function U(u) which is defined by: <Z c > and σ c (2) for u 1 < u < (u 1 + 1)/2 are denoted as: To obtain the expression for (u i + 1)/2 < u <1, in Equation (29) we should change the interval of the first integral to (1 − u, u − u i ) and that of the second integral to (0, 1 − u), and change the sign of the variable η in the integrands. We estimated the δν L of the employed LD to be 2.0 MHz by fitting the theoretical curve defined by Equation (28) with the measured mean reflectogram, as shown by an orange curve in Figure 13b.
We calculated the S/N as a function of u from 30 reflectograms which are shown in Figure 13a and plotted the result in Figure 14a. The distribution had rapidly varying components, but the overall profile had a peak at the center of the pumping and decayed as the distance from the center increased. We calculated σ c (2) and σ c (4) as a function of u numerically by using Equation (29) and Equations (A67) to (A74) in Appendix E, respectively. Then we plotted the two distributions of S/N, which we obtained by substituting σ c (2) and σ c (2) + σ c (4) into the denominator of Equation (18), as shown by the red and blue lines, respectively, in Figure 14a. Since the calculated values differed from the measured values by a factor of up to 2, it was clear that we could obtain the approximate S/N value by using the theoretical expression even when the length of the optical fiber approached the coherence length of the pump light and the resultant S/N was degraded. However, they still had peaks on both sides of the center of pumping and there was a noticeable discrepancy between the calculated and measured reflection profiles. Then we measured the reflectograms from the 40 m long PM fiber and plotted the S/N distribution together with those which we calculated using σ c (2) and σ c (2) + σ c (4) in Figure 14b. The best S/N value was further degraded to 8, and the calculated distribution had still the same two peaks. In each figure, the S/N distribution shown in black was obtained from distributions acquired by sweeping the tunable laser 30 times. The ratio <Zc > /√σc (2) as a function of u was calculated and plotted with a red line. The correction term σc (4) is expressed by the sum of the eight terms A to H, which were represented by Equations (A67) to (A74) in Appendix E, respectively. The ratio <Zc > /√(σc (2) + σc (2) ) as a function of u was calculated and plotted with a blue line. The employed parameters were (a) Δrms = 3.74×10 6 rad/s, τe = 115 ns, uc = 0.545. (b) Δrms = 5.65×10 6 rad/s, τe = 402 ns, uc = 0.515.
Since the correction term σc (4) did not lead to a fundamental solution for suppressing the side peaks, we concluded that the discrepancy was not caused by the approximation of the factor exp{− iΔ(t)(τc−τ)} in the integrand of <Ic(t)> as the power series expansion. In deriving the expressions for σc (2) and σc (2) + σc (4) , we made one more approximation such that where we theoretically showed that the approximation error was of the order of π 2 fuδνrmsτe 2 /2. When we tested the 40 m long PM fiber, we had fu = 25 kHz, δνrms = Δrms/2π = 0.9 MHz and τe = 402 ns, and with these values we estimated the approximation error to be around 0.018 rad, and this meant that the approximation was still valid. Therefore, the origin of the noticeable discrepancy that we observed resulted not from the error due to the two kinds of approximations that we employed but from the simple method we introduced as an effect of the finite coherence time.
Returning to the basics of statistical averaging, therefore, we should calculate the mean value and variance of the absolute square of the Fourier inverse transform with two different kinds of statistically independent random processes such as Gaussian noise due to technical noise and AM and phase noise inherent in the laser diode [25,26]. If our reflectometry technique is to be applied to mid and long-range distributed strain sensing, the employed optical fiber will range in length from hundreds of meters to several kilometers. As described in our experiment, S/N will be degraded both by frequency fluctuations due to technical noise and by the finite coherence time unless we use an excellent DFB LD as the pump light source. Once we succeed in deriving the theoretical formula for S/N including the finite coherence time, we will be able to obtain the detailed specifications for the pump light source needed to achieve a high S/N even when we test an optical fiber several kilometers in length.

Conclusions
Signal-dependent speckle-like noise has been a serious factor in Brillouin-grating based coherent FMCW reflectometry. In this paper we theoretically and experimentally showed that the noise was generated by the frequency fluctuations of the pump light from the DFB LD. By assuming that the frequency of the pump light was modulated by the technical noise from the employed current source, we derived theoretical formulas for the mean value and the variance of the fluctuating power of the Stokes light, which contained the second and fourth-order moments of the frequency fluctuations. We adopted six experimental conditions for the reflectogram measurement where we used two current sources with different technical noises and set the center of pumping to the center and both sides of a 1.35 m long optical fiber under test. Under each condition we numerically Since the correction term σ c (4) did not lead to a fundamental solution for suppressing the side peaks, we concluded that the discrepancy was not caused by the approximation of the factor exp{− i∆(t)(τ c −τ)} in the integrand of <I c (t)> as the power series expansion. In deriving the expressions for σ c (2) and σ c (2) + σ c (4) , we made one more approximation such that φ(t − τ/2) − φ(t − τ c + τ/2) ≈ ∆(t)(τ c − τ), where we theoretically showed that the approximation error was of the order of π 2 f u δν rms τ e 2 /2. When we tested the 40 m long PM fiber, we had f u = 25 kHz, δν rms = ∆ rms /2π = 0.9 MHz and τ e = 402 ns, and with these values we estimated the approximation error to be around 0.018 rad, and this meant that the approximation was still valid. Therefore, the origin of the noticeable discrepancy that we observed resulted not from the error due to the two kinds of approximations that we employed but from the simple method we introduced as an effect of the finite coherence time.
Returning to the basics of statistical averaging, therefore, we should calculate the mean value and variance of the absolute square of the Fourier inverse transform with two different kinds of statistically independent random processes such as Gaussian noise due to technical noise and AM and phase noise inherent in the laser diode [25,26]. If our reflectometry technique is to be applied to mid and long-range distributed strain sensing, the employed optical fiber will range in length from hundreds of meters to several kilometers. As described in our experiment, S/N will be degraded both by frequency fluctuations due to technical noise and by the finite coherence time unless we use an excellent DFB LD as the pump light source. Once we succeed in deriving the theoretical formula for S/N including the finite coherence time, we will be able to obtain the detailed specifications for the pump light source needed to achieve a high S/N even when we test an optical fiber several kilometers in length.

Conclusions
Signal-dependent speckle-like noise has been a serious factor in Brillouin-grating based coherent FMCW reflectometry. In this paper we theoretically and experimentally showed that the noise was generated by the frequency fluctuations of the pump light from the DFB LD. By assuming that the frequency of the pump light was modulated by the technical noise from the employed current source, we derived theoretical formulas for the mean value and the variance of the fluctuating power of the Stokes light, which contained the second and fourth-order moments of the frequency fluctuations. We adopted six experimental conditions for the reflectogram measurement where we used two current sources with different technical noises and set the center of pumping to the center and Under each condition we numerically calculated the mean value and variance along the optical fiber using the data for the power spectral density of the frequency fluctuations, and the resulting S/N distribution agreed with the measured distribution within 20% in most parts of the fiber. Although the reflectometry was proposed originally for the short-range diagnosis of miniaturized optical waveguides, our success shows that the reflectometry has a great potential for extending its application area to the middle and long-range fiber-optic distributed strain sensing. τ e )/2 < τ < τ e are obtained by changing between ξ 1 and ξ 2 and by changing the sign of the variable χ in the respective expressions for τ i <τ < (τ i + τ e )/2. Since the integrands of A and B, which are represented by Equations (A31) and (A32), respectively, have the same forms as the first term in Equation (16), we can straightforwardly express the results for A and B as: Considering that C and D, which are expressed by Equation (A33) and (A34), respectively, contain the same function r(τ + χ/β) r (τ − χ/β) in their integrands as the second term in Equation (16), we find that: should be taken in the domain which is defined by -ξ 1 ≤ χ 1 , χ 2 , χ 1 + χ 2 ≤ξ 2 .
In this appendix we derive the expressions for the terms A to H, which all contain the same factor |r 0 | 4 . The correction term σ (4) is the sum of these terms so that the variance of σ (2) + σ (4) is also proportional to |r 0 | 4 , meaning that the standard deviation is proportional to |r 0 | 2 . Considering that the mean signal level is also proportional to|r 0 | 2 as clearly seen from Equation (12), S/N, which is defined by the ratio of the mean signal level to the standard deviation, either no longer includes the factor r 0 , or it is independent of r 0 . From this independence, we redefine the individual terms A to H as the ones that are obtained by letting |r 0 | = 1 in their original expressions.
The correction term σ c (4) including the coherence effect are obtained by replacing r(τ) with r(τ)V(τ c − τ) in the constituent terms represented by Equations (A31) to (A38). Assuming that the coherence function V(τ) defined by Equation (25) is a real-valued function, the constituent terms A to H represented by Equations (A59) to (A66) should be changed to the following expressions for u i < u < (u i + 1)/2 using the function U(u) which is defined by U(u) = V(τ e u):