Beschreibung

Verfahren zur Signalverarbeitung für Messsignale eines Wirbeldurchflussmessaufnehmers

[0001] Die Erfindung betrifft ein Verfahren zur Signalverarbeitung für Messsignale eines Wirbeldurchflussmessaufnehmers .

[0002] Wirbeldurchflussmessaufnehmern werden in der industriellen Messtechnik zur Messung von Volumenströmen eingesetzt.

[0003] Sie arbeiten nach dem Prinzip der Karman' sehen Wirbelstraße. Dabei wird ein Medium, dessen Volumendurchfluss gemessen werden soll, durch ein Messrohr geleitet, in dem ein Staukörper eingesetzt ist. Hinter dem Staukörper bilden sich abwechselnd beidseitig Wirbel mit entgegen gesetztem Drehsinn. Diese Wirbel erzeugen jeweils einen lokalen Unterdruck. Die Druckschwankungen werden von dem Messaufnehmer erfasst und in elektrische Messsignale umgewandelt. Die Wirbel bilden sich innerhalb der zulässigen Einsatzgrenzen des Messaufnehmers sehr regelmäßig aus. Die Anzahl der pro Zeiteinheit erzeugten Wirbel verhält sich proportional zum Volumendurchfluss.

[0004] Zur Erfassung der Druckschwankungen eignen sich z.B. kapazitive Sensoren, wie die in der Europäischen Patentschrift EP-Bl 0 229 933 beschrieben sind.

[0005] Die vom Messaufnehmer abgeleiteten elektrischen Messsignale weisen im Healfall einen sinusförmigen Verlauf auf. Hieraus wird heutzutage die Anzahl der pro Zeiteinheit erzeugten Wirbel anhand von Nulldurchgängen des elektrischen Signals bestimmt. Dabei wird z.B. die Zeit von einem Nulldurchgang zum nächsten gemessen, oder es wird z.B. in einem Zeitintervall die Zeit vom ersten zum letzten Nulldurchgang und die Anzahl der in dem Zeitintervall auftretenden Nulldurchgänge bestimmt.

[0006] Unter realen Messbedingungen kann der Fall eintreten, dass dem im Healfall sinusförmigen Signal Störungen überlagert sind, die dazu führen können, dass ein zusätzlicher Nulldurchgang auftritt, bzw. dass ein Nulldurchgang verschwindet. Solche Messfehler sind nur sehr schwer zu erkennen.

[0007] Ein weiterer Nachteil derartiger Verfahren besteht darin, dass zusätzlich im Signal vorhandene Informationen, wie z.B. dessen Amplitude oder Oberschwingungen, ungenutzt bleiben.

[0008] Diese Informationen lassen sich nutzen, indem beispielsweise mittels einer schnellen Fouriertransformation (FFT) eine Frequenzbestimmung durchgeführt wird. Derartige Verfahren erfordern jedoch, dass das Messsignal mit einer hohen Abtastrate

digitalisiert wird und eine große Anzahl von Abtastwerten zumindest vorübergehend gespeichert wird. Außerdem erfordert diese Form der Frequenzbestimmung eine große Anzahl von mathematischen Operationen. [0009] Wirbeldurchflussmesser sind heute bereits als so genannte 2-Draht Geräte erhältlich. Derartige 2-Draht Geräte weisen zwei Anschlussleitungen auf, über die sowohl die Energieversorgung des Geräts als auch die übertragung von Messergebnissen erfolgt. In der Industrie hat sich ein Standard für diese 2-Draht Geräte durchgesetzt, gemäß dem die Energieversorgung über eine 24 Volt Spannungsquelle erfolgt und die 2-Draht Geräte einen durch die Anschlussleitungen fließenden Strom in Abhängigkeit von den zu übertragenden Messergebnissen auf Werte zwischen 4 mA und 20 niA regeln. Entsprechend steht diesen Messaufnehmern nur eine geringe elektrische Leistung zur Verfügung. Hierdurch ist die erzielbare Rechenleistung eng begrenzt.

[0010] Es ist eine Aufgabe der Erfindung ein Verfahren zur Signalverarbeitung für Messsignale eines Wirbeldurchflussmessaufnehmers mit geringem Bedarf an Rechenleistung und Speicherplatz anzugeben. [0011] Hierzu besteht die Erfindung in einem Verfahren zur Signalverarbeitung für

