2. System configuration and operation principle
Figure 1
shows the proposed test method, which uses a DFB-LD, a photodiode (PD), and an electrical spectrum analyzer (ESA). The DFB-LD periodically emits a chirped test signal with lightwave frequencies of f
1 and f
2 into a distributed amplifier or transmission fiber of several tens to hundreds of kilometers under test. The double-Rayleigh-scattered powers of the pulsed test signal of f
1 and f
2 overlap the direct signal of f
1 and f
2 at the PD. The direct signal acts as the pulsed local light in a heterodyne receiver. The ESA detects the electrical power of the beat current between the double-Rayleigh-scattered power and the direct signal at the frequency of |f
1-f
2|.
Fig. 1 Proposed test method for MPI measurement using chirped test signal with self-heterodyne detection
The detected electrical power of the beat current corresponds to the optical power product of the total forward-scattered signal (
) and direct signal (
) as
. Multipath interference is defined as the optical power ratio of the total forward-scattered signal to the direct signal:
. Therefore, MPI is simply obtained using these relations as
where
is a coefficient to be determined with the calibration procedure. Calibration was carried out using Rayleigh backscattered light of the test signal from a 1.55-μm dispersion- shifted single-mode fiber (SMF/DS) of 25.3 km in length.
The advantage of our method is that the beat spectrum ideally concentrates at two frequencies: zero and |f1-f2|. The frequency |f1-f2| can be appropriately determined to be in the frequency band of low noise. Consequently, highly sensitive MPI measurement becomes possible.
3. Synthesis of chirped test signal
To obtain the chirped test signal, basic characteristics of the DFB-LD need to be clarified. Chirped characteristics are measured using the self-heterodyne fiber loop interferometer shown in
Fig. 2
. A pulsed reference frequency is sampled by SW1 with a 0.2-μs window one time every 5.12 ms (256 times x 20 μs pulse period) from the 50 kHz repetition rate chirped stream and stored in a loss compensated fiber loop with an EDFA. The repetition rate of 50 kHz is selected by the pulse period of nearly twice time of the fiber loop round trip time of 10.1 μs. SW3 selects a pulsed reference frequency of N circuit from the fiber loop outputs. SW2 is used to initialize the loss compensated fiber loop. The beat frequency between the chirped pulse stream and the reference pulse via the fiber loop after circulation is measured using an ESA. The lightwave frequency of the reference pulse into the PDs is shifted by −240 MHz due to acousto-optical modulators used as SW1 and SW3. Because of the time difference between the delay and the half period of the pulse, the temporal beating position (sample timing) in the pulse period shifts according to the number of circuits. By reconstructing the ESA output, the temporal chirp characteristics over a pulse period can be evaluated. This is similar to a sampling oscilloscope operation. Owing to the self-heterodyne configuration and the ESA with a resolution bandwidth of 300 kHz, highly stable chirp evaluation of less than 1 MHz becomes possible.
Fig. 2 Chirp characteristics measurement using self-heterodyne fiber loop interferometer
Figure 3
shows the chirp characteristics of the output emitted from a directly modulated 1558- nm DFB-LD at a repetition rate of 50 kHz by a rectangular pulse stream. The threshold current of the DFB-LD was 13.4 mA and the optical power at an operating current of 40 mA was 6.5 mW. The temporal chirp response was clearly obtained from the beat frequency in the
Fig. 3 (b).
Fig. 3 Chirp characteristics of output lightwave emitted from directly modulated DFB-LD with 50-kHz rectangular pulse stream
Chirped characteristics arise from the temperature and carrier density change of the LD active layer caused by the electrical modulation signal [
7H. Shalom, A. Zadok, M. Tur, P. J. Legg, W. D. Cornwell, and I. Andonovic, “On the various time constants of wavelength changes of a DFB laser under direct modulation,” IEEE J. Quantum Electron. 34(10), 1816–1822 (1998). [CrossRef]
-
8A. Yariv, Optical Electronics in Modern Communications fifth edition (Oxford University Press, New York, USA, 1997), Chap. 15.
]. This process is modeled using the following equation:
where
Δf(t) is the frequency shift from the average (
):
,
ΔJQ is the heat generated from the modulation current
ΔILD,
z(f) is the frequency domain thermal impedance,
F−1 is the inverse Fourier transform, and
kf is a coefficient.
The parameters in
Eq. (2) are extracted by measuring the temporal response as a function of the bias current and pulse amplitude. From the peak-to-peak chirp response, an empirical formula for
ΔJQ and
ΔILD relation is expressed as
where
a and
b are coefficients, and
Ibias is the bias current. Thermal impedance of the analytical expression obtained from a two-dimensional diffusion equation of a cylindrical coordinate is used:
where
ω is the angular frequency,
j is
,
K0 and
K1 are modified Bessel functions. The time constant
τd is extracted by fitting the slope of the temporal chirp response. The coefficient
kz is obtained from the peak-to-peak chirp response and the time constant. The coefficient
kf in the second term of
Eq. (2) is extracted from the frequency jump observed at the alternate instant of the modulation rectangular pulse polarity. The experimentally obtained parameters in
Eqs. (2)–
(4) are summarized as
a, 353 (MHz/mA);
b, 3 (MHz/mA);
kf, 26 (MHz/mA);
kz, –0.58; and
τd, 5(μs).
The first and second terms of
Eq. (2) are the dominant and fractional roles in determining the chirp characteristics, respectively. To periodically synthesize the appropriate chirped test signal with lightwave frequencies of f
1 and f
2, the DFB-LD modulation current
ΔILD is calculated from
Eq. (2) according to the calculation algorithm for DSP shown in
Fig. 4
. Analytical calculation is difficult due to chirp model nonlinearity. An iterative process of about five steps is used. In blocks 2 and 4, calculation with the chirp model is carried out. First, inverse calculation from tentative target pattern to drive current is carried out with the dominant term only at block 2. Second, the effect of the fractional term is taken into account at block 4. Blocks 1, 3, and 5 are pre and post calculation blocks for the iterative process. In block 3, low pass filters and a limiter block, the operating range of the drive current is restricted appropriately. By monitoring the convergence, the driving signal is obtained from block 3. For convergence of the iterative process, an appropriate cutoff frequency is important. In our case, about 5MHz was used.
Fig. 4 Outline of drive signal calculation algorithm
The tentative target chirp patterns for the first n = 1 and successive nth iterations are expressed as
The inverse calculation in block 2 is carried out using fast Fourier transform (FFT) and inverse FFT (IFFT) as follows:
The conversion of
ΔJQ to
ΔILD is carried out by solving
Eq. (3).
Figure 5(a)
shows a temporal target chirp pattern for synthesis. The peak-to-peak frequency deviation |f
1-f
2| was 300 MHz, the transition time from f
1 to f
2 was 0.9 μs, and the repetition rate was 50 kHz. We selected the repetition rate of 50 kHz for the following two reasons. First, the systematic error caused by the low repetition rate test pulse in the MPI measurement is expected to be less than 0.5 dB for a test fiber length that exceeds 2km by the 50 kHz repetition rate test signal [
6K. Aida, T. Okada, and Y. Hinako, “Multipath interference test method for distributed amplifier using self-heterodyne technique,” IEICE Trans. ELECTRON , E 90-C(1), 18–24 (2007).
]. Second, the self-heterodyne fiber loop interferometer for measuring chirp characteristics can operate around 50 kHz.
Figure 5(b) shows a synthesized drive signal at a 40-mA bias current attained after iterations of five loops. The frequency variation from a temporal target chirp pattern is expected to be less than 0.2 MHz, except the transition time from the simulation. At the transition time, large frequencies over shoots by 10 MHz are anticipated.
Figure 5(c) shows the measured optical power waveform emitted from the DFB-LD modulated by the synthesized drive signal, and
Fig. 5(d) shows the measured temporal chirp pattern obtained by the synthesized drive signal. The peak-to-peak frequency deviation of 300 MHz was attained as expected. The frequencies, f
1 and f
2, exhibited variations of about +/− 1 MHz, which were larger than those in the simulation(0.2MHz), except the transition time of 1 μs. This discrepancy may be due to parameters errors used in the simulation. To overcome the frequency variation of about +/−1 MHz, the resolution bandwidth (RBW) of the ESA is adjusted to the widest value of 3 MHz for detecting the beat power with the frequency variation around |f
1-f
2|, although the minimum ESA RBW is 300 kHz. Then the frequency variations are less than the ESA RBW.
Figure 5(e) shows that the beat spectrum can be easily detected with the ESA for the test signal passing through a 25.3-km SMF/DS. A large spectrum peak was observed at 300 MHz. Thus, the drive signal was successfully synthesized for enabling highly sensitive MPI measurement.
Fig. 5 Synthesized voltage waveform for DFB-LD direct modulation and attained test signal characteristics
4. MPI measurements of single-mode fibers
The transmitter and receiver pair calibration for the MPI measurement was carried out using a Mach-Zehnder interferometer of two paths, which was described by Aida et al. [
6K. Aida, T. Okada, and Y. Hinako, “Multipath interference test method for distributed amplifier using self-heterodyne technique,” IEICE Trans. ELECTRON , E 90-C(1), 18–24 (2007).
]. One path contained an attenuator for adjusting the direct signal power to the receiver. The other path contained a Rayleigh backscattered light generator consisted of a three-port circulator and a far-end-terminated 25.3-km SMF/DS, and provided sufficient optical power to compensate the insertion loss of the optical components. The light from each path can be measured using an optical power meter and can be adjusted. Using this procedure, a calibrated MPI test signal is generated in the desired value. Precise calibrations can be carried out in the receiver saturated power region as well as liner operation. In
Fig. 6
, the MPI values are calibrated at several reference points in the range of −74 to −40 dB.
Fig. 6 MPI measurement of SMF and SMF/DS
Multipath interference measurements were carried out using various lengths of a standard SMF and SMF/DS. The coupled power into the photodiode of the receiver was adjusted between −7 and −3 dBm, appropriately.
Figure 6 shows the relationship between MPI and fiber length, plotted for the SMF (●) and SMF/DS (▲). The solid lines in
Fig. 6 are theoretical MPIs for both fibers, calculated from the following equation:
The parameters in
Eq. (7) for the fiber under test were obtained from the backscattering and loss measurements as follows:
The measured MPI values agreed well with the theoretical ones for both fibers of over 10 km in length. The discrepancies between the measured MPI values and the theoretical ones over 10 km length were less than +/− 0.3 dB. A relatively large discrepancy of 1.4 dB was observed for the SMF of 3 km length. The discrete reflections from connectors at the fiber both ends were arisen, although SPC (about −50 dB) connectors were used. Because the inherent MPI of the 3-km fiber was calculated to be as small as −79.2 dB, discrete reflections strongly affected MPI measurement. This test method improved the minimum detectable MPI to as low as −78 dB, compared with that of −50 dB with conventional test methods [
1C. R. S. Fludger and R. J. Mears, “Electrical measurements of multipath interference in distributed Raman amplifiers,” J. Lightwave Technol. 19(4), 536–545 (2001). [CrossRef]
,
3S. Faralli and F. Di Pasquale, “Impact of double Rayleigh scattering noise in distributed higher order Raman pumping schemes,” IEEE Photon. Technol. Lett. 15(6), 804–806 (2003). [CrossRef]
] and −70 dB with a method that uses an analog equalizing circuit for emitting a chirped test signal from a DFB-LD [
6K. Aida, T. Okada, and Y. Hinako, “Multipath interference test method for distributed amplifier using self-heterodyne technique,” IEICE Trans. ELECTRON , E 90-C(1), 18–24 (2007).
].
The proposed method can estimate the MPI value in the rage of −78 to −57 dB as shown in
Fig. 6, because of the synthesized chirped test signal with the self-heterodyne detection and their calibrations. The estimated MPI values agreed well with the theoretical ones calculated by
Eq. (7) for both fibers of over 10 km in length. This means that a fiber itself can be used as a reference MPI generator for low range of MPI value. But the proposed method is more effective to measure the wide range of MPI value.
Figure 7
shows the relationship between the MPI of fiber lines containing discrete reflections at both ends plotted for 3.0- and 40.1-km SMFs. The reflection coefficient at the far end of the fiber varied from −47 to −26.5 dB. The near end reflection coefficient was fixed at −46 dB. The solid lines in
Fig. 7 are theoretical MPIs calculated from the following equation:
Fig. 7 MPI measurement of fiber lines with discrete reflections
The measured MPI values agreed well with the theoretical ones for both fiber lines with discrete reflections.