WO2020188888A1

WO2020188888A1 - Vibration detection optical fiber sensor and vibration detection method

Info

Publication number: WO2020188888A1
Application number: PCT/JP2019/045789
Authority: WO
Inventors: 祥宏神田
Original assignee: 沖電気工業株式会社
Priority date: 2019-03-18
Filing date: 2019-11-22
Publication date: 2020-09-24
Also published as: JP2020153704A; JP7238507B2

Abstract

In the present invention, a light intensity acquisition unit acquires the intensity of output light resulting from the backscattering of probe light in an optical fiber. An OTDR waveform acquisition means sequentially acquires and stores OTDR waveforms that are the waveforms of the output light for each optical pulse sent from the light intensity acquisition unit. A moving average acquisition means acquires and stores a kth moving average waveform by averaging kth to (k + M - 1)th temporally continuous OTDR waveforms. A differential waveform acquisition means acquires a differential waveform p'(x) for a (k + M + N)th OTDR waveform and the kth moving average waveform. A position-specific statistical information acquisition means acquires a sample average μ(x) and a standard deviation σ(x) for each position of kth to (k + i - 1)th (where i is an integer greater than or equal to 2) differential waveforms. A degree-of-abnormality acquisition means calculates the degrees of abnormality at each position. An abnormality determination means determines whether there has been abnormal vibration on the basis of the values of the degrees of abnormality, and when there has been abnormal vibration, acquires the position at which the abnormal vibration occurred.

Description

Vibration detection Optical fiber sensor and vibration detection method

The present disclosure is a vibration detection optical fiber sensor and vibration that can be used to detect illegal intrusion of people in large facilities such as power plants and factories, and to detect cracks in large civil engineering structures such as bridges and roads. Regarding the detection method.

When a light wave is input to an optical fiber, backscattered light is generated as it propagates. The backscattered light generated at each position in the longitudinal direction of the optical fiber is observed with a delay of the time required for the input light wave (hereinafter referred to as input light) to reciprocate from the input end of the optical fiber to the position where the backscattered light is generated. Will be done. When there is a breaking point or a point that greatly attenuates the propagating light in the longitudinal direction of the optical fiber, the light intensity of the backscattered light corresponding to the position changes. This principle is used to detect a break point of a communication optical fiber, and is known as a time domain reflection measurement (OTDR: Optical Time Domain Reflectometry).

In OTDR, vibration transmitted to an optical fiber can be detected simply by replacing the laser used to generate an optical pulse as input light with a laser having a small frequency drift and a narrow line width (for example, J.C. Juarez, See W. Maier, K. N. Choi, and H. F. Taylor, “Distributed Fiber-Optic Intrusion Sensor System,” IEEE JLT, vol. 23, No. 6, June 2005, pp. 2081-2087).

In OTDR using such a laser, the waveform of the light intensity of the backscattered light is the strong interference of the coherent backscattered light generated at multiple scattering centers in the optical fiber while the optical pulse propagates through the optical fiber. Observed as a result.

When vibration is applied to the optical fiber, the refractive index and birefringence of the optical fiber change only at that position. As a result, the relative phase difference between the plurality of coherent backscattered lights from the plurality of scattering centers at the position where the vibration is applied changes. In the superposition of coherent waves, the waveform changes when the relative phase difference between the waveforms changes. Therefore, in the observed OTDR waveform, only the waveform at the time corresponding to the position where the vibration is applied changes. By utilizing this phenomenon, the occurrence and position of vibration can be detected by calculating the difference between the observed OTDR waveform and the OTDR waveform observed at a previous time.

This method is called phase sensitive OTDR or φ-OTDR. In the phase sensitive OTDR, the difference from the OTDR waveform at the previous time is obtained. At this time, even if there is no abnormal vibration, if the amplitude fluctuation of the OTDR waveform is large, a large noise will be generated. This causes the presence or absence of abnormal vibration to be overlooked or false detection for the purpose of detecting abnormal vibration. In order to reduce this noise, methods called moving average (moving averaging) and moving differential (moving differential) are known (for example, Y. Lu, T. Zhu, L. Chen, and X. Bao, “ Distributed Vibration Sensor Based on Coherent Detection of Phase-OTDR, ”IEEE JLT, vol. 28, No. 22, Nov. 15, 2010, pp.3243-3249).