Messsignale eines Wirbeldurchflussmessaufnehmers zur Messung eines Durchflusses eines Medium durch ein Messrohr,

[0012] - der einen im Messrohr angeordneten Staukörper, und [0013] - einen Sensor zur Erfassung von im Bereich des Staukörpers auftretenden [0014] Druckschwankungen und zur Umwandlung dieser Druckschwankungen in [0015] ein elektrisches Messsignal aufweist, bei dem

[0016] - zumindest ein Teil des Messsignals abgetastet und digitalisiert wird, [0017] - eine Autokorrelation des digitalisierten Messsignals berechnet wird, und [0018] - der Durchfluss anhand mindestens einer Eigenschaft der Autokorrelation [0019] abgeleitet wird.

[0020] Gemäß einer Ausgestaltung des Verfahrens ist die Eigenschaft eine Frequenz oder eine Periode der Autokorrelation, und die Frequenz oder die Periode der Autokorrelation wird gleich einer Frequenz oder einer Periode des Messsignals gesetzt.

[0021] Gemäß einer ersten Weiterbildung wird eine Nullstelle oder ein Minimum der Autokorrelation bestimmt, und die Eigenschaft der Autokorrelation wird anhand der Lage der Nullstelle oder des Minimums bestimmt wird. [0022] Gemäß einer Weiterbildung wird die Lage des Minimums durch eine Anpassung

einer Parabel an Punkte der Autokorrelation bestimmt wird. [0023] Gemäß einer anderen Weiterbildung ist die Eigenschaft der Autokorrelation die

Frequenz oder Periode derselben und wird durch eine Anpassung einer Kosinus

Funktion an die Autokorrelation ermittelt. [0024] Gemäß einer Ausgestaltung durchläuft das Messsignal vor der Bildung der

Autokorrelation einen adaptiven Filter. [0025] Gemäß einer Weiterbildung der letztgenannten Ausgestaltung ist der adaptive Filter ein adaptiver Bandpassfilter. [0026] Gemäß einer Weiterbildung der letztgenannten Ausgestaltung wird eine Frequenz der Autokorrelation bestimmt, und ein Frequenzbereich in dem der Bandpassfilter durchlässig ist anhand der Frequenz der Autokorrelation eingestellt. [0027] Gemäß einer anderen Weiterbildung der letztgenannten Ausgestaltung wird als adaptiver Bandpassfilter ein Adaptive Line Enhancer verwendet. [0028] Gemäß einer anderen Weiterbildung der letztgenannten Ausgestaltung ist der adaptive Filter ein Notch-Filter. [0029] Gemäß einer Weiterbildung der letztgenannten Weiterbildung führt der adaptive

Notch-Filter eine Filterfunktion, die ein Restsignal liefert, das zur Optimierung des

Filters dient, und eine hierzu komplementäre Filterfunktion aus, die das Messsignal liefert, anhand dessen die Autokorrelation bestimmt wird. [0030] Weiter besteht die Erfindung in einem Wirbeldurchflussmessaufnehmers zur

Messung eines Durchflusses eines Mediums durch ein Messrohr, mit [0031] - einen im Messrohr angeordneten Staukörper und [0032] - einen Sensor zur Erfassung von im Bereich des Staukörpers [0033] auftretenden Druckschwankungen und zur Umwandlung dieser [0034] Druckschwankungen in ein elektrisches Messsignal aufweist, [0035] - einem Analog-Digital Wandler, der dazu dient, zumindest ein Teil des [0036] Messsignals abzutasten und zu digitalisieren, und [0037] - einer Signalverarbeitung, die dazu dient, eine Autokorrelation des [0038] digitalisierten Messsignals zu berechnen und den Durchfluss [0039] anhand mindestens einer Eigenschaft der Autokorrelation abzuleiten.

[0040] Ein Vorteil der Erfindung besteht darin, dass durch die Autokorrelation nicht nur einzelne Messpunkte des Messsignals eingehen, sondern Verlauf und Amplituden des

Messsignals mit berücksichtigt werden. Hierdurch wird eine höhere Messgenauigkeit erzielt. [0041] Ein weiterer Vorteil besteht darin, das dem Messsignal überlagertes weißes

Rauschen sich in der Autokorrelation nur zur Korrelationszeit Null auswirkt, und damit sehr leicht ausgeblendet werden kann.

