Distributed Optical Fiber Vibration Measurement Device and Distributed Optical Fiber Vibration Measurement Method

Abstract

A distributed optical fiber vibration measurement device includes a phase constant difference computation unit causing the first backscattered light generated at the points of a plurality of optical fibers under test that are integrated to interfere with another light to obtain two AC components and determining a phase constant difference from the two AC components, a phase distribution data creation unit comparing amplitudes between the two AC components obtained by causing the first backscattered light generated at the points of the optical fibers under test to interfere with the other light and selecting, for each of the points, phase data regarding one of the two AC components having a larger amplitude to create phase distribution data with the phase data, and a vibration measurement unit identifying an optical path length difference between any two points of the optical fiber under test to measure vibration of the optical fiber under test.

Claims

1. A distributed optical fiber vibration measurement device for causing a first test light pulse to enter a plurality of optical fibers under test that are integrated to generate first backscattered light at points of the plurality of optical fibers under test and heterodyne-detecting the first backscattered light to measure vibration applied to positions of an optical fiber under test of the plurality of optical fibers under test, the distributed optical fiber vibration measurement device comprising: a phase constant difference computation unit configured to cause the first backscattered light generated at the points of the plurality of optical fibers under test to interfere with another light having an optical frequency different from an optical frequency of the first test light pulse to obtain two AC components and determine a phase constant difference from the two AC components; a phase distribution data creation unit configured to compare amplitudes between the two AC components obtained by causing the first backscattered light generated at the points of the plurality of optical fibers under test to interfere with the other light and select, for each of the points, phase data regarding an AC component of the two AC components having a larger amplitude to create phase distribution data with the phase data arranged over time; and a vibration measurement unit configured to identify, from the phase constant difference and the phase distribution data, an optical path length difference between any two points of the optical fiber under test to measure vibration of the optical fiber under test.

2. The distributed optical fiber vibration measurement device according to claim 1, wherein the other light is second backscattered light generated by causing a second test light pulse having an optical frequency different from the optical frequency of the first test light pulse to enter the plurality of optical fibers under test at a predetermined interval from the first test light pulse.

3. A distributed optical fiber vibration measurement device for causing a first test light pulse and a second test light pulse having an optical frequency different from an optical frequency of the first test light pulse to enter an optical fiber under test to generate first backscattered light and second backscattered light at points of the optical fiber under test and heterodyne-detecting the first backscattered light and the second backscattered light to measure vibration applied to positions of the optical fiber under test, the distributed optical fiber vibration measurement device comprising: a phase constant difference computation unit configured to cause the first backscattered light and the second backscattered light to interfere respectively with another light having an optical frequency different from optical frequencies of the first test light pulse and the second test light pulse to obtain two AC components and determine a phase constant difference from the two AC components; a phase distribution data creation unit configured to compare amplitudes between the two AC components obtained by causing the first backscattered light and the second backscattered light to interfere with the other light respectively and select, for each of the points, phase data regarding an AC component of the two AC components having a larger amplitude to create phase distribution data with the phase data arranged over time; and a vibration measurement unit configured to identify, from the phase constant difference and the phase distribution data, an optical path length difference between any two points of the optical fiber under test to measure vibration of the optical fiber under test.

4. The distributed optical fiber vibration measurement device according to claim 3, wherein the other light includes third backscattered light and fourth backscattered light generated by causing a third test light pulse and a fourth test light pulse having optical frequencies different from the optical frequencies of the first test light pulse and the second test light pulse to enter the optical fiber under test at predetermined intervals from the first test light pulse and the second test light pulse.

5. A distributed optical fiber vibration measurement method of causing a first test light pulse to enter a plurality of optical fibers under test that are integrated to generate first backscattered light at points of the plurality of optical fibers under test and heterodyne-detecting the first backscattered light to measure vibration applied to positions of an optical fiber under test of the plurality of optical fibers under test, the distributed optical fiber vibration measurement method comprising: causing the first backscattered light generated at the points of the plurality of optical fibers under test to interfere with another light having an optical frequency different from an optical frequency of the first test light pulse to obtain two AC components and determining a phase constant difference from the two AC components; comparing amplitudes between the two AC components obtained by causing the first backscattered light generated at the points of the plurality of optical fibers under test to interfere with the other light and selecting, for each of the points, phase data regarding an AC component of the two AC components having a larger amplitude to create phase distribution data with the phase data arranged over time; and identifying, from the phase constant difference and the phase distribution data, an optical path length difference between any two points of the optical fiber under test to measure vibration of the optical fiber under test.

6. A distributed optical fiber vibration measurement method of causing a first test light pulse and a second test light pulse having an optical frequency different from an optical frequency of the first test light pulse to enter an optical fiber under test to generate first backscattered light and second backscattered light at points of the optical fiber under test and heterodyne-detecting the first backscattered light and the second backscattered light to measure vibration applied to positions of the optical fiber under test, the distributed optical fiber vibration measurement method comprising: causing the first backscattered light and the second backscattered light to interfere respectively with another light having an optical frequency different from optical frequencies of the first test light pulse and the second test light pulse to obtain two AC components and determining a phase constant difference from the two AC components; comparing amplitudes between the two AC components obtained by causing the first backscattered light and the second backscattered light to interfere with the other light respectively and selecting, for each of the points, phase data regarding an AC component of the two AC components having a larger amplitude to create phase distribution data with the phase data arranged over time; and identifying, from the phase constant difference and the phase distribution data, an optical path length difference between any two points of the optical fiber under test to measure vibration of the optical fiber under test.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0022] FIG. 1 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a first embodiment of the present disclosure.

[0023] FIG. 2 is a diagram illustrating a test light pulse delivery sequence according to the first embodiment of the present disclosure.

[0024] FIG. 3 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a first embodiment of the present disclosure.

[0025] FIG. 4 is a diagram illustrating a test light pulse delivery sequence according to the first embodiment of the present disclosure.

[0026] FIG. 5 is a diagram illustrating a calculation processing flow in a calculation processing device according to the first embodiment of the present disclosure.

[0027] FIG. 6 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a second embodiment of the present disclosure.

[0028] FIG. 7 is a diagram illustrating a test light pulse delivery sequence according to the second embodiment of the present disclosure.

[0029] FIG. 8 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a third embodiment of the present disclosure.

[0030] FIG. 9 is a diagram illustrating a test light pulse delivery sequence according to the third embodiment of the present disclosure.

[0031] FIG. 10 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a fourth embodiment of the present disclosure.

[0032] FIG. 11 is a diagram illustrating a test light pulse delivery sequence according to the fourth embodiment of the present disclosure.

DESCRIPTION OF EMBODIMENTS