The moving average is a method of reducing noise when the difference is calculated by obtaining the moving average for each position of a plurality of OTDR waveforms observed with the passage of time. The movement difference is the difference from the OTDR waveform observed at a slightly distant time, instead of calculating the difference from the OTDR waveforms adjacent to each other on the time axis when calculating the difference between the OTDR waveforms observed at different times. Is a method of calculating. When the duration of the vibration is sufficiently long compared to the time interval for acquiring the OTDR waveform, a large difference value can be obtained as compared with the method of calculating the difference between the OTDR waveforms adjacent to each other on the time axis.

However, even if the moving average and the moving difference are used for the OTDR waveform, the ratio of the values with and without abnormal vibration is about several times. As long as the presence or absence of vibration is evaluated only by the difference, the magnitude of the difference when there is abnormal vibration is, in principle, limited to the difference between the maximum value and the minimum value of the possible amplitude of the OTDR waveform at the maximum. As described above, since the ratio of the values when there is abnormal vibration and when there is no abnormal vibration is small, the operation of the vibration detection optical fiber sensor is unstable. Further, while the phase-sensitive OTDR can detect even a slight vibration, the degree of noise depends on the environment in which the optical fiber is laid. Therefore, it is difficult to distinguish the presence or absence of an abnormality in an environment affected by wind or the like.

The present disclosure statistically evaluates changes in the OTDR waveform by utilizing the nature of noise in addition to the use of moving averages and moving differences in OTDRs, such as phase-sensitive OTDRs, to determine the presence or absence of abnormal vibrations. Provided are a vibration detection optical fiber sensor and a vibration detection method for clearly discriminating.

The first aspect of the present disclosure is a vibration detection optical fiber sensor, which comprises a light source unit, an optical fiber, a light intensity acquisition unit, an OTDR waveform acquisition unit, a moving average acquisition unit, a difference waveform acquisition unit, and the like. It includes a position-specific statistical information acquisition unit, an abnormality degree acquisition unit, and an abnormality determination unit.

The light source unit generates an optical pulse as probe light. The probe light is input to the optical fiber, and the light intensity acquisition unit acquires the intensity of the output light in which the probe light is backscattered by the optical fiber.

The OTDR waveform acquisition unit sequentially acquires and stores the OTDR waveform, which is the waveform of the output light for each optical pulse sent from the light intensity acquisition unit. The moving average acquisition unit acquires the kth moving average waveform by averaging the temporally continuous OTDR waveforms of the kth (k is an integer of 1 or more) to the k + M-1 (M is an integer of 2 or more). And remember. The difference waveform acquisition unit acquires the moving average waveform of the k + N-1 (N is an integer of 2 or more) and the difference waveform p'(x) of the kth moving average waveform. The position-specific statistical information acquisition unit acquires the sample mean μ (x) and the standard deviation σ (x) of each position for the difference waveforms of the kth to k + i-1 (i is an integer of 2 or more). The abnormality degree acquisition unit acquires the abnormality degree waveform indicating the abnormality degree by calculating the abnormality degree a (p') of each position using the following equation (1).

The abnormality determination unit acquires the presence or absence of abnormal vibration and the position of the abnormal vibration from the value of the degree of abnormality.

The second aspect of the present disclosure is an exemplary embodiment in the above aspect, and the light source unit may be configured to include a narrow line width laser, a function generator, and an intensity modulator. The narrow line width laser is a light source that generates a laser beam, and the line width may be 10 kHz or less. Even if the function generator generates an electrical pulse at a constant frequency. The intensity modulator may generate an optical pulse by converting the laser beam into an optical pulse with an electric pulse.

A third aspect of the present disclosure is a vibration detection method, which includes an optical pulse generation process, a light intensity acquisition process, an OTDR waveform acquisition process, a moving average acquisition process, a difference waveform acquisition process, and a position-based method. It includes a statistical information acquisition process, an abnormality degree acquisition process, and an abnormality determination process.

