WO2014064771A1 - 異常検知システム及び異常検知方法 - Google Patents

異常検知システム及び異常検知方法 Download PDF

Info

Publication number
WO2014064771A1
WO2014064771A1 PCT/JP2012/077354 JP2012077354W WO2014064771A1 WO 2014064771 A1 WO2014064771 A1 WO 2014064771A1 JP 2012077354 W JP2012077354 W JP 2012077354W WO 2014064771 A1 WO2014064771 A1 WO 2014064771A1
Authority
WO
WIPO (PCT)
Prior art keywords
intensity distribution
optical fiber
abnormality detection
fir filter
light
Prior art date
Application number
PCT/JP2012/077354
Other languages
English (en)
French (fr)
Inventor
宇野　和史
武井　文雄
丈夫 笠嶋
恭子 只木
石鍋　稔
Original Assignee
富士通株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 富士通株式会社 filed Critical 富士通株式会社
Priority to EP12887297.5A priority Critical patent/EP2913647B1/en
Priority to PCT/JP2012/077354 priority patent/WO2014064771A1/ja
Priority to JP2014543051A priority patent/JP5880728B2/ja
Publication of WO2014064771A1 publication Critical patent/WO2014064771A1/ja
Priority to US14/694,440 priority patent/US9347803B2/en
Priority to US15/149,508 priority patent/US9528860B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
    • G01D5/35354—Sensor working in reflection
    • G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection

Definitions

  • the present invention relates to an abnormality detection system and an abnormality detection method.
  • an abnormality detection system including a temperature distribution measuring device (Distributed Temperature Sensor: DTS) that uses an optical fiber as a temperature sensor may be employed.
  • DTS distributed Temperature Sensor
  • an optical fiber is laid around a pipe or tank, and the end of the optical fiber is connected to a temperature distribution measuring device. Then, laser light is irradiated into the optical fiber from the temperature distribution measuring device, and the Raman scattered light generated in the optical fiber is detected by the temperature distribution measuring device to obtain the temperature of the pipe, tank, etc. Determine if there is an abnormality.
  • An object of the present invention is to provide an abnormality detection system and an abnormality detection method capable of detecting abnormalities occurring in facilities such as chemical plants, refineries and thermal power plants at an early stage.
  • an optical fiber is connected to one end side and the other end side of the optical fiber, and light is incident on the optical fiber from the one end side, and a first intensity distribution of backscattered light
  • the backscattered light detection unit that enters the optical fiber from the other end side to acquire the second intensity distribution of the backscattered light, and the first acquired by the backscattered light detection unit
  • an abnormality detection system including a data processing unit that calculates for each position in the length direction of the fiber and determines the presence or absence of abnormality from the result.
  • light is incident on one end side of the optical fiber by a backscattered light detection unit to obtain a first intensity distribution of backscattered light
  • the backscattered light detection unit A step of acquiring light on the other end side of the optical fiber to obtain a second intensity distribution of backscattered light, and a first FIR (Finite Impulse Response) filter on the first intensity distribution by the data processing unit.
  • a first FIR Finite Impulse Response
  • an anomaly occurring in a facility such as a chemical plant, a refinery factory, or a thermal power plant can be detected at an initial stage.
  • FIG. 1 is a view showing a state in which an optical fiber is wound with a constant tension around a portion where a branch pipe is welded to a main pipe.
  • FIG. 2 is a diagram showing the results of examining transmission loss when a moderate bending is applied to an optical fiber, when a slightly strong bending is applied, and when a strong bending is applied.
  • FIG. 3 is a diagram specifically illustrating moderate bending, slightly strong bending, and strong bending.
  • FIG. 4 is a diagram showing a result of differentiating the intensity distribution of the return light shown in FIG. 2 twice.
  • FIG. 5 is a diagram showing the intensity distribution of Raman scattered light when bending stress is applied to the optical fiber, and the result obtained by differentiating the intensity distribution twice.
  • FIG. 1 is a view showing a state in which an optical fiber is wound with a constant tension around a portion where a branch pipe is welded to a main pipe.
  • FIG. 2 is a diagram showing the results of examining transmission loss when a moderate
  • FIG. 6 is a view showing both the intensity distribution of the return light and the temperature distribution obtained by the temperature distribution measuring apparatus.
  • FIG. 7 is a block diagram illustrating the anomaly detection system according to the embodiment.
  • FIG. 8A is a view showing a state where one end side of the optical fiber and the beam splitter are optically connected
  • FIG. 8B is an optical view between the other end side of the optical fiber and the beam splitter. It is a figure which shows the state connected.
  • FIG. 9 is a flowchart (part 1) illustrating the abnormality detection method.
  • FIG. 10 is a flowchart (part 2) illustrating the abnormality detection method.
  • FIG. 11 is a diagram showing a first intensity distribution (NTS1) and a second intensity distribution (NTS2).
  • FIGS. 12A and 12B are diagrams illustrating an example of the FIR filter.
  • FIG. 13 shows the result of convolution of the FIR filter shown in FIG. 12A with the first intensity distribution (NTS1) and the FIR filter shown in FIG. 12B with the second intensity distribution (NTS2).
  • FIG. 14 is a graph obtained by calculating the product of the FIRNTS1 value and the FIRNTS2 value for each position in the length direction of the optical fiber.
  • FIGS. 15A and 15B are enlarged views of a part of FIG.
  • FIG. 16 is a diagram for explaining a second method for calculating the transmission loss.
  • FIG. 17 is a diagram illustrating a display example of transmission loss.
  • FIG. 18A is a diagram showing a result obtained by differentiating the first intensity distribution NTS1 of FIG. 11 twice
  • FIG. 18B is a diagram showing a part of FIG.
  • FIG. 19 is a diagram showing the results of obtaining the peak height (peak value) from FIGS. 14 and 15.
  • FIG. 20 is a diagram showing the relationship between the peak height and the average loss.
  • FIG. 21 is a diagram showing the result of obtaining the average loss using the calibration function F ( ⁇ P).
  • FIG. 22 is a diagram showing the positions of the half widths of the peaks P4 to P7 and the positions of LF and LR.
  • FIG. 23 is a diagram showing linear approximation of the first intensity distribution NTS1 and the second intensity distribution NTS2.
  • FIG. 24 is a diagram illustrating a result of obtaining transmission loss (linear approximate loss) at the positions of peaks P5 and P6.
  • FIG. 25 is a diagram showing availability 1.
  • FIG. 26 is a diagram showing availability 2.
  • FIG. 27 is a diagram illustrating availability 3.
  • an anomaly is detected by utilizing the fact that the transmission loss of the optical fiber changes due to the application of bending stress.
  • FIG. 1 is a view showing a state in which an optical fiber 13 is wound with a constant tension around a portion where a branch pipe 12 is welded to a main pipe 11.
  • the optical fiber 13 is partially fixed to the branch wiring 12 with a tape 14.
  • the flow of the solution or gas in the main pipe 11 and the branch pipe 12 changes, and the temperature of the main pipe 11 and the branch pipe 12 changes. Due to this temperature change, the main pipe 11 and the branch pipe 12 expand or contract, and the bending stress or tensile stress applied to the optical fiber 13 changes.
  • the optical fiber 13 has a large transmission loss when a bending stress or a tensile stress exceeding a certain level is applied. Therefore, for example, it is possible to determine whether there is an abnormality by comparing the transmission loss at the time of past operation or stop and the current transmission loss.
  • the horizontal axis is the distance (position in the length direction of the optical fiber), the vertical axis is the intensity of the light, and the optical fiber is moderately bent. It is a figure which shows the result of having investigated the transmission loss when applying a strong bending.
  • the moderate bending is the bending shown in FIG. 3A (bending radius is about 10 mm), and the slightly stronger bending is slightly stronger than FIG. 3A (see FIG. 3B).
  • the strong bend is a slightly stronger bend than that shown in FIG. 3B (see FIG. 3C).
  • the intensity of the return light is normalized with reference to the light quantity at a position in the length direction of the optical fiber of 0 m.
  • FIG. 2 shows that transmission loss corresponding to the degree of bending occurs at a position of about 340 m in the length direction of the optical fiber.
  • the optical fiber is moderately bent during normal operation and a certain amount of transmission loss occurs at a specific position in the length direction of the optical fiber. In this case, if the transmission loss of the optical fiber changes rapidly, it can be determined that some abnormality has occurred.
  • Patent Document 1 describes that the intensity distribution of return light is differentiated twice for the purpose of measuring the position and connection loss of an optical fiber with high accuracy. It is conceivable to use this method for detecting the presence or absence of abnormality.
  • FIG. 4 is a diagram showing the result of differentiating the return light intensity distribution shown in FIG. 2 twice.
  • a one-dot chain line in FIG. 4 indicates a range of 3 ⁇ ( ⁇ is a standard deviation).
  • the threshold value in order to remove noise components, it is necessary to set the threshold value to about 3 ⁇ .
  • the threshold value is set to 3 ⁇ , the peak level of the middle bending portion is slightly higher than the noise level, and it cannot be said that the detection reliability is high.
  • the optical pulse detector (Optical Time Domain Reflectometer: OTDR) used in Patent Document 1 uses Rayleigh scattered light, but Raman scattered light (Stokes light and anti-Stokes light) used in a temperature distribution measuring device (DTS). Similar results can be obtained using. If the temperature distribution measuring device is used, the temperature distribution measurement and the abnormality detection can be performed simultaneously.
  • the horizontal axis indicates the distance (position in the length direction of the optical fiber), the vertical axis indicates the amount of light (right side) and the threshold (left side), and the Raman scattered light when bending stress is applied. It is a figure which shows collectively intensity distribution and the result of having differentiated the intensity distribution twice.
  • NTS indicates a measurement result using Stokes light
  • TTS indicates a measurement result using anti-Stokes light.
  • the broken line in FIG. 5 indicates the range of 3 ⁇ .
  • the abnormality cannot be detected until the transmission loss becomes large to some extent, that is, the abnormality cannot be detected in the initial stage. .
  • the delay in detecting abnormalities leads to the expansion of accidents, so it is desirable to detect abnormalities at an early stage.
  • ⁇ Temperature is another reason why abnormalities cannot be detected in the initial stage. Since the intensity of the Raman scattered light changes with temperature, the intensity of the return light is related to the stress applied to the optical fiber and the temperature.
  • the horizontal axis represents distance (position in the length direction of the optical fiber), and the vertical axis represents temperature (left side) and return light intensity (right side). It is the figure which showed collectively the temperature distribution obtained by this.
  • the intensity of the return light is normalized based on the light amount at a position of 0 m in the length direction of the optical fiber.
  • bending stress is applied at a position of about 345 m in the length direction of the optical fiber.
  • FIG. 6 shows that when the transmission radius is increased by reducing the bending radius of the optical fiber, the temperature measurement error increases.
  • FIG. 6 there is no problem because the temperature difference before and after the position where the bending stress is applied is due to the difference in the optical fiber laying path.
  • the stress applied to the optical fiber laid in the rack changes for some reason, thereby changing the intensity of the return light, resulting in an accurate temperature. May not be detected.
  • the temperature is detected higher than the actual temperature and a case where the temperature is detected lower than the actual temperature.
  • an abnormality detection system that can detect an abnormality occurring in a facility such as a chemical plant, a refinery, or a thermal power plant at an early stage will be described.
  • FIG. 7 is a block diagram illustrating the anomaly detection system according to the embodiment.
  • the anomaly detection system includes a loop type light detection device 20 and a data processing device 30 that processes data output from the light detection device 20.
  • the loop-type photodetector 20 is an example of a backscattered light detector
  • the data processor 30 is an example of a data processor.
  • the loop type light detection device 20 includes a laser light source 21, a beam splitter 22, a transmission path switch 23, a light detection circuit unit 24, and a calculation unit 25, and is connected to an optical fiber 26 for use. Both ends of the optical fiber 26 are connected to the transmission path switch 23, and are laid around the pipes 11 and 12 as shown in FIG. 1 and partially fixed to the pipes 11 and 12 with a tape 14 or the like.
  • Laser light with a predetermined pulse width is output from the laser light source 21 at a constant cycle.
  • the laser light passes through the beam splitter 22 and enters the optical fiber 26 via the transmission path switch 23.
  • the transmission path switch 23 switches the laser light transmission path at a constant cycle. That is, the transmission path switch 23 is in a state where one end side of the optical fiber 26 and the beam splitter 22 are optically connected (see FIG. 8A), and the other end side of the optical fiber 26 and the beam splitter 22. Are alternately switched to the state in which they are optically connected (see FIG. 8B).
  • a part of the light that has entered the optical fiber 26 is backscattered by molecules constituting the optical fiber 26.
  • the backscattered light returns through the optical fiber 26 and reaches the beam splitter 22 through the transmission path switch 23. Then, the light is reflected by the beam splitter 22 and reaches the light detection circuit unit 24.
  • the photodetection circuit unit 24 is provided with a filter (not shown) that separates light of a predetermined wavelength and a light receiving element (not shown) that receives light of a predetermined wavelength separated by the filter. A signal corresponding to the intensity of the light received by the light receiving element is output from the light detection circuit unit 24.
  • the calculation unit 25 includes a computer.
  • the arithmetic unit 25 stores changes with time of the signal output from the light detection circuit unit 24 and outputs the data to the data processing device 30.
  • the data processing device 30 includes a computer.
  • the data processing device 30 processes the data output from the light detection device 20 to determine whether there is an abnormality as will be described later, and executes a preset process such as generating an alarm when it is determined as abnormal. To do.
  • An optical pulse detector that uses Rayleigh scattered light may be used as the loop-type photodetecting device 20, and a temperature distribution measuring device that uses Raman scattered light (Stokes light and anti-Stokes light) ( DTS) may be used.
  • OTDR optical pulse detector
  • DTS temperature distribution measuring device that uses Raman scattered light
  • a temperature distribution measuring device is used as the light detection device 20, the temperature distribution can be measured together with the detection of the abnormality.
  • the inventors of the present application have proposed a temperature measurement method in which the temperature measurement value at other measurement points is corrected using a transfer function on the basis of the temperature at a specific measurement point (Patent Document 3, etc.). According to this method, the temperature of the measurement points set at intervals of 10 cm to several tens of cm along the length direction of the optical fiber can be detected with high accuracy.
  • step S11 an initialization flag “1” is set at a predetermined position in the length direction of the optical fiber 26.
  • step S11 an initialization flag “1” is set at a position where a transmission loss is predicted as in the connection unit. Also, when a certain amount of stress is applied in advance to a specific position of the optical fiber 26, the initialization flag “1” is set at that position.
  • step S12 the data processing device 30 acquires the intensity distribution in the length direction of the optical fiber obtained when the one end side of the optical fiber 26 is irradiated with the laser light (see FIG. 8A). , The first intensity distribution. Further, the data processing device 30 acquires the intensity distribution in the length direction of the optical fiber obtained when the other end side of the optical fiber 26 is irradiated with the laser light (see FIG. 8B), and the second Intensity distribution.
  • the horizontal axis indicates the distance (position in the length direction of the optical fiber 26), and the vertical axis indicates the intensity (light quantity) of the return light (backscattered light) detected by the light detection device 20.
  • NTS1 1 intensity distribution
  • NS2 2nd intensity distribution
  • the horizontal axis in FIG. 11 indicates a predetermined position on the laser light source 21 side as a reference position in a state where the optical detection device 20 and one end side of the optical fiber 26 are optically connected as shown in FIG. (0m) is the distance from the reference position. Further, the transmission loss at the positions indicated by A1 and A2 in FIG. 11 is caused by a connector that optically connects between the light detection device 20 and one end side and the other end side of the optical fiber 26.
  • the intensity distribution (NTS1) of the backscattered light when the laser light is irradiated to one end of the optical fiber 26 has a higher intensity of the backscattered light closer to the reference position (0 m position).
  • the intensity of backscattered light decreases as the distance from the reference position increases.
  • the intensity distribution (NTS2) of the backscattered light when the other end side of the optical fiber 26 is irradiated with the laser light has a lower intensity of the backscattered light closer to the reference position. The farther away from the position, the higher the intensity of the backscattered light.
  • step S13 the process proceeds to step S13, and the data processing device 30 applies a differential FIR (Finite Impulse Response) filter to the first intensity distribution (NTS1) and the second intensity distribution (NTS2). Make it work.
  • FIR Finite Impulse Response
  • a differential FIR filter is a filter having both characteristics of a differential and a low-pass filter, unlike a unit step type differential filter.
  • An example of the FIR filter used in the present embodiment is shown in FIGS.
  • FIG. 12A shows an FIR filter used when acting on the first intensity distribution (NTS1) from the reference position side
  • FIG. 12B shows the second intensity distribution (NTS2) from the reference position side. It is an FIR filter used when acting.
  • the FIR filter of FIG. 12B has a positive / negative symmetrical characteristic with respect to the FIR filter of FIG.
  • the FIR filter of Fig.12 (a) can be used.
  • FIG. 13 convolves the FIR filter shown in FIG. 12A with the first intensity distribution (NTS1), and converts the FIR filter shown in FIG. 12B into the second intensity distribution (NTS2). It is a figure which shows the result of convolution.
  • the result of convolving the FIR filter shown in FIG. 12A with the first intensity distribution (NTS1) is shown as FIRNTS1, and the FIR filter shown in FIG. 12B is changed into the second intensity distribution (NTS2).
  • the result of convolution is shown as FIRNTS2.
  • the transmission loss is emphasized by convolving the first FIR filter and the second FIR filter with the first intensity distribution and the second intensity distribution, respectively.
  • FIRNTS1 since one end side of the optical fiber 26 is set as the reference position, FIRNTS1 has high sensitivity on the reference position side, and FIRNTS2 has high sensitivity on the side opposite to the reference position.
  • step S14 the data processing device 30 calculates the product of the FIRNTS1 value and the FIRNTS2 value for each position L in the length direction of the optical fiber 26.
  • FIG. 14 is a graph in which the product of the FIRNTS1 value and the FIRNTS2 value is calculated for each position in the length direction of the optical fiber 26 and graphed.
  • 15A is an enlarged view of the range of 50 to 100 m in FIG. 14, and
  • FIG. 15B is an enlarged view of the range of 1800 to 2100 m in FIG.
  • step S15 the data processing device 30 compares the value of the product of the FIRNTS1 value and the FIRNTS2 value with the set value for each position L.
  • the product of the FIRNTS1 value and the FIRNTS2 value is larger than the set value, it is determined that there is an abnormality, and an abnormality flag “1” is set at the position L.
  • the set value may be set to about 3 ⁇ , for example. However, when a default loss described later is registered at the position L, the set value is (default loss + error range).
  • step S16 the data processing device 30 determines whether or not there is a position where the abnormality flag “1” is set. If it is determined that there is a position for which the abnormality flag “1” is set, the process proceeds to step S17. If it is determined that there is no position for which the abnormality flag “1” is set, the process returns to step S12 to continue the processing. .
  • step S17 the data processing device 30 determines whether or not the initialization flag “1” is set at the position where the abnormality flag “1” is set. When it is determined that the initialization flag “1” is set at the position where the abnormality flag “1” is set, the process proceeds to step S18, and when it is determined that the initialization flag “1” is not set. The process proceeds to step S21.
  • the data processing apparatus 30 registers the FIRNTS1 value and the FIRNTS2 value at the position L where the initialization flag “1” is set as default loss, and stores the default loss in the data processing apparatus 30.
  • the FIRNTS1 value at the position L is described as FIRNTS1 (L)
  • the FIRNTS2 value at the position L is described as FIRNTS2 (L).
  • step S19 the data processing device 30 resets the initialization flag “1”. And it returns to step S12 and repeats the process mentioned above.
  • step S15 since the initialization flag has already been reset, if there is a part where the abnormality flag “1” is set in step S15, the process proceeds from step S17 to step S21.
  • step S21 the data processing device 30 generates an alarm to notify the abnormality, for example. Thereafter, the process proceeds to step S22, and the data processing device 30 performs the loss amount quantification calculation by the first method.
  • the data processing device 30 calculates the loss amount by, for example, the method described below.
  • the data processing device 30 calculates the abnormality detection signal P (L) by the following equation (1) for each position.
  • the abnormality detection signal P (L) FIRNTS1 (L) ⁇ FIRNTS2 (L) (1)
  • the abnormality detection signal P (L) shows a peak waveform.
  • the effective peak height ⁇ P is The following equation (2) is obtained.
  • ⁇ P P (Lalert) ⁇ Pave (2)
  • the transmission loss Loss1 (L) is calculated by the following equation (3) with the calibration function being F ( ⁇ P).
  • step S23 the data processing device 30 attempts to quantify the loss amount by the second method.
  • Another loss point may be a position of an adjacent peak, or a position where a default loss is registered such as a connection portion of an optical fiber. Then, the front and rear of the position where the abnormality is detected are linearly approximated by the least square median method (LMedS) or the like, and the loss (dB) is obtained from the difference between the intercepts.
  • LMedS least square median method
  • the abnormality detection signal P (L) is calculated by the following equation (4) for each position.
  • the abnormality detection signal P (L) FIRNTS1 (L) ⁇ FIRNTS2 (L) (4)
  • the abnormality detection signal P (L) shows a peak waveform.
  • the effective peak height ⁇ P is The following equation (5) is obtained.
  • the positions LF and LR at both ends of the half-value width of the peak of interest are LF N and LR N, and the half-value width of the peak on the reference position side is more than that.
  • the two positions at both ends are LF N-1 and LR N-1 .
  • Two positions at both ends of the half-value width of the peak on the opposite side are denoted as LF N + 1 and LR N + 1 .
  • first intensity distribution (NTS1) and the second intensity distribution (NTS2) are linearly approximated between LR N-1 and LF N.
  • first intensity distribution (NTS1) and the second intensity distribution (NTS2) are linearly approximated between LR N and LF N + 1 .
  • ⁇ PNTS1 (Lalert) 2 ⁇ abs (PNTS11 (Lalert) ⁇ PNTS12 (Lalert)) / (PNTS11 (Lalert) + PNTS12 (Lalert)) (6)
  • ⁇ PNTS2 (Lalert) 2 ⁇ abs (PNTS21 (Lalert) ⁇ PNTS22 (Lalert)) / (PNTS21 (Lalert) + PNTS22 (Lalert)) (7)
  • the transmission loss Loss2 (L) is calculated by the following equation (8).
  • LF N and LR N may be set at positions slightly away from the peak center rather than the half-width position.
  • the amount of light in the case where a loss has occurred exhibits a kink-like characteristic. Therefore, if the area showing the kink-like characteristic is removed, the residual with respect to the approximate straight line data becomes smaller.
  • step S24 the data processing device 30 determines whether or not the transmission loss Loss2 (L) is obtained by the second method.
  • the process proceeds to step S25.
  • step S25 the transmission loss Loss1 (L) obtained by the first method and the transmission loss Loss2 (L) obtained by the second method are compared, and the accuracy A (%) is obtained by the following equation (9).
  • A 10 B ⁇ 100 (9)
  • B ( ⁇ abs (Loss1 (L) ⁇ Loss2 (L)) ⁇ 10).
  • the unit of transmission loss Loss1 (L) and transmission loss Loss2 (L) is dB.
  • the data processing device 30 displays, for example, the position where the abnormality has occurred, the average value of the transmission loss Loss1 (L) and the transmission loss Loss2 (L), and the accuracy A on the display.
  • FIG. 17 shows a display example of transmission loss.
  • “loss by peak” is the transmission loss Loss1 (L)
  • “loss by linear approximation” is the transmission loss Loss2 (L).
  • step S24 the process proceeds from step S24 to step S26. Then, the data processing device 30 displays, for example, the position where the abnormality has occurred and the transmission loss Loss1 (L) on the display.
  • the abnormality detection system can accurately detect a minute change in stress applied to an optical fiber. Thereby, the occurrence of an abnormality in a facility such as a chemical plant, a refinery factory, or a thermal power plant can be detected at an early stage, and the occurrence of an accident or the expansion of the accident can be prevented.
  • FIG. 13 is a diagram obtained by convolving the FIR filters of FIGS. 12A and 12B with the first intensity distribution and the second intensity distribution shown in FIG.
  • the FIR filter of FIGS. 12A and 12B is convoluted with the intensity distribution of FIG. 11, the difference in the amount of light before and after the portion where the stress is applied is emphasized.
  • FIG. 14 shows the result of calculating the product of the FIRNTS1 value and the FIRNTS2 value at each position in the length direction of the optical fiber.
  • the product of the FIRNTS1 value and the FIRNTS2 value is calculated at each position in the length direction of the optical fiber. By calculating, the difference in the amount of light before and after the portion where the stress is applied is further emphasized.
  • FIG. 18 (a) is a diagram showing a result obtained by differentiating the first intensity distribution NTS1 of FIG. 11 twice
  • FIG. 18 (b) is a diagram showing an enlarged part of FIG. 18 (a).
  • FIG. 19 shows the peak height (peak value) obtained from FIGS.
  • the average loss in FIG. 19 is obtained by calculating a difference in light quantity visually from an intensity distribution diagram and converting it to a transmission loss.
  • the peak of P4 has no significant difference with respect to the noise level, and it is difficult to detect this peak.
  • the peaks of P1 and P7 can be reliably detected. That is, it can be said that the method according to the present embodiment has detection sensitivity at least twice as high as that of the method of differentiating the intensity distribution of the return light twice.
  • FIG. 20 is a diagram showing the relationship between the peak height and the average loss based on FIG. From FIG. 20, the calibration function F ( ⁇ P) can be defined as the following equation (8).
  • the calibration function F ( ⁇ P) 0.02694 ⁇ ln ( ⁇ P) ⁇ 0.17762 (8)
  • the calibration function F ( ⁇ P) is basically the same.
  • FIG. 21 shows the result of calculating the average loss (unit: dB) using this calibration function F ( ⁇ P). From FIG. 21, there is no significant difference between the transmission loss obtained from the peak value by the method disclosed in the present embodiment and the transmission loss obtained visually from the intensity distribution of the return light. Therefore, it was verified that the transmission loss obtained by the method disclosed in the embodiment is highly reliable.
  • FIG. 22 shows the positions of the half widths of the peaks P4 to P7 in the length direction of the optical fiber (the first half width position and the second half width position), and the positions of LF and LR.
  • the least square method or the least square median method (LMedS) or the like paying attention to the peak P5, by the least square method or the least square median method (LMedS) or the like, as shown in FIG. 23, the left side of the LF in the first intensity distribution NTS1 and the second intensity distribution NTS2, LR On the right side, between LF and LR, linear approximation is performed. Then, from these straight lines, the difference in the amount of light at the position P5 is determined, and the transmission loss is calculated.
  • the least square method or the least square median method (LMedS) or the like is used to linearly approximate the left side of LF, the right side of LR, and between LF and LR. The amount of light at the position is determined and the transmission loss is calculated.
  • FIG. 24 shows the result of determining the transmission loss (linear approximate loss) at the positions of the peaks P5 and P6 in this way.
  • FIG. 24 also shows the average loss (see FIG. 21) visually determined from the intensity distribution of the return light.
  • FIG. 25 shows a case where the present invention is applied to a system that detects an abnormality in the laid state of an optical fiber when air conditioning is performed by laying the optical fiber in a server rack of a data center.
  • the optical fiber 26 is pulled out from an optical fiber cassette (not shown) installed under the floor, and the optical fiber 26 is laid on the server rack 41 using the jig 43.
  • the transmission loss is reduced.
  • the reduction in the transmission loss not only causes the problem that the temperature detection accuracy is lowered as described above, but also leads to a reduction in the life of the optical fiber 26.
  • the technology disclosed in the embodiment it is possible to detect in real time a decrease in transmission loss due to inadequate winding around the jig 43.
  • the information is notified to the administrator, and the administrator notifies the operator of the location of the installation error and the correction instruction, so that the optical fiber 26 can be correctly installed again.
  • FIGS. 26A to 26D show examples in which the abnormality detection system described in the embodiment is applied to abnormality detection of a pipe connection portion as shown in FIG.
  • FIG. 26 (a) it is assumed that a high-temperature liquid or gas flows in the main pipe 51 during plant operation.
  • the pipe 51 expands during plant operation, and the pipe 51 contracts when the plant is stopped.
  • the data processor 30 stores the loss amount of the optical fiber 26 during plant operation and shutdown.
  • reference numeral 54 denotes a tape for fixing the optical fiber 26
  • reference numeral 55 denotes a heat insulating material and a protective tube arranged around the main pipe 51.
  • the normal operation is performed when the next operation is performed as shown in FIG. 26 (c). Rather than the branch pipe 52 is pushed out. Then, as shown in FIG. 26D, when the next stop state is reached, the optical fiber 26 is pulled without returning the pushed branch pipe 52, and an abnormality is detected by the abnormality detection system.
  • FIG. 27 shows an example in which the abnormality detection system described in the embodiment is applied to the cultivation of high-quality fruits and theft prevention in a greenhouse.
  • a temperature distribution measuring device DTS
  • the temperature in the house is managed based on the measurement result.
  • the temperature distribution measuring device is also used as the light detection device 20 in FIG. 7 and connected to the data processing device 30 to be used for detecting an abnormality.
  • the abnormality detection system detects an abnormality
  • the administrator is notified of the occurrence of the abnormality along with the lighting of the patrol lamp and the start of the alarm buzzer. As a result, enormous damage can be suppressed.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

