Propagation time measuring device

Abstract

To provide a technique that can accurately determine a propagation time.SOLUTION: There is provided a device that determines the propagation time of an acoustic signal through mutual correlation analysis between a transmission signal and a received signal, and the device determines the propagation time through mutual correlation analysis between the transmission signal and the received signal before removal of reverberation, removes signals subsequent to a timing based on the determined propagation time from the received signal as reverberation, and determines the propagation time again through mutual correlation analysis between the transmission signal and the received signal after the removal of reverberation.SELECTED DRAWING: Figure 3

Description

[Technical field]

[0001] The present invention relates to a technique for measuring the flight time of an acoustic signal. [Background technology]

[0002] Devices that can measure the flow rate and velocity of the fluid flowing through a pipe non-destructively based on the propagation time of an acoustic signal that propagates through the pipe using a sensor attached to the outside of the pipe have been put to practical use. This type of device generally uses ultrasonic waves as the acoustic signal, and is called an "ultrasonic flowmeter" or "ultrasonic flowmeter."

[0003] For example, Patent Document 1 discloses an apparatus that uses a pair of ultrasonic transducers arranged on the upstream and downstream sides of a pipe to determine the flow rate of a fluid based on the propagation time difference between ultrasonic waves propagated in the forward direction of the fluid flow and ultrasonic waves propagated in the reverse direction. The apparatus of Patent Document 1 calculates the propagation time difference by calculating the cross-correlation between the received signal of the upstream ultrasonic transducer and the received signal of the downstream ultrasonic transducer. [Prior art documents] [Patent documents]

[0004] [Patent Document 1] JP 2017-187310 A Summary of the Invention [Problem to be solved by the invention]

[0005] However, since a transducer that receives an acoustic signal such as an ultrasonic wave cannot immediately stop when the received acoustic signal disappears, the electrical signal (received signal) output from the transducer contains reverberation. If the received signal contains reverberation, the reverberation becomes noise, making it impossible to accurately calculate the propagation time, and therefore the flow velocity and flow rate of the fluid flowing through the pipe.

[0006] The present invention has been made in consideration of the above-mentioned circumstances, and has an object to provide a technique capable of determining the propagation time with high accuracy. [Means for solving the problem]

[0007] In order to achieve the above object, the present invention employs the following configuration.

[0008] A first aspect of the present invention provides a propagation time measuring device comprising: a plurality of transducers arranged at different positions relative to a pipe through which a fluid flows, the plurality of transducers including at least a first transducer that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer that receives the acoustic signal transmitted from the first transducer and propagated through the fluid in the pipe and converts the acoustic signal into a reception signal as an electric signal; and a signal processing unit that removes reverberation from the reception signal and determines a propagation time of the acoustic signal from the first transducer to the second transducer by cross-correlation analysis between the transmission signal and the reception signal after the reverberation has been removed, the signal processing unit determines the propagation time of the acoustic signal from the first transducer to the second transducer by cross-correlation analysis between the transmission signal and the reception signal before the reverberation has been removed, removes signals after a timing based on the determined propagation time as the reverberation, and determines again the propagation time of the acoustic signal from the first transducer to the second transducer by cross-correlation analysis between the transmission signal and the reception signal after the reverberation has been removed.

[0009] Reverberation in the received signal affects the desired propagation time. The timing at which reverberation occurs can be determined with high accuracy based on the propagation time determined using the received signal before reverberation is removed. In the above-described configuration, the signal after the timing thus determined is removed as reverberation from the received signal, and the propagation time is determined again using the received signal after reverberation has been removed. In this manner, reverberation can be removed with high accuracy, and the propagation time can be determined with high accuracy. Furthermore, since unnecessary signal values ​​(reverberation signal values) are not used in calculating the propagation time, it is possible to shorten the processing time required to calculate the propagation time and reduce the power consumption required to calculate the propagation time.

[0010] It is considered that the receiving transducer receives an acoustic signal having the same time length as the time length of the transmission signal from the time when the propagation time has elapsed, and the signal of the reception signal after the reception of the acoustic signal is considered to be reverberation. Therefore, in the first aspect, the timing may be the time when the total time of the propagation time calculated by cross-correlation analysis between the transmission signal and the reception signal before the reverberation is removed and the time length of the transmission signal has elapsed from the time when the transmission signal is input to the first transducer. In this way, it is possible to more accurately remove reverberation and more accurately calculate the propagation time.

[0011] A second aspect of the present invention provides a propagation time measuring device comprising: a plurality of transducers arranged at different positions relative to a pipe through which a fluid flows, the plurality of transducers including at least a first transducer that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer that receives the acoustic signal transmitted from the first transducer and propagated through the fluid in the pipe and converts it into a reception signal as an electric signal; and a signal processing unit that removes reverberation from the reception signal and determines a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmission signal and the reception signal after the reverberation has been removed, wherein the signal processing unit approximately obtains an envelope of a waveform of the reception signal from the reception signal before the reverberation has been removed, determines a threshold based on the envelope, and removes, as the reverberation, a signal after the timing at which the signal value converges below the determined threshold.

[0012] The waveform of the received signal has a relatively large envelope. The reverberation signal value is relatively small. Therefore, the timing at which reverberation occurs can be accurately determined based on the envelope. In the above-mentioned configuration, the envelope is approximately determined from the received signal, and the signal after the timing at which the signal value converges below a threshold based on the envelope is removed as reverberation from the received signal. Then, the propagation time is determined using the received signal after the reverberation has been removed. In this way, the reverberation can be accurately removed, and the propagation time can be accurately determined. Furthermore, since unnecessary signal values ​​(reverberation signal values) are not used in the calculation of the propagation time, the processing time required for the calculation of the propagation time can be shortened, and the power consumption required for the calculation of the propagation time can be reduced.

[0013] In a second aspect, the signal processor may determine, as the threshold, a value obtained by reducing a peak of the envelope by a predetermined factor. In the second aspect, the signal processor may determine, as the threshold, a predetermined number of discrete values ​​counted from the larger one of a plurality of discrete values ​​constituting the envelope.

[0014] A third aspect of the present invention provides a system including a plurality of transducers arranged at different positions relative to a pipe through which a fluid flows, the plurality of transducers including at least a first transducer that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer that receives the acoustic signal transmitted from the first transducer and propagated through the fluid in the pipe and converts the acoustic signal into a reception signal as an electric signal, and a signal processing unit that removes reverberation from the reception signal, and obtains a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmission signal and the reception signal after the reverberation has been removed, the signal processing unit determining a propagation time of the acoustic signal from the first transducer to the second transducer after a timing at which a signal value converges to or below a predetermined number of peaks counted from the larger side among a plurality of peaks shown in the reception signal before the reverberation is removed. The present invention provides a propagation time measuring device that is characterized by removing the above-mentioned signal as the reverberation.

