Method of ultrasonic flow measurement and device for the implementation of the method
09719822 · 2017-08-01
Assignee
Inventors
Cpc classification
G01F1/667
PHYSICS
G01F1/66
PHYSICS
G01F1/668
PHYSICS
International classification
G01F1/66
PHYSICS
Abstract
A method of ultrasonic clamp-on flow measurement according to the transit time difference method and apparatus for the implementation of the method, wherein the electromechanical transducer element of at least one of the two acoustic transducers is comprised of at least two array elements and a correction factor is determined by comparing the transit times between the acoustic transducers while using different array elements.
Claims
1. Method of ultrasonic clamp-on flow measurement using two acoustic transducers with transducer elements installed on a measuring pipe, where at least one transducer element is comprised of at least two array elements, the method comprising, in an analysis phase carrying out the following steps: a. measurement of a transit time t1 of an acoustic wave traveling between a transducer element of a first acoustic transducer and a first array element of a second acoustic transducer, b. measurement of a transit time t2 of an acoustic wave traveling between the transducer element of the first acoustic transducer and a second array element of the second acoustic transducer, and c. calculation of a time difference Δtc between the transit times t1 and t2, and determining a calibration factor Kaf using the calculated time difference Δtc, and in a subsequent operating phase, using the calibration factor Kaf for measurement of a volumetric flow rate Q in order to compensate for the acoustic influence of the measuring pipe, wherein a position of the first array element of the second acoustic transducer is offset from a position of the second array element of the second acoustic transducer in a flow direction of fluid in the measuring pipe, and wherein: a calculation of a path length difference Δlc is carried out based on the time difference Δtc multiplied by a known sound speed ca in a transducer wedge, resulting in an angle change of the angle of incidence Δalpha with arcsin(Δlc/Δs), where Δs is a distance between the first and second array elements, from which the calibration factor Kaf is calculated as ca/sin(alpha+Δalpha), or the calibration factor Kaf is calculated on the basis of the time difference Δtc according to the equation Kaf=Ka/(1+Ka*cos(alpha)*Δtc/Δs), where Ka is the ratio of a sound speed ca in a transducer wedge and the sine of the angle of incidence alpha in the transducer wedge and Δs is a distance between the first and second array elements, or the calibration factor Kaf is determined from the calculated time difference Δtc using an approximation function fc(Δtc) and the approximation function fc(Δtc) is calculated from a number N of value pairs (Kaf_i, Δtc_i), where the value pairs (Kaf_i, Δtc_i) are obtained before the flow measurement by measuring Kaf_i and Δtc_i on each of a number N of pipes i, on which the acoustic transducers are installed for customary ultrasonic clamp-on flow measurement, where the measurement of Kaf_i is carried out by measuring the transit times tx1 and tx2, while the second acoustic transducer is located at the position x1 or x2, where the transit times are the transit times between a transducer element of the first acoustic transducer and a transducer element of the second acoustic transducer, where the points x1 and x2 lie symmetrically with respect to a point x0, and Kaf_i is calculated as Kafxi=(x2−x1)/(tx2−tx1), and the measurement of Δtc_i is carried out according to steps a through c, while the second acoustic transducer (2) is located at the position x0.
2. Method of ultrasonic clamp-on flow measurement using two acoustic transducers with transducer elements installed on a measuring pipe, where at least one transducer element is comprised of at least two array elements, the method comprising, in an analysis phase, calculating a cross-correlation function of at least two received signals, where a first received signal is a signal emitted by the transducer element of a first acoustic transducer and received and digitized by a first array element of a second acoustic transducer, and a second received signal is a signal emitted by the transducer element of the first acoustic transducer and received and digitized by a second array element of the second acoustic transducer, calculating a temporal position of a maximum of the cross-correlation function, the temporal position being the time difference Δtc, and determining a calibration factor Kaf using a calculated time difference Δtc, and in a subsequent operating phase, using the calibration factor Kaf for measurement of a volumetric flow rate Q in order to compensate for the acoustic influence of the measuring pipe, wherein a position of the first array element of the second acoustic transducer is offset from a position of the second array element of the second acoustic transducer in a flow direction of fluid in the measuring pipe, and wherein: a calculation of a path length difference Δlc is carried out based on the time difference Δtc multiplied by a known sound speed ca in a transducer wedge, resulting in an angle change of the angle of incidence Δalpha with arcsin(Δlc/Δs), where Δs is a distance between the first and second array elements, from which the calibration factor Kaf is calculated as ca/sin(alpha+Δalpha), or the calibration factor Kaf is calculated on the basis of the time difference Δtc according to the equation Kaf=Ka/(1+Ka*cos(alpha)*Δtc/Δs), where Ka is the ratio of a sound speed ca in a transducer wedge and the sine of the angle of incidence alpha in the transducer wedge and Δs is a distance between the first and second array elements, or the calibration factor Kaf is determined from the calculated time difference Δtc using an approximation function fc(Δtc) and the approximation function fc(Δtc) is calculated from a number N of value pairs (Kaf_i, Δtc_i), where the value pairs (Kaf_i, Δtc_i) are obtained before the flow measurement by measuring Kaf_i and Δtc_i on each of a number N of pipes i, on which the acoustic transducers are installed for customary ultrasonic clamp-on flow measurement, where the measurement of Kaf_i is carried out by measuring the transit times tx1 and tx2, while the second acoustic transducer is located at the position x1 or x2, where the transit times are the transit times between a transducer element of the first acoustic transducer and a transducer element of the second acoustic transducer, where the points x1 and x2 lie symmetrically with respect to a point x0, and Kaf_i is calculated as Kafxi=(x2−x1)/(tx2−tx1), and the measurement of Δtc_i is carried out while the second acoustic transducer (2) is located at the position x0.
3. Apparatus for a method of ultrasonic clamp-on flow measurement using first and second acoustic transducers with transducer elements, which are installed on a measuring pipe, where at least the second transducer element is comprised of at least first and second array elements, the apparatus comprising an emitting unit, at least one array multiplexer with at least two inputs, a direction switch, a receiving amplifier, a processing unit, a control unit, and a calculation unit, wherein either the first array element or the second array element or both array elements of the second acoustic transducer are connected to a channel of the direction switch via the array multiplexer, and a second channel of the direction switch is connected to the first transducer element, and the direction switch with two remaining terminals thereof is connected to the emitting unit and the receiving amplifier, and the receiving amplifier is connected to the processing unit, the processing unit configured to determine a time difference Δtc and transmits it to the calculation unit and the control unit is configured to control the direction switch and the array multiplexer as well as the receiving amplifier, the processing unit and the calculation unit, in order to realize an analysis phase as well as an operating phase using the same first and second acoustic transducers without changing position thereof on the measuring pipe, where during the analysis phase the direction switch is configured to be switched so that the emitting unit is connected to terminals of the first acoustic transducer and the transducer element thereof, and the receiving amplifier is connected to the array multiplexer, and the array multiplexer is configured to be controlled so that either only the first array element or the second array element of the second transducer element is connected to the receiving amplifier, and during the operating phase the array multiplexer is configured to be switched by the control unit so that the array elements of the second transducer are connected in parallel.
4. The apparatus according to claim 3, wherein the processing unit is comprised of an analog-to-digital converter, a buffer memory and a digital correlator.
5. The apparatus according to claim 4, wherein the correlator comprises a digital signal processor.
6. The apparatus according to claim 3, wherein the processing unit comprises a time measuring unit, a buffer memory and a difference generator.
7. The apparatus according to anyone of claims 3 to 6 in that, wherein the first transducer element is comprised of at least two array elements and the array multiplexer connects either the first array element or the second array element or both array elements of the first transducer element with a channel of the direction switch.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) In the following, the design examples of the invention are explained in detail with the help of the illustrations. The figures show the following:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
DETAILED DESCRIPTION OF THE INVENTION
(10) The arrangement shown in
(11) The idea behind the invention is to measure this deviation from parallelism of the wave front to the transducer element of the acoustic transducer receiving the sound wave from the transit time differences which occur along the geometry of the receiving transducer element. In order to achieve that, it is suggested to design the transducer element of at least one of the two acoustic transducers as an array consisting of at least two array elements.
(12)
delta_tc=t2−t1 Eq. (5)
(13) During an analysis phase, the two transit times t1 and t2 are therefore measured and the time difference delta_tc is calculated to obtain a measured value of the time difference delta_tc. The measured value of the time difference is used to determine a calibration factor Kaf. During the subsequent operating phase, the flow measurement is carried out. The calibration factor determined during the analysis phase is used to compensate for the acoustic influence of the measuring pipe. To achieve this, the transducer constant Ka in Eq. (4) is replaced by Kaf.
(14) During the operating phase, the two array elements 4a and 4b of the transducer element 4 are electrically connected in parallel. Preferably, the two combined array elements of the transducer element 4 form a transducer element of the same size as the transducer element 3. Due to the parallel connection, the transducer element 4 has about the same acoustic properties as the transducer element 3.
(15) One advantageous embodiment of the invention is obtained by determining the time difference delta_tc directly from the received signals without a previous measurement of the transit times t1 and t2. To achieve this, the cross correlation function of the received signals s1 and s2 is calculated, where s1 is the signal emitted by the transducer element 3 of the first acoustic transducer 1 and received and digitized by the array element 4a of the second acoustic transducer 2, and s2 is the signal emitted by the transducer element 3 of the first acoustic transducer 1 and received and digitized by the array element 4b of the second acoustic transducer 2. If the two signals s1 and s2 have the same transit time, the maximum of the cross-correlation function over time is at the point t=0. A difference of the transit times of the two signals manifests itself in a displacement of the maximum of the cross-correlation function by exactly this difference. The time difference delta_tc can therefore be determined by determining the maximum of the cross-correlation function of the signals s1 and s2.
(16) Another advantageous embodiment of the invention is described in the following. The path length difference delta_tc is calculated as the time difference delta_tc multiplied by the sound speed ca in the transducer wedge.
Delta_tc=ca*delta_tc Eq. (6)
(17) The sound speed ca in the transducer wedge can be assumed as known. The change of the angle delta_alpha is therefore
Delta_alpha=arcsin(delta_lc/delta_s) Eq. (7)
(18) Here, delta_s is the distance between the array elements as defined in
Kaf=ca/sin(alpha+delta_alpha) Eq. (8)
(19) Another possible embodiment of the invention results from the fact that the equation for the calculation of the calibration factor Kaf on the basis of the time difference delta_tc is derived from the following considerations.
(20) The imagined displacement of the array element 4a by the distance delta_x as shown in
Kafx=delta_x/delta_tx Eq. (9)
(21) In order to avoid the physical displacement, the transit time t2 is instead measured and used to calculate the transit time t20 between the transducer element 3 and the array element 4 at the imagined position. The transit time t20 is obtained by from t2 by subtracting the transit time ta along the distance delta_i through the transducer wedge.
t20=t2−ta Eq. (10)
(22) Therefore, the time difference delta_tx is according to Eq. (5):
delta_tx=t1−t20=t1−(t2−ta)=delta_tc+ta Eq. (11)
(23) Therefore, Eq. (9) yields
Kafx=delta_x/(delta_tc+ta) Eq. (12)
(24) The transit time to results from the distance delta_i and the sound speed ca in the transducer wedge as
ta=delta_l/ca Eq. (13)
(25) Delta_l, delta_s and delta_x form a right triangle. Therefore, delta_x can be replaced by delta_s using the angle alpha. Therefore:
delta_x=delta_s/cos(alpha) Eq. (14)
(26) In Eq. (13), delta_l is replaced by the product of tan(alpha) and delta_s:
ta=tan(alpha)*delta_s/ca. Eq. (15)
(27) Inserting Eq. (14) and Eq. (15) in Eq. (12) yields:
Kafx=delta_s/(cos(alpha)(delta_tc+tan(alpha)*delta_s/ca)) Eq. (16)
(28) With Ka according to Eq. (2), Eq. (16) can be simplified to
Kafx=Ka/(1+Ka*cos(alpha)*delta_tc/delta_s) Eq. (17)
(29) Therefore, the measurement of the calibration factor Kafx according to Eq. (17) replaces the spatial displacement delta_x with an equivalent measurement of the time difference delta_tc using the parameters of the acoustic transducer. Therefore, Eq. (17) provides a formula for determining the calibration factor Kaf:
Kaf=Ka/(1+Ka*cos(alpha)*delta_tc/delta_s) Eq. (18)
(30) Eq. (8) and Eq. (18) are alternative ways to describe the calibration factor Kaf. The difference between Kaf and the transducer constant Ka, which results only from the properties of the transducer, is caused by the effect of the pipe wall. A special case is where there is no pipe wall influence at all, meaning that delta_tc=0. In this case, Eq. (8) and Eq. (18) yield Kaf=Ka.
(31) Eq. (8) yields the same result as Eq. (18). Namely, Eq. (8) can be converted to Eq. (18) by approximating sin(alpha+delta_alpha) using the first term of the Taylor series and substituting arcsin(delta_lc/delta_s) with the argument delta_lc/delta_s.
(32) Another possible embodiment of the invention is obtained by experimentally determining the relationship between the time difference delta_tc and the calibration factor Kaf. By doing this, Kaf can be determined using the method described in DE 102004031274 B4. This method is applied in a laboratory to a sufficient choice of measuring pipes which are possible in practical applications. For this, the acoustic transducers are positioned on pipe 5 as shown in Error! Reference source not found. Then, the transit times tx1 and tx2 at the positions x1 and x2 of the acoustic transducer 2 are measured, where the points x1 and x2 lie symmetrically with respect to a point x0. The transit times tx1 and tx2 are the transit times between the transducer element 3 of the acoustic transducer 1 and the transducer element 4 of the acoustic transducer 2. The two array elements 4a and 4b of the transducer element 4 are electrically connected in parallel.
(33) As described in DE 102004031274 B4, the calibration factor Kaf, which already contains the influence of the pipe wall, is calculated as
Kaf=(x2−x1)/(tx2−tx1). Eq. (19)
(34) Afterwards, the acoustic transducer 2 is placed at position x0 and the time difference delta_tc is measured. This yields a value pair Kaf_i and delta_tc_i for the pipe i used in the experiment. The experiment is repeated for a number N of additional pipe segments which are selected in such way that the extent of their influence on the calibration factor varies. This yields a number N of value pairs (Kaf_i, delta_tc_i). Applying a common method, these value pairs are used to calculate an approximation function fc(delta_tc) suitable to approximate Kaf from delta_tc.
(35) During the operating phase, this can be used to calculate the calibration factor Kaf from the measured value delta_tc:
Kaf=fc(delta_tc) Eq. (20)
(36) The function fc serves to cover as much as possible of the variation range of Kaf which can occur in practice. To achieve this, the variation range of the properties of the pipes used in the experiment has to correspond to that of the pipes present in the practical application of the flowmeter. For example, if the acoustic transducers are suitable for steel and stainless steel pipes with pipe wall thicknesses between 4 and 8 mm, one possibility would be to conduct the experiment on pipes with the pipe wall thicknesses 4, 6 and 8 mm for both steel and stainless steel. The data base can be extended by conducting the experiment using different steel or stainless steel grades with varying elastic properties.
(37) One advantageous embodiment of the invention results from the fact that the transducer element 4 shown in
(38) The number of array elements can be continued in the manner described. With the increasing number of array elements, the transducer element used during the analysis phase becomes more similar to the interconnection of all array elements used for the flow measurement. On the other hand, the time difference delta_tc, which causes a certain deviation of the angle of the wave front decreases, meaning that it can be measured less exactly.
(39) The properties of the two acoustic transducers used in a flow measurement arrangement according to the time difference principle should be as similar as possible in order to minimize the zero offset error. One advantageous embodiment of the invention results from the fact that, during the operating phase, the two array elements 4a and 4b are connected in parallel and, by doing this, the size of the transducer element of the acoustic transducer 2 is made largely identical to that of the acoustic transducer 1. Furthermore, it is possible to realize the transducer element 3 of the acoustic transducer 1 as an array. One such arrangement is shown in
(40) Using the arrangement shown in
(41) Other possible embodiments of the invention are achieved, as previously described, by increasing the number of array elements.
(42) The measurement of the time difference delta_tc is only intended to measure the influence of the deviation of the wave front from parallelism to the receiving transducer element. However, if the fluid in the pipe is flowing, the transit time of the sound is also influenced by the flow. If the flow velocity is constant during the measurement of t1 and t2, the flow has no effect on the difference. However, in practice, this is never achieved in an ideal manner. In practice, the flow is mostly in a turbulent state. In this case, the flow velocity randomly fluctuates around a mean value. The influence of such random fluctuations in case of a constant mean value of the flow velocity can be eliminated by averaging several successive measured values of the time difference. However, if the mean value of the flow velocity is not constant, this procedure will not help. For example, a time-linear increase of the flow velocity results in that the flow velocity has a systematically greater influence on the transit time t1 than on the transit time t2. The pause times between the individual measurements should be kept as short as possible. This procedure can also be applied to the present invention. Obviously, the effect of the change of flow is smaller if the time intervals between the measurements of t1 and t2 are shorter.
(43) For the implementation of the described method, a device consisting of an emitting unit S, at least one array multiplexer AMUX, a direction switch DMUX, a receiving amplifier V, a processing unit SPU and a control unit CTRL and a calculation unit CALC is proposed (
(44) During the analysis phase, the direction switch DMUX is switched in such way that the emitting unit S is connected to the terminals of the acoustic transducer 1 and its corresponding transducer element, and the receiving amplifier V is connected to the array multiplexer AMUX. The array multiplexer AMUX is controlled in such way either only the array element 4a or the array element 4b of the acoustic transducer 2 is connected to the receiving amplifier V.
(45) The processing unit SPU connected to the receiving amplifier V determines the time difference delta_tc to be measured. The processing unit SPU shown in
(46) The calibration factor Kaf is calculated by the calculation unit CALC from the time difference delta_tc determined in the correlator according to Eq. 8 or Eq. 18.
(47) In another design variation (
(48) During the operating phase, the array multiplexer AMUX is switched by the control unit CTRL in such way that the array elements 4a and 4b are connected in parallel. The direction switch DMUX alternately operates the acoustic transducer 1 as an emitter and the acoustic transducer 2 as a receiver, and the acoustic transducer 2 as an emitter and the acoustic transducer 1 as a receiver. The processing unit SPU is now used to measure the transit time difference Δt caused by the flow and the transit time in the fluid t1, and uses them to calculate the flow rate according to Eq. 4 in the calculation unit CALC, where the transducer constant Ka is substituted by the calibration factor Kaf determined during the analysis phase.
(49) In another design variation (
(50) In another design variation, at least one of the transducer elements 3 and 4 is made up of more than two array elements. In this case, the corresponding array multiplexer of the devices according to