Measurement pickup for determining the mass flow rate of a liquid

Abstract

Disclosed is a measurement pickup for determining the mass flow rate of a liquid comprising: at least one measurement tube for carrying the liquid having an inlet-side end section and an outlet-side end section; a support body on which the measurement tube is mounted an exciter for exciting vibrations of the measurement tube; at least one vibration sensor for detecting vibrations of the measurement tube; an operating and evaluation circuit for driving the exciter, for receiving the signals of the vibration sensor, and for ascertaining a measured value representing the mass flow rate. The operating and evaluation circuit comprises an adaptive low-pass filter for filtering the sequence of measured values representing the mass flow rate, wherein the low-pass filter has at least one filter parameter dependent on at least one adaptive controlled variable that, for its part, is dependent on the gas content of the test medium.

Claims

1. A measurement pickup for determining the mass flow rate of a liquid measured medium, the measurement pickup comprising: a measurement tube for carrying the liquid measured medium, wherein the measurement tube has an inlet-side end section and an outlet-side end section; a support body on which the measurement tube is mounted in a region of the inlet-side end section and in a region of the outlet-side end section; an exciter for exciting vibrations of the measurement tube in at least one vibration mode; at least two vibration sensors for detecting vibrations of the measurement tube and for outputting vibration-dependent signals; and an operating and evaluation circuit for driving the exciter, for receiving the vibration-dependent signals of the at least two vibration sensors, and for calculating a sequence of mass flow rate measured values as a function of a sequence of time-discrete phase difference values between the vibration-dependent signals of the at least two vibration sensors, wherein the operating and evaluation circuit includes an adaptive low-pass filter for filtering the sequence of the mass flow rate measured values, wherein the adaptive low-pass filter has a filter parameter dependent on a gas content of the liquid measured medium, and wherein the operating and evaluation circuit is designed to calculate the filter parameter as a function of a gas-content-dependent fluctuation of a natural frequency of the measurement tube vibration.

2. The measurement pickup according to claim 1, wherein the filter parameter is independent of the mass flow rate.

3. The measurement pickup according to claim 1, wherein the filter parameter depends on a signal representing the level of the gas content.

4. The measurement pickup according to claim 1, wherein the filter parameter includes a weighting function that weights a current, preliminary mass flow rate measured value in determining a filtered mass flow rate measured value to be output, wherein the weighting function depends on a current weighting time constant τ that is a function of the gas content of the liquid measured medium, and wherein the weighting function is a function that decreases monotonously with an increasing weighting time constant τ.