【課題】化学プラント、製油工場及び火力発電所等の施設で発生する異常を初期段階で検知できる異常検知システム及び異常検知方法を提供する。 【解決手段】異常検知システムは、光ファイバ２６と、後方散乱光検出部２０と、データ処理部３０とを有する。後方散乱光検出部２０は、光ファイバ２６の一端側及び他端側に接続され、一端側から光ファイバ２６に光を入射して後方散乱光の第１の強度分布を取得し、他端側から光ファイバ２６に光を入射して後方散乱光の第２の強度分布を取得する。データ処理部３０は、後方散乱光検出部３０で取得した第１の強度分布に対し第１のＦＩＲフィルタを作用させて得た値と、第２の強度分布に対し第２のＦＩＲフィルタを作用させて得た値との積を、光ファイバ２６の長さ方向の各位置毎に計算し、その結果から異常の有無を判定する。

Description

異常検知システム及び異常検知方法
　本発明は、異常検知システム及び異常検知方法に関する。
　化学プラント、製油工場及び火力発電所等のように大量の可燃物、爆発物又は危険物を取り扱う施設では、配管やタンクの腐食や減肉を早期に検知して、重大事故を未然に防ぐことが重要である。
　そのために、光ファイバを温度センサとして使用する温度分布測定装置（Distributed Temperature Sensor：ＤＴＳ）を備えた異常検知システムが採用されることがある。
　この種の異常検知システムでは、例えば配管やタンクの周囲に光ファイバを敷設し、光ファイバの端部を温度分布測定装置に接続する。そして、温度分布測定装置から光ファイバ内にレーザ光を照射し、光ファイバ内で発生するラマン散乱光を温度分布測定装置で検出して配管やタンク等の温度を取得し、その結果を基に異常の有無を判定する。
　化学プラント、製油工場及び火力発電所等の施設では、異常検知の遅れが重大事故につながることがあるため、異常の発生をより早い段階で検知できるシステムが望まれる。
特開平９－１８４２８号公報 特開平２－１２３３０４号公報 国際公開第２０１０／１２５７１２号
　化学プラント、製油工場及び火力発電所等の施設で発生する異常を初期段階で検知できる異常検知システム及び異常検知方法を提供することを目的とする。
　開示の技術の一観点によれば、光ファイバと、前記光ファイバの一端側及び他端側に接続され、前記一端側から前記光ファイバに光を入射して後方散乱光の第１の強度分布を取得し、前記他端側から前記光ファイバに光を入射して後方散乱光の第２の強度分布を取得する後方散乱光検出部と、前記後方散乱光検出部で取得した前記第１の強度分布に対し第１のＦＩＲ（Finite Impulse Response）フィルタを作用させて得た値と、前記第２の強度分布に対し第２のＦＩＲフィルタを作用させて得た値との積を、前記光ファイバの長さ方向の各位置毎に計算し、その結果から異常の有無を判定するデータ処理部とを有する異常検知システムが提供される。
　開示の技術の他の一観点によれば、後方散乱光検出部により光ファイバの一端側に光を入射して後方散乱光の第１の強度分布を取得し、前記後方散乱光検出部により前記光ファイバの他端側に光を入射して後方散乱光の第２の強度分布を取得する工程と、データ処理部により、前記第１の強度分布に第１のＦＩＲ（Finite Impulse Response）フィルタを作用させて得た値と、前記第２の強度分布に第２のＦＩＲフィルタを作用させて得た値との積を、前記光ファイバの長さ方向の各位置毎に計算する工程とを有する異常検知方法が提供される。
　上記一観点に係る異常検知システム及び異常検知方法によれば、化学プラント、製油工場及び火力発電所等の施設で発生する異常を初期段階で検知できる。
図１は、主配管に枝配管が溶接された部分に光ファイバを一定の張力で巻き付けた状態を示す図である。 図２は、光ファイバに中程度の曲げを加えたとき、やや強い曲げを加えたとき、及び強い曲げを加えたときの伝送損失を調べた結果を示す図である。 図３は、中程度の曲げ、やや強い曲げ及び強い曲げを具体的に示す図である。 図４は、図２に示す戻り光の強度分布を２回微分した結果を示す図である。 図５は、光ファイバに曲げ応力を加えたときのラマン散乱光の強度分布と、その強度分布を２回微分した結果とを併せて示す図である。 図６は、戻り光の強度分布と温度分布測定装置により得られる温度分布とを併せて示した図である。 図７は、実施形態に係る異常検知システムを示すブロック図である。 図８（ａ）は光ファイバの一端側とビームスプリッタとの間が光学的に接続した状態を示す図、図８（ｂ）は光ファイバの他端側とビームスプリッタとの間が光学的に接続した状態を示す図である。 図９は、異常検知方法を示すフローチャート（その１）である。 図１０は、異常検知方法を示すフローチャート（その２）である。 図１１は、第１の強度分布（ＮＴＳ１）及び第２の強度分布（ＮＴＳ２）を示す図である。 図１２（ａ），（ｂ）は、ＦＩＲフィルタの一例を示す図である。 図１３は、図１２（ａ）に示すＦＩＲフィルタを第１の強度分布（ＮＴＳ１）に畳み込みし、図１２（ｂ）に示すＦＩＲフィルタを第２の強度分布（ＮＴＳ２）に畳み込みした結果を示す図である。 図１４は、FIRNTS1値とFIRNTS2値との積を、光ファイバの長さ方向の各位置毎に計算してグラフ化した図である。 図１５（ａ），（ｂ）は、図１４の一部を拡大して示す図である。 図１６は、伝送損失を計算する第２の方法を説明する図である。 図１７は、伝送損失の表示例を示す図である。 図１８（ａ）は図１１の第１の強度分布ＮＴＳ１を２回微分した結果を示す図、図１８（ｂ）は図１８（ａ）の一部を拡大して示す図である。 図１９は、図１４，図１５から、ピーク高さ（ピーク値）を求めた結果を示す図である。 図２０は、ピーク高さと平均損失との関係を表した図である。 図２１は、校正関数Ｆ(ΔＰ)を用いて平均損失を求めた結果を示す図である。 図２２は、ピークＰ４～Ｐ７の半値幅の位置、並びにＬＦ及びＬＲの位置を示す図である。 図２３は、第１の強度分布ＮＴＳ１及び第２の強度分布ＮＴＳ２の直線近似を示す図である。 図２４は、ピークＰ５，Ｐ６の位置における伝送損失（直線近似損失）を求めた結果を示す図である。 図２５は、利用可能性１を示す図である。 図２６は、利用可能性２を示す図である。 図２７は、利用可能性３を示す図である。
　以下、実施形態について説明する前に、実施形態の理解を容易にするための予備的事項について説明する。
　実施形態に係る異常検知システムでは、曲げ応力の印加により光ファイバの伝送損失が変化することを利用して、異常を検知する。
　図１は、主配管１１に枝配管１２が溶接された部分に光ファイバ１３を一定の張力で巻き付けた状態を示す図である。光ファイバ１３は、テープ１４により枝配線１２に部分的に固定している。
　プラントの運転や停止にともなって、主配管１１及び枝配管１２内の溶液又はガスの流れが変化し、主配管１１及び枝配管１２の温度が変化する。この温度変化により主配管１１及び枝配管１２が膨張又は収縮し、光ファイバ１３に加わる曲げ応力や引張り応力が変化する。
　光ファイバ１３は、ある程度以上の曲げ応力又は引張り応力が加えられると、伝送損失が大きくなる。そこで、例えば過去の運転時又は停止時の伝送損失と現在の伝送損失とを比較することで、異常の有無を判定することが可能になる。
　図２は、横軸に距離（光ファイバの長さ方向の位置）をとり、縦軸に戻り光の強度をとって、光ファイバに中程度の曲げを加えたとき、やや強い曲げを加えたとき、及び強い曲げを加えたときの伝送損失を調べた結果を示す図である。
　なお、中程度の曲げとは図３（ａ）に示す程度の曲げ（曲げ半径１０ｍｍ程度）であり、やや強い曲げとは図３（ａ）よりも若干強い曲げ（図３（ｂ）参照）であり、強い曲げとは図３（ｂ）よりも若干強い曲げ（図３（ｃ）参照）である。また、図２では、戻り光の強度を、光ファイバの長さ方向０ｍの位置の光量を基準に規格化している。
　図２から、光ファイバの長さ方向の約３４０ｍの位置で、曲げの程度に応じた伝送損失が発生していることがわかる。
　例えば、通常運転時には光ファイバに中程度の曲げが加わっており、光ファイバの長さ方向の特定の位置で一定量の伝送損失が発生しているものとする。この場合、光ファイバの伝送損失が急激に変化したら、何らかの異常が発生したものと判定することができる。
　図２からわかるように、戻り光の強度は光ファイバの長さ方向の位置に応じて変化するため、単に戻り光の強度から異常の有無を判定することはできない。従って、異常の検知を自動化するためには、戻り光の強度の変化を検出することが重要になる。
　特許文献１には、光ファイバの接続部の位置及び接続損失を高精度で測定することを目的とし、戻り光の強度分布を２回微分することが記載されている。この手法を、異常の有無の検知に利用することが考えられる。
　図４は、図２に示す戻り光の強度分布を２回微分した結果を示す図である。図４中の一点鎖線は、３σ（σは標準偏差）の範囲を示している。
　この図４に示すように、戻り光の強度分布を２回微分すると戻り光の強度変化が強調され、伝送損失及び伝送損失が生じた位置を比較的精度よく検出することができる。
　但し、ノイズ成分を除去するためにはしきい値を３σ程度に設定することが必要になる。しきい値を３σに設定すると、中程度の曲げの部分のピークレベルはノイズレベルよりも若干高い程度であり、検出の信頼性が高いとはいえない。
　図２から、やや強い曲げが加えられた位置の前後の戻り光の光量差は約２．４％であることがわかる。つまり、この方法では、曲げ応力や引張り応力により２．４％以上の光量の変化がないと、異常の発生を精度よく検知することはできない。
　特許文献１で使用している光パルス検出器（Optical Time Domain Reflectmeter：ＯＴＤＲ）ではレイリー散乱光を使用するが、温度分布測定装置（ＤＴＳ）で使用するラマン散乱光（ストークス光及び反ストークス光）を使用しても、同様の結果が得られる。温度分布測定装置を使用すると、温度分布の測定と異常の検知とを同時に行うことができる。
　図５は、横軸に距離（光ファイバの長さ方向の位置）をとり、縦軸に光量（右側）及びしきい値（左側）をとって、曲げ応力を加えたときのラマン散乱光の強度分布と、その強度分布を２回微分した結果とを併せて示す図である。なお、図５中、ＮＴＳはストークス光による測定結果を示し、ＴＴＳは反ストークス光による測定結果を示している。また、図５中の破線は３σの範囲を示している。
　この図５から、戻り光を２回微分する方法では、ストークス光及び反ストークス光を使用しても、光量の差が２％程度ないと異常を精度よく検出することができないことがわかる。
　このように、戻り光の強度分布を２回微分して異常の検出を行う方法では、伝送損失がある程度大きくなるまで異常を検知することができない、すなわち異常を初期段階で検知できないという問題がある。工場や化学プラントなどの施設では、異常の検知の遅れが事故拡大につながるため、異常を初期段階で検知することが望まれる。
　異常を初期段階で検知できない別の理由に、温度の影響がある。ラマン散乱光の強度は温度により変化するため、戻り光の強度は光ファイバに加えられる応力と温度とに関係する。
　図６は、横軸に距離（光ファイバの長さ方向の位置）をとり、縦軸に温度（左側）及び戻り光の強度（右側）をとって、戻り光の強度分布と温度分布測定装置により得られる温度分布とを併せて示した図である。
　なお、図６において、戻り光の強度は、光ファイバの長さ方向の０ｍの位置の光量を基準に規格化している。また、ここでは光ファイバの長さ方向の約３４５ｍの位置に曲げ応力を印加している。
　図６から、光ファイバの曲げ半径を小さくして伝送損失を大きくすると、温度の測定誤差が増大することがわかる。なお、図６において、曲げ応力を印加した位置の前後の温度差は、光ファイバの敷設経路の違いによるものなので問題ない。
　例えば、ハウジング型データセンターにおいて光ファイバを用いて温度監視を行っているときに、ラックに敷設した光ファイバに加わる応力が何らかの理由により変化し、それにより戻り光の強度が変化して正確な温度を検知できないことがある。この場合、実際の温度よりも高く検知されるときと、実際の温度よりも低く検知されるときとがある。
　実際の温度よりも高く検知されるときには、異常がないにもかかわらず異常があると判定されてしまう。逆に、実際の温度よりも低く検知されるときには、許容上限温度を超えているにもかかわらず異常がないと判定されてしまう。
　以下の実施形態では、化学プラント、製油工場及び火力発電所等の施設で発生する異常を初期段階で検知できる異常検知システムについて説明する。
　（実施形態）
　図７は、実施形態に係る異常検知システムを示すブロック図である。
　本実施形態に係る異常検知システムは、ループ型光検出装置２０と、光検出装置２０から出力されるデータを処理するデータ処理装置３０とを有する。ループ型光検出装置２０は後方散乱光検出部の一例であり、データ処理装置３０はデータ処理部の一例である。
　ループ型光検出装置２０は、レーザ光源２１と、ビームスプリッタ２２と、伝送経路切替器２３と、光検出回路部２４と、演算部２５とを有し、光ファイバ２６に接続して使用する。光ファイバ２６は、その両端が伝送経路切替器２３に接続され、例えば図１のように配管１１，１２の周囲に敷設されて、テープ１４等により配管１１，１２に部分的に固定される。
　レーザ光源２１からは、所定のパルス幅のレーザ光が一定の周期で出力される。このレーザ光は、ビームスプリッタ２２を透過し、伝送経路切替器２３を介して光ファイバ２６内に進入する。
　伝送経路切替器２３は、一定の周期でレーザ光の伝送経路を切替える。すなわち、伝送経路切替器２３は、光ファイバ２６の一端側とビームスプリッタ２２との間が光学的に接続した状態（図８（ａ）参照）と、光ファイバ２６の他端側とビームスプリッタ２２との間が光学的に接続した状態（図８（ｂ）参照）とを交互に切替える。
　光ファイバ２６内に進入した光の一部は、光ファイバ２６を構成する分子により後方散乱される。後方散乱光は光ファイバ２６を戻り、伝送経路切替器２３を通ってビームスプリッタ２２に到達する。そして、ビームスプリッタ２２により反射されて、光検出回路部２４に到達する。
　光検出回路部２４には所定の波長の光を分離するフィルタ（図示せず）と、フィルタにより分離された所定の波長の光を受光する受光素子（図示せず）とが設けられている。そして、光検出回路部２４からは、受光素子で受光した光の強度に応じた信号が出力される。
　演算部２５は、コンピュータを含んで構成される。この演算部２５は、光検出回路部２４から出力される信号の経時的変化を記憶し、それらのデータをデータ処理装置３０に出力する。
　データ処理装置３０は、コンピュータを含んで構成される。そして、データ処理装置３０は、後述するように光検出装置２０から出力されるデータを処理して異常の有無を判定し、異常と判定したときには警報を発生するなどの予め設定された処理を実行する。
　ループ型光検出装置２０として、レイリー（Rayleigh）散乱光を使用する光パルス検出器（ＯＴＤＲ）を使用してもよく、ラマン散乱光（ストークス光及び反ストークス光）を使用する温度分布測定装置（ＤＴＳ）を使用してもよい。光検出装置２０として温度分布測定装置を使用した場合は、異常の検知とともに、温度分布の測定を行うこともできる。
　なお、本願発明者らは、特定の測定ポイントの温度を基準にし、伝達関数を用いてその他の測定ポイントの温度測定値を補正する温度測定方法を提案している（特許文献３等）。この方法によれば、光ファイバの長さ方向に沿って１０ｃｍ～数１０ｃｍの間隔で設定された測定ポイントの温度を精度よく検出することができる。
　以下、本実施形態係る異常検知システムにおける異常検知方法について、図９，図１０に示すフローチャートを参照して説明する。
　まず、ステップＳ１１において、光ファイバ２６の長さ方向の所定の位置に、初期化フラグ“１”を設定する。
　例えばコネクタや融着によって２本の光ファイバ２６を光学的に接続した場合、接続部では必然的に伝送損失が発生する。従って、接続部である程度の伝送損失が発生していても、異常ではない。ステップＳ１１では、接続部のように伝送損失が予測される位置に、初期化フラグ“１”を設定する。また、光ファイバ２６の特定の位置に予めある程度の応力を印加しておく場合も、その位置に初期化フラグ“１”を設定する。
　次に、ステップＳ１２において、データ処理装置３０は、光ファイバ２６の一端側にレーザ光を照射したとき（図８（ａ）参照）に得られる光ファイバの長さ方向における強度分布を取得して、第１の強度分布とする。また、データ処理装置３０は、光ファイバ２６の他端側にレーザ光を照射したとき（図８（ｂ）参照）に得られる光ファイバの長さ方向における強度分布を取得して、第２の強度分布とする。
　図１１は、横軸に距離（光ファイバ２６の長さ方向の位置）をとり、縦軸に光検出装置２０で検出される戻り光（後方散乱光）の強度（光量）をとって、第１の強度分布（ＮＴＳ１）及び第２の強度分布（ＮＴＳ２）を示す図である。
　なお、図１１の横軸は、図８（ａ）のように光検出装置２０と光ファイバ２６の一端側とが光学的に接続された状態において、レーザ光源２１側の所定の位置を基準位置（０ｍ）とし、基準位置からの距離を示している。また、図１１中にＡ１，Ａ２で示す位置における伝送損失は、光検出装置２０と光ファイバ２６の一端側及び他端側との間を光学的に接続するコネクタで発生したものである。
　図８（ａ）のように光ファイバ２６の一端側にレーザ光を照射したときの後方散乱光の強度分布（ＮＴＳ１）は、基準位置（０ｍの位置）に近いほうが後方散乱光の強度が高く、基準位置から離れるほど後方散乱光の強度が低下する。また、図８（ｂ）のように光ファイバ２６の他端側にレーザ光を照射したときの後方散乱光の強度分布（ＮＴＳ２）は、基準位置から近いほうが後方散乱光の強度が低く、基準位置から離れるほど後方散乱光の強度が高くなる。
　次に、ステップＳ１３に移行し、データ処理装置３０は、第１の強度分布（ＮＴＳ１）及び第２の強度分布（ＮＴＳ２）に対し、微分系のＦＩＲ（Finite Impulse Response：有限インパルス応答）フィルタを作用させる。
　微分系のＦＩＲフィルタとは、ユニットステップ型の微分フィルタと異なり、微分とローパスフィルタの特性を併せもつフィルタである。本実施形態で使用するＦＩＲフィルタの一例を、図１２（ａ），（ｂ）に示す。
　図１２（ａ）は第１の強度分布（ＮＴＳ１）に対し基準位置側から作用させるときに用いるＦＩＲフィルタであり、図１２（ｂ）は第２の強度分布（ＮＴＳ２）に対し基準位置側から作用させるときに用いるＦＩＲフィルタである。図１２（ｂ）のＦＩＲフィルタは、図１２（ａ）のＦＩＲフィルタに対し、正負対称の特性を有する。なお、第２の強度分布（ＮＴＳ２）に対し基準位置と反対の側から作用させる場合は、図１２（ａ）のＦＩＲフィルタを使用することができる。
　図１３は、図１２（ａ）に示すＦＩＲフィルタを、第１の強度分布（ＮＴＳ１）に畳み込み（convolution）し、図１２（ｂ）に示すＦＩＲフィルタを、第２の強度分布（ＮＴＳ２）に畳み込みした結果を示す図である。図１３では、図１２（ａ）に示すＦＩＲフィルタを第１の強度分布（ＮＴＳ１）に畳み込みした結果をFIRNTS1で示し、図１２（ｂ）に示すＦＩＲフィルタを第２の強度分布（ＮＴＳ２）に畳み込みした結果をFIRNTS2で示している。
　この図１３に示すように、第１の強度分布及び第２の強度分布にそれぞれ第１のＦＩＲフィルタ及び第２のＦＩＲフィルタを畳み込みすることにより、伝送損失が強調される。この場合、光ファイバ２６の一端側を基準位置にしているので、FIRNTS1は基準位置側の感度が高く、FIRNTS2は基準位置と反対側の感度が高い。
　次に、ステップＳ１４に移行し、データ処理装置３０は、光ファイバ２６の長さ方向の各位置Ｌ毎に、FIRNTS1値とFIRNTS2値との積を計算する。
　図１４は、FIRNTS1値とFIRNTS2値との積を、光ファイバ２６の長さ方向の各位置毎に計算してグラフ化した図である。また、図１５（ａ）は図１４の５０ｍ～１００ｍの範囲を拡大して示す図であり、図１５（ｂ）は図１４の１８００ｍ～２１００ｍの範囲を拡大して示す図である。
　これらの図１４，図１５からわかるように、FIRNTS1値とFIRNTS2値との積を演算することにより、ノイズ成分が減少し、伝送損失がより一層強調される。この例では、Ｐ１～Ｐ７の７つのピークがあり、それらの位置で伝送損失が発生していることがわかる。
　次に、ステップＳ１５に移行し、データ処理装置３０は、各位置Ｌ毎に、FIRNTS1値とFIRNTS2値との積の値と設定値とを比較する。そして、FIRNTS1値とFIRNTS2値との積の値が設定値よりも大きいときには異常ありと判定し、その位置Ｌに異常フラグ“１”を設定する。
　設定値は、例えば３σ程度に設定すればよい。但し、位置Ｌに後述するデフォルト損失が登録されている場合は、設定値を（デフォルト損失＋誤差範囲）とする。
　なお、ステップＳ１１で初期化フラグ“１”が設定された部分、すなわち光ファイバの接続部や予めある程度の応力が印加されている部分では、設定値を超える伝送損失が生じる。このため、１回目のループでは、それらの部分に異常フラグ“１”が設定される。
　次に、ステップＳ１６に移行し、データ処理装置３０は、異常フラグ“１”が設定されている位置があるか否かを判定する。異常フラグ“１”が設定されている位置があると判定した場合はステップＳ１７に進み、異常フラグ“１”が設定されている位置がないと判定した場合はステップＳ１２に戻って処理を継続する。
　ステップＳ１６からステップＳ１７に移行した場合、データ処理装置３０は、異常フラグ“１”が設定されている位置に初期化フラグ“１”が設定されているか否かを判定する。異常フラグ“１”が設定されている位置に初期化フラグ“１”が設定されていると判定した場合はステップＳ１８に移行し、初期化フラグ“１”が設定されていないと判定した場合はステップＳ２１に移行する。
　ステップＳ１７からステップＳ１８に移行した場合、データ処理装置３０は、初期化フラグ“１”が設定されている位置ＬのFIRNTS1値及びFIRNTS2値をデフォルト損失として登録し、データ処理装置３０内に記憶する。以下、位置ＬにおけるFIRNTS1値をFIRNTS1(Ｌ)と記載し、位置ＬにおけるFIRNTS2値をFIRNTS2(Ｌ)と記載する。
　その後、ステップＳ１９に移行して、データ処理装置３０は初期化フラグ“１”をリセットする。そして、ステップＳ１２に戻り、上述した処理を繰り返す。
　２回目以降のループでは、初期化フラグが既にリセットされているため、ステップＳ１５で異常フラグ“１”が設定された個所があれば、ステップＳ１７からステップＳ２１に移行する。
　ステップＳ２１に移行すると、データ処理装置３０は、例えば警報を発生して異常を知らせる。その後、ステップＳ２２に移行し、データ処理装置３０は、第１の方法により損失量の定量化計算を行う。
　第１の方法では、データの規格化を行ってベースラインに対するピーク高さから構成テーブルを基に損失量（ｄＢ）を求める。データ処理装置３０は、例えば以下に説明する方法により、損失量を計算する。
　まず、データ処理装置３０は、各位置毎に、異常検知信号Ｐ(Ｌ)を下記（１）式で算出する。
　Ｐ(Ｌ)＝FIRNTS1(Ｌ)・FIRNTS2(Ｌ)  …（１）
　異常が発生した位置では、異常検知信号Ｐ(Ｌ)がピーク波形を示す。そのピークの最大高さの位置をLalert（Ｌ＝Lalert）とし、Lalert近傍のピーク波形を除いた異常検知信号Ｐ(Ｌ)の平均値をＰaveとしたときに、実効的なピーク高さΔＰは下記（２）式のようになる。
　ΔＰ＝Ｐ(Lalert)－Ｐave  …（２）
　ここで、伝送損失Ｌoss1(Ｌ)は、校正関数をＦ(ΔＰ)として、下記（３）式により計算する。
　Ｌoss1(Ｌ)＝－１０・log（１－Ｆ(ΔＰ))　…（３）
　但し、Ｆ(ΔＰ)は、ａ，ｂを定数としたときに、下記（３）式であらわされる。
　Ｆ(ΔＰ)＝ａ・ln(ΔＰ)－ｂ　…（３）
　このようにして第１の方法で損失量Ｌoss1(Ｌ)を求めた後、ステップＳ２３に移行する。そして、ステップＳ２３において、データ処理装置３０は、第２の方法による損失量の定量化を試みる。
　第２の方法では、注目する損失箇所の両端に別の損失箇所が存在することが前提となる。別の損失箇所は、隣接するピークの位置でもよく、光ファイバの接続部のようなデフォルト損失が登録された位置でもよい。そして、異常を検知した位置の前後を最小２乗メディアン法（ＬＭｅｄＳ）等で直線近似して、その切片の差から損失（ｄＢ）を求める。
　ここでも、第１の方法と同様に、各位置毎に、異常検知信号Ｐ(Ｌ)を下記（４）式により計算する。
　Ｐ(Ｌ)＝FIRNTS1(Ｌ)・FIRNTS2(Ｌ)  …（４）
　異常が発生した位置では、異常検知信号Ｐ(Ｌ)がピーク波形を示す。そのピークの最大高さの位置をLalert（Ｌ＝Lalert）とし、Lalert近傍のピーク波形を除いた異常検知信号Ｐ(Ｌ)の平均値をＰaveとしたときに、実効的なピーク高さΔＰは下記（５）式のようになる。
　ΔＰ＝Ｐ(Lalert)－Ｐave  …（５）
　その後、ピークの半値幅の両端の位置ＬＦ，ＬＲを求める。
　例えば図１６に示すように、注目するピーク（図１６の例ではＰ６）の半値幅の両端の位置ＬＦ，ＬＲを、ＬＦN，ＬＲNとし、それよりも基準位置側のピークの半値幅の両端の２つの位置をＬＦN-1，ＬＲN-1とする。また、逆側のピークの半値幅の両端の２つの位置をＬＦN+1，ＬＲN+1とする。
　そして、ＬＲN-1とＬＦN間で、第１の強度分布（ＮＴＳ１）及び第２の強度分布（ＮＴＳ２）を直線近似する。同様に、ＬＲNとＬＦN+1との間で第１の強度分布（ＮＴＳ１）及び第２の強度分布（ＮＴＳ２）を直線近似する。
　このようにして得た４つの直線のＬ＝Lalertの位置における値を、PNTS11(Lalert)、PNTS12(Lalert)、PNTS21(Lalert)、PNTS22(Lalert)とする。この場合、Lalertの位置における各直線の規格化された差ΔPNTS1(Lalert)及びΔPNTS2(Lalert)は、下記（６）式、（７）式のようになる。
　ΔPNTS1(Lalert)＝２・abs(PNTS11(Lalert)－PNTS12(Lalert))／(PNTS11(Lalert)＋PNTS12(Lalert))  …（６）
　ΔPNTS2(Lalert)＝２・abs(PNTS21(Lalert)－PNTS22(Lalert))／(PNTS21(Lalert)＋PNTS22(Lalert))  …（７）
　これらの値を用いて、伝送損失Ｌoss2(Ｌ)を、下記（８）式により計算する。
　Ｌoss2(Lalert)＝－１０・log(１－（ΔPNTS1（Lalert)＋ΔPNTS2(Lalert))／２)　…（８）
　ＬＦN及びＬＲNは、半値幅の位置よりも若干ピーク中央から離れる位置に設定してもよい。一般的に、損失が発生している場合の光量はキンク的な特性を示すため、そのキンク的な特性を示す領域を外したほうが、近似した直線のデータに対する残差が小さくなるためである。
　その後、ステップＳ２４に移行し、データ処理装置３０は、第２の方法により伝送損失Ｌoss2(Ｌ)が得られたか否かを判定する。第２の方法により伝送損失Ｌoss2(Ｌ)が得られた場合は、ステップＳ２５に移行する。そして、ステップＳ２５において、第１の方法により得られた伝送損失Ｌoss1(Ｌ)及び第２の方法で得られた伝送損失Ｌoss2(Ｌ)を比較し、下記（９）式で確度Ａ（％）を計算する。
　Ａ＝１０B×１００　…（９）
　但し、Ｂ＝（－abs(Ｌoss1(Ｌ)－Ｌoss2(Ｌ)）÷１０）である。また、伝送損失Ｌoss1(Ｌ)及び伝送損失Ｌoss2(Ｌ)の単位はｄＢである。
　そして、データ処理装置３０は、例えばディスプレイに、異常が発生した位置と、伝送損失Ｌoss1(Ｌ)及び伝送損失Ｌoss2(Ｌ)の平均値と、確度Ａとを表示する。図１７に伝送損失の表示例を示す。図１７において、“損失ｂｙピーク”は伝送損失Ｌoss1(Ｌ)であり、“損失ｂｙ直線近似”は伝送損失Ｌoss2(Ｌ)である。
　一方、第２の方法により伝送損失Ｌoss2(Ｌ)が得られていなければ、ステップＳ２４からステップＳ２６に移行する。そして、データ処理装置３０は、例えばディスプレイに、異常が発生した位置と、伝送損失Ｌoss1(Ｌ)とを表示する。
　本実施形態に係る異常検知システムは、光ファイバに印加される応力の微小な変化を精度よく検知できる。これにより、化学プラント、製油工場及び火力発電所等の施設における異常の発生を初期段階で検知でき、事故の発生又は事故の拡大を防止できる。
　（考察）
　図１３は、図１１に示す第１の強度分布及び第２の強度分布に対し、図１２（ａ），（ｂ）のＦＩＲフィルタを畳み込みして得た図である。図１１の強度分布に対し、図１２（ａ），（ｂ）のＦＩＲフィルタを畳み込みすると、応力が印加された個所の前後における光量の差が強調される。
　図１４は、光ファイバの長さ方向の各位置でFIRNTS1値とFIRNTS2値との積を計算した結果を示している。
　第１の強度分布に含まれるノイズ成分と第２の強度分布に含まれるノイズ成分の位相が同じとは限らないので、光ファイバの長さ方向の各位置でFIRNTS1値とFIRNTS2値との積を計算することにより、応力が印加された個所の前後における光量の差がより一層強調される。
　図１８（ａ）は図１１の第１の強度分布ＮＴＳ１を２回微分した結果を示す図、図１８（ｂ）は図１８（ａ）の一部を拡大して示す図である。
　図１８（ａ），（ｂ）に示すように、強度分布を２回微分する方法では、図８（ａ）中に一点鎖線で示すノイズレベルを３σとしたときに、ノイズレベルに対し有意の差があるピークが１つしか検出されない。これに対し、実施形態に開示した方法では、図１４，図１５に示すように、７つのピークＰ１～Ｐ７を検出することができる。
　図１４，図１５から、ピーク高さ（ピーク値）を求めたものが図１９である。なお、図１９中の平均損失は、強度分布の図から目視により光量差を求めて、伝送損失に換算したものである。
　戻り光の強度分布を２回微分する方法では、Ｐ４のピークはノイズレベルに対し有意の差がなく、このピークを検出することは困難である。それに対し、実施形態に係る方法では、Ｐ１，Ｐ７のピークを確実に検出することができる。すなわち、本実施形態に係る方法は、戻り光の強度分布を２回微分する方法に比べて、少なくとも２倍以上の検出感度を有しているということができる。
　図２０は、図１９を基に、ピーク高さと平均損失との関係を表した図である。この図２０から校正関数Ｆ(ΔＰ)を、下記（８）式のように定義することができる。
　Ｆ(ΔＰ)＝0.02694・ln(ΔＰ)－0.17762  …（８）
　光ファイバ２６の種類及び光検出装置２０の種類が同じシステムでは、校正関数Ｆ(ΔＰ)は基本的に同じである。光ファイバ２６の種類及び光検出装置２０の種類が異なる場合は、実験により予め校正関数Ｆ(ΔＰ)を求めておくことが必要になる。
　校正関数Ｆ(ΔＰ)がln(ΔＰ)を含むのは、以下の理由による。すなわち、光ファイバの長さ方向のある位置における光量をＰow1としたときに、それよりもＬだけ離れた位置における光量Ｐow2は、損失係数をαとしたときに、Ｐow２＝Ｐow１・exp（－α・Ｌ）であらわされる。従って、損失係数αは、α＝(ln(Ｐow１÷Ｐow２))÷Ｌとなる。校正関数Ｆ(ΔＰ)は損失係数αに関係するため、校正関数Ｆ(ΔＰ)にはln(ΔＰ)が含まれる。
　この校正関数Ｆ(ΔＰ)を用いて平均損失（単位：ｄＢ）を求めた結果を図２１に示す。この図２１から、本実施形態に開示した方法によりピーク値から求めた伝送損失と、戻り光の強度分布から目視で求めた伝送損失との間に大きな差は認められない。従って、実施形態に開示した方法により求めた伝送損失の信頼性が高いということが検証できた。
　以下、前述の第２の方法について、より詳細に説明する。
　例えば、ＬＦ，ＬＲとして、ピーク中央から１ｍ外側の位置を採用するものとする。図２２に、光ファイバの長さ方向におけるピークＰ４～Ｐ７の半値幅の位置（前半値幅位置及び後半値幅位置）、ＬＦ及びＬＲの位置を示す。
　ここで、ピークＰ５に着目し、最小２乗法又は最小２乗メディアン法（ＬＭｅｄＳ）等により、図２３に示すように第１の強度分布ＮＴＳ１及び第２の強度分布ＮＴＳ２におけるＬＦよりも左側、ＬＲよりも右側、ＬＦとＬＲとの間を、それぞれ直線近似する。そして、これらの直線から、Ｐ５の位置の光量の差を決定し、伝送損失を算出する。これと同様に、ピークＰ６に着目し、最小２乗法又は最小２乗メディアン法（ＬＭｅｄＳ）等により、ＬＦよりも左側、ＬＲよりも右側、ＬＦとＬＲとの間をそれぞれ直線近似して、Ｐ６の位置の光量を決定し、伝送損失を計算する。
　このようにして各ピークＰ５，Ｐ６の位置における伝送損失（直線近似損失）を求めた結果を、図２４に示す。また、図２４には、戻り光の強度分布から目視で求めた平均損失（図２１参照）を併せて示している。
　この図２４からわかるように、戻り光の強度分布から目視で求めた平均損失と第２の方法で求めた伝送損失とは、ほぼ一致している。
　以下、開示した技術の利用可能性について説明する。
　（利用可能性１）
　図２５は、データセンターのサーバラックに光ファイバを敷設して空調管理を行う際に、光ファイバの敷設状態の異常を検知するシステムに適用した事例を示している。
　データセンター内にサーバラック４１が新設された場合、床下に設置された光ファイバカセット（図示せず）から光ファイバ２６を引き出し、治具４３を用いてサーバラック４１に光ファイバ２６を敷設する。
　作業者は正しく敷設したつもりでも、例えば光ファイバ２６の治具４３への巻き付けに不備があると、伝送損失が低下する。伝送損失の低下は、前述したように温度の検出精度が低下するという問題だけでなく、光ファイバ２６の寿命の低下にもつながる。
　実施形態に開示した技術を使用すると、治具４３への巻き付け不備による伝送損失の低下をリアルタイムに検知できる。その情報が管理者に通知され、管理者から作業者に敷設ミスの場所と修正指示とを通知することにより、光ファイバ２６を正しく敷設し直すことができる。
　（利用可能性２）
　図２６（ａ）～（ｄ）は、実施形態で説明した異常検知システムを図１に示したような配管接続部の異常検知に適用した事例を示している。
　図２６（ａ）のように、プラント運用中は、主配管５１内を高温の液体又はガスが流れているものとする。この場合、プラント運用中は配管５１が膨張し、プラントが停止しているときには配管５１が収縮する。データ処理装置３０には、プラント運用中及び停止中における光ファイバ２６の損失量を記憶しておく。
　なお、図２６（ａ）～（ｄ）において、５４は光ファイバ２６を固定するテープであり、５５は主配管５１の周囲に配置された断熱材及び保護管を示している。
　図２６（ｂ）のように、例えばプラントの停止にともなって主配管５１と枝配管５２との接合溶接部に金属疲労が発生した場合、図２６（ｃ）のように次に運用したときには通常よりも枝配管５２が外に押し出される。そして、図２６（ｄ）のように、次に停止状態になると、押し出された枝配管５２が戻りきらずに光ファイバ２６が引っ張られ、異常検知システムにより異常が検知される。
　このように、プラント等において配管接続部の異常を初期段階で検知することにより、重大事故の発生が回避される。
　（利用可能性３）
　図２７は、実施形態で説明した異常検知システムを、ビニールハウス内での高級果物の栽培及び盗難防止に適用した事例を示す。
　ここでは、クラウンメロンのハウス栽培において、土中温度、周囲環境温度、及び果物の温度を、温度分布測定装置（ＤＴＳ）により測定し、その測定結果に基づいてハウス内の温度を管理しているものとする。また、ここでは、温度分布測定装置を図７の光検出装置２０としても使用し、データ処理装置３０に接続して異常の検知にも使用しているものとする。
　例えば盗難者がメロンを盗難しようとする場合、メロンに巻きつけられた光ファイバ２６を解こうとする。盗難者が慎重に作業をすれば光ファイバ２６を切断することは無いが、光ファイバ２６を解こうとすると微小伝送損失の発生は避けられない。これにより、異常検知システムにおいて、異常を検知することができる。
　異常検知システムが異常を検知すると、パトランプの点灯や警報ブザーの始動と共に管理者に異常の発生が通知される。これにより、甚大な被害を抑えることができる。
                                                                                

Claims (20)

  1. 　光ファイバと、
    　前記光ファイバの一端側及び他端側に接続され、前記一端側から前記光ファイバに光を入射して後方散乱光の第１の強度分布を取得し、前記他端側から前記光ファイバに光を入射して後方散乱光の第２の強度分布を取得する後方散乱光検出部と、
    　前記後方散乱光検出部で取得した前記第１の強度分布に対し第１のＦＩＲ（Finite Impulse Response）フィルタを作用させて得た値と、前記第２の強度分布に対し第２のＦＩＲフィルタを作用させて得た値との積を、前記光ファイバの長さ方向の各位置毎に計算し、その結果から異常の有無を判定するデータ処理部と
    　を有することを特徴とする異常検知システム。
  2. 　前記第１のＦＩＲフィルタ及び前記第２のＦＩＲフィルタは、微分とローパスフィルタとの特性を併せもつことを特徴とする請求項１に記載の異常検知システム。
  3. 　前記第１のＦＩＲフィルタ及び前記第２のＦＩＲは、正負対称の特性であることを特徴とする請求項２に記載の異常検知システム。
  4. 　前記第２のＦＩＲフィルタは、前記第１のＦＩＲフィルタと同一の特性を有することを特徴とする請求項２に記載の異常検知システム。
  5. 　前記データ処理部は、前記第１の強度分布に前記第１のＦＩＲフィルタを作用させて得た値と、前記第２の強度分布に前記第２のＦＩＲフィルタを作用させて得た値との積を計算し、それにより得られたグラフのピークから、異常が発生した位置と伝送損失とを計算することを特徴とする請求項１又は２に記載の異常検知システム。
  6. 　前記グラフのピーク高さをΔＰとしたときに伝送損失は関数Ｆ(ΔＰ)であらわすことができ、前記データ処理部は前記ピーク位置Ｌにおけるデシベル表示の伝送損失Ｌoss1(Ｌ)を、Ｌoss1(Ｌ)＝－１０・log（１－Ｆ(ΔＰ))により算出することを特徴とする請求項５に記載の異常検知システム。
  7. 　前記データ処理部は、前記ピークの前後における前記第１の強度分布及び第２の強度分布をそれぞれ直線近似し、前記ピーク位置における後方散乱光の光量変化量から伝送損失を算出することを特徴とする請求項５に記載の異常検知システム。
  8. 　前記後方散乱光検出部が、ラマン散乱光を検出して前記光ファイバの長さ方向の温度分布を取得可能であることを特徴とする請求項１乃至７のいずれか１項に記載の異常検知システム。
  9. 　前記光ファイバが、配管に沿って敷設されていることを特徴とする請求項１乃至８のいずれか１項に記載の異常検知システム。
  10. 　前記光ファイバが、電子機器に敷設されていることを特徴とする請求項１乃至８のいずれか１項に記載の異常検知システム。
  11. 　前記光ファイバが、植物の周囲に敷設されていることを特徴とする請求項１乃至８のいずれか１項に記載の異常検知システム。
  12. 　後方散乱光検出部により光ファイバの一端側に光を入射して後方散乱光の第１の強度分布を取得し、前記後方散乱光検出部により前記光ファイバの他端側に光を入射して後方散乱光の第２の強度分布を取得する工程と、
    　データ処理部により、前記第１の強度分布に第１のＦＩＲ（Finite Impulse Response）フィルタを作用させて得た値と、前記第２の強度分布に第２のＦＩＲフィルタを作用させて得た値との積を、前記光ファイバの長さ方向の各位置毎に計算する工程と
    　を有することを特徴とする異常検知方法。
  13. 　前記第１のＦＩＲフィルタ及び前記第２のＦＩＲフィルタは、微分とローパスフィルタとの特性を併せもつことを特徴とする請求項１２に記載の異常検知方法。
  14. 　前記第１のＦＩＲフィルタ及び前記第２のＦＩＲは、正負対称の特性であることを特徴とする請求項１３に記載の異常検知方法。
  15. 　前記第２のＦＩＲフィルタは、前記第１のＦＩＲフィルタと同一の特性を有することを特徴とする請求項１３に記載の異常検知方法。
  16. 　前記データ処理部は、前記第１の強度分布に前記第１のＦＩＲフィルタを作用させて得た値と、前記第２の強度分布に前記第２のＦＩＲフィルタを作用させて得た値との積を計算し、それにより得られたグラフのピークから、異常が発生した位置と伝送損失とを計算することを特徴とする請求項１２又は１３に記載の異常検知方法。
  17. 　前記グラフのピーク高さをΔＰとしたときに伝送損失は関数Ｆ(ΔＰ)であらわすことができ、前記データ処理部は前記ピーク位置Ｌにおけるデシベル表示の伝送損失Ｌoss1(Ｌ)を、Ｌoss1(Ｌ)＝－１０・log（１－Ｆ(ΔＰ))により算出することを特徴とする請求項１６に記載の異常検知方法。
  18. 　前記データ処理部は、前記ピークの前後における前記第１の強度分布及び前記第２の強度分布をそれぞれ直線近似し、前記ピーク位置における後方散乱光の光量変化量から伝送損失を算出することを特徴とする請求項１６に記載の異常検知方法。
  19. 　前記データ処理部は、前記第１の強度分布に前記第１のＦＩＲフィルタを作用させて得た値と前記第２の強度分布に前記第２のＦＩＲフィルタを作用させて得た値との積の値が設定値を超えた位置があっても、当該位置が予め伝送損失が予測される位置の場合は異常とは判定しないことを特徴とする請求項１２に記載の異常検知方法。
  20. 　前記後方散乱光検出部として、温度分布測定装置を使用することを特徴とする請求項１２に記載の異常検知方法。
                                                                                    
PCT/JP2012/077354 2012-10-23 2012-10-23 異常検知システム及び異常検知方法 WO2014064771A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EP12887297.5A EP2913647B1 (en) 2012-10-23 2012-10-23 Abnormality detecting system and abnormality detecting method
PCT/JP2012/077354 WO2014064771A1 (ja) 2012-10-23 2012-10-23 異常検知システム及び異常検知方法
JP2014543051A JP5880728B2 (ja) 2012-10-23 2012-10-23 異常検知システム及び異常検知方法
US14/694,440 US9347803B2 (en) 2012-10-23 2015-04-23 Abnormality detection system and abnormality detection method
US15/149,508 US9528860B2 (en) 2012-10-23 2016-05-09 Abnormality detection system and abnormality detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/077354 WO2014064771A1 (ja) 2012-10-23 2012-10-23 異常検知システム及び異常検知方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/694,440 Continuation US9347803B2 (en) 2012-10-23 2015-04-23 Abnormality detection system and abnormality detection method

Publications (1)

Publication Number Publication Date
WO2014064771A1 true WO2014064771A1 (ja) 2014-05-01

Family

ID=50544170

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/077354 WO2014064771A1 (ja) 2012-10-23 2012-10-23 異常検知システム及び異常検知方法

Country Status (4)

Country Link
US (2) US9347803B2 (ja)
EP (1) EP2913647B1 (ja)
JP (1) JP5880728B2 (ja)
WO (1) WO2014064771A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015222225A (ja) * 2014-05-23 2015-12-10 横河電機株式会社 分布型温度測定装置
EP2977734A1 (en) * 2014-07-16 2016-01-27 Yokogawa Electric Corporation Optical fiber temperature distribution measuring device
WO2016181541A1 (ja) * 2015-05-13 2016-11-17 富士通株式会社 温度測定装置、温度測定方法および温度測定プログラム
DE102015110528A1 (de) * 2015-06-30 2017-01-05 Aiq Dienstleistungen Ug (Haftungsbeschränkt) Filtern von verteilten Erfassungsdaten

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5880728B2 (ja) * 2012-10-23 2016-03-09 富士通株式会社 異常検知システム及び異常検知方法
EP2913645B1 (en) * 2012-10-23 2017-03-22 Fujitsu Limited Abnormality detecting system and abnormality detecting method
US9512711B2 (en) * 2014-02-24 2016-12-06 Halliburton Energy Services, Inc. Portable attachment of fiber optic sensing loop
US11441926B2 (en) 2018-01-10 2022-09-13 Hifi Engineering Inc. Method and system for detecting events in a conduit

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02123304A (ja) 1988-10-31 1990-05-10 Fujikura Ltd 火災検知用光ファイバ
JPH0918428A (ja) 1995-06-30 1997-01-17 Ando Electric Co Ltd 光ファイバの試験方法
JP2584478B2 (ja) * 1988-03-14 1997-02-26 東京電力株式会社 光ファイバ後方散乱光の受信信号処理方法
JPH09304536A (ja) * 1996-05-16 1997-11-28 Hitachi Ltd 光ファイバ温度・放射線計測装置
JPH10117424A (ja) * 1996-08-23 1998-05-06 Sumitomo Electric Ind Ltd 電力ケーブルの埋設深度測定方法
JPH11183290A (ja) * 1997-12-24 1999-07-09 Toshiba Corp 光ファイバ計測システム
WO2010125712A1 (ja) 2009-05-01 2010-11-04 富士通株式会社 温度測定システム及び温度測定方法
JP2011232138A (ja) * 2010-04-27 2011-11-17 Neubrex Co Ltd 分布型光ファイバセンサ

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4588255A (en) * 1982-06-21 1986-05-13 The Board Of Trustees Of The Leland Stanford Junior University Optical guided wave signal processor for matrix-vector multiplication and filtering
JPH04332835A (ja) 1991-05-08 1992-11-19 Agency Of Ind Science & Technol 分布温度データの補正処理方法
JP2962452B2 (ja) * 1993-08-16 1999-10-12 川崎製鉄株式会社 低温設備の異常検出方法
US5642445A (en) * 1995-10-31 1997-06-24 The United States Of America As Represented By The Secretary Of The Navy System for determining an interior or exterior acoustic noise level of an enclosed structure and noise reduction device incorporating such system
WO2004104536A1 (en) * 2003-05-23 2004-12-02 Sensor Highway Limited Distributed temperature sensing system with remote reference coil
GB0420238D0 (en) * 2004-09-11 2004-10-13 Univ Strathclyde Fibre optic sensing system
JP4929949B2 (ja) * 2006-09-26 2012-05-09 住友電気工業株式会社 光ファイバ分布型センサ及び光ファイバ分布型検知方法
JPWO2008084637A1 (ja) * 2007-01-10 2010-04-30 住友電気工業株式会社 対象物の温度分布測定方法及びセンサユニット
JP5664658B2 (ja) * 2010-10-29 2015-02-04 富士通株式会社 温度測定システム及び温度測定方法
US8818199B2 (en) * 2012-02-06 2014-08-26 Adtran, Inc. Correlation systems and methods with error compensation
EP2913645B1 (en) * 2012-10-23 2017-03-22 Fujitsu Limited Abnormality detecting system and abnormality detecting method
JP5880728B2 (ja) * 2012-10-23 2016-03-09 富士通株式会社 異常検知システム及び異常検知方法
EP2913646B1 (en) * 2012-10-26 2018-04-11 Fujitsu Limited Temperature measuring system and abnormality detecting method
US10316643B2 (en) * 2013-10-24 2019-06-11 Baker Hughes, A Ge Company, Llc High resolution distributed temperature sensing for downhole monitoring

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2584478B2 (ja) * 1988-03-14 1997-02-26 東京電力株式会社 光ファイバ後方散乱光の受信信号処理方法
JPH02123304A (ja) 1988-10-31 1990-05-10 Fujikura Ltd 火災検知用光ファイバ
JPH0918428A (ja) 1995-06-30 1997-01-17 Ando Electric Co Ltd 光ファイバの試験方法
JPH09304536A (ja) * 1996-05-16 1997-11-28 Hitachi Ltd 光ファイバ温度・放射線計測装置
JPH10117424A (ja) * 1996-08-23 1998-05-06 Sumitomo Electric Ind Ltd 電力ケーブルの埋設深度測定方法
JPH11183290A (ja) * 1997-12-24 1999-07-09 Toshiba Corp 光ファイバ計測システム
WO2010125712A1 (ja) 2009-05-01 2010-11-04 富士通株式会社 温度測定システム及び温度測定方法
JP2011232138A (ja) * 2010-04-27 2011-11-17 Neubrex Co Ltd 分布型光ファイバセンサ

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2913647A4 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015222225A (ja) * 2014-05-23 2015-12-10 横河電機株式会社 分布型温度測定装置
EP2977734A1 (en) * 2014-07-16 2016-01-27 Yokogawa Electric Corporation Optical fiber temperature distribution measuring device
JP2016023955A (ja) * 2014-07-16 2016-02-08 横河電機株式会社 光ファイバ温度分布測定装置
US9835503B2 (en) 2014-07-16 2017-12-05 Yokogawa Electric Corporation Optical fiber temperature distribution measuring device
WO2016181541A1 (ja) * 2015-05-13 2016-11-17 富士通株式会社 温度測定装置、温度測定方法および温度測定プログラム
JPWO2016181541A1 (ja) * 2015-05-13 2018-03-01 富士通株式会社 温度測定装置、温度測定方法および温度測定プログラム
US10712211B2 (en) 2015-05-13 2020-07-14 Fujitsu Limited Temperature measurement device, temperature measurement method, and computer-readable non-transitory medium
AU2019202787B2 (en) * 2015-05-13 2020-10-08 Fujitsu Limited Temperature measurement device, temperature measurement method, and temperature measurement program
DE102015110528A1 (de) * 2015-06-30 2017-01-05 Aiq Dienstleistungen Ug (Haftungsbeschränkt) Filtern von verteilten Erfassungsdaten
DE102015110528B4 (de) 2015-06-30 2017-02-09 Aiq Dienstleistungen Ug (Haftungsbeschränkt) Filtern von verteilten Erfassungsdaten

Also Published As

Publication number Publication date
JP5880728B2 (ja) 2016-03-09
JPWO2014064771A1 (ja) 2016-09-05
US20160252463A1 (en) 2016-09-01
US9347803B2 (en) 2016-05-24
US9528860B2 (en) 2016-12-27
EP2913647A4 (en) 2015-11-18
EP2913647B1 (en) 2017-05-03
US20150241251A1 (en) 2015-08-27
EP2913647A1 (en) 2015-09-02

Similar Documents

Publication Publication Date Title
JP5880728B2 (ja) 異常検知システム及び異常検知方法
JP5867618B2 (ja) 異常検知システム及び異常検知方法
US8988669B2 (en) Power monitor for optical fiber using background scattering
JP6483806B2 (ja) ファイバーオプティクスを用いて行うパイプラインの完全性の監視
JP5975108B2 (ja) 温度測定システム及び異常検知方法
US10088342B2 (en) Abnormality detection system and abnormality detection method
CA2567551A1 (en) Fibre optic sensor method and apparatus
WO2019012616A1 (ja) 温度測定装置、温度測定方法および温度測定プログラム
JP6257539B2 (ja) 光ファイバ水位計測装置及び水位計測方法
JP2962452B2 (ja) 低温設備の異常検出方法
JP2016121958A (ja) 漏水検知システム、漏水検知装置及び漏水検知方法
US10401573B2 (en) Affixing fiber optic sensing element to an apparatus
WO2013112730A1 (en) Monitoring power combines
KR101634679B1 (ko) 광 센서 선형 감지기를 이용한 화재감지 시스템
JP5008029B2 (ja) 温度分布測定システム及び温度分布測定方法
JPH04276518A (ja) 光ファイバ形分布計測装置
CN117607626A (zh) 电缆局放和电缆温度在线一体化监测方法
JP2011209225A (ja) 温度分布校正機能を有する温度分布測定装置及び温度分布校正方法
Heggdal Evaluation of fiber for fault detection in a Direct Electrical Heating System
FR2873505A1 (fr) Systeme de detection perimetrique par analyse interferometrique

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12887297

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014543051

Country of ref document: JP

Kind code of ref document: A

REEP Request for entry into the european phase

Ref document number: 2012887297

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2012887297

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE