Method and apparatus for adaptive line enhancement in Coriolis mass flow meter measurement

• US 5,555,190 A
• Filed: 07/12/1995
• Issued: 09/10/1996
• Est. Priority Date: 07/12/1995
• Status: Expired due to Term

• Current Assignee: Micro Motion Incorporated (Emerson Electric Company)
• Original Assignee: Micro Motion Incorporated (Emerson Electric Company)
• Inventors: Derby, Howard V., Rajan, Seeraman, Bose, Tamal
• Primary Examiner(s): Trammell, James P.
• Application Number:
  US08/501,411
• Time in Patent Office: 426 Days
• Field of Search: 364/509, 364/510, 364/572, 73/861.38, 73/861.37, 73/861.03, 73/861, 324/76, 324/77-85, 324/601
• US Class Current: 702/45
• CPC Class Codes: G01F 1/8431 electronic circuits
• CPC Class Codes: G01F 1/8436 signal processing
• CPC Class Codes: G01F 1/8477 with multiple measuring conduits
• CPC Class Codes: G01F 15/024 involving digital counting

Abstract

An apparatus and method for determining frequency and phase relationships of vibrating flow tubes in a Coriolis mass flow meter. Adaptive line enhancement (ALE) techniques and apparatus are used in a digital signal processing (DSP) device to accurately determine frequency and phase relationships of the vibrating flow tube and to thereby more accurately determine mass flow rate of a material flowing through the mass flow meter. In a first embodiment, an adaptive notch filter is used to enhance the signal from each corresponding sensor signal on the vibrating flow tubes. In a second embodiment, a plurality of adaptive notch filters are cascaded to enhance the signal from each corresponding sensor signal. In both embodiments an anti-aliasing decimation filter associated with each sensor signal reduces the computational complexity by reducing the number of samples from a fixed frequency A/D sampling device associated with each sensor signal. Computational adjustments are performed to compensate for spectral leakage between the fixed sampling frequency and the variable fundamental frequency of the vibrating flow tubes. Despite this added computational complexity, the present invention is simpler than prior designs and provides better noise immunity due to the adaptive notch filtration. Heuristics are applied to the weight adaptation algorithms of the notch filters to improve convergence of the digital filters and to reduce the possibility of instability of the filters interfering with mass flow measurements.

Claims

1. An apparatus for measuring mass flow rate of a material in a Coriolis mass flow meter having a flow tube and having first and second sensors associated with said flow tube for generating output signals indicative of the oscillatory movement of said flow tube, said apparatus comprising:

analog to digital converter means for periodically sampling said sensor output signals and for converting sampled sensor output signals to digital form to generate a sequence of discrete sampled values representative of said output signals, including any undesirable components, of each of said first and second sensors;

digital notch filtration means, responsive to the generation of said sequence of discrete sampled values, for generating a sequence of discrete enhanced values, each discrete enhanced value corresponding to a sample of said sequence of discrete sampled values with said undesirable components removed;

phase value determination means, responsive to the generation of said sequence of discrete enhanced values, for generating the phase values of the oscillatory movement of said flow tube indicated by said sequence of discrete enhanced values;

phase difference means, responsive to the generation of said phase values, for determining a phase difference between the output signals of said first and second sensors; and

mass flow measurement means, responsive to the determination of phase difference, for determining a mass flow rate value of the material flowing through the flow tube.

2. The apparatus of claim 1 further comprising:

notch adaptation means, cooperative with said digital notch filtration means, for altering filter parameters of said digital notch filtration means to affect the notch capability to reject undesirable components of the output signals of said first and second sensors.

3. The apparatus of claim 2 wherein said filter parameters comprise variable polynomial coefficients applied to said discrete sampled values to enhance said discrete sampled values, wherein said variable polynomial coefficients determine the center frequency of the notch of said digital notch filtration means.

4. The apparatus of claim 3 wherein said notch adaptation means further comprises weight adaptation means for adjusting said variable polynomial coefficients applied to said discrete sampled values in enhancing said discrete sampled values to alter the center frequency of said digital notch filtration means.

5. The apparatus of claim 4 further comprising:

means for detecting that the ratio of a discrete enhanced value to a corresponding noise signal has fallen below a predetermined threshold value and for generating a fault signal responsive to said detection; and

means responsive to generation of said fault signal for adjusting said variable polynomial coefficients.

6. The apparatus of claim 4 further comprising:

stability test means for detecting instability in the adjustment of said variable polynomial coefficients and for generating an instability signal to indicate detection of the instability; and

means responsive to generation of said instability signal to adjust said variables polynomial coefficients to reduce said instability.

7. The apparatus of claim 2 wherein said filter parameters comprise a variable debiasing parameter applied to said discrete sampled values in enhancing said discrete sampled values, wherein said variable debiasing parameter determines the frequency spectrum width of the notch of said digital notch filtration means.

8. The apparatus of claim 7 wherein said notch adaptation means further comprises weight adaptation means for adjusting said variable debiasing parameter applied to said discrete sampled values in enhancing said discrete sampled values.