[0033] Embodiments of the present disclosure will be described in detail below with reference to the drawings.

[0034] A distributed optical fiber vibration measurement device described in the embodiments causes a first test light pulse obtained by frequency-shifting light from a light source to enter each of an integrated plurality of test optical fibers. The distributed optical fiber vibration measurement device causes first backscattered light generated at a plurality of points of the plurality of optical fibers under test to interfere with another light for heterodyne detection to measure vibration applied to each of the positions of the optical fiber under test. The distributed optical fiber vibration measurement device includes a phase constant difference computation unit, a phase distribution data creation unit, and vibration measurement unit. The phase constant difference computation unit causes the first backscattered light generated at the plurality of points of the plurality of optical fibers under test to interfere with the other light, and determines a phase constant difference from two AC components obtained. The phase distribution data creation unit causes the first backscattered light generated at the plurality of points of the plurality of optical fibers under test to interfere with the other light and compares the amplitudes between two AC components obtained. Then, the phase distribution data creation unit selects phase data regarding the AC component with the larger amplitude for each of the plurality of points and arranges the phase data over time to create phase distribution data. The vibration measurement unit identifies, from the phase constant difference and the phase distribution data, an optical path length difference between any two points of each of the optical fibers under test to measure vibration of the optical fiber under test.

[0035] In addition, the distributed optical fiber vibration measurement device described in the embodiments causes a first test light pulse and a second test light pulse to enter optical fibers under test, the first test light pulse being obtained by subjecting light from a light source to a first frequency shift and the second test light pulse being obtained by subjecting light from the light source to a second frequency shift. The distributed optical fiber vibration measurement device causes first backscattered light and second backscattered light generated at a plurality of points of the optical fibers under test to interfere with another light for heterodyne detection, and measures vibration applied to a plurality of positions of the optical fibers under test. The distributed optical fiber vibration measurement device includes a phase constant difference computation unit, a phase distribution data creation unit, and a vibration measurement unit. The phase constant difference computation unit causes each of the first backscattered light and the second backscattered light to interfere with another light to obtain two AC components, and determines a phase constant difference from the two AC components. The distribution data creation unit compares the amplitudes between two AC components obtained by causing each of the first backscattered light and the second backscattered light to interfere with the other light. The distribution data creation unit selects phase data regarding the AC component with the larger amplitude for each of the plurality of points and arranges the phase data over time to create phase distribution data. The vibration measurement unit identifies, from the phase constant difference and the phase distribution data, an optical path length difference between any two points of each of the optical fibers under test to measure vibration of the optical fiber under test.

[0036] According to the distributed optical fiber vibration measurement device, by selecting and employing one of the two types of signals having a higher signal intensity as illustrated in FIG. 2, vibration detection can be performed by using a signal with a high signal-to-noise ratio, thus reducing the influence of fading noise. In a method for utilizing local light and backscattered light of a test light pulse resulting from a frequency shift to heterodyne-detect interference between the local light and the backscattered light, an accurate and sensitive distributed optical fiber vibration measurement can be achieved with fading noise reduced. Alternatively, in a method for heterodyne-detecting interference between rays of backscattered light by using a plurality of test light pulses frequency-shifted from each other, accurate and sensitive distributed optical fiber vibration measurement can be achieved with fading noise reduced. Furthermore, the distributed optical fiber vibration measurement device can contribute to reducing noise caused by optical frequency drift, which is involved in fading noise and which causes noise when low-frequency vibration is measured.

First Embodiment

[0037] FIG. 3 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device according to a first embodiment. The distributed optical fiber vibration measurement device of the present embodiment branches light from a single light source into rays of light used as test light and local light, multiplexes backscattered light with the local light for heterodyne detection, and acquires, from the optical detection unit, AC components of two electrical signals. Here, the backscattered light is obtained by causing test light to simultaneously enter two different optical fibers under test.

[0038] Output light from a light source 1 emitting coherent light is branched into two lines by a branching element 2. One ray of the branched light is used as local light, whereas the other is incident on an optical frequency control unit 3 as test light. The optical frequency control unit 3 applies a frequency shift Δω.sub.A between the test light and the local light.

[0039] The test light from the optical frequency control unit 3 is converted into a single pulse by the light pulsing unit 5, and the pulse is caused to enter two optical fibers under test 10 via an optical amplifier 6 and a branching element 8. The light pulsing unit 5 is driven by a drive signal from a pulse signal generation unit 7. Backscattered light generated in the optical fibers under test 10 is multiplexed with the local light fed from the branching element 11 via an optical circulator 9, and the resultant light is converted into electrical signals by balanced optical reception units 12 and 13. The balanced optical reception units 12 and 13 can multiplex the backscattered light with the local light for heterodyne detection to acquire AC components of the two electrical signals.

[0040] Outputs from the balanced optical reception units 12 and 13 are partially branched, the partial outputs are mixed by a mixer 18, and from the resultant signal, high frequency components are removed via a filter 17. Then, the signal is converted to a digital signal by a digitization unit 14, and the digital signal is processed by a calculation processing unit 15. Additionally, output signals from the balanced optical reception units 12 and 13 are input to the digitization unit 14 similarly to the input to the mixer 18, and converted into digital signals by the digitization unit 14, and the digital signals are processed by the calculation processing unit 15.

[0041] The results of the signal processing by the calculation processing unit 15 are displayed on the display unit 16 as appropriate.

[0042] Here, the two optical fibers under test need to be integrated in such a manner as to be vibrated in the same manner. The integration refers to a state in which the two optical fibers under test are physically bound to be vibrated in the same manner. Specifically, the optical fibers under test may be bundled in the same coating resin as in an optical fiber tape, or may be tightly formed into a cable, or may have a multiple waveguide structure as in a multi-core optical fiber, for example.

[0043] Here, it is sufficient that the optical frequency control unit 3 is specifically an external modulator having a function for changing the optical frequency in accordance with an electrical signal from a driving unit 4 used as a drive source, and as is well known, the function is provided in a phase modulator or an amplitude modulator using an acoustic optical switch or LiNbO.sub.3, or an Single Side Band Suppressed Carrier (SSB-SC) modulator.

[0044] It is sufficient that the light pulsing unit 5 has a function for converting test light into a pulse, and as is well known, the function is provided in a high-speed variable attenuator, a semiconductor optical amplifier, an acoustic optical switch, and the like. Note that, for the light pulsing unit 5, an acoustic optical switch can be used to provide the function of the optical frequency control unit 3, or a semiconductor optical amplifier can be used to provide the function of the optical amplifier 6.