[0042] Die Erfindung und weitere Vorteile werden nun anhand der Figuren der Zeichnung, in denen fünf Ausführungsbeispiele dargestellt sind, näher erläutert; gleiche Teile sind in den Figuren mit gleichen Bezugszeichen versehen.

[0043] Fig. 1 zeigt ein Wirbeldurchflussmessgerät;

[0044] Fig. 2 zeigt einen Messaufnehmer eines Wirbeldurchflussmessgeräts;

[0045] Fig. 3 zeigt schematisch, eine Ausbildung von Wirbeln hinter einem

[0046] Staukörper;

[0047] Fig. 4 zeigt ein Messsignal als Funktion der Zeit;

[0048] Fig. 5 zeigt eine Autokorrelation des Messsignals von Fig. 4 als Funktion der

[0049] Korrelationszeit;

[0050] Fig. 6 zeigt eine Autokorrelation eines sinusförmigen Messsignals dem

[0051] ausschließlich weißes Rauschen überlagert;

[0052] Fig. 7 zeigt einen ein Minimum enthaltenden Abschnitt der Autokorrelation

[0053] und eine an Punkte in der Umgebung des Minimums angepasste

[0054] Parabel;

[0055] Fig. 8 zeigt eine an einen Sensor des Messaufnehmers angeschlossene

[0056] Schaltung mit einem adaptiven Filter; und

[0057] Fig. 9 zeigt eine an einen Sensor des Messaufnehmers angeschlossene

[0058] Schaltung mit einem Notch^ilter.

[0059] Fig. 1 zeigt ein Beispiel eines Wirbeldurchflussmessgeräts, und Fig. 2 ein Beispiel eines zugehörigen Messaufnehmers. Dieser arbeitet nach dem Prinzip der Karman' sehen Wirbelstraße. Ein Medium, dessen Volumendurchfluss gemessen werden soll, wird durch ein Messrohr 1 geleitet. In einem in Fig. 1 durch einen Pfeil dargestellten Strömungspfad des Mediums ist ein Staukörper 3 eingebracht. Hinter dem Staukörper 3 bilden sich abwechselnd beidseitig Wirbel mit entgegen gesetztem Drehsinn. Dies ist in Fig. 3 schematisch dargestellt. Die Wirbel erzeugen jeweils einen lokalen Unterdruck. Die Druckschwankungen werden von dem Messaufnehmer erfasst und in elektrische Messsignale umgewandelt. Die Wirbel bilden sich innerhalb der zulässigen Einsatzgrenzen des Messaufnehmers sehr regelmäßig aus. Die Anzahl der pro Zeiteinheit erzeugten Wirbel verhält sich proportional zum Volumendurchfluss. Die Frequenz f der Messsignale ist proportional zur Fließgeschwindigkeit des Mediums im Messrohr 1, die wiederum proportional zum Volumendurchfluss ist.

[0060] Es ist ein Sensor 5 zur Erfassung von im Bereich des Staukörpers 3, z.B. am oder

hinter dem Staukörper 3, auftretenden Druckschwankungen und zur Umwandlung dieser Druckschwankungen in ein elektrisches Messsignal U _s vorgesehen. Hierzu eignen sich Drucksensoren, wie z.B. die eingangs erwähnten kapazitiven Sensoren. Der Sensor 5 ist in Strömungsrichtung hinter dem Staukörper 3 im Messrohr 1 angeordnet. Die vom Sensor 5 abgeleiteten elektrischen Messsignale U _s weisen im Healfall einen sinusförmigen Verlauf auf. Fig. 4 zeigt ein solches Messsignal U _s als Funktion der Zeit t. Die Frequenz f des Messsignals U _s ist proportional zur Fließgeschwindigkeit des Mediums im Messrohr 1 und damit zum Durchfluss.

[0061] Erfindungsgemäß wird zumindest ein Teil des Messsignals U _s abgetastet und digitalisiert. Dabei genügt es einen Signalzug, der nur eine geringe Anzahl Schwingungen, z.B. 10 Schwingungen, enthält, mit hinreichender Abtastrate, z.B. 10 Punkte pro Schwingung, abzutasten und zu digitalisieren. Die einzelnen Abtastpunkte sind in Fig. 4 eingezeichnet und stehen im Anschluss an die Digitalisierung als Messpunkte [UsCt ₁ ); tj zur Verfügung, wobei U ₈ (I ₁ ) der zum Abtastzeitpunkt ζ aufgenommene Wert des Messsignals U _s und i ein ganzzahliger laufender Index mit i = 1 .. M ist. Zwischen zwei aufeinander folgenden Abtastzeitpunkten ζ und t ₁₊₁ liegt jeweils ein durch die Abtastrate vorgegebener identischer Zeitraum δt.

