Patent Description:
With an optical fiber sensing technology, it is possible to measure and detect various physical and chemical quantities such as temperature, strain, displacement, vibration, and pressure, and there are advantages such as a small diameter, low weight, flexibility, durability, and electromagnetic noise tolerance that electric sensors do not have. This can also be used for a distribution type sensor capable of measurement at a specific or freely selected position along an optical fiber. Various physical quantities of a measurement target can be measured from transmitted light, reflected light, backscattered light, and the like obtained by causing probe light to pass through a sensing fiber disposed in the measurement target. It is also possible to measure, in a distribution manner, of a temporal change (dynamic strain) of strain applied to the sensing fiber due to mechanical vibration of the measurement target at different positions in a longitudinal direction of the sensing fiber.

As a vibration distribution measurement technology, optical frequency domain reflectometry (OFDR) disclosed in NPL <NUM> (hereinafter, related art) is known. In OFDR, frequency swept light from a frequency swept light source is input to a sensing fiber as probe light. A beat signal between backscattered light returning from the sensing fiber and local light split from the frequency swept light is taken. A beat frequency of the beat signal corresponds to a distance to scatterers on the sensing fiber. A relationship between the beat frequency and the distance is used to measure a waveform of the backscattered light at each distance from an incidence end of the sensing fiber, that is, a distribution waveform of the backscattered light. A freely selected section of the distribution waveform is cropped, and a section of the sensing fiber corresponding thereto is defined as a strain sensor section. The Fourier transform (a square of an absolute value) of the distribution waveform in the strain sensor section indicates an optical spectrum of the backscattered light in the strain sensor section.

In OFDR, frequency sweep is repeated in a light source and probe light is repeatedly input to the sensing fiber, making it possible to capture a vibration phenomenon itself of an observation target, that is, a time waveform of vibration. For dynamic strain measurement, in the related art, frequency swept probe light is incident on the sensing fiber at a frequency higher than a frequency of vibration that can occur in the sensing fiber. Frequency modulation caused by vibration to light propagating through the sensing fiber causes a frequency offset in a beat signal observed by an optical receiver. In OFDR, the distance of the scatterers to be measured fluctuates due to the frequency offset caused by vibration that is a measurement target. Thus, the related art has proposed a scheme to compensate for the fluctuation of the distance measured through digital signal processing and accurately measuring the vibration at a designated position from a probe light incident point.

Further prior art documents are <CIT> showing a method and device for measuring difference in propagation delay time between spatial channels, <CIT> showing a dispersion correction in optical frequency-domain reflectometry and <CIT> showing a vibration distribution measuring system, vibration waveform analysis method, vibration waveform analyzing device, and analyzing program.

A vibration distribution measurement device that accurately measures vibration at a designated position without using signal processing to compensate for distance fluctuation is provided. An index of tolerance of vibration distribution measurement to distance fluctuation (distance offset), which further simplifies the vibration distribution measurement device, is also clarified.

An embodiment of the present disclosure is a device and a method according to the respective claims.

A vibration distribution measurement device that accurately measures dynamic strain without increasing a load of signal processing is provided.

A vibration distribution measurement device of the present disclosure provides measurement of a vibration distribution with a reduced load of signal processing on a sensing fiber. The vibration distribution measurement device of the present disclosure is based on OFDR and can measure dynamic strain at a freely selected position on the sensing fiber. The dynamic strain that is a measurement target means a phenomenon in which the amount of strain fluctuates with time, and includes vibration of which the amount of strain periodically fluctuates in a sinusoidal manner and aperiodic vibration of which the amount of strain randomly changes with time. That is, this is a concept including vibration including one or a plurality of frequency components. This also includes detection of one change in amount of strain, for example, when an object is placed on the sensing fiber and pressure changes at that moment.

As will be described below, dynamic strain is measured by detecting a spectrum change (shift) with respect to an optical spectrum obtained in reference measurement. Thus, for the dynamic strain in the vibration distribution measurement device of the present disclosure, all phenomena in which the amount of strain (intensity) changes with time are detection targets. Examples of a use pattern include measuring cable vibration due to wind using a sensing fiber to measure wind pressure, and detecting and reproducing sound generated around a sensing fiber for a function of a sound collecting microphone. The following vibration distribution measurement device and method of the present disclosure can be applied not only to the above-described use patterns but also to various fields.

In the related art, fluctuation of a measurement distance due to a frequency offset of a beat signal caused by vibration that is a measurement target is corrected through digital signal processing. In each sweep of an optical frequency, a distance offset at a designated position on a sensing fiber is determined and strain at a corrected distance is measured, but signal processing for determining the distance offset requires time and arithmetic power. When the vibration is measured in a distribution manner at many measurement points of the entire sensing fiber or when dynamic strain is measured over a long period of time, a load of signal processing in a vibration distribution measurement device increases. In view of the problem of the related art, the present disclosure provides a device that accurately measures vibration at a designated position without using digital signal processing to compensate for distance fluctuation.

<FIG> is a diagram illustrating a configuration of a vibration distribution measurement device according to the present disclosure. A vibration waveform analysis device <NUM> broadly divided into an optical circuit unit that supplies the frequency swept light and receives signal light, a sensing fiber <NUM>, and an optical reception and analysis unit that converts the signal light from the sensing fiber <NUM> into a beat signal and analyzes the beat signal. The optical circuit unit includes a frequency swept light source <NUM>, an optical splitter <NUM>, an optical circulator <NUM>, and an optical combiner <NUM>. The frequency swept light source <NUM> is, for example, a laser capable of the frequency sweep at high speed, and receives a trigger signal from a trigger source <NUM> and repeatedly outputs the frequency swept light at a predetermined sweep period. The frequency swept light is split into two by the optical splitter <NUM>, and one becomes probe light Lprobe that is supplied to the optical circulator <NUM>, and the other becomes the local light Llocal that is supplied to the optical combiner <NUM>.

The probe light Lprobe is incident on the sensing fiber <NUM> via the optical circulator <NUM> which is a directional coupling element, backscattered by scatterers in the sensing fiber <NUM>, and returned to the incidence side. The signal light Lsignal in which backscattered light beams Lbs from respective positions of the sensing fiber <NUM> are superimposed exits the sensing fiber <NUM> and propagates to the optical combiner <NUM> via the optical circulator <NUM>. Thus, the optical circuit unit repeatedly supplies the frequency swept light to the sensing fiber, receives backscattered light from the sensing fiber, and combines local light, which is a part of the frequency swept light, with the backscattered light.