[0045] FIG. 4 illustrates a test light pulse delivery sequence according to the present example. The optical frequency of the light from the light source 1 is ω. A test light pulse with an optical frequency ω+Δω.sub.A resulting from a frequency shift by the optical frequency control unit 3 is repeatedly pulsed with a pulse width W by the light pulsing unit 5, the pulse is amplified by the optical amplifier 6, and the amplified pulse is branched by the branching element 8 into two pulses, which are then delivered to the respective optical fibers under test 10.

[0046] The repetition time of the test light pulse is set equal to or longer than the round-trip propagation delay time of the maximum length of the optical fiber under test 10.

[0047] Here, AC components I.sub.1(t) and I.sub.2(t) of beat electrical signals output from the two balanced optical reception unit 12 and 13 are respectively given by (Expression 3) and (Expression 4) below, with respect to the output time t and the optical frequency ω of the coherent light from the light source, the frequency shift Δω.sub.A, and the pulse width W of the test light pulse.

[00003]$\begin{array}{cc}\mathrm{Math}.3& \\ I_1\left(t\right)\propto 2E_L{E}_S\underset{i=1}{\overset{N}{.Math.}}{a}_i\cos \left\{\Delta {\omega }_{A}t-\Delta {\omega }_A{\tau }_i+\theta \left(t-{\tau }_i\right)-\theta \left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right\}.Math.\mathrm{rect}\left(\frac{t-{\tau }_i}{W}\right)=2E_L{E}_S{R}_t\left(t\right)\cos \left({\mathrm{\Delta \omega }}_{A}t+{\theta }_1+{\mathrm{\Delta \theta }}_L\left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right)& \left(\mathrm{Expression}3\right)\\ \mathrm{Math}.4& \\ I_2\left(t\right)\propto 2E_L{E}_S\underset{i=1}{\overset{N}{.Math.}}{b}_i\cos \left\{\Delta {\omega }_{A}t-{\mathrm{\Delta \omega }}_A{\tau }_i+\theta \left(t-{\tau }_i\right)-\theta \left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right\}.Math.\mathrm{rect}\left(\frac{t-{\tau }_i}{W}\right)=2E_L{E}_S{R}_2\left(t\right)\cos \left({\mathrm{\Lambda \omega }}_{A}t+{\theta }_2+{\mathrm{\Lambda \theta }}_L\left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right)& \left(\mathrm{Expression}4\right)\end{array}$

[0048] Here, E.sub.L and E.sub.S are respectively electric field amplitudes of the local light and the test light, N is the number of scatterers present in the length direction of each of the optical fibers under test, i (1≤i≤N) is the number of each scatterer, a.sub.i and b.sub.i are each the reflectance of the scatterer, τ.sub.i is the round-trip propagation delay time from an incident end of the optical fiber under test to the scatterer i, θ(t) is an initial phase of the test light, and θ.sub.FUT(τ) is a phase change caused by vibration applied to the optical fiber under test and is an object to be measured. R.sub.1(t) and R.sub.2(t) are backscattering coefficient distributions of the respective optical fibers under test, and Δθ.sub.L(t) means an initial phase difference of light between the time of output of the test light pulse and the time of interference with the local light, i.e., phase noise from the laser, and is the same between the two optical fibers under test. According to assumed conditions, θ.sub.FUT(τ) is similarly the same between the two optical fibers under test. θ.sub.1 and θ.sub.2 are each a phase constant attributed to fading, and conditions for interference caused by fading, in other words, scatterers forming the optical fiber, vary among the optical fibers under test, and thus the phase constants differ from each other.

[0049] The AC components of (Expression 3) and (Expression 4) are converted into digital signals by the digitization unit 14, and the digital signals are input to the calculation processing unit 15. Note that in the present embodiment, the frequency of the AC component is Δω.sub.A, and the two balanced optical reception units 12 and 13 and the digitization unit 14 require a band of 2Δω.sub.A or higher on the basis of the Nyquist theorem.

[0050] Before processing by the digitization unit 14 at the subsequent stage of the two balanced optical reception unit 12 and 13, the electrical signals resulting from branching are synthesized by the mixer 18. A synthetic AC signal which is the synthesized AC signal, is given by the product of the AC components, and is as represented by (Equation 5).

Math. 5

I.sub.1(t)×I.sub.2(t)=4E.sub.L.sup.2E.sub.S.sup.2R.sub.1(t)R.sub.2(t)[cos{2(Δω.sub.xt+Δθ.sub.L(t)+θ.sub.FUT(t))+θ.sub.1+θ.sub.2}+cos(ϕ.sub.1,2)]  (Equation 5)

[0051] Here, ϕ.sub.1,2=θ.sub.1−θ.sub.2. As seen in (Equation 5), the spectrum of the synthetic AC signal includes a signal component centered at a double frequency of the frequency shift Δω.sub.A applied by the optical frequency control unit 3 and a signal component with an amplitude modulated into the vicinity of a direct current according to the phase constant difference |θ.sub.1−θ.sub.2|. The DC component can be easily separated from the signal component centered at the double frequency of Δω.sub.A by using a low-pass filter (filter 17).

[0052] Next, assuming that i (i=1 or 2) is the number of each AC component, an amplitude A.sub.i(t) and a phase P.sub.i(t) for each of the two AC components represented in (Expression 3) and (Expression 4) are given from the AC components by (Equation 6) and (Equation 7) below.

[00004]$\begin{array}{cc}\mathrm{Math}.6& \\ A_i\left(t\right)=\sqrt{{I_i\left(t\right)}^2+H{\left\{I_i\left(t\right)\right\}}^2}& \left(\mathrm{Equation}6\right)\\ \mathrm{Math}.7& \\ P_i\left(t\right)=\arctan \left[H\left\{I_i\left(t\right)\right\}/I_i\left(t\right)\right]-{\mathrm{\Lambda \omega }}_{A}t& \left(\mathrm{Equation}7\right)\end{array}$

[0053] Here, H{I(t)} means Hilbert transform of a signal I(t) and is an orthogonal signal obtained by shifting the phase of I(t) by 90 degrees. Similarly, by performing Hilbert transform on the DC component of the synthetic AC signal, the phase constant difference ϕ.sub.1,2 can be recognized from the DC component.

[0054] Note that, in the above description, all of the amplitude and phase are acquired from the electrical signals represented by (Expression 3) to (Equation 5) by performing Hilbert transform in digital signal processing but that, for example, an optical 90° hybrid may be used instead of the balanced optical reception units 12 and 13. Alternatively, an electrical 90° hybrid may be located following the balanced optical reception units 12 and 13, and the amplitude and phase may be extracted by an analog signal circuit.