[0062] Erfindungsgemäß wird eine Autokorrelation AK(T) als Funktion der

Korrelationszeit T des digitalisierten Messsignals U _s (t) berechnet, und der Durchfluss anhand einer Eigenschaft der Autokorrelation AK(T) abgeleitet.

[0063] Die Autokorrelation AK(T) kann beispielsweise gemäß der folgenden Berechnungsvorschrift:

[0064]

[0065] AK(T) = ∑ ^N U ₈ (I ₁ ) U ₈ (I ₁ + T)

[0066] i =1

[0067] bestimmt werden. Die Messpunkte U ₈ (I ₁ ) wurden zu diskreten Zeitpunkten ζ abgetastet. Entsprechend lässt sich die Autokorrelation AK(T) für diskrete Korrelationszeiten T _k berechnen, mit T _k = k δt wobei k= 0, 1, 2 ... K ist und die Bedingung N + K < M erfüllt ist.

[0068] Die Autokorrelationsfunktion eines periodischen Signals ist ebenfalls periodisch und weist die gleiche Frequenz f und die gleiche Periode P auf, wie das Ausgangssignal, fct das Ausgangssignal sinusförmig, z.B. in der Form U(t) = U ₀ sin (2 π f t), so ist die zugehörige Autokorrelationsfunktion eine Kosinusfunktion AK u(T) = A ₀ cos (2 π f T) mit der gleichen Frequenz f.

[0069] Die Autokorrelationsfunktion eines unkorrelierten Signals S, z.B. von weißem

Rauschen weist lediglich für die Korrelationszeit T= 0 einen von Null verschiedenen Wert auf. Für alle anderen Korrelationszeiten T ≠ 0 ist die zugehörige Autokorrelation AKs(T) = 0.

[0070] Im Healfall des in Fig. 4 dargestellten sinusförmigen Messsignals U _s (t) gilt somit:

[0071]

[0072] AK(T) = E ^N U ₈ Ct ₁ ) UsCt ₁ + T) = C cos(2π f T)

[0073] i =1

[0074] wobei C eine Konstante und

[0075] f gleich der Frequenz des Messsignals U _s (t) ist.

[0076] Die zugehörige Autokorrelationsfunktion AK(T) ist in Fig. 5 dargestellt.

[0077] Besteht das Messsignal aus einem sinusförmigen Signal, dem ausschließlich weißes Rauschen überlagert ist, so gilt die zuletzt genannte Beziehung nur für Korrelationszeiten T ≠ 0. Bei T = O weist die Autokorrelation AK(T) einen Wert C ₀ auf der größer als C ist. Eine entsprechende Autokorrelationsfunktion ist in Fig. 6 dargestellt. Weißes Rauschen wirkt sich in der Autokorrelation nur auf den Wert bei T=O aus. Alle anderen Bereiche der Autokorrelation sind rauschfrei.

[0078] Es wird eine Eigenschaft der Autokorrelation AK(T) bestimmt und der Durchfluss anhand dieser Eigenschaft abgeleitet. Vorzugsweise ist die Eigenschaft der Autokorrelation AK(T) eine Frequenz f oder eine Periode P der Autokorrelation AK(T). Wie oben erläutert, ist die Frequenz f der Autokorrelation, bzw. deren Periode P im wesentlichen gleich der Frequenz f bzw. der Periode P des Messsignals U _s . Letztere ist ein Maß für den Durchfluss und kann, wie oben beschrieben zur Bestimmung des Durchflusses eingesetzt werden. Vorzugsweise wird dabei die Frequenz f oder die Periode P der Autokorrelation AK(T) gleich der Frequenz f oder der Periode P des Messsignals U _s gesetzt. Dies bietet den Vorteil, dass sich Frequenz f oder Periode P des Messsignals U _s anhand der Autokorrelation AK(T) sehr viel genauer bestimmen lassen, als dies durch eine Auswertung des ursprünglichen Messsignals U _s> z.B. durch eine Erfassung der Nulldurchgänge des Messsignals U _s möglich wäre. Bei einer Erfassung der Nulldurchgänge des Messsignals U _s kann sich jede Störung auf die Anzahl der Nulldurchgänge auswirken. Demgegenüber wirkt sich bei der Autokorrelation weißes Rauschen nur lokal bei der Korrelationszeit T = O aus. Ein weiterer Vorteil der Autokorrelation AK(T) besteht darin, dass deren Genauigkeit nicht von der Genauigkeit eines einzelnen Messpunktes abhängt, sondern Messpunkte eines Signalzuges eingehen, und damit Informationen über Amplitude und Verlauf des genutzten Signalzuges ausgenutzt werden. Entsprechend ist die Autokorrelation AK(T)

