PL197438B1 - Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa - Google Patents

Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa

Info

Publication number
PL197438B1
PL197438B1 PL352214A PL35221400A PL197438B1 PL 197438 B1 PL197438 B1 PL 197438B1 PL 352214 A PL352214 A PL 352214A PL 35221400 A PL35221400 A PL 35221400A PL 197438 B1 PL197438 B1 PL 197438B1
Authority
PL
Poland
Prior art keywords
signals
frequency
signal
flow tube
samples
Prior art date
Application number
PL352214A
Other languages
English (en)
Other versions
PL352214A1 (en
Inventor
Denis Henrot
Original Assignee
Micro Motion
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=23352280&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=PL197438(B1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by Micro Motion filed Critical Micro Motion
Publication of PL352214A1 publication Critical patent/PL352214A1/xx
Publication of PL197438B1 publication Critical patent/PL197438B1/pl

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F15/00Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
    • G01F15/02Compensating or correcting for variations in pressure, density or temperature
    • G01F15/022Compensating or correcting for variations in pressure, density or temperature using electrical means
    • G01F15/024Compensating or correcting for variations in pressure, density or temperature using electrical means involving digital counting
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78Direct mass flowmeters
    • G01F1/80Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84Coriolis or gyroscopic mass flowmeters
    • G01F1/8409Coriolis or gyroscopic mass flowmeters constructional details
    • G01F1/8413Coriolis or gyroscopic mass flowmeters constructional details means for influencing the flowmeter's motional or vibrational behaviour, e.g., conduit support or fixing means, or conduit attachments
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78Direct mass flowmeters
    • G01F1/80Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84Coriolis or gyroscopic mass flowmeters
    • G01F1/8409Coriolis or gyroscopic mass flowmeters constructional details
    • G01F1/8436Coriolis or gyroscopic mass flowmeters constructional details signal processing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78Direct mass flowmeters
    • G01F1/80Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84Coriolis or gyroscopic mass flowmeters
    • G01F1/845Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits
    • G01F1/8468Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits
    • G01F1/8472Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits having curved measuring conduits, i.e. whereby the measuring conduits' curved center line lies within a plane
    • G01F1/8477Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits having curved measuring conduits, i.e. whereby the measuring conduits' curved center line lies within a plane with multiple measuring conduits
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Measuring Volume Flow (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

1. Sposób okre slania w la sciwo sci materia lu przep lywaj a- cego przez rur e przep lywow a przep lywomierza typu Corioli- sa, polegaj acy na tym, ze wprawia si e przewód w drgania i mierzy si e drgania przewodu za pomoc a pierwszego czujni- ka przesuni ecia i drugiego czujnika przesuni ecia jednocze- snie wytwarzaj ac sygna l cyfrowy, a nast epnie dostarcza si e próbki sygna lów otrzymywanych z pierwszego czujnika prze- suni ecia i drugiego czujnika przesuni ecia z pierwsz a pr edko- sci a próbkowania do cyfrowego procesora sygna lowego, po czym zmienia si e pr edko sc próbkowania z pierwszej pr edko- sci próbkowania do po sredniej pr edko sci próbkowania za po- moc a pierwszej decymacji tych próbek, znamienny tym, ze dokonuje si e demodulacji (403) próbek o po sredniej pr ed- ko sci próbkowania i zmienia si e pr edko sc próbkowania demo- dulowanych próbek do wymaganej pr edko sci próbkowania poprzez przeprowadzenie drugiej decymacji (404) próbek sygna lów, a nast epnie z tych próbek decymowanych sygna- lów wyznacza si e cz estotliwo sc (405) drga n rury przep lywo- wej (103A, 103B) przy pierwszym czujniku przesuni ecia (105) i przy drugim czujniku przesuni ecia (105'), za s na podstawie tych cz estotliwo sci wyznacza si e w la sciwo sci materia lu prze- p lywaj acego przez rur e przep lywow a (103A, 103B). PL PL PL PL

Description

(12) OPIS PATENTOWY (19) PL (21) Numer zgłoszema: 352214 (22) Date zgtoszema: 12·06·2000 (86) Data i numer zgłoszenia międzynarodowego:
12.06.2000, PCT/US00/16061 (87) Data i numer publikacji zgłoszenia międzynarodowego:
04.01.2001, WO01/01083 PCT Gazette nr 01/01 (11) 197438 113) B1 (51 ) Int.Cl.
G01F 1/84 (2006.01)
Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa
(30) Pierwszeństwo: 28.06.1999,US,09/344,840 (73) Uprawniony z patentu: MICRO MOTION, INC.,Boulder,US
(43) Zgłoszenie ogłoszono: 11.08.2003 BUP 16/03 (72) Twórca(y) wynalazku: Denis Henrot,Louisville,US
(45) O udzieleniu patentu ogłoszono: 31.03.2008 WUP 03/08 (74) Pełnomocnik: Słomińska-Dziubek Anna, POLSERVICE, Kancelaria Rzeczników Patentowych Sp. z o.o.
(57) 1. Sposób określania właściwości materiałuprzepływającego przez rurę przepływową przepływomierza typu Coriolisa, polegający na tym, że wprawia się przewód w drgania i mierzy się drgania przewodu za pomocą pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia jednocześnie wytwarzając sygnał cyfrowy, a następnie dostarcza się próbki sygnałów otrzymywanych z pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia z pierwszą prędkością próbkowania do cyfrowego procesora sygnałowego, po czym zmienia się prędkość próbkowania z pierwszej prędkości próbkowania do pośredniej prędkości próbkowania za pomocą pierwszej decymacji tych próbek, znamienny tym, że dokonuje się demodulacji (403) próbek o pośredniej prędkości próbkowania i zmienia się prędkość próbkowania demodulowanych próbek do wymaganej prędkości próbkowania poprzez przeprowadzenie drugiej decymacji (404) próbek sygnałów, a następnie z tych próbek decymowanych sygnałów wyznacza się częstotliwość (405) drgań rury przepływowej (103A, 103B) przy pierwszym czujniku przesunięcia (105) i przy drugim czujniku przesunięcia (105'), zaś na podstawie tych częstotliwości wyznacza się właściwości materiału przepływającego przez rurę przepływową (103A, 103B).
PL 197 438 B1
Opis wynalazku
Przedmiotem wynalazku jest sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa.
Znane jest stosowanie przepływomierzy masowych wykorzystujących zjawisko Coriolisa do pomiaru przepływu masowego i innych informacji o materiałach przepływających przez kanał. Przykładowe takie przepływomierze opisano w patentach USA o numerach 4,109,534 z 29 sierpnia 1978, 4,491,025 z 1 stycznia 1985, i Re 31,450 z 11 lutego 1982, wszystkie na nazwisko J. E. Smith i in. Te przepływomierze są wyposażone w jedną lub więcej rur przepływowych o kształcie prostym lub zakrzywionym. Każdy kształt rury przepływowej w przepływomierzu ma pewien zespół naturalnych modów drgań, które mogą być typu prostego zginania, skręcania lub sprzężone. Każda rura przepływowa jest pobudzana do drgań w rezonansie z jednym z tych naturalnych modów. Materiał przepływa do wnętrza miernika przepływu z dołączonego przewodu po stronie dopływowej przepływomierza, jest kierowany przez jedną lub więcej rur przepływowych, i opuszcza przepływomierz przez stronę odpływową. Naturalne mody drgań układu wypełnionego materiałem drgającym są częściowo wyznaczone przez łączną masę rur przepływowych i materiału płynącego wewnątrz rur przepływowych.
Przy braku przepływu przez przepływomierz, wszystkie punkty wzdłuż rury przepływowej oscylują z identyczną fazą pod wpływem przyłożonej siły napędowej, a ewentualny niewielki wstępny stały offset fazowy może zostać skorygowany. W miarę, jak rozpoczyna się przepływ płynu, siły Coriolisa powodują, że każdy punkt wzdłuż rury przepływowej ma inną fazę. Faza po stronie dopływowej rury przepływowej opóźnia się względem wzbudnika, natomiast faza po stronie odpływowej wyprzedza fazę wzbudnika. Na rurze przepływowej rozmieszczone są czujniki przesunięcia wytwarzające sygnały sinusoidalne reprezentujące ruch rury przepływowej. Sygnały wyprowadzane z czujników przesunięcia są przetwarzane w celu wyznaczenia różnicy fazowej między czujnikami przesunięcia. Różnica faz między dwoma sygnałami czujników jest proporcjonalna do strumienia masowego przepływu materiału w rurze.
Przepływomierz typu Coriolisa zaopatrzony jest w nadajnik generujący sygnał wzbudzający, napędzający wzbudnik i wyznaczający strumień przepływu masowego i inne właściwości materiału na podstawie sygnałów otrzymanych z czujników przesunięcia. Konwencjonalny nadajnik jest wykonany z układów analogowych, które są przeznaczone do generowania sygnału wzbudzającego i do detekcji sygnałów z czujników przesunięcia. Nadajniki analogowe zostały zoptymalizowane w ciągu lat i stały się stosunkowo tanimi w wytwarzaniu. Pożądane jest, zatem, opracowanie przepływomierza typu Coriolisa wykorzystującego nadajnik konwencjonalny.
Pewnym problemem jest to, że nadajnik konwencjonalny musi pracować z sygnałami w wąskim paśmie częstotliwości roboczych. Ten zakres częstotliwości roboczych zawiera się zwykle między 20 Hz a 200 kHz. Ogranicza to konstruktorów do tego wąskiego pasma częstotliwości roboczych. Ponadto, wąski zakres częstotliwości roboczych uniemożliwia wykorzystanie konwencjonalnego nadajnika w niektórych przepływomierzach, jak na przykład przepływomierzu z prostą rurą przepływową, która pracuje w wyższym zakresie częstotliwości wynoszącym 300 Hz - 800 Hz. Przepływomierze z prostą rurą przepływową działające w zakresie 300 Hz - 800 Hz mają tendencję do wykazywania mniejszej czułości na efekty Coriolisa wykorzystywane do pomiaru strumienia przepływu masowego. Zatem do obliczania przepływu masowego potrzebne są bardziej precyzyjne pomiary różnicy faz między czujnikami.
W celu wykorzystania jednego typu nadajnika w kilku konstrukcjach przepływomierzy Coriolisa, pracujących na kilku różnych częstotliwościach, wytwórcy przepływomierzy Coriolisa stwierdzają, że pożądane byłoby wykorzystanie cyfrowego procesora sygnałowego do generowania sygnałów wzbudzających, i do przetwarzania sygnałów otrzymywanych z czujników przesunięcia. Cyfrowy procesor sygnałowy jest pożądany, ponieważ przez digitalizację sygnałów z czujników przesunięcia jako sygnałów odbieranych przez nadajnik można uniknąć większych wymagań narzucanych na rozdzielczość i dokładność analogowych części elektronicznych ze strony przepływomierzy działających na wyższych częstotliwościach, na przykład konstrukcji z prostą rurą przepływową. Ponadto, możliwe jest zmodyfikowanie instrukcji do procesów sygnalizacyjnych stosowane w cyfrowym procesorze sygnałowym, odpowiednio do pracy na kilku różnych częstotliwościach.
Jednakowoż cyfrowe procesory sygnałowe mają kilka wad w porównaniu z konwencjonalnymi nadajnikami na układach analogowych. Pierwszy problem polega na tym, że cyfrowy procesor sygnałowy jest kosztowniejszy w produkcji, ponieważ jego układy są bardziej złożone. Po drugie, cyfrowy procesor sygnałowy wymaga płytki montażowej o większej powierzchni, co może stwarzać problemy, kiedy podstawowe znaczenie w konstrukcji przepływomierza ma zajmowana przestrzeń. Po trzecie,
PL 197 438 B1 cyfrowe procesory sygnałowe wymagają większej mocy zasilania, niż układy analogowe. Zużycie mocy stanowi problem zwłaszcza, kiedy procesor musi pracować z maksymalną prędkością zegara dla wykonania wszystkich obliczeń potrzebnych przy przetwarzaniu sygnałów i generowaniu wyników pomiaru właściwości materiału, jak na przykład przepływu masowego. Z tych wszystkich powodów, w technice występuje zapotrzebowanie na cyfrowy procesor sygnałowy nadający się do adaptowania do kilku konstrukcji przepływomierzy, niekosztowny w produkcji i pozwalający zmniejszyć ilość mocy obliczeniowej potrzebnej do wykonywania niezbędnych obliczeń.
W patencie USA nr 5,555,190 opisano urządzenie i sposób wyznaczania zależności częstotliwościowych i fazowych drgających rur przepływowych w przepływomierzu masowym Coriolisa. Do dokładnego określania zależności częstotliwościowych i fazowych rury przepływowej i dzięki temu do dokładniejszego wyznaczenia przepływu masowego materiału płynącego przez przepływomierz masowy w cyfrowym urządzeniu przetwarzania sygnałów wykorzystuje się adaptacyjne metody i urządzenia. Do sygnałów od czujników na drgających rurach przepływowych przyporządkowane filtry decymacyjne zapobiegające zjawisku zawijania widma. Filtry zmniejszają liczbę próbek z analogowo-cyfrowego urządzenia próbkującego pracującego ze stałą częstotliwością, przyporządkowanego do każdego sygnału czujnikowego.
W patencie USA nr 5,734,112 opisano sposób wyznaczania ciśnienia w przepływomierzu masowym działającym na zasadzie zjawiska Coriolisa. Rury przepływowe przepływomierza typu Coriolisa są pobudzane do drgań zarówno w trybie zginania, jak i w trybie skręcania. Stosunek częstotliwości podstawowych, na których rury przepływowe drgają w każdym z tych dwóch trybów drgań jest proporcjonalny do ciśnienia w wewnątrz rury przepływowej. W korzystnej odmianie wykonania metodą tworzenia sumy/różnicy wstępnie wyizolowuje się nałożone wzajemnie sinusoidy reprezentujące częstotliwości podstawowe dwóch trybów drgań. Następnie do szybkiego oszacowania częstotliwości podstawowych każdego z dwóch trybów drgań wykorzystuje się szybkie sprzężone filtry cyfrowe. Szacowane częstotliwości są następnie wykorzystywane przez łańcuchy filtracyjne obejmujące cyfrowe filtry wycinające i filtry pasmowo-przepustowe, jak również rekursywne sposoby cyfrowej filtracji maksymalnego prawdopodobieństwa do uwydatnienia szacowanych wartości podstawowych częstotliwości trybu zginania i trybu skręcania. Uwydatnione wartości podstawowych szacowanych częstotliwości trybu zginania i trybu skręcania wykorzystuje się do wyznaczenia ciśnienia wewnątrz rur przepływowych w funkcji stosunku tych dwóch częstotliwości, jak również do wypośrodkowania łańcuchów filtrów wycinających i filtrów pasmowych wykorzystywanych do uwydatniania wartości częstotliwości trybu zginania obu kanałów czujników drgań, dla obliczeń strumienia przepływu. Tak wyznaczone ciśnienie może być wtedy wykorzystywane do korygowania obliczeń przepływu masowego lub do innych celów związanych z pomiarem ciśnienia jako takim.
Według wynalazku, sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa, polegający na tym, że wprawia się przewód w drgania i mierzy się drgania przewodu za pomocą pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia jednocześnie wytwarzając sygnał cyfrowy, a następnie dostarcza się próbki sygnałów otrzymywanych z pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia z pierwszą prędkością próbkowania do cyfrowego procesora sygnałowego, po czym zmienia się prędkość próbkowania z pierwszej prędkości próbkowania do pośredniej prędkości próbkowania za pomocą pierwszej decymacji tych próbek, charakteryzuje się tym, że dokonuje się demodulacji próbek o pośredniej prędkości próbkowania i zmienia się prędkość próbkowania demodulowanych próbek do wymaganej prędkości próbkowania poprzez przeprowadzenie drugiej decymacji próbek sygnałów, a następnie z tych próbek decymowanych sygnałów wyznacza się częstotliwość drgań rury przepływowej przy pierwszym czujniku przesunięcia i przy drugim czujniku przesunięcia, zaś na podstawie tych częstotliwości wyznacza się właściwości materiału przepływającego przez rurę przepływową.
Korzystnie jako częstotliwości sygnałów wyznacza się znormalizowane częstotliwości sygnałów.
Demodulowania próbek dokonuje się obliczając znormalizowaną pulsację znormalizowanej częstotliwości sygnałów i obliczając iloczyn skalarny pulsacji znormalizowanej i sygnałów pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia, z przeniesieniem sygnałów na częstotliwość środkową.
Znormalizowaną częstotliwość sygnałów wyznacza się poprzez przeprowadzenie demultipleksowania sygnałów na składowe współfazowe i składowe kwadraturowe, a następnie całkuje się składowe współfazowe i składowe kwadraturowe i kolejno dokonuje się multipleksowania składowych współfazowych i składowych kwadraturowych wytwarzając sygnały całkowane cyfrowo, po czym wyzna4
PL 197 438 B1 cza się stosunek amplitudy sygnałów do amplitudy sygnałów całkowanych cyfrowo otrzymując znormalizowane częstotliwości sygnałów.
Po dokonaniu całkowania i przed dokonaniem multipleksowania całkowane składowe kwadraturowe przesyła się do kompensatora.
Po dokonaniu całkowania i przed dokonaniem multipleksowania przesyła się całkowane składowe współfazowych do kompensatora.
Częstotliwość drgań rury przepływowej wyznacza się poprzez odczytanie znormalizowanej częstotliwości sygnałów, a następnie dokonanie modulacji znormalizowanej częstotliwości sygnałów i przeprowadzenie zespolonej demodulacji sygnałów z wykorzystaniem zmodulowanej częstotliwości znormalizowanej.
Wyznaczając częstotliwości drgań rury przepływowej dokonuje się decymacji zdemodulowanych sygnałów i dokonuje się zespolonej korelacji sygnałów z obu czujników przesunięcia wyznaczając różnicę fazową między sygnałami.
Wyznaczając częstotliwości drgań rury przepływowej mierzy się różnicę fazy między sygnałami z pierwszego czujnika przesunięcia, a sygnałami z drugiego czujnika przesunięcia.
Jako właściwość materiału przepływającego przez rurę przepływową wyznacza się jego natężenie przepływu masowego materiału.
Jako właściwość materiału przepływającego przez rurę przepływową wyznacza się jego gęstość.
Według niniejszego wynalazku proponuje się przechowywanie procesów zapamiętanych w pamięci i wykonywanych przez procesor w celu przetwarzania sygnałów odbieranych z przetworników przesunięcia znajdujących się na rurze przepływowej. Procesy według wynalazku dają wiele korzyści, umożliwiające wykorzystywanie jednego typu cyfrowego procesora sygnałowego w wielu typach przepływomierzy Coriolisa.
Pierwszą zaletą procesów według niniejszego wynalazku jest to, że procesy te nie powodują utraty dokładności pomimo wykorzystywania arytmetyki stałoprzecinkowej zamiast zmiennoprzecinkowej. Drugą zaletą procesów według niniejszego wynalazku jest to, że procesy mogą być implementowane na wielu tanich cyfrowych procesorach sygnałowych, o niewielkiej mocy. Instrukcje dla procesów według niniejszego wynalazku są dostatecznie małe, aby mogły rezydować w wewnętrznej pamięci cyfrowego procesora sygnałowego, co eliminuje konieczność stosowania pamięci zewnętrznej o szybkim dostępie, która zwiększa koszt, zużycie mocy i miejsca zajmowane dla nadajnika na płycie. Procesy mają niewielką liczbę struktur obliczeniowych, co poprawia przenośność procesów między niedrogimi prostymi procesorami.
Trzecią zaletą jest to, że wymagania obliczeniowe procesów są zminimalizowane. To zmniejszenie wymagań obliczeniowych umożliwia pracę cyfrowego procesora sygnałowego z częstotliwością cyklu zegarowego mniejszą, niż maksymalna częstotliwość cyklu zegarowego procesora, co zmniejsza zużycie mocy procesora.
Nadajnik, który realizuje procesy według niniejszego wynalazku zawiera wymienione poniżej elementy elektroniczne. Sygnały analogowe z przetworników przesunięcia odbierane są przez przetwornik analogowo-cyfrowy („A/D” - Analog to Digital convertor). Sygnały cyfrowe po konwersji podawane są do standardowego procesora cyfrowego. Procesor cyfrowy jest to jednostka przetwarzająca, która wykonuje czytelne dla maszyny polecenia, które są przechowywane w pamięci dołączonej do procesora za pośrednictwem magistrali. Typowy procesor cyfrowy ma pamięć stałą (ROM), która przechowuje polecenia potrzebne do wykonywanych procesów, na przykład procesów według niniejszego wynalazku. Procesor jest połączony również z pamięcią o dostępie swobodnym, która przechowuje polecenia dla procesu aktualnie realizowanego, i dane potrzebne do realizacji procesu. Procesor może również generować sygnały dla przepływomierza typu Coriolisa. W celu podania sygnału wzbudzającego do układu wzbudzającego, procesor cyfrowy może być połączony za pośrednictwem przetwornika cyfrowo-analogowego (D/A - Digital to Analog), który odbiera sygnały cyfrowe od procesora i podaje sygnały analogowe do układu wzbudzającego.
Procesy według niniejszego wynalazku realizują funkcje wyznaczania częstotliwości sygnałów otrzymywanych z czujników przesunięcia, jak również różnicy faz między sygnałami. Najpierw próbki sygnałów są odbierane z czujników przesunięcia, z pierwszą prędkością próbkowania. Prędkość próbkowania jest to ilość danych wejściowych otrzymywanych z przetworników przesunięcia, która charakteryzuje sygnały z przetworników przesunięcia. Realizowana jest pierwsza decymacja próbek, z pierwszej prędkości próbkowania do pożądanej prędkości próbkowania, decymacja jest to po prostu konwersja z pierwszej liczby próbek na mniejszą liczbę próbek. Decymację wykonuje się dla zwiększenia
PL 197 438 B1 rozróżniania próbkowanych sygnałów z zapewnieniem większej dokładności obliczeń częstotliwości sygnałowej każdego sygnału. Następnie przy pośredniej prędkości próbkowania odbywa się demodulacja. Realizowana jest następnie druga decymacja próbek, do pożądanej prędkości próbkowania. Na podstawie próbek sygnałów, przy pożądanej prędkości próbkowania, wyznaczana jest częstotliwość drgań rur przepływowych przy pierwszym czujniku przesunięcia i przy drugim czujniku przesunięcia.
W celu wykorzystywania tych samych procesów z różnymi przepływomierzami, o różnych częstotliwościach, można realizować również następujące kroki. Oblicza się oszacowanie częstotliwości drgań przepływomierza. Ta oszacowana częstotliwość jest następnie wykorzystywana do demodulowania sygnałów z każdego przetwornika przesunięcia na pewną składową l i pewną składową Q. Składowa l i składowa Q służą następnie do przemieszczania sygnałów do pewnej częstotliwości środkowej, jeżeli częstotliwość robocza sygnałów jest większa od pewnej częstotliwości przejściowej. Po demodulacji, sygnały mogą być drugi raz poddawane decymacji dla powtórnego poprawienia rozdzielczości sygnałów.
Następnie jest wyizolowywana i precyzyjnie mierzona częstotliwość dominująca sygnałów. Następnie wylicza się przesunięcie do częstotliwości zerowej, dla zarówno składowej l, jak i składowej Q sygnałów. Wtedy każda składowa może być ponowna poddana decymacji dla poprawienia dokładności pomiaru. Możliwe jest przy tym, przez zastosowanie odpowiedniej filtracji dolnoprzepustowej, zwężenie w potrzebnym stopniu pasma każdego sygnału. Realizuje się następnie zespoloną korelację, która zapewnia wyznaczenie różnicy fazowej między sygnałami.
Przedmiot wynalazku jest objaśniony w przykładach realizacji w oparciu o rysunek, na którym fig. 1 przedstawia przepływomierz typu Coriolisa zaopatrzony w nadajnik cyfrowy realizujący z wieloma prędkościami procesy przetwarzania sygnału z przetwornika przesunięcia, według niniejszego wynalazku; fig. 2 - schemat blokowy cyfrowego nadajnika sygnału; fig. 3 - sieć działań operacji wykonywanych przez nadajnik cyfrowy; fig. 4 - sieć działań procesu generowania danych na podstawie sygnałów otrzymywanych z czujników przesunięcia; fig. 5 - proces do realizacji decymacji próbek sygnału z przetwornika przesunięcia; fig. 6 - proces obliczania szacowanej częstotliwości sygnałów otrzymywanych z przetworników przesunięcia; fig. 7 - proces do realizacji dwuwartościowego wyboru częstotliwości odbieranych sygnałów; fig. 8 - proces demodulacji odbieranych sygnałów; a fig. 9 - sposób wyznaczania z odbieranych sygnałów danych o wibracjach rury przepływowej.
Przepływomierz typu Coriolisa - fig. 1
Figura 1 przedstawia przepływomierz 5 typu Coriolisa zawierający zespół mierzący 10 Coriolisa i nadajnik 20. Nadajnik 20 jest dołączony do zespołu mierzącego 10 za pośrednictwem przewodów 100 dla podawania ścieżką 26 danych o gęstości, prędkości przepływu masowego, prędkości przepływu objętościowego, temperatury, sumarycznego przepływu masowego i przyrostu gęstości. Opisano konstrukcję przepływomierza typu Coriolisa, jakkolwiek dla specjalisty jest oczywiste, że niniejszy wynalazek może być wykorzystywany w połączeniu z dowolnym urządzeniem zawierającym wibracyjną rurę przepływową do pomiaru właściwości materiału. Drugim przykładem takiego urządzenia jest densytometr z wibrującą rurą przepływową, który nie ma dodatkowych możliwości pomiarowych masowego przepływomierza typu Coriolisa.
Zespół mierzący 10 zawiera parę kołnierzy 101 i 101', rurę rozgałęźna 102 i rury przepływowe 103A i 103B. Wzbudnik 104 i czujniki 105 i 105' przesunięcia są połączone z rurami przepływowymi 103A-B. Paski usztywniające 106 i 106' służą do wyznaczenia osi W i W, wokół której oscyluje każda rura przepływowa.
Po wstawieniu przepływomierza 10 w układ (nie przedstawiony) rurociągu, który przenosi materiał przetwarzany podlegający pomiarowi, materiał wchodzi w zespół mierzący 10 przez kołnierz 101. przechodzi przez rozgałęzienie 102, gdzie jest kierowany do wejścia rur przepływowych 103A i 103B. przepływa przez rury przepływowe 103A i 103B i z powrotem do rozgałęzienia 102, skąd przez kołnierz 101' wychodzi z zespołu mierzącego 10.
Rury przepływowe 103A i 103B są dobrane i odpowiednio zamontowane do rozgałęzienia 102 tak, aby miały w zasadzie taki sam rozkład mas, momentów bezwładności i modułów sprężystości wokół osi zginania, odpowiednio W-W i W-W. Rury przepływowe rozchodzą się na zewnątrz od rozgałęzienia w zasadzie równolegle.
Rury przepływowe 103A - 103B są wzbudzane wzbudnikiem 104 w przeciwnych kierunkach wokół ich odpowiednich osi zginania W i W, co jest nazywane pierwszym niewspółfazowym modem zginania. Wzbudnik 104 może zawierać jedną lub więcej znanych struktur, jak na przykład magnes zamocowany do rury przepływowej 103A ze znajdującą się naprzeciwko niego cewką zamocowaną do
PL 197 438 B1 rury przepływowej 103B, przez którą przepływa prąd przemienny do wzbudzania drgań obu rur przepływowych. Odpowiedni sygnał wzbudzający podawany jest przez elektroniczne układy mierzące 20, za pośrednictwem doprowadzenia 110, do wzbudnika 104.
Nadajnik 20 otrzymuje sygnały prędkościowe, lewy i prawy, pojawiające się na przewodach, odpowiednio 111 i 111'. Nadajnik 20 wytwarza sygnał wzbudzający pojawiający się na doprowadzeniu 110 i powodujący nadawanie przez wzbudnik 104 drgań rurom przepływowym 103A i 103B. Nadajnik 20 przetwarza sygnały prędkościowe, lewy i prawy, dla obliczenia masowego strumienia przepływu i gęstości materiału przechodzącego przez zespół mierzący 10. Ta informacja zostaje podana na ścieżkę 26.
Specjaliście jest wiadome, że przepływomierz 5 Coriolisa jest w budowie bardzo podobny do densytometru z drgającą rurą przepływową. Densytometry z drgającą rurą przepływową również wykorzystują rurę wibrującą, przez którą przepływa płyn, lub w przypadku densytometru z próbką, wewnątrz której utrzymywany jest płyn. Densytometry z rurą drgającą również wykorzystują układ wzbudzający do pobudzania rury przepływowej do drgań. Densytometry z drgającą rurą przepływową zwykle wykorzystują tylko jeden sygnał sprzężenia zwrotnego, ponieważ pomiar gęstości wymaga tylko pomiaru częstotliwości, a pomiar fazy nie jest konieczny. Opisy wynalazku zawarte w niniejszym dokumencie mają zastosowanie również do densytometrów z drgającą rurą przepływową.
Cyfrowy nadajnik 20 - fig. 2
Figura 2 przedstawia części składowe nadajnika cyfrowego 20, ścieżki 111 i 111' transmitują sygnały prędkościowe, lewy i prawy, z zespołu 10 przepływomierza do nadajnika 20. Sygnały prędkościowe są odbierane przez przetwornik analogowo-cyfrowy (A/D) 203 w zespole 20 układów elektronicznych. Przetwornik analogowo-cyfrowy 203 dokonuje konwersji sygnałów prędkościowych, lewego i prawego, nadających się do wykorzystywania przez jednostkę przetwarzającą 201, i transmituje sygnały cyfrowe ścieżką 210 - 210'. Jakkolwiek przedstawiony w postaci oddzielnych części składowych, przetwornik analogowo-cyfrowy 203 może być pojedynczym przetwornikiem cyfrowo-analogowym, jak na przykład układ scalony, który zawiera 2 przetworniki, tak że sygnały od obu przetworników przesunięcia mogą być przetwarzane jednocześnie. Sygnały cyfrowe przenoszone są ścieżkami 210 - 210' do procesora 201. Dla specjalisty jest oczywiste, że do procesora 201 może być dołączona dowolna liczba czujników przesunięcia i innych czujników, jak na przykład oporowy czujnik temperatury RTD, do wyznaczania temperatury rury przepływowej.
Sygnały wzbudnika są transmitowane ścieżką 21.2, która dostarcza sygnały do przetwornika cyfrowo-analogowego (D/A) 202. Przetwornik cyfrowo-analogowy 202 jest zwykłym przetwornikiem cyfrowo-analogowym i może być oddzielnym przetwornikiem cyfrowo-analogowym lub przetwornikiem scalonym w układzie kodeka stereofonicznego. Sygnały analogowe z przetwornika cyfrowo-analogowego 202 są transmitowane do obwodu wzbudzającego 290 za pośrednictwem ścieżki 291. Układ wzbudzający 291 przy tym podaje sygnał wzbudzający do wzbudnika 104 za pośrednictwem ścieżki 110. Ścieżka 26 przenosi sygnały do elementów wejściowych i wyjściowych (nie pokazane), które umożliwiają odbieranie i przekazywanie przez nadajnik 20 danych od, i do operatora.
Jednostką przetwarzającą 201 jest mikroprocesor, procesor lub grupa procesorów, która odczytuje polecenia z pamięci i realizuje polecenia wykonując różne funkcje przepływomierza. Realizowane funkcje obejmują, choć nie tylko, obliczanie przepływu masowego materiału, obliczanie przepływu objętościowego materiału, i obliczanie gęstości materiału, które mogą być przechowywane w postaci poleceń w pamięci stałej (ROM) 220. Procesor 201 realizuje operacje w pamięci ROM 220 za pośrednictwem ścieżki 221. Dane, jak również polecenia do realizacji różnych funkcji są przechowywane w pamięci 230 o dostępie swobodnym (RAM). Operacje odczytu i zapisu w pamięci RAM 230 procesor 201 realizuje za pośrednictwem ścieżki 231.
Przegląd operacji realizowanych przez nadajnik 20 - fig. 3.
Na fig. 3 przedstawiono przegląd funkcji realizowanych przez cyfrowy nadajnik 20 pracujący z przepływomierzem 5. Proces 300 rozpoczyna się w kroku 301 od generowania przez nadajnik 20 sygnału wzbudzającego. Sygnał wzbudzający jest wtedy podawany do wzbudnika 104 za pośrednictwem ścieżki 110. W kroku 302, nadajnik cyfrowy 20 odbiera sygnały z przetworników przesunięcia 105 i 105' spowodowane drganiami rur przepływowych 103A-B przy przechodzeniu materiału przez te rury przepływowe 103A-B. Dane o sygnałach, takie jak częstotliwość i różnica fazowa między sygnałami, są generowane przez nadajnik cyfrowy 20 w kroku 303. Następnie w kroku 304 obliczane są informacje o właściwościach materiału przepływającego przez rury przepływowe 103A i 103B. na przykład o strumieniu przepływu masowego, gęstości i strumieniu przepływu objętościoPL 197 438 B1 wego. Proces 300 następnie jest powtarzany przez cały czas pracy przepływomierza 5 Coriolisa w przewodzie rurowym.
Proces generowania danych o sygnałach przetworników przesunięcia według niniejszego wynalazku - fig. 4
Figura 4 ilustruje proces 400, który jest procesem generacji danych, takich jak częstotliwość sygnałów odbieranych z przetworników przesunięcia 105 i 105', które mierzą oscylacje rur przepływowych 103A-B w przepływomierzu 5 Coriolisa. Proces 400 stwarza kilka zalet dla pracy cyfrowego nadajnika 20. Pierwsza zaleta procesu polega na tym, że nie ma strat dokładności pomimo wykorzystywania arytmetyki stałoprzecinkowej zamiast arytmetyki zmiennoprzecinkowej. Umożliwia to stosowanie taniego procesora małej mocy. Druga zaleta polega na tym, że wymagania dotyczące pamięci dla procesu 400 są niewielkie na tyle, że może on rezydować w pamięci wewnętrznej procesora, co eliminuje konieczność dysponowania szybką magistralą między procesorem a pamięcią zewnętrzną. Wymagania dotyczące przetwarzania są zredukowane dzięki procesowi 400, który umożliwia pracę procesora z częstotliwością cyklu zegarowego sygnałów znacznie mniejszą od jego maksymalnej częstotliwości cyklu zegarowego.
Proces 400 rozpoczyna się w kroku 401 decymacja prędkości próbkowania sygnałów odbieranych z przetworników przesunięcia, z pierwszej prędkości próbkowania na drugą, mniejszą prędkość próbkowania. W korzystnej odmianie wykonania, sygnały są poddawane decymacji z pierwszej prędkości próbkowania wynoszącej 48 kHz do drugiej prędkości próbkowania wynoszącej 4 kHz. Decymacja prędkości próbkowania zwiększa rozdzielczość sygnałów, co podnosi dokładność obliczeń, jak to opisano poniżej na fig. 5. W korzystnej odmianie wykonania, zmniejszenie prędkości próbkowania z 48 kHz do 4 kHz zwiększa rozdzielczość próbki z B bitów do B + 1,79 bitów.
W kroku 402 z próbkowanych sygnałów oblicza się oszacowanie częstotliwości sygnału. Proces korzystny w przypadku obliczania szacowanej częstotliwości sygnału przedstawiono na fig. 6. Oszacowaną częstotliwość sygnału następnie wykorzystuje się do demodulowania odbieranych sygnałów w etapie 403. Proces demodulowania sygnałów cyfrowych podano na fig. 7 i 8. Druga decymacja próbkowanych sygnałów jest realizowana w kroku 404. Druga decymacja redukuje próbki sygnału z drugiej prędkości próbkowania do trzeciej prędkości próbkowania, z^i^k^^^j^<c rozdzielczość próbkowanego sygnału. W korzystnej odmianie wykonania, redukcja następuje z prędkości próbkowania wynoszącej 4 kHz do 800 Hz, co zwiększa rozdzielczość do B + 2,95 bitów. Ta redukcja jest wykonywana w taki sam sposób, jak decymacja w kroku 401.
Po drugiej decymacji w kroku 404, dokonuje się obliczeń na podstawie odebranych sygnałów. Wykonuje się to jako część procesu przedstawionego na fig. 9. Po usunięciu szumu, w kroku 405 wyznaczana jest częstotliwość sygnału pochodzącego z każdego przetwornika przesunięcia. W kroku 406 wyznacza się różnicę fazową między sygnałem z pierwszego przetwornika przesunięcia a sygnałem z drugiego przetwornika przesunięcia. Następnie w kroku 407 wyznacza się amplitudę każdego sygnału. Proces 400 następnie jest albo powtarzany przez cały czas pracy przepływomierza, albo proces 400 się kończy.
Proces decymacji prędkości próbkowania sygnałów z przetworników przesunięcia - fig. 5
Figura 5 przedstawia proces decymacji prędkości próbek odbieranych z przetworników przesunięcia. Ten sam proces wykorzystuje się do decymacji wykonywanej w każdym z kroków 401, 404 i w procesie wyznaczania częstotliwości sygnału. W każdym z tych kroków, proces 500 wykonuje się dla sygnałów z każdego przetwornika przesunięcia oddzielnie. Różnicą między decymacją wykonywaną w każdym kroku jest długość wejściowych wektorów danych, jak to opisano poniżej.
Decymacja opisana w procesie 500 jest zaimplementowana z wykorzystaniem blokowego sposobu przetwarzania, przy rozmiarze wektora wejściowego równym stosunkowi decymacji. Stosunek decymacji jest to wartość zmniejszenia częstotliwości próbkowania w wyniku decymacji. Wykorzystanie tego blokowego sposobu przetwarzania jest korzystne pod względem operacyjnym, ponieważ proces wymaga powtarzania dopiero z wyjściową, a nie wejściową prędkością danych. Zasada opóźnionej rekursywnej filtracji blokowej polega na tym, że reprezentacją zmiennej stanu sygnałów jest:
Xk+1 — A*Xk+Buk; yk — C*Xk+D*uk;
gdzie:
A,B,C,D — macierze reprezentujące stan systemu; Xk — N+1 wektor stanu w chwili czasowej k;
PL 197 438 B1 uk = dana wejściowa; a yk = dana wyjściowa reprezentująca sygnał po decymacji. Na podstawie indukcji jest oczywiste, że:
Xk+m ( Am Am-1B A m,B m-2 ... B
yk C D 0 ... 0
= CA CB D ... 0
yk+M~ CAm., CAm_,B CAm.2B ... D
Uk u
k+l
Przy decymacji pewnego sygnału z współczynnikiem M, zostaje zachowana tylko co M-ta próbka. Zatem, z powyższej macierzy można wyeliminować wszystkie wiesze z wyjątkiem ostatniego, otrzymując następujące równanie:
(x Ί -*k+M (A A-B -to
K^k+M-l) c CAM-2P Ca.U-?3
Na podstawie powyższego jest oczywiste, że liczba operacji sumowania/mnożenia w przypadku jednej rekursji powyższego równania wynosi:
Nmac = (N+1)*(N+M) gdzie:
Nmac = liczba operacji akumulowania/mnożenia;
N = rząd macierzy A; a
M = rozmiar bloku, który jest równy współczynnikowi decymacji procesu. Zatem obciążenie obliczeniami w przypadku wykonywania decymacji wynosi
Rsvd = Fout*N MAC gdzie:
Rsvd = obciążenie procesora obliczeniami, a
Fout = reprezentuje współczynnik wydajności wyjściowej filtru.
Do realizacji decymacji potrzebne są następujące pamięci:
pamięć do przechowywania każdego współczynnika filtru, która może być tylko pamięcią stałą (ROM);
pamięć do przechowywania wektora stanu Xk filtru, która musi być pamięcią do zapisu i odczytu (RAM); i pamięć wejściowego bloku buforowego współczynnika filtru (odczytowo-zapisowa).
Na fig. 5 przedstawiono proces decymacji z wykorzystaniem powyższego sposobu przetwarzania blokowego. Proces 500 rozpoczyna się w kroku 501 od odebrania m próbek do bufora i utworzenia bloku wejściowego. Blok wejściowy następnie w kroku 502 mnoży się przez wektor stanu. Wyniki wychodzące co m-tą próbkę następnie w kroku 503 są wyprowadzane do wykorzystania w innych obliczeniach. Po kroku 503 proces 500 kończy się.
Proces szacowania częstotliwości odbieranych sygnałów - fig. 6
Proces 600 jest procesem szacowania częstotliwości odbieranych sygnałów, w celu demodulowania sygnałów w następnym kroku procesu. Zanim będzie możliwość zdemodulowania sygnałów proces 600 musi się skończyć w kroku 403. Następną demodulację opisano poniżej i przedstawiono na fig. 7.
Na fig. 6 przedstawiono proces 600 szacowania częstotliwości sygnałów. Proces 600 wykonuje się na dowolnym z odebranych sygnałów. Proces 600 zaczyna się w kroku 601 od demultipleksowania próbkowanego sygnału cyfrowego do składowej współfazowej (l) i składowej kwadraturowej (Q). Następnie w kroku 602 odbywa się cyfrowe całkowanie składowej l i składowej Q sygnału. W kroku 603 na podstawie scałkowanego sygnału wylicza się sygnał kompensacji. Następnie w kroku 604 składowe sygnału są multipleksowane, z wygenerowaniem sygnału całkowanego cyfrowo. Następnie w kroku 605 oblicza się stosunek sygnału oryginalnego do sygnału scałkowanego cyfrowo. Stosunek
PL 197 438 B1 daje pewne oszacowanie częstotliwości sygnału, które można wykorzystać do demodulowania sygnałów w procesie 700.
W procesie 600 do oszacowania częstotliwości wykorzystuje się filtry o stałych współczynnikach. Ponieważ nie jest potrzebny żaden algorytm rekursywny, to proces 600 jest zawsze zbieżny. Poza tym, proces 600 nadąża za szybkimi zmianami częstotliwości. Szacowana częstotliwość na końcu procesu 600 jest dana następującym równaniem:
Fest = (oest/2n)x(Fx/12) gdzie:
Fest = częstotliwość szacowana; ω est = pulsacja znormalizowana; a Fk = częstotliwość próbek.
Proces dla dwustanowego selektora częstotliwości - fig. 8
Proces 800 przedstawiony na fig. 8 jest opcjonalnym procesem, który może być realizowany między oszacowaniem częstotliwości a demodulacją sygnału. Dwustanowy selektor częstotliwości jest potrzebny do wyznaczenia częstotliwości będącej przedmiotem zainteresowania. Proces 1000 (patrz fig. 9), który dokładnie mierzy częstotliwość sygnału wykazuje błąd systematyczny estymacji i wolniejszą zbieżność w zakresie częstotliwości normalizowanej:
|Fo| <0,5o|F0|>0,45 gdzie:
Fo = normalizowana częstotliwość sygnału.
Z tego równania wynika, że proces wyznaczania częstotliwości nie jest dokładny, kiedy prędkość próbkowania wynosi 4 kHz a częstotliwości mierzonego sygnału są rzędu tylko 20 Hz. Proces 800 zapobiega tej sytuacji umożliwiając wykorzystywania przetwarzania 1000 w dużym paśmie częstotliwości. Odbywa się to przez stwierdzenie, w następujący sposób, czy proces ma działać w trybie częstotliwości górnych, czy dolnych. Proces 800 rozpoczyna się w kroku 801 przez stwierdzenie, czy szacowana częstotliwość jest mniejsza, czy większa od częstotliwości odniesienia. W korzystnej odmianie wykonania wybiera się częstotliwość odniesienia równą 250 Hz. Częstotliwość 250 Hz wybrano jako częstotliwość między normalnymi częstotliwościami roboczymi konwencjonalnych rur przepływowych i prostych rur przepływowych.
Jeżeli rzeczywista oszacowana częstotliwość jest mniejsza od częstotliwości odniesienia, to w kroku 803 następuje zwrócenie jako częstotliwości szacowanej zera. Jeżeli rzeczywista częstotliwość jest większa od częstotliwości odniesienia, to częstotliwość szacowaną oblicza się w kroku 802 jako równą rzeczywistej oszacowanej częstotliwości zmniejszonej o 120 Hz.
Proces do demodulowania sygnałów odbieranych - fig. 7
Figura 7 przedstawia proces 700 do demodulowania odbieranych sygnałów cyfrowych. Proces wykorzystuje oszacowaną częstotliwość, która została albo wyliczona w procesie 600 (fig. 6) albo w procesie 800 (fig. 8). Proces 700 zaczyna się w kroku 701 obliczeniem znormalizowanej pulsacji, która wyraża się następującym równaniem:
ωά = 2^(12Fd)/Fs gdzie:
ad = pulsacja znormalizowana;
Fd = częstotliwość szacowana; a
Fs = częstotliwość próbek.
W kroku 702 oblicza się rzeczywisty współczynnik „skręcenia”, zgodnie z następującym równaniem: Wd = cos(o dk) i χβ = Αβ cos(o dk + φβ), gdzie:
β= sygnał odebrany od jednego z czujników pobierających.
Iloczyn skalarny „współczynnika skręcenia” i rzeczywistego odebranego sygnału oblicza się w kroku 702 z następującego równania:
yp = WkxP(K) = (Afc/2) {(^(((ω+ω^+φ) + cos((o-rod)k+φ)}
PL 197 438 B1
Należy zauważyć, że jeżeli proces 800 daje w wyniku tryb częstotliwościowy niski, to ustala się taki tryb, w którym szacowana częstotliwość jest równa zeru, y^k) = X(k)· W przeciwnym przypadku zmodulowany sygnał wyjściowy ma przedstawione poniżej dwie składowe:
F± = (Fs/12){(o ± Od)/2n}
Jednakowoż można temu zaradzić, przez opisaną poniżej decymację podwójną. Sygnałem będącym przedmiotem zainteresowania jest pierwsza składowa, której odpowiada znak w powyższym równaniu. Drugi sygnał, odpowiadający znakowi „+” w tym równaniu zostanie wyfiltrowany w następnym procesie decymacji, w kroku 905 procesu 900 (patrz fig. 9).
Proces generowania danych na podstawie odbieranych sygnałów - fig. 9
Figura 9 przedstawia proces 900 generowania danych o sygnałach odbieranych z czujników, proces 900 rozpoczyna się w kroku 901 od obliczenia pewnej adaptacji parametru filtru wycinającego, który oblicza się w sposób następujący. Wiadomo, że parametr filtru wycinającego jest parametrem adaptacyjnym reprezentowanym przez:
= 1 + a.,z 1 + z 2 (7 1 + om/' + cr/ 2
Gdzie α < 1 jest parametrem zbieżności korygującym szerokość pasma filtru, a a1 jest parametrem szukanym po całkowitej zaadaptowaniu. Załóżmy I aj <2, i wtedy zaznaczmy, że:
a1 = -2 cos(o)
Punkty zerowe charakterystyki H(z) dane są równaniem:
z± = exp(± j o) gdzie z = punkty zerowe; j = a stała; a ω = sygnał.
Zatem bieguny sygnałów wyrażone są w następującym równaniu:
ρ± = α exp(±jo)
W kroku 902, „a” jest obliczane dla każdego sygnału. Wartość „a” jest obliczana z wykorzystaniem jednego z wielu algorytmów konwencjonalnych, jak RLS, RML lub SGN. Minimalizuje to energię sygnału.
W kroku 902 wyznacza się częstotliwość każdego sygnału. W celu wyznaczenia częstotliwości sygnału, dokonuje się wyznaczenia częstotliwości normalizowanej w odniesieniu do sygnału po decymacji, z wykorzystaniem następującego równania:
Fo = (1/2n)arccos(-a-|/2);
gdzie:
Fo = częstotliwość normalizowana w odniesieniu do sygnału decymowanego; a a1 jest aktualną wartością adaptacyjną parametru filtru wycinającego.
Można następnie, w kroku 9o2, wyznaczyć częstotliwość sygnałów przez pomnożenie częstotliwości normalizowanej przez częstotliwość decymacji (Fd = Fo x Fs). W kroku 9o3 na sygnałach dokonuje się modulacji kwadraturowej. Modulacja kwadraturowa jest wykonywana przez takie wybranie sygnału demodulacyjnego, jak pokazane poniżej, gdzie sygnały dominującej częstotliwości są przesunięte do zera. W procesie 9o3 oblicza się iloczyn skalarny sygnału demodulacyjnego i sygnałów odbieranych. Sygnał modulacyjny jest reprezentowany w sposób następujący:
o o = 2n Fo gdzie:
o o = pulsacja sygnału modulującego; a
Fo = częstotliwość normalizowanego sygnału obliczonego powyżej w kroku 9o2.
Jak wspomniano powyżej, odbierane sygnały można przedstawić jako:
x.<k) = A cos (ook + φβ;
gdzie β = każdy sygnał z czujnika 1o5 -1o5' przesunięcia.
PL 197 438 B1
Z powyższych równań wynika, że wynikiem demodulacji kwadraturowej jest:
zp(k) = Wixp(k) = A/2 (θχρ(ίφβ) + exp(-j(2o ok + φβ))}.
W celu dalszego zwiększenia rozdzielczości sygnału, dokonuje się decymacji w kroku 904. W korzystnej odmianie wykonania ta decymacja jest decymacją 40-krotną, która jest wykonywana na składowych zarówno l, jak i Q, odbieranych sygnałów. Decymacja ta redukuje wynik zespolonej demodulacji kwadraturowej do:
zp(k) = (a/2) exp(ju).
Po wykonaniu decymacji, w kroku 905 otrzymuje się różnicę fazową. W przykładowej odmianie wykonania różnicę fazową oblicza się w sposób następujący. Pierwszy z sygnałów, odbierany albo z lewego, albo z prawego czujnika przesunięcia, jest sprzężony, zgodnie z następującym równaniem:
z rpo(k) (Arp0/2)exp(-j®rp0)
Następnie ten sygnał jest mnożony z drugim sygnałem z otrzymaniem zespolonej korelacji między sygnałami przetworników przesunięcia, zgodnie z następującym równaniem:
- A2 q(k) = z1s,(k)z2nd(k) = — expjMw2nd))
Zatem następnie, w kroku 906 wyznaczana jest różnica fazowa, i jest dana następującym równaniem:
Φ(Μ = arg(q(k)) = Φ15 - φ2nd
Różnicę fazową następnie można wykorzystać do obliczania strumienia przepływu masowego i innych właściwości materiału przepływającego przez materiał.
Powyżej zamieszczono opis nadajnika cyfrowego 20 do przepływomierza 5 Coriolisa oraz procesy wyznaczania danych dotyczących sygnałów odbieranych przez nadajnik 20.
Należy oczekiwać, że możliwe jest zaprojektowanie alternatywnych cyfrowych procesorów sygnałowych i procesów, stanowiących naruszenie niniejszego wynalazku, przedstawionego w załączonych zastrzeżeniach, albo literalnie albo na zasadzie rozwiązań równoważnych.

Claims (11)

  1. Zastrzeżenia patentowe
    1. Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa, polegający na tym, że wprawia się przewód w drgania i mierzy się drgania przewodu za pomocą pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia jednocześnie wytwarzając sygnał cyfrowy, a następnie dostarcza się próbki sygnałów otrzymywanych z pierwszego czujnika przesunięcia i drugiego czujnika przesunięcia z pierwszą prędkością próbkowania do cyfrowego procesora sygnałowego, po czym zmienia się prędkość próbkowania z pierwszej prędkości próbkowania do pośredniej prędkości próbkowania za pomocą pierwszej decymacji tych próbek, znamienny tym, że dokonuje się demodulacji (403) próbek o pośredniej prędkości próbkowania i zmienia się prędkość próbkowania demodulowanych próbek do wymaganej prędkości próbkowania poprzez przeprowadzenie drugiej decymacji (404) próbek sygnałów, a następnie z tych próbek decymowanych sygnałów wyznacza się częstotliwość (405) drgań rury przepływowej (103A, 103B) przy pierwszym czujniku przesunięcia (105) i przy drugim czujniku przesunięcia (105'), zaś na podstawie tych częstotliwości wyznacza się właściwości materiału przepływającego przez rurę przepływową (103A, 103B).
  2. 2. Sposób według zastrz. 1, znamienny tym, że jako częstotliwości sygnałów wyznacza się znormalizowane częstotliwości sygnałów.
  3. 3. Sposób według zastrz. 2, znamienny tym, że demodulowania próbek dokonuje się obliczając (701) znormalizowaną pulsację znormalizowanej częstotliwości sygnałów i obliczając (703) iloczyn skalarny pulsacji znormalizowanej i sygnałów pierwszego czujnika przesunięcia (105) i drugiego czujnika przesunięcia (105'), z przeniesieniem sygnałów na częstotliwość środkową.
  4. 4. Sposób według zastrz. 2, znamienny tym, że znormalizowaną częstotliwość sygnałów wyznacza się poprzez przeprowadzenie demultipleksowania (601) sygnałów na składowe współfazowe i składowe kwadraturowe, a następnie całkuje się (602) składowe współfazowe i składowe kwadraturo12
    PL 197 438 B1 we i kolejno dokonuje się multipleksowania (604) składowych współfazowych i składowych kwadraturowych wytwarzając sygnały całkowane cyfrowo, po czym wyznacza się (605) stosunek amplitudy sygnałów do amplitudy sygnałów całkowanych cyfrowo otrzymując znormalizowane częstotliwości sygnałów.
  5. 5. Sposób według zastrz. 4, znamienny tym, że po dokonaniu całkowania (602) i przed dokonaniem multipleksowania (604) całkowane składowe kwadraturowe przesyła się (603) do kompensatora.
  6. 6. Sposób według zas^z. 4, znam ienny tym i że po dokonaniu caakowania (602) i preed dokonaniem multipleksowania (604) przesyła się (603) całkowane składowe współfazowych do kompensatora.
  7. 7. Sposób weeług ο^ιζ. (, znamienny tym i że cczsto-llwośś drggńiury prgzeruwowe- ((00A, 103B) wyznacza się poprzez odczytanie (900) znormalizowanej częstotliwości sygnałów, a następnie dokonanie modulacji (902) znormalizowanej częstotliwości sygnałów i przeprowadzenie zespolonej demodulacji (903) sygnałów z wykorzystaniem zmodulowanej częstotliwości znormalizowanej.
  8. 8. Sposób według zassic 7, znamienny tym, że wyznaczając często-llwości drgań rury przepływowej (103A, 103B) dokonuje się decymacji (904) zdemodulowanych sygnałów i dokonuje się (905) zespolonej korelacji sygnałów z obu czujników przesunięcia (105, 105') wyznaczając różnicę fazową między sygnałami.
  9. 9. Sposób według zas^z. 1, znamienny tym. że wyznaczając często-llwości drgań rury przepływowej (103A, 103B) mierzy się (910) różnicę fazy między sygnałami z pierwszego czujnika przesunięcia (105), a sygnałami z drugiego czujnika przesunięcia (105').
  10. 10. Sposób według 1, znamienny tym, że jako wł^^^^ii^c^^tó materiału pr^^^f^pł^i^^j^c^^ę^o przez rurę przepływową (103A, 103B) wyznacza się jego natężenie przepływu masowego materiału.
  11. 11. Sposób według zas^z. 1, znamienny tym, że jako właściwość matenału przepPywającego przez rurę przepływową (103A, 103B) wyznacza się jego gęstość.
PL352214A 1999-06-28 2000-06-12 Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa PL197438B1 (pl)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/344,840 US6505131B1 (en) 1999-06-28 1999-06-28 Multi-rate digital signal processor for signals from pick-offs on a vibrating conduit
PCT/US2000/016061 WO2001001083A1 (en) 1999-06-28 2000-06-12 Multi-rate digital signal processor for vibrating conduit sensor signals

Publications (2)

Publication Number Publication Date
PL352214A1 PL352214A1 (en) 2003-08-11
PL197438B1 true PL197438B1 (pl) 2008-03-31

Family

ID=23352280

Family Applications (1)

Application Number Title Priority Date Filing Date
PL352214A PL197438B1 (pl) 1999-06-28 2000-06-12 Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa

Country Status (18)

Country Link
US (1) US6505131B1 (pl)
EP (1) EP1190222B1 (pl)
JP (1) JP3713236B2 (pl)
KR (1) KR100472994B1 (pl)
CN (1) CN1199034C (pl)
AR (1) AR024467A1 (pl)
AT (1) ATE316237T1 (pl)
AU (1) AU766488B2 (pl)
BR (1) BR0012075B8 (pl)
CA (1) CA2378027C (pl)
DE (1) DE60025607T2 (pl)
DK (1) DK1190222T3 (pl)
HK (1) HK1047789B (pl)
MX (1) MXPA02000126A (pl)
MY (1) MY125425A (pl)
PL (1) PL197438B1 (pl)
RU (1) RU2234683C2 (pl)
WO (1) WO2001001083A1 (pl)

Families Citing this family (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8447534B2 (en) 1997-11-26 2013-05-21 Invensys Systems, Inc. Digital flowmeter
US6311136B1 (en) * 1997-11-26 2001-10-30 Invensys Systems, Inc. Digital flowmeter
US7404336B2 (en) 2000-03-23 2008-07-29 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
US7784360B2 (en) 1999-11-22 2010-08-31 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
US20030216874A1 (en) 2002-03-29 2003-11-20 Henry Manus P. Drive techniques for a digital flowmeter
US8467986B2 (en) * 1997-11-26 2013-06-18 Invensys Systems, Inc. Drive techniques for a digital flowmeter
US7124646B2 (en) * 1997-11-26 2006-10-24 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
DE01918944T1 (de) * 2000-03-23 2004-10-21 Invensys Systems, Inc., Foxboro Korrektur für eine zweiphasenströmung in einem digitalen durchflussmesser
DE10045537A1 (de) * 2000-09-13 2002-03-21 Flowtec Ag Mess- und Betriebsschaltung für einen Coriolis-Massedurchflußmesser
US6694279B2 (en) * 2001-02-16 2004-02-17 Micro Motion, Inc. Methods, apparatus, and computer program products for determining structural motion using mode selective filtering
GB2391313B (en) * 2002-07-26 2006-01-18 Abb Ltd Deriving measurements from periodic data
US7059199B2 (en) * 2003-02-10 2006-06-13 Invensys Systems, Inc. Multiphase Coriolis flowmeter
US7188534B2 (en) * 2003-02-10 2007-03-13 Invensys Systems, Inc. Multi-phase coriolis flowmeter
US7013740B2 (en) * 2003-05-05 2006-03-21 Invensys Systems, Inc. Two-phase steam measurement system
DE10322851A1 (de) * 2003-05-19 2004-12-16 Endress + Hauser Flowtec Ag, Reinach Coriolis-Durchflußmeßgerät
GB2402222B (en) * 2003-05-30 2007-04-25 Abb Ltd Phase measurement in measuring device
US7072775B2 (en) * 2003-06-26 2006-07-04 Invensys Systems, Inc. Viscosity-corrected flowmeter
US7065455B2 (en) * 2003-08-13 2006-06-20 Invensys Systems, Inc. Correcting frequency in flowtube measurements
KR101206350B1 (ko) * 2004-12-29 2012-12-04 마이크로 모우션, 인코포레이티드 유량계를 위한 고속 주파수 및 위상 평가
EP1889012B1 (en) * 2005-05-20 2019-10-30 Micro Motion, Inc. Meter electronics and methods for determining void fraction of gas
EP1889010B1 (en) * 2005-05-20 2016-12-28 Micro Motion, Inc. Meter electronics and methods for rapidly determinig a mass fraction of a multi-phase fluid from a coriolis flow meter signal
KR20130022429A (ko) * 2005-05-27 2013-03-06 마이크로 모우션, 인코포레이티드 코리올리 유량계를 통해 유동하는 물질의 비균일성을 신속하게 탐지하기 위한 방법 및 계측 전자장치
US8781759B2 (en) * 2005-08-18 2014-07-15 Micro Motion, Inc. Meter electronics and methods for processing sensor signals for a multi-phase flow material in a flowmeter
BRPI0520555B1 (pt) * 2005-09-19 2017-11-14 Micro Motion, Inc. Electronic meter apparatus and method for determining a rigidity parameter of a flow meter
RU2376556C1 (ru) * 2005-09-20 2009-12-20 Майкро Моушн, Инк. Измерительная электроника и способы для генерирования сигнала возбуждения для вибрационного расходомера
DK1949047T3 (en) * 2005-10-18 2015-04-20 Micro Motion Inc Electronic device for measurement and the methods for determining a difference in phase between a first sensor signal and a second sensor signal of a flow meter
US7617055B2 (en) 2006-08-28 2009-11-10 Invensys Systems, Inc. Wet gas measurement
US8290630B2 (en) * 2008-09-30 2012-10-16 Rockwell Automation Technologies, Inc. Condition monitoring parameter normalization system and method
JP4436884B1 (ja) * 2009-02-06 2010-03-24 株式会社オーバル 信号処理方法、信号処理装置、およびコリオリ流量計
EP2430406B1 (en) 2009-05-11 2021-01-13 Micro Motion, Inc. A flow meter including a balanced reference member and its method
US20120041695A1 (en) 2010-08-16 2012-02-16 Csi Technology, Inc. Integrated vibration measurement and analysis system
US9625295B2 (en) * 2011-08-08 2017-04-18 Hach Company Doppler measurement instrument and doppler measurement processing method
BR112015026826B1 (pt) 2013-04-23 2020-11-10 Micro Motion, Inc método de gerar um sinal de acionamento para um sensor vibratório, e, sensor vibratório
WO2017143577A1 (en) * 2016-02-26 2017-08-31 Micro Motion, Inc. Meter electronics for two or more meter assemblies
EP3420319B1 (en) * 2016-02-26 2023-05-10 Micro Motion, Inc. Communicating with two or more hosts
CN114112020A (zh) * 2016-02-26 2022-03-01 高准公司 确定振动传感器类型
CN107131920A (zh) * 2016-02-26 2017-09-05 高准公司 用于计量电子器件的子板
EP3433586B1 (en) * 2016-03-25 2021-06-23 Micro Motion, Inc. Method for maximizing flowmeter turndown and related apparatus
AU2020220614B2 (en) 2019-02-12 2022-07-28 Paltech Process for moulding polymeric foam core sandwich articles
EP4095492B1 (de) 2021-05-26 2023-09-27 Krohne Messtechnik GmbH Verfahren zum betreiben eines coriolis-massedurchflussmessgeräts und entsprechendes coriolis-massedurchflussmessgerät

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1219887A (en) 1969-02-10 1971-01-20 Robert Baron Jacobs Fluid flow measurement systems
FR2321217A1 (fr) * 1975-08-13 1977-03-11 Cit Alcatel Dispositif de traitement d'un signal echantillonne
US4996871A (en) 1989-06-02 1991-03-05 Micro Motion, Inc. Coriolis densimeter having substantially increased noise immunity
HU216462B (hu) * 1991-08-01 1999-06-28 Micro Motion Inc. Eljárás és áramlásmérő tömegáramlásisebesség-információ meghatározására
DE4316297C1 (de) * 1993-05-14 1994-04-07 Fraunhofer Ges Forschung Frequenzanalyseverfahren
US5361036A (en) * 1993-08-12 1994-11-01 Rockwell International Corporation Complex digital demodulator employing Chebychev-approximation derived synthetic sinusoid generation
US5741980A (en) * 1994-11-02 1998-04-21 Foster-Miller, Inc. Flow analysis system and method
US5555190A (en) * 1995-07-12 1996-09-10 Micro Motion, Inc. Method and apparatus for adaptive line enhancement in Coriolis mass flow meter measurement
US5734190A (en) * 1996-03-11 1998-03-31 Eastman Kodak Company Imager having a plurality of cylindrical lenses
US5926096A (en) * 1996-03-11 1999-07-20 The Foxboro Company Method and apparatus for correcting for performance degrading factors in a coriolis-type mass flowmeter
US5734112A (en) 1996-08-14 1998-03-31 Micro Motion, Inc. Method and apparatus for measuring pressure in a coriolis mass flowmeter
EP0866319B1 (de) 1997-03-21 1999-11-24 Endress + Hauser Flowtec AG Mess- und Betriebsschaltung eines Coriolis-Massedurchflussmessers
JPH11118848A (ja) * 1997-10-14 1999-04-30 Advantest Corp スペクトラムアナライザ

Also Published As

Publication number Publication date
AR024467A1 (es) 2002-10-02
US6505131B1 (en) 2003-01-07
ATE316237T1 (de) 2006-02-15
WO2001001083A1 (en) 2001-01-04
KR20020019932A (ko) 2002-03-13
JP3713236B2 (ja) 2005-11-09
BR0012075B8 (pt) 2013-02-19
CN1199034C (zh) 2005-04-27
CA2378027A1 (en) 2001-01-04
AU766488B2 (en) 2003-10-16
EP1190222A1 (en) 2002-03-27
PL352214A1 (en) 2003-08-11
HK1047789A1 (en) 2003-03-07
HK1047789B (zh) 2005-12-09
DE60025607T2 (de) 2006-08-03
EP1190222B1 (en) 2006-01-18
DK1190222T3 (da) 2006-05-22
MXPA02000126A (es) 2002-07-30
KR100472994B1 (ko) 2005-03-09
JP2003503690A (ja) 2003-01-28
RU2234683C2 (ru) 2004-08-20
MY125425A (en) 2006-07-31
DE60025607D1 (de) 2006-04-06
AU5606100A (en) 2001-01-31
BR0012075A (pt) 2002-03-26
CA2378027C (en) 2006-12-19
BR0012075B1 (pt) 2013-01-08
CN1371471A (zh) 2002-09-25

Similar Documents

Publication Publication Date Title
PL197438B1 (pl) Sposób określania właściwości materiału przepływającego przez rurę przepływową przepływomierza typu Coriolisa
JP4960967B2 (ja) 流量計の第1のセンサ信号と第2のセンサ信号との間の位相差を決定するための計器電子装置及び方法
KR100342180B1 (ko) 코리올리 질량유량계의 압력측정방법과 장치
AU2005322423B2 (en) High speed frequency and phase estimation for flow meters
JP4373781B2 (ja) コリオリ流量計における流管と流体との特性の決定
JP2005502041A5 (pl)
RU2328703C2 (ru) Программируемые электронные средства расходомера кориолиса, способ эксплуатации электронных средств и считываемый процессором носитель программного продукта
WO2000004347A1 (en) Improved vibrating conduit process parameter sensors, operating methods and computer program products utilizing complex modal estimation