The optical reception and analysis unit includes a balanced optical receiver <NUM>, an A/D converter <NUM>, an acquired data storage (memory) <NUM>, and an analysis unit <NUM>. The signal light Lsignal and the local light Llocal from the optical splitter <NUM> are combined by the optical combiner <NUM>, and the beat signal Sbeat between Lsignal and Llocal is output as an electrical signal by the balanced optical receiver <NUM>. The beat signal Sbeat is sampled by an A/D converter <NUM> and measured as a digital signal in synchronization with the trigger signal to the frequency swept light source <NUM>. The measured beat signal at each time is stored in the acquired data storage <NUM>, and the analysis unit <NUM> performs measurement of the vibration distribution and analysis of a vibration waveform. With the vibration distribution measurement device illustrated in <FIG>, a backscattered light spectrum in a freely selected section of the sensing fiber <NUM> can be measured, and the vibration distribution and the vibration waveform can be measured and analyzed.

An operation of the optical reception and analysis unit is summarized as follows. A beat signal is generated from the local light and the backscattered light, and Fourier transform is executed for the beat signal to obtain a backscattered light waveform. Next, a Fourier transform is executed in a freely selected spectrum analysis section of the backscattered light waveform to obtain an optical spectrum, and the amount of strain of the dynamic strain in a section of the sensing fiber corresponding to the spectrum analysis section is obtained based on the amount of spectrum shift between a spectrum obtained in the reference measurement and a spectrum obtained from the backscattered light for one sweep of the frequency swept light.

The "dynamic strain" is a temporal change in amount of strain in a sensing fiber section and is obtained from the amount of strain for each sweep of the frequency swept light that is repeatedly supplied.

The entire vibration distribution measurement device <NUM> in <FIG> can be regarded as a vibration distribution measurement system. Further, the analysis of the vibration time waveform of the dynamic strain can be performed as arithmetic (calculation) processing in an analysis unit <NUM> in which, for example, a central processing unit (CPU), a digital signal processor (DSP), or the like is used. The present disclosure also has an aspect as an analysis program that implements a method including processing steps to be described below. Hereinafter, the fluctuation of the measurement distance due to the frequency offset of the beat signal caused by a vibration that is a measurement target, which is a common problem in the related art and the present disclosure, will be described.

<FIG> is a model diagram illustrating a generation situation of the backscattered light in the sensing fiber assumed in the vibration distribution measurement system. In the sensing fiber model <NUM> of <FIG>, the probe light Lprobe incident from an input end of the measurement system at a left end (not illustrated) propagates in a first dynamic strain section <NUM> of the sensing fiber <NUM>. The probe light Lprobe that has passed through the first dynamic strain section <NUM> is backscattered by the scatterers constituting a second the dynamic strain section <NUM>, and propagates in a reverse direction through the first dynamic strain section <NUM> again, and returns to an input end of the sensing fiber. A situation in which the backscattered light Lbs is emitted from the input end as the signal light Lsignal and input to the optical circulator <NUM> in <FIG> is shown.

The emitted signal light Lsignal is subjected to phase modulation, as described below, due to the mechanical vibration in the first dynamic strain section <NUM>. According to the related art, the phase modulation to which the signal light Lsignal is subjected due to the mechanical vibration of the first dynamic strain section <NUM> will be described as follows.

<FIG> is a diagram illustrating an example of a waveform of phase modulation due to the mechanical vibration in the dynamic strain section. In <FIG>, a time section corresponding to one period Pmod of a phase modulation waveform generated by the mechanical vibration in the first dynamic strain section <NUM> illustrated in <FIG> is shown. In a time section of Pmod, N beat signal measurements are performed in measurement periods (Mes1 to MesN), and each of the measurement periods has a maximum of one sweep period of the frequency swept light.

A round-trip time for the probe light Lprobe to propagate in the first dynamic strain section <NUM>, be backscattered by the scatterer in the second the dynamic strain section <NUM>, and propagate in the first dynamic strain section <NUM> again, as illustrated in <FIG>, is considered. When this round-trip time is sufficiently shorter than the period Pmod of the phase modulation, a phase modulation waveform on an outward path and a phase modulation waveform on a return path can be regarded as being the same. Thus, the amount of phase modulation to which the light propagating back and forth in the first dynamic strain section <NUM> is subjected is twice that of the phase modulation in the outward path. An electrical field waveform Esignal(t) in one measurement of the signal light Lsignal propagating back and forth is expressed by Relationship (<NUM>) below. <NUM>] <MAT>.

Here, ν<NUM> indicates an initial frequency in the frequency sweep, γ indicates a frequency sweep speed, and ϕ(t) indicates a waveform of phase modulation to which the light propagating back and forth is subjected. An instantaneous frequency νinst. (t) of the electrical field waveform expressed by Relationship (<NUM>) is expressed by Equation (<NUM>) below. <NUM>] <MAT>.

In the instantaneous frequency of Equation (<NUM>), a third term on the right side corresponds to frequency modulation due to the dynamic strain. As illustrated in <FIG>, when the measurement time (Mes1, ···, MesN) of the beat signal is sufficiently shorter than a minimum period (the period Pmod in <FIG>) among frequency modulations (phase modulations) generated due to a plurality of dynamic strains applied to the sensing fiber, the phase modulation according to the dynamic strain can be approximated as linear phase modulation. Because the measurement time of the beat signal is a time section of each of "Mes1" and "Mes2" in <FIG>, when Mes1 and Mes2 << Pmod is satisfied, the phase modulation according to the dynamic strain is a frequency modulation having the amount of frequency modulation with a constant value. The amount of frequency modulation with a constant value corresponds to a slope of an arrow of the "instantaneous frequency νinst. " in <FIG>. Thus, when Equation (<NUM>) is further approximated, the instantaneous frequency νinst. (t) of the electrical field waveform is expressed as Equation (<NUM>). <NUM>] <MAT>.

In the instantaneous frequency of Equation (<NUM>), νoffset indicates a frequency offset that is given to the signal light as a linearly approximated phase modulation. In the OFDR, the beat signal between the local light and the signal light is measured, and the beat frequency of the beat signal Sbeat is associated with a distance from an incidence end of the sensing fiber to the strain section that is a target. A beat frequency fbeat between the signal light having the frequency offset given by Equation (<NUM>) and the local light is obtained by Equation (<NUM>) below. <NUM>] <MAT>.

In Equation (<NUM>), z is the distance, c is speed of light in the fiber, and 2z/c is a delay time between the local light and the signal light, which is the backscattered light from the second the dynamic strain section. Further, zoffset indicates a distance offset that is caused by the frequency offset νoffset.