In the optical pulse generation process, an optical pulse is generated as probe light. In the light intensity acquisition process, the intensity of the output light in which the probe light is backscattered by the optical fiber is acquired. In the OTDR waveform acquisition process, the OTDR waveform, which is the waveform of the output light for each optical pulse acquired in the light intensity acquisition process, is sequentially acquired and stored. In the moving average acquisition process, the kth moving average waveform is acquired by averaging the temporally continuous kth (k is an integer of 1 or more) to k + M-1 (M is an integer of 2 or more) OTDR waveforms. And remember. In the difference waveform acquisition process, the moving average waveform of the k + N-1 (N is an integer of 2 or more) and the difference waveform p'(x) of the kth moving average waveform are acquired. In the position-specific statistical information acquisition process, the sample mean μ (x) and standard deviation σ (x) of each position are acquired for the difference waveforms of the kth to k + i-1 (i is an integer of 2 or more). In the abnormality degree acquisition process, an abnormality degree waveform indicating the abnormality degree is acquired by calculating the abnormality degree a (p') at each position using the above equation (1). In the abnormality determination process, the presence or absence of abnormal vibration and the position of the abnormal vibration, if any, are acquired from the value of the degree of abnormality.

A fourth aspect of the present disclosure is an exemplary embodiment in the third aspect, wherein the optical pulse generation process includes a process of generating a laser beam using a narrow line width laser having a line width of 10 kHz or less and a constant frequency. A process of generating an electric pulse and a process of generating an optical pulse by converting a laser beam into an electric pulse with an electric pulse may be provided.

According to the above aspect, the vibration detection optical fiber sensor and the vibration detection method of the present disclosure not only take the difference waveform between the OTDR waveform and the moving average waveform, but also obtain the sample average and standard deviation of each position of the difference waveform. The degree of abnormality is obtained by using, and the presence or absence of abnormal vibration is determined from the value of this degree of abnormality. As a result, the ratio of the peak values when there is an abnormality and when there is no abnormality becomes larger than the ratio of the peak values in the difference waveform, and the presence or absence of the abnormality becomes clearer. Furthermore, it is possible to clearly judge the detection of distributed abnormal vibrations that suppress the dependence of sensitivity on the distance from the input end of the optical fiber.

It is a block diagram which shows the schematic structure of the vibration detection optical fiber sensor. It is a schematic diagram which shows the concept of determination of anomalous vibration. It is a figure which shows the result of the characteristic test.

Hereinafter, exemplary embodiments of the present disclosure will be described with reference to the drawings, but the shape, size, and arrangement of each component are only schematically shown to the extent that the present disclosure can be understood. .. Further, although a preferable configuration example of the present disclosure will be described below, numerical conditions and the like are merely suitable examples. Therefore, the present disclosure is not limited to the following exemplary embodiments, and many modifications or modifications can be made to achieve the effects of the present disclosure without departing from the scope of the configuration of the present disclosure.

(Vibration detection optical fiber sensor)
An exemplary embodiment of the vibration detection optical fiber sensor will be described with reference to FIG. FIG. 1 is a block diagram showing a schematic configuration of a vibration detection optical fiber sensor.

This vibration detection optical fiber sensor is configured to include a light source unit 10, an optical fiber 20, an optical circulator 30, and a measurement unit 40. This vibration detection optical fiber sensor is used for OTDR.

The light source unit 10 periodically generates an optical pulse as probe light. The spatial resolution of the vibration detection optical fiber sensor depends on the width of this optical pulse. Further, the measurement distance of the vibration detection optical fiber sensor depends on the frequency of the optical pulse. The optical pulse takes 5 ns to propagate 1 m through the optical fiber 20. When observing backscattered light, it takes a round trip time of forward propagation and reverse propagation, so a delay of 10 ns per 1 m occurs. For example, when the pulse width is 100 ns and the frequency is 5 kHz, the spatial resolution is 10 m and the maximum measurement distance is 20 km.

The light source unit 10 includes, for example, a laser light source 12, an intensity modulator 14, a function generator 16, and an optical amplifier 18.

The laser light source 12 generates laser light as continuous light in the communication wavelength band. As the laser light source 12, it is preferable to use a so-called narrow line width laser having a line width of 10 kHz or less. When a narrow line width laser is used as the laser light source 12, this vibration detection optical fiber sensor can be used for a phase sensitive OTDR. The wavelength of the laser beam may be arbitrary, but it is preferable to use a standard single-mode optical fiber with a low loss of 1550 nm. The laser beam generated by the laser light source 12 is sent to the intensity modulator 14.