[0015] The received signal shows multiple peaks that gradually increase and then decrease. The reverberation signal value is relatively small. Therefore, the timing at which reverberation occurs can be accurately determined based on the multiple peaks. In the above-mentioned configuration, a predetermined number of peaks from the larger side among the multiple peaks is set as a threshold, and the signal from the received signal after the signal value converges to or below the threshold is removed as reverberation, and the propagation time is determined using the received signal after the reverberation has been removed. In this way, the reverberation can be accurately removed, and the propagation time can be accurately determined. Furthermore, since unnecessary signal values ​​(reverberation signal values) are not used in the calculation of the propagation time, the processing time required for the calculation of the propagation time can be shortened and the power consumption required for the calculation of the propagation time can be reduced.

[0016] In each of the first, second, and third aspects, the signal processing unit may further include a storage unit that stores information about the timing in advance, and the signal processing unit may remove the reverberation based on the information stored in advance by the storage unit. The information about the timing may be information indicating the timing itself, or may not be information about the timing itself. For example, the information about the timing may be information indicating a propagation time calculated using a received signal before reverberation is removed, or may be information indicating a threshold based on an envelope or a plurality of peaks.

[0017] In each of the first, second and third aspects, the signal processing unit may acquire the information about the timing using a plurality of reception signals obtained by performing transmission from the first transducer and reception by the second transducer a plurality of times. In this way, the timing at which reverberation occurs can be determined with higher accuracy. As a result, reverberation can be removed with higher accuracy, and the propagation time can be determined with higher accuracy.

[0018] The first and second oscillators may be arranged on each of the first, second and third side faces to face each other with the piping in between. The first and second oscillators may be arranged on each of the first, second and third side faces to be at different positions in a longitudinal direction of the piping.

[0019] In each of the first, second, and third aspects, the transmission signal may be input to the second transducer, and a switching unit may be configured to switch the first transducer to output the reception signal after receiving the acoustic signal transmitted from the second transducer, and the signal processing unit may further remove reverberation from the reception signal output from the first transducer, and calculate a propagation time of the acoustic signal from the second transducer to the first transducer by performing a cross-correlation analysis between the transmission signal input to the second transducer and the reception signal output from the first transducer after the reverberation has been removed. In this way, the propagation time when the acoustic signal propagates from the upstream side to the downstream side and the propagation time when the acoustic signal propagates from the downstream side to the upstream side can be accurately calculated for the same propagation path.

[0020] The signal processing unit may determine a flow velocity and / or a flow rate of the fluid in the pipe based on a difference between a propagation time of the acoustic signal from the first transducer to the second transducer and a propagation time of the acoustic signal from the second transducer to the first transducer, thereby making it possible to obtain information about the fluid in the pipe with high accuracy.

[0021] The present invention may be understood as a propagation time measuring device having at least a part of the above configuration, or as a flow velocity measuring device, a flow rate measuring device, a flow meter, a flow rate sensor, etc. The present invention may also be understood as a propagation time measuring method, a flow velocity measuring method, or a flow rate measuring method including at least a part of the above processing, or as a program for realizing such a method, or a recording medium on which such a program is non-temporarily recorded. Each of the above may be combined with each other as much as possible to constitute the present invention. Effect of the Invention

[0022] According to the present invention, the propagation time can be determined with high accuracy. [Brief description of the drawings]

[0023] [Figure 1] FIG. 1 is a diagram illustrating a schematic configuration of a propagation time measuring device. [Diagram 2] FIG. 2 is a cross-sectional view showing an example of installation of a transducer in a pipe. [Diagram 3] FIG. 3 is a flowchart showing the flow of the measurement operation of the propagation time measuring device in the first embodiment. [Figure 4] FIG. 4A is a diagram showing an example of a transmitted signal, FIG. 4B is a diagram showing an example of an acoustic signal based on the transmitted signal of FIG. 4A, FIG. 4C is a diagram showing an example of a received signal, and FIG. 4D is a diagram showing an example of a cross-correlation function. [Diagram 5] FIG. 5A is a diagram showing an example of a cross-correlation function, and FIG. 5B is a diagram showing an example of a Hilbert transform. [Figure 6] FIG. 6 is a diagram illustrating dereverberation in the first embodiment. [Figure 7] FIG. 7 is a flowchart showing the flow of the measurement operation of the propagation time measuring device in the second embodiment. [Figure 8] FIG. 8 is a diagram illustrating dereverberation in the second embodiment. [Figure 9] FIG. 9 is a flowchart showing the flow of the measurement operation of the propagation time measuring device in the third embodiment. [Figure 10] FIG. 10 is a diagram illustrating dereverberation in the third embodiment. [Figure 11] FIG. 11 is a flowchart showing the flow of the measurement operation of the propagation time measuring device in the fourth embodiment. [Figure 12] FIG. 12 is a diagram illustrating dereverberation in the fourth embodiment. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0024] <Application Examples> An application example of the propagation time measuring device will be described with reference to FIG.

[0025] The propagation time measuring device 1 includes two or more transducers 101, and receives an acoustic signal transmitted from one of the transducers (e.g., 101a) by another transducer (e.g., 101b), and calculates the time (propagation time) required for the acoustic signal to propagate along a path between the two transducers. Since the transducers 101 are disposed at different positions relative to the pipe 120, the acoustic signal propagating between the two transducers 101 passes through (crosses) the inside of the pipe 120. Therefore, the propagation time of the acoustic signal is not constant, but changes depending on the state of the fluid 121 flowing inside the pipe 120 (e.g., flow velocity, flow rate, presence of air bubbles or foreign matter, etc.). Therefore, by utilizing the propagation time measured by the propagation time measuring device 1, the state of the fluid 121 inside the pipe 120 can be measured nondestructively.

[0026] The fluid 121 may be a liquid or a gas as long as it is a material capable of propagating an acoustic signal. The acoustic signal is typically an ultrasonic wave, but may also include sound waves in the audible range.

[0027] The propagation time measuring device 1 uses cross-correlation analysis to calculate the propagation time. For example, if the driving electrical signal given to the transmitting transducer 101 is called the "transmitted signal" and the electrical signal output from the receiving transducer 101 is called the "received signal," the propagation time measuring device 1 calculates the cross-correlation function between the transmitted signal and the received signal. Then, based on the maximum peak position of the cross-correlation function, the propagation time measuring device 1 determines the lag (time delay) of the received signal relative to the transmitted signal. This lag corresponds to the propagation time of the acoustic signal from the transmitting transducer 101 to the receiving transducer 101. do.