<FIG> is a diagram illustrating an influence of the frequency modulation on the signal light in a time domain. In <FIG>, a state of a temporal change in the local light Llocal, the signal light Lsignal (backscattered light), and the beat frequency fbeat of the beat signal Sbeat between the local light and the signal light for two measurements of the optical frequency sweep (measurement periods Mes1 and Mes2) is shown. In each measurement, a delay of the local light Llocal and the signal light Lsignal is τ, and attention is paid only to backscattering from one dynamic strain section (strain sensor) which is a measurement target corresponding to a delay τ. Thus, it should be noted that the signal light input to the optical receiver is a signal in which backscattered light beams with different delay times from all the dynamic strain sections of the sensing fiber have been superimposed. The beat signal output from the optical receiver also includes different beat frequency components from all the dynamic strain sections.

As illustrated in <FIG>, the beat frequency of the beat signal Sbeat between the signal light Lsignal (backscattered light) and the local light Llocal in the dynamic strain section corresponding to the delay τ fluctuates as fbeat1 and fbeat2 for each measurement due to the frequency offset νoffset fluctuating as νoffset1 and νoffset2 caused by the dynamic strain. A difference between fbeat1 in the beat signal <NUM> and fbeat2 in the beat signal <NUM> means that the measurement distance fluctuates in the OFDR in which the beat frequency is associated with the distance from the incidence end of the fiber to the dynamic strain section.

<FIG> is a diagram illustrating an influence of the frequency modulation on the signal light as a spectrum change in a beat frequency domain. As illustrated in <FIG>, a distance offset zoffset assigned to each measurement changes because of the frequency offset νoffset due to the dynamic strain with respect to a peak position of the spectrum S<NUM> in a case in which there is no frequency offset. As a result, the distance to the dynamic strain section (strain sensor), which is a measurement target, is also measured with fluctuation. In dynamic strain distribution measurement using OFDR, a temporal change in the spectrum of the strain sensor due to mechanical vibration is measured at a freely selected distance on the sensing fiber. In Equation (<NUM>) above, it can be understood that the distance offset zoffset changes in each measurement due to mechanical vibration, which is the dynamic strain, and a distance to the strain sensor to be measured changes in each measurement. Measurement accuracy of the vibration distribution measurement device cannot be improved due to the fluctuation of the measurement distance.

In the related art (NPL <NUM>), for a change in the distance offset zoffset in each measurement in which the above-described frequency sweep is repeated, the amount of distance offset is computed in each measurement, and the strain sensor section is moved by the computed amount of distance offset so that the spectrum of the backscattered light is calculated. More specifically, in an electrical field E(τn) of the backscattered light obtained from the beat signal, a "window section for analyzing the optical frequency response" is shifted by the computed amount of distance offset. However, compensation arithmetic with the computation of the amount of distance offset or the shift calculation in the strain sensor section requires a high arithmetic power in the analysis unit of the vibration distribution measurement device. When the vibration is measured in a distribution manner at many measurement points of the entire sensing fiber or when the dynamic strain is measured over a long period of time, a load of signal processing in the vibration distribution measurement device is increased.

The inventors found that vibration can be analyzed with high accuracy using the vibration distribution measurement device illustrated in <FIG> without performing distance offset compensation arithmetic processing as in the related art by newly adding an additional condition in analysis of dynamic strain of OFDR.

Because OFDR uses the frequency swept light as probe light as disclosed in NPL <NUM>, it is possible to measure a spectrum S(v) of Rayleigh backscattered light at a freely selected position on the sensing fiber. <NUM>] <MAT>.

In the spectrum S(ν) of Relationship (<NUM>) above, v indicates an optical frequency, E(τn) indicates an electrical field of backscattered light from the delay τn (distance zn = cτn/<NUM>), and c indicates a speed of light in the optical fiber, N<NUM> indicates a position at which the optical frequency response is analyzed, and N<NUM> indicates a length for analysis of the optical frequency response.

As discussed in Relationship (<NUM>) and Equations (<NUM>) to (<NUM>) described above, the spectrum of Rayleigh backscattered light linearly shifts depending on the amount of strain in a longitudinal direction of the sensing fiber. In the vibration distribution measurement device of the embodiment, reference measurement for computing the spectrum shift is first performed under a condition that a period of vibration to which the sensing fiber is subjected is sufficiently longer than a period of the frequency sweep of the probe light, for measurement of the dynamic strain in a certain sensor section (not covered by the claimed invention). Next, a cross-correlation peak between a spectrum Sref(ν) obtained in this reference measurement and a spectrum Ssig(ν) obtained in the n-th measurement is computed so that the amount of spectrum shift (amount of strain) is analyzed.

However, in the OFDR in which the beat frequency is associated with the distance as described above, when the dynamic strain (vibration) is applied to the sensing fiber, a section in which the strain is analyzed changes through optical frequency modulation caused by the vibration. Thus, for the spectrum Sref(ν) obtained in the reference measurement and the spectrum Ssig(ν) obtained in the n-th measurement, a "deviation" corresponding to the distance offset occurs, and different fiber sections are assigned. <NUM>] <MAT>
[Math. <NUM>] <MAT>.

In Relationship (<NUM>) and Equation (<NUM>) above, Nd indicates a distance deviation amount (delay deviation amount) due to frequency modulation caused by the vibration, νoffset indicates the amount of frequency modulation caused by the vibration, γ indicates a frequency sweep speed of the frequency swept light, and T indicates the measurement time of the beat signal. The distance deviation amount Nd is obtained through signal processing for estimating a beat frequency offset from a cross-correlation between the waveform of the backscattered light in the reference measurement and the waveform of the backscattered light in each measurement, as in NPL <NUM>, for example.

<FIG> is a diagram illustrating the distance deviation of a measurement section resulting from frequency modulation due to vibration. <FIG> illustrates an electrical field waveform of the backscattered light after the beat signal is subjected to Fourier transform, and a horizontal axis corresponds to the beat frequency. As already described, in the OFDR, the beat frequency corresponds to the distance from the incidence end of the sensing fiber. An upper waveform in <FIG> shows the electrical field waveform of the backscattered light in the reference measurement. Sample points in a case in which there is no distance offset (distance deviation) in the window section for analyzing the optical frequency response corresponding to the length N<NUM> are shown. A lower waveform in <FIG> is the electrical field waveform of the backscattered light in the n-th measurement and illustrates sample points in a state in which the distance offset Nd occur. Further, when Fourier transform is performed on the window section of N<NUM>, for example, the spectra Sref(ν) and Ssig(ν) of the window section can be obtained.