The function generator 16 generates a rectangular electric pulse. This electric pulse is sent to the intensity modulator 14. The electric pulse generated by the function generator 16 has, for example, a pulse width of 100 nsec and a repetition frequency of 5 kHz. The output of the function generator 16 is also sent to the analog-to-digital (A / D) converter 44, which will be described later, and is used as a trigger signal.

The intensity modulator 14 converts the laser beam into an optical pulse with an electric pulse to generate an optical pulse. This optical pulse is sent to the optical amplifier 18. The pulse width and frequency of the optical pulse generated by the intensity modulator 14 are both the same as the electric pulse generated by the function generator 16. In this example, the optical pulse has a pulse width of 100 nsec and a repetition frequency of 5 kHz.

The optical pulse generated by the intensity modulator 14 is amplified by the optical amplifier 18 and then sent to the optical fiber 20 as probe light via the optical circulator 30.

The probe light sent to the optical fiber 20 propagates in the optical fiber 20, and backscattered light is generated along with the propagation of the probe light. This backscattered light is sent to the measuring unit 40 as output light via the optical circulator 30.

The measurement unit 40 includes a photodetector 42 and an A / D converter 44 as a light intensity acquisition unit, and a calculation unit 50. The output light input from the optical fiber 20 to the measuring unit 40 is sent to the photodetector 42.

The photodetector 42 can be composed of, for example, a photodiode (PD). The photodetector 42 converts the output light into an electric signal by square detection and sends it to the A / D converter 44. Here, the light intensity of the backscattered light is small. Therefore, it is desirable that the photodetector 42 has a sensitivity to receive light even at about -30 dBm.

The A / D converter 44 converts the electric signal received from the photodetector 42 into a digital signal. Here, the sampling frequency of the A / D converter 44 may be large enough to sample the OTDR waveform, and about 200 MHz is sufficient.

The digital signal obtained by the A / D converter 44 is sent to the arithmetic unit 50.

The arithmetic unit 50 statistically evaluates changes in the OTDR waveform using a digital signal, detects the presence or absence of abnormal vibration, and identifies the position of vibration. As the arithmetic unit 50, for example, a commercially available personal computer (PC) in which a program for detecting the presence or absence of abnormal vibration and specifying the position of vibration is installed can be used. Here, as an example, the arithmetic unit 50 will be described as being configured to include a CPU (Central Processing Unit) 60, a RAM (Random Access Memory) 52, a ROM (Read Only Memory) 54, and a storage unit 56. The CPU 60 realizes each functional unit described later by executing a program stored in the ROM 54. The processing result of each functional unit is temporarily stored in the RAM 52.

The OTDR waveform is periodically sent from the A / D converter 44 to the arithmetic unit 50. The OTDR waveform is represented by, for example, taking x indicating the distance from the input end of the optical fiber 20 on the horizontal axis and taking the signal strength p (x, k) at the distance x on the vertical axis. Here, k is an integer of 1 or more and corresponds to the number of the optical pulse input to the optical fiber 20. The number k of this optical pulse indicates the elapsed time. The end of the optical fiber 20 on the side connected to the optical circulator 30 is used as the input end.

The OTDR waveform acquisition unit 62 sequentially stores, that is, stores the OTDR waveforms of the kth to k + M-1, which are OTDR waveforms for M pieces (M is an integer of 2 or more) that are continuous in time, in the storage unit 56. After M OTDR waveforms are stored, the oldest OTDR waveform is deleted and a new OTDR waveform is stored in the storage unit 56.

Each time the moving average acquisition unit 64 receives the OTDR waveform from the A / D converter 44, the moving average of the kth moving average is obtained by averaging the M OTDR waveforms of the kth to k + M-1 that are continuous in time. Get the waveform. The kth moving average waveform is stored in the storage unit 56. The time-consecutive N moving average waveforms (N is an integer of 2 or more) acquired by the moving average acquisition unit 64 are sequentially stored in the storage unit 56. After N moving average waveforms are stored, the oldest moving average waveform is deleted and a new moving average waveform is stored in the storage unit 56.