[0055] Additionally, the synthetic AC signal represented by (Equation 5) may be synthesized by digital signal processing instead of mixing in the analog signal circuit. In this case, mixer 18 and filter 17 in FIG. 3 can be omitted.

[0056] Furthermore, the above description involves the use of the DC component of the synthetic AC signal. On the other hand, in a case where the signal component centered at the double frequency of the frequency shift Δω.sub.A is extracted by the high pass filter, Δω.sub.A can be monitored through this signal component. By applying the monitored frequency shift Δω.sub.A to the processing of (Equation 7) instead of assuming the fixed frequency shift Δω.sub.A, the monitored Δω.sub.A can be used for compensation for Δω.sub.A in a case described below. The case described below refers to a case where frequency drift is present in the optical frequency control unit 3, with Δω.sub.A changing slowly. In this case, it should be noted that the band of the digitization unit 14 needs to be 4Δω.sub.A or more in order to correctly acquire the synthetic AC signal in (Equation 5).

[0057] FIG. 5 illustrates a calculation processing flow of the calculation processing unit of the distributed optical fiber vibration measurement device according to the first embodiment.

[0058] After the digitization unit 14 samples the electrical signals (S301), the processing in (Expression 3) to (Equation 7) results in creation of amplitude and phase distribution data for each AC component (S302) and creation, from the synthetic AC component, of distribution data regarding the phase constant difference ϕ.sub.n (S304).

[0059] After the processing in S302, the distributed optical fiber vibration measurement device compares the amplitudes of the two AC components, selects the AC component with the larger amplitude, and arranges the phase data over time (S303). In S303, due to fading, which of the AC components is selected at each point of the optical fiber under test is random, but selection of the AC component with the larger amplitude results in creation of phase data from a signal with a higher signal-to-noise ratio.

[0060] After the processing in the S303 and S304, the distributed optical fiber vibration measurement device computes the phase change amount D.sub.n(t) at the same point of the optical fiber under test, i.e., at the same round-trip propagation delay time (S305). Specifically, the distributed optical fiber vibration measurement device computes a difference from the phase obtained during the preceding incidence of a test light pulse. However, in this computation, in a case where the AC component corresponding to the original data for the phase obtained during the (n−1)-th incidence differs from the AC component of the phase obtained during the n-th incidence, the phase constant difference computed in S304 is present between the phases. Thus, the distributed optical fiber vibration measurement device further executes processing for canceling the phase constant difference from the phase change amount. Note that FIG. 5 illustrates equations for a case where ϕ.sub.1,2=θ.sub.1−θ.sub.2 is used as the phase constant difference ϕ.sub.n but that, in a case where ϕ.sub.2,1=θ.sub.2−θ.sub.1 is used as the phase constant difference ϕ.sub.n, a sign for addition or subtraction of ϕ.sub.n may be reversed. Furthermore, in S306, from the phase change amount D.sub.n(t) thus obtained, computation is made of the difference between any two adjacent points of the optical fiber under test, i.e., between two amounts of change in phase D.sub.n(t) having a very small round-trip propagation delay time difference Δt, with a phase difference change ΔD.sub.n(t) in any section eventually output. The phase difference change ΔD.sub.n(t) obtained by causing the repeated incidence of the test light pulse corresponds to a temporal change in optical path length difference between any two adjacent points of the optical fiber under test, that is, a temporal change in the amount of strain, and ΔD.sub.n(t) corresponds to vibration.

[0061] Note that, in the present embodiment, the processing in S306 is executed after the processing in S305 as the calculation processing step but that the processing in S305 may be executed after the processing in S306. Specifically, in this method, after the phase data with a high signal-to-noise ratio is arranged in S303, the phase difference between any two points is computed in S306, and a temporal change in phase difference is computed in S305. This can be easily understood because, in steps S305 and S306, the values for the phase and the phase constant difference are simply subtracted.

[0062] As described above, by acquiring two different AC components from the two optical fibers under test and determining and compensating for the phase constant difference from the synthetic AC component, data with a higher signal-to-noise ratio can be selected, enabling accurate vibration measurement. Here, the phase constant difference is a phase constant difference occurring when the AC component used as the original signal from which the phase is computed is different from the signal from which the preceding phase is computed. Note that the description of the present example relates to the two optical fibers under test, but with a reception configuration for acquiring a plurality of electrical signals, two or more optical fibers under test can be used to implement the present embodiment on the assumption that the number of combinations of synthetic AC signals is similarly increased.

Second Embodiment

[0063] FIG. 6 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a second embodiment. The distributed optical fiber vibration measurement device of the present embodiment uses an optical detection unit to heterodyne-detect interference between rays of backscattered light obtained by causing test light to enter optical fibers under test, and performs this detection on two different optical fibers under test to acquire AC components of two electrical signals. Here, the above-described distributed optical fiber vibration measurement device causes two pulses having different optical frequencies and a time difference to enter the optical fibers under test as the above-described test light. The distributed optical fiber vibration measurement device of the present embodiment does not require the local light used in the first embodiment. In other words, the present embodiment differs from the first embodiment in a unit for acquiring AC components of two electrical signals, and only differences in configuration and signal processing involved in the above-described change will be described below.

[0064] Output light from the light source 1 emitting coherent light is incident on the optical frequency control unit 3 as test light. The optical frequency control unit 3 modulates the optical frequency in two stages at predetermined time intervals. The optical frequency difference at this time is Δω.sub.A. For the test light, the light pulsing unit 5 forms two pulses having a first optical frequency ω and a second optical frequency ω+Δω.sub.A, and causes the pulses to enter the two optical fibers under test 10 via the optical amplifier 6 and the branching element 8.

[0065] Rays of backscattered light generated in the optical fibers under test 10 travel through the optical circulators 9 and are pre-amplified by an optical amplifier 21 and an optical amplifier 22, and the rays of light then travel through optical filters 23 and 24 for removing spontaneous emissions from the optical amplifiers 21 and 22 and are converted into electrical signals by optical detection units 25 and 26. The optical detection units 25 and 26 multiplexes two rays of backscattered light with different optical frequencies for heterodyne detection, and the distributed optical fiber vibration measurement device can acquire AC components of two electrical signals from the optical detection units 25 and 26.

[0066] Two outputs from the optical detection units 25 and 26 are partially branched, and the partial outputs are input to a mixer 28, in which the partial outputs are multiplexed. Then, high frequency components are removed from the multiplexing results via a filter 27, and the resultant signal is converted into digital signal by the digitization unit 14. Then, the digital signal is subjected to signal processing by the calculation processing unit 15. Additionally, similarly to being branched and input to the mixer 28, the two output signals from the optical detection units 25 and 26 are input to the digitization unit 14, and the signals are converted into digital signals by the digitization unit 14. The digital signals are subjected to signal processing in the calculation processing unit 15.