[0028] When the signal waveform of the transmission signal is sufficiently preserved in the received signal, a clear peak appears in the cross-correlation function, and the lag between the two signals (i.e., the propagation time) can be obtained with high accuracy. However, since the receiving transducer 101 cannot immediately stop even when the acoustic signal to be received disappears, the received signal contains reverberation. If the received signal contains reverberation, the reverberation becomes noise, and the cross-correlation function causes a peak to appear at a time position different from the actual propagation time, and the propagation time cannot be obtained with high accuracy. As a result, the flow velocity and flow rate of the fluid 121 flowing through the pipe 120 cannot be obtained with high accuracy. Recently, there has been an increasing need for a device capable of measuring a minute flow rate of fluid, and in order to realize this, a performance capable of measuring the propagation time with an accuracy of, for example, nanosecond order to picosecond order is required.

[0029] Therefore, the propagation time measuring device 1 removes reverberation from the received signal, and calculates the propagation time using the received signal after the reverberation has been removed. Specifically, the propagation time measuring device 1 removes reverberation by any of the following methods 1 to 3. Note that the removal of reverberation may be a process that excludes the signal value of the reverberation from the cross-correlation analysis (calculation of the cross-correlation function), and may be, for example, a process that sets (cuts) the signal value corresponding to the reverberation to zero, or a process that makes the signal value of the reverberation distinguishable from other signal values.

[0030] (Method 1) In method 1, the propagation time measurement device 1 calculates the propagation time by performing cross-correlation analysis between the transmitted signal and the received signal before reverberation is removed, and removes the signal after the timing based on the calculated propagation time as reverberation.The propagation time measurement device 1 then calculates the propagation time again by performing cross-correlation analysis between the transmitted signal and the received signal after reverberation has been removed.

[0031] Reverberation contained in the received signal affects the calculated propagation time. Therefore, the timing at which reverberation occurs can be calculated with high accuracy based on the propagation time calculated using the received signal before reverberation is removed. In method 1, the signal from the received signal after the timing calculated in this way is removed as reverberation, and the propagation time is calculated again using the received signal after reverberation has been removed. In this way, reverberation can be accurately removed, the correct peak (the peak at the position corresponding to the actual propagation time) can be accurately detected from the cross-correlation function, and the propagation time can be accurately calculated.

[0032] (Method 2) In method 2, the propagation time measuring device 1 approximately calculates the envelope of the waveform of the received signal from the received signal before reverberation is removed, determines a threshold based on the calculated envelope, and removes as reverberation the signal from the point in time when the signal value converges below the determined threshold.

[0033] The waveform of the received signal has a relatively large envelope. The reverberation signal value is relatively small. Therefore, the timing at which reverberation occurs can be determined accurately based on the envelope. In method 2, the envelope is approximately determined from the received signal, and the signal from the received signal after the signal value converges below a threshold based on the envelope is removed as reverberation. The flight time is then determined using the received signal after reverberation has been removed. In this way, reverberation can be accurately removed, the correct peak can be accurately detected from the cross-correlation function, and the flight time can be accurately determined.

[0034] (Method 3) In the third method, the propagation time measuring device 1 measures the propagation time of the signal after the timing at which the signal value converges to or below a predetermined number of peaks from the larger side among the multiple peaks shown in the received signal before the reverberation is removed. is removed as reverberation.

[0035] The received signal shows multiple peaks that gradually increase and then decrease. The reverberation signal value is relatively small. Therefore, the timing at which reverberation occurs can be determined accurately based on the multiple peaks. In method 3, a predetermined number of peaks from the larger side among the multiple peaks is set as a threshold, and the signal from the received signal after the signal value converges to below the threshold is removed as reverberation, and the flight time is determined using the received signal after reverberation has been removed. In this way, reverberation can be accurately removed, the correct peak can be accurately detected from the cross-correlation function, and the flight time can be accurately determined.

[0036] Furthermore, in any of the methods 1 to 3, unnecessary signal values ​​(reverberation signal values) are not used in the calculation of the propagation time (specifically, the calculation of the cross-correlation function), so that it is possible to shorten the processing time required for the calculation of the propagation time and reduce the power consumption required for the calculation of the propagation time.

[0037] First Embodiment (Device configuration) A specific configuration of the propagation time measuring device 1 will be described with reference to Figures 1 and 2. Figure 1 is a block diagram that shows a schematic configuration of the propagation time measuring device 1, and Figure 2 is a cross-sectional view showing an example of installing a transducer in a pipe. The propagation time measuring device 1 of this embodiment is a device for non-destructively measuring the flow velocity and flow rate of a fluid 121 flowing in a pipe 120, and is also called an ultrasonic flowmeter or ultrasonic flow sensor.

[0038] The propagation time measuring device 1 has a device body 100 and a plurality of transducers 101. The device body 100 and each transducer 101 are connected by a cable. In this embodiment, two transducers 101 are provided: a first transducer 101a arranged on the upstream side in the longitudinal direction of the pipe 120, and a second transducer 101b arranged on the downstream side of the first transducer 101a. Hereinafter, when it is necessary to distinguish between the two transducers, they will be referred to as "first transducer 101a" and "second transducer 101b", and when the description is common to both, they will simply be referred to as "transducer 101". Note that the number of transducers 101 is not limited to two, and three or more transducers 101 may be provided.

[0039] The vibrator 101 is a device that converts between an electric signal and an acoustic signal, and is also called a transducer. For example, a piezoelectric element that converts between voltage and force by the piezoelectric effect can be used as the vibrator 101. As shown in FIG. 2, each vibrator 101 is embedded in a clamp 30 made of resin. When the clamp 30 clamps the pipe 120, the two vibrators 101a and 101b face each other across the pipe 120, and each vibrator 101a and 101b is installed so that a straight line connecting the two vibrators 101a and 101b forms a predetermined angle θ with respect to the axis of the pipe 120. By adopting such a clamp structure, the vibrator 101 can be easily attached to an appropriate position on the existing pipe 120 (without modifying the pipe 120). In order to make the clamp 30 and the pipe 120 tightly contact each other and to match the acoustic impedance, grease, gel, or the like may be applied between the clamp 30 and the pipe 120. The angle θ is called the propagation angle of the acoustic signal. The propagation angle θ can be any value, but when using the transit time difference method described later, it is recommended to set the angle θ in the range of 0 degrees < θ < 90 degrees, preferably 20 degrees < θ < 60 degrees.