9. The apparatus of claim 8 further comprising:

means for detecting that the ratio of a discrete enhanced value to a corresponding noise signal has fallen below a predetermined threshold value and for generating a fault signal responsive to the detection; and

means responsive to said fault signal for adjusting said variable debiasing parameter.

10. The apparatus of claim 1 wherein said phase difference means further comprises:

windowing means for defining a window comprising a plurality of sequential ones of said discrete enhanced values; and

Goertzel filtration means for decimating the discrete enhanced values in said window to determine a phase value for said window.

11. The apparatus of claim 10 further comprising Hanning windowing means for weighting each of said plurality of discrete enhanced values in said window wherein the weights are determined as:

space="preserve" listing-type="equation">h(k)=(1/2) (1-cos (2π

k/(2N-1)))

where;

N is half the number of discrete enhanced values in the window, and

k is the index of the value to which the weight h(k) is applied.

14. The apparatus of claim 1 wherein said digital notch filtration means further comprises:

first digital notch filtration means, responsive to the generation of said sequence of discrete sampled values, for generating an intermediate sequence of discrete values, each intermediate discrete value corresponding to a sample of said sequence of discrete sampled values with signals representative of said undesirable components partially removed; and

second digital notch filtration means, responsive to the generation of said intermediate sequence of discrete values, for generating said sequence of discrete enhanced values, each discrete enhanced value corresponding to an intermediate discrete value of said intermediate sequence of discrete values with signals representative of said undesirable components removed.

15. The apparatus of claim 14 further comprising:

first notch adaptation means, cooperative with said first digital notch filtration means, for altering filter parameters of said first digital notch filtration means which are determinative of the characterization of signals as representative of noise; and

second notch adaptation means, cooperative with said second digital notch filtration means, for altering filter parameters of said second digital notch filtration means which are determinative of the characterization of signals as representative of noise.

16. The apparatus of claim 15 wherein said filter parameters comprise variable polynomial coefficients applied to said discrete sampled values to enhance said discrete sampled values.

17. The apparatus of claim 16 wherein said first notch adaptation means further comprises weight adaptation means for adjusting said variable polynomial coefficients applied to said discrete sampled values to generate said intermediate sequence of discrete values, and wherein said second notch adaptation means further comprises weight adaptation means for adjusting variable polynomial coefficients applied to said intermediate sequence of discrete values to generate said sequence of discrete enhanced values.

12. The apparatus of claim I wherein said phase difference means further comprises:

windowing means for defining a plurality of windows, each of said plurality of windows comprising a plurality of sequential ones of said discrete enhanced values; and

Goertzel filtration means for decimating said plurality of discrete enhanced values in each of said plurality of windows to determine a phase value for said each of said plurality of windows, wherein each of said windows comprises an equal number of said discrete enhanced values and wherein each of said windows is offset from an earlier window by an equal number of said discrete enhanced values.

13. The apparatus of claim 12 further comprising Hanning windowing means for weighting each of said plurality of discrete enhanced values in each of said windows wherein the weights are determined as:

space="preserve" listing-type="equation">h(k)=(1/2) (1-cos (2π

k/(2N-1)))

where;

N is half the number of discrete enhanced values in the window, and

k is the index of the value to which the weight h(k) is applied.

18. In a Coriolis mass flow meter having a flow tube and having first and second sensors associated with the flow tube for generating output signals indicative of the oscillatory movement of the flow tube, a method for measuring the mass flow rate of a material flowing through said flow tube of said flow meter comprising the steps of:

periodically converting analog output signals generated by the first and second sensors into digital form to generate a sequence of discrete sampled values representative of said output signals, including any undesirable components, of each of said first and second sensors;

applying said sequence of discrete sampled values to digital notch filtration means to generate a sequence of discrete enhanced values, each discrete enhanced value corresponding to a sample of said sequence of discrete sampled values with signals representative of noise removed;

applying said sequence of discrete enhanced values to phase value determination means to determine phase values of the oscillatory movement of the flow tube indicated by said sequence of said discrete enhanced values;

applying said phase values to phase difference computation means to determine a phase difference between the output signals of said first and second sensors; and

determining the mass flow rate of the material flowing through said flow meter responsive to the determination of phase difference.

19. The method of claim 18 further comprising:

altering filter parameters of said digital notch filtration means to adjust the digital notch filtration means to compensate for changes in the frequency of oscillations of the flow tube.

20. The method of claim 19 wherein said filter parameters comprise variable polynomial coefficients applied to said discrete sampled values to enhance said discrete sampled values, wherein said variable polynomial coefficients determine the center frequency of the notch of said digital notch filtration means.

21. The method of claim 20 wherein the altering step further comprises adjusting said variable polynomial coefficients applied to said discrete sampled values in enhancing said discrete sampled values to alter the center frequency of said digital notch filtration means.

22. The method of claim 21 further comprising:

determining a ratio of a discrete enhanced value to a corresponding noise signal;