5. The measurement pickup according to claim 4, wherein the weighting function is a function of the reciprocal value of the weighting time constant τ, including of the type:
(1−exp(−(t/τ)), where t is a period between two preliminary measured values.

6. The measurement pickup according to claim 4, wherein the operating and evaluation circuit is designed: to first ascertain a time constant trend value τ.sub.s as a function of current signals dependent on the gas content of the liquid measured medium, to compare a last valid weighting time constant τ.sub.i-1 with the time constant trend value τ.sub.s, to increase the value of the weighting time constant τ in comparison to the last valid weighting time constant τ.sub.i-1 when the last valid weighting time constant τ.sub.i-1 is less than the time constant trend value τ.sub.s, and to reduce the value of the weighting time constant τ in comparison to the last valid weighting time constant τ.sub.i-1 when the last valid weighting time constant τ.sub.i-1 is greater than the time constant trend value τ.sub.s.

7. The measurement pickup according to claim 6, wherein the following applies to the case of a weighting time constant to be increased: (τ.sub.i−τ.sub.i-1)/t=a, where t is a time interval between two consecutive measured values, and a is a constant which is not more than 8.

8. The measurement pickup according to claim 6, wherein the following applies to the case of time constants to be reduced: τ.sub.i/τ.sub.i-1=b, where 0.9<b<1.

9. The measurement pickup according to claim 6, wherein the operating and evaluation circuit is further designed to ascertain the time constant trend value τ.sub.s as a function of a fluctuation of a natural frequency of the measurement tube vibration normalized with the square of the natural frequency of the measurement tube vibration.

10. The measurement pickup according to claim 9, wherein the operating and evaluation circuit is further designed to ascertain the time constant trend value τ.sub.s as follows: ${\tau }_{s,i}=X.Math.k_{\tau }.Math.\left(\frac{.Math.f_i-f_{i-1}.Math.}{f_i^2}\right)+\left(1-k_{\tau }\right).Math.{\tau }_{s,i-1},$ where X is a constant, f.sub.i is a current natural frequency of a measurement tube vibration, and k.sub.τ, is a constant weighting factor.

11. The measurement pickup according to claim 9, wherein the operating and evaluation circuit is further designed to ascertain the time constant trend value τ.sub.s as follows: ${\tau }_{s,i}=X.Math.\left(\frac{{\mathrm{FFD}}_i}{{\mathrm{FD}}_i^2}\right),\mathrm{where}\text{:}$ ${\mathrm{FFD}}_i=k_{\mathrm{FF}}.Math..Math.f_i-f_{i-1}.Math.+\left(1-k_{\mathrm{FF}}\right).Math.{\mathrm{FFD}}_{i-1},\mathrm{and}\text{:}$ ${\mathrm{FD}}_i=k_F.Math.f_i+\left(1-k_F\right).Math.{\mathrm{FD}}_{i-1},$ and where X is a constant, f.sub.i is a current natural frequency of a measurement tube vibration, and k.sub.FF, and k.sub.F are constant weighting factors.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The invention is now described in further detail on the basis of an exemplary embodiment shown in the drawings.

(2) They show:

(3) FIG. 1a shows an illustration of an exemplary embodiment of a mass flow measuring device having two measurement tubes in side view with a partially removed housing;

(4) FIG. 1b shows the mass flow measuring device of FIG. 1 in a perspective view with a partially removed housing;

(5) FIG. 2a shows mass flow measured values as a function of the gas content in the unfiltered state, filtered with a conventional filter and filtered with the adaptive filter of the mass flow measurement pickup according to the present disclosure;

(6) FIG. 2b shows an illustration enlarged in comparison to FIG. 2a of mass flow measured values as a function of the gas content.

(7) FIG. 3 shows measured values of a mass flow measuring device according to the invention in the filtered and unfiltered state.

DETAILED DESCRIPTION

(8) FIGS. 1a and 1b show, by way of example, a mass flow measuring device 2 with a longitudinal axis L, which mass flow measuring device operates according to the Coriolis principle and is designed according to the present invention. The method according to the invention can in particular be carried out by this mass flow measuring device 2. The mass flow measuring device 2 has two vibratorily mounted measurement tubes A and B, which respectively have an arcuate shape and extend in parallel to one another. The mass flow measuring device 2 can be inserted into a pipeline (not shown) in such a way that the fluid flowing in the pipeline flows through the two measurement tubes A and B. Flow dividers and flow collectors 4, 6 are respectively provided on the inlet side and on the outlet side of the measurement tubes A and B.

(9) An exciter 8 is arranged between the two measurement tubes A and B. In the present embodiment, the exciter 8 is positioned at a reversal point of the arc formed by the two measurement tubes A and B in each case. The two measurement tubes A and B can be periodically deflected in relation to one another by the exciter 8 so that they execute bending vibrations. The two measurement tubes A and B are also mechanically coupled to one another on the inlet side and on the outlet side by corresponding coupling elements 10, 12.

(10) Two vibration sensors 14, 16 respectively extend between the two measurement tubes A and B at an inlet-side and at an outlet-side section thereof. In the present embodiment, the two vibration sensors 14, 16 can each detect the change in distance between the two measurement tubes A, B, i.e., their combined amplitude. The vibration sensors 14, 16 each output a sensor voltage as a function of the vibrations of the measurement tubes A and B. This is initially an analog measurement signal in each case. Excitation of the exciter 8 by application of a corresponding excitation voltage and processing and evaluation of the analog measurement signals provided by the vibration sensors 14, 16 takes place by a correspondingly designed operating and evaluation circuit 18, which is shown in FIGS. 1 and 2 only schematically by a box.

(11) A phase difference Δφ(t) of the vibrations of the two measurement tubes A, B between the two measuring points formed by the vibration sensors 14, 16 is to be determined in order to ascertain measured values of the mass flow rate {dot over (m)}(t). The mass flow rate of the fluid flowing in the pipeline can be determined in particular, as is known to the person skilled in the art, from this phase difference Δφ(t) as follows:
{dot over (m)}(t):=k.Math.tan(Δφ(t)/2)/(2π.Math.f),

(12) where k is a constant specific to the respective measuring device, and f describes the current vibration frequency of the measurement tubes.

(13) In order to determine a sequence Δφ.sub.i=Δφ(t.sub.i) of discrete values of the phase difference Δφ(t), the operating and evaluation circuit 18 respectively determines, from the analog measurement signals provided by the vibration sensors 14, 16, the phase information φ.sub.1(t.sub.i), φ.sub.2(t.sub.i) of the vibration of the two measurement tubes A, B at the measuring points formed by the two vibration sensors 14, 16. For this purpose, in particular the analog measurement signal, which is respectively provided by the vibration sensors 14, 16 and is generally formed by a sensor voltage corresponding to the vibrations, can be sampled at a high sampling frequency, for example at 40 kHz. As a result, a time-discrete, one-dimensional measurement signal having a temporal sequence of measured values is obtained. As a rule, the individual measured values of the sequence are also quantized in order to enable digital processing. Furthermore, in the present embodiment, each measurement signal is converted into an analytical signal consisting of a real part R(t.sub.i) and an imaginary part I(t.sub.i). For this purpose, two filters having a phase difference of 90° may, for example, be used in parallel in a known manner. Furthermore, the present embodiment provides that the data rate of the analytical signal is reduced. This can be done, for example, via corresponding decimation stages. From the analytic signal, as is familiar to the person skilled in the art, (time-dependent) amplitude information A.sub.1(t.sub.i), A.sub.2(t.sub.i) and (time-dependent) phase information φ.sub.1(t.sub.i), φ.sub.2(t.sub.i) of the vibration of the two measurement tubes A, B can respectively be obtained at the respective measuring point. The phase difference Δφ.sub.i=Δφ(t.sub.i) between the two measuring points formed by the vibration sensors 14, 16 can be obtained by forming the difference of the phase information φ.sub.1(t.sub.i), φ.sub.2(t.sub.i) of the vibrations detected by the two vibration sensors 14, 16. Since the signal processing in the present embodiment takes place largely digitally, the respective processed measurement signals are time-discrete so that specific instants t.sub.i are referred to in each case.

(14) The operating and evaluation circuit as a rule evaluates the amplitude information A.sub.1(t.sub.i), A.sub.2(t.sub.i) and the phase information φ.sub.1(t.sub.i), φ.sub.2(t.sub.i). In particular, the excitation of the measurement tubes A, B by the exciter 8 is controlled in each case as a function of these variables.

(15) The operating and evaluation circuit is further designed to determine, from the sequence of time-discrete phase difference values Δφ.sub.i, a corresponding sequence of measured values of the mass flow rate through the pipeline {dot over (m)}.sub.i as follows:
{dot over (m)}.sub.i:=k.Math.tan(Δφ.sub.i/2)/(2π.Math.f.sub.i),

(16) where {dot over (m)}.sub.i:=m(t.sub.i), Δφ.sub.i=Δφ(t.sub.i), and f.sub.i=f(t.sub.i).

(17) Other signal processing or digital or analog signal processing alternative to the above-explained signal processing, such as amplification, zero point compensation, etc., can also be performed in order to provide the sequence of measured values of the mass flow rate {dot over (m)}.sub.i. This can take place both before and after the signal processing described below with an adaptive filter.

(18) The sequence of measured values of the mass flow rate {dot over (m)}.sub.i, or measured values for short, is, in particular with fluctuating gas content of a measured medium flowing in the pipeline, subject to differently strong fluctuations, as shown in FIGS. 2a and 2b, wherein FIG. 2b shows a section with higher resolution. The figures first show unfiltered mass flow measured values. Depending on the gas content, the measured values fluctuate to such an extent that a specific measured value can hardly be identified. The narrower, light-colored band in the middle shows the associated measured value profile on the basis of a conventional low-pass filtering of the measured values. In this respect, it becomes clear, in particular on the basis of the section in FIG. 2b, that these filtered data also have considerable fluctuations. The dark line extending within the light-colored band shows the result of the measured value processing by means of the adaptive filter of the mass flow measuring device according to the invention. It can be seen that the signal has significantly lower fluctuations and can nevertheless follow the actual measured value profile at least as well as the signal processed with a constant filter.

(19) The effective operating principle of the adaptive filter is made possible by controlling a filter parameter as a function of variables which are independent of the phase difference measurement signal, as explained in more detail below.

(20) For operation, the adaptive filter requires a current time constant τ.sub.i:=τ(t.sub.i), which in the present exemplary embodiment is ascertained as a function of the fluctuation of the vibration frequency of the measurement tube vibrations. Wherein the fluctuation, for its part, is also detected via a filter in order to avoid excessive changes in the time constant.

(21) In this case, a current time constant trend value τ.sub.s,i is first ascertained, to which a last valid time constant is compared. If the time constant trend value is greater than the last valid time constant τ.sub.i-1 by more than a tolerance value, a currently valid time constant τ.sub.i with a value greater than that of the last valid time constant is ascertained. If the time constant trend value is smaller than the last valid time constant τ.sub.i-1, a currently valid time constant τ.sub.i with a value smaller than that of the last valid time constant is ascertained.

(22) In this exemplary embodiment, the current time constant trend value τ.sub.s,i is ascertained on the basis of the fluctuation of the frequency of the measurement tube vibrations, which is normalized with the square of this frequency. Specifically, the filter uses the following equation:

(23) ${\tau }_{s,i}=k_{\tau s}.Math.X.Math.\left(\frac{.Math.f_i-f_{i-1}.Math.}{f_i^2}\right)+\left(1-k_{\tau s}\right).Math.{\tau }_{s,i-1}$

(24) The weighting factor k.sub.τs with which a current frequency fluctuation contributes to the current time constant trend value is a few % and can be ascertained as:
k.sub.τs=(1−exp(−t.sub.s/τ.sub.k)),

(25) where t.sub.s is the measurement interval, that is to say the time between two frequency measured values, for example 20 ms. τ.sub.K is a defined time constant which is a multiple of t.sub.s, for example 0.5 s<τ.sub.k<2 s, in particular τ.sub.k=1 s.

(26) The factor X is a constant having a value between 10 and 200, for example 50.

(27) When the thus-ascertained time constant trend value τ.sub.s,i is greater than the sum of the last valid weighting time constant τ.sub.i-1 and a tolerance value, which, for example, corresponds to half a measurement interval t.sub.s, the weighting time constant is increased, for example by addition of a constant time increment, which corresponds to a measurement interval t.sub.s, for example 20 ms.

(28) If, on the other hand, the time constant trend value τ.sub.s,i is less than the last valid weighting time constant τ.sub.i-1, the latter is multiplied by a factor b<1 in order to ascertain a current weighting time constant τ.sub.i, i.e., τ.sub.i=b.Math.τ.sub.i-1, where: 0.9<b<1, in particular 0.95<b<1, preferably 0.975<b<1.

(29) The combination of increases in the weighting time constant by addition of a moderate increment with a reduction by multiplication by a factor <1 has the effect on the one hand that the weighting time constant increases relatively quickly with occurring gas contents starting from small values, but that it then, in contrast to a sequence of multiplications, does not increase exponentially without bounds, and that on the other hand the weighting time constant decreases exponentially as the gas content decreases so that the measured value ascertainment can be adapted quickly to the more favorable conditions. An example of a jump response of the weighting time constant to a jump of the time constant trend value is shown in FIG. 3. Herein, the dotted line shows a jump of the time constant trend value τ.sub.s,i of, for example, 1 s to 10 s and back, while the solid line illustrates the reaction of the weighting time constant τ.sub.i.

(30) The weighting time constant τ.sub.i contributes to the ascertainment of a weighting factor k.sub.m(τ.sub.i) for ascertaining a filtered measured value as follows:

(31) k.sub.m(τ.sub.i)=(1−exp(−(t.sub.s/τ.sub.i)), where t.sub.s is again the measurement interval. A current, filtered measured value {dot over (m)}.sub.f,i={dot over (m)}.sub.f (t.sub.i) is ascertained by a weighted addition of the currently unfiltered measurement signal {dot over (m)}.sub.i={dot over (m)}(t.sub.i) to the previous filtered measurement signal {dot over (m)}.sub.f,i-1={dot over (m)}.sub.f(t.sub.i-1), i.e.:
{dot over (m)}.sub.f,i=k.sub.m(τ.sub.i).Math.{dot over (m)}.sub.i+(1−k.sub.m(τ.sub.i)){dot over (m)}.sub.f,i-1.

(32) As shown in FIGS. 2a and 2b, this leads to a significantly improved damping of the fluctuations of the mass flow measured values.

(33) If the frequency fluctuations are used to ascertain the time constant trend value, the current frequency fluctuations and their normalizing frequency value contribute, as described above, with a weighting factor k.sub.τs to the time constant trend value, which is the function of a time constant.

(34) Since the relative fluctuations of the fluctuations of the frequency are significantly greater than the relative fluctuations of the frequency, it may be expedient to use a larger time constant for the fluctuations than for the normalizing frequency. As shown in the following embodiment.

(35) ${\tau }_{s,i}=X.Math.\left(\frac{{\mathrm{FFD}}_i}{{\mathrm{FD}}_i^2}\right),\mathrm{where}\text{:}$${\mathrm{FFD}}_i=k_{\mathrm{FF}}.Math..Math.f_i-f_{i-1}.Math.+\left(1-k_{k_{\mathrm{FF}}}\right).Math.{\mathrm{FFD}}_{i-1},\mathrm{and}\text{:}$${\mathrm{FD}}_i=k_F.Math.f_i+\left(1-k_F\right).Math.{\mathrm{FD}}_{i-1},$

(36) where X is a constant, f.sub.i is a current natural frequency of a measurement tube vibration, and k.sub.FF, k.sub.F, are weighting factors, in particular constant weighting factors, where:
k.sub.FF,<k.sub.F.

(37) This ensures that a faster adaptation of the adaptive filter to an actually changing vibration frequency of the measurement tubes takes place.

(38) Instead of the fluctuations of the vibration frequency of the measurement tubes, the damping of the vibration amplitudes A.sub.1(t.sub.i), A.sub.2(t.sub.i) at the current resonant frequency or the fluctuation of this damping can also be used. The damping D.sub.i=D(t.sub.i) results, for example, from a ratio of a current excitation signal E.sub.i=E(t.sub.i) for exciting the vibrations and, for example, the sum of the amplitudes at resonance, i.e.:
D.sub.i=E.sub.i/(A.sub.1(t.sub.i)+A.sub.2(t.sub.i))

(39) A normalized current fluctuation of the FD.sub.i damping can be ascertained as follows:
FDi=(Di−Di−1)/Di−1

(40) Based on these variables, a time constant trend value, with which the filter is to be controlled as described above, can be determined with appropriate normalizations in accordance with the methods described in connection with the frequency fluctuation.