[0067] FIG. 7 illustrates a test light pulse delivery sequence according to the present embodiment. Two test light pulses are repeatedly delivered with a predetermined time difference t.sub.d, the test light pulses having the optical frequency ω of the light from the light source 1 and the optical frequency ω+Δω.sub.A resulting from a frequency shift by the optical frequency control unit 3. The two test light pulses have the same pulse width W. The repetition time of the test light pulse is set equal to or longer than the round-trip propagation delay time of the maximum length of the optical fiber under test 10.

[0068] Here, AC components I.sub.1(t) and I.sub.2(t) of beat electrical signals output from the optical detection units 25 and 26 are respectively given by (Expression 8) and (Expression 9) below, with respect to the output time t and the optical frequency ω of the coherent light from the light source 1, the frequency shift Δω.sub.A, and the pulse width W of the test light pulse.

[00005]$\begin{array}{cc}\mathrm{Math}.8& \\ I_1\left(t\right)\propto 2E_S^2\underset{i=1}{\overset{N}{.Math.}}\underset{j=1}{\overset{N}{.Math.}}{a}_i{a}_j\cos \left\{{\mathrm{\Delta \omega }}_{A}t+\omega {\tau }_i-\left(\omega +{\mathrm{\Delta \omega }}_A\right){\tau }_j+\theta \left(t-{\tau }_i\right)-\theta \left(t-{\tau }_j\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right\}.Math.\mathrm{rect}\left(\frac{t-{\tau }_i}{W}\right).Math.\mathrm{rect}\left(\frac{t-{\tau }_j}{W}\right)=2E_S^2{R}_1\left(t\right)\cos \left({\mathrm{\Delta \omega }}_{A}t+{\theta }_1+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right)& \left(\mathrm{Expression}8\right)\\ \mathrm{Math}.9& \\ I_2\left(t\right)\propto 2E_S^2\underset{i=1}{\overset{N}{.Math.}}\underset{j=1}{\overset{N}{.Math.}}{b}_i{b}_j\cos \left\{{\mathrm{\delta \omega }}_{A}t+{\mathrm{\omega \tau }}_i-\left(\omega +\Delta {\omega }_A\right){\tau }_j+\theta \left(t-{\tau }_i\right)-\theta \left(t-{\tau }_j\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right\}.Math.\mathrm{rect}\left(\frac{t-{\tau }_i}{W}\right).Math.\mathrm{rect}\left(\frac{t-{\tau }_j}{W}\right)=2E_S^2{R}_2\left(t\right)\cos \left({\mathrm{\Delta \omega }}_{A}t+{\theta }_2+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right)& \left(\mathrm{Expression}9\right)\end{array}$

[0069] Here, E.sub.S is the electric field amplitude of the test light, N is the number of scatterers present in the length direction of each of the optical fibers under test, i (1≤i≤N) is the number of each scatterer, a.sub.i, b.sub.i is the reflectance of the scatterer i, τ.sub.i is the round-trip propagation delay time from the incident end of the optical fiber under test to the scatterer i, θ(t) is the initial phase of the test light, and θ.sub.FUT(τ) is a phase change caused by vibration applied to the optical fiber under test and is an object to be measured. R.sub.1(t) and R.sub.2(t) are the backscattering coefficient distributions of the respective optical fibers under test, and according to the assumed conditions, θ.sub.FUT(τ) is the same between two optical fibers under test. Here, when the time difference between the two test light pulses incident on the respective optical fibers under test is sufficiently shorter than the coherence length of the light source, θ(t−τ.sub.i)−θ(t−τ.sub.j), which means phase noise from the laser, is very small and negligible. θ.sub.1 and θ.sub.2 are each a phase constant attributed to fading, and conditions for interference caused by fading, in other words, scatterers forming the optical fiber, vary among the optical fibers under test, and thus the phase constants differ from each other.

[0070] As seen in (Expression 8) and (Expression 9), the form of the expression for each AC component is substantially the same as that in the first embodiment except that, with the phase term focused on, the term meaning phase noise from the laser is substantially negligible at this stage.

[0071] In the stage subsequent to the optical detection units 25 and 26, the electrical signals resulting from branching are synthesized in the mixer 28, and then high frequency components are removed from the signal by the filter 27. The synthetic AC signal resulting from synthesis in the mixer 28 is given by the product of the AC components, and is as described below.

Math. 10

I.sub.1(t)×I.sub.2(t)=4E.sub.L.sup.2R.sub.1(t)R.sub.2(t)[cos{2(Δω.sub.At+θ.sub.out(τ))+θ.sub.1+θ.sub.2}+cos(ϕ.sub.1,2)]  (Equation 10)

[0072] The synthetic AC signal can also be treated as in the first embodiment except that, with the phase term focused on, the term meaning phase noise from the laser is substantially negligible at this stage. Accordingly, the calculation processing flow in the calculation processing unit of the present embodiment is the same as that in the first embodiment except that, with reference to FIG. 5, the calculation step in S306 is unnecessary. This is because, in the present embodiment, the time difference is set for the two test light pulses with different frequencies incident on the respective optical fibers under test, and thus a phase change in any section of each optical fiber under test corresponding to the time difference is treated as the signal θ.sub.FUT(τ) in each of the AC components computed in (Expression 8) and (Expression 9).

[0073] As described above, accurate vibration measurement can be implemented by the present embodiment, which differs from the first embodiment in that the distributed optical fiber vibration measurement device causes test light pulses with different frequencies to be incident with a time difference instead of the local light and that the interference between the resultant plurality of rays of backscattered light is measured as an AC component.

Third Embodiment

[0074] FIG. 8 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device of a third embodiment. The distributed optical fiber vibration measurement device of the present embodiment branches light from the light source into rays of light used as the test light and the local light, and shifts the frequency of the test light with respect to the frequency of the local light by a plurality of different values to generate frequency-shifted test light at predetermined time intervals (frequency-multiplexed test light). Furthermore, the distributed optical fiber vibration measurement device of the present embodiment causes the test light to enter one optical fiber under test to obtain backscattered light, and multiplexes the backscattered light with the local light for heterodyne detection. Furthermore, the distributed optical fiber vibration measurement device of the present embodiment performs digitization processing on the electrical signal obtained and then uses digital signal processing to perform frequency separation to acquire two AC components. The present embodiment differs from the first embodiment and the second embodiment in that instead of using two optical fibers under test or two balanced optical reception units (or two optical detection units), the present embodiment executes digital signal processing to separate frequency-multiplexed test light into rays of light. Only differences in configuration and signal processing involved in the above-described change will be described below.

[0075] Output light from the light source 1 emitting coherent light is branched into two lines by the branching element 2. One ray of the branched light is used as local light, whereas the other is incident on an optical frequency control unit 3 as test light. The optical frequency control unit 3 modulates the optical frequency in two stages (shifts the frequency by different values) at predetermined time intervals. For the test light, the light pulsing unit 5 forms two test light pulses including a first test light pulse having an optical frequency ω+ω.sub.1 and a second test light pulse having an optical frequency ω+ω.sub.2, and causes the test light pulses to enter the optical fiber under test 10 via the optical amplifier 6.

[0076] The backscattered light generated in the optical fiber under test 10 travels through the optical circulator 9 and is converted into an electrical signal by the balanced optical reception unit 12. The electrical signal is converted into a digital signal by the digitization unit 14, and the digital signal is subjected to signal processing by the calculation processing unit 15. The balanced optical reception unit 12 causes test light (first test light pulse and second test light pulse) to enter one optical fiber under test to obtain backscattered light, and multiplexes the backscattered light with the local light for heterodyne detection to converts the backscattered light into an electrical signal.

[0077] FIG. 9 illustrates a test light pulse delivery sequence according to the present embodiment. Light from the light source 1 is frequency-shifted by the optical frequency control unit 3, and two test light pulses having the optical frequencies ω+ω.sub.1 and ω+ω.sub.2 are repeatedly delivered with a predetermined time difference t.sub.d. The two test light pulses (third test light pulse and fourth test light pulse) has the same pulse width W. The repetition time of the test light pulse is set equal to or longer than the round-trip propagation delay time of the maximum length of the optical fiber under test 10.

[0078] Here, the beat electrical signal output from the balanced optical reception unit 12 is digitized by the digitization unit 14, and then in the calculation processing unit 15, the resultant signal is separated into two signals having beat frequencies ω1 and ω2 by using a digital filter. The AC components I.sub.1(t) and I.sub.2(t) of the two signals resulting from the separation are respectively given by (Expression 11) and (Expression 12) below, with respect to the output time t and the optical frequency ω of the coherent light from the light source 1, frequency shifts ω.sub.1 and ω.sub.2, and the pulse width W of the test light pulse.

[00006]$\begin{array}{cc}\mathrm{Math}.11& \\ I_1\left(t\right)\propto 2E_L{E}_S\underset{i=1}{\overset{N}{.Math.}}{a}_i\cos \left\{{\omega }_{1}t-{\omega }_1{\tau }_i+\theta \left(t-{\tau }_i\right)-\theta \left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right\}.Math.\mathrm{rect}\left(\frac{t-{\tau }_i}{W}\right)=2E_L{E}_S{R}_1\left(t\right)\cos \left({\omega }_{1}t+{\theta }_1+{\mathrm{\Delta \theta }}_L\left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right)& \left(\mathrm{Expression}11\right)\\ \mathrm{Math}.12& \\ I_2\left(t\right)\propto 2E_L{E}_S\underset{i=1}{\overset{N}{.Math.}}{a}_i\cos \left\{{\omega }_2\left(t-t_d\right)-{\omega }_2{\tau }_i+\theta \left(t-t_d-{\tau }_i\right)-\theta \left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau -t_d\right)\right\}.Math.\mathrm{rect}\left(\frac{t-t_d-{\tau }_i}{W}\right)-2E_L{E}_S{R}_2\left(t-t_d\right)\cos \left({\omega }_{2}t+{\theta }_2+{\mathrm{\Delta \theta }}_L^i\left(t\right)+{\theta }_{\mathrm{FUT}}\left(\tau -t_d\right)\right)& \left(\mathrm{Expression}12\right)\end{array}$

[0079] Here, E.sub.L and E.sub.S are respectively the electric field amplitudes of the local light and the test light, N is the number of scatterers present in the length direction of each of the optical fibers under test, i (1≤i≤N) is the number of each scatterer, a.sub.i is the reflectance of the scatterer i, τ.sub.i is the round-trip propagation delay time from the incident end of the optical fiber under test to the scatterer i, θ(t) is the initial phase of the test light, and θ.sub.FUT(τ) is a phase change caused by vibration applied to the optical fiber under test and is an object to be measured. R.sub.1(t) and R.sub.2(t) are respectively the backscattering coefficient distributions of the optical fibers under test at ω.sub.1 and ω.sub.2. Δθ.sub.L(t) is the initial phase difference of light between the time of output of the test light pulse and the time of interference with the local light, meaning phase noise from the laser, and is eventually negligible as in the first and second embodiments.

[0080] Here, the two AC components according to the present embodiment mean the results of the interference between the backscattered light resulting from the first test light pulse and local light and the results of the interference between the backscattered light resulting from the second test light pulse, that is incident with the time difference t.sub.d, and local light. Accordingly, the above-described two AC components are signals containing the influences of vibrations differing from each other in time by t.sub.d, i.e., θ.sub.FUT(τ) and θ.sub.FUT(τ−t.sub.d), respectively. θ.sub.1 and θ.sub.2 are each a phase constant attributed to fading, and are different constants because the first test light pulse and the second test light pulse have different optical frequencies.

[0081] Now, assuming that T is a measurable vibration frequency, in a case where t.sub.d<<1/T holds true, i.e., in a case where the time difference t.sub.d is very small with respect to the vibration period, it is possible to assume that θ.sub.FUT(τ)≈θ.sub.FUT(τ−t.sub.d). As a specific example, assuming that the measurable maximum vibration frequency is 20 kHz, which is the upper limit of the human audible range, and t.sub.d=200 ns, the time difference t.sub.d corresponds to one-250th of a vibration period of 50 μs, and thus even the maximum value is approximately 0.03 rad in terms of the phase change amount, which is small enough to be negligible.

[0082] The AC components indicated in (Expression 11) and (Expression 12) are synthesized by digital signal processing in the calculation processing unit 15. This synthetic AC signal is corrected with the time shift of the AC component I.sub.2(t) by the incident time difference t.sub.d taken into account, and is further modified as follows, by using θ.sub.FUT(τ)≈θ.sub.FUT(τ−t.sub.d).

Math. 13

I.sub.1(t)×I.sub.2(t+t.sub.d)=4E.sub.L.sup.2E.sub.S.sup.2R.sub.1(t)R.sub.2(t)[cos{(ω.sub.1+Ω.sub.2)t+2(Δθ.sub.L(t)+θ.sub.OUT(t))+θ.sub.1+θ.sub.2}+cos((ω.sub.1−ω.sub.2)t+ϕ.sub.1,2)]   (Equation 13)

[0083] (Equation 13) indicates that the synthetic AC signal represented by this equation includes a frequency component of the sum ω.sub.1+ω.sub.2 of the optical frequency shifts applied to the first test light pulse and the second test light pulse, and a frequency component of the difference ω.sub.1−ω.sub.2 between the optical frequency shifts applied to the first test light pulse and the second test light pulse. By using a low-pass filter, the frequency component of the difference ω.sub.1−ω.sub.2 is extracted. As is the case with the first embodiment and the second embodiment, the distributed optical fiber vibration measurement device can extract the phase term using an orthogonal component shifted by 90 degrees and execute subtraction of the known optical frequency shifts ω.sub.1 and ω.sub.2 applied to the first test light pulse and the second test light pulse, to learn the phase constant difference ϕ.sub.1,2. Subsequently, the distributed optical fiber vibration measurement device can execute processing as is the case with the first embodiment.

[0084] Here, when a reception band is designed, given that the signal band of the test light pulse is 2/W, the condition for the design needs to be such that, for the optical frequencies of the first test light pulse and the second test light pulse, |ω.sub.1−ω.sub.2|≥1/W. Additionally, the band for the balanced optical reception unit 12 and the digitization unit 14 needs to be 2(ω.sub.1+ω.sub.2) or more on the basis of the Nyquist theorem.

[0085] As described above, in the present embodiment, the distributed optical fiber vibration measurement device produces frequency-multiplexed test light, causes rays of test light with different frequencies to simultaneously enter one optical fiber under test to obtain backscattered light, and multiplexes the backscattered light with the local light for heterodyne detection. In addition, the distributed optical fiber vibration measurement device performs frequency separation using digital signal processing executed by one optical detection unit (balanced optical reception unit 12), acquiring two AC signals. This enables the present embodiment, which is different from the first embodiment, to implement accurate vibration measurement.

[0086] Note that, in the present embodiment as well, under the condition that the time difference t.sub.d is very small compared to the vibration period, more than two frequencies may be multiplexed and that, in this case, assuming that the number of synthetic AC signals combined is also increased, multiplexing of even three or more frequencies can be implemented.

[0087] The present embodiment can be implemented without the use of two or more optical fibers under test or optical detection units, allowing the configuration to be simplified compared to the first embodiment and the second embodiment, thus enabling an economical device to be provided.

Fourth Embodiment

[0088] FIG. 10 is a block configuration diagram illustrating a distributed optical fiber vibration measurement device according to a fourth embodiment. The distributed optical fiber vibration measurement device of the present embodiment uses test light frequency-coded at predetermined time intervals, causes the test light to enter the optical fiber under test to obtain rays of backscattered light, uses the optical detection unit to receive interference between the rays of backscattered light, and subjects the resultant signal to digitization processing and then to digital signal processing to frequency-separate the signal into any two AC components. A difference from the third embodiment is that, instead of the signals corresponding to the first test light pulse and the second test light pulse of the third embodiment, test light pulses each frequency-multiplexed and coded are used. Specifically, the present embodiment uses a unit for acquiring the AC components of two electrical signals which differs from the corresponding unit in the third embodiment, and only differences in configuration and signal processing involved in the above-described change will be described below.

[0089] Output light from the light source 1 emitting coherent light is incident on the optical frequency control unit 3 as test light. The optical frequency control unit 3 modulates the optical frequency in four stages at predetermined time intervals. For the test light, the light pulsing unit 5 forms one test light pulse (one test light pulse in which a third test light pulse, a fourth test light pulse, a fifth test light pulse, and a sixth test light pulse are continuously formed) including the frequency component modulated in four stages, and is incident on the optical fiber under test 10 via the optical amplifier 6. Backscattered light generated in the optical fiber under test 10 travels through the optical circulator 9 and is pre-amplified by an optical amplifier 30, and the light then travels through an optical filter 31 for removing spontaneous emission from the optical amplifier 30 and is converted into an electrical signal by an optical detection unit 32. The electrical signal is converted into a digital signal by the digitization unit 14, and the digital signal is subjected to signal processing by the calculation processing unit 15.

[0090] FIG. 11 is a diagram illustrating a test light pulse delivery sequence according to the present embodiment. The light from the light source 1 is modulated in four stages ω.sub.A, ω.sub.B, ω.sub.A+ω.sub.1, and ω.sub.B+ω.sub.2 at predetermined time intervals W. The repetition time of the test light pulse is set equal to or longer than the round-trip propagation delay time of the maximum length of the optical fiber under test 10.

[0091] Here, given that the signal band of each frequency component of the test light pulse is 2/W, the condition for the design needs to be such that, for the reception band, |ω.sub.1−ω.sub.2|≥1/W. Additionally, the band for the optical detection unit 32 and the digitization unit 14 needs to be 2(ω.sub.1+ω.sub.2) or more on the basis of the Nyquist theorem. Furthermore, the difference between OA and OB needs to be larger than the reception band for the optical detection unit 32 and the digitization unit 14. With the design as described above, the optical detection unit 32 detects the beat frequency components ω.sub.1 (corresponding to the difference frequency component between ω.sub.A and ω.sub.A+ω.sub.1) and ω.sub.2 (corresponding to the difference frequency component between ω.sub.B and ω.sub.B+ω.sub.2).

[0092] Here, the beat electrical signal is output from the optical detection unit 32, and is digitized by the digitization unit 14. Subsequently, the calculation processing unit 15 employs a digital filter to separate the beat electrical signal into two signals having the beat frequencies ω.sub.1 and ω.sub.2. The AC components I.sub.1(t) and I.sub.2(t) of the two signals resulting from the separation at this time are respectively given by (Expression 14) and (Expression 15) below, with respect to the output time t of the coherent light from the light source.

[00007]$\begin{array}{cc}\mathrm{Math}.14& \\ I_1\left(t\right)\propto 2E_S^2\underset{i=1}{\overset{N}{.Math.}}\underset{j=1}{\overset{N}{.Math.}}{a}_i{a}_j\cos \left\{{\omega }_{1}t+{\omega }_A{\tau }_i-\left({\omega }_1+{\omega }_A\right){\tau }_j+\theta \left(t-{\tau }_i\right)-\theta \left(t-{\tau }_j\right)+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right\}.Math.\mathrm{rect}\left(\frac{t-{\tau }_i}{W}\right).Math.\mathrm{rect}\left(\frac{t-{\tau }_j}{W}\right)=2E_S^2{R}_1\left(t\right)\cos \left({\omega }_{1}t+{\theta }_1+{\theta }_{\mathrm{FUT}}\left(\tau \right)\right)& \left(\mathrm{Expression}14\right)\\ \mathrm{Math}.15& \\ I_2\left(t\right)\propto 2E_S^2\underset{i=1}{\overset{N}{.Math.}}\underset{j=1}{\overset{N}{.Math.}}{a}_i{a}_j\cos \left\{{\omega }_2\left(t-W\right)+{\omega }_B{\tau }_i-\left({\omega }_2+{\omega }_B\right){\tau }_j+\theta \left(t-W-{\tau }_i\right)-\theta \left(t-W-{\tau }_j\right)+{\theta }_{\mathrm{FUT}}\left(\tau -W\right)\right\}.Math.\mathrm{rect}\left(\frac{t-W-{\tau }_i}{W}\right).Math.\mathrm{rect}\left(\frac{t-W-{\tau }_j}{W}\right)=2E_S^2{R}_2\left(t\right)\cos \left({\omega }_{2}t+{\theta }_2+{\theta }_{\mathrm{FUT}}\left(\tau -W\right)\right)& \left(\mathrm{Expression}15\right)\end{array}$

[0093] Here, E.sub.S is the electric field amplitude of the test light, N is the number of scatterers present in the length direction of each of the optical fibers under test, i (1≤i≤N) is the number of each scatterer, a.sub.i is the reflectance of the scatterer i, τ.sub.i is the round-trip propagation delay time from the incident end of the optical fiber under test to the scatterer i, θ(t) is the initial phase of the test light, and θ.sub.FUT(τ) is a phase change caused by vibration applied to the optical fiber under test and is an object to be measured. R.sub.1(t) and R.sub.2(t) are the backscattering coefficient distributions of the respective optical fibers under test at OA and OB. Here, when the time difference W between the frequency components OA and OB incident on the optical fiber under test is sufficiently shorter than the coherence length of the light source, θ(t−τ.sub.i)−θ(t−τ.sub.j) and θ(t−W−τ.sub.i)−θ(t−W−τ.sub.j), which mean phase noise from the laser, are very small and negligible. θ.sub.1 and θ.sub.2 are each a phase constant attributed to fading, and differ from each other due to the different optical frequencies of the two test light pulses.

[0094] Now, assuming that T is a measurable vibration frequency, in a case where W<<1/T holds true, i.e., in a case where the time difference W is very small with respect to the vibration period, it is possible to assume that θ.sub.FUT(τ)≈θ.sub.FUT(τ−W). As a specific example, assuming that the measurable maximum vibration frequency is 20 kHz, which is the upper limit of the human audible range, and W=200 ns, the time difference W corresponds to one-250th of a vibration period of 50 μs, and thus even the maximum value is approximately 0.03 rad in terms of the phase change amount, which is small enough to be negligible.

[0095] The AC components indicated in (Expression 14) and (Expression 15) are synthesized by digital signal processing in the calculation processing unit 15. This synthetic AC signal is corrected with the time shift of the AC component i.sub.2(t) by the incident time difference W taken into account, and is further modified as represented in (Equation 16) below, by using θ.sub.FUT(τ)≈θ.sub.FUT(τ−W).

Math. 16

I.sub.1(t)×I.sub.2(t+W)=2E.sub.S.sup.AR.sub.1(t)R.sub.2(t)[cos{(ω.sub.1+ω.sub.2)t+2θ.sub.FUT(τ)+θ.sub.1+θ.sub.2}cos((ω.sub.1−ω.sub.2)t+ϕ.sub.1,2)]  (Equation 16)

[0096] (Equation 16) indicates that the synthetic AC signal represented by this equation includes a frequency component of the sum ω.sup.1+ω.sub.2 of the optical frequency shifts applied to the two test light pulses, and a frequency component of the difference ω.sub.1−ω.sub.2 between the optical frequency shifts applied to the two test light pulses. By using a low-pass filter for the above-described signals, the frequency component of the difference ω.sub.1−ω.sub.2 is extracted. As is the case with the first embodiment and the second embodiment, the distributed optical fiber vibration measurement device can extract the phase term using an orthogonal component shifted by 90 degrees and execute subtraction of the known optical frequency shifts ω.sub.1 and ω.sub.2 to learn the phase constant difference ϕ.sub.1,2. Subsequently, the distributed optical fiber vibration measurement device can execute processing as is the case with the second embodiment.

[0097] In spite of the use of different acquisition units, the distributed optical fiber vibration measurement device described in the first to fourth embodiments can learn, from the synthetic AC component, the phase constant difference attributed to fading by acquiring the two AC components under a condition as described below. This means that, even in a case where the distributed optical fiber vibration measurement device compares the amplitudes of the two AC components and selects the AC component with a higher signal-to-noise ratio, thus changing the selected AC component itself, a difference in phase constant difference can be corrected when a change in the time direction is computed. Thus, the two AC components can be used to implement accurate vibration measurement that is less susceptible to fading noise. Here, the above-described condition described below is the condition that the optical frequencies of the test light pulses or the optical fibers under test or both vary. Additionally, the above-described change in the time direction refers to a temporal change in the phase of the AC component selected for each repeated incidence of the test light pulse.

[0098] Additionally, in the embodiment in which the phase is computed from the interference of backscattered light, if optical frequency drift in the light source or the optical frequency control unit is present, the optical frequency drift becomes a direct factor that increases phase errors in a process described below. However, by taking a measure for monitoring the beat frequency without removing high frequency components of the synthetic AC component as in the first embodiment, the optical frequency drift can be corrected in other embodiments. Here, the above-described process described below is a process in which the distributed optical fiber vibration measurement device computes the phase from the AC component or the synthetic AC component, and a process of computing a phase change with respect to time.

[0099] Note that, in the first embodiment to the fourth embodiment, the calculation processing unit 15 can be implemented by a computer and a program and that the program can be recorded in a recording medium or provided through a network.

[0100] In addition, various disclosures can be formed by appropriately combining a plurality of constituent elements disclosed in the above embodiments. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, components of different embodiments may be appropriately combined.

REFERENCE SIGNS LIST

[0101] 1 Light source [0102] 2 Branching element [0103] 3 Optical frequency control unit [0104] 4 Driving unit [0105] 5 Light pulsing unit [0106] 6 Optical amplifier [0107] 7 Pulse signal generation unit [0108] 8 Branching element [0109] 9 Optical circulator [0110] 10 Optical fiber under test [0111] 11 Branching element [0112] 12, 13 Balanced optical reception unit [0113] 14 Digitization unit [0114] 15 Calculation processing unit [0115] 16 Display unit [0116] 17, 27 Filter [0117] 18, 28 Mixer [0118] 23, 24 Optical filter [0119] 25, 26 Optical Detection unit [0120] 21, 22, 30 Optical Amplifier [0121] 31 Optical filter [0122] 32 Optical detection unit