zur Bestimmung der Frequenz f bzw. der Periode P des Messsignals U _s sehr viel besser geeignet als das Messsignal U _s selbst.

[0079] Unter der Annahme, dass das Messsignal U _s (t) ein sinusförmiges Signal ist, dem kein oder ausschließlich weißes Rauschen überlagert ist, könnte die Frequenz f der Autokorrelationsfunktion AK(T) und damit die Frequenz des Sensorsignals bereits anhand eines einzigen Wertes der Autokorrelation für T ≠ 0 berechnet werden. Diese Vorgehensweise hat jedoch bei realen Messsignalen keine hochgenauen reproduzierbaren Ergebnisse geliefert. Eine Ursache hierfür liegt darin, dass reale Messsignale neben einem sinusförmigen Signal und einem weißen Rauschen zusätzliche Störsignale enthalten. Erfindungsgemäß wird daher eine Eigenschaft des der Autokorrelation AK(T) zur Bestimmung des Durchflusses herangezogen. Dabei genügt es, die Berechnung der im wesentlichen kosinusförmigen Autokorrelation auf markante Abschnitte derselben zu begrenzen. Markante Abschnitte liegen beispielsweise im Bereich von Nulldurchgängen, sowie von Minima oder Maxima der Autokorrelation .

[0080] Gemäß einer ersten Variante wird der Verlauf der Autokorrelation AK(T) im

Bereich des ersten Nulldurchgangs herangezogen. Hierzu wird die Autokorrelation AK(T) beispielsweise solange für ansteigende Korrelationszeiten T, mit T = k δt mit k = 0, 1, 2 .. berechnet, bis mindestens ein erster negativer Wert auftritt, d.h. AK (x δt ) < 0, wobei x eine natürliche Zahl zwischen 1 und K ist. Anhand dieser Werte wird eine erste Nullstelle NSl der Autokorrelation AK(T) bestimmt. Die Lage der ersten Nullstelle NS 1 kann beispielsweise durch eine einfache Interpolation ermittelt werden, bei der durch den Punkt (AK(x δt ) < 0; x δt ) und den vorangehenden Punkt (AK((x-l) δt ) > 0; (x-1) δt ), bei dem die Autokorrelation einen positiven Wert aufweist, eine Gerade gelegt wird, deren Nulldurchgang dann gleich der ersten Nullstelle NS 1 gesetzt wird. Die erste Nullstelle NS 1 ergibt sich hiernach gemäß folgender Berechnungsvorschrift:

[0081] NS 1 = ( x - AK (x δt)/[AK (x δt)- AK((x- 1) δt]) δt

[0082] Die erste Nullstelle NS 1 entspricht einer Viertelperiode der im wesentlichen kosinusförmigen Autokorrelation AK(T). Eine volle Periode P der Autokorrelation AK(T) hat dementsprechend eine Dauer die gleich dem vierfachen der ersten Nullstelle NS 1 ist. Es gilt: P = 4 NS 1. Dementsprechend ergibt sich eine Frequenz f für die Autokorrelation AK(T) mit f = 1/ (4 NSl). Diese Frequenz f , gibt die Frequenz f des Messsignals U _s (t) sehr genau wieder. Die Frequenz f des Messsignals U _s (t) wird folglich gleich der Frequenz f der Autokorrelation AK(T) gesetzt, die anhand der Lage

der ersten Nullstelle NSl bestimmt wurde. Das gleiche gilt entsprechend für die Periode P.

[0083] Der gesuchte Durchfluss ist proportional zur Frequenz f des Messsignals U _s und lässt sich somit bereits anhand der ersten Nullstelle NS 1 der Autokorrelation AK(T) bestimmen. Analog kann der Durchfluss natürlich anhand der zugehörigen Periode P bestimmt werden.

[0084] Ein Vorteil dieses Verfahrens besteht darin, dass nur eine sehr geringe Anzahl der Werte der Autokorrelation benötigt wird und die erforderlichen Rechenoperation sehr einfache Operationen sind, die nur wenige Rechenschritte benötigen. Entsprechend wird für dieses Verfahren nur sehr geringer Speicherplatz und eine geringe Rechenleistung benötigt.

[0085] Gemäß einer zweiten Variante wird der Verlauf der Autokorrelation AK(T) im

Bereich eines Minimums, vorzugsweise des ersten Minimums Ml herangezogen. Das erste Minimum Ml kann auf verschiedene Weise bestimmt werden. Beispielsweise wird die Autokorrelation AK(T) solange für ansteigende Werte von T, mit T = k δt mit k = 0, 1, 2 .. berechnet, bis mindestens ein Wert AK(y δt) der Autokorrelation AK(T) ermittelt wurde, der größer als ein unmittelbar vorangegangener Wert AK((y-l) δt) ist. Hieraus ergibt sich die ungefähre Lage des Minimums Ml. Die genaue Position des Minimums Ml, d.h. die Korrelationszeit T _Min , bei der das Minimum Ml auftritt, kann dann anhand einer Ableitung der Autokorrelation AK(T) in diesem Bereich bestimmt werden. Hierbei wird ausgenutzt, das eine Ableitung einer Funktion, dort, wo die Funktion ein Minimum aufweist, gleich null ist. Es kann folglich eine Nullstelle der Ableitung bestimmt und gleich dem gesuchten ersten Minimum Ml gesetzt werden.

[0086] Alternativ kann, z.B. auf die oben beschriebene Weise, die erste Nullstelle NSl bestimmt werden und die Lage des Minimums Ml anhand der Lage der ersten Nullstelle NSl abgeschätzt werden. Dabei wird ausgenutzt, dass das Minimum Ml bei einer Korrelationszeit T _Mm zu erwarten ist, die gleich dem doppelten der Korrelationszeit der ersten Nullstelle NSl ist. Zusätzlich wird vorzugsweise eine Ableitung der Autokorrelation in dem Bereich berechnet, in dem das Minimum Ml zu erwarten ist und ein Nulldurchgang der Ableitung ermittelt. Die Nullstelle der Ableitung entspricht der Position des ersten Minimums Ml.

[0087] Vorzugsweise wird die genaue Lage des Minimums Ml durch eine Anpassung einer Parabel p(T) an Punkte der Autokorrelation im Bereich des Minimums bestimmt. Dies kann, wie in Fig. 7 dargestellt, anhand von mindestens drei Punkten der Autokorrelation AK(T) erfolgen. Als Punkte eignen sich z.B. die oben angeführten

[AK(yδt); yδt], [AK((y-l)δt); (y-l)δt] und der diesen vorangehende Punkt [AK((y-2)δt); (y-2)δt]. Die Parabel p(T) ist in Fig. 7 als gestrichelte Linie eingezeichnet.

[0088] Ebenso kann der Nulldurchgang der Ableitung der Autokorrelation AK(T) bei der Auswahl der Punkte herangezogen werden. Dabei wird beispielsweise derjenige Punkt der Autokorrelation, der den geringsten Abstand zu der Nullstelle der Ableitung aufweist, sowie die beiden rechts und links davon liegenden unmittelbar dazu benachbarten Punkte verwendet.

[0089] An diese drei Punkte wird, wie in Fig. 7 dargestellt eine Parabel angepasst. Hierzu kann der folgende Ansatz verwendet werden:

[0090] p(T) = P ₀ (T - T _M111 ) ² + _Pl

[0091] wobei p ₀ , Pi und T _Mm Koeffizienten sind, die anhand der drei Punkte der

Autokorrelation AK(T) bestimmt werden. Die Korrelationszeit T _Min entspricht der Position des Minimums der Parabel p(T) und gibt die Korrelationszeit des Minimums der Autokorrelation AK(T) sehr genau wieder. Entsprechend wird sie gleich der Korrelationszeit T _Min des Minimums Ml des Autokorrelation AK(T) gesetzt. Hierdurch ist das Minimum auch dann sehr genau bestimmbar, wenn nur wenige Punkte der Autokorrelation zur Verfügung stehen.

[0092] Die Parabel kann natürlich auch an mehr als drei Punkte der Autokorrelation angepasst werden, z.B. anhand der Methode der Minimierung der Abstandsquadrate. Dies kann insb. einen Genauigkeitsgewinn bringen, wenn das Messsignal mit einer hohen Abtastrate digitalisiert wird und somit deutlich mehr als die oben angegebenen 10 Punkte pro Schwingung vorliegen.

[0093] Die Korrelationszeit T _Min des ersten Minimums Ml entspricht einer halben Periode Vi P. Die Frequenz f der Autokorrelation AK(T) und damit des Messsignals U _s ergibt sich somit zu f = 1/(2 T _Min ).

[0094] Beide oben beschriebenen Verfahren kommen mit einer sehr geringen Anzahl von Punkten der Autokorrelation aus. Die Anzahl lässt sich sogar noch weiter verringern, wenn vorab ein Schätzwert für die Frequenz f, bzw. für die Periode P ermittelt wird. Ein solcher Schätzwert kann z.B. ein messgerätspezifischer Wert sein oder aber ein in einer vorangegangenen Messung ermittelter Wert für eine dieser Größen sein. Liegt ein solcher Schätzwert vor, so kann der Bereich, in dem die erste Nullstelle NSl bzw. das erste Minimum Ml zu erwarten ist eingegrenzt werden und die Berechnung der Autokorrelation auf diese Bereiche begrenzt werden.

[0095] Die Genauigkeit, mit der der Durchfluss bestimmbar ist, lässt sich weiter erhöhen,

indem eine Kosinus Funktion K(T) = cos (2 π f _c T) an die Autokorrelation AK(T) angepasst wird. Hierzu wird vorzugsweise mindestens eine volle Periode der Autokorrelation AK(T) berechnet und es wird die erhaltene Autokorrelation normiert. Vorzugsweise wird hierbei der Wert der Autokorrelation bei der Autokorrelationszeit T= 0 ausgeklammert, da sich in diesem Wert das weiße Rauschen niederschlägt. Die Anpassung kann z.B. erfolgen, indem eine Summe der Abstandsquadrate J zwischen den Werten der Kosinus Funktion und den zugehörigen normierten Werten ak(T J der Autokorrelation minimiert werden, wobei die Frequenz f _c der Kosinusfunktion als Anpassungsparameter dient.

[0096] Die Summe der Abstandsquadrate kann beispielsweise gemäß folgender Berechnungsvorschrift bestimmt werden:

[0097] J(f _c ) =∑ ^L (ak(l δt) - cos (2 π f _c 1 δt)) ²

[0098] 1 =1

[0099] worin 1 ein von 1 bis L laufender ganzzahliger Index,

[0100] ak(l δt) die normierten Werte der Autokorrelation, und

[0101] f _c der Anpassungsparamenter sind.

[0102] Anschließend wird diejenige Frequenz f _c , bei der eine optimale Anpassung vorliegt, gleich der Frequenz f der Autokorrelation AK(T) und gleich der Frequenz f des Messsignals U _s gesetzt, anhand derer dann aufgrund der bestehenden Proportionalitäten der Durchfluss ermittelt wird. Die optimale Anpassung entspricht in dem genannten Beispiel dem Minimum der Summe der Abstandsquadrate J(f _c ).

[0103] In realen Messsituationen kann das Messsignal neben dem im wesentlichen sinusförmigen gesuchten Signal und einem eventuell vorhandenem weißen Rauschen zusätzliche Störsignale, z.B. nicht stationäres Rauschen, überlagert sein, und es kann Amplitudenschwankungen aufweisen. Auch in diesen Situationen zeigen die beschriebenen Verfahren sehr genaue Messergebnisse.

[0104] Die Messgenauigkeit lässt sich noch weiter verbessern, indem das Messsignal U _s vor der Bildung der Autokorrelation einen Filter durchläuft und die Autokorrelation AK(T) anhand des gefilterten Messsignals bestimmt wird.

[0105] Fig. 8 zeigt eine entsprechende Schaltung, die den Sensor 5 und einen daran angeschlossenen Filter 7 aufweist. Die gefilterten Messsignale sind einem Analog- Digital- Wandler 9 zugeführt, der die digitalen gefilterten Messsignale einer Signalverarbeitungseinheit 11, z.B. einem Mikroprozessor, zuführt. Die Signalverarbeitungseinheit 11 bestimmt auf die oben beschriebene Weise die Autokorrelation und deren Eigenschaft, anhand derer nachfolgend der Durchfluss

abgeleitet wird. Dies kann ebenfalls mittels der Signalverarbeitungseinheit 11 erfolgen.

[0106] Der Filter 7, den das Messsignal U _s vor der Bildung der Autokorrelation AK(T) durchläuft ist z.B. ein adaptiver Bandpassfilter. Vorzugsweise wird die Frequenz f der Autokorrelation AK(T) bestimmt und ein Frequenzbereich, in dem der Bandpassfilter durchlässig ist, wird anhand der Frequenz f der Autokorrelation eingestellt. Hierbei wird ausgenutzt, dass die Frequenz f der Autokorrelation AK(T) im wesentlichen gleich der Frequenz f des Messsignals U _s ist. Der Filter 7 wird dabei kontinuierlich an die Frequenz f des Messsignals U _s angepasst, indem der nur für einen Frequenzbereich in der unmittelbaren Umgebung der Frequenz f des Messsignals U _s durchlässig gestellt wird. Entsprechend werden im Messsignal enthaltene Störsignale anderer Frequenzen unmittelbar herausgefiltert.

[0107] Alternativ kann der Durchlassbereich des Bandpassfilters anhand einer unmittelbar anhand des Messsignals abgeleiteten Frequenz des Messsignals eingestellt werden. Hierzu können beispielsweise auf herkömmliche Weise die Nulldurchgänge des Messsignals gezählt und anhand deren Anzahl pro Zeiteinheit, die Frequenz bestimmt werden.

[0108] Anstelle von oder zusätzlich zu dem dargestellten analogen Filter 7 kann auch ein digitaler Filter 13 eingesetzt werden. Dieser ist dann, wie in Fig. 7 gestrichelt eingezeichnet, hinter dem Analog Digital Filter 9 anzuordnen. Hier eignen sich z.B. digitale adaptive Bandpassfilter, die analog zu dem oben im Zusammenhang mit dem analogen adaptiven Bandpassfilter Gesagten eingesetzt werden.

[0109] Ein weiteres bevorzugtes Ausführungsbeispiel für einen digitalen Filter 13 sind so genannte Adaptive Line Enhancer (ALE). Derartige Filter stellen sich automatisch auf die Frequenz des Messsignals ein und bewirken eine Dämpfung des Untergrundrauschens .

[0110] Alternativ oder zusätzlich kann, insb. bei hohen Abtastraten, ein Mittelwertfilter 15 eingesetzt werden. Dies ist in Fig. 8 ebenfalls dargestellt. Mittelwertfilter bewirken eine Glättung des Messsignals, indem sie eine Mittelung über mehrere aufeinander folgende Messwerte ausführen.

[0111] Fig. 9 zeigt eine weitere Schaltung, mittels der das Messsignal vor der Bildung der Autokorrelation einer Filterung unterzogen wird. Als Filter ist hier ein adaptiver Notch Filter 17 vorgesehen. Notch Filter sind hochselektive Filter für einen schmalen Sperrbereich. Ein durch eine Filterfunktion H(z) beschriebener Notch Filter 17 liefert an dessen Ausgang ein Restsignal R, das alle außerhalb des Sperrbereichs liegenden Anteile des Eingangssignals enthält. üblicherweise wird dieses Restsignal R zur

Optimierung des Notch Filters 17 eingesetzt. Die Optimierung des Notch Filters 17 lässt sich z.B. wie in Fig. 9 schematisch dargestellt durch eine Anpassung des Filters 17 dahingehend bewirken, das das Restsignal R minimiert wird.

[0112] Erfindungsgemäß führt der Notch Filter 17, zusätzlich die zu der Filterfunktion H(z) komplementäre Filterfunktion 1 - H(z) aus, die aus dem eingehenden Sensorsignal das eigentliche Messsignal U _s herausgefiltert. Dieses herausgefilterte Messsignal U _s wird dann, wie oben beschrieben weiter verarbeitet, indem es der Signalverarbeitungseinheit 11 zugeführt wird, dessen Autokorrelation AK(T) bestimmt und anhand mindestens einer Eigenschaft der Autokorrelation AK(T) der Durchfluss bestimmt wird.

[0113]

Tabelle 1

1 Messrohr

3 Staukörper

5 Sensor

7 Filter

9 Analog Digital Wandler

11 Signalverarbeitungseinheit

13 digitaler Filter

15 Mittelwertfilter

17 Notch Filter

[0114]