The difference waveform acquisition unit 66 obtains the kth difference waveform p'(x) by taking the difference between the k + N-1 moving average waveform and the kth moving average waveform stored in the storage unit 56. get. In addition, this difference may be referred to as a moving difference. The kth moving average waveform is stored in the storage unit 56. The time-continuous i moving average waveforms (i is an integer of 2 or more) acquired by the moving average acquisition unit 64 are sequentially stored in the storage unit 56. After i moving average waveforms are stored, the oldest moving average waveform is deleted and a new moving average waveform is stored in the storage unit 56.

The position-specific statistical information acquisition unit 68 acquires the sample average μ (x) and the standard deviation σ (x) of each position for the difference waveforms of the kth to k + i-1.

Here, the sample average μ (x) and the standard deviation σ (x) can be obtained, for example, at intervals of about the spatial resolution of this vibration detection optical fiber sensor. Alternatively, it may be the interval between measurement points adjacent to each other in the longitudinal direction of the optical fiber 20 when viewed from the digital data acquired by the A / D converter 44.

The abnormality degree acquisition unit 70 acquires a waveform indicating the abnormality degree (hereinafter, also referred to as an abnormality degree waveform) by calculating the abnormality degree at each position using the above equation (1).

After that, the abnormality determination unit 72 determines the presence or absence of abnormal vibration from the degree of abnormality. The threshold value used to determine the presence or absence of the degree of abnormality may be set to an appropriate value by monitoring the influence of the environment of the installation location.

The number M of waveforms when calculating the moving average is determined, for example, based on the line width of the laser light source 12. When the line width of the laser light source 12 is 500 Hz, it is preferable that the average time is 2 msec (= 1/500 Hz) or less. When the frequency of the optical pulse is 5 kHz, for example, M can be set to 20. In this case, the moving average waveform becomes the average for a time of 4 msec.

When the position-specific statistical information acquisition unit 68 acquires the sample average μ (x) and the standard deviation σ (x), it is preferable to perform the calculation without abnormal vibration. Therefore, for example, when the abnormality determination unit 72 determines that there is an abnormality vibration, the position-specific statistical information acquisition unit 68 excludes the sample corresponding to the position, and the sample average μ (x) and the standard deviation. It can be configured to acquire σ (x). Further, when the abnormality determination unit 72 determines that there is an abnormal vibration, the moving average acquisition unit 64 may be configured to acquire the moving average waveform excluding the corresponding OTDR waveform.

The time difference between the OTDR waveform and the moving average waveform when acquiring the difference waveform should be set to be equal to or longer than the period of abnormal vibration. For example, when the frequency of the optical pulse is 5 kHz, the period of the optical pulse is 0.2 msec (= 1/5 kHz). On the other hand, the frequency of vibration generated by hitting the optical fiber is about several hundred Hz. In this case, since the period is about several msec, the difference between temporally adjacent OTDR waveforms is small. Therefore, it is preferable to set the time difference to about 4 msec.

The number i of the difference waveforms when the position-specific statistical information acquisition unit 68 acquires the sample mean μ (x) and the standard deviation σ (x) of each position is the number M of the waveforms when calculating the moving average. Similarly, for example, it is determined based on the line width of the laser light source 12. When the line width of the laser light source 12 is 500 Hz, it is preferable that the average time is 2 msec (= 1/500 Hz) or less. When the frequency of the optical pulse is 5 kHz, i can be set to 20, for example.

Here, with reference to FIG. 2, the concept of determining abnormal vibration used in the present disclosure will be described. FIG. 2 is a schematic diagram for explaining the concept of determining abnormal vibration. FIG. 2 (1) shows the difference waveform, and FIG. 2 (2) shows the degree of abnormality obtained based on the difference waveform shown in FIG. 2 (1). In (1) of FIG. 2, as an example, i = 4 and four difference waveforms from the first to the fourth are shown.

The three axes (1) in FIG. 2 show the distance from the input end of the optical fiber 20, the elapsed time, and the light intensity of the difference waveform. Further, in FIG. 2 (2), the horizontal axis shows the distance from the input end of the optical fiber, and the vertical axis shows the degree of abnormality.

As shown in (1) of FIG. 2, when a disturbance different depending on the position is applied to the optical fiber 20, the value of the difference waveform is obtained in the section where the disturbance is applied even if there is no abnormal vibration. Becomes larger.