[0040] The device main body 100 mainly comprises a control circuit 102, a D / A converter 103, an A / D converter 104, a switch 105, and an output unit 106. The control circuit 102 is a circuit that controls each part of the propagation time measuring device 1, processes signals, performs calculations, etc. The D / A converter 103 is a circuit that performs D / A conversion and signal amplification based on a transmission signal (digital data) input from the control circuit 102, and outputs a transmission signal (analog signal) of a predetermined voltage to the transducer 101. The A / D converter 104 is a circuit that performs A / D conversion on the reception signal (analog signal) input from the transducer 101 at a predetermined sampling period and outputs the reception signal (digital data) to the control circuit 102. The switch 105 is a switch that switches the connection relationship between the D / A converter 103 and the A / D converter 104 and the first transducer 101a and the second transducer 101b. The transducer 101 connected to the D / A converter 103 is the transmitting side, and the transducer 101 connected to the A / D converter 104 is the receiving side. The output unit 106 is a device that outputs information such as the results of signal processing and arithmetic processing by the control circuit 102, and is, for example, a display. Note that the device main body 100 may be provided with an input unit (for example, a button, a touch panel, etc.) for a user to operate, or a communication circuit (for example, a WiFi module, etc.) that transmits information to an external device (for example, an external computer or server, etc.).

[0041] As shown in FIG. 1, the control circuit 102 includes a transmission signal generating unit 110, a signal processing unit 111, and a storage unit 112. The transmission signal generating unit 110 has a function of generating data of a transmission signal used in measurement and outputting the data to the D / A converter 103. The signal processing unit 111 has a function of calculating the propagation time of an acoustic signal based on a transmission signal and a reception signal, and further calculating the flow velocity and / or flow rate of a fluid from the propagation time. The signal processing unit 111 also has a function of removing reverberation from a reception signal. In this embodiment, the signal processing unit 111 removes reverberation by the above-mentioned method 1. The storage unit 112 stores waveform data that defines the waveform of a transmission signal. The transmission signal generating unit 110 reads the waveform data from the storage unit 112 and generates data of a transmission signal.

[0042] The control circuit 102 may be configured by a computer having, for example, a CPU (processor), a RAM, a non-volatile storage device (for example, a ROM, a flash memory, a hard disk, etc.), an I / O, etc. In this case, the CPU expands a program stored in the storage device into the RAM and executes the program, thereby providing the functions of the transmission signal generation unit 110 and the signal processing unit 111. The form of the computer does not matter. For example, it may be a personal computer, an embedded computer, a smartphone, a tablet terminal, etc. Alternatively, all or a part of the functions provided by the control circuit 102 may be configured by a circuit such as an ASIC or an FPGA. Alternatively, the control circuit 102 may perform the processing described later in cooperation with another computer by utilizing distributed computing or cloud computing technology.

[0043] The material, size, and shape of the pipe 120 are not important. For example, metal pipes or resin pipes may be used. The size of the pipe 120 may be a standard size defined by JIS or ANSI, or may be a unique size. The method of the present embodiment has an advantage that it can measure a minute flow rate with high accuracy, and is particularly suitable for measuring small pipes such as 1 / 8 inch pipes (outer diameter: 3.18 mm, inner diameter: 1.59 mm), 1 / 4 inch pipes (outer diameter: 6.35 mm, inner diameter: 3.97 mm), and 1 / 2 inch pipes (outer diameter: 12.70 mm, inner diameter: 9.53 mm). The pipes are not limited to straight pipes, and may be bent pipes or curved pipes, and the cross-sectional shape of the pipes is also arbitrary.

[0044] (Measurement operation) The flow of the measurement operation of the propagation time measuring device 1 will be described with reference to the flowchart of FIG.

[0045] In step S100, the transmission signal generating unit 110 of the control circuit 102 reads from the storage unit 112 waveform data of the transmission signal.

[0046] In step S101, the control circuit 102 controls the switch 105 to connect the first oscillator 101a to the D / A converter 103 and connect the second oscillator 101b to the A / D converter 104. As a result, the first oscillator 101a becomes the transmitting side and the second oscillator 101b becomes the receiving side.

[0047] In step S102, the transmission signal generating unit 110 generates a transmission signal based on the waveform data read in step S100, and outputs the transmission signal to the D / A converter 103. In addition, this transmission signal is temporarily stored in a RAM (work memory) for use in a cross-correlation analysis at a later stage.

[0048] In step S103, the transmission signal that has been D / A converted and amplified by the D / A converter 103 is input to the first transducer 101a, and an acoustic signal based on the transmission signal is transmitted from the first transducer 101a. FIG. 4A is an example of a transmission signal, with the horizontal axis representing time and the vertical axis representing signal value. FIG. 4B is an example of an acoustic signal based on the transmission signal of FIG. 4A, with the horizontal axis representing time and the vertical axis representing sound pressure. Here, the reason why the amplitude of the acoustic signal is not constant even though the amplitude of the transmission signal is constant is due to the frequency characteristics of the transducer. The acoustic signal reaches the second transducer 101b via the clamp 30, the pipe 120, and the fluid 121.

[0049] In step S104, the second transducer 101b converts the received acoustic signal into a reception signal and outputs it to the A / D converter 104. FIG. 4C shows an example of a reception signal, with the horizontal axis representing time and the vertical axis representing signal value. In FIG. 4C, the vertical axis is enlarged compared to FIG. 4A. Since the acoustic signal attenuates during the propagation process, the amplitude (voltage) of the reception signal is on the order of about 1 / 100 to 1 / 1000 of the transmission signal. For example, the amplitude of the transmission signal in FIG. 4A is about 30V, whereas the amplitude of the reception signal in FIG. 4C is about 10mV. Also, as shown in FIG. 4C, the reception signal contains various noises including reverberation. The reception signal A / D converted by the A / D converter 104 is taken into the control circuit 102 and temporarily stored in a RAM (work memory).

[0050] In step S105, the signal processing unit 111 reads the transmission signal and the reception signal from the RAM, and calculates the cross-correlation function between the two signals. Fig. 4D shows an example of the cross-correlation function, with only the vicinity of the maximum peak enlarged. The horizontal axis is time (amount of time shift), and the vertical axis is the value of cross-correlation normalized so that the height of the maximum peak is 1. Note that the cross-correlation function is a known technique, so a detailed explanation will be omitted here.