As disclosed in NPL <NUM> or NPL <NUM>, the electrical field E(τn) of the backscattered light has a distribution in a longitudinal direction of the optical fiber in a freely selected section (a section from τ<NUM> to τN) that follows a Gaussian distribution. Thus, for the electrical field E(τn) of the backscattered light, when a distance deviation occurs due to frequency modulation caused by the vibration and fiber sections (sensor sections) assigned as the spectrum analysis section in the respective measurements completely differ, a correlation between the spectra disappears. Thus, there is no correlation between the Sref(ν) obtained in the reference measurement and the spectrum Ssig(ν) obtained in the n-th measurement, and the spectrum shift due to the dynamic strain cannot be computed.

However, when a spectrum analysis length N<NUM> of the electrical field E(τn) of the backscattered light is larger than the delay deviation Nd due to frequency modulation caused by the dynamic strain, a spectrum correlation is maintained by a common part illustrated in <FIG>. When a cross-correlation peak of two spectra of Sref(ν) and Ssig(ν) given in Relationship (<NUM>) is analyzed, the following equation is obtained. <NUM>] <MAT>.

In Equation (<NUM>), a peak signal-to-noise ratio (PSNR) indicates a noise ratio of the cross-correlation peak. In Equation (<NUM>), PSNR indicates that, when the spectrum analysis length N<NUM> is larger than the distance deviation Nd due to frequency modulation, this is not buried in the noise level and a cross-correlation has a peak. Thus, PSNR is an index indicating tolerance of a measurement method itself to a vibration, which is an observation target, in the vibration distribution measurement, and is one item indicating vibration measurement performance when the vibration distribution measurement using OFDR is performed.

<FIG> and <FIG> are diagrams illustrating a measurement example of a vibration distribution of the sensing fiber in which vibration is applied to two different positions through a comparison between the device of the related art and the device of the present disclosure. <FIG> is a diagram illustrating a configuration of a sensing fiber a vibration distribution of which has been analyzed. <FIG> is a diagram illustrating a comparison of vibration distribution analysis results between the vibration distribution measurement device of the related art and the vibration distribution measurement device of the present disclosure.

Referring to <FIG>, a total length of the sensing fiber that is a measurement target is <NUM>, and there is no vibration over <NUM> from an input end on the OFDR device side at a left end of <FIG> (between <NUM> to <NUM>). Subsequently, vibration at a frequency of <NUM> is applied in a section of a length of <NUM> (<NUM> to <NUM>). In the next section of a length of <NUM> (<NUM> to <NUM>), there is no vibration again, and subsequently, in a section of a length of <NUM> (<NUM> to <NUM>), vibration with a frequency of <NUM> is applied. Finally, there is no vibration in a section of a length of <NUM> (<NUM> to <NUM>). An end of the sensing fiber is an open end in an angled physical contact (APC) polished state.

Probe light of frequency swept light at a repetition frequency of <NUM> was repeatedly incident on the sensing fiber having the configuration in <FIG> <NUM> times from the vibration distribution measurement device illustrated in <FIG>, a beat signal was acquired from the backscattered light, and a distribution of the dynamic strain was analyzed. A starting light frequency of the probe light was <NUM> THz, and the frequency sweep speed was <NUM>/ms.

<FIG> illustrates results of analyzing the vibration distribution under the above-described measurement conditions in comparison with a vibration distribution measurement device of the related art in <FIG> and the vibration distribution measurement device of the present disclosure in <FIG>. In each case, a vertical axis indicates a length of the sensing fiber, which corresponds to a length of approximately <NUM>. A horizontal axis indicates time, and although color is not displayed in <FIG>, a spatial distribution and temporal fluctuation of an intensity of the strain are shown by shade of a region represented by a coordinate plane of the distance on the vertical axis and the time on the horizontal axis. In an original diagram of a graph of <FIG>, the amount of strain (µε) is shown by color at a right end, but in <FIG>, the amount is shown by only shade for reference. Further, time positions of a + peak and a - peak of the strain of the vibration at <NUM> corresponding to the distance of <NUM> to <NUM> in <FIG> are outlined.

In the case of the related art of (a), the length of the spectrum analysis section illustrated in <FIG> is <NUM>, and distance deviations Nd given by a vibration source of <NUM> and <NUM> under the above-described measurement conditions is <NUM> and <NUM>, respectively. Compared with the length of the spectrum analysis section of <NUM>, <NUM> of the distance deviation Nd given by the vibration source of <NUM> occupies <NUM>% or more, and it cannot be said that the length of the spectrum analysis section is sufficiently longer than the distance deviation Nd. In <FIG>, shade corresponding to the vibration at a frequency of <NUM> can be seen along the horizontal axis up to a distance of about <NUM>, but when the distance exceeds <NUM>, noise becomes predominant and a temporal fluctuation of the vibration at a frequency of <NUM> becomes unclear. A state in which a deviation occurs in an analysis section of the electrical field waveform of the backscattered light and PSNR in Equation (<NUM>) is degraded is shown. In the vibration distribution measurement device of the related art, a distance deviation amount is caused by frequency modulation due to the dynamic strain and the longer a distance to an observation point of the sensing fiber, the more the amount of frequency modulation is accumulated. The measurement is disabled from a position at which the accumulated amount of the frequency modulation, that is, the distance deviation Nd exceeds N<NUM>, as illustrated in <FIG>.

On the other hand, in the measurement of the vibration distribution measurement device of the present disclosure of <FIG>, the length of the spectrum analysis section is <NUM>, and a condition that the spectrum analysis section is sufficiently longer than the distance deviation Nd (<NUM> and <NUM>) given by the vibration source of <NUM> and <NUM> is satisfied. Two vibrations illustrated in <FIG> are clearly measured in an entire region of the sensing fiber having <NUM>. That is, shade corresponding to the vibration at a frequency of <NUM> (strain intensity increases or decreases <NUM> times in <NUM> seconds, which corresponds to <NUM>) can be observed in a range corresponding to <NUM> to <NUM> on a vertical axis. Further, shade corresponding to vibration at a frequency of <NUM> (strain intensity increases or decreases <NUM> times, which corresponds to <NUM>) can be observed in a range corresponding to <NUM> to <NUM> on the vertical axis. Thus, by making the spectrum analysis length N<NUM> of the electrical field E(τn) of the backscattered light larger than the delay deviation Nd due to the frequency modulation caused by the dynamic strain in the present embodiment, the dynamic strain can be measured without PSNR deterioration.

<FIG> is a flowchart illustrating a procedure for designating a reference waveform in a vibration distribution analysis method that can be performed by the vibration distribution measurement device of the present disclosure. This is a step of designating the reference waveform, which is a procedure before the measurement of the dynamic strain. The reference measurement is based on a spectrum in a state in which there is no dynamic strain in the sensing fiber. However, a waveform in a vibrated state can be used as the reference waveform as long as the waveform is a waveform measured with a signal measurement time sufficiently shorter than a minimum period of the dynamic strain applied to the sensing fiber (not covered by the claimed invention). That is, Mes1 and Mes2 << Pmod is established, as described with reference to <FIG>. In such a case, for example, a distribution waveform of an initial measurement result may be used as a reference waveform as it is.

