WO2019240231A1 - 漏えい調査装置、漏えい調査方法、および記録媒体 - Google Patents

漏えい調査装置、漏えい調査方法、および記録媒体 Download PDF

Info

Publication number
WO2019240231A1
WO2019240231A1 PCT/JP2019/023546 JP2019023546W WO2019240231A1 WO 2019240231 A1 WO2019240231 A1 WO 2019240231A1 JP 2019023546 W JP2019023546 W JP 2019023546W WO 2019240231 A1 WO2019240231 A1 WO 2019240231A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
measurement time
peak
cross
correlation function
Prior art date
Legal status (The legal status 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 status listed.)
Ceased
Application number
PCT/JP2019/023546
Other languages
English (en)
French (fr)
Inventor
茂樹 篠田
裕文 井上
純一郎 又賀
孝寛 久村
菊池　克
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
NEC Corp
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 NEC Corp filed Critical NEC Corp
Priority to EP19818550.6A priority Critical patent/EP3809107A4/en
Priority to US17/056,463 priority patent/US11402290B2/en
Priority to JP2020525661A priority patent/JP7052869B2/ja
Publication of WO2019240231A1 publication Critical patent/WO2019240231A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/24—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic or ultrasonic vibrations
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F17D5/02—Preventing, monitoring, or locating loss
    • F17D5/06—Preventing, monitoring, or locating loss using electric or acoustic means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/20—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using special tracer materials, e.g. dye, fluorescent material, radioactive material
    • G01M3/22—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using special tracer materials, e.g. dye, fluorescent material, radioactive material for pipes, cables or tubes; for pipe joints or seals; for valves; for welds; for containers, e.g. radiators
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/24—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic or ultrasonic vibrations
    • G01M3/243—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic or ultrasonic vibrations for pipes

Definitions

  • the present invention relates to a technique for investigating leakage of piping.
  • Patent Documents 1 and 2 disclose such a technique.
  • the leak investigation device is installed in two places where two sensors are installed at two places on the pipe line or at two points on the valve plugs connected to the pipe line.
  • the vibration waveform from the sound source propagating in the pipe is detected by two sensors.
  • the leakage investigation apparatus calculates a cross-correlation function between the two detected vibration waveforms.
  • the cross-correlation function is used to calculate a time difference when two sensors detect similar vibration waveforms.
  • the leakage investigation device specifies the position of the sound source from the time difference thus calculated.
  • the sound source may be a disturbance due to the use of water, not leakage. However, it is difficult to accurately determine the type of sound source only from the vibration waveform. Therefore, the leakage investigation apparatus related to the related technology repeatedly detects the vibration waveform a plurality of times using two sensors. In all detections, if a sound source is specified at the same position, it is determined that the sound source is leaking.
  • An object of the present invention is to provide a leak investigation device and the like that can more efficiently determine whether there is a leak.
  • a leakage investigation apparatus uses a measurement time determination unit that determines a measurement time and at least two sensors installed in a pipe. Vibration measurement means for measuring a vibration waveform, cross-correlation function calculation means for calculating a cross-correlation function of the measured vibration waveform, and peak detection for detecting a peak of the cross-correlation function at least twice during the measurement time And a leakage determining means for determining that there is a leakage when the peak is repeatedly detected for the measurement time.
  • a leakage investigation method determines a measurement time, and measures a vibration waveform for the measurement time using at least two sensors installed in a pipe. Calculating a cross-correlation function of the measured vibration waveform, detecting a peak of the cross-correlation function twice or more during the measurement time, and if the peak is repeatedly detected for the measurement time, leakage is detected. It includes determining that there is.
  • a recording medium measures a vibration waveform for the measurement time by determining the measurement time and using at least two sensors installed in a pipe. Calculating a cross-correlation function of the measured vibration waveform, detecting a peak of the cross-correlation function twice or more during the measurement time, and repeatedly detecting the peak for the measurement time.
  • a program for causing a computer to determine that there is a leak is stored.
  • the presence or absence of leakage can be determined more efficiently.
  • FIG. 1 is a diagram schematically showing a configuration of a leakage investigation system according to Embodiment 1.
  • FIG. It is a block diagram which shows the structure of the data processor with which the leakage investigation system concerning Embodiment 1 was equipped. It is a histogram figure of the usage time of water which the data processor concerning Embodiment 1 generates.
  • 3 is a flowchart showing a flow of operation of a data processor provided in the leakage investigation system according to the first embodiment. An example of data of two vibration waveforms measured by two measuring instruments is shown. An example of the cross-correlation function calculated based on the data of two vibration waveforms is shown. It is a heat map which shows the time transition of a cross correlation function.
  • FIG. 10 is a flowchart showing a flow of operation of the data processor according to the second embodiment. It is a block diagram which shows the structure of the leak investigation apparatus concerning Embodiment 3.
  • FIG. 10 is a diagram illustrating a hardware configuration of an information processing apparatus according to a fourth embodiment.
  • FIG. 1 is a diagram schematically showing a configuration of a leakage investigation system 100 according to the present embodiment.
  • the leakage investigation system 100 includes a pipe 1, a valve plug 2, two measuring instruments 3 a and 3 b, and a data processor 4.
  • Leakage investigation system 100 may include three or more measuring instruments.
  • Piping 1 is used for transporting fluids such as water.
  • the pipe 1 is provided with a plurality of valve plugs 2. However, FIG. 1 shows only two valve plugs 2. As shown in FIG. 1, the valve plug 2 is provided with two measuring instruments 3a and 3b.
  • the pipe 1 may be a part of the water supply main, or may be a water supply main or drainage pipe connected from the water supply main to a factory, facility, home, or the like.
  • the two measuring instruments 3a and 3b are devices for measuring a vibration waveform generated in the pipe 1.
  • the vibration waveform data (see FIG. 5) measured by the two measuring instruments 3 a and 3 b is transmitted to the data processor 4.
  • the data processor 4 determines the presence or absence of leakage in the pipe 1 based on the vibration waveform data received from the two measuring instruments 3a and 3b, and displays the determination result on the display unit 20.
  • the data processor 4 also determines the type of sound source. Specifically, the type of sound source is leakage and disturbance.
  • the disturbance includes, for example, the use of water.
  • the data processor 4 is an example of a leakage investigation device.
  • the measuring instruments 3a and 3b and the data processor 4 may be provided in the same device or may be provided in different devices.
  • the method of transmitting and receiving data between the measuring instruments 3a and 3b and the data processor 4 is not limited by means such as wired, wireless, and recording media.
  • the measuring instruments 3a and 3b and the data processor 4 may transmit and receive data at all times, or after a certain amount of data is accumulated in the measuring instruments 3a and 3b, the accumulated data is collectively collected in the data processor 4. May be sent to.
  • the measuring instruments 3a and 3b may detect the vibration waveform by any principle.
  • the measuring instruments 3a and 3b may include one of a vibration sensor, a water pressure sensor, and a hydrophone.
  • the measuring instruments 3a and 3b should just have a sensitivity about the frequency band of the vibration waveform which propagates the inside of the piping 1 at least.
  • the fluid flowing in the pipe 1 is water
  • the type of fluid is not particularly limited.
  • the fluid may be gas, petroleum, or other liquid or gas.
  • FIG. 2 is a diagram illustrating a configuration of the data processor 4 provided in the leakage investigation system 100.
  • the data processor 4 includes a measurement data storage unit 10, a usage time database (DB) 11, a histogram function calculation unit 12, a measurement time determination unit 13, a vibration measurement unit 14, and a cross correlation function.
  • a calculation unit 15, an integration time determination unit 16, a peak detection unit 17, a continuity determination unit 18, a leakage determination unit 19, and a display unit 20 are provided.
  • the data processor 4 may include an output unit for outputting data to an external device (not shown) instead of the display unit 20.
  • the measurement data storage unit 10 stores vibration waveform data measured by the measuring devices 3a and 3b by the vibration measuring unit 14, respectively.
  • the measurement data storage unit 10 may further store temporary processing data generated by each unit of the data processor 4 (excluding the measurement data storage unit 10 and the usage time DB 11) and result information.
  • the usage time DB 11 includes pre-recorded statistical data of the usage time of each water.
  • the water use time means the time during which water is continuously flowing.
  • the histogram function calculation unit 12 acquires statistical data of water usage time from the usage time DB 11. And the histogram function calculation part 12 calculates the histogram which shows the relationship between the usage time of water and frequency, ie, the frequency
  • FIG. 3 shows an example of a histogram generated by the histogram function calculation unit 12 from statistical data indicating the relationship between the usage time and the number of uses of water.
  • the horizontal axis represents usage time
  • the vertical axis represents frequency, that is, the number of times of use corresponding to water usage time.
  • the usage count of water decreases as the usage time increases. For example, although the usage time of water in daily life depends on the application, it is usually several seconds to several tens of seconds, and is rarely about several hours at the longest. Water usage in everyday life rarely lasts for a day or more.
  • the measurement time determination unit 13 receives a histogram from the histogram function calculation unit 12.
  • the measurement time determination unit 13 refers to the histogram generated by the histogram function calculation unit 12 and determines, as the measurement time, a time such that the probability that the sound source is likely to be a disturbance is less than a predetermined value. .
  • the measurement time determination unit 13 notifies the vibration measurement unit 14 of information indicating the determined measurement time.
  • the measurement end time t m is, since t m ⁇ t u, becomes the average value + 2 [sigma] or more or mean + 3 [sigma] or more.
  • ⁇ represents a standard deviation when it is assumed that the histogram shown in FIG. 3 follows a normal distribution.
  • the method for determining the measurement time t m is not particularly limited.
  • the measurement end time coincides with the measurement time.
  • the measurement time is also represented by the symbol t m in the same way as the measurement end time.
  • the probability that the sound source will be misidentified is called the misclassification rate.
  • the measurement time determination unit 13 may set the measurement time t m by multiplying the minimum measurement time tu by an arbitrary safety margin coefficient s (> 1). Thereby, the misclassification rate of the sound source can be further reduced.
  • the vibration measurement unit 14 receives information indicating the measurement time from the measurement time determination unit 13.
  • the vibration measuring unit 14 measures the vibration waveform of the pipe 1 using the two measuring instruments 3a and 3b. Then, the vibration measurement unit 14 accumulates the measured vibration waveform data in the data storage unit 10.
  • the vibration measurement unit 14 acquires vibration waveform data for the measurement time t m from the data storage unit 10 and transmits it to the cross-correlation function calculation unit 15.
  • Cross-correlation function calculating unit 15 the vibration measurement unit 14, receives the data of the measuring time t m fraction of the vibration waveform.
  • Cross-correlation function calculating unit 15 using the data of the vibration waveform of the measured time t m fraction received from the vibration measuring unit 14 calculates a cross-correlation function of the measured time t m min. Specifically, the cross-correlation function calculation unit 15 determines a certain time (“window time”), and slides the window time from the start time to the end time of the measurement time t m , A cross-correlation function is calculated. Details of the window time will be described later.
  • the cross-correlation function indicates the magnitude of the correlation between the vibration waveforms measured by the two measuring instruments 3a and 3b.
  • the cross-correlation function calculation unit 15 transmits information indicating the calculated cross-correlation function for the measurement time t m to the leakage determination unit 19. As described below, leakage determination unit 19, based on the cross-correlation function of the measured time t m min, to determine whether there is a sound source to the pipe 1. When the discrimination result that the sound source is present in the pipe 1 is received from the leak discrimination unit 19, the cross-correlation function calculation unit 15 transmits information indicating the calculated cross-correlation function for the measurement time t m to the integration time determination unit 16. Send to.
  • the integration time determination unit 16 receives information indicating the cross-correlation function for the measurement time t m from the cross-correlation function calculation unit 15.
  • the cross-correlation function of the measured time t m min peak by the sound source appears.
  • the peak value of the cross-correlation function obtained by the cross-correlation function calculation unit 15 is not sufficiently large with respect to noise (that is, the S / N is low) such that it can be clearly distinguished from noise.
  • the peak detector 17 described later cannot detect this peak.
  • the integration time determination unit 16 determines the time for accumulating the value of the cross-correlation function (hereinafter sometimes referred to as integration time).
  • integration time the calculation for accumulating the value of the cross-correlation function for a certain time.
  • the integration time determination unit 16 first determines the minimum integration time ⁇ T based on the following Equation 2.
  • the minimum integration time ⁇ T is the minimum value of the integration time.
  • the integration time determination unit 16 calculates the minimum integration time ⁇ T according to Equation 2 shown below.
  • Equation 2 psr is the maximum value (hereinafter simply referred to as the maximum value) of the cross-correlation function accumulated over the measurement time T (same as t m in FIG. 3 described later), and psr th is the mutual value detected as a peak. Represents the threshold value of the correlation function. When there is a peak in the integrated cross-correlation function, the maximum value psr is the value of that peak. On the other hand, when there is no peak in the integrated cross correlation function, the maximum value psr is due to noise.
  • the peak threshold value psr th is determined in advance by S / N which is a signal-to-noise ratio from the sound source.
  • the minimum integrated time ⁇ T increases as the maximum value psr decreases and as the measurement time T increases. For example, when the measurement time T is 30 minutes, psr is 20, and psr th is 6, the minimum integration time ⁇ T is 2.7 minutes from Equation 2.
  • the integration time determination unit 16 determines the integration time so as to be longer than the calculated minimum integration time ⁇ T and shorter than the measurement time (hereinafter referred to as t m ).
  • Integrated time determination section 16 the determined information indicating the accumulated time, together with information indicating a cross-correlation function of the measured time t m min, and transmits to the peak detector 17.
  • Peak detector 17, measurement time (hereinafter referred to as t m) repeatedly detecting the peak of the cross correlation function in.
  • the peak detector 17, N ( ⁇ 2) times during the measurement time t m, for detecting the peak of the cross correlation function.
  • the peak detection unit 17 determines immediately after the start of the measurement time t m as the first sampling start time.
  • the peak detector 17 integrates the cross-correlation function values from the first sampling start time over the first integration time.
  • the peak detection unit 17 calculates a maximum value exceeding the threshold value psr th from the integrated cross-correlation function in the first sampling period (that is, from the first sampling start time to after the first integration time). Detect as a peak.
  • the peak detection unit 17 determines a time after a predetermined time from the first sampling start time in the measurement time t m as the second sampling start time.
  • the predetermined time is a time obtained by dividing a time obtained by subtracting a time from the measurement start time to the first sampling start time from the measurement time t m by N ⁇ 1.
  • the peak detection unit 17 integrates the value of the cross-correlation function from the second sampling start time over the second integration time.
  • the predetermined time may be the same as the integration time, but may be shorter or longer than the integration time. That is, the predetermined time is independent of the accumulated time.
  • the peak detection unit 17 detects the number of the sampling period (or sampling start time) in which the peak of the integrated cross-correlation function was detected during the measurement time t m and the number of the sampling period (or sampling start time) in which the peak could not be detected. Is transmitted to the continuity determination unit 18.
  • the continuity determination unit 18 determines the number of sampling periods (or sampling start times) at which the peaks of the accumulated cross-correlation function were detected during the measurement time t m from the peak detection unit 17 and the sampling period during which no peak was detected. Information indicating the number (or sampling start time) is received.
  • the peak detecting section 17 determines whether the detected repeating peak during the measurement time t m. Specifically, the continuity determination section 18, during the sampling period of the peak detector 17 does not detect a peak during the measurement time t m to (sampling start time until after the integration time), and peaks were severed discriminated To do. On the other hand, the peak detector 17 is all sampling periods in the measurement time t m, when can detect the peak, respectively, the continuity determination section 18, during the measurement time t m, the peak of the cross correlation function is continued Is determined. Thus, the continuity determination section 18, during the measurement time t m, to determine whether the peak of the cross correlation function is continued. The continuity determination unit 18 transmits the determination result to the leakage determination unit 19.
  • the leakage determination unit 19 receives information indicating the cross-correlation function for the measurement time t m from the cross-correlation function calculation unit 15. When the maximum value psr of the cross-correlation function accumulated over the measurement time t m is equal to or less than the threshold value psr th , the leakage determination unit 19 determines that there is no sound source in the pipe 1. On the other hand, when the maximum value psr of the cross-correlation function accumulated over the measurement time t m exceeds the threshold value psr th , the leakage determination unit 19 determines that there is a sound source in the pipe 1.
  • the leakage determination unit 19 transmits a determination result that there is a sound source to the cross-correlation function calculation unit 15.
  • the leakage determination unit 19 receives a determination result regarding whether or not the peak of the cross-correlation function continues during the measurement time t m from the continuity determination unit 18.
  • the leakage determination unit 19 determines the presence or absence of leakage based on the continuity determination result by the continuity determination unit 18. More specifically, when the peak detection unit 17 repeatedly detects a peak over the measurement time t m , the leakage determination unit 19 determines that the sound source is leakage.
  • the leakage determination unit 19 determines that the sound source is a disturbance, that is, is not a leakage.
  • the leakage determination unit 19 causes the display unit 20 to display the determination result.
  • FIG. 4 is a flowchart showing the operation of the data processor 4 provided in the leakage investigation system 100.
  • the histogram function calculation unit 12 generates the above-described histogram (see FIG. 3) using the statistical data of the usage time of water acquired from the usage time DB 11 (S001).
  • the measurement time determination unit 13 refers to the histogram function calculated by the histogram function calculation unit 12 and determines the measurement time t m of the vibration waveform so that the misclassification rate of the sound source is equal to or less than a predetermined value. (S002). Measuring time determining unit 13, the determined information indicating the measured time t m and notifies the vibration measurement section 14.
  • the vibration measurement unit 14 receives a notification of information indicating the measurement time t m from the measurement time determination unit 13.
  • the vibration measurement unit 14 acquires vibration waveform data for the measurement time t m measured by the two measuring instruments 3a and 3b (see FIG. 1) from the data storage unit 10 (S003). Thereafter, the vibration measurement unit 14 outputs vibration waveform data for the measurement time t m to the cross-correlation function calculation unit 15.
  • the cross-correlation function calculation unit 15 calculates a cross-correlation function using the data received from the vibration measurement unit 14 (S004). More specifically, the cross-correlation function calculation unit 15 calculates a cross-correlation function based on vibration waveform data for a certain time (hereinafter referred to as window time) obtained from the two measuring instruments 3a and 3b. To do.
  • the window time is determined based on the distance between the two measuring instruments 3a and 3b, for example.
  • the cross-correlation function calculation unit 15 sets a time obtained by dividing the distance between the two measuring instruments 3a and 3b by the speed at which vibration propagates through the pipe 1 as the window time.
  • Cross-correlation function calculating unit 15 in between the start time of the measured time t m to the end time, while sliding the window time, and calculates the cross-correlation function.
  • the cross-correlation function of the measured time t m min is obtained.
  • the cross-correlation function calculation unit 15 transmits information indicating the cross-correlation function for the measurement time t m to the leakage determination unit 19.
  • the leakage determination unit 19 receives information indicating the cross-correlation function for the measurement time t m from the cross-correlation function calculation unit 15.
  • the leakage determination unit 19 integrates the cross-correlation function over the measurement time t m and extracts the above-described maximum value psr from the integrated cross-correlation function. When the maximum value psr does not exceed the threshold value psr th capable of detecting the peak of the cross-correlation function, there is no peak in the cross-correlation function during the measurement time t m . If there is no peak in the cross-correlation function integrated over the measurement time t m (No in S005), the leakage determination unit 19 determines that there is no sound source (S009A).
  • the leakage determination unit 19 notifies the cross-correlation function determination unit 15 of the determination result that there is a sound source.
  • the cross-correlation function calculation unit 15 transmits information indicating the cross-correlation function for the measurement time t m to the integration time determination unit 16.
  • the integration time determination unit 16 receives information indicating the cross-correlation function for the measurement time t m from the cross-correlation function calculation unit 15. As described above, the integration time determination unit 16 calculates the minimum integration time ⁇ T, and further determines an integration time longer than the minimum integration time ⁇ T and shorter than the measurement time t m (S006).
  • Peak detector 17, every predetermined time during the measurement time t m, for detecting the peak of the cross correlation function is integrated over the integration time (S007).
  • the predetermined time may be the same as or longer than the integration time, or may be shorter than the integration time.
  • the predetermined time is equal to or longer than the integration time, the previous integration time and the next integration time do not overlap at all (when there is no overlap).
  • the predetermined time is shorter than the accumulated time, some of them overlap (when there is an overlap).
  • the continuity determination unit 18 determines whether or not the peak of the cross-correlation function accumulated by the peak detection unit 17 has been detected in the repeated processing from the first to the Nth ( ⁇ 2) times. In other words, the continuity determination unit 18 determines whether or not the peak of the cross-correlation function exists without any interruption during the measurement time t m (S008). Note during the measurement time t m, even when the peak detector 17 continues to detect the peak, when the position of the peak varies depending on the sampling period, these peaks with different sound sources. Thus, the continuity determination section 18, during the measurement time t m, the peak detector 17 to determine whether detected repeatedly a peak at the same position.
  • the leakage determination unit 19 determines that the sound source is leakage (S009B).
  • the leakage determination unit 19 determines that the sound source is a disturbance. That is, peak detector 17, during the measurement time t m, at least once, if not detected peaks, leakage determination unit 19, the sound source is determined to be in an external disturbance.
  • the leakage determination unit 19 causes the display unit 20 to display the determination result regarding the presence / absence of leakage and the type of sound source.
  • the output unit may output the determination result by the leakage determination unit 19 to a device such as a monitor, or output it as data that can be processed by a computer. May be. Thus, the operation of the data processor 4 is completed.
  • FIG. 5 shows an example of data of two vibration waveforms measured by the vibration measuring unit 14 using the measuring instruments 3a and 3b.
  • the vertical axis of the graph shown in FIG. 5 is acceleration (unit is arbitrary), and the horizontal axis is time (unit: [s]).
  • FIG. 5 shows a characteristic correlation.
  • FIG. 6 shows the position of the sound source obtained from the cross-correlation function calculated by the cross-correlation function calculation unit 15.
  • the vertical axis of the graph shown in FIG. 6 is a value (unit is arbitrary) indicating the strength of correlation between two vibration waveforms, and the horizontal axis is the distance from the measuring instrument 3a or 3b to the sound source (unit: [m]). ).
  • the distance from the measuring instrument 3 a or 3 b to the sound source is obtained by converting the time difference between the same vibration waveforms in the two graphs shown in FIG. 5 using the speed of sound in the pipe 1.
  • the distance L1 (L + ⁇ L) / 2 from the measuring instrument 3a to the sound source
  • the value of the cross-correlation function suddenly increases at a position of a distance of about 250 m from the measuring instrument 3a or 3b. That is, there is a sound source at a distance of about 250 m from the measuring instrument 3a or 3b.
  • the cross-correlation function shown in FIG. 6 indicates the position of the sound source in a certain sampling period.
  • the continuity determination unit 18 determines whether or not a peak is repeatedly detected every predetermined time over the measurement time t m. There is a need.
  • FIG. 8 is a graph showing the relationship between the distance L (vertical axis) from the measuring instrument 3a or measuring instrument 3b to the sound source and the sampling period (horizontal axis) in which the peak of the cross-correlation function is detected. Also, the black circles in the figure indicate the sampling period during which the peak detector 17 detects the peak.
  • the leakage determination unit 19 determines that the sound source is leakage.
  • the peak due to a sound source in the position L2 is not detected for some period of measurement time t m. Therefore, the leakage determination unit 19 determines that the sound source is a disturbance due to the use of water or the like.
  • the leakage determination unit 19 can accurately determine the sound source.
  • FIG. 9 is a graph showing the relationship between the distance L (vertical axis) from the measuring instrument 3a or the measuring instrument 3b to the sound source and the sampling period (horizontal axis) in which the peak of the cross-correlation function is detected.
  • the leakage determination unit 19 determines that the sound source at the position L1 is leakage. Also relates to a sound source in the position L2, the measurement time t m, the peak is continuously being detected. Therefore, the leakage determination unit 19 erroneously determines that the sound source at the position L2 is also leakage. The reason is that the measurement time t m is too short.
  • the histogram function calculation unit 12 calculates a histogram function (see FIG. 3) from the statistical data of the water usage time. Then, the minimum measurement time tu is calculated so that the probability that the sound source is likely to be a disturbance due to the use of water is equal to or less than a predetermined value. The measurement time t m is determined so as to be longer than the minimum measurement time t u . For this reason, the leakage determination unit 19 is not likely to erroneously determine that the disturbance is leakage, and thus can correctly determine the type of the sound source.
  • FIG. 10 is a graph showing the relationship between the distance L (vertical axis) from the measuring instrument 3a or measuring instrument 3b to the sound source and the sampling period (horizontal axis) in which the peak of the cross-correlation function is detected.
  • L vertical axis
  • sampling period horizontal axis
  • the leakage determination unit 19 erroneously determines that the sound source is a disturbance. The reason is that the integration time is too short.
  • the vibration waveform (signal) from the sound source and noise are superimposed on the data generated by the measuring instruments 3a and 3b. If the integration time is too short, the signal is buried in noise, so the peak detector 17 may not be able to correctly detect the peak of the cross correlation function. Specifically, when the integration time is shorter than the above-described minimum integration time ⁇ T, the peak detection unit 17 cannot correctly detect the peak of the cross correlation function, and as a result, a situation as shown in FIG. 10 may occur. There is.
  • the above-described threshold psr th is determined based on the S / N ratio so that the vibration waveform (signal) is not buried in noise.
  • the minimum integration time ⁇ T is calculated based on the above-described equation 2. If the integration time is equal to or greater than the minimum integration time ⁇ T, the signal is not buried in noise, so that the peak detector 17 can correctly detect the peak of the cross-correlation function. Therefore, there is a low possibility that the leakage determination unit 19 erroneously determines that the leakage is a disturbance.
  • FIG. 11 is a block diagram showing the configuration of the data processor 4a provided in the leakage investigation system according to this embodiment.
  • the data processor 4a according to the present embodiment replaces the measurement time determination unit 13 included in the data processor 4 according to the first embodiment with a discrimination accuracy rate function calculation unit 21 and a measurement time.
  • a determination unit 22, a survey efficiency function calculation unit 23, an economic efficiency function calculation unit 24, and a survey plan database (DB) 25 are further provided.
  • Other configurations of the data processor 4a are the same as those of the data processor 4 according to the first embodiment.
  • the data processor 4a is an example of a leakage investigation device.
  • FIG. 12 shows a graph representing a function of the discrimination accuracy rate with respect to time and a graph representing a function of investigation efficiency with respect to time.
  • the discrimination accuracy rate indicates the probability that the data processor 4a can correctly discriminate the type of sound source as described in the first embodiment.
  • determination correct answer rate function calculation part 21 acquires the statistical data of the usage time of water from usage time DB11. And the discrimination
  • the measurement time determination unit 22 determines the measurement lower limit time t L using the discrimination accuracy rate function calculated by the discrimination accuracy rate function calculation unit 21. Specifically, the measurement time determination unit 22 determines the measurement lower limit time t L so that the discrimination accuracy rate exceeds the first predetermined value ⁇ .
  • the first predetermined value ⁇ may be determined by a user, for example.
  • the investigation efficiency function calculation unit 23 acquires information indicating a leakage investigation plan from the investigation plan DB 25.
  • the information indicating the leakage investigation plan includes information indicating the total length L (km) of the pipe 1, the scheduled investigation period Y (year), and the annual measurement days D (day / year).
  • the investigation efficiency is calculated as L / Y / D (km / day).
  • the survey efficiency decreases as the time, that is, the number of measurement days increases.
  • Measurement time determining unit 22 uses the survey efficiency function investigation efficiency function calculating unit 23 calculates and determines the measurement limit time t U. Specifically, the measurement time determining unit 22, as research efficiency does not fall below a second predetermined value beta, determines a measurement limit time t U.
  • the second predetermined value ⁇ may also be determined by the user, for example.
  • the measurement time t m may satisfy the relationship t L ⁇ t m ⁇ t U.
  • the measurement time determination unit 22 determines the measurement time t m from the viewpoint of maximizing economic efficiency, as will be described below.
  • the economic efficiency function calculation unit 24 calculates an economic efficiency function using the discrimination correct rate function and the survey efficiency function. For example, the economic efficiency function calculation unit 24 may calculate the economic efficiency function shown in FIG. 13 by weighting and adding the discrimination correct rate function and the survey efficiency function shown in FIG.
  • FIG. 13 shows an example of the economic efficiency function calculated by the economic efficiency function calculation unit 24.
  • economic efficiency function between times t L and the time t U, having a greater maximum value than either of the values in the value and the time t U at time t L. This is because, as described above, the discrimination accuracy rate increases monotonously with time, while the investigation efficiency decreases monotonously with time.
  • FIG. 14 is a flowchart showing an operation flow of the data processor 4a.
  • the processing executed by the other members of the data processor 4a is the same as the processing described in the first embodiment (see FIG. 4).
  • determination correct rate function calculation part 21 acquires the statistical data of the usage time of water from usage time DB11. And the discrimination
  • the measurement time determination unit 22 determines the measurement lower limit time t L with reference to the discrimination correct answer rate function so that the discrimination correct answer rate is equal to or greater than a predetermined value (S0022).
  • the investigation efficiency function calculation unit 23 refers to the investigation plan DB 25 and calculates an investigation efficiency function indicating the investigation efficiency of the leakage investigation (S0023).
  • Measurement time determining unit 22 refers to the determination accuracy rate function, so that the above value survey efficiency predetermined, to determine the measurement limit time t U (S0024).
  • the economic efficiency function calculation unit 24 calculates an economic efficiency function from the discrimination correct rate function and the leakage investigation efficiency function (S0025).
  • the measurement time determining unit 22 refers to the economic efficiency function and determines the measurement time t m (FIG. 13) so that the economic efficiency is maximized between the measurement lower limit time t L and the measurement upper limit time t U. (S0026).
  • the economic efficiency function is calculated from the discrimination correct rate function and the survey efficiency function. Then, the measurement time is determined so as to maximize the economic efficiency. Therefore, the presence or absence of leakage can be determined more efficiently from the viewpoint of economic efficiency.
  • FIG. 15 is a block diagram showing a configuration of a leakage investigation apparatus 4b according to the present embodiment.
  • the leakage investigation device 4 b includes a measurement time determination unit 31, a vibration measurement unit 32, a cross correlation function calculation unit 33, a peak detection unit 34, and a leakage determination unit 35.
  • the measurement time determination unit 31 determines the measurement time.
  • the vibration measuring unit 32 measures the vibration waveform with respect to the measurement time using at least two sensors installed in the pipe.
  • the cross-correlation function calculation unit 33 calculates a cross-correlation function of the measured vibration waveform.
  • the peak detector 34 detects the peak of the cross-correlation function for each integration time shorter than the measurement time.
  • the leakage determination unit 35 determines that there is a leakage when a peak is detected for each integration time with respect to the measurement time.
  • the vibration waveform is measured with respect to the measurement time using at least two sensors installed in the pipe.
  • the peak of the cross-correlation function of the measured vibration waveform is detected.
  • the peak of the cross-correlation function repeats for the measurement time, it is determined that there is a leak. Therefore, the presence or absence of leakage can be determined efficiently by only one measurement.
  • each component of each device represents a functional unit block. A part or all of each component of each device is realized by an arbitrary combination of an information processing device 900 and a program as shown in FIG. 16, for example.
  • FIG. 16 is a block diagram illustrating an example of a hardware configuration of the information processing apparatus 900 that realizes each component of each apparatus.
  • the information processing apparatus 900 includes the following configuration as an example.
  • CPU Central Processing Unit
  • ROM Read Only Memory
  • RAM Random Access Memory
  • a program 904 loaded into the RAM 903
  • a storage device 905 that stores the program 904
  • a drive device 907 that reads and writes the recording medium 906
  • a communication interface 908 connected to the communication network 909 ⁇ Input / output interface 910 for inputting / outputting data -Bus 911 connecting each component
  • Each component of each device in each embodiment is realized by the CPU 901 acquiring and executing a program 904 that realizes these functions.
  • a program 904 for realizing the function of each component of each device is stored in advance in, for example, the storage device 905 or the ROM 902, and the CPU 901 loads the RAM 903 and executes it as necessary.
  • the program 904 may be supplied to the CPU 901 via the communication network 909, or may be stored in the recording medium 906 in advance, and the drive device 907 may read the program and supply it to the CPU 901.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Examining Or Testing Airtightness (AREA)
  • Pipeline Systems (AREA)

Abstract

より効率的に漏えいの有無を判別する。振動計測部（３２）は、計測時間を決定する計測時間決定部（３１）と、配管に設置された少なくとも２つのセンサを用いて、計測時間について、振動波形を計測し、相互相関関数算出部（３３）は、計測された振動波形の相互相関関数を算出し、ピーク検出部（３４）は、計測時間中に２回以上、相互相関関数のピークを検出し、漏えい判別部（３５）は、計測時間について、ピークが繰り返し検出された場合、漏えいがあると判別する。

Description

漏えい調査装置、漏えい調査方法、および記録媒体
　本発明は、配管の漏えいを調査するための技術に関する。
　地中に埋設された配管に生じた孔またはヒビなどの空隙からの漏えいを、非掘削かつ非侵襲で調査する技術が存在する。例えば、特許文献１および２には、そのような技術が開示されている。そのような技術の一例では、漏えい調査装置は、配管の管路上の二か所か、あるいは管路と接続する弁栓類の二か所に、二つのセンサを設置して、設置された二つのセンサによって、配管内を伝搬する音源からの振動波形を検出する。
　そして、漏えい調査装置は、検出した二つの振動波形の相互相関関数（cross-correlation function）を算出する。相互相関関数は、二つのセンサが類似する振動波形を検出した時刻の差を算出するために使用される。漏えい調査装置は、このように算出される時刻の差から、音源の位置を特定する。
　音源は、漏えいではなく、水の使用等による外乱である可能性もある。しかしながら、振動波形のみから、音源の種類を正確に判別することは困難である。そこで、関連する技術に係わる漏えい調査装置は、二つのセンサを用いて、振動波形を複数回、繰り返して検出する。そして、すべての検出において、同じ位置で音源が特定された場合、音源は漏えいであると判別する。
特開平１１－２０１８５８号公報 国際公開第２０１４／０５０９２３号
　しかしながら、上述した関連する技術では、音源の種類を判別するために、時間帯や日付を変えて、振動波形を何度も検出する必要がある。そのため、調査効率が低いという問題がある。
　本発明の目的は、より効率的に漏えいの有無を判別することが可能な漏えい調査装置等を提供することにある。
　上記の課題を解決するために、本発明の一態様に係わる漏えい調査装置は、計測時間を決定する計測時間決定手段と、配管に設置された少なくとも２つのセンサを用いて、前記計測時間について、振動波形を計測する振動計測手段と、計測された前記振動波形の相互相関関数を算出する相互相関関数算出手段と、前記計測時間中に２回以上、前記相互相関関数のピークを検出するピーク検出手段と、前記計測時間について、前記ピークが繰り返し検出された場合、漏えいがあると判別する漏えい判別手段とを備えている。
　上記の課題を解決するために、本発明の一態様に係わる漏えい調査方法は、計測時間を決定し、配管に設置された少なくとも２つのセンサを用いて、前記計測時間について、振動波形を計測し、計測された前記振動波形の相互相関関数を算出し、前記計測時間中に２回以上、前記相互相関関数のピークを検出し、前記計測時間について、前記ピークが繰り返し検出された場合、漏えいがあると判別することを含む。
　上記の課題を解決するために、本発明の一態様に係わる記録媒体は、計測時間を決定することと、配管に設置された少なくとも２つのセンサを用いて、前記計測時間について、振動波形を計測し、計測された前記振動波形の相互相関関数を算出することと、前記計測時間中に２回以上、前記相互相関関数のピークを検出することと、前記計測時間について、前記ピークが繰り返し検出された場合、漏えいがあると判別することとをコンピュータに実行させるためのプログラムを格納している。
　本発明の一態様によれば、より効率的に漏えいの有無を判別することができる。
実施形態１に係わる漏えい調査システムの構成を概略的に示す図である。 実施形態１に係わる漏えい調査システムが備えたデータ処理器の構成を示すブロック図である。 実施形態１に係わるデータ処理器が生成する、水の使用時間のヒストグラム図である。 実施形態１に係わる漏えい調査システムが備えたデータ処理器の動作の流れを示すフローチャートである。 ２つの計測器によって計測される２つの振動波形のデータの一例を示す。 ２つの振動波形のデータに基づいて算出される相互相関関数の一例を示す。 相互相関関数の時間遷移を示すヒートマップである。 音源の種類の判別方法の一実施例を示す図である。 音源の種類の判別方法の比較例を示す図である。 音源の種類の判別方法の他の比較例を示す図である。 実施形態２に係わるデータ処理器の構成を示すブロック図である。 実施形態２に係わるデータ処理器が算出する判別正解率関数および調査効率関数の一例を示すグラフである。 実施形態２に係わるデータ処理器が算出する経済効率関数の一例を示すグラフである。 実施形態２に係わるデータ処理器の動作の流れを示すフローチャートである。 実施形態３に係わる漏えい調査装置の構成を示すブロック図である。 実施形態４に係わる情報処理装置のハードウェア構成を示す図である。
　以下、本発明の実施の形態について、図面を用いて説明する。尚、すべての図面において、同様な構成要素には同様の符号を付し、適宜説明を省略する。
　〔実施形態１〕
　（漏えい調査システム１００）
　図１は、本実施形態に係わる漏えい調査システム１００の構成を概略的に示す図である。図１に示すように、漏えい調査システム１００は、配管１と、弁栓２と、２つの計測器３ａ、３ｂと、データ処理器４とを含む。なお、漏えい調査システム１００は、３つ以上の計測器を備えていてもよい。
　配管１は、水などの流体を搬送するために使用される。配管１には、複数の弁栓２が設けられている。ただし、図１は２つの弁栓２のみを示している。図１に示すように、弁栓２には、２つの計測器３ａ、３ｂが設置されている。配管１は、給水本管の一部であってもよいし、給水本管から工場、施設、自宅などへ接続する給水本管または排水管であってもよい。
　２つの計測器３ａ、３ｂは、配管１において発生する振動波形を計測するための装置である。２つの計測器３ａ、３ｂが計測した振動波形のデータ（図５参照）は、データ処理器４に送信される。
　データ処理器４は、２つの計測器３ａ、３ｂから受信した振動波形のデータに基づいて、配管１における漏えいの有無を判別し、その判別結果を表示部２０に表示する。また、データ処理器４は、音源の種類も判別する。音源の種類とは、具体的には、漏えいおよび外乱のことである。外乱には、例えば、水の使用によるものが含まれる。データ処理器４は、漏えい調査装置の一例である。
　計測器３ａ、３ｂおよびデータ処理器４は、同一の機器に備えられていてもよいし、別々の機器に備えられていてもよい。
　計測器３ａ、３ｂとデータ処理器４との間におけるデータの送受信方法は、有線や無線、記録媒体等の手段によって限定されない。計測器３ａ、３ｂおよびデータ処理器４は、常時にデータを送受信してもよいし、一定量のデータが計測器３ａ、３ｂに蓄積された後、蓄積されたデータが一括でデータ処理器４に送信されてもよい。
　図１に示すように、配管１に漏えい孔５がある場合、漏えい孔５からの漏えいを音源とする振動が発生する。そのため、振動波形が配管１内を伝搬する。２つの計測器３ａ、３ｂは、同期して、この振動波形を計測する。
　計測器３ａ、３ｂは、どのような原理で振動波形を検出してもよい。例えば、計測器３ａ、３ｂは、振動センサ、水圧センサ、およびハイドロフォンのうち一つを備えていてもよい。また、計測器３ａ、３ｂは、少なくとも、配管１内を伝搬する振動波形の周波数帯域について、感度を有していればよい。
　以下では、配管１内を流れる流体が水である場合について説明する。しかしながら、流体の種類は特に限定されない。例えば、流体はガス、石油、またはその他の液体あるいは気体であってもよい。
　（データ処理器４）
　図２は、漏えい調査システム１００が備えたデータ処理器４の構成を示す図である。図２に示すように、データ処理器４は、計測データ記憶部１０、使用時間データベース（ＤＢ；Data Base）１１、ヒストグラム関数算出部１２、計測時間決定部１３、振動計測部１４、相互相関関数算出部１５、積算時間決定部１６、ピーク検出部１７、連続性判別部１８、漏えい判別部１９、および表示部２０を備えている。あるいは、データ処理器４は、表示部２０の代わりに、図示しない外部機器にデータを出力するための出力部を備えていてもよい。
　計測データ記憶部１０は、振動計測部１４によって、計測器３ａ、３ｂがそれぞれ計測した振動波形のデータを格納される。計測データ記憶部１０は、データ処理器４の各部（計測データ記憶部１０および使用時間ＤＢ１１を除く）が生成する一時処理データ、および結果情報をさらに記憶していてもよい。
　使用時間ＤＢ１１は、予め記録された、１回ごとの水の使用時間の統計データを含む。ここでは、水の使用時間とは、水が継続して流されている時間を意味する。
　ヒストグラム関数算出部１２は、使用時間ＤＢ１１から、水の使用時間の統計データを取得する。そして、ヒストグラム関数算出部１２は、取得したデータを用いて、水の使用時間と、頻度すなわち水の使用回数との関係を示すヒストグラムを算出する。ヒストグラム算出部１２は、算出したヒストグラムを計測時間決定部１３へ送信する。
　図３は、水の使用時間と使用回数との関係を示す統計データからヒストグラム関数算出部１２が生成するヒストグラムの一例を示す。図３に示すグラフでは、横軸が使用時間を表し、縦軸が頻度すなわち水の使用時間に対応する使用回数を表す。図３に示すグラフから分かるように、水の使用時間がある程度大きいならば、使用時間が大きくなるほど、水の使用回数が小さくなる。例えば、日常生活における水の使用時間は、用途にもよるが、通常は数秒から数十秒であり、長くても数時間程度になることはまれである。日常生活における水の使用時間は、１日あるいはそれ以上継続することはほとんどない。
　計測時間決定部１３は、ヒストグラム関数算出部１２からヒストグラムを受信する。計測時間決定部１３は、ヒストグラム関数算出部１２が生成したヒストグラムを参照して、音源が外乱である蓋然性を表す確率が、予め決めた所定値以下になるような時間を、計測時間として決定する。計測時間決定部１３は、決定した計測時間を示す情報を、振動計測部１４へ通知する。
　図３は、音源が外乱である蓋然性を示す確率が、予め決めた所定値以下になるような使用時間をｔｕで示す。計測時間決定部１３は、ｔｍ≧ｔｕを満たすように、計測終了時刻ｔｍを決定する。これにより、音源の誤判別、すなわち音源が外乱である場合に、音源は漏えいであると判別すること、およびその逆、が生じるリスクを低減することができる。例えば、ヒストグラムが正規分布に従うと仮定した場合、ｔｕは、ヒストグラムの平均値＋２σまたは平均値＋３σであってもよい。この場合、計測終了時刻ｔｍは、ｔｍ≧ｔｕなので、平均値＋２σ以上または平均値＋３σ以上になる。ここで、σは、図３に示すヒストグラムが正規分布に従うと仮定した場合の標準偏差を表す。しかしながら、計測時間ｔｍの決定方法は特に限定されない。
　計測開始時刻を基準（ゼロ）とすれば、計測終了時刻は、計測時間と一致する。以下では、計測時間も計測終了時刻と同じように符号ｔｍで表す。
　音源の誤判別が生じる確率を誤判別率と呼ぶ。例えば、計測時間決定部１３が、水の統計的な使用時間の９５％を含むように、最小計測時間ｔｕを決定した場合、音源の誤判別率は５％となる。なぜならば、使用時間が最小計測時間ｔｕよりも長い可能性は５％である。したがって、最小計測時間ｔｕよりも長い計測時間にわたって、音源が発生する音が継続する場合、音源が漏えいであることの確からしさは、９５％（＝１００％－５％）である。
　あるいは、計測時間決定部１３は、以下の式１に示すように、最小計測時間ｔｕに、任意の安全余裕係数ｓ（＞１）を掛けたものを、計測時間ｔｍとしてもよい。これにより、音源の誤判別率をさらに低下させることができる。
　（式１）
ｔｍ＝ｓ・ｔｕ
　振動計測部１４は、計測時間決定部１３から、計測時間を示す情報を受信する。振動計測部１４は、２つの計測器３ａ、３ｂを用いて、配管１の振動波形を計測する。そして、振動計測部１４は、計測した振動波形のデータを、データ記憶部１０に蓄積する。振動計測部１４は、計測時間ｔｍ分の振動波形のデータを、データ記憶部１０から取得して、相互相関関数算出部１５へ送信する。
　相互相関関数算出部１５は、振動計測部１４から、計測時間ｔｍ分の振動波形のデータを受信する。相互相関関数算出部１５は、振動計測部１４から受信した計測時間ｔｍ分の振動波形のデータを用いて、計測時間ｔｍ分の相互相関関数を算出する。具体的には、相互相関関数算出部１５は、一定の時間（「窓時間」）を決定し、計測時間ｔｍの開始時刻から終了時刻までの間で、上記の窓時間をスライドさせながら、相互相関関数を算出する。なお窓時間の詳細は後述する。
　相互相関関数は、２つの計測器３ａ、３ｂがそれぞれ計測した振動波形の相関の大きさを示す。相互相関関数算出部１５は、算出した計測時間ｔｍ分の相互相関関数を示す情報を、漏えい判別部１９へ送信する。後述するように、漏えい判別部１９は、計測時間ｔｍ分の相互相関関数に基づいて、配管１に音源があるかどうかを判別する。
　また、漏えい判別部１９から、配管１に音源ありという判別結果を受信した場合、相互相関関数算出部１５は、算出した計測時間ｔｍ分の相互相関関数を示す情報を、積算時間決定部１６へ送信する。
　積算時間決定部１６は、相互相関関数算出部１５から、計測時間ｔｍ分の相互相関関数を示す情報を受信する。配管１において、漏えい又は水の使用がある場合、計測時間ｔｍ分の相互相関関数には、音源によるピークが現れる。相互相関関数算出部１５によって得られた相互相関関数のピークの値が、ノイズと明確に区別できるほど、ノイズに対して十分に大きくない（すなわちＳ／Ｎが低い）場合がある。その結果、後述するピーク検出部１７は、このピークを検出することができない。しかし、ある時間分の相互相関関数の値を累積すれば、相互相関数のピークは、ノイズに対して、徐々に大きくなり、目立ってくる。そこで、積算時間決定部１６は、相互相関関数の値を累積する時間（以下、積算時間と呼ぶ場合がある）を決定する。ここでは、ある時間分の相互相関関数の値を累積する計算を積算と呼ぶ。
　具体的には、積算時間決定部１６は、まず、以下の式２に基づいて、最小積算時間ΔＴを決定する。最小積算時間ΔTは、積算時間の最小値である。まず、積算時間決定部１６は、以下に示す式２にしたがって、最小積算時間ΔＴを算出する。
　（式２）
Figure JPOXMLDOC01-appb-I000001
　式２において、ｐｓｒは、計測時間Ｔ（後述する図３のｔｍと同じ）にわたって積算した相互相関関数の最大値（以下、単に最大値と呼ぶ）、ｐｓｒｔｈは、ピークとして検出される相互相関関数の閾値、をそれぞれ表す。積算した相互相関関数にピークがある場合、最大値ｐｓｒは、そのピークの値になる。一方、積算した相互相関関数にピークがない場合、最大値ｐｓｒはノイズによるものである。ピーク閾値ｐｓｒｔｈは、音源からのシグナル対ノイズの比であるＳ／Ｎによって、予め定められる。式２によれば、最大値ｐｓｒが小さいほど、また計測時間Ｔが長いほど、最小積算時間ΔＴは大きくなる。例えば、計測時間Ｔが３０分であり、ｐｓｒが２０、ｐｓｒｔｈが６である場合、式２から、最小積算時間ΔＴは、２．７分となる。
　積算時間決定部１６は、このように算出した最小積算時間ΔＴよりも長く、かつ計測時間（以下、ｔｍと記載する）より短くなるように、積算時間を決定する。積算時間決定部１６は、決定した積算時間を示す情報を、計測時間ｔｍ分の相互相関関数を示す情報とともに、ピーク検出部１７へ送信する。
　ピーク検出部１７は、積算時間決定部１６から、積算時間を示す情報を、計測時間ｔｍ分の相互相関関数を示す情報を受信する。ピーク検出部１７は、計測時間（以下、ｔｍと記載する）中に繰り返し相互相関関数のピークを検出する。具体的には、ピーク検出部１７は、計測時間ｔｍ中にＮ（≧２）回、相互相関関数のピークを検出する。
　例えば、ピーク検出部１７は、計測時間ｔｍの開始直後を、第１回目のサンプリング開始時刻として決定する。ピーク検出部１７は、第１回目のサンプリング開始時刻から、第１回目の積算時間にわたって、相互相関関数の値を積算する。そして、ピーク検出部１７は、積算した相互相関関数から、閾値ｐｓｒｔｈを超える最大値を、第１回目のサンプリング期間（すなわち第１回目のサンプリング開始時刻から第１回目の積算時間後まで）におけるピークとして検出する。
　次に、ピーク検出部１７は、計測時間ｔｍ中の第１回目のサンプリング開始時刻から、所定の時間だけ後の時刻を、第２回目のサンプリング開始時刻として決定する。例えば、所定の時間は、計測開始から第１回目のサンプリング開始時刻までの時間を計測時間ｔｍから差し引いた時間を、Ｎ－１で割った時間である。
　ピーク検出部１７は、同様に、第２回目のサンプリング開始時刻から、第２回目の積算時間にわたって、相互相関関数の値を積算する。ピーク検出部１７は、積算した相互相関関数から、閾値ｐｓｒｔｈを超える最大値を、第２回目のサンプリング期間（すなわち第２回目のサンプリング開始時刻から第２回目の積算時間後まで）におけるピークとして検出する。
　このようにして、ピーク検出部１７は、計測時間ｔｍ分の相互相関関数のうち、第ｎ（ｎ＝１・・・，Ｎ）回目のサンプリング開始時刻から積算時間分の相互相関関数を抽出し、抽出した相互相関関数の値を積算時間にわたって積算する処理、および、積算した相互相関関数のピークを検出する処理を、計測時間ｔｍ分の相互相関関数に対し、Ｎ（≧２）回繰り返し行う。
　なお、所定の時間は、積算時間と同一であってもよいが、積算時間より短くあるいは長くてもよい。すなわち、所定の時間は、積算時間とは無関係である。
　ピーク検出部１７は、計測時間ｔｍ中、積算した相互相関関数のピークを検出したサンプリング期間の番号（あるいはサンプリング開始時刻）と、ピークを検出できなかったサンプリング期間の番号（あるいはサンプリング開始時刻）とを示す情報を、連続性判別部１８へ送信する。
　連続性判別部１８は、ピーク検出部１７から、計測時間ｔｍ中、積算した相互相関関数のピークを検出したサンプリング期間の番号（あるいはサンプリング開始時刻）と、ピークを検出できなかったサンプリング期間の番号（あるいはサンプリング開始時刻）とを示す情報を受信する。連続性判別部１８から受信した情報に基づき、連続性判別部１８は、ピーク検出部１７が計測時間ｔｍ中に繰り返しピークを検出しているかどうかを判別する。
　具体的には、連続性判別部１８は、ピーク検出部１７が計測時間ｔｍ中にピークを検出できなかったサンプリング期間中（サンプリング開始時刻から積算時間後まで）に、ピークが断絶したと判別する。一方、ピーク検出部１７が計測時間ｔｍ中の全てのサンプリング期間において、それぞれピークを検出できた場合、連続性判別部１８は、計測時間ｔｍ中、相互相関関数のピークが継続していると判別する。これにより、連続性判別部１８は、計測時間ｔｍ中、相互相関関数のピークが継続しているかどうかを判別する。連続性判別部１８は、その判別結果を漏えい判別部１９へ送信する。
　漏えい判別部１９は、相互相関関数算出部１５から、計測時間ｔｍ分の相互相関関数を示す情報を受信する。計測時間ｔｍにわたって積算した相互相関関数の最大値ｐｓｒが閾値ｐｓｒｔｈ以下である場合、漏えい判別部１９は、配管１に音源なしと判別する。一方、計測時間ｔｍにわたって積算した相互相関関数の最大値ｐｓｒが閾値ｐｓｒｔｈを超える場合、漏えい判別部１９は、配管１に音源ありと判別する。この場合、漏えい判別部１９は、音源ありという判別結果を、相互相関関数算出部１５へ送信する。
　また漏えい判別部１９は、連続性判別部１８から、計測時間ｔｍ中、相互相関関数のピークが継続しているかどうかに関しての判別結果を受信する。漏えい判別部１９は、連続性判別部１８による連続性の判別結果に基づいて、漏えいの有無を判別する。より詳細には、計測時間ｔｍにわたって、ピーク検出部１７が繰り返しピークを検出した場合、漏えい判別部１９は、音源は漏えいであると判別する。一方、計測時間ｔｍ中に少なくとも１度、ピーク検出部１７がピークを検出していない場合、漏えい判別部１９は、音源は外乱である、すなわち漏えいではないと判別する。漏えい判別部１９は、その判別結果を、表示部２０に表示させる。
　（動作フロー）
　図４は、漏えい調査システム１００が備えたデータ処理器４の動作を示すフローチャートである。図４に示すように、まず、ヒストグラム関数算出部１２は、使用時間ＤＢ１１から取得した水の使用時間の統計データを用いて、上述したヒストグラム（図３参照）を生成する（Ｓ００１）。
　計測時間決定部１３は、ヒストグラム関数算出部１２が算出したヒストグラム関数を参照して、音源の誤判別率が予め定められた所定値以下になるように、振動波形の計測時間ｔｍを決定する（Ｓ００２）。計測時間決定部１３は、決定した計測時間ｔｍを示す情報を、振動計測部１４へ通知する。
　振動計測部１４は、計測時間決定部１３から、計測時間ｔｍを示す情報の通知を受信する。振動計測部１４は、データ記憶部１０から、２つの計測器３ａ、３ｂ（図１参照）により計測された計測時間ｔｍ分の振動波形のデータを取得する（Ｓ００３）。その後、振動計測部１４は、計測時間ｔｍ分の振動波形のデータを、相互相関関数算出部１５へ出力する。
　相互相関関数算出部１５は、振動計測部１４から受信したデータを用いて、相互相関関数を算出する（Ｓ００４）。より詳細には、相互相関関数算出部１５は、２つの計測器３ａ、３ｂから得られた一定の時間（以下、窓時間と呼ぶ）分の振動波形のデータに基づいて、相互相関関数を算出する。窓時間は、例えば、２つの計測器３ａ、３ｂの距離に基づいて決められる。相互相関関数算出部１５は、２つの計測器３ａ、３ｂの距離を、配管１中を振動が伝播する速度で割って得られる時間を、上記の窓時間とする。
　相互相関関数算出部１５は、計測時間ｔｍの開始時刻から終了時刻までの間で、上記の窓時間をスライドさせながら、相互相関関数を算出する。これにより、計測時間ｔｍ分の相互相関関数が得られる。相互相関関数算出部１５は、計測時間ｔｍ分の相互相関関数を示す情報を、漏えい判別部１９へ送信する。
　漏えい判別部１９は、相互相関関数算出部１５から、計測時間ｔｍ分の相互相関関数を示す情報を受信する。漏えい判別部１９は、計測時間ｔｍにわたって、相互相関関数を積算し、積算した相互相関関数から、上述した最大値ｐｓｒを抽出する。最大値ｐｓｒが、相互相関関数のピークを検出可能な閾値ｐｓｒｔｈを超えない場合、計測時間ｔｍ中の相互相関関数には、ピークが１つも存在しない。計測時間ｔｍにわたって積算した相互相関関数にピークが存在しない場合（Ｓ００５でＮｏ）、漏えい判別部１９は、音源なしと判別する（Ｓ００９Ａ）。
　一方、最大値ｐｓｒが閾値ｐｓｒｔｈを超える場合、計測時間ｔｍ中、相互相関関数にピークが存在する。相互相関関数にピークが存在する場合（Ｓ００５でＹｅｓ）、漏えい判別部１９は、音源ありという判別結果を、相互相関関数判別部１５へ通知する。この場合、相互相関関数算出部１５は、計測時間ｔｍ分の相互相関関数を示す情報を、積算時間決定部１６へ送信する。
　積算時間決定部１６は、相互相関関数算出部１５から、計測時間ｔｍ分の相互相関関数を示す情報を受信する。積算時間決定部１６は、上述したように、最小積算時間ΔＴを算出し、最小積算時間ΔＴよりも長く、かつ計測時間ｔｍより短い積算時間をさらに決定する（Ｓ００６）。
　ピーク検出部１７は、計測時間ｔｍ中の所定の時間ごとに、積算時間にわたって積算された相互相関関数のピークを検出する（Ｓ００７）。本実施形態１では、ピーク検出部１７は、積算時間算出部１６から受信した計測時間ｔｍ分の相互相関関数のうち、第ｎ（ｎ＝１・・・，Ｎ）回目のサンプリング開始時刻から積算時間分の相互相関関数の値を積算する処理と、積算した相互相関関数のピークを検出する処理とを、Ｎ回繰り返し実行する。
　上述したように、所定の時間は、積算時間と同じ長さかあるいは積算時間よりも長くてもよいし、積算時間より短くてよい。所定の時間が積算時間と同じ長さかあるいは積算時間よりも長い場合、前の積算時間と次の積算時間とがまったく重ならない（オーバーラップなしの場合）。一方、所定の時間が積算時間よりも短い場合、それらの一部が重なる（オーバーラップありの場合）。
　連続性判別部１８は、第１回から第Ｎ（≧２）回までの繰り返しの処理において、ピーク検出部１７が積算した相互相関関数のピークを検出したかどうかを判別する。換言すれば、連続性判別部１８は、計測時間ｔｍ中、相互相関関数のピークが１度も断続することなく存在するのかどうかを判別する（Ｓ００８）。
　なお計測時間ｔｍ中、ピーク検出部１７がピークを検出し続けた場合であっても、ピークの位置がサンプリング期間によって異なる場合、これらのピークは異なる音源による。したがって、連続性判別部１８は、計測時間ｔｍ中、ピーク検出部１７が同じ位置でピークを繰り返し検出したかどうかを判別する。
　計測時間ｔｍ中、相互相関関数のピークが継続的に繰り返す場合（Ｓ００８でＹｅｓ）、漏えい判別部１９は、音源が漏えいであると判別する（Ｓ００９Ｂ）。
　一方、計測時間ｔｍ中、相互相関関数のピークが継続していない場合（Ｓ００８でＮｏ）、漏えい判別部１９は、音源が外乱である、すなわち漏えいではないと判別する（Ｓ００９Ｃ）。また、計測時間ｔｍ中、相互相関関数のピークが断続的に繰り返す場合も、音源が外乱であると、漏えい判別部１９は判別する。つまり、ピーク検出部１７が、計測時間ｔｍ中に、少なくとも１回、ピークを検出しなかった場合、漏えい判別部１９は、音源が外乱であると判別する。
　漏えい判別部１９は、漏えいの有無、および音源の種類についての判別結果を、表示部２０に表示させる。あるいは、データ処理器４が出力部を備えている場合、出力部は、漏えい判別部１９による判別結果を、モニタ等のデバイスに出力してもよいし、あるいはコンピュータが処理可能なデータとして出力してもよい。以上で、データ処理器４の動作は終了する。
　（相互相関関数）
　図５は、振動計測部１４が計測器３ａ、３ｂを用いて計測する２つの振動波形のデータの一例を示す。図５に示すグラフの縦軸は加速度（単位は任意）であり、横軸は時間（単位：［ｓ］）である。漏えいに基づく振動波形、または、水の使用等による外乱に基づく振動波形のうち少なくとも一方が存在する場合、図５に示す振動波形のデータは、特徴的な相関を示す。
　図６は、相互相関関数算出部１５が算出する相互相関関数から得られる音源の位置を示す。図６に示すグラフの縦軸は、２つの振動波形の相関の強さを表す値（単位は任意）であり、横軸は、計測器３ａまたは３ｂから音源までの距離（単位：［ｍ］）である。図６に示すグラフにおいて、計測器３ａまたは３ｂから音源までの距離は、配管１内での音速を用いて、図５に示す２つのグラフにおける同じ振動波形の時間差を変換することによって得られる。具体的には、計測器３ａから音源までの距離Ｌ１と、計測器３ｂから音源までの距離Ｌ２との差ΔＬ（＝Ｌ１－Ｌ２）は、（管内の音速）×（同じ振動波形の時間差）によって得られる。そして、計測器３ａから計測器３ｂまでの距離をＬ（＝Ｌ１＋Ｌ２）とすると、計測器３ａから音源までの距離Ｌ１＝（Ｌ＋ΔＬ）／２、計測器３ｂから音源までの距離Ｌ２＝（Ｌ－ΔＬ）／２が得られる。
　図６では、計測器３ａまたは３ｂからの距離約２５０ｍの位置において、相互相関関数の値が突然大きくなる。すなわち、計測器３ａまたは３ｂからの距離約２５０ｍの位置に音源がある。図６に示す相互相関関数は、あるサンプリング期間における音源の位置を示している。
　図７は、相互相関関数の時間遷移を示すヒートマップである。図７に示すグラフは、縦軸が距離（単位：［ｍ］）であり、横軸が時間（単位：［ｓ］）である。図７では、計測器３ａまたは３ｂから距離約５５０ｍの位置に、最大値が存在する。したがって、計測器３ａまたは３ｂから距離約５５０ｍの位置には、漏えいに基づく音源、または、外乱に基づく音源が存在する。この音源が、漏えいに基づくか、または、外乱に基づくかを判別するためには、計測時間ｔｍにわたって、所定の時間ごとにピークが繰り返し検出されたかどうかを、連続性判別部１８が判別する必要がある。
　（実施例）
　図８を参照して、音源の種類の判別方法の一実施例を説明する。図８は、計測器３ａまたは計測器３ｂから音源までの距離Ｌ（縦軸）と、相互相関関数のピークが検出されたサンプリング期間（横軸）との関係を示すグラフである。また、図中の黒丸は、ピーク検出部１７がピークを検出したサンプリング期間を示す。Ｌ１およびＬ２で示す位置には、それぞれ音源がある。ここでは、位置Ｌ１にある音源は漏えいであり、位置Ｌ２にある音源は、水の使用等による外乱である。
　図８に示すように、位置Ｌ１にある音源によるピークは、所定の時間ごとに、計測時間ｔｍにわたって検出されている。したがって、漏えい判別部１９は、音源が漏えいであると判別する。一方、位置Ｌ２にある音源によるピークは、計測時間ｔｍの一部の期間では検出されていない。そのため、漏えい判別部１９は、音源が、水の使用等による外乱であると判別する。
　この例のように、積算時間および計測時間がどちらも適切である場合、漏えい判別部１９は、音源を正確に判別することができる。
　（比較例１；計測時間が短すぎる場合）
　図９を参照して、計測時間ｔｍが適切ではない場合に起こり得る状況の一例を説明する。図９は、計測器３ａまたは計測器３ｂから音源までの距離Ｌ（縦軸）と、相互相関関数のピークが検出されたサンプリング期間（横軸）との関係を示すグラフである。ここでは、位置Ｌ＝Ｌ１にある音源は漏えいであり、位置Ｌ＝Ｌ２にある音源は、水の使用による外乱である。
　図９に示すように、位置Ｌ１にある音源に関し、計測時間ｔｍについて、連続してピークが検出されている。そのため、漏えい判別部１９は、位置Ｌ１にある音源が漏えいであると判別する。また、位置Ｌ２にある音源に関し、計測時間ｔｍについて、連続してピークが検出されている。そのため、漏えい判別部１９は、位置Ｌ２にある音源も漏えいであると誤判別する。その理由は、計測時間ｔｍが短すぎるためである。
　図３を参照すると、計測時間ｔｍが最小計測時間ｔｕよりも短い場合、図９に示すような状況が生じる可能性がある。
　本実施形態では、水の使用時間の統計データから、ヒストグラム関数算出部１２がヒストグラム関数（図３参照）を算出する。そして、音源が水の使用による外乱である蓋然性を表す確率が、予め決めた所定値以下になるように、最小計測時間ｔｕが算出される。計測時間ｔｍは、最小計測時間ｔｕよりも長くなるように決定される。そのため、漏えい判別部１９は、外乱を漏えいであると誤判別する可能性は低いので、音源の種類を正しく判別することができる。
　（比較例２；積算時間が短すぎる場合）
　図１０を参照して、積算時間が適切ではない場合に起こり得る状況の一例を説明する。図１０は、計測器３ａまたは計測器３ｂから音源までの距離Ｌ（縦軸）と、相互相関関数のピークが検出されたサンプリング期間（横軸）との関係を示すグラフである。ここでは、位置Ｌ＝Ｌ１にある音源は漏えいであり、位置Ｌ＝Ｌ２にある音源は、水の使用による外乱である。
　図１０に示すように、位置Ｌ１にある音源に関し、計測時間ｔｍの一部の期間について、ピークが検出されていない。そのため、漏えい判別部１９は、音源が外乱であると誤判別する。その理由は、積算時間が短すぎるためである。
　計測器３ａ、３ｂが生成するデータには、音源からの振動波形（シグナル）とノイズとが重畳している。積算時間が短すぎる場合、シグナルがノイズに埋没するため、ピーク検出部１７は、相互相関関数のピークを正しく検出できない可能性がある。具体的には、積算時間が、上述した最小積算時間ΔＴよりも短い場合、ピーク検出部１７は相互相関関数のピークを正しく検出できず、その結果、図１０に示すような状況が生じる可能性がある。
　本実施形態では、振動波形（シグナル）がノイズに埋没しないように、Ｓ／Ｎ比に基づいて、上述した閾値ｐｓｒｔｈを定めている。最小積算時間ΔＴは、上述した式２に基づいて算出される。積算時間が最小積算時間ΔＴ以上であれば、シグナルがノイズに埋没することはないため、ピーク検出部１７は、相互相関関数のピークを正しく検出することができる。したがって、漏えい判別部１９は、漏えいを外乱であると誤判別する可能性は低い。
　（本実施形態の効果）
　本実施形態の構成によれば、計測時間および積算時間が適切であれば、相互相関関数のピークを正しく検出することができる。そして、計測時間について、相互相関関数のピークが繰り返すかどうかに基づいて、音源の種類を判別する。したがって、１回の計測であっても、漏えいの有無を、効率的に判別することができる。
　〔実施形態２〕
　前記実施形態では、水の使用時間のヒストグラムに基づいて、計測時間を決定する構成を説明した。本実施形態では、経済効率の観点から、計測時間を決定する構成を説明する。なお、本実施形態に係わる漏えい調査システムの構成は、前記実施形態１に係わる漏えい調査システム１００（図１参照）と同じである。
　（データ処理器４ａ）
　図１１は、本実施形態に係わる漏えい調査システムが備えたデータ処理器４ａの構成を示すブロック図である。図１１に示すように、本実施形態に係わるデータ処理器４ａは、前記実施形態１に係わるデータ処理器４が備えた計測時間決定部１３に代えて、判別正解率関数算出部２１、計測時間決定部２２、調査効率関数算出部２３、経済効率関数算出部２４、および調査計画データベース（ＤＢ）２５をさらに備えている。データ処理器４ａのその他の構成は、前記実施形態１に係わるデータ処理器４と同じである。データ処理器４ａは、漏えい調査装置の一例である。
　（判別正解率と調査効率）
　図１２を参照して、判別正解率および調査効率について説明する。図１２は、時間に対する判別正解率の関数を表すグラフと、時間に対する調査効率の関数を表すグラフとを示す。判別正解率は、前記実施形態１で説明したように、データ処理器４ａが音源の種類を正しく判別できる確率を示す。
　判別正解率関数算出部２１は、使用時間ＤＢ１１から、水の使用時間の統計データを取得する。そして、判別正解率関数算出部２１は、使用時間ＤＢ１１から取得したデータを用いて、時間に対する判別正解率の高さを示す判別正解率関数を算出する。図１２に示すように、判別正解率は、時間とともに増加する。
　計測時間決定部２２は、判別正解率関数算出部２１が算出した判別正解率関数を用いて、計測下限時間ｔＬを決定する。具体的には、計測時間決定部２２は、判別正解率が第１の所定値αを超えるように、計測下限時間ｔＬを決定する。第１の所定値αは、例えば、ユーザによって決定されてよい。
　調査効率関数算出部２３は、調査計画ＤＢ２５から、漏えいの調査計画を示す情報を取得する。例えば、漏えいの調査計画を示す情報には、配管１の総延長Ｌ（ｋｍ）、調査予定期間Ｙ（年）、および年間の計測日数Ｄ（日／年）を示す各情報が含まれる。この場合、調査効率は、Ｌ／Ｙ／Ｄ（ｋｍ／日）と算出される。図１２に示すように、調査効率は、時間すなわち計測日数が長くなるほど低下する。
　計測時間決定部２２は、調査効率関数算出部２３が算出した調査効率関数を用いて、計測上限時間ｔＵを決定する。具体的には、計測時間決定部２２は、調査効率が第２の所定値βを下回らないように、計測上限時間ｔＵを決定する。第２の所定値βも、例えば、ユーザによって決定されてよい。
　図１２から分かるように、判別正解率を高くするために、計測時間ｔｍを長くするほど、調査効率が低下する。一方、調査効率を高くするために、計測時間ｔｍを短くするほど、判別正解率が低下する。すなわち、判別正解率と調査効率とはトレードオフの関係にある。計測時間ｔｍは、ｔＬ≦ｔｍ≦ｔＵの関係を満たせばよい。本実施形態では、計測時間決定部２２は、以下で説明するように、経済効率を最大化する観点で、計測時間ｔｍを決定する。
　（経済効率）
　経済効率関数算出部２４は、判別正解率関数および調査効率関数を用いて、経済効率関数を算出する。例えば、経済効率関数算出部２４は、図１２に示す判別正解率関数と調査効率関数とを、それぞれ重み付けして足し合わせることによって、図１３に示す経済効率関数を算出してもよい。
　図１３は、経済効率関数算出部２４が算出する経済効率関数の一例を示す。図１３に示すように、経済効率関数は、時刻ｔＬと時刻ｔＵとの間で、時刻ｔＬにおける値および時刻ｔＵにおける値のどちらよりも大きい極大値を有する。上述したように、判別正解率は時間とともに単調に増大する一方、調査効率は時間とともに単調に低下するからである。
　（動作フロー）
　図１４を参照して、本実施形態に係わるデータ処理器４ａの動作を説明する。図１４は、データ処理器４ａの動作の流れを示すフローチャートである。本実施形態では、判別正解率関数算出部２１、計測時間決定部２２、および調査効率関数算出部２３が係わる処理のみを説明する。データ処理器４ａの他の部材が実行する処理は、前記実施形態１で説明した処理（図４参照）と同じである。
　図１４に示すように、判別正解率関数算出部２１は、使用時間ＤＢ１１から、水の使用時間の統計データを取得する。そして、判別正解率関数算出部２１は、使用時間ＤＢ１１から取得したデータを用いて、音源の種類を正しく判別できる確率を示す判別正解率関数を算出する（Ｓ００２１）。
　計測時間決定部２２は、判別正解率関数を参照して、判別正解率が予め決めた値以上になるように、計測下限時間ｔＬを決定する（Ｓ００２２）。
　調査効率関数算出部２３は、調査計画ＤＢ２５を参照して、漏えい調査の調査効率を示す調査効率関数を算出する（Ｓ００２３）。
　計測時間決定部２２は、判別正解率関数を参照して、調査効率が予め決めた値以上になるように、計測上限時間ｔＵを決定する（Ｓ００２４）。
　経済効率関数算出部２４は、判別正解率関数および漏えい調査効率関数から、経済効率関数を算出する（Ｓ００２５）。
　計測時間決定部２２は、経済効率関数を参照して、計測下限時間ｔＬと計測上限時間ｔＵとの間で、経済効率が最大になるように、計測時間ｔｍ（図１３）を決定する（Ｓ００２６）。
　（本実施形態の効果）
　本実施形態の構成によれば、判別正解率関数および調査効率関数から、経済効率関数を算出する。そして、経済効率を最大化するように、計測時間を決定する。そのため、経済効率の観点で、より効率的に漏えいの有無を判別することができる。
　〔実施形態３〕
　（漏えい調査装置４ｂ）
　図１５は、本実施形態に係わる漏えい調査装置４ｂの構成を示すブロック図である。図１５に示すように、漏えい調査装置４ｂは、計測時間決定部３１と、振動計測部３２と、相互相関関数算出部３３と、ピーク検出部３４と、漏えい判別部３５とを備えている。
　計測時間決定部３１は、計測時間を決定する。
　振動計測部３２は、配管に設置された少なくとも２つのセンサを用いて、計測時間について、振動波形を計測する。
　相互相関関数算出部３３は、計測された振動波形の相互相関関数を算出する。
　ピーク検出部３４は、計測時間よりも短い積算時間ごとに、相互相関関数のピークを検出する。
　漏えい判別部３５は、計測時間について、積算時間ごとにピークが検出された場合、漏えいがあると判別する。
　（本実施形態の効果）
　本実施形態の構成によれば、配管に設置された少なくとも２つのセンサを用いて、計測時間について、振動波形を計測する。計測された振動波形の相互相関関数のピークを検出する。計測時間について、相互相関関数のピークが繰り返す場合、漏えいがあると判別する。したがって、１回の計測のみで、効率的に漏えいの有無を判別することができる。
　〔実施形態４〕
　（ハードウェア構成について）
　本開示の各実施形態において、各装置の各構成要素は、機能単位のブロックを示している。各装置の各構成要素の一部又は全部は、例えば図１６に示すような情報処理装置９００とプログラムとの任意の組み合わせにより実現される。図１６は、各装置の各構成要素を実現する情報処理装置９００のハードウェア構成の一例を示すブロック図である。
　図１６に示すように、情報処理装置９００は、一例として、以下のような構成を含む。
　　・ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）９０１
　　・ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）９０２
　　・ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）９０３
　　・ＲＡＭ９０３にロードされるプログラム９０４
　　・プログラム９０４を格納する記憶装置９０５
　　・記録媒体９０６の読み書きを行うドライブ装置９０７
　　・通信ネットワーク９０９と接続する通信インターフェース９０８
　　・データの入出力を行う入出力インターフェース９１０
　　・各構成要素を接続するバス９１１
　各実施形態における各装置の各構成要素は、これらの機能を実現するプログラム９０４をＣＰＵ９０１が取得して実行することで実現される。各装置の各構成要素の機能を実現するプログラム９０４は、例えば、予め記憶装置９０５やＲＯＭ９０２に格納されており、必要に応じてＣＰＵ９０１がＲＡＭ９０３にロードして実行される。なお、プログラム９０４は、通信ネットワーク９０９を介してＣＰＵ９０１に供給されてもよいし、予め記録媒体９０６に格納されており、ドライブ装置９０７が当該プログラムを読み出してＣＰＵ９０１に供給してもよい。
　以上、実施形態（及び実施例）を参照して本発明を説明したが、本発明は上記実施形態（及び実施例）に限定されるものではない。上記実施形態（及び実施例）の構成や詳細には、本発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
　この出願は、２０１８年６月１５日に出願された日本出願特願２０１８－１１４２４６を基礎とする優先権を主張し、その開示の全てをここに取り込む。
　　４　　データ処理器
　４ａ　　データ処理器
　４ｂ　　漏えい調査装置
　１２　　ヒストグラム関数算出部
　１３　　計測時間決定部
　１４　　振動計測部
　１５　　相互相関関数算出部
　１６　　積算時間決定部
　１７　　ピーク検出部
　１９　　漏えい判別部
　２０　　表示部
　２１　　判別正解率関数算出部
　２２　　計測時間決定部
　２３　　調査効率関数算出部
　２４　　経済効率関数算出部
　１００　漏えい調査システム

Claims (7)

  1. 　計測時間を決定する計測時間決定手段と、
    　配管に設置された少なくとも２つのセンサを用いて、前記計測時間について、振動波形を計測する振動計測手段と、
    　計測された前記振動波形の相互相関関数を算出する相互相関関数算出手段と、
    　前記計測時間中に２回以上、前記相互相関関数のピークを検出するピーク検出手段と、
    　前記計測時間について、前記ピークが繰り返し検出された場合、漏えいがあると判別する漏えい判別手段と
    を備えた漏えい調査装置。
  2. 　水の使用時間を示す情報から、１回の水の使用時間と頻度との関係を示すヒストグラム関数を算出するヒストグラム関数算出手段をさらに備え、
    　前記計測時間決定手段は、前記ヒストグラム関数に基づいて、音源が水の使用による外乱である蓋然性を表す確率が、予め決めた所定値以下になるような時間を、前記計測時間として決定する
    ことを特徴とする請求項１に記載の漏えい調査装置。
  3. 　前記ピーク検出手段による、前記ピークの検出結果に基づいて、前記計測時間について前記ピークが繰り返すかどうかを判別する連続性判別手段をさらに備え、
    　前記計測時間について、前記ピークが繰り返すと前記連続性判別手段が判別した場合、前記漏えい判別手段は、前記振動波形の音源は漏えいであると判別し、
    　前記計測時間について、前記ピークが繰り返さないと前記連続性判別手段が判別した場合、前記漏えい判別手段は、前記振動波形の音源は漏えいではないと判別する
    ことを特徴とする請求項１または２に記載の漏えい調査装置。
  4. 　前記振動波形のＳ／Ｎ比に基づいて、前記振動波形の相互相関関数の値を累積する積算時間を決定する積算時間決定手段をさらに備え、
    　前記ピーク検出手段は、前記積算時間にわたって積算した前記相互相関関数のピークを検出する
    ことを特徴とする請求項１から３のいずれか１項に記載の漏えい調査装置。
  5. 　水の使用時間の統計データを用いて、前記振動波形の音源の種類を正しく判別できる確率と時間との関係を示す判別正解率関数を算出する判別正解率関数算出手段と、
    　調査計画の情報を用いて、漏えい調査の効率と時間との関係を示す調査効率関数を算出する算出手段と、
    　前記判別正解率関数および前記漏えい調査効率関数から、時間に対する漏えい調査の経済効率を示す経済効率関数を算出する経済効率関数算出手段と、をさらに備え、
    　前記計測時間決定手段は、前記経済効率関数に基づいて、漏えい調査の経済効率が最大になるように、前記計測時間を決定する
    　ことを特徴とする請求項１から４のいずれか１項に記載の漏えい調査装置。
  6. 　計測時間を決定し、
    　配管に設置された少なくとも２つのセンサを用いて、前記計測時間について、振動波形を計測し、
    　計測された前記振動波形の相互相関関数を算出し、
    　前記計測時間中に２回以上、前記相互相関関数のピークを検出し、
    　前記計測時間について、前記ピークが繰り返し検出された場合、漏えいがあると判別する
    ことを含む漏えい調査方法。
  7. 　計測時間を決定することと、
    　配管に設置された少なくとも２つのセンサを用いて、前記計測時間について、振動波形を計測し、
    　計測された前記振動波形の相互相関関数を算出することと、
    　前記計測時間中に２回以上、前記相互相関関数のピークを検出することと、
    　前記計測時間について、前記ピークが繰り返し検出された場合、漏えいがあると判別することと
    をコンピュータに実行させるためのプログラムを格納した、一時的でない記録媒体。
PCT/JP2019/023546 2018-06-15 2019-06-13 漏えい調査装置、漏えい調査方法、および記録媒体 Ceased WO2019240231A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP19818550.6A EP3809107A4 (en) 2018-06-15 2019-06-13 LEAK INSPECTION DEVICE, LEAK INSPECTION PROCEDURE AND RECORDING MEDIA
US17/056,463 US11402290B2 (en) 2018-06-15 2019-06-13 Leakage inspection device, leakage inspection method, and storage medium
JP2020525661A JP7052869B2 (ja) 2018-06-15 2019-06-13 漏えい調査装置、漏えい調査方法、およびプログラム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2018-114246 2018-06-15
JP2018114246 2018-06-15

Publications (1)

Publication Number Publication Date
WO2019240231A1 true WO2019240231A1 (ja) 2019-12-19

Family

ID=68842310

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/023546 Ceased WO2019240231A1 (ja) 2018-06-15 2019-06-13 漏えい調査装置、漏えい調査方法、および記録媒体

Country Status (4)

Country Link
US (1) US11402290B2 (ja)
EP (1) EP3809107A4 (ja)
JP (1) JP7052869B2 (ja)
WO (1) WO2019240231A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3892976B1 (fr) * 2020-04-09 2023-03-08 Sagemcom Energy & Telecom SAS Procede de detection et de localisation d'une fuite de fluide
CN116658835A (zh) * 2023-05-22 2023-08-29 山东普赛通信科技股份有限公司 基于主动式声波测量的管网泄漏检测定位方法及系统

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111414955B (zh) * 2020-03-17 2023-08-25 昆仑数智科技有限责任公司 石油钻井井漏溢流智能检测方法、装置以及电子设备
CN113606504B (zh) * 2021-08-20 2022-07-22 合肥工业大学 供液管道泄漏量声学检测方法
CN114166214B (zh) * 2021-11-12 2023-06-06 海南大学 基于sta/lta算法的管道连接器检测方法和系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11201858A (ja) 1998-01-19 1999-07-30 Osaka Gas Co Ltd 相関法を用いる管路系の振動測定方法
US20080314122A1 (en) * 2007-06-21 2008-12-25 National Research Council Of Canada Monitoring of leakage in wastewater force mains and other pipes carrying fluid under pressure
KR20130064403A (ko) * 2011-12-08 2013-06-18 한국원자력연구원 매설배관 누설 탐지용 상호상관함수기법의 정확도 향상을 위한 기계 잡음 제거 방법
WO2014050923A1 (ja) 2012-09-28 2014-04-03 日本電気株式会社 漏洩判定方法、漏洩判定システム、及びプログラム
WO2015146082A1 (ja) * 2014-03-26 2015-10-01 日本電気株式会社 漏洩検知装置、漏洩検知方法、およびプログラムを格納する記録媒体
JP2016024023A (ja) * 2014-07-18 2016-02-08 積水化学工業株式会社 異常音発生位置の特定方法
WO2017094846A1 (ja) * 2015-12-03 2017-06-08 日本電気株式会社 装置、方法及び記録媒体
JP2018044799A (ja) * 2016-09-12 2018-03-22 積水化学工業株式会社 分析データの作成方法、ノイズフィルター、漏水有無判定方法、漏水場所の特定方法、振動データの計測方法
JP2018114246A (ja) 2017-01-20 2018-07-26 花王株式会社 肌状態評価方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006003311A (ja) * 2004-06-21 2006-01-05 Mitsubishi Electric Corp 異常箇所検出装置
GB0914463D0 (en) * 2009-08-19 2009-09-30 Sev Trent Water Ltd Leak detector
JP6413741B2 (ja) * 2014-12-16 2018-10-31 日本電気株式会社 振動発生源推定装置、方法およびプログラム

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11201858A (ja) 1998-01-19 1999-07-30 Osaka Gas Co Ltd 相関法を用いる管路系の振動測定方法
US20080314122A1 (en) * 2007-06-21 2008-12-25 National Research Council Of Canada Monitoring of leakage in wastewater force mains and other pipes carrying fluid under pressure
KR20130064403A (ko) * 2011-12-08 2013-06-18 한국원자력연구원 매설배관 누설 탐지용 상호상관함수기법의 정확도 향상을 위한 기계 잡음 제거 방법
WO2014050923A1 (ja) 2012-09-28 2014-04-03 日本電気株式会社 漏洩判定方法、漏洩判定システム、及びプログラム
WO2015146082A1 (ja) * 2014-03-26 2015-10-01 日本電気株式会社 漏洩検知装置、漏洩検知方法、およびプログラムを格納する記録媒体
JP2016024023A (ja) * 2014-07-18 2016-02-08 積水化学工業株式会社 異常音発生位置の特定方法
WO2017094846A1 (ja) * 2015-12-03 2017-06-08 日本電気株式会社 装置、方法及び記録媒体
JP2018044799A (ja) * 2016-09-12 2018-03-22 積水化学工業株式会社 分析データの作成方法、ノイズフィルター、漏水有無判定方法、漏水場所の特定方法、振動データの計測方法
JP2018114246A (ja) 2017-01-20 2018-07-26 花王株式会社 肌状態評価方法

Non-Patent Citations (1)

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

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3892976B1 (fr) * 2020-04-09 2023-03-08 Sagemcom Energy & Telecom SAS Procede de detection et de localisation d'une fuite de fluide
CN116658835A (zh) * 2023-05-22 2023-08-29 山东普赛通信科技股份有限公司 基于主动式声波测量的管网泄漏检测定位方法及系统

Also Published As

Publication number Publication date
JPWO2019240231A1 (ja) 2021-06-03
JP7052869B2 (ja) 2022-04-12
EP3809107A1 (en) 2021-04-21
US20210116323A1 (en) 2021-04-22
EP3809107A4 (en) 2021-08-04
US11402290B2 (en) 2022-08-02

Similar Documents

Publication Publication Date Title
JP7052869B2 (ja) 漏えい調査装置、漏えい調査方法、およびプログラム
US7729876B2 (en) Diagnostic device for use in process control system
EP1007931B1 (en) Improved method for detecting leaks in pipelines
US10704982B2 (en) Sensor recording analysis apparatus and method
CN100573381C (zh) 用于过程控制系统的诊断设备
KR101876730B1 (ko) 상수도 누수탐지 감시 시스템
US10184611B2 (en) Detecting fluid properties of a multiphase flow in a condensate drain
US7357034B1 (en) Dynamic transient pressure detection system
GB2549420A (en) Pipe inspection system, pipe inspection device, pipe inspection method, and recording medium
EP3594653A1 (en) Diagnosis cost output device, diagnosis cost output method, and computer-readable recording medium
US20190137044A1 (en) Leakage position analyzing system, leakage position analyzing method, leakage position analyzing device, and computer-readable recording medium
WO2024177834A1 (en) Pipe material detection using acoustical wave propagation
CN101501469A (zh) 用于诊断在压力测量变换器中脉冲导管的阻塞的方法及压力测量变换器
US11029181B2 (en) Vortex flowmeter with flow instability detection
EP3591364A1 (en) Measurement time specification device, sensing device, measurement time specification method, and computer-readable recording medium
US20250085184A1 (en) Probability-based pipeline leak monitoring tool
JP2019109194A (ja) 流量計測装置
JP2008203000A (ja) 詰まり診断装置及び詰まり診断方法
JP7095743B2 (ja) 流体漏洩診断装置、流体漏洩診断システム、流体漏洩診断方法、及び、流体漏洩診断プログラム
JP2024129574A (ja) データ処理装置、データ処理方法、およびデータ処理プログラム
Jin Investigation on Parameters Affecting the Performance of Negative Pressure Wave Leak Detection Systems
EA048867B1 (ru) Способ и интегрированная система непрерывного контроля трубопровода с текучей средой под давлением
Jara-Arriagada et al. Field validation of a pressure transient source localization method in a water distribution network
JP2023125843A (ja) 漏水検知装置、漏水検知方法、漏水監視システム、およびプログラム
US20210164859A1 (en) Analyzing device, analysis method, and storage medium

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: 19818550

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020525661

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2019818550

Country of ref document: EP