[0051] In step S106, the signal processing unit 111 obtains the time position of the maximum peak (apex) in the cross-correlation function calculated in step S105 as the propagation time T1 of the acoustic signal from the first transducer 101a to the second transducer 101b. Here, the propagation time T1 is shifted from the actual propagation time due to the influence of reverberation of the received signal, etc.

[0052] Due to the constraints of digital signal processing, the cross-correlation function calculated in step S105 is composed of discrete data. Therefore, as shown in FIG. 5A, the points at which the data of the cross-correlation function is obtained (points indicated by black circles) do not necessarily coincide with the positions of the peaks (apexes). Therefore, in step S106, the signal processor 111 may estimate the position of the maximum peak after approximately determining the shape of the maximum peak and its surrounding waveform (maximum peak waveform) from the discrete data of the cross-correlation function. For example, as shown in FIG. 5B, the signal processor 111 may convert data near the maximum peak in the cross-correlation function into phase data by a Hilbert transform. Then, the converted data may be linearly approximated, and the zero crossing point (the position where the phase becomes zero) of the approximated line may be regarded as the position of the maximum peak. Alternatively, the signal processor 111 may estimate the shape of the maximum peak waveform by interpolating data near the maximum peak of the cross-correlation function by polynomial approximation, and determine the position of the apex. By performing such processing, it is possible to determine the position at which the cross-correlation is maximized, i.e., the propagation time of the acoustic signal, with a resolution higher than the sampling period of the A / D conversion.

[0053] In step S107, the signal processing unit 111 extracts the signal from the received signal read in step S105 after the timing based on the propagation time T1 calculated in step S106. The reverberation is removed from the received signal as reverberation. Here, it is considered that the receiving transducer receives an acoustic signal with the same time length as the time length of the transmission signal from the timing when the propagation time has elapsed, and the signal after the reception of the acoustic signal is completed is considered to be reverberation. For this reason, as shown in FIG. 6, the signal processing unit 111 removes, from the received signal, the signal after the timing when the sum of the propagation time T1 and the time length T0 of the transmission signal has elapsed from the timing when the transmission signal is input to the first transducer 101a. In this way, the reverberation can be removed with high accuracy.

[0054] In step S108, the signal processing unit 111 calculates a cross-correlation function between the transmission signal read in step S105 and the reception signal from which reverberation has been removed in step S107.

[0055] In step S109, the signal processing unit 111 obtains the time position of the maximum peak (apex) in the cross-correlation function calculated in step S108 as the propagation time T1' of the acoustic signal from the first transducer 101a to the second transducer 101b. Since reverberation has been removed from the received signal, the propagation time T1' that is close to the actual propagation time is obtained.

[0056] In step S110, the control circuit 102 controls the switch 105 to connect the D / A converter 103 to the second transducer 101b and connect the A / D converter 104 to the first transducer 101a. That is, the transducers on the transmitting side and the receiving side are switched. The processes in the following steps S111 to S118 are the same as those in steps S102 to S109 (however, the "first transducer 101a" is replaced with the "second transducer 101b" and the "second transducer 101b" is replaced with the "first transducer 101a"). In step S115, the propagation time T2 of the acoustic signal from the second transducer 101b to the first transducer 101a is obtained from the cross-correlation function based on the received signal before reverberation is removed. In step S118, a propagation time T2' of the acoustic signal from the second transducer 101b to the first transducer 101a is obtained from a cross-correlation function based on the received signal after reverberation has been removed.

[0057] By the above processing, the propagation time T1' of the acoustic signal from the first transducer 101a to the second transducer 101b and the propagation time T2' of the acoustic signal from the second transducer 101b to the first transducer 101a are obtained. When the fluid 121 in the pipe 120 is flowing, a time difference occurs between the propagation time T1' and the propagation time T2' according to the flow speed of the fluid 121. Therefore, the flow speed and the flow rate of the fluid 121 can be calculated using the propagation time T1' and the propagation time T2'. Since the reverberation is removed and the propagation time T1' and the propagation time T2' are obtained with high accuracy, the flow speed and the flow rate of the fluid 121 can also be obtained with high accuracy.

[0058] In step S119, the signal processing unit 111 calculates the flow velocity V of the fluid 121 by the following formula.

number

[0059] Here, V is the flow velocity of the fluid, L is the propagation path length inside the pipe, θ is the propagation angle, Tab is the propagation time T1' from the upstream transducer to the downstream transducer, Tba is the propagation time T2' from the downstream transducer to the upstream transducer, and To is the propagation time of the parts other than the fluid. Note that the propagation time To of the parts other than the fluid is, for example, the time it takes for the acoustic signal to propagate through the clamp 30 and the pipe 120, and can be obtained in advance by experiment or simulation if the specifications of the pipe 120 (inner diameter, outer diameter, material, etc.) are known.

[0060] In step S120, the signal processing unit 111 calculates the flow rate Q of the fluid by the following formula: .

number

[0061] Here, Q is the flow rate of the fluid, V is the flow velocity of the fluid, and A is the cross-sectional area inside the pipe. The cross-sectional area A is assumed to be known.

[0062] In step S121, the signal processing unit 111 outputs the processing results (for example, the propagation time, the flow velocity, the flow rate, etc.) to the output unit .

[0063] (Advantages of this embodiment) According to the configuration of the present embodiment described above, the timing at which reverberation occurs is accurately determined based on the propagation time determined using the received signal before reverberation is removed. Then, the signal after the timing thus determined is removed from the received signal as reverberation, and the propagation time is again determined using the received signal after reverberation is removed. In this way, reverberation can be accurately removed, and the propagation time can be accurately determined. Therefore, the present invention can be applied to situations requiring high accuracy, such as measurement of minute flow rates. Furthermore, since unnecessary signal values ​​(reverberation signal values) are not used in the calculation of the propagation time, the processing time required for the calculation of the propagation time can be shortened, and the power consumption required for the calculation of the propagation time can be reduced.

[0064] <Second embodiment> The propagation time measurement device 1 according to the second embodiment of the present invention removes reverberation by the above-mentioned method 2. Since the basic configuration is similar to that of the first embodiment, the following description will focus on the configuration parts that differ from the first embodiment.

[0065] 7 is a flowchart showing the flow of the measurement operation of the propagation time measurement device 1 according to the second embodiment. The same steps as in the flowchart of the first embodiment (FIG. 3) are denoted by the same step numbers.

[0066] The processes in steps S100 to S104 are the same as those in the first embodiment. Then, in step S200, the signal processing unit 111 reads the received signal from the RAM, and approximately calculates the envelope of the waveform of the received signal, as shown in Fig. 8. The envelope can be calculated using a known technique, and therefore a detailed description thereof will be omitted here.