On the other hand, in the vibration detection optical fiber sensor of the present disclosure, the position-specific statistical information acquisition unit 68 acquires the sample average μ (x) and standard deviation σ (x) of each position, and the abnormality determination unit 72 causes an abnormality. Judge the presence or absence of abnormal vibration from the degree. By normalizing with the sample mean μ (x) and standard deviation σ (x) at each position in this way, as shown in FIG. 2 (2), the degree of anomaly in the section where the disturbance is always applied is small. And its impact is underestimated. Therefore, the position dependence of the disturbance can be suppressed.

(Characteristic test)
The characteristic test of this vibration detection optical fiber sensor will be described. Here, the length of the optical fiber 20 is set to 18.2 km. The optical fiber at a position of 15.6 to 15.8 km from the input end side of the optical fiber 20 was vibrated and abnormal vibration was applied.

The laser light source 12 was a narrow line width laser having a line width of 500 Hz. The pulse width and frequency of the optical pulse were set to 100 nsec and 5 kHz, respectively. Further, the number M of the OTDR waveforms used when calculating the moving average was set to 20, and the time difference when acquiring the difference waveform was set to 4 msec. This time difference of 4 msec is equivalent to 20 light pulses when converted to the number N of optical pulses. Further, the number i of the difference waveforms when acquiring the sample mean μ (x) and the standard deviation σ (x) was set to 20.

That is, in this test, the number M of waveforms used when calculating the moving average, the number N of optical pulses converted from the time difference when acquiring the difference waveform, the sample mean μ (x) and the standard deviation σ ( The number i of the difference waveforms when acquiring x) is equal to 20. The storage unit stores M (= 20) OTDR waveforms, N (= 20) moving average waveforms, and i (= 20) moving difference waveforms.

(1) and (2) of FIG. 3 are diagrams showing the results of the characteristic test. (1) of FIG. 3 is a diagram showing a movement difference, in which the horizontal axis shows the distance [unit: km] from the input end of the optical fiber 20, and the vertical axis shows the signal level [unit: mV]. It is shown. FIG. 3 (2) is a diagram showing the degree of abnormality, in which the horizontal axis shows the distance [unit: km] from the input end of the optical fiber 20 and the vertical axis shows the degree of abnormality.

As shown in (2) of FIG. 3, in the range of 0 to 15.6 km and the range of 15.8 to 18.2 km, the signal level of the movement difference is the same as that no abnormal vibration is applied. The signal level varies with distance.

On the other hand, as shown in FIG. 3 (2), in the vibration detection optical fiber sensor of the present disclosure, the degree of abnormality is calculated using the above equation (1), and as a result, the range is 0 to 15.6 km. In the range of 15.8 to 18.2 km, the degree of anomaly is almost 0, and it can be clearly distinguished from the degree of anomaly in the section of 15.6 to 15.8 km to which the abnormal vibration is applied.

In this way, this vibration detection optical fiber sensor statistically evaluates the fluctuation of the OTDR waveform at each position in the vicinity of each time with the passage of time. Thereby, the detection of the distributed abnormal vibration that suppresses the dependence of the sensitivity on the distance from the input end of the optical fiber 20 can be clearly determined. Furthermore, even if the external environment changes slowly, the effect of suppressing the influence on the sensitivity of vibration detection can be expected.

(Other exemplary embodiments)
Here, an example of measuring the intensity of backscattered light by using a PD as a light intensity acquisition unit has been described, but the present invention is not limited to this. Y. Lu, T. Zhu, L. Chen, and X. Bao, “Distributed Vibration Sensor Based on Coherent Detection of Phase-OTDR,” IEEE JLT, vol. 28, No. 22, Coherent detection may be performed by interfering the backscattered light with the laser light generated by the laser light source, as described in Nov. 15, 2010, pp.3243-3249. Alternatively, a 90-degree light hybrid receiver may be used to observe the phase change of the backscattered light. Heterodyne detection may also be performed on the backscattered light.

The entire disclosure of Japanese Patent Application No. 2019-050121 filed on March 18, 2019 is incorporated herein by reference in its entirety.

All documents, patent applications, and technical standards described herein are to the same extent as if the individual documents, patent applications, and technical standards were specifically and individually stated to be incorporated by reference. Incorporated herein by reference.

Claims

A light source that generates an optical pulse as probe light,
The optical fiber to which the probe light is input and
A light intensity acquisition unit that acquires the intensity of the output light in which the probe light is backscattered by the optical fiber.
An OTDR waveform acquisition unit that sequentially acquires and stores an OTDR waveform that is a waveform of output light for each optical pulse sent from the light intensity acquisition unit.
A moving average waveform of the kth is acquired and stored by averaging the OTDR waveforms of the kth (k is an integer of 1 or more) to the k + M-1 (M is an integer of 2 or more) that are continuous in time. Average acquisition department and
A moving average waveform of the k + N-1 (N is an integer of 2 or more), a difference waveform acquisition unit for acquiring the difference waveform p'(x) of the kth moving average waveform, and
For the difference waveforms of the kth to k + i-1 (i is an integer of 2 or more), a position-specific statistical information acquisition unit that acquires the sample mean μ (x) and standard deviation σ (x) at each position, and
An abnormality degree acquisition unit that acquires an abnormality degree waveform indicating the abnormality degree by calculating the abnormality degree a (p') at each position using the following equation (1).
From the value of the degree of abnormality, the presence or absence of abnormal vibration, and if there is abnormal vibration, an abnormality determination unit that acquires the position, and
A vibration detection fiber optic sensor.
When the position-specific statistical information acquisition unit determines that there is abnormal vibration in the abnormality determination unit, the sample mean μ (x) and standard deviation σ (x) are acquired except for the sample corresponding to the position. The vibration detection optical fiber sensor according to claim 1.
The vibration detection optical fiber sensor according to claim 1, wherein the moving average acquisition unit acquires a moving average waveform excluding the corresponding OTDR waveform when the abnormality determination unit determines that there is an abnormal vibration.
The light source unit
A light source that generates laser light, a narrow line width laser with a line width of 10 kHz or less, and
A function generator that generates electrical pulses at a constant frequency,
An intensity modulator that generates the optical pulse by converting the laser beam into an optical pulse with the electric pulse.
The vibration detection optical fiber sensor according to any one of claims 1 to 3.
An optical pulse generation process that generates an optical pulse as probe light,
The light intensity acquisition process for acquiring the intensity of the output light in which the probe light is backscattered by the optical fiber, and
An OTDR waveform acquisition process in which OTDR waveforms, which are waveforms of output light for each optical pulse acquired in the light intensity acquisition process, are sequentially acquired and stored, and
A moving average that acquires and stores the kth moving average waveform by averaging the temporally continuous kth (k is an integer of 1 or more) to k + M-1 (M is an integer of 2 or more) OTDR waveforms. Acquisition process and
The difference waveform acquisition process for acquiring the moving average waveform of the k + N-1 (N is an integer of 2 or more) and the difference waveform p'(x) of the kth moving average waveform.
For the difference waveforms of the kth to k + i-1 (i is an integer of 2 or more), the position-specific statistical information acquisition process for acquiring the sample mean μ (x) and standard deviation σ (x) at each position, and
Anomalous degree acquisition process for acquiring anomalous degree waveform indicating anomaly degree by calculating anomaly degree a (p') at each position using the following equation (1), and
From the value of the degree of abnormality, the presence or absence of abnormal vibration, and the abnormality determination process of acquiring the position if there is abnormal vibration,
A vibration detection method.
In the position-specific statistical information acquisition process, when it is determined that there is abnormal vibration in the abnormality determination process, the sample mean μ (x) and standard deviation σ (x) are acquired except for the sample corresponding to the position. The vibration detection method according to claim 5.
The vibration detection method according to claim 5, wherein in the moving average acquisition process, when it is determined that there is an abnormal vibration in the abnormality determination process, the moving average waveform is acquired by excluding the corresponding OTDR waveform.
The light pulse generation process is
The process of generating laser light using a narrow line width laser with a line width of 10 kHz or less,
The process of generating electrical pulses at a constant frequency,
The process of generating the optical pulse by converting the laser beam into an optical pulse with the electric pulse, and
The vibration detection method according to any one of claims 5 to 7, further comprising.