Step (<NUM>-<NUM>): Measurement of the distribution waveform of the backscattered light is performed (Fourier transform of the beat signal). Here, the optical frequency response (beat signal) of the entire sensing fiber is measured.

Step (<NUM>-<NUM>): The spectrum analysis section is designated. Here, the spectrum analysis length N<NUM> of the electrical field E(τn) of the backscattered light is set. Here, the length N<NUM> of the spectrum analysis section of the electrical field E(τn) of the backscattered light is set to be longer than the delay deviation Nd due to the frequency modulation caused by the dynamic strain.

Step (<NUM>-<NUM>): Measurement of the backscattered light spectrum in the designated spectrum analysis section is performed. Here, the optical frequency response in the designated spectrum analysis section is subjected to Fourier transform so that the spectrum at each distance corresponding to the spectrum analysis section is analyzed. It is possible to extract a freely selected spectrum analysis section of the backscattered light waveform and obtain a spectrum of the backscattered light in the freely selected section of the sensing fiber.

A difference from the related art is that, in step (<NUM>-<NUM>), as a new analysis condition, the spectrum analysis length N<NUM> is limited with respect to the delay deviation Nd.

<FIG> is a flowchart illustrating an analysis procedure of a time waveform of the dynamic strain in a vibration distribution analysis method that can be performed by the vibration distribution measurement device of the present disclosure.

In this analysis procedure, first, the following steps (<NUM>-<NUM>) to (<NUM>-<NUM>) are repeated N times at a measurement time i = <NUM> to N-<NUM>. N in <FIG> is the number of repeated measurements of OFDR (see <FIG>).

(<NUM>-<NUM>) Measurement of the distribution waveform of the backscattered light at the measurement time i is performed (Fourier transform of the beat signal).

(<NUM>-<NUM>) The backscattered light spectrum at measurement time i is calculated.

After the above iterative step ends, a temporal change in the spectrum of the strain sensor (spectrum analysis section) is analyzed as step (<NUM>-<NUM>).

Next, the following steps (<NUM>-<NUM>) to (<NUM>-<NUM>) are repeated N times at i = <NUM> to N-<NUM>.

(<NUM>-<NUM>) A cross-correlation between the spectrum and the reference spectrum at the measurement time i is obtained.

(<NUM>-<NUM>) A spectrum shift at the measurement time i is computed. Because the optical spectrum in a freely selected section is spectrum shifted depending on the amount of strain, the amount of spectrum shift with respect to the optical spectrum obtained in the reference measurement is analyzed every N measurements of i = <NUM> to N-<NUM>.

Finally, as step (<NUM>-<NUM>), the time waveform of the dynamic strain is obtained in a range of time i = <NUM> to N-<NUM> from results of the spectrum shifts of steps (<NUM>-<NUM>) to (<NUM>-<NUM>) and the procedure is ended.

The Fourier transform arithmetic processing and the calculation processing of the spectrum and the amount of spectrum shift in the above-described steps of <FIG> and <FIG> can be performed after pieces of beat signal data at a certain number of times are accumulated by the acquired data storage <NUM> of the vibration distribution measurement device in <FIG> and a series of measurements ends. It should be noted that it is not always necessary to perform the above-described arithmetic processing in real time in synchronization with the application of actual frequency swept light to the sensing fiber.

Thus, a method of measuring a dynamic strain of a sensing fiber of the present disclosure in the device illustrated in <FIG> includes repeatedly supplying frequency swept lights to the sensing fiber, receiving backscattered light from the sensing fiber, combining local light, which is a part of a frequency swept light of the frequency swept lights, with the backscattered light to generate a beat signal, executing Fourier transform on the beat signal to obtain a backscattered light waveform, executing Fourier transform in a freely selected spectrum analysis section of the backscattered light waveform to obtain an optical spectrum, and obtaining, based on the amount of spectrum shift between a spectrum obtained in reference measurement and a spectrum obtained from the backscattered light for one sweep of the frequency swept lights, a strain amount of the dynamic strain in a section of the sensing fiber corresponding to the spectrum analysis section, in which a length of the spectrum analysis section is set to be larger than a distance deviation amount Nd of the spectrum analysis section caused by frequency modulation due to the dynamic strain.

The respective steps of <FIG> and <FIG> are executed by a program on a computer constituting the analysis unit <NUM> in <FIG>, but the analysis unit <NUM> in <FIG> and optionally the acquired data storage <NUM> can be disposed at places away from the components <NUM> to <NUM> in <FIG> and connected to a network.

The vibration distribution measurement system and the vibration distribution measurement device according to the embodiment of the present invention can be implemented by a computer and a program that execute the above analysis method, and the program can be recorded on a recording medium or provided through a network.

Although the processing for calculating the spectrum shift is included in the above-described embodiment, a more detailed condition for setting the length N of the spectrum analysis section with respect to the distance deviation amount Nd will be clarified in the next embodiment.

In the vibration distribution measurement device of the first embodiment described above, the dynamic strain in the sensing fiber can be performed with high accuracy under the condition that the length of the spectrum analysis section for the backscattered light is larger than the delay deviation Nd due to the frequency modulation caused by the dynamic strain. However, a reference waveform is obtained (<FIG>), the delay deviation Nd caused by a vibration phenomenon is estimated as the distance offset (steps <NUM>-<NUM> to <NUM>-<NUM>), the spectrum analysis section is shifted by Nd, and the dynamic strain is calculated (steps <NUM>-<NUM>). Calculations related to these distance offsets require processing time.

The inventors further examined a method capable of further simplifying the arithmetic processing in the vibration distribution measurement device and shortening the processing time. A relationship between the length of the spectrum analysis section and the delay deviation Nd, which was focused on in the vibration distribution measurement device of the first embodiment, was examined in further consideration of statistical properties of the Rayleigh backscattered light, and a concept of "tolerance to the distance offset" was further proposed. In the vibration distribution measurement device of the second embodiment to be described below, a concept of "tolerance" will be described in detail, and a more simplified vibration distribution measurement device and method configured based on the tolerance will be described. In the following description, the delay deviation (delay deviation) Nd is referred to as the distance offset Nd, but these are used as having the same meaning. It should also be noted that the spectrum shift on the optical frequency in the spectrum analysis of the first embodiment corresponds to the distance offset Nd converted into a distance on the sensing fiber.

In measurement of the backscattered light in OFDR, an influence of occurrence of the distance offset (delay deviation) in a measurement section that is a target due to the vibration of the measurement target is ascertained more quantitatively. Thus, the properties of Rayleigh backscattered light were statistically analyzed so that "tolerance" of the dynamic strain measurement to the distance offset Nd was obtained. A theoretical value of the "tolerance" was described by two parameters including the "length of spectrum analysis section" and the "distance offset Nd", as a probability that a vibration can be measured correctly even when there is the distance offset.

<FIG> is a flowchart for obtaining a theoretical value of tolerance of the dynamic strain measurement to the distance offset Nd. In the flowchart of <FIG>, an electrical field waveform of the Rayleigh backscattered light is generated by a different combination of "spectrum analysis section length N" and the "distance offset Nd" by simulation. In an optical spectrum corresponding to the generated electrical field waveform, a probability of a correlation peak value between the reference measurement and the trial measurement exceeding a correlation noise level was obtained. In the simulation, trial measurement in which a correlation peak was calculated for about <NUM> combinations as different combinations of N and Nd was repeated. With the trial measurement, a probability that the correlation peak can be detected was obtained with N and Nd as parameters.

<FIG> is a diagram schematically illustrating processing in each of the steps of simulation regarding the tolerance of the dynamic strain measurement. Hereinafter, a simulation method will be described with reference to the steps of <FIG> and the schematic diagrams of descriptions of processing operations of <FIG> alternately. As will be described below, <FIG> corresponds to step <NUM>-<NUM> in <FIG>, <FIG> corresponds to step <NUM>-<NUM> in <FIG>, and <FIG> corresponds to step <NUM>-<NUM> in <FIG>.

Referring back to <FIG>, simulation is started in step <NUM>-<NUM>, and N and Nd for one trial measurement are set in step <NUM>-<NUM>. In next step <NUM>-<NUM>, a reference waveform and a trial measurement waveform are generated. That is, an electrical field waveforms E(τ) for each of the reference measurement and the trial measurement are generated, as illustrated in <FIG>. Here, the electrical field waveform was generated as a random variable whose average value of an amplitude is <NUM> and follows a normal distribution. The distance offset Nd was set for a certain length N of the spectrum analysis section in each of the reference measurement and the trial measurement. The distance offset Nd was randomly set either before or after the spectrum analysis section for the reference measurement and the trial measurement. Thus, there are Nd non-common parts between the two generated waveforms.

Then, in step <NUM>-<NUM> of <FIG>, each of the electrical field waveforms in the reference measurement and the trial measurement is spectrally analyzed and the spectrum S(v) is obtained. That is, the spectrum of the reference measurement and the spectrum of the trial measurement are obtained, as illustrated in <FIG>. When no distance offset occurs, the two spectra mostly match, as illustrated in <FIG>. When the distance offset occurs, spectra in different sections of the sensing fiber are analyzed and spectrum shapes do not match.

In step <NUM>-<NUM> of <FIG>, a spectrum correlation is analyzed for each of the spectra in the reference measurement and the trial measurement, and the presence or absence of a correlation peak is determined. That is, an optical frequency shift Δν is obtained from the spectrum of the reference measurement and the spectrum of the trial measurement, as illustrated in <FIG>. A correlation value calculated here exceeds the noise level as illustrated in <FIG> when there is a common part between an electrical field waveform of the reference measurement and an electrical field waveform of the trial measurement, and a peak of the correction value appears. A case in which a correlation peak level exceeds the noise level is defined as a state in which vibration can be measured correctly. When the correlation peak level exceeds the noise level, the amount of strain can be measured correctly and thus, it is possible to obtain the amount of strain applied to the sensing fiber without shifting the spectrum analysis section according to the amount of distance offset. In one trial measurement, a determination is made as to whether the correlation peak level exceeds the noise level, the trial measurement is repeated (NTrial times) for one combination of N and Nd, and a determination is made as to whether the vibration measurement is correctly performed. Simulation of the probability distribution in which vibration can be measured correctly without the reference measurement in an N-Nd space can be performed through repetition of steps <NUM>-<NUM> to <NUM>-<NUM>.

Referring back to <FIG>, after the presence or absence of the correlation peak is determined in step <NUM>-<NUM>, the processing returns to step <NUM>-<NUM>, and steps <NUM>-<NUM> to <NUM>-<NUM> are repeated NTrial times. When repetition of the total number of trials NTrial ends, a probability distribution in which vibration can be correctly measured at a certain point in the N-Nd space without the reference measurement is calculated in step <NUM>-<NUM>. Specifically, the number of times NSuccess at which the correlation peak can be detected is divided by the number of trials NTrial so that the probability is computed. After the probability is calculated, Nd or N is changed in the next step <NUM>-<NUM> and the processing returns to step <NUM>-<NUM>. In step <NUM>-<NUM>, a combination of Nd and N is determined in a range and particle size that are targets, and the probability is obtained at all points in the N-Nd coordinate space. After all combinations of Nd and N are set and the probability is obtained, the simulation ends in step <NUM>-<NUM>.

As described with reference to <FIG> and <FIG>, it is possible to obtain the probability distribution in which vibration can be measured correctly without the reference measurement in the combination of N and Nd. When the distance offset Nd assumed in the measurement target can be known in advance before measurement, the length N of the spectrum analysis section is simply selected appropriately based on this probability distribution, making it possible to measure the vibration correctly without the reference measurement.

<FIG> is a diagram illustrating a probability distribution with which vibration (dynamic strain) can be correctly measured without the reference measurement using N and Nd as parameters. <FIG> is a diagram illustrating the probability distribution for the combination of N and Nd in shades. A horizontal axis indicates the length N of the spectrum analysis section, and a vertical axis indicates the distance offset Nd. It is difficult to see because the original color diagram is converted into monochrome shade, but a distribution is such that a probability value P gradually increases from <NUM> to <NUM> from an upper left corner to a lower right corner of a graph. <FIG> is a diagram illustrating a relationship between the distance offset Nd and a probability P with which vibration can be measured correctly without the reference measurement, with the length N of the spectrum analysis section as a parameter.

Referring to <FIG>, it can be seen that, when the length N of the spectrum analysis section is larger, the correlation peak can be detected with a probability P = <NUM> for a longer distance offset Nd. That is, when N is set to be larger, the tolerance to distance offset increases, and the probability with which the vibration (the dynamic strain) can be measured correctly without the reference measurement is increased. More specifically, referring to <FIG>, in the measurement of the dynamic strain that is a target, for example, when the assumed distance offset is Nd = <NUM>, the probability P that can be read from a curve of N = <NUM> is equal to or higher than <NUM> and is substantially close to <NUM>. <FIG> illustrates that, when Nd = <NUM> is assumed in the measurement of the target and N is set to <NUM> or more, the vibration can be correctly measured without the reference measurement with a probability of <NUM>% or more. When Nd having different values are assumed, the length N of the spectrum analysis section in which the vibration can be correctly measured with a predetermined probability value P without the reference measurement can be determined according to each figure of <FIG>.

A distribution of the probability value P in the N-Nd coordinate space illustrated in <FIG> is obtained by simulations described in <FIG> and <FIG>, but indicates tolerance to the distance offset Nd in the dynamic strain measurement, that is, the vibration distribution measurement. With only two parameters, N and Nd, it is possible to know limit conditions under which vibration can be measured correctly without the reference measurement, and limit conditions can be used for design of the tolerance to the distance offset Nd in the vibration distribution measurement. When the length N of the spectrum analysis section is set so that the reference measurement is unnecessary, this means that the computation of the distance offset becomes unnecessary. Next, when N is selected based on a theoretical probability value P in the N-Nd coordinate space illustrated in <FIG>, the validity as an index of tolerance to Nd of the probability distribution obtained by simulation is verified.

<FIG> is a diagram illustrating a configuration of a vibration system used for verification of the validation as a tolerance index of theoretical probability value in the N-Nd coordinate space. A total length from a signal input end of the OFDR device to the end of the sensing fiber, which is a measurement target, is <NUM> and a vibration at a frequency of <NUM> is applied in a section (<NUM> to <NUM>) of <NUM> from the input end. In the next section (<NUM> to <NUM>) of the length of <NUM>, there is no vibration. The end of the sensing fiber is an open end in an angled physical contact (APC) polished state. In the vibration system of <FIG>, the distance offset assumed in a vibration at a frequency of <NUM> in a first section is <NUM> (Nd = <NUM>). Nd = <NUM> is a dimensionless amount obtained by dividing the distance offset of <NUM> by a spatial resolution of OFDR of <NUM>.

<FIG> is a diagram illustrating results of vibration measurement when two different N values have been selected. The vibration measurement results when the lengths N of the spectrum analysis sections are set to <NUM> and <NUM>, respectively, for the vibration system having the configuration illustrated in <FIG> are shown. It should be noted that any of the measurements is made "without distance offset arithmetic processing". <FIG> show time (ms) on a horizontal axis and a length of the sensing fiber on a vertical axis, are diagrams that are the same as <FIG>, and show a vibrated state detected by the sensing fiber having a length of <NUM> in <FIG>. Although color is not displayed in <FIG>, a spatial distribution and temporal fluctuation of an intensity of the strain are shown by shade of a region represented by a coordinate plane of the time on the horizontal axis and the distance on the vertical axis. In a graph of <FIG>, the amount of strain (µε) is shown by color at a right end, but the amount is shown by only shade for reference. Further, time positions of a + peak and a - peak of the strain of the vibration at <NUM> (a period is <NUM>) in a range of a distance of <NUM> to <NUM> in <FIG> are outlined.

<FIG> shows a vibrated state measured when the length N of the spectrum analysis section is set to <NUM>. The vibration at <NUM> is observed in a section of distance <NUM> to <NUM>, but a peak position fluctuates in a time axis direction, noise is observed, and the vibration is unclear. Here, referring to a simulation theoretical value in <FIG>, a probability P of positions of Nd = <NUM> and N = <NUM> is about <NUM>, and corresponds to a state in which the probability with which vibration can be measured correctly without the reference measurement is <NUM>.

On the other hand, <FIG> shows a vibrated state measured when the length N of the spectrum analysis section is set to <NUM>. A vibration at <NUM> in a section of a distance of <NUM> to <NUM> is shown very clearly, and there are no fluctuation and noise in the time axis direction. Referring back to <FIG>, a probability P of positions of Nd = <NUM> and N = <NUM> being outside a graph is substantially <NUM>, and the probability with which vibration can be measured correctly without the reference measurement is <NUM>, which corresponds to a state in which the vibration measurement can be performed <NUM>% correctly. It can be understood from <FIG> that two measurement results in which N has been selected based on a theoretical probability distribution in the N-Nd coordinate space illustrated in <FIG> correspond to a probability value P of a peak of a correlation value exceeding a noise level, and theoretical probability distribution in the N-Nd coordinate space shows the tolerance to the distance offset Nd.

<FIG> show a probability of being analyzed as an "vibration-free state" as a measurement result over a section in which the length of the sensing fiber in the vibration-free state is <NUM> to <NUM>. A horizontal axis indicates time (ms), a vertical axis indicates the probability of being analyzed as a "vibration-free state", and the probability of being analyzed as no vibration should be originally <NUM> in a section of <NUM> to <NUM>.

Here, referring to <FIG> in which the length N of the spectrum analysis section is set to <NUM>, the probability of being analyzed as a "vibration-free state" periodically decreases to nearly <NUM>. This corresponds to the fact that various levels of strain (µε) are erroneously observed periodically and noise appears in the corresponding section of <NUM> to <NUM> of <FIG>. In <FIG>, it can be seen that a measured value and a theoretical value are shown together, and when N = <NUM> is set, false vibration is observed even in the "vibration-free state".

On the other hand, referring to <FIG> in which the length N of the spectrum analysis section is set to <NUM>, the probability of being analyzed as a "vibration-free state" does not fall below <NUM>, and it can be said that the probability is substantially always close to <NUM>. That is, a section in which the length of the sensing fiber in a vibration-free state is <NUM> to <NUM> is correctly analyzed as the "vibration-free state" with a probability of substantially <NUM>%. This corresponds to the fact that a strain level is displayed as <NUM> in a state in which there is substantially no noise in a corresponding section of <NUM> to <NUM> in <FIG>.

It can be understood from <FIG> that the two measurement results in which N has been selected based on theoretical probability distribution in the N-Nd coordinate space illustrated in <FIG> correspond to the probability value P of the peak of the correlation value exceeding the noise level, and theoretical probability distribution in the N-Nd coordinate space shows the tolerance to the distance offset Nd. Thus, when the length N of the spectrum analysis section is selected depending on theoretical probability distribution in the N-Nd coordinate space obtained in <FIG>, the dynamic strain can be measured correctly even when arithmetic processing of the offset distance is not performed without the reference measurement. The length N of the spectrum analysis section in this case is selected so that a probability of the correlation peak exceeding the noise level has a desired value or more close to <NUM>, when the correlation measurement is performed on the spectrum of the reference measurement with respect to the value of the assumed distance offset Nd.

It was found that the vibration can be measured correctly without the reference measurement and arithmetic processing of the distance offset or the like by using the probability distribution of the correlation value peak level exceeding the noise level as an index of tolerance using two parameters N and Nd shown in <FIG> and selecting an appropriate N. N at which the probability of the correlation peak exceeding the noise level when the correlation measurement is performed on the spectrum of the reference measurement becomes about <NUM> is selected with respect to the assumed distance offset Nd based on the theoretical probability in the N-Nd coordinate space, so that vibration measurement with high tolerance to the distance offset Nd can be achieved.

A value of the above-described probability P can be set to a desired value close to <NUM> depending on accuracy required for the vibration measurement device or a vibration phenomenon of the measurement target. For example, in an environment in which high accuracy is required for a measured value, the length N of the spectrum analysis section is determined with respect to the assumed distance offset Nd, with the probability value P having a desired value being <NUM>. Further, in an environment in which high accuracy is not required for a measured value, N can be determined from a wider range with respect to an assumed Nd, the probability value P having a desired value being <NUM>.

As described above, the length N of the spectrum analysis section in which there is tolerance to the distance offset Nd is set based on theoretical probability in the N-Nd coordinate space, making it possible to greatly simplify the configurations of the vibration distribution measurement device and method of the present embodiment as compared to the first embodiment. In the flowchart of <FIG>, a value in which there is the tolerance to the assumed distance offset Nd is selected for the length N of the spectrum analysis section.

Thus, the strain measurement in the vibration distribution measurement device and method of the second embodiment can be performed by the following steps.

Step III above corresponds to step <NUM>-<NUM> in <FIG>, step IV corresponds to step <NUM>-<NUM> in <FIG>, step V corresponds to repetition of steps <NUM>-<NUM> and <NUM>-<NUM> in <FIG>, and step VI corresponds to repetition of steps <NUM>-<NUM>, <NUM>-<NUM> and <NUM>-<NUM>, and step <NUM>-<NUM> in <FIG>. Further, steps I to VI above are substantially the same as N repeated static strain measurements.

Thus, a method of measuring a dynamic strain of the present disclosure includes repeatedly supplying frequency swept lights to a sensing fiber, receiving backscattered light from the sensing fiber, combining local light, which is a part of a frequency swept light of the frequency swept lights, with the backscattered light to generate a beat signal, obtaining a backscattered light waveform from the beat signal, obtaining an optical spectrum in a freely selected spectrum analysis section of the backscattered light waveform, and obtaining a strain amount of the dynamic strain in a section of the sensing fiber corresponding to the spectrum analysis section, in which a length N of the spectrum analysis section is set with respect to an assumed distance deviation amount Nd caused by the dynamic strain such that a probability of a correlation peak exceeding a noise level has a desired value when correlation measurement between a spectrum of reference measurement and the optical spectrum is performed, the desired value being approximate to <NUM>.

Further, a configuration of the vibration distribution measurement device is not changed from the configuration of <FIG>, and steps I to VI above is performed as processing operations different from the first embodiment in the analysis unit <NUM>. Thus, a device that measures the dynamic strain of the present disclosure includes an optical circuit unit that repeatedly supplies frequency swept lights to a sensing fiber, receive backscattered light from the sensing fiber, and combine local light, which is a part of a frequency swept light of the frequency swept lights, with the backscattered light, and an optical reception and analysis unit that generates a beat signal from the local light and the backscattered light, obtain a backscattered light waveform from the beat signal, obtain an optical spectrum in a freely selected spectrum analysis section of the backscattered light waveform, and obtain a strain amount of the dynamic strain in a section of the sensing fiber corresponding to the spectrum analysis section, in which a length N of the spectrum analysis section is set with respect to an assumed distance deviation amount Nd caused by the dynamic strain such that a probability of a correlation peak exceeding a noise level has a desired value when correlation measurement between a spectrum of reference measurement and the optical spectrum is performed, the desired value being approximate to <NUM>.

As described above, the series of steps obtained by simplifying the flowcharts of <FIG> and <FIG> of the present embodiment does not include arithmetic operation processing of the relevant distance offset. Thus, this is nothing but processing for measuring static strain in the OFDR. This becomes the same state as when an amount of static strain applied to the sensing fiber, which does not fluctuate with time, is measured by OFDR instead of the dynamic strain such as a vibration phenomenon. Thus, it becomes possible to measure "the dynamic strain" which is a fluctuation of temporal strain by simply repeating measurement of "static strain" periodically. When the dynamic strain is steady, once the maximum value of the distance offset Nd is actually measured and N is determined to be a preparatory step described in step <NUM> described above, a distance offset processing step can be omitted and thus the measurement time can be greatly shortened. That is, a device and a measuring method that accurately measure vibration at a designated position without using digital signal processing to compensate for distance fluctuation become possible.

The probability distribution of the correlation value peak level exceeding the noise level, in which the two parameters N and Nd shown in <FIG> have been used, shows an example in which N is set to <NUM> or more when Nd = <NUM> is assumed in the measurement of the target. In the N-Nd coordinate space, this can be specified as a region (boundary) in which the probability value P becomes about <NUM>. That is, the region is a region in which P(N, Nd) = <NUM> in the probability value P given in <FIG>.

As described in detail above, the vibration distribution measurement device of the present disclosure provides a device and a measuring method that accurately measure vibration at a designated position without using digital signal processing to compensate for distance fluctuation.

Claim 1:
A device (<NUM>) configured to measure dynamic strain of a sensing fiber (<NUM>), the device comprising:
an optical circuit unit (<NUM>, <NUM>, <NUM>, <NUM>) configured to repeatedly supply frequency swept lights to the sensing fiber,
receive backscattered light from the sensing fiber, and
combine local light, the local light being a part of a frequency swept light of the frequency swept lights, with the backscattered light; and
an optical reception and analysis unit (<NUM>, <NUM>, <NUM>) configured to generate a beat signal from the local light and the backscattered light,
execute Fourier transform on the beat signal to obtain a backscattered light waveform,
execute Fourier transform in a freely selected spectrum analysis section of the backscattered light waveform to obtain an optical spectrum, and
obtain, based on the amount of spectrum shift between a spectrum obtained in reference measurement and a spectrum obtained from the backscattered light for one sweep of the frequency swept lights, a strain amount of the dynamic strain in a section of the sensing fiber corresponding to the spectrum analysis section,
characterized in that
a length of the spectrum analysis section is set to be larger than a distance deviation amount Nd of the spectrum analysis section caused by frequency modulation due to the dynamic strain, and
the reference measurement is based on a spectrum in a state where there is no dynamic strain in the sensing fiber.