[0067] In step S201, the signal processor 111 determines a threshold value obtained by reducing the peak of the envelope obtained in step S200 by a predetermined factor, and detects the time Ta until the signal value converges to the determined threshold value or less from the received signal read in step S200. The predetermined factor is not particularly limited, but is, for example, 1 / 2.

[0068] In step S202, the signal processing unit 111 removes, as reverberation, the signal after the time Ta detected in step S201 from the received signal read in step S200, as shown in FIG.

[0069] In step S203, the signal processing unit 111 reads the transmission signal from the RAM, and calculates a cross-correlation function between the read transmission signal and the reception signal from which the reverberation has been removed in step S202.

[0070] In step S204, the signal processing unit 111 calculates the time position of the maximum peak (apex) in the cross-correlation function calculated in step S203 from the first oscillator 101a to the second oscillator 101b. The propagation time of the acoustic signal from the receiving point to the receiving point is calculated as T1'. Since reverberation has been removed from the received signal, the propagation time T1' that is close to the actual propagation time is calculated.

[0071] The processes of steps S110 to S113 are the same as those of the first embodiment. The processes of steps S210 to S214 are the same as those of steps S200 to S204 (however, the "first transducer 101a" is replaced with the "second transducer 101b" and the "second transducer 101b" is replaced with the "first transducer 101a"). In step S211, the signal processing unit 111 determines a value obtained by reducing the peak of the envelope obtained in step S210 by a predetermined factor as a threshold value, and detects a time Tb from the reception signal read in step S210 until the signal value converges to the determined threshold value or less. In step S214, a propagation time T2' of the acoustic signal from the second transducer 101b to the first transducer 101a is obtained from a cross-correlation function based on the reception signal after the reverberation is removed. The processes of steps S119 to S121 are the same as those of the first embodiment.

[0072] According to the configuration of the present embodiment described above, the envelope of the waveform of the received signal before reverberation is removed is obtained from the received signal, and a value obtained by reducing the peak of the obtained envelope by a predetermined factor is determined as a threshold value. Then, the signal after the timing when the signal value converges below the determined threshold value is removed as reverberation from the received signal, and the propagation time is calculated using the received signal after reverberation is removed. In this way, reverberation can be accurately removed and the propagation time can be accurately calculated. Therefore, the present invention can be applied to situations requiring high accuracy, such as measurement of minute flow rates. Furthermore, since unnecessary signal values ​​(reverberation signal values) are not used in the calculation of the propagation time, the processing time required for the calculation of the propagation time can be shortened and the power consumption required for the calculation of the propagation time can be reduced.

[0073] <Third embodiment> The propagation time measurement device 1 according to the third embodiment of the present invention eliminates reverberation by the above-mentioned method 2. However, the specific method for eliminating reverberation differs from that of the second embodiment. Since the basic configuration is the same as that of the above embodiment, the following description will focus on the components that differ from the above embodiment.

[0074] 9 is a flowchart showing the flow of the measurement operation of the propagation time measurement device 1 according to the third embodiment. The same steps as in the flowchart of the second embodiment (FIG. 7) are denoted by the same step numbers.

[0075] The processes in steps S100 to S104 and step S200 are the same as those in the second embodiment. The envelope data obtained in step S200 is discrete data. The time intervals of the discrete values ​​constituting the envelope are not particularly limited, but are equal to the time intervals of the peaks indicated by the received signal, specifically, the time intervals of the peaks indicated by the transmitted signal.

[0076] In step S301, the signal processor 111 determines a predetermined number of discrete values ​​from the larger side among the multiple discrete values ​​constituting the envelope obtained in step S200 as a threshold value. Then, the signal processor 111 detects the time Ta from the received signal read in step S200 until the signal value converges to the determined threshold value or less. The predetermined number is not particularly limited, but is, for example, 10.

[0077] In step S302, the signal processing unit 111 removes, as reverberation, the signal after the time Ta detected in step S301 from the received signal read in step S200, as shown in FIG.

[0078] In step S303, the signal processing unit 111 reads the transmission signal from the RAM, and calculates a cross-correlation function between the read transmission signal and the reception signal from which reverberation has been removed in step S302.

[0079] In step S304, the signal processing unit 111 obtains the time position of the maximum peak (apex) in the cross-correlation function calculated in step S303 as the propagation time T1' of the acoustic signal from the first transducer 101a to the second transducer 101b. Since reverberation has been removed from the received signal, the propagation time T1' close to the actual propagation time is obtained.

[0080] The processes of steps S110 to S113 and step S210 are the same as those of the second embodiment. The processes of steps S311 to S314 are the same as those of steps S301 to S304 (however, the "first transducer 101a" is replaced with the "second transducer 101b" and the "second transducer 101b" is replaced with the "first transducer 101a"). In step S311, the signal processing unit 111 determines, as a threshold, a predetermined number of discrete values ​​from the larger side among the multiple discrete values ​​constituting the envelope obtained in step S210. Then, the signal processing unit 111 detects a time Tb from the reception signal read in step S210 until the signal value converges to the determined threshold or less. In step S314, a propagation time T2' of the acoustic signal from the second transducer 101b to the first transducer 101a is obtained from a cross-correlation function based on the reception signal after reverberation is removed. The processes in steps S119 to S121 are the same as those in the second embodiment.

[0081] According to the configuration of the present embodiment described above, the envelope of the waveform of the received signal before reverberation is removed is obtained from the received signal, and a predetermined number of discrete values ​​constituting the obtained envelope are determined as a threshold value. Then, the signal after the timing at which the signal value converges below the determined threshold value is removed as reverberation from the received signal, and the propagation time is calculated using the received signal after reverberation is removed. In this way, reverberation can be accurately removed, and the propagation time can be accurately calculated. Therefore, the present invention can be applied to situations requiring high accuracy, such as measurement of minute flow rates. Furthermore, since unnecessary signal values ​​(reverberation signal values) are not used in the calculation of the propagation time, the processing time required for the calculation of the propagation time can be shortened, and the power consumption required for the calculation of the propagation time can be reduced.

[0082] <Fourth embodiment> The propagation time measurement device 1 according to the fourth embodiment of the present invention removes reverberation by the above-mentioned method 3. Since the basic configuration is similar to that of the above embodiment, the following description will focus on configuration parts that differ from the above embodiment.

[0083] 11 is a flowchart showing the flow of the measurement operation of the propagation time measurement device 1 according to the fourth embodiment. The same steps as in the flowchart of the first embodiment (FIG. 3) are denoted by the same step numbers.