determining whether said ratio has fallen below a predetermined threshold value;

generating a fault signal responsive to the determination that said ratio has fallen below said predetermined threshold value; and

adjusting said variable polynomial coefficients responsive to generation of said fault signal.

23. The method of claim 21 further comprising:

determining whether said variable polynomial coefficients are outside an acceptable range of stable values;

generating an instability signal responsive to a determination that said variable polynomial coefficients are unstable; and

adjusting said variable polynomial coefficients responsive to generation of said instability signal to reduce said instability.

24. The method of claim 19 wherein said filter parameters comprise a variable debiasing parameter applied to said discrete sampled values in enhancing said discrete sampled values, wherein said variable debiasing parameter determines the frequency spectrum width of the notch of said digital notch filtration means.

25. The method of claim 24 wherein said step of altering said parameters further comprises adjusting said variable debiasing parameter applied to said discrete sampled values for enhancing said discrete sampled values.

26. The method of claim 25 further comprising:

determining a ratio of a discrete enhanced value to a corresponding noise signal;

determining whether said ratio has fallen below a predetermined threshold value;

generating a fault signal responsive to the determination that said ratio has fallen below said predetermined threshold value; and

adjusting said variable debiasing parameter responsive to generation of said fault signal.

27. The method of claim 18 wherein the application of said phase values to said phase difference computation means further comprises:

defining a window comprising a plurality of sequential ones of said discrete enhanced values; and

decimating the discrete enhanced values in said window through a Goertzel filtration to determine a phase value for said window.

28. The method of claim 27 further comprising:

determining a Hanning window weight for weighting each of said plurality of discrete enhanced values in said window wherein the weights are determined as;

space="preserve" listing-type="equation">h(k)=(1/2)(1-cos (2π

k/(2N-1)))

where;

N is half the number of discrete enhanced values in the window, and

k is the index of the value to which the weight h(k) is applied.

29. The method of claim 18 wherein the application of said phase values to said phase difference computation means further comprises:

defining a plurality of windows, each of said plurality of windows comprising a plurality of sequential ones of said discrete enhanced values; and

decimating the discrete enhanced values in each of said plurality of windows through a Goertzel filter to determine a phase value for said each of said plurality of windows, wherein each of said windows comprises an equal number of said discrete enhanced values and wherein each of said windows is offset from an earlier window by an equal number of said discrete enhanced values.

30. The method of claim 29 further comprising:

determining a Hanning window weight for weighting each of said plurality of discrete enhanced values in each of said plurality of windows wherein the weights are determined as;

space="preserve" listing-type="equation">h(k)=(1/2) (1-cos (2π

k/(2N-1)))

where;

N is half the number of discrete enhanced values in the window, and

k is the index of the value to which the weight h(k) is applied.

31. The method of claim 18 wherein the step of applying said discrete sampled values to said digital notch filtration means further comprises:

operating a first digital notch filtration means to generate an intermediate sequence of discrete values, each intermediate discrete value corresponding to a sample of said sequence of discrete sampled values with signals representative of noise partially removed; and

operating a second digital notch filtration means responsive to the generation of said intermediate sequence of discrete values to generate a sequence of discrete enhanced values, each discrete enhanced value corresponding to an intermediate discrete value of said intermediate sequence of discrete values with signals representative of noise removed.

32. The method of claim 31 further comprising:

operating first notch adaption means to alter filter parameters of the first digital notch filtration means to adjust the first digital notch filtration means to compensate for changes in the frequency of oscillations of the flow tube; and

operating second notch adaption means to alter filter parameters of the second digital notch filtration means to adjust the second digital notch filtration means to compensate for changes in the frequency of oscillations of the flow tube.

33. The method of claim 32 wherein said filter parameters comprise variable polynomial coefficients applied to said discrete sampled values to enhance said discrete sampled values.

34. The method of claim 33 wherein said step of altering said first notch adaptation means further comprises .operating weight adaptation means for adjusting variable polynomial coefficients applied to said discrete sampled values to generate said intermediate sequence of discrete values and wherein said step of operating said second notch adaptation means further comprises operating weight adaptation means for adjusting variable polynomial coefficients applied to said intermediate sequence of discrete values to generate said sequence of discrete enhanced values.

35. An apparatus for measuring the mass flow rate of a material in a Coriolis mass flow meter having a flow tube and having first and second sensors associated with the flow tube for generating output signals indicative of the oscillatory movement of said flow tube, said apparatus comprising:

analog to digital converter means for periodically sampling sensor output signals at a fixed rate and for converting sampled sensor output signals to digital form to generate a sequence of discrete sampled values representative of said output signals of each of said first and second sensors;

digital filtration means, responsive to the generation of said sequence of discrete sampled values, for generating a sequence of discrete enhanced values;

phase difference means, responsive to the generation of said sequence of discrete enhanced values, for determining a phase difference between said output signals of said first and second sensors; and

mass flow measurement means, responsive to the determination of said phase difference, for determining a mass flow rate value of the material flowing through said flow tube.