[0084] The processing of steps S100 to S104 is the same as that of the first embodiment. Then, in step S400, the signal processing unit 111 reads the received signal from the RAM, and among the multiple peaks shown by the received signal, determines a predetermined number of peaks from the larger side as a threshold value. Then, the signal processing unit 111 detects the time Ta from the received signal until the signal value converges to the determined threshold value or less. The predetermined number is not particularly limited, but is, for example, 10.

[0085] In step S401, the signal processing unit 111 removes, as reverberation, the signal after the time Ta detected in step S400 from the received signal read in step S400, as shown in FIG.

[0086] In step S402, the signal processing unit 111 reads the transmission signal from the RAM, and calculates the cross-correlation function between the read transmission signal and the received signal from which the reverberation has been removed in step S401.

[0087] In step S403, the signal processing unit 111 obtains the time position of the maximum peak (apex) in the cross-correlation function calculated in step S402 as the propagation time T1' of the acoustic signal from the first transducer 101a to the second transducer 101b. Since reverberation has been removed from the received signal, the propagation time T1' close to the actual propagation time is obtained.

[0088] The processing of steps S110 to S113 is the same as that of the first embodiment. The processing of steps S410 to S413 is the same as that of steps S400 to S403 (however, the "first transducer 101a" is replaced with the "second transducer 101b" and the "second transducer 101b" is replaced with the "first transducer 101a"). In step S410, the signal processing unit 111 reads the reception signal (the reception signal acquired in step S113) from the RAM, and determines the predetermined number of peaks from the larger side among the multiple peaks shown by the reception signal as the threshold value. Then, the signal processing unit 111 detects the time Tb from the reception signal until the signal value converges to the determined threshold value or less. In step S413, the propagation time T2' of the acoustic signal from the second transducer 101b to the first transducer 101a is obtained from the cross-correlation function based on the reception signal after the reverberation is removed. The processes in steps S119 to S121 are the same as those in the first embodiment.

[0089] According to the configuration of the present embodiment described above, among the multiple peaks shown by the received signal before the reverberation is removed, the peak that is a predetermined number from the larger side is determined as the threshold. Then, the signal after the timing at which the signal value converges below the determined threshold is removed from the received signal as reverberation, and the propagation time is calculated using the received signal after the reverberation is removed. In this way, the reverberation can be accurately removed, and the propagation time can be accurately calculated. Therefore, it is possible to apply the present embodiment to situations that require high accuracy, such as measuring a minute flow rate. Furthermore, since unnecessary signal values ​​(signal values ​​of reverberation) are not used in calculating the propagation time, it is possible to shorten the processing time required to calculate the propagation time and reduce the power consumption required to calculate the propagation time. Note that, depending on the sampling period (resolution) of the received signal, a peak that should be shown in the received signal may not be shown. In such a case, it is expected that the first to third embodiments can obtain a more accurate propagation time, but the processing of calculating a tentative propagation time and calculating an envelope is increased.

[0090] <Other> The above embodiment merely exemplifies the configuration of the present invention. The present invention is not limited to the above specific embodiment, and various modifications are possible within the scope of the technical concept. For example, in the device of the above embodiment, the propagation time of the acoustic signal is measured, and then the measured value of the propagation time is used to calculate the flow velocity and flow rate of the fluid, but the calculation of the flow velocity and flow rate is not essential. The propagation time measuring device may simply perform a process of measuring the propagation time (at least one of the propagation time T1' and the propagation time T2'). In addition, if the propagation time is simply measured, the propagation angle θ may be 90 degrees. In addition, in the above embodiment, a clamp-on type device that can be attached so as to sandwich the piping is exemplified, but a piping built-in type device configuration may be adopted. In addition, the number of transducers may be three or more, and a transducer pair used for propagating the acoustic signal from the upstream side to the downstream side and a transducer pair used for propagating the acoustic signal from the downstream side to the upstream side may be separated. The transmitting transducer and the receiving transducer need only be positioned at different positions relative to the piping, and do not have to be positioned opposite each other across the piping, nor do they have to be positioned at different positions in the longitudinal direction of the piping.

[0091] The storage unit 112 may store information regarding the timing at which reverberation occurs in advance, and the signal processing unit 111 may remove the reverberation based on the information stored in the storage unit 112. The information regarding the timing at which reverberation occurs may be information indicating the timing itself (time T1+T0 or time T2+T0 in the first embodiment, times Ta and Tb in the second to fourth embodiments, etc.), or may not be information indicating the timing itself. For example, the information regarding the timing at which reverberation occurs may be a propagation time calculated using a received signal including reverberation (propagation time T1 or propagation time T2 in the first embodiment, etc.). The information on the timing of reverberation may be information indicating a propagation time T1′ (the propagation time T2) based on an envelope or a plurality of peaks. The information on the timing of reverberation may be information indicating a threshold value based on an envelope or a plurality of peaks (the threshold value used in the second to fourth embodiments). The information on the timing of reverberation is acquired, for example, in a pre-measurement by a method similar to that described in the first to fourth embodiments. In the main measurement, the signal processing unit 111 reads the information from the storage unit 112 instead of obtaining information on the timing of reverberation. For example, when the propagation time T1′ is obtained by the method of the first embodiment, the processing of steps S101 to S106 in FIG. 3 is performed in the pre-measurement, and the signal processing unit 111 stores the obtained propagation time T1 in the storage unit 112. In the main measurement, the processing of steps S101 to S104 is performed, the signal processing unit 111 reads the propagation time T1 from the storage unit 112, and the processing of steps S107 to S109 is performed.

[0092] Furthermore, when measuring the propagation time T1' or the propagation time T2', the pair of transducers may transmit and receive an acoustic signal multiple times. The signal processor 111 may obtain information about the timing at which reverberation occurs by using multiple reception signals obtained by multiple transmissions and receptions. Since the resolution is apparently increased by synthesizing multiple reception signals, the timing at which reverberation occurs can be obtained more accurately by using multiple reception signals. For example, the propagation time T1, the propagation time T2, the envelope, the peak of the reception signal, and the like can be obtained more accurately. As a result, reverberation can be removed more accurately, and the propagation time can be obtained more accurately.

[0093] In the second to fourth embodiments, the time Ta for removing reverberation from the reception signal obtained by the second transducer 101b (corresponding to the acoustic signal from the first transducer 101a to the second transducer 101b) and the time Tb for removing reverberation from the reception signal obtained by the first transducer 101a (corresponding to the acoustic signal from the second transducer 101b to the first transducer 101a) are determined individually, but the present invention is not limited to this. For example, only one of the times Ta and Tb may be detected and the detected time may be used as the other of the times Ta and Tb. This can reduce the processing load and processing time. For example, if only the time Ta is detected and the time Ta is also used as the time Tb, the processes of steps S210 and S211 in FIG. 7 (detection of the envelope, determination of the threshold, and detection of the time Tb), the processes of steps S210 and S311 in FIG. 9 (detection of the envelope, determination of the threshold, and detection of the time Tb), and the process of step S410 in FIG. 11 (determination of the threshold and detection of the time Tb) are unnecessary.

[0094] Similarly, in the first embodiment, one of the propagation times T1, T2 may be used as the other of the propagation times T1, T2. Also, in the second to fourth embodiments, one of the two thresholds for determining the two times Ta, Tb may be used as the other of the two thresholds. This also makes it possible to reduce the processing load and processing time.

[0095] <Appendix 1> a plurality of transducers (101a, 101b) arranged at different positions relative to a pipe (120) through which a fluid (121) flows, the plurality of transducers (101a, 101b) including at least a first transducer (101a) that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer (101b) that receives the acoustic signal transmitted from the first transducer (101a) and propagated through the fluid (121) in the pipe (120) and converts it into a reception signal as an electric signal; a signal processing unit (111) that removes reverberation from the received signal and obtains a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmitted signal and the received signal after the reverberation has been removed; Equipped with The signal processing unit (111) determining a propagation time of the acoustic signal from the first transducer (101a) to the second transducer (101b) by a cross-correlation analysis between the transmission signal and the reception signal before the reverberation is removed; A signal after a timing based on the obtained propagation time is removed as the reverberation. A propagation time of the acoustic signal from the first transducer (101a) to the second transducer (101b) is again calculated by a cross-correlation analysis between the transmitted signal and the received signal after the reverberation has been removed. A propagation time measuring device (1).

[0096] <Appendix 2> a plurality of transducers (101a, 101b) arranged at different positions relative to a pipe (120) through which a fluid (121) flows, the plurality of transducers (101a, 101b) including at least a first transducer (101a) that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer (101b) that receives the acoustic signal transmitted from the first transducer (101a) and propagated through the fluid (121) in the pipe (120) and converts it into a reception signal as an electric signal; a signal processing unit (111) that removes reverberation from the received signal and obtains a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmitted signal and the received signal after the reverberation has been removed; Equipped with The signal processing unit (111) approximately calculates an envelope of a waveform of the received signal from the received signal before the reverberation is removed, determines a threshold based on the envelope, and removes the signal as the reverberation after the timing when the signal value converges to or below the determined threshold. A propagation time measuring device (1).

[0097] <Appendix 3> a plurality of transducers (101a, 101b) arranged at different positions relative to a pipe (120) through which a fluid (121) flows, the plurality of transducers (101a, 101b) including at least a first transducer (101a) that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer (101b) that receives the acoustic signal transmitted from the first transducer (101a) and propagated through the fluid (121) in the pipe (120) and converts it into a reception signal as an electric signal; a signal processing unit (111) that removes reverberation from the received signal and obtains a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmitted signal and the received signal after the reverberation has been removed; Equipped with The signal processing unit (111) removes, as the reverberation, a signal after a timing at which a signal value converges to or below a predetermined number of peaks from the larger side among a plurality of peaks shown by the received signal before the reverberation is removed. A propagation time measuring device (1). [Explanation of symbols]

[0098] 1: Propagation time measuring device 30: Clamp 100: Device body 101: transducer 101a: first transducer 101b: second transducer 102: Control circuit 103: D / A converter 104: A / D converter 105: Switch 106: Output 110: Transmission signal generating unit 111: Signal processing unit 112: Storage unit 120: Piping 121: Fluid

Claims

1. A plurality of transducers arranged at different positions relative to a pipe through which a fluid flows, the plurality of transducers including at least a first transducer that converts a transmission signal as an electric signal into an acoustic signal, and a second transducer that receives the acoustic signal transmitted from the first transducer and propagated through the fluid in the pipe and converts it into a reception signal as an electric signal; a signal processing unit that removes reverberation from the received signal and obtains a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmitted signal and the received signal after the reverberation has been removed; Equipped with The signal processing unit includes: determining a propagation time of the acoustic signal from the first transducer to the second transducer by a cross-correlation analysis between the transmission signal and the dereverberated reception signal; A signal after a timing based on the obtained propagation time is removed as the reverberation. A cross-correlation analysis between the transmitted signal and the dereverberated received signal is used to again determine a propagation time of the acoustic signal from the first transducer to the second transducer. A propagation time measuring device comprising:

2. The timing is a time when the total time of the propagation time calculated by a cross-correlation analysis between the transmission signal and the reception signal before the reverberation is removed and the time length of the transmission signal has elapsed from the time when the transmission signal is input to the first transducer.

2. The propagation time measuring device according to claim 1 .

3. A storage unit that stores information regarding the timing in advance, 3. The propagation time measuring device according to claim 1, wherein the signal processing unit removes the reverberation based on the information previously stored in the storage unit.

4. The signal processing unit acquires the information regarding the timing using a plurality of reception signals obtained by performing transmission from the first oscillator and reception by the second oscillator a plurality of times.

4. The propagation time measuring device according to claim 1, wherein the propagation time measuring device is a propagation time measuring device.

5. The first transducer and the second transducer are disposed so as to face each other across the pipe.

5. The propagation time measuring device according to claim 1, wherein the propagation time measuring device is a propagation time measuring device.

6. The first transducer and the second transducer are disposed at different positions in the longitudinal direction of the pipe.

6. The propagation time measuring device according to claim 1, wherein the propagation time measuring device is a propagation time measuring device.

7. a switching unit that switches the first transducer so that the transmission signal is input to the second transducer and the first transducer that has received the acoustic signal transmitted from the second transducer outputs the reception signal, The signal processing unit further removes reverberation from the reception signal output from the first transducer, and obtains a propagation time of the acoustic signal from the second transducer to the first transducer by a cross-correlation analysis between the transmission signal input to the second transducer and the reception signal output from the first transducer from which the reverberation has been removed.

7. The propagation time measuring device according to claim 1, wherein the propagation time measuring device is a propagation time measuring device.

8. The signal processing unit determines a flow velocity and / or a flow rate of the fluid in the pipe based on a difference between a propagation time of the acoustic signal from the first transducer to the second transducer and a propagation time of the acoustic signal from the second transducer to the first transducer.

8. A propagation time measuring device according to claim 7.

Images