JP2000018995A - ディジタル流量計 - Google Patents

ディジタル流量計

Info

Publication number
JP2000018995A
JP2000018995A JP10378592A JP37859298A JP2000018995A JP 2000018995 A JP2000018995 A JP 2000018995A JP 10378592 A JP10378592 A JP 10378592A JP 37859298 A JP37859298 A JP 37859298A JP 2000018995 A JP2000018995 A JP 2000018995A
Authority
JP
Japan
Prior art keywords
flow meter
signal
sensor
conduit
digital
Prior art date
Legal status (The legal status 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 status listed.)
Withdrawn
Application number
JP10378592A
Other languages
English (en)
Inventor
P Manus Henry
ピー. ヘンリー マニュス
W David Clarke
ダブリュ. クラーク ディビッド
H James Vignos
エイチ. ヴィグノス ジェイムズ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Schneider Electric Systems USA Inc
Original Assignee
Foxboro Co
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=26746870&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=JP2000018995(A) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by Foxboro Co filed Critical Foxboro Co
Publication of JP2000018995A publication Critical patent/JP2000018995A/ja
Withdrawn legal-status Critical Current

Links

Classifications

    • 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/8404Coriolis or gyroscopic mass flowmeters details of flowmeter manufacturing methods
    • 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/8431Coriolis or gyroscopic mass flowmeters constructional details electronic circuits
    • 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/8481Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits having loop-shaped measuring conduits, e.g. the measuring conduits form a loop with a crossing point
    • G01F1/8486Coriolis or gyroscopic mass flowmeters arrangements of measuring means, e.g., of measuring conduits vibrating measuring conduits having loop-shaped measuring conduits, e.g. the measuring conduits form a loop with a crossing point with multiple measuring conduits
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F25/00Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
    • G01F25/10Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
    • G01N9/002Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity using variation of the resonant frequency of an element vibrating in contact with the material submitted to analysis
    • G01N2009/006Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity using variation of the resonant frequency of an element vibrating in contact with the material submitted to analysis vibrating tube, tuning fork

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Signal Processing (AREA)
  • Measuring Volume Flow (AREA)

Abstract

(57)【要約】 (修正有) 【課題】従来のアナログの方法に対して、より大なる反
応性、精度、順応性を有し、正確且つ精巧な制御アルゴ
リズムの適用が可能なディジタル質量流量計を提供す
る。 【解決手段】 ディジタル流量計100は、振動可能な
導管と、導管に接続されて導管に運動を与えるように操
作可能な駆動回路115と、導管に接続されて導管の運
動を検出するように操作可能なセンサ110とを含む。
デジタルコントローラ105が、駆動回路とセンサとの
間に接続される。デジタルコントローラ105は、回路
を含み、この回路は、センサからセンサ信号を受け取
り、デジタル信号処理を使用してセンサ信号に基づいた
駆動信号を生成し、駆動信号を駆動回路へ出力し、セン
サからの信号に基づいて導管を流れる材料の特性の測定
値を生成する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、ディジタル流量計
に関する。
【0002】
【背景技術】流量計は、導管を移動する材料に関する情
報を提供する。例えば、質量流量計は、導管を移動する
材料の直接表示を行う。同様に、密度流量計、すなわち
比重計は、導管の中を流れる材料の密度の表示を行う。
質量流量計も、材料の密度の表示を行うこともある。
【0003】コリオリ型の質量流量計は、周知のコリオ
リ効果に基づいており、回転導管を流れる材料は、コリ
オリの力による作用を受けて加速されて放射状に移動す
る塊になる。多くのコリオリ型の質量流量計は、導管の
長手方向と直交する回動軸を中心に導管をサイン波状に
振動させることによって、コリオリ力を誘導する。かか
る質量流量計において、移動流体質量によって受けるコ
リオリの反応力は、導管そのものに移動して、回転面内
のコリオリ力のベクトルの方向に、導管の偏差、すなわ
ち偏位として発現する。
【0004】エネルギは、周期的な力を印加する駆動機
構によって導管に供給されて、導管を振動させる。駆動
機構の1つのタイプは、印加電圧と比例した力を与える
電気機械的式駆動回路である。振動流量計において、印
加電圧は周期的であり、一般にサイン波状である。入力
電圧の周期は、導管の運動が導管の振動の共振モードと
整合するように選択される。これによって、振動を維持
するために必要なエネルギが減少する。振動流量計は、
導管の振動に対応する瞬間的な周波数および位相情報を
担持するセンサ信号が増幅されて、電気機械的式駆動回
路を使用する導管にフィードバックされる、フィードバ
ックループを使用することがある。
【0005】
【発明の概要】本発明は、制御・測定システムを使用し
て導管の振動を制御するとともに質量流および密度の測
定値を生成する、ディジタル質量流量計などのディジタ
ル流量計を提供する。導管に接続されたセンサは、制御
・測定システムに信号を供給する。制御・測定システム
は、信号を処理して質量流の測定値を生成し、ディジタ
ル信号処理を使用して、導管を駆動する信号を生成す
る。次に、駆動信号は、導管の振動を誘起する力に変換
される。
【0006】ディジタル質量流量計は、従来の、アナロ
グの方法に勝る多くの効果を提供する。制御の点から、
ディジタル処理技術を使用すると、従来のアナログの方
法に対して、より大なる反応性、精度、順応性を提供す
る、正確且つ精巧な制御アルゴリズムの適用が可能にな
る。ディジタル制御システムによって、導管の振動を制
御する際、負の利得を使用することができる。このよう
に、導管振動について位相が180度ずれている駆動信
号が、振動の振幅を減らすために使用することができ
る。これの実際の関連性は、特に減衰での突然の低下が
振動振幅の望ましくない増加を生成する、大きく且つ可
変の減衰状況においては、特に重要である。可変減衰状
況の1つの例は、空気混和が導管の中を流れる材料に生
じる時である。
【0007】振動の大きさが、ユーザ制御により変更で
きる一定の設定値に制御されるとき、負のフィードバッ
クを行う能力も重要である。負のフィードバックによっ
て、振動の設定値の減少は、設定値の増加と同様に高速
で実行できる。対照的に、正のフィードバックのみに頼
るアナログ計器は、利得をゼロに設定して、振幅を減ら
した設定値に下げるために、システムの減衰を待たなけ
ればならない。
【0008】測定の点から、ディジタル質量流量計は、
高い情報帯域幅を提供できる。例えば、ディジタル測定
システムは、18ビットの精度で且つ55kHzのサン
プリングレートで作動しているアナログーディジタルコ
ンバータを使用する。ディジタル測定システムも、デー
タをフィルタにかけて処理する精巧なアルゴリズムを使
用し、センサからの生のデータで開始して最終的な測定
データまで継続する。これによって、例えば1周期につ
き5ナノ秒の位相精度等の、極めて高い精度が実現され
る。生のセンサデータで始まるディジタル処理も、既存
の測定技術への拡張も斟酌して、時間変化する振幅、周
波数、ゼロオフセットを検出して補償すること等によっ
て、現実的な状況での性能を改善する。
【0009】制御および測定の改良は、相互に作用して
更に改善される。例えば、振動振幅の制御は、振幅測定
値の質に依存する。通常の条件下において、ディジタル
質量流量計は、所望の設定値の20ppm内に振動を維
持する。同様に、制御が改善されると、測定値に正の利
益をもたらす。振動の安定性が増加すると、一定の振幅
(すなわち、一定の設定値)を必要としない計器に対し
てさえも、測定値の質を改善する。例えば、改善された
安定性によって、測定値の計算に使用される仮定は、広
範囲の条件に対して有効になる。
【0010】ディジタル質量流量計によって、測定およ
び制御プロセスを有する完全に新しい機能(例えば診断
法)の統合が可能になる。例えば、空気混和プロセスの
有無を検出するアルゴリズムは、空気混和が検出される
場合、測定及び制御の両方で生じる補償動作によって実
行される。ディジタル質量流量計のさらなる効果は、使
用されるハードウェアの限られた量から生じる。そし
て、これによって、製造時や現場での組立、デバッグ、
修理を行うために、計器が簡単になる。改善された性能
のための現場での迅速な修理や、機械部品(例えばルー
プ、フランジ、センサ、駆動回路)の摩耗の補償は、可
能である。何となれば、計器は、容易に交換可能な規格
のハードウェア部品を使用しており、ソフトウェアの変
更を比較的容易に行うことができるからである。更に、
診断、測定、制御の統合は、ハードウェアの簡単さとソ
フトウェアにおいて実行される機能のレベルとによっ
て、容易になる。例えば低電力部品や、性能が改善され
た部品などの、新しい機能は、制御システム全体の主要
な再設計をせずに、統合される。
【0011】1の概念において、ディジタル流量計は、
振動可能な導管と、導管に接続されて導管に運動を与え
るように動作可能な駆動回路と、導管に接続されて導管
の運動を検出するように動作可能なセンサとを含む。駆
動回路とセンサの間に接続された制御・測定システム
は、センサからセンサ信号を受け取り、ディジタル信号
処理を使用してセンサ信号に基づいた駆動信号を生成
し、駆動信号を駆動回路に供給し、センサからの信号に
基づいて導管を流れる材料の特性の測定値を生成する電
気回路を含む。
【0012】実施例は、以下の特徴を含む。計器は、導
管に接続されて導管の運動を検出するように動作可能な
第2のセンサを含む。この場合、制御・測定システム
は、第2のセンサに接続されて、第2のセンサから第2
のセンサ信号を受け取り、第1および第2のセンサ信号
に基づいて駆動信号を生成し、第1および第2のセンサ
信号に基づいて導管の中を流れる材料の特性の測定値を
生成する。制御・測定システムは、第1および第2のセ
ンサ信号をディジタル的に合成して、センサ信号の組合
せに基づいて駆動信号を生成する。
【0013】制御・測定システムは、2つの駆動回路に
対して別々の駆動信号を生成してもよい。駆動信号は、
例えば異なる周波数や振幅を有することがある。ディジ
タル流量計は、駆動回路に供給される電流を測定する回
路を含むことがある。回路は、駆動回路と直列の抵抗器
と、抵抗器と並列であるとともに抵抗器の電圧を測定す
るように構成されたアナログ−ディジタルコンバータと
を含み、測定電圧をディジタル値に変換して、ディジタ
ル値を制御・測定システムに供給する。
【0014】ディジタル流量計には、導管の入口で第1
の圧力を測定するために接続された第1の圧力センサ
と、導管の出口で第2の圧力を測定するために接続され
た第2の圧力センサとを含む。アナログ−ディジタルコ
ンバータが、接続されて、第1の圧力センサ及び第2の
圧力センサにより生成された信号をディジタル値に変換
し、ディジタル値を制御・測定システムに供給するよう
に構成されている。温度センサを、導管の入口及び出口
での温度を測定するために接続することができる。
【0015】制御・測定システムは、第1のセンサ信号
の周波数を評価し、第1のセンサ信号を使用して位相差
を計算し、更に計算された位相差を使用して測定値を生
成することによって、特性の測定値を生成する。制御・
測定システムは、一方のセンサ信号の振幅を調節するこ
とによって、センサ信号の振幅差を補償する。例えば、
制御・測定システムは、一方のセンサ信号の振幅にセン
サ信号の振幅の比を乗算する。
【0016】センサ信号が周期的である場合、制御・測
定システムは、1セット内でセンサ信号を処理する。各
セットは、周期的センサ信号の全周期に対するデータを
含み、連続したセットは、周期的センサ信号の周期を重
畳するデータを含む。制御・測定システムは、前の周期
の周波数を使用して、周期の終点を評価する。制御・測
定システムは、周期が更なる処理に値するかどうかを判
別するために1周期のデータを分析する。例えば、シス
テムは、周期のデータがデータとして予測された動作に
一致しないとき、周期はさらなる処理に値しないと判別
する。ここで、予測された動作は、前の周期のパラメー
タに基づいている。一実施例において、システムは、周
期に対する周波数が前の周期に対する周波数とは閾値以
上異なるとき、周期はさらなる処理に値しないと判別す
る。システムは、周期に対する周波数が前の周期に対す
る周波数と等しければ、周期内のある点での値を、生じ
るであろう値に比較することによって、周波数が異なる
か否かを判別する。
【0017】制御・測定システムは、センサ信号のゼロ
クロスを検出し、ゼロクロスの間のサンプルを数えるこ
とによって、センサ信号の周波数を測定する。システム
は、反復曲線あてはめ技術を使用して、センサ信号の周
波数を測定する。制御・測定システムは、フーリエ解析
を使用して、センサ信号の振幅を測定し、駆動信号を生
成する際の測定された振幅を使用する。
【0018】制御・測定システムは、各センサ信号に対
する位相偏差を測定し、位相偏差を比較することによっ
て位相差を測定する。システムは、フーリエ解析を使用
して、位相差を測定することもできる。制御・測定シス
テムは、各センサ信号に対する、周波数、振幅、位相偏
差を測定し、センサ信号の周波数の平均に対して位相偏
差を拡大・縮小する。制御・測定システムは、複数の方
法を使用して、位相差を計算し、1つの方法の結果を計
算された位相差として選択する。
【0019】制御・測定システムは、センサ信号を合成
して合成信号を生成し、合成された信号に基づいて駆動
信号を生成する。例えば、制御・測定システムは、セン
サ信号を合計して、合成信号を生成し、利得を合成信号
に適用することによって駆動信号を生成する。一般に、
制御・測定システムは、第1のモードの信号生成を使用
して駆動信号を生成することによって、導管の運動を開
始し、第2のモードの信号生成を使用して駆動信号を生
成することによって、導管の運動を維持する。第4のモ
ードの信号生成は、導管振動の所望の初期周波数等の、
所望の特性を有する周期的信号の合成であり、第2のモ
ードの信号生成は、センサ信号を含むフィードバックル
ープを使用する。
【0020】他の実施例において、第1のモードの信号
生成は、センサ信号を含むフィードバックループの使用
を含み、第2のモードの信号生成は、所望の特性を有す
る周期信号の合成を含む。例えば、制御・測定システム
は、大なる利得を合成信号にに適用して導管の運動を開
始させ、運動の開始後の駆動信号として、センサ信号の
位相及び周波数に基づいた位相および周波数を有する周
期的信号を生成することによって、駆動信号を生成す
る。合成された信号の所望の特性は、センサ信号の周波
数および位相に相当する周波数および位相である。
【0021】制御・測定システムは、順応性のある周期
的駆動信号を生成する。例えば、計器は、制御・測定シ
ステムと駆動回路との間に接続された正及び負の直流電
源を含み、制御・測定システムは、センサ信号に基づい
た位相および周波数を有する間隔で電流源のオン・オフ
を切り替えることによって、駆動信号を生成する。制御
・測定システムは、センサ信号の位相及び周波数に相当
する位相及び周波数等の、センサ信号の特性に相当する
特性を有するサイン波を合成することによって、駆動信
号を生成する。
【0022】制御・測定システムは、センサ信号の特性
に基づいて駆動信号を生成する時に使用する利得をディ
ジタル的に生成する。例えば、制御・測定システムは、
センサ信号の振幅に基づいて、利得をディジタル的に生
成する。駆動回路は、導管を振動運動させるように動作
可能である。制御・測定システムは、導管振動の振幅を
調整するPI制御アルゴリズムをディジタル的に実行す
る。制御・測定システムは、導管の振動振幅をユーザが
制御する値に維持するために、センサ信号に基づいて駆
動信号をディジタル的に生成する。これに賛成して、制
御・測定システムは、振幅がユーザ制御の値を越えると
きは、負の駆動信号を生成し、これによって、駆動回路
は、導管の運動に対抗し、振動振幅がユーザ制御の値未
満のときは、正の駆動信号を生成し、これによって、駆
動回路は、導管を運動させる。
【0023】制御・測定システムは、センサ信号に基づ
いて利得信号を生成するコントローラと、コントローラ
に接続されて利得信号を受け取るとともにこの利得信号
に基づいて駆動信号を生成する乗算ディジタル−アナロ
グコンバータとを含む。ディジタル流量計が、導管に接
続されて導管の運動を検出するように動作可能な第2の
センサを含むとき、制御・測定システムは、測定値を生
成するコントローラと、第1のセンサとコントローラと
の間に接続されて第1のディジタルセンサ信号をコント
ローラに出力する第1のアナログ−ディジタルコンバー
タと、第2のセンサとコントローラとの間に接続されて
第2のディジタルセンサ信号をコントローラに出力する
第2のアナログ−ディジタルコンバータとを含む。コン
トローラは、ディジタルセンサ信号を合成して合成信号
を生成するとともに、第1および第2のディジタルセン
サ信号に基づいて利得信号を生成する。制御・測定シス
テムは、合成信号および利得信号をコントローラから受
け取るために接続された乗算ディジタル−アナログコン
バータを含み、合成信号および利得信号の積として駆動
信号を生成する。
【0024】制御・測定システムは、センサ信号に負の
利得を選択的に印加して、導管の運動を減らす。制御・
測定システムは、センサ信号のゼロオフセットを補償で
きる。ゼロオフセットは、利得変化に起因する成分と、
利得の非線形性に起因する成分とを含み、制御・測定シ
ステムは、別々に2つの成分を補償できる。制御・測定
システムは、一つ以上の訂正因子を生成するとともに、
訂正因子を使用してセンサ信号を修正することによっ
て、ゼロオフセットを補償する。
【0025】制御・測定システムは、第1および第2の
センサ信号に対する位相偏差を計算できる。位相偏差
は、センサ信号の基本周波数に相当するセンサ信号の成
分に対するゼロ位相の点と、センサ信号のゼロクロス点
との間の差として定義される。制御・測定システムは、
計算された位相偏差を合成して位相差を生成する。制御
・測定システムは、第1のセンサ信号の周波数を評価
し、第1のセンサ信号の周波数とは異なる第2のセンサ
信号の周波数を評価し、更に、評価された周波数を使用
してセンサ信号間の位相差を計算することによって、特
性の測定値を生成する。
【0026】センサが速度センサであるときに、制御・
測定システムは、センサ信号の周波数、振幅、位相を評
価し、速度センサおよび絶対位置センサ間の性能の差を
考慮するために、評価された周波数、振幅、位相を訂正
する。振動の見かけの振幅(すなわち、センサが速度セ
ンサのときは振動の速度)を制御する代わりに、システ
ムは、センサ信号を評価された周波数で割ることによっ
て、真の振幅を制御する。この訂正は、振幅の制御およ
びノイズの減少を改善する。
【0027】制御・測定システムは、センサ信号の第1
のパラメータを評価し、第2のパラメータの変化率を測
定し、測定された変化率に基づいて評価された第1のパ
ラメータを訂正する。例えば、システムは、導管の振動
の周波数や振幅の測定された変化率に基づいて、センサ
信号の評価された周波数、振幅、位相を訂正する。シス
テムは、各センサ信号に対して別々の訂正を実行でき
る。
【0028】ディジタル流量計は、質量流量計であって
もよく、導管を流れる材料の特性は質量流量となる。デ
ィジタル流量計は、比重計であってもよく、導管を流れ
る材料の特性は、材料の密度となる。制御・測定システ
ムは、初期の質量流量を測定し、導管を流れる材料の見
かけの密度を測定し、見かけの密度を材料の周知の密度
と比較して密度差を測定し、密度差に基づいて初期の質
量流量を調節して調節された質量流量を生成することに
よって、導管の空気混和の影響の原因を説明する。シス
テムは、調節された質量流量を調節して減衰の影響を考
慮することによって、導管の空気混和の影響の原因を更
に説明する。導管の空気混和の影響の原因を更に説明す
るために、システムは、第1および第2のセンサ信号の
振幅の差に基づいて、調節された質量流量を調節する。
【0029】振動可能な導管は、2つの平行な平面ルー
プを含む。センサおよび駆動回路は、ループの間に接続
される。計器は、単一のワイヤ対にのみ電力を受け取る
電力回路を含む。電力回路は、ディジタル制御・測定シ
ステムと、駆動回路とに電力を提供し、ディジタル制御
・測定システムは、単一のワイヤ対にのみ導管を流れる
材料の特性の測定値を伝えるように動作可能である。電
力回路は、ディジタル制御・測定システムに電力を提供
する定出力回路と、2つのワイヤから過剰の電力によっ
て充電される駆動コンデンサとを含む。ディジタル制御
・測定システムは、駆動コンデンサを放電して駆動回路
に給電し、さらに、駆動コンデンサの充電レベルをモニ
タして、コンデンサの充電レベルが閾値レベルに達した
後、駆動コンデンサを放電する。ディジタル制御・測定
システムは、周期的に駆動コンデンサを放電して、1対
のワイヤでの双方向通信を実行する。
【0030】制御・測定システムは、周期的信号の期間
において第1のデータセットを集めて、第1のデータセ
ットを処理して駆動信号と測定値とを生成する。システ
ムは、同時に、第1のデータセットを処理しながら、セ
ンサ信号の次の期間において第2のデータセットを集め
る。第1のデータセットに相当する期間は、第2のデー
タセットに相当する期間と重なることがある。
【0031】制御・測定システムは、駆動信号を制御し
てセンサ信号の振幅を一定の設定値に維持し、駆動信号
が第1の閾値レベルを越えたとき、一定の設定値を減ら
し、駆動信号が第2の閾値レベル未満であり、且つ一定
の設定値が最大許容値未満であるときは、一定の設定値
を増やす。第1の閾値レベルは、最大許容駆動信号の9
5%以下である。
【0032】制御・測定システムは、測定値に対する不
確定性分析を実行する。この場合、制御・測定システム
は、測定値と不確定性分析の結果とを制御システムに送
る。制御・測定システムは、駆動信号の位相を調整する
ためにディジタル処理を使用して、センサと駆動回路と
の間に接続された部品とセンサとに関連する時間遅延を
補償する。
【0033】
【好ましい実施例の記載】本発明の特徴および効果は、
図面を含む以下の記載と請求項とから明らかにする。図
1を参照すると、ディジタル質量流量計100は、ディ
ジタルコントローラ105、一つ以上の運動センサ11
0、一つ以上の駆動回路115、導管120(以下、流
管とも称す)、温度センサ125を含む。ディジタルコ
ントローラ105は、例えば、プロセッサ、フィールド
・プログラマブルゲートアレイ、ASIC、他のプログ
ラマブル論理アレイやゲートアレイ、プロセッサコアを
有するプログラマブル論理回路等の一つ以上を使用して
実行される。ディジタルコントローラは、運動センサ1
10から入力される信号に少なくとも基づいて、導管1
20を流れる質量流の測定値を生成する。ディジタルコ
ントローラも、駆動回路115を制御して、導管120
に運動を誘起する。この運動は、運動センサ110によ
り検出される。
【0034】導管120を流れる質量流は、駆動回路1
15によって供給される駆動力に応答して、導管に誘起
される運動に対応する。特に、質量流は、導管の温度と
同様に、運動の位相および周波数に対応する。ディジタ
ル質量流量計も、導管の中を流れる材料の密度の測定値
を提供する。密度は、運動の周波数および導管の温度に
対応する。記載されている技術の多くは、質量流の測定
よりはむしろ密度の測定を行う比重計に応用できる。
【0035】導管の温度は、温度センサ125を使用し
て測定され、その剛性および寸法等の導管の特性に影響
を及ぼす。ディジタルコントローラは、これらの温度作
用を補償する。ディジタルコントローラ105の温度
は、例えば、ディジタルコントローラの動作周波数に影
響を及ぼす。一般に、コントローラ温度の影響は、十分
に小さく、無視できると考えられている。しかしなが
ら、場合によっては、ディジタルコントローラは、固体
装置を使用して、コントローラの温度を測定して、コン
トローラの温度の影響を補償する。 A.機械的設計 1の実施例において、図2および図3に図示するよう
に、導管120は、小領域が導管に対して空間を除去し
たりまたは確保するパイプライン(図示せず)に挿入さ
れるように設計されている。導管120は、パイプライ
ンとの接続用の装着フランジ12と、パイプラインと垂
直に向きが定められた2つの平行な平坦ループ18,2
0を支持する中心マニホールドブロック16とを含む。
電磁気駆動回路46とセンサ48とは、ループ18、2
0の各端部の間に取付けられる。2つの駆動回路46の
各々は、図1の駆動回路115に相当するが、2つのセ
ンサ48の各々は、図1のセンサ120に相当する。
【0036】ループの反対側の端部の駆動回路46は、
大きさが等しく且つ符号が反対の電流(すなわち、位相
が180度ずれている電流)が給電されて、ループ1
8、20の直線領域26が同一平面の垂直二等分線56
を中心に回転され、これは、点P(図3)で、管を横切
る。駆動回路に供給される給電電流を繰り返し反転する
(例えばサイン波状に制御する)と、直線領域26の各
々は、線56−56、すなわちループの対称軸を中心
に、水平面内で弓ひも形を掃引する振動運動を被る。下
方の丸いターン部38、40でのループの横方向全体の
移動は、小さく、直径が1インチ(2.54cm)のパ
イプの2フィート(60.96cm)の長い直線領域2
6に対して1/16のインチ(1.6mm)のオーダで
ある。振動の周波数は、大抵約80〜90Hzである。 B.導管の運動 ループ18,20の直線領域の運動を、図4、図5、図
6の3つのモードに示す。図5に示す駆動モードにおい
て、ループは、対応する点Pを中心に180度位相がず
れて駆動され、故に、2つのループは、同期をとって、
反対に(inthe opposite sense)
に回転する。従って、AおよびCなどの、それぞれの端
部は、周期的に一緒になり、離れて行く。
【0037】図5に示す駆動運動は、図4に示すコリオ
リモードの運動を含む。これは、ループの間で反対方向
であり、互いにわずかに向かうように(または離れるよ
うに)直線領域26を移動させる。コリオリの効果は、
mvWに直接対応する。ここで、mはループの断面にあ
る材料の質量であり、vは質量が移動する速度(容積流
量)であり、Wはループの角速度(W=Wsinω
t)であり、mvは質量流量である。2つの直線領域が
サイン波状に駆動されて、各速度がサイン波状に変化す
るときに、コリオリの効果は最も大きい。これらの条件
の下に、コリオリの効果は、駆動信号に対して位相が9
0度ずれている。
【0038】図6は、同じ方向にループを変位させる望
ましくない共通モード運動を示す。このタイプの運動
は、図2および図3の実施例のパイプラインの軸方向の
振動によって生じることがある。何となれば、ループが
パイプラインに対して垂直であるからである。図5に示
されるタイプの振動は、非対称モードと呼ばれ、図4の
コリオリモードは、対称モードと呼ばれる。非対称のモ
ードの振動の固有周波数は、足のねじれの弾性力に依存
する。通常、図2および図3に示す形の導管に対する非
対称モードの共振周波数は、対称モードの共振周波数よ
りも高い。質量流測定値のノイズ感度を減らすために、
ある質量流量に対するコリオリ力を最大にすることが望
ましい。上記したように、ループは共振周波数で駆動さ
れ、コリオリ力は、ループが振動している周波数(すな
わち、ループの角速度)に直接対応している。それゆえ
に、ループは、非対称モードにおいて駆動され、これは
共振振動数をより高くする傾向がある。
【0039】他の実施例は、導管の設計が異なってい
る。例えば、単一のループや、直線チューブ領域を、導
管として使用してもよい。 C.電子工学的な設計 ディジタルコントローラ105は、ループの互いに対向
する端部に位置するセンサ48(すなわち、運動センサ
110)によって生じた信号を処理することによって、
質量流量を測定する。各センサによって、生じる信号
は、ループがセンサの隣りに位置する駆動回路によって
駆動される相対速度に相当する成分と、ループに誘起さ
れるコリオリ力に起因するループの相対速度に相当する
成分とを含む。ループは、非対称モードで駆動され、そ
の結果、駆動速度に相当するセンサ信号の成分は、振幅
が等しいが、符号が反対になる。その結果生じたコリオ
リ力は、対称モードにあり、故に、コリオリの速度に相
当するセンサ信号の成分は、振幅が等しく、符号も同じ
である。このように、信号の差を計算すると、コリオリ
の速度成分を相殺し、駆動速度と比例する差になる。同
様に、信号を合計すると、駆動速度成分を相殺し、コリ
オリ速度と比例する和になる。これは、その結果、コリ
オリ力と比例する。次に、この和は、質量流量を測定す
るために使用される。
【0040】1.アナログ制御システム ディジタル質量流量計100は、従来のアナログ質量流
量計に対してかなりの効果を奏する。後の議論のため
に、図7に、従来の質量流量計のアナログ制御システム
400を示す。センサ48の各々は、電圧信号を生成
し、信号VA0はセンサ48aにより生成され、信号V
BOはセンサ48bにより生成される。VA0、VB0
は、センサの位置での相対的なループの速度に相当す
る。処理の前に、信号VA0、VB0は、対応する入力
増幅器405、410で増幅されて、信号VA1、V
B1を生成する。増幅器およびセンサでのアンバランス
を補正するために、入力増幅器410は、同期復調器4
15及び積分器420を含むフィードバックループから
来るバランス信号によって制御される可変利得を有す
る。
【0041】増幅器405の出力部で、信号VA1は、
次式の形式を取る。
【0042】
【数1】 VA1=Vsinωt+Vcosωt そして、増幅器410の出力部で、信号VB1は、次式
の形式を取る。
【0043】
【数2】 VB1=−Vsinωt+Vcosωt 但し、VおよびVは、それぞれ駆動電圧およびコリ
オリ電圧であり、ωは駆動モードの角周波数である。電
圧VA1、VB1は、演算増幅器425によって差が計
算されて、VDRVを生成する。
【0044】
【数3】 VDRV=VA1−VB1=2Vsinωt 但し、VDRVは駆動運動に相当して、駆動回路に電力
を供給するために使用される。駆動回路に電力を供給す
ることに加えて、VDRVは、ゼロクロス検出器430
に供給される。ゼロクロス検出器430は、VDRV
周波数に相当する周波数を有する出力方形波FDRV
生成する。FDRVは、ディジタル位相同期ループ回路
435への入力として使用される。FDRVは、プロセ
ッサ440にも供給される。
【0045】電圧VA1、VB1は、演算増幅器445
により合計されて、以下に示すVCORを生成する。
【0046】
【数4】 VCOR=VA1+VB1=2Vcosωt 但し、VCORは、誘起されたコリオリ運動に対応す
る。VCORは、同期復調器450に供給され、同期復
調器450は、ゲーティング信号Qと周波数が異なり且
つ位相も異なるVCORの成分を拒絶することによっ
て、質量に正比例する出力電圧Vを生成する。位相同
期ループ回路435は、Qを生成する。このQは、V
DRVと同じ周波数(ω)を有し且つVDRVと位相が
90度異なる(すなわち、VCORと位相が一致する)
直交基準信号である。従って、同期復調器450は、ω
以外の周波数を拒絶し、故に、Vは、ωでVCOR
振幅に相当する。この振幅は、導管の質量に正比例す
る。
【0047】Vは、電圧・周波数コンバータ455に
供給される。電圧・周波数コンバータ455は、V
振幅に相当する周波数を有する方形波信号Fを生成す
る。次に、プロセッサ440は、FをFDRVで割っ
て質量流量の測定値を生成する。ディジタル位相同期ル
ープ回路435も、基準信号Iを生成する。この基準信
号Iは、VDRVと位相が同じであり、フィードバック
ループ制御増幅器410において、同期復調器415の
ゲート制御を行うために使用される。対応する入力信号
の駆動成分が乗算された入力増幅器405、410の利
得が等しいとき、演算増幅器445での加算動作は、信
号VCORにゼロ駆動成分を生成する(すなわち、V
DRVと同じ位相信号が無い)。入力増幅器405、4
10の利得が等しくないとき、駆動成分はVCORに存
在する。この駆動成分は、同期復調器415によって抽
出されて、積分器420によって積分されて、入力増幅
器410の利得を補正する誤差電圧を生成する。利得が
あまりにも高いときや、あまりにも低いとき、同期復調
器415は、出力電圧を生成し、この出力電圧により、
積分器によって、利得を修正する誤差電圧を変化させ
る。利得が所望の値に達するとき、同期変調器の出力
は、ゼロに行き、誤差電圧は、変化をやめて、利得を所
望の値で維持する。
【0048】2.ディジタル制御システム 図8は、ディジタル質量流量計100の一実施例500
のブロック図を示す。このディジタル質量流量計は、導
管120、駆動回路46、図2および図3のセンサ48
を、ディジタルコントローラ505とともに含む。セン
サ48からのアナログ信号は、アナログ−ディジタ
ル(”A/D”)コンバータ510によってディジタル
信号に変換されて、コントローラ505に供給される。
A/Dコンバータは、別々のコンバータとして、また
は、単一のコンバータの別々のチャネルとして構成して
も良い。
【0049】ディジタル−アナログ(”D/A”)コン
バータ515は、駆動回路46を駆動するために、コン
トローラ505からのディジタル制御信号をアナログ信
号に変換する。各駆動回路に対して別々の駆動信号を使
用することは、多数の効果を有する。例えば、システム
は、診断処理に対して、対称駆動モードと非対称駆動モ
ードとの間で容易に切り替わることができる。他の実施
例において、コンバータ515によって生成された信号
は、駆動回路46に供給される前に、増幅器によって増
幅される。また他の実施例において、単一のDAコンバ
ータを、両方の駆動回路に供給する駆動信号を生成する
ために使用してもよく、駆動信号は、駆動回路の一方に
供給される前に反転されて、導管120を対称モードで
駆動する。
【0050】高精度抵抗器520および増幅器525
は、各駆動回路46に供給される電流を測定するために
使用される。A/Dコンバータ530は、測定された電
流をディジタル信号に変換して、ディジタル信号をコン
トローラ505に供給する。コントローラ505は、駆
動信号を生成する際、測定された電流を使用する。温度
センサ535および圧力センサ540は、それぞれ、導
管の入口545及び出口550での温度と圧力とを測定
する。A/Dコンバータ555は、測定された値をディ
ジタル信号に変換して、ディジタル信号をコントローラ
505に供給する。コントローラ505は、様々な方法
で測定された値を使用する。例えば、圧力測定値の差
は、導管の背圧を測定するために使用される。導管の剛
性は、背圧にとともに変化するので、コントローラは、
測定された背圧に基づいて導管の剛性を説明する。
【0051】追加の温度センサ560は、A/Dコンバ
ータによって使用される水晶発振器565の温度を測定
する。A/Dコンバータ570は、この温度測定値を、
コントローラ505による使用のためにディジタル信号
に変換する。A/Dコンバータの入出力関係は、コンバ
ータの動作周波数に応じて変化し、動作周波数は、水晶
発振器の温度に応じて変化する。従って、コントローラ
は、温度の測定値を使用して、A/Dコンバータによっ
て供給されるデータや、システム較正のデータを調節す
る。
【0052】図8の実施例において、ディジタルコント
ローラ505は、図9に示すプロシージャ600に従っ
てA/Dコンバータ510によって生成されるディジタ
ル化されたセンサ信号を処理して、駆動回路46に供給
される駆動信号及び質量流の測定値を生成する。最初
に、コントローラは、センサからデータを収集する(ス
テップ605)。このデータを使用して、コントローラ
は、センサ信号の周波数を測定し(ステップ610)、
センサ信号からのゼロオフセットを除去し(ステップ6
15)、センサ信号の振幅および位相を測定する(ステ
ップ620、ステップ625)。コントローラは、これ
らの計算された値を使用して、駆動信号を生成し(ステ
ップ630)、質量流の測定値および他の測定値を生成
する(ステップ635)。駆動信号および測定値を生成
した後、コントローラは、新しいデータのセットを収集
して、このプロシージャを繰り返す。プロシージャ60
0の行程は、順番に、または並行に実行されたり、また
は順序を変化させて実行することができる。
【0053】周波数、ゼロオフセット、振幅、位相の関
係のため、1の推定値が、他を計算するときに使用され
る。これは、精度を改善するために、計算の繰り返しに
至る。例えば、センサ信号のゼロオフセットを測定する
際に使用される初期周波数決定は、オフセット・除去さ
れたセンサ信号を使用して修正される。更に、適切な場
合、周期に対して生成される値を、次のサイクルに対す
る開始推定値として使用できる。
【0054】図8は、ディジタル流量計に含まれる一般
的なハードウェアを示す。 a.データ収集 議論を簡単にするために、2つのセンサからのディジタ
ル化された信号を信号SV、SVとし、信号SV
は、センサ48aから出て、信号SVは、センサ48b
から出るものとする。新しいデータが絶えず生成される
にもかかわらず、計算は、両方のセンサの1つの全周期
に相当するデータに基づいていると仮定する。充分なデ
ータのバッファリングによって、この条件は、データを
処理する平均時間がデータの収集に要する時間よりも短
い限り、真である。周期のために実行される作業は、周
期が完了したことを決める行程と、周期の周波数(また
は、SVおよびSVの周波数)を計算する行程と、
SVおよびSVの振幅を計算する行程と、SV
よびSVの位相差を計算する行程とを含む。複数の実
施例において、これらの計算は、次に対するスタートと
して前の周期の終点を使用して、周期毎に繰り返され
る。他の実施例において、周期は、180度だけ、また
は他の角度(例えば90度)だけ重なり、故に、周期
は、それに先行して続く周期内に含まれる。
【0055】図10および図11は、信号SV、SV
からのサンプルデータの2つのベクトルを示し、それ
ぞれ、(sv1_in)および(sv2_in)と名付
ける。各ベクトルの第1のサンプリング点は、公知であ
り、ベクトルによって表されるサイン波のゼロクロスに
相当する。(sv1 in)に対し、第1のサンプリン
グ点は、負の値から正の値へのゼロクロスであり、(s
v2 in)に対し、第1のサンプリング点は、正の値
から負の値へのゼロクロスである。
【0056】周期に対する実際の開始点(すなわち、実
際のゼロクロス)は、サンプリング点とは正確にはまず
一致しない。このため、初期サンプリング点(star
sample SV1およびstart sampl
SV2)は、周期の開始直前に生じるサンプリング
点である。第1のサンプリング点と実際の開始との間の
差を説明するために、アプローチとして、周期が実際に
開始する次のサンプルと開始サンプルとの間の位置(s
tart offset SV1またはstart
ffset SV2)を使用する。
【0057】信号SVとSVとの間に位相偏差が存
在するので、(sv1 in)および(sv2 in)
は、異なるサンプリング点で開始する。サンプルレート
および位相差の両方が高い場合、(sv1 in)のス
タートと(sv2 in)のスタートとの間にはサンプ
ルの差が存在する。この差は、位相のオフセットの生の
推定値を提供し、計算された位相オフセットに対するチ
ェックとして使用できる。これを以下に説明する。例え
ば、55kHzでサンプリングするとき、1のサンプル
は、およそ0.5度の位相シフトに相当し、1のサンプ
ルは、およそ800のサンプル点に相当する。
【0058】コントローラが、例えば和(A+B)およ
び差(A−B)等の関数を使用して、BがAと同じ周波
数を持つように重みが付けられている時、さらなる変数
(例えば(satrt sample sum)および
(start offset sum)は、各関数に対す
る期間の開始を追跡する。和および差の関数は、SV
およびSVの中間に位相偏差を有する。
【0059】一実施例において、センサからのデータを
記憶するために使用されるデータ構造は、各センサ毎に
円形のリストであり、周期における最大数のサンプルの
少なくとも二倍の容量を有する。このデータ構造によっ
て、処理は、現在の周期に対するデータに対して実行さ
れ、一方、割り込みや他の技術が、次の周期に対するデ
ータをリストに加えるために使用される。
【0060】処理は、全周期に相当するデータに対して
実行されて、サイン波に基づいた方法を使用するときの
誤差を避ける。従って、周期に対するデータを組み立て
る際の第1の作業は、周期が始まり、終わるところを測
定することである。重畳しない周期が使用されるとき、
周期の開始は、前の周期の端部として識別される。重畳
する周期が使用されて、周期が180度重畳するとき、
周期の開始は、前の周期の中間点として、または前の周
期に先行する周期の終点として、識別される。
【0061】周期の終了点が、前の周期のパラメータに
基づいて、更にパラメータが周期から周期へと所定の量
以上変化しないという仮定の下に、最初に評価される。
例えば、5%が、最後の周期の値からの最大許容変化と
して使用され、これは、55kHzのサンプリングレー
トで、連続する周期に対して振幅または周波数の5%の
増加または減少の繰り返しが、1秒で5000%近くの
変化になるので、合理的である。
【0062】5%を振幅および周波数の最大許容増加と
して指定して、連続する周期において5度の最大位相変
化を許容することによって、信号SVに対する周期の
終了点での上限に対する従来の推定値が、以下のように
測定される。
【0063】
【数5】 但し、(start sample SV1)が(sv
in)の第1のサンプルであり、(サンプル レー
ト)はサンプリングレートであり、(est freq)
は前の周期からの周波数である。信号SVに対する周
期の終了端の上限(end sample SV2)
は、同じように測定される。
【0064】周期の終了端が識別されたあと、周期が処
理に値するか否かについての簡単なチェックが行われ
る。例えば、導管が失速したり、または、センサ波形が
著しくゆがめられたとき、周期は、処理に値しない。適
切な周期のみを処理すると、演算をかなり減らすことが
できる。周期の適性を判別する1つの方法は、周期のあ
る点を調べて予測された動作を確かめることである。上
記したように、最後の周期の振幅および周波数によっ
て、現在の周期に対応する値の有効な開始評価が得られ
る。これらの値を使用して、周期の30度、150度、
210度、330度に相当する点が調べられる。振幅お
よび周波数が前の周期の振幅および周波数に正確に匹敵
する場合、これらの点は、それぞれ(est amp/
2)、(est amp/2)、(−est amp/
2)、(−est amp/2)に相当する値を有す
る。但し、(est amp)は、信号の評価振幅(す
なわち、前の周期からの振幅)である。振幅および周波
数において5%の変化を許容することによって、不等式
が各四分の一周期毎に生成される。30度の点に対し
て、不等式は以下のように表せる。
【0065】
【数6】 他の点に対する不等式は、同じ形式を有し、(est
amp SV1)項の符号及びオフセット項(x/36
0)は、適切な値を有する。これらの不等式は、導管が
合理的な方法で振動していることをチェックするために
使用される。
【0066】測定値処理は、ベクトル(sv1 in)
(スタート:終了)、(sv2 in)(スタート:終
了)に生じる。但し、
【0067】
【数7】 となる。信号に対するスタート点および終了点間の差
は、信号の周波数を表す。
【0068】b.周波数測定 ディスクリートにサンプリングされた純粋なサインの周
波数は、パルス幅の間の遷移を検出し(すなわち、正ま
たは負のゼロクロスを検出することによって)、各期間
でのサンプル数を数えることによって、計算される。こ
の方法を使用することによって、例えば55kHzで8
2.2Hzサイン波をサンプリングすると、0.15%
の最大誤差で周波数を評価する。より高い精度は、例え
ば(start offset SV1)および(st
art offset SV2)を使用して、ゼロクロ
スが実際に起きたときのサンプルの一部を評価すること
によって得られる。ランダムノイズおよびゼロオフセッ
トによって、この方法の精度は減少する。
【0069】図12および図13に示すように、周波数
測定のより精密な方法は、サイン波の2乗の補間を使用
する。この方法によって、2次関数(quadrati
function)は、2乗したサイン波の最小点と
一致するようにはめ込まれ、2次関数のゼロは、周波数
を測定するために使用される。もし、次式を満たす場
合、
【0070】
【数8】 (但し、svが時刻tでのセンサ電圧であり、Aが振
動振幅、xが時刻tでのラジアン角(すなわち、x
=2πft)であり、δがゼロオフセット、εは分布
N(0,1)を有する確率変数であり、σはノイズの変
化である)、二乗関数は、次式で与えられる。
【0071】
【数9】 但し、xが2πに近いとき、sin(x)およびs
in(x)は、x0t=x−2π、およびx0t
としてそれぞれ近似される。従って、2πに近いx
の値に対して、aは、次式で近似できる。
【0072】
【数10】 これは、ノイズを加えた純粋な二次式(最小値はゼロ以
外δ=0を仮定する)であり、ノイズの大きさはσ及び
δに依存している。線形補間を使用してもよい。
【0073】この曲線あてはめ技術に関連した誤差の源
は、ランダムノイズ、ゼロオフセット、真の二次式から
の偏差である。曲線のあてはめは、ランダムノイズのレ
ベルに非常に敏感である。センサ電圧のゼロオフセット
は、正弦2乗関数のノイズの大きさを増やし、ゼロオフ
セット除去の重要性を示している(後述)。最小値から
移動すると、純粋なサイン波の2乗は、完全な2次式で
はない。最も重大な抽出項は、四次である。対照的に、
線形補間にとって最も重大な抽出項は、三次である。
【0074】この曲線あてはめ技術に関連する自由度
は、データ点のどれくらい、またはどのデータ点が使用
されるかに関係がある。最小数は3であるが、最小二乗
法の使用により、(より大きい演算費用で)それ以上を
使用することができる。かかる適合は、ランダムノイズ
に影響されない。図12は、2次近似が最小点から約2
0度離れたところまでは良いことを示している。最小点
からさらに離れたデータを使用すると、ランダムノイズ
の影響は減少するが、正弦2乗関数での2次以外の項
(すなわち、四次以上)に起因する誤差が増える。
【0075】図14に、曲線あてはめ技術を実行するた
めのプロシージャ900を示す。第1行程として、コン
トローラは変数を初期化する(ステップ905)。これ
らの変数は、(エンド ポイント),すなわちゼロクロ
スの最良推定値、(ep int)、すなわちエンド
ポイントに最も近い整数値、s〔0..i〕、すなわち
全サンプル点のセット、z[k]、すなわち(エンド
ポイント)に最も近いサンプル点の2乗、z〔0..n
−1〕、すなわち(エンド ポイント)を計算するため
に使用される2乗されたサンプル点のセット、n、すな
わち(エンド ポイント(n=2k+1))を計算するた
めに使用されるサンプル点の数、ステップ 長、すなわ
ちzの連続する値の間のsのサンプル数、反復 カウン
ト、すなわちコントローラが実行する繰り返しのカウン
トを含む。
【0076】次に、コントローラは、「エンド ポイン
ト」の第1の推定値を生成する(ステップ910)。コ
ントローラは、前の周期からの評価された周波数に基づ
いて評価されたゼロクロス点を計算し、評価されたゼロ
クロス点の周囲を(前後に)検索して真のクロスポイン
トに最も近いものを見つけることによって、この推定値
を生成する。次に、コントローラは、真のクロス点近傍
のサンプルのより小さな振幅を有するサンプル点に等し
い「エンド ポイント」を設定する。
【0077】次に、コントローラは、曲線あてはめのた
めの点の数、nを設定する(ステップ915)。コント
ローラは、11kHzのサンプルレートに対して5に等
しく、44kHzのサンプルレートに対して21に等し
いnを設定する。次に、コントローラは、「反復 カウ
ント」を0に設定し(ステップ920)、「反復 カウ
ント」を増加してプロシージャの繰り返しを開始する
(ステップ925)。
【0078】プロシージャの反復部分の第1行程とし
て、コントローラは、「反復 カウント」の値に基づい
て「ステップ 長」を選択する(ステップ930)。コ
ントローラは、「反復 カウント」がそれぞれ1、2、
3のいずれかに等しいかどうかに応じて、6、3、1に
等しい「ステップ 長」を設定する。次に、コントロー
ラは、「エンド ポイント」および0.5の合計の整数
部として「ep int」を測定し(ステップ93
5)、zアレイを満たす(ステップ940)。例えば、
nが5に等しいとき、z
〔0〕=s〔ep int−2
ステップ 長〕、z〔1〕=s〔ep int−ス
テップ 長〕、z〔2〕=s〔ep int〕、z
〔3〕=s〔ep int+ステップ 長〕、z
〔4〕=s〔ep int+2ステップ 長〕であ
る。
【0079】次に、コントローラは、例えばザビツキィ
・ゴレイ(Savitzky−Golay)フィルタ等
のフィルタを使用して、z〔k−1〕、z〔k〕、z
〔k+1〕の平滑な値を計算する(ステップ945)。
ザビツキィ・ゴレイフィルタは、ニューメリカルレシピ
(Numerical Recipes in C)第
650頁乃至第655頁(第2版、ケンブリッジ大学出
版、1995年)にプレス(Press)等により記載
されている。次に、コントローラは、二次式をz〔k−
1〕、z〔k〕、z〔k+1〕のあてはめ(ステップ9
50)、二次式(z)および対応する位置(x)の
最小値を計算する(ステップ955)。
【0080】xが(k−1)および(k+1)に相当
する点の間にある場合(ステップ960)、コントロー
ラは、xに等しい「エンド ポイント」を設定する
(ステップ965)。その後で、「反復 カウント」が
3未満の場合(ステップ970)、コントローラは、
「反復 カウント」を+1増加して(ステップ92
5)、プロシージャの反復部分を繰り返す。
【0081】xが(k−1)および(k+1)に相当
する点の間にある場合(ステップ960)、または「反
カウント」が、3に等しい場合(ステップ97
0)、コントローラは、プロシージャの反復部分を出
る。次に、コントローラは、「エンド ポイント」と周
期の開始点との間の差に基づいて周波数を計算する。そ
して、これは公知である(ステップ975)。
【0082】本質において、プロシージャ900によっ
て、コントローラは、各試みのより短い「ステップ
長」を使用して、「エンド ポイント」に戻る3つの試
みを実行する。試みに対する結果である最小値が、曲線
のあてはめに使用される点の外側に来る(すなわち、補
間より外挿になっている)場合、これは、前の推定値ま
たは新しい推定値のいずれかが貧弱であることと、ステ
ップサイズの減少が許可されないこととを示している。
【0083】プロシージャ900は、センサによって生
成される少なくとも3つの異なるサイン波に適用され
る。これらは、信号SV、SV、二つの重み付き合
計を含む。さらに、ゼロオフセットが除去されたと仮定
すると、これらの信号に対して生成された周波数の推定
値は、互いに独立している。これは、各々の誤差が独立
しているので、明らかに、信号SV、SVにあては
まる。しかし、信号SV、SV間での質量流および
対応する位相差が、周波数の計算に対して十分に大き
く、場合毎に異なるサンプルに基づいている限り、これ
は、重み付き和に対してあてはまる。これが正しいと
き、周波数推定値のランダム誤差も、互いに独立すべき
である。
【0084】3つの独立した周波数の推定値が構成され
て、改善された推定値を提供する。この合成された推定
値は、単に3つの周波数推定値の平均値である。 c.ゼロオフセットの補償 コリオリ送信機での重要な誤差の源は、各センサ電圧の
ゼロオフセットである。ゼロオフセットは、前置増幅回
路およびアナログ−ディジタルコンバータのドリフトに
よってセンサ電圧信号にもたらされる。ゼロオフセット
効果は、差動回路の使用による正及び負電圧に対する前
置増幅利得のわずかな差によって悪化することがある。
各誤差源は、送信機の間で変化し、送信機温度ととも
に、素子の摩耗による時間の経過とともに変化する。
【0085】コントローラにより使用されるゼロオフセ
ットの補償技術の例を、以下に詳細に記載する。一般
に、コントローラは、周波数評価及び積分法を使用し
て、センサ信号の各々のゼロオフセットを測定する。次
に、コントローラは、それらの信号からゼロオフセット
を除去する。信号SV、SVからゼロオフセットを
除去した後、コントローラは、それらの信号の周波数を
再び計算して、周波数の改善された周波数評価を行う。
【0086】d.振幅測定 振幅は、様々な可能性のある用途を有する。これらは、
フィードバックを介した導管振動の調整と、駆動回路波
形を合成するときのセンサ電圧の寄与の釣り合いと、位
相測定値に対する和および差の計算と、測定値を補正す
るための振幅変化率の計算とを含む。
【0087】1つの実施例において、コントローラは、
信号SV、SVの評価振幅を使用して、信号S
、SVの和及び差、和及び差の積を計算する。和
および差を測定する前に、コントローラは、一方の信号
を補正して、2つのセンサの利得間の差を考慮する。例
えば、コントローラは、信号SVの振幅に対する信号
SVの振幅の比に基づいて、信号SVに対してデー
タを補償するので、両方の信号が同じ振幅を持つ。
【0088】コントローラは、計算された和に基づい
て、周波数のさらなる推定値を生成する。これの推定値
は、以前の周波数の推定値と共に平均化されて、信号の
周波数の精製された推定値を生成したり、または以前の
推定値と交換したりする。コントローラは、フーリエに
基づく技術に応じて振幅を計算して、高調波の影響を除
去する。(ゼロクロス技術を使用して識別されるよう
な)期間Tに亘るセンサ電圧x(t)は、オフセットと
一連の高調波の項によって、以下のように表される。
【0089】
【数11】 この式によって、ゼロでないオフセットaは、ゼロで
ない余弦項aになる。対象の振幅は、基本成分の振幅
(すなわち、周波数ωでの振幅)であるが、高調波成分
の振幅(すなわち、周波数kω、但し、kは1より大で
ある)をモニタすることは、診断目的のための値になり
える。a、bの値は、次式で計算される。
【0090】
【数12】
【0091】
【数13】 各高調波の振幅は、次式から得られる。
【0092】
【数14】 積分は、2次式の補正を有するシンプソンの方法を使用
して計算される(後述)。この方法の主要な演算は、純
粋な正弦およ余弦関数を計算している。
【0093】e.位相測定 コントローラは、多数の方法を使用して、信号SV
SVの間の位相差を計算する。例えば、コントローラ
は、次式によって、t=0での開始時間に対する各高調
波の位相偏差を測定する。
【0094】
【数15】 位相偏差は、周期の開始(すなわち、ゼロクロス点)と
周波数ωのSV(t)の成分に対するゼロ位相の点との
間の差として、単一波形内で解釈される。位相偏差は、
波形全体に対する平均であるので、周期の中間点からの
位相偏差として使用される。理想的には、振動のゼロオ
フセットおよび一定振幅の無い状態で、位相偏差は周期
毎にゼロとなるべきである。コントローラは、同じ時間
に対する各センサ電圧の位相偏差を比較することによっ
て、位相差を測定する。
【0095】振幅および位相は、高調波の影響を除去す
るフーリエ法を使用して生成される。この方法は、導管
の両端部が同じ周波数で振動していることを仮定しない
という効果を有する。方法の第1行程として、周波数の
推定値は、ゼロクロスを使用して生成され、周期の開始
および終了の間の時間を測定する。周波数の線形変化が
仮定される場合、この推定値は、期間にわたる時間平均
周波数に等しい。評価され、仮定した周期の時間−不変
量、周波数ωを使用して、コントローラは、次式を計算
する。
【0096】
【数16】 但し、SV(t)は、センサ電圧波形(すなわち、SV
(t)、またはSV(t))である。次に、コント
ローラは、振幅および位相の推定値を判別する。
【0097】
【数17】 次に、コントローラは、各センサ信号の平均周波数及び
位相が波形全体を表しているものと仮定して、位相差を
計算する。これらの周波数はSV、SVに対して異
なるので、対応する位相は、平均周波数に拡大・縮小さ
れる。更に、位相は、同じ開始点(すなわち、SV
周期の中間点)へシフトされる。拡大・縮小された後、
減算されて位相差を求める。
【0098】
【数18】 但し、hは、サンプル長であり、中間点は、サンプルの
項において次式として定義される。
【0099】
【数19】 一般に、位相および振幅は、2つのセンサに対して同じ
時間フレームでは計算されない。流量がゼロであると
き、2つの周期の中間点は一致する。しかし、高い流量
で発散するので、計算は、時間内に一致しないサンプル
セットに基づいている。これは、変化する質量流の条件
下で、位相ノイズの増加につながる。十分な流量で、
(360度から)4度の位相シフトは、SV、SV
データセットの99%のサンプルだけが一致しているこ
とを意味する。空気混和の条件下では、はるかに大きな
位相シフトが観察され、これは、重複の低流量にさえ至
ることがある。
【0100】図15に、この問題に触れた修正された方
法1000を示す。最初に、コントローラは、SV
SVデータセットの周波数(f、f)と中間点
(d、d)とを見つける(ステップ1005)。最
後の周期から周波数に線形シフトを仮定して、コントロ
ーラは、SVの中間点(f2m1)でSVの周波数
と、SVの中間点(f1m2)でSVの周波数を計
算する(ステップ1010)。
【0101】次に、コントローラは、f1m2、f
2m1の周波数を仮定して、中間点m、mを用いで
新しいデータセット(d1m2、d2m1)の開始点と
終了点とを計算する(ステップ1015)。これらの終
了点が、必ずしもゼロクロス点と一致するわけではな
い。しかし、これは、フーリエに基づく計算に対する必
要条件ではない。次に、コントローラは、セットd
2m1の位相および振幅のフーリエ計算を実行する。
この位相差計算は、上記に概説した(ステップ102
0)。dおよびd2m1の中間点は同一であるから、
(スケール−シフト SV)は常にゼロであり、無視
できる。コントローラは、データセットd、d1m2
に対してこれらの計算を繰り返す(ステップ102
5)。次に、コントローラは、測定値を生成するため
に、計算された振幅および位相差の平均値を生成する
(ステップ1030)。中間点m、mの間に充分な
分離が存在するとき、コントローラは、二つの結果のセ
ットを使用して、位相および振幅の変化率の局所的な推
定値を提供する。
【0102】コントローラは、SV、SVの差を計
算し、計算した差を二乗し、その結果を積分することか
らなる差動増幅方法を使用する。他の方法によれば、コ
ントローラは、サイン波を合成し、信号SV、SV
の差をサイン波に乗算し、結果を積分する。コントロー
ラは、信号SV、SVの積も積分し、これは周波数
(2f)を有するサイン波となり(但し、fは信号SV
、SVの平均周波数である)、または、コントロー
ラは、積を二乗したり、その結果を積分したりする。コ
ントローラは、積のサイン波と同等のコサイン波を合成
し、合成したコサイン波に積のサイン波を乗算して、次
にコントローラが積分する周波数(4f)のサイン波を
生成する。コントローラは、これらの方法のうちの複数
を使用して、別々の位相測定値を生成し、次に、最終的
な位相測定値として別々の測定値の平均値を計算する。
【0103】差動・増幅法は、次式で始まる。
【0104】
【数20】 但し、φは、センサ間の位相差である。基礎的な三角恒
等式(trigonometricidentitie
s)が、次式の如く、信号間の和(Sum)および差
(Diff)を定義するために使用される。
【0105】
【数21】 これらの関数は、それぞれ、2Acos(φ/2)、
2Asin(φ/2)の振幅を有する。コントローラ
は、SV、SVに対するデータから(Sum)及び
(Diff)に対するデータセットを計算し、次に上記
方法のうちのいくつかを使用して、それらのデータセッ
トによって表される信号の振幅を計算する。次に、コン
トローラは、計算された振幅を使用して、位相差φを計
算する。
【0106】または、位相差は、次式にて定義される関
数(Prod)を使用して計算される。
【0107】
【数22】 これは、振幅(Asinφ)、周波数(2f)の関数
である。(Prod)は、サンプル毎に生成され、φ
は、結果としてのサイン波の振幅から計算される。
【0108】特に、位相の計算は、前の計算(すなわ
ち、SV、SVの周波数および振幅の計算)の精度
に依存している。コントローラは、位相の別々の(完全
に独立していなければ)推定値を提供する。これは、合
成されて、改善された推定値を与える。 F.駆動信号生成 コントローラは、利得を信号SV、SVの差に適用
することによって、駆動信号を生成する。コントローラ
は、正の利得(結果として正のフィードバック)または
負の利得(結果として負のフィードバック)のいずれか
を適用する。
【0109】一般に、導管のQは、十分に高く、故に、
導管は一定のディスクリートな周波数のみで共振する。
例えば、導管に対する最低共振周波数は、駆動周波数に
拘わらず、プロセス流体の密度に依存して、65Hzと
95Hzとの間である。このように、導管を共振周波数
で駆動してサイクル毎のエネルギ損失を最小にすること
は、望ましい。センサ電圧を駆動回路にフィードバック
すると、駆動周波数は共振周波数に移ることができる。
【0110】フィードバックを使用して駆動信号を生成
する他の実施例として、上記により測定された周波数及
び位相を有する純粋なサイン波が、合成されて駆動回路
に送られる。この方法は、共振周波数の高調波等の、不
要な高周波成分を除去する効果がある。この方法によっ
て、アナログ−ディジタルコンバータ、処理、ディジタ
ル−アナログ変換器によって導入される時間遅延を補償
でき、確実に、駆動信号の位相がセンサ信号の位相の中
間点に相当させる。この補償は、システム成分の時間遅
延を測定し、時間遅延相当する位相シフトを導入するこ
とによって、行われる。
【0111】導管を駆動する他の方法は、方形波パルス
を使用する。これは、別の合成方法であり、この方法
は、一定の(正及び負の)直流電源が所定の間隔でオン
・オフが切り替えられて必要なエネルギを提供する。切
り替えは、センサ電圧位相と同期している。この方法
は、ディジタル−アナログ変換器を必要としないので、
有効である。
【0112】一般に、導管の振幅は、開始時に所望の値
を高速で得るべきであるが、測定値機能を高速で行うた
めに、計器を傷つける可能性のあるかなりのオーバーシ
ュートの無い状態でそうすべきである。所望高速スター
トアップは、高利得を設定することによって成し遂げら
れ、故に、導管の高いQ及びランダムノイズの存在は、
導管の運動を始動させるためには十分である。1の実施
例において、高利得および正のフィードバックが、導管
の運動を始動させるために使用される。安定した動作が
得られると、システムは、駆動信号を生成するために、
合成方法に切り替わる。
【0113】図16乃至図27を参照すると、合成方法
は、高利得がそうできないときに導管運動を始動させる
ためにも使用される。例えば、センサ電圧の直流電圧オ
フセットがランダムノイズよりかなりより大きい場合、
高利得の適用は、振動運動を誘起しない。この状態を、
図16乃至図19に示す。図16乃至図19において、
高利得が、およそ0.3秒で適用される。図16および
図17に示すように、高利得を適用すると、一方の駆動
信号が、大なる正の値(図16)を仮定し、他方が大な
る負の値(図17)を仮定する。駆動信号の振幅は、セ
ンサ信号のノイズとともに変化する(図18および図1
9)。しかし、増幅されたノイズは、振動を誘起するた
めに駆動信号の符号を変化させるには不充分である。
【0114】図20乃至図23は、複数の周期における
方形波の印加によって、振動の高速スタートアップが生
じることができる様子を示している。直径が2インチの
導管の振動は、およそ2秒の間に安定する。導管振動の
安定は、図20および図21に示すように、駆動信号の
振幅の減少により示される。図24乃至図27は、1イ
ンチ導管がおよそ0.5秒間に安定する様子を示す。
【0115】方形波も、動作中に使用されて、導管振動
の問題を補正する。例えば、ある状況において、流量計
の導管は、例えば1.5kHzのオーダの周波数等の、
共振周波数の高調波で振動し始めることは公知である。
かかる高周波振動が検出されるときに、所望の周波数を
有する方形波は、導管振動を共振周波数に戻すために使
用される。
【0116】g.測定値の生成 コントローラは、アナログコントローラによって使用さ
れる方法と同様な方法で、質量流測定値をディジタル的
に生成する。コントローラは、密度等の他の測定値も生
成できる。一実施例において、コントローラは、2つの
センサ信号の間の位相差(phase diff)、導
管の振動周波数(freq)、プロセス温度(tem
p)基づいて、質量流を計算する。
【0117】
【数23】 但し、Tcは、較正温度であり、(MF−MF)は
較正プロシージャにおいて計算される較正定数であり、
(noneu mf)は非工学ユニットでの質量流であ
る。
【0118】コントローラは、導管の振動周波数および
プロセス温度に基づいて密度を計算する。
【0119】
【数24】 但し、(D−D)は較正プロシージャの間に生成さ
れる較正定数である。 D.積分法 多数の積分法が有効であり、様々な技術が、異なるレベ
ルの計算力を必要とし、異なるレベルの精度を提供して
いる。記載した実施例において、シンプソンの方法の変
形が使用される。基礎的な技術は、次式の如く表され
る。
【0120】
【数25】 但し、tはサンプルkでの時間であり、yは対応す
る関数値であり、hは行程長である。この規則は、奇数
個のデータポイント(すなわち、少なくとも3つの点)
で任意のデータベクトルに対して繰り返し適用され、デ
ータポイントの立方スプラインをあてはめて積分するこ
とと等価である。データポイントの数が偶然にも偶数の
場合、いわゆる第(3/8)番目規則が、間隔の1の終
了端に適用される。
【0121】
【数26】 以前に説明したように、各周期は、サンプリング点から
若干ずれたところ(例えば、start offset
SV1)で始まり、または終了する。積分法の精度
は、これらのオフセットを考慮することによって、かな
り改善される。例えば、半周期サイン波の積分におい
て、その結果の一貫した過小評価を避けるために、部分
サンプルに相当する領域を計算に含む必要がある。
【0122】2つのタイプの関数、すなわち、サイン関
数またはサイン二乗関数が、記載された計算において積
分される。両者は、終了点が起こるところでは、ゼロ近
傍に容易に近似される。終了点で、サイン波は、近似的
に線形であり、サイン二乗関数は、近似的に2次式であ
る。これら2つのタイプの関数の点から、3つの異なる
積分方法が評価された。これらは、終了点を補正しない
シンプソンの方法、線形の終了部を補正するシンプソン
の方法、2次式を補正するシンプソンの方法である。
【0123】積分法は、任意のアナログ−ディジタル切
り捨て誤差もシミュレーションすることなく、純粋なサ
イン関数およびサイン2乗関数をサンプリングすること
によって、テストされた。積分が計算され、その結果
は、信号の真の振幅と比較された。これらの計算の誤差
の唯一の源は、積分法に起因する。得られた結果を表A
および表Bに示す。
【0124】
【表1】
【0125】
【表2】 サイン関数に対し、線形補正を有するシンプソンの方法
は、最小の標準偏差で偏りが無く、一方、修正のないシ
ンプソンの方法は、負の誤差に偏り、2次補正を有する
シンプソンの方法は、比較的高い標準偏差を有した。サ
イン2乗関数に対して、誤差は一般に減少し、2次補正
は、最高の結果を提供した。これらの評価に基づいて、
サイン関数を積分するときは線形補正が使用され、サイ
ン2乗関数を積分するときは、2次補正が使用される。 E.同期変調法 図28は、センサ信号を処理する他のプロシージャ14
00を示す。プロシージャ1400は、同期変調に基づ
いている。これは、例えば、デニス(Denys)等に
よる、「同期性の損失に対するフランス未来防衛計画に
対する電圧位相の測定」(電力供給についてのIEEE
論文集、7(1)、62−69、1992)、ベゴビッ
ク(Begovic)等による「高調波が存在する場合
の電力網の周波数追跡」(電力供給についてのIEEE
論文集、8(2)、480−486、1993)に記載
されている。
【0126】第1に、コントローラは、システムの公称
動作周波数の初期推定値を生成する(ステップ140
5)。コントローラは、この公称周波数から、信号x
〔k〕(例えば、SV)の周波数の偏差を測定しよう
とする。
【0127】
【数27】 但し、Aは、信号のサイン波部分の振幅であり、ω
公称周波数(例えば、88Hz)であり、Δωは、公称
周波数からの偏差であり、hはサンプリングの間隔であ
り、φは位相シフトであり、ε(k)は、加えられたノ
イズおよび高調波に相当する。
【0128】この測定値を生成するために、コントロー
ラは、公称周波数で振動する2つの信号を合成する(ス
テップ1410)。信号は、0または(π/2)だけ位
相シフトされ、単一の振幅を有する。コントローラは、
これらの信号の各々に元の信号を乗算して、信号y
を生成する(ステップ1415)。
【0129】
【数28】 但し、y、yの第1項は高周波(例えば176H
z)成分であり、第2項は低周波(例えば0Hz)成分
である。次に、コントローラは、ローパスフィルタを使
用して、高周波成分を除去する(ステップ1420)。
【0130】
【数29】 但し、ε〔k〕およびε〔k〕は、元の信号からフ
ィルタ処理されたノイズを表す。コントローラは、これ
らの信号を合成して、u〔k〕を生成する(ステップ1
425)。
【0131】
【数30】 これは、周波数偏差に関する不可欠な情報を担持する。
示すように、u〔k〕は、u〔k〕の実数成分を表
し、u〔k〕は虚数成分を表す。
【0132】コントローラは、u〔k〕の実数および虚
数成分を使用して、周波数偏差、Δfを計算する(ステ
ップ1430)。次に、コントローラは、周波数偏差を
公称周波数に加算して(ステップ1435)、実際のΔ
f、fを与える。
【0133】
【数31】 コントローラは、u〔k〕の実数および虚数成分を使用
して、元の信号の振幅を測定する。特に、コントローラ
は、次式として振幅を測定する(ステップ1440)。
【0134】
【数32】 次に、コントローラは、2つのセンサ信号の位相差を測
定する(ステップ1445)。後述するローパスフィル
タの適用後に残っているノイズ(ε〔k〕およびε
〔k〕)を無視できると仮定すれば、y’〔k〕およ
びy’〔k〕のノイズフリーバージョン(y
〔k〕およびy 〔k〕)は、次式で表される。
【0135】
【数33】 これらの信号を乗算すると、次式となる。
【0136】
【数34】 0Hz近傍のカットオフ周波数を有するローパスフィル
タによって信号をフィルタ処理すると、不必要な成分が
除去されて、次式に示すものが残る。
【0137】
【数35】 この式から、位相差が、次のように計算される。
【0138】
【数36】 このプロシージャは、プロシージャがこの周波数から偏
差のみを測定するように、動作周波数が最初に評価され
る精度に依存する。よい評価が得られる場合、狭帯域フ
ィルタを使用でき、このことはプロシージャを非常に正
確なものにしている。典型的な流量計に対して、動作周
波数は、約95Hz(空)および82Hz(充填)であ
る。半分の範囲の第1の近似(88Hz)が使用され、
これによって、ローパスフィルタは13Hzのカットオ
フが可能になる。非常に小さいカットオフ周波数はサイ
ン波の振幅を減らすので、カットオフ周波数の選択は注
意しなければならない。
【0139】測定値の精度は、使用されるフィルタ処理
特性に依存する。不感帯でのフィルタの減衰は、高調波
の除去量を測定するが、カットオフ周波数が小さくなる
と、ノイズ除去を改善する。 F.PI制御を有する計器 図29および図30は、他の技術を使用して駆動回路に
供給する信号を生成するコントローラ1505を有する
計器1500を示す。アナログ−ディジタルコンバータ
1510は、センサ48からの信号をディジタル化し
て、ディジタル化された信号をコントローラ1505に
出力する。コントローラ1505は、ディジタル化され
た信号を使用して、駆動回路毎の利得を計算する。利得
は、導管に所望の振動を生成するために最適である。利
得は、正または負のいずれでも良い。次に、コントロー
ラ1505は、利得を、乗算ディジタル−アナログ変換
器1515に供給する。他の実施例において、直列に配
置した2つ以上の乗算ディジタル−アナログ変換器を、
単一のより感度のある乗算ディジタル−アナログ変換器
を実行するために使用してもよい。
【0140】コントローラ1505は、ディジタル化さ
れたセンサ信号を使用して、駆動信号も生成する。コン
トローラ1505は、これらの駆動信号をディジタル−
アナログ変換器1520に供給する。ディジタル−アナ
ログ変換器1520は、信号をアナログ信号に変換し、
このアナログ信号は、乗算ディジタル−アナログ変換器
1515に供給される。
【0141】乗算ディジタル−アナログ変換器1515
は、アナログ信号にコントローラ1505からの利得を
掛けて、導管を駆動する信号を生成する。次に、増幅器
1525は、これらの信号を増幅して、駆動回路46に
供給する。乗算ディジタル−アナログ変換器によって実
行される乗算を、コントローラ1505に実行させるこ
とによって、同様の結果が得られる。このとき、乗算デ
ィジタル−アナログ変換器は、標準のディジタル−アナ
ログ変換器に置換しても良い。
【0142】図30に、制御方法をより詳細に示す。コ
ントローラ1505内で、ディジタル化されたセンサ信
号が振幅検出器1550に供給される。振幅検出器15
50は、例えば上記の技術を使用して、導管の運動の振
幅の測定値a(t)を測定する。次に、加算器1555
は、振幅a(t)および所望の振幅aを使用して、次
式に示す誤差e(t)を計算する。
【0143】e(t)=a−a(t) 誤差e(t)は、比例積分(proportional
−integral;”PI”)制御ブロック1560
によって使用されて、利得K(t)を生成する。この
利得に、センサ信号の差が掛けられて、駆動信号を生成
する。PI制御ブロックは、変化している状態に対して
高速反応ができるようにする。振幅検出器1550、加
算器1555、PI制御ブロック1560は、コントロ
ーラ1505によって処理されるソフトウェアとして、
または分離回路として実行される。
【0144】1.制御プロシージャ 計器1500は、図31に示すプロシージャ1600に
応じて動作する。まず最初に、コントローラは、センサ
から、ディジタル化されたデータを受け取る(ステップ
1605)。その後で、プロシージャ1600は、並行
する3つのブランチ、すなわち、測定値ブランチ161
0、駆動信号生成ブランチ1615、利得生成ブランチ
1620を含む。
【0145】測定値ブランチ1610において、ディジ
タル化されたセンサデータが使用されて、前述のよう
に、振幅、周波数、位相の測定値を生成する(ステップ
1625)。次に、これらの測定値が使用され、質量流
量を計算し(ステップ1630)更に他の変数を計算す
る。一般に、コントローラ1505は、測定値ブランチ
1610を実行する。
【0146】駆動信号生成ブランチ1615において、
2つのセンサからのディジタル化された信号は、差が計
算されて、利得を掛ける信号を生成し、駆動信号を生成
する(ステップ1635)。上述のように、この差の生
成動作は、コントローラ1505によって実行される。
一般に、差の生成動作は、センサ信号間の振幅差を説明
する重み付き差を生成する。
【0147】利得生成ブランチ1620において、比例
・積分制御ブロックを使用して、利得が計算される。上
述したように、導管の運動の振幅a(t)は、測定され
て(ステップ1640)、所望の振幅aから引かれて
(ステップ1645)、誤差e(t)を計算する。別々
の行程として説明したが、振幅a(t)の生成は、測定
値生成ステップ1625の振幅の生成に相当する。最後
に、PI制御ブロックは、誤差e(t)を使用して、利
得を計算する(ステップ1650)。
【0148】計算された利得に、差信号が掛けられて、
駆動回路に供給する駆動信号を生成する(ステップ16
55)。上述のように、この乗算動作は、乗算DAコン
バータによって実行され、またはコントローラによって
実行される。 2.PI制御ブロック PI制御ブロックの目的は、振幅aを有する純粋なサ
イン波振動を導管に維持することである。導管の動作
は、簡単な質量−バネシステムとして次式に表すように
モデル化される。
【0149】
【数37】 但L、xは、平衡からの質量の変位及び時間の関数であ
り、ωは固有周波数であり、ζは減衰因子であり、小
さい(例えば0.001)と仮定される。出力y(t)
および入力i(t)の関数としてのこの運動方程式の解
は、電気回路と類似しており、この中で、供給電流i
(s)と検出された出力電圧y(s)との間の伝達関数
は、次式で表される。
【0150】
【数38】 導管に所望の振動を生成するために、利得K(t)を
有する正のフィードバックループが、「低速」外側ルー
プによって自動的に調節される。すなわち、
【0151】
【数39】 となる。システムが、「2−時間−スケール(two−
time−scales)」特性を有すると仮定する
と、これは、K(t)の変化が十分に低速なために、
上記で得られたxの方程式の解が、一定の減衰を仮定す
ることによって得られることを意味する。
【0152】ゼロ定常状態誤差を与える2−項PI制御
ブロックは、次式として表される。
【0153】
【数40】 但し、誤差e(t)(すなわち、a−a(t))は、
PI制御ブロックへの入力であり、K、Kは定数で
ある。1実施例において、a=10によって、K
0.02およびK=0.0005のコントローラの定
数は、振動が高速に形成される反応を提供する。しか
し、このPI制御ブロックは非線形であり、これは、結
果として、設計および動作の困難につながる。
【0154】振動振幅の動作の線形モデルは、x(t)
がAejwtに等しいと仮定することによって導かれ
る。そして、これは、次式になる。
【0155】
【数41】 これらの式をループの振動の式に代入し、実数項と虚数
項に分けると次式として表される。
【0156】
【数42】 A(t)も、次式で表すことができる。
【0157】
【数43】 この式の解を次式に示す。
【0158】
【数44】 a(t)を(logA(t))に等しいものとして定義
することによって変数を変形すると、A(t)の式は、
次式として表すことができる。
【0159】
【数45】 但し、Kは、明らかに時間に依存する。ラプラス変換
を行うと、以下に示す式になる。
【0160】
【数46】 それは、図32に示すように、伝達関数とみなせる。こ
の図は、全てのK及びaに対して線形であるので、で
コントローラの設計に対して特に大切であり、なお、唯
一の仮定は「2−時間スケール」特性である。閉ループ
の特性は、この仮定に関して強力であり、故に、実際に
達成しうる高速応答を、容易に設計することができる。
【0161】図32から、項(ζω)は、コントロー
ラにより除去されることが必要である「ロード外乱」
(すなわち、kK/2は、定数となるa(t)用のω
と等しくなければならない)である。ゼロ定常状態誤
差にたいして、これは、外側ループコントローラが積分
器(すなわち、非常に大きい利得)を有しなければなら
ないことを意味する。このように、適切なPIコントロ
ーラ(Cs)は、(K(1+1/sT))であると
仮定される。但し、Tは、定数である。比例項は、安
定性のために必要である。しかしながら、項(ζω
は、安定性やコントローラ設計に影響を与えず、これ
は、代わりに次式に示す開ループ伝達閏数に基づいてい
る。
【0162】
【数47】 を変化させる根の軌跡を、図33に示す。小さいK
に対して、ゆっくりと不十分に減衰する根がある。K
が増加するにつれて、根は、コントローラ利得が
(K)=8/(kT)となる点Pで実数になる。特に、
この理論は、Tの選択時にいかなる制限も置かない点
に注意すべきである。それゆえに、反応は、原則とし
て、KおよびTの適切な選択によって必要とされる
のと同様に高速で、精確に減衰される。
【0163】極点が、点Pで純粋に実数であるにもかか
わらず、これは、閉ループ行程反応にオーバーシュート
がないことを意味しない。これは、所望の値、a、誤
差eの間の伝達開数を調べることによって、容易に分か
る。
【0164】
【数48】 但し、Pは、2次の多項式である。ステップ入力(a
(s)=α/s)によっては、反応は、αp’(t)
と書くことができる。但し、p(t)は、1/p
(s)の逆変換であり、(aexp(−λt)+
exp(−λt))に等しい。信号p(t)は増
加し、次に、ゼロに減衰する。故に、e(t)はp’と
比例し、符号を変えなければならず、a(t)のオーバ
ーシュートを意味する。設定点aは、予めフィルタ処
理されて、次式に示す疑似設定値a を与える。
【0165】
【数49】 ここで、Tは、周知のコントローラパラメータであ
る。このプレフィルタを用いて、真のコントローラの極
点は、オーバシュートフリーの行程反応を提供すべきで
ある。この特徴は、オーバーシュート(例えば、機械的
妨害や部品の過剰ストレス)に対して物理的な制約があ
るので、有効である。
【0166】図33の根軌跡は、唯一の力学が、内側ル
ープの利得・振幅対数の伝達関数(図31)と、外側ル
ープのPIコントローラC(s)とである(すなわち、
振幅の対数(a=logA)が即座に測定される)こと
を仮定する。しかし、Aは、成長したりまたは減衰して
いる可能性のある振幅であり、故に、一般に、根本にあ
る正弦曲線を考慮せずに測定することはできない。上記
方法に加えて、Aを測定する方法は、複数存在する。準
安定状態での使用に適したものもある。例えば、サイン
波信号s(t)=sin(ωt+φ)が測定波形y
(t)=A(t)sin(ωt+φ1)にロックする
位相同期ループが、使用される。このように、振幅の測
定値、a=logAは、これらの信号を(適切な保護お
よびフィルタを用いて)掛けることによって得られる。
この方法は、恐らく定常状態近傍では適切であるが、ロ
ックする前の開始状態に対しては不十分である。
【0167】他の方法は、ピーク追跡器を使用し、この
ピーク追跡器は、ピーク追従アルゴリズムがコントロー
ラにおいて実行されるゼロクロス検出器を含む。しか
し、ゼロクロス方法は、ノイズの影響を受けやすい。更
に、ピーク追跡器からの結果は、半周期毎にのみ利用で
き、故に、コントローラの更新に対してサンプル間隔を
命令する。
【0168】最後に、AM検出器を使用できる。サイン
波をy(t)=Asinωtとすると、Aの推定値
が、
【0169】
【数50】 から得られる。但し、F{}は、1のDC利得を有する
適宜のローパスフィルタである。AM検出器は、最も簡
単な方法である。さらに、それは、いかなる特定周波数
の振動の存在も想定せず、それゆえに、スタートアップ
状態の間利用できる。これは、結果として生じる振動の
スペクトルに影響を及ぼす高調波のリークが、内側ルー
へと存在するという欠点に苦しむ。更に、フィルタは、
外側ループに余分な力学を加え、故に、反応速度とスペ
クトルの純度との間で妥協することが必要である。特
に、フィルタの効果は、最高のTの選択を強制するこ
とである。
【0170】abs(y)に対するフーリエ級数は、次
式で知られている。
【0171】
【数51】 このように、出力は、(π/2)だけ拡大・縮小され
て、正しい直流出力Aを与え、(偶数次の)高調波項
(acos2kωt)は、フィルタ処理される必要
がある。全フィルタが、DC成分を通過せしめて、他の
周波数の全てを減らす必要があるので、2ω以下のカ
ットオフを有する「レンガ−壁」フィルタは充分であ
る。しかし、フィルタの力学は、閉ループの動作に影響
を及ぼす。フィルタの共通選択は、バターワース(Bu
tterworth)形式にある。例えば、設計区切点
周波数ωを有する三次のローパスフィルタを次式で示
す。
【0172】
【数52】 設計周波数で、反応は3dB下がり、2ωで、−18
dB(0.12)であり、4ωで、−36dB(0.
015)下がる。高次バターワースフィルタは、鋭いロ
ールオフ(roll−off)を有するが、大部分の極
点は、複雑で、制御ループの根の軌跡に悪影響を及ぼす
ことがある。 G.ゼロオフセット補償 上記したように、ゼロオフセットは、前置増幅回路のド
リフトによって、さらにアナログ−ディジタルコンバー
タによって、センサ電圧信号にもたらされる。差動回路
の使用による正及び負電圧に対する前置増幅利得の若干
の差は、ゼロオフセット効果を悪化させることがある。
誤差は、送信機温度および部品の疲労によって、送信機
の間で変化する。
【0173】音質(すなわち、比較的低コスト)アナロ
グ−ディジタルコンバータが、経済的な理由により使用
される。これらの装置は、高い優先順位としてのDCオ
フセットおよび振幅安定性を伴って設計されない。図3
4乃至図37は、オフセットと正や負の利得とが、この
種のコンバータ(AD1879コンバータ)に対するチ
ップ動作温度によってどのように変化するかについてを
示す。図示された傾向の再現性は乏しく、この傾向に基
づいた温度補償を考慮に入れたとしても、残留ゼロオフ
セットおよび正・負利得は、残りのものと一致しない。
【0174】2つのセンサ電圧のゼロクロス点の時差を
使用して位相を計算する場合、DCオフセットは位相誤
差になる。この影響を、図38乃至図40に示す。各グ
ラフは、真の位相偏差がゼロ(すなわち、ゼロフロー)
のとき、ディジタル送信機によって測定されるような計
算された位相偏差を示す。図38は、正のゼロクロスで
開始する全周期に基づいて計算される位相を示す。平均
値は、0.00627度である。
【0175】図39は、負のゼロクロスで開始する、計
算された位相を示す。平均値は、0.0109度であ
る。図40は、半周期毎に計算された位相を示す。図4
0は、図38および図39からのデータにインターリー
ブする。平均位相(−0.00234)は、図38およ
び図39のものよりもゼロに近いが、信号の標準偏差
は、およそ6倍高い。
【0176】フーリエ法に基づいたものなどの、精巧な
位相測定値技術は、DCオフセットに不感である。しか
し、これらの技術が使用されるときでも、ゼロオフセッ
トを除去することが望ましい。何となれば、データは、
ゼロクロス点によって表現された全周期パケットにおい
て処理されるからである。これによって、例えば、見か
けの位相および周波数の振幅変調の効果の分析が簡単に
なる。更に、正及び負電圧間の利得の不整合は、誤差を
測定値法に含ませることになる。
【0177】位相検出のゼロクロス技術は、ゼロオフセ
ットおよび利得不整合誤差の効果と、その結果である除
去を説明するために使用される。図41および図42
は、ゼロフローを有する位相の長期ドリフトを示す。各
点は、有効なデータの1分間の平均値を表す。図41は
平均の位相を示し、図42は、位相の標準偏差を示す。
数時間にわたり、ドリフトは大きい。このように、毎日
計器をゼロに合わせた(大抵の用途において過剰な保守
条件と考えられている)としても、それでも、かなりの
位相ドリフトが存在する。
【0178】1.補償技術 電圧オフセットおよび利得不整合を扱う技術は、ディジ
タル送信機の計算能力を使用して、ゼロフロー条件を必
要としない。この技術は、周期毎に一組の計算を使用す
る。これは、合理的な期間(例えば10,000周期)
にわたり平均化されるときに、主たる変化の領域(例え
ば設定点の変化、空気混和の発生)を排除することで、
所望のゼロオフセットおよび利得不整合補償に収束す
る。
【0179】3次以上の高調波の存在を仮定すると、セ
ンサ電圧SV(t)に対する所望の波形は、次式の形式
を取る。
【0180】
【数53】 但し、Aは、基本周波数成分の振幅であり、A−A
は、3次高調波成分の振幅を示す。しかし、実際に、
実際の波形は、負および正の利得G、Gの間の不整
合とゼロオフセットZ(これは、ゼロ近傍の値を有す
る)とで混合されている。普遍性の損失なしで、G
1に等しく、Gは、次式によって与えられると、仮定
することができる。
【0181】 但し、εは利得不整合を表す。この技術は、振幅A
および周波数ωは一定であると仮定する。これは正し
い。何となれば、Zおよびεの推定値は、多数の周
期(例えば、約1分の動作で起きて、インターリーブさ
れた10000の周期)に対する平均に基づいているか
らである。この技術を使用するとき、コントローラは、
周波数及び振幅の相当量の変化の存在を調べて、分析の
有効性を保証する。高次の高調波が存在すると、特定の
高調波に対して位相および振幅情報を抽出するフーリエ
法の使用に至る。これは、SV(t)の積分と、変調サ
インまたはコサイン関数との乗算とを伴う。
【0182】ゼロオフセットは、関数の形式と同様に、
積分限界に影響を与える。ゼロオフセットが存在するの
で、振幅および位相の計算のための出発点は、周期的波
形SV(t)のゼロ位相点にない。ゼロオフセットZ
に対して、対応する位相偏差は、近似的に次式で表され
る。
【0183】
【数54】 小さい位相、φZ0=(−Z/A)に対して、対応
する時間遅延によって、
【0184】
【数55】 となる。積分は、拡大・縮小されて、制限値(すなわ
ち、Zおよびε方法ゼロとして)は、関係する高調
波の振幅に等しくなる。対象の第1の2つの積分は、次
式で表される。
【0185】
【数56】 これらの積分は、実際には、センサ電圧データの通常の
フーリエ解析の間に計算されているものを表す。下付き
文字「1」は、1次高調波を表す。
【0186】
【数57】 N及びPは、それぞれ、負または正の半周期を示し、s
およびcは、それぞれサインまたはコサインのいずれの
変調関数が使用されたを示している。
【0187】厳密に言えば、中央のゼロクロス点と、そ
れに対応する積分制限値とは、(π/ω+tZ0)より
は(π/ω−tZ0)によって与えられる。しかし、正
確なゼロクロス点よりは正確な中点を使用すると、簡単
な分析と、より優れた数値の動作(主にゼロクロス点の
位置の誤差による)とになる。正確な中点の使用によっ
て導入される唯一の誤差は、上記積分の各々の小領域に
誤った利得((1+ε)の代わりに1,その逆もあ
り)が掛けられることである。しかし、これらの誤差
は、Z εのオーダであり、無視できると考えられ
ている。
【0188】コンピュータの代数学を使用し、Zおよ
びεが小さいと仮定すると、積分に対する1次推定値
は、次式として導かれる。
【0189】
【数58】 和、差、積分とその推定値との比を含む有効且つ関数す
る関数が判別される。積分の和は、次式として表すこと
ができる。
【0190】
【数59】 一方、推定値の和は、等しい。
【0191】
【数60】 同様に、積分の差は、次式として表される。
【0192】
【数61】 一方、推定値の差は、次式で表される。
【0193】
【数62】 最後に、積分の比は、次式で表される。
【0194】
【数63】 一方で、推定値の比は、次式で表される。
【0195】
【数64】 対応するコサイン積分は、次式として定義される。
【0196】
【数65】 推定値、和は、次式で表される。
【0197】
【数66】 第2高調波積分は、次式で表される。
【0198】
【数67】 推定値は次式で表される。
【0199】
【数68】 和は次式で表される。
【0200】
【数69】 積分は、周期ごとに数値が計算される。以下に記載する
ように、さまざまな振幅、ゼロオフセット、利得に関し
て積分の値を評価する式は、再配置されて、計算された
積分に基づいて、ゼロオフセットおよび利得項の推定値
を与える。 2. 評価式の精度を、例で示す。基礎的な積分の各々に対し
て、3つの値が、提供される。すなわち、終了時の補正
を有するシンプソンの方法を使用した、積分の「真の」
値(ロンベルグ(Romberg)積分を使用したマス
キャド(Mathcad)内で計算される)、評価方程
式を使用している値、シミュレーションモードで作動す
るディジタル送信機によって計算される値である。
【0201】このように、例えば、次式により計算され
たIに対する値は、0.101353である。
【0202】
【数70】 一方、次式として計算された推定値(I
1Ps est)は、0.101358である。
【0203】
【数71】 シミュレーションモードのディジタル送信機を使用して
計算された値は、0.101340である。これらの計
算は、表Cに示すパラメータ値を使用する。
【0204】
【表3】 これらのパラメータ値を使用したことによる正確な評価
値及びシミュレーション結果を、表Dに示す。
【0205】
【表4】 このように、少なくとも選択された特定の値に対して、
1次方程式によって与えられた推定値は、極めて正確で
ある。Zおよびεがゼロに近づくにつれて、推定値
とシミュレーションとの両方の誤差は、ゼロに近づく。 3.実施例 積分に対する1次の推定値は、高調波の振幅、ゼロオフ
セット、利得不整合に関して、連立非線形方程式を定義
する。方程式が非線形であるので、正確な解は、容易に
利用できない。しかし、補正の繰り返しによる近似は、
限られた演算オーバーヘッドを合理的に収束する。
【0206】導管に固有の比が、A〜Aに対して仮
定される。このように、振幅A−Aの全でを計算す
る試みは、実行されない。その代わりに、AおよびA
だけが、上記で定義した積分方程式を使用して推定さ
れる。相対振幅の経験に基づいて、Aは、(A
2)として近似され、Aは、(A/10)として近
似される。ゼロオフセットの補償法は、図43に示すプ
ロシージャ2200に従って実行される。周期毎に、コ
ントローラは、積分I1Ps、11Ns、I1Pc、I
1Nc、I2Ps、I2Ns、関係がある関数sum
1s、ratio1s、sum1c、sum2sを計算
する(ステップ2205)。これは、周波数、振幅、位
相を測定するために使用される従来のフーリエ計算より
は、追加する必要な計算が少なくなる。
【0207】10000周期毎に、コントローラは、従
来の変化率評価法を使用して、センサ電圧振幅Aの傾
きをチェックする(ステップ2210)。振幅が一定で
あれば(ステップ2215)、コントローラは、ゼロオ
フセットおよび利得不整合用の計算に進む。このチェッ
クは、周波数安定性に対するテストまでのばすことがで
きる。
【0208】計算を実行するために、コントローラは、
最後の10000周期に対する関数(例えば、sum
1s)に対する平均値を生成する。次に、コントローラ
は、ゼロオフセットおよび利得不整合の最初の推定値を
作る(ステップ2225)。
【0209】
【数72】 これらの値を使用して、コントローラは、逆利得因子
(k)と振幅因子(amp factor)を計算する
(ステップ2230)。
【0210】
【数73】 コントローラは、逆利得因子および振幅因子を使用し
て、振幅の第1の推定値を作る(ステップ2235)。
【0211】
【数74】 次に、コントローラは、以下の計算によって推定値を改
善し、必要に応じて繰り返す(ステップ2240)。
【0212】
【数75】 コントローラは、標準の技術を使用して、Zおよびε
の値の収束を検査する。実際は、第1の繰り返しの後
の補正は小さく、経験上、3回繰り返せば十分である。
【0213】最後に、コントローラは、生のデータを調
整して、Zおよびεを除去する(ステップ224
5)。次に、コントローラは、プロシージャを繰り返
す。一旦、ゼロオフセットおよび利得不整合が生のデー
タから除去されると、Zおよびεに対する次の値を
生成する際に使用される関数(すなわち、sum1s
は、補正されたデータに基づいている。従って、Z
よびεに対するこれらの次の値は、残留ゼロオフセッ
トおよび利得不整合を反映して、以前に生成した値とと
もに合計されて、実際のゼロオフセットおよび利得不整
合を生成する。生のデータを調節する1つの方法におい
て、コントローラは、調整パラメータ(例えば、S1
off、S2 off)を生成する。この調整パラメー
タは、センサからのアナログ信号をディジタルデータに
変換する際に使用される。
【0214】図44乃至図46、図47および図48
に、プロシージャ2200を使用して得られた結果を示
す。短期間の動作を、図44乃至図46に示す。これ
は、スタートアップ後の5分間に得られた連続的な位相
推定値を示し、プロシージャ用の時間が、出力に影響を
及ぼし始めることができる。位相は、正のゼロクロス、
負のゼロクロス、その両方に基づいて示される。
【0215】正および負の平均値の差は、インターリー
ブされたデータセットの平均ゼロオフセットの減少と対
応して、20分の1に減少する。対応する標準偏差は、
およそ6分の1に減少した。長期間の動作を、図47お
よび図48に示す。最初の大きなゼロオフセットは、高
速で補正され、次に、位相偏差は、長時間にわたってゼ
ロ近傍に維持される。第1のわずかな値を排除した、平
均位相偏差は、6.14e−6であり、これは、プロシ
ージャが、電圧オフセットおよび利得不均衡の変化の補
償に成功していることを強力に示唆している。
【0216】ディジタルコリオリ計器に対するZ、ε
に対する値の多くは、信号SVに対しては、Z
(−7.923e−4)、ε=(−1.754
−5)でSあり、信号SVに対しては、Z=(−
8.038e−4)、ε=(+6.93e−4)であ
る。 H.力学的解析 一般に、コリオリ計器に対する従来の測定値計算は、導
管の各側面での振動の周波数および振幅が一定であり、
導管の各側面での周波数がいわゆる共振周波数と同じで
あり等しい、と仮定している。位相は、一般に導管の各
側面で別々に測定されず、2つの側面間の位相差は、測
定プロセスの間一定であると仮定している。デジタル計
器を使用した、半周期ごとの周波数、位相、振幅の正確
な測定値は、パラメータ値が約数秒間に平均化されると
きに、これらの仮定は単に有効であることを証明する。
100Hzまたはより高周波数でみると、これらのパラ
メータは、かなりの変化を示す。例えば、通常の動作
中、SVの周波数および振幅の値は、対応するSV
の値と負の相関を強く示す。従って、従来の測定値アル
ゴリズムは、これらのダイナミックな変化に対するノイ
ズの原因にさらされる。ノイズは、測定値計算速度が増
加するにつれて、より大きくなる。他のノイズの項が、
例えば流管力学、非線形性力学などの物理的な要因(例
えば振幅によって変化する流管の剛性)、または絶対位
置データより速度データを提供するセンサ電圧の動的な
結果によって導かれる。
【0217】記載の技術は、デジタル計器の高精度を引
き出して、導管の力学的な動作をモニタし且つ補償し
て、質量流および密度等のより正確なプロセス変数を提
供するために、ノイズを減らす。これは、周波数、位
相、振幅の変化率、流管力学、現実的な力学物理等のか
かる影響をモニタし且つ補償することによって、行われ
る。各側面での周波数が同じと仮定しない位相差計算
は、すでに記載した。他の補償技術を以下に記載する。
【0218】力学的な影響に対するモニタリングおよび
補償は、各センサレベルで生じて、位相、周波数、振
幅、または他のパラメータの補正済みの推定値を提供す
る。更なる補償は、導管のレベルで起きることがある。
ここでは、例えば位相差および平均周波数の計算の間
に、両方のセンサからのデータが合成される。これらの
2つのレベルは、広範囲の補償を行うために一緒に使用
してもよい。
【0219】このように、流量計による瞬間的な質量流
および密度の測定値は、流量計の動作の力学的な影響を
モデル化して考慮することによって、改良される。一般
に、コリオリ流量計の位相ノイズの80%以上が、測定
すべきプロセス条件よりも、流管力学(「鳴動」と呼ば
れることがある)を原因とすることがある。力学モデル
を適用すると、位相ノイズは4分の1から10分の1に
減少し、流量測定の性能の大幅な改善につながる。単一
のモデルは、流量および振動振幅の全てに有効である。
一般に、演算要件は、無視できる。
【0220】力学的分析は、互いに隔離されたセンサ信
号の各々で実行される。これは、導管の2つの側面間の
力学的相互作用のモデル化を回避する、または少なくと
も遅らせる。これは、各センサでの力学よりもはるかに
複雑になる傾向がある。また、個々のセンサ信号の分析
は、導管の2つの側面にプロセス流体から異なる力が作
用している空気混和やバッチスタートアップ等の状況に
おいて成功する傾向がある。
【0221】一般に、力学的分析は、これらのパラメー
タに対して計算された値の時間で変化する振幅、周波
数、位相の影響を考慮する。周波数および振幅は各セン
サ電圧の毎に容易に定義されるが、位相は、慣習上、セ
ンサ電圧間の差の点から定義される。しかし、フーリエ
解析を使用する場合、各センサに対する位相は、周期の
中間点と平均180度の位相点との差に関して定義され
る。
【0222】3つのタイプの力学的効果は、測定値誤差
と、いわゆる「フィードバック」効果と、「速度」効果
とである。振幅および位相を計算するアルゴリズムは、
周波数、振幅、位相が対象の時間間隔に亘り一定である
ことを仮定しているので、測定値誤差が生じる。測定値
アルゴリズムの特性は、これらのパラメータの変化を補
正することによって改良される。
【0223】フィードバック効果は、振動の一定振幅を
維持するために、エネルギを導管に供給して導管からの
エネルギ損失を埋め合わせるために生じる。エネルギを
導管に加える必要性は、振動振幅が所望の設定値からず
れ始めたあとにのみ、認識される。この結果、振動して
いる導管に対する運動方程式の減衰項は、ゼロではな
く、代わりに、ゼロ近傍で絶えず揺れている。導管の固
有周波数は変化しないが、これは、振幅の僅かな変化と
関連するゼロクロスのシフト(すなわち、位相変化)に
よってあいまいにされている。
【0224】センサ電圧は導管速度を観測するが、導管
の位置を表すものとして解析されるので、速度効果が生
じる。この結果は、振幅の変化率が見かけの周波数およ
び位相に影響することである。たとえ、これらのパラメ
ータの本当の値が一定であったとしてもである。 1.振幅変調に対するセンサ−レベル補償 力学的効果を補正する1つの方法は、センサ信号の振幅
をモニタし、振幅の変化に基づいて調整をする。力学的
効果を分析するために、位相、周波数、振幅の推定値
は、周期毎の各センサ電圧に対して測定される、と仮定
する。図49に示すように、計算は、完全であるが重複
する周期に基づいている。各周期は、ゼロクロス点で、
前の周期の中間点で始まる。正の周期は、最初のゼロク
ロスの直後に正の電圧で始まり、負の周期は負の電圧で
始まる。このように、周期nは正であり、周期(n−
1)および(n+1)は負である。ゼロオフセット補正
が実行されてゼロオフセットが無視できると仮定する。
また、高次の高調波が存在すると仮定する。
【0225】振幅、周波数、位相の線形変化を仮定す
る。この仮定の下に、1の周期の間の各パラメータの平
均値は、周期の中点でのパラメータの瞬間的な値に等し
い。周期は180度だけ重なるので、1の周期に対する
平均値は、次の周期の開始値と等しい。例えば、周期n
は、時間0から(2π/ω)までである。振幅、周波
数、位相の平均値は、中間点、(π/ω)での瞬間的な
値に等しく、この中間点は、周期(n+1)の開始点で
ある。周期(n+1)は、時間(π/ω)から(3π/
ω)までである。もちろん、これらのタイミングは近似
である。何となれば、ωも時間とともに変化するからで
ある。
【0226】a.力学的効果の補償プロシージャ コントローラは、図50に示すプロシージャ2600に
応じて力学的効果を考慮する。最初に、コントローラ
は、ゼロクロスを使用して、上記のように、周期の開始
および終了間の時間を測定することによって、周波数推
定値を生成する(ステップ2605)。周波数が線形に
変化すると仮定すれば、この推定値は、この期間におけ
る時間の平均周波数と等しい。
【0227】次に、コントローラは、推定された周波数
を使用して、上記のフーリエ法を使用しつつ、振幅およ
び位相の第1の推定値を生成する(ステップ261
0)。上記したように、この方法は、高調波の影響を除
去する。位相は、短一の波形において、周期(すなわち
ゼロクロス点)の開始と、位相偏差として表される、周
波数ωのSV(t)の成分に対するゼロ位相の点との差
になるように、解釈される。位相偏差は、全波形に対す
る平均であるので、周期の中間点からの位相偏差として
使用される。理想的には、振動のゼロオフセットおよび
一定振幅の無い状態で、位相偏差は、周期ごとにゼロで
あるべきである。しかし、実際には、それは、高レベル
の振動を呈し、質量流を補正して振幅の力学的変化を説
明する、優れた基礎を提供する。
【0228】次に、コントローラは、位相差を計算する
(ステップ2615)。多数の位相差の定義が存在し得
るが、この解析は、各センサ信号の平均位相及び周波数
が全波形を表しているものと仮定する。これらの周波数
はSVおよびSVに対して異なるので、対応する位
相は、平均の周波数に拡大・縮小される。更に、位相
は、同じ開始点(すなわち、SVの周期の中心点)へ
シフトされる。拡大・縮小した後に、これらが引かれ
て、位相差を提供する。
【0229】次に、コントローラは、周期nに対する振
幅の変化率を測定する(ステップ2620)。
【0230】
【数76】 この計算は、周期(n+1)からの振幅が、周期nの変
化率を計算するときに利用できると仮定する。生の振幅
計算が行われたあとの1の周期において補正が行われた
場合、これは可能である。変化率の正確な推定値を有す
ることの効果と、良好な測定値の補正とは、補正された
測定値が、一実施例においては5ミリ秒のオーダである
場合の、遅れよりも重要である。最新で生成された情報
は、常に導管の制御のため(すなわち、駆動信号の生成
のため)に使用される。
【0231】必要に応じて、振幅の補正が適用された
後、(下記に記載するように)変化率の補正された推定
値を計算してもよい。これは、振幅および変化率の最良
値に対する収束への反復になる。 b.フィードバックおよび速度効果に対する周波数補償 上記したように、フィードバックループの力学的側面
は、設定値を中心とする振幅の若干の偏位により、位相
に、時間変化するシフトを導く。これは、測定された周
波数に結果としてなり、これは、導管の固有周波数とは
異なるが、ゼロクロスに基づいている。速度センサが使
用される場合、位相にさらなるシフトが生じる。この追
加のシフトは、導管の位置振幅の変化とも関係してい
る。力学的分析は、これらの影響をモニタして補償す
る。従って、コントローラは、計算された振幅変化率を
使用して、周波数の推定値を補正する(ステップ262
5)。
【0232】導管の振動振幅を一定に維持するために使
用されるフィードバックループにおける振動導管の位置
は、次式として表される。
【0233】
【数77】 但し、θ(t)は、フィードバック効果によって生じる
位相遅延である。振動している導管の機械的Qは、大抵
1000のオーダである。これは、振幅および位相の若
干の偏差を意味する。これらの条件の下に、θ(t)
は、次式で与えられる。
【0234】
【数78】 各センサが速度を測定するので、
【0235】
【数79】 となる。但し、γ(t)は速度効果によって生じる位相
遅延であり、次式で与えられる。
【0236】
【数80】 導管の機械的Qは、大抵1000のオーダであり、それ
ゆえ、振幅および位相の変化が小さいから、次式の如く
仮定することは、理にかなっている。
【0237】
【数81】 これは、SV(t)に対する式を次式の如く簡単に表す
ことができることを意味する。
【0238】
【数82】 同じ理由により、速度オフセットの位相遅延に対する式
は、次式に示すように簡単になる。
【0239】
【数83】 フィードバックおよび速度効果の位相遅延を合計する
と、全位相遅延が得られる。
【0240】
【数84】 そして、SV(t)に対する式は、次式として表され
る。
【0241】
【数85】 これにより、実際の振動周波数は、振動の固有周波数と
は見分けることが可能である。前者を観測するが、後者
は密度計算に役立つ。合理的な時間長に対して、適切な
振幅制御を仮定すると、これらの2つの周波数の平均
は、同じである(何となれば、振幅の平均変化率はゼロ
でなければならないからである)。しかし、改善された
瞬間的な密度測定のために、力学的効果に対して実際の
振動周波数を補償して固有周波数を得ることは、望まし
い。これは、瞬間の密度が高速で時間とともに変化す
る、曝気された流体を扱う際に、特に有効である。
【0242】周期nに対して観測された見かけの周波数
は、周期(n−1)、(n+1)の中間点で起きている
ゼロクロスによって表される。速度変化による位相遅延
は、周期の見かけの開始および終了に影響する。
【0243】
【数86】 この分析に基づいて、次式に示す積分された誤差項(e
rror sum)を使用して修正が行われる。
【0244】
【数87】 但し、起動時の(error sum)の値(すなわ
ち、周期ゼロでの値)は、次式で与えられる。
【0245】
【数88】 これらの式は(1/8π)の値を有する定数項を含む
が、実際のデータが、(1/8π)の項がより適してい
ることを示している。この矛盾は、更なる解析により解
かれるモデル化されていない力学に起因する。
【0246】上記計算は、真の振動振幅、Aが有効であ
ることを仮定している。しかし、実際問題として、セン
サ電圧SVのみが観察される。このセンサ電圧は、次式
の如く表すことができる。
【0247】
【数89】 この式の振幅、amp SV(t)は、次式で表され
る。
【0248】
【数90】 この振幅の変化率は、次式で表される。
【0249】
【数91】 故に、以下の評価式を使用できる。
【0250】
【数92】 c.フィードバックおよび速度効果の周波数補償の適用 図51乃至図62は、プロシージャ2600の適用が、
直径が1インチ(2.54cm)の導管を有する計器か
らの真のデータに対して、固有周波数の推定値、さらに
プロセス密度をどのように改善するかについてを示す。
図の各々は、10000のサンプルを示し、これらは1
分間に亘り集められたものである。
【0251】図51および図52は、振幅の設定値に対
するランダムな変化が適用された時から取られた、SV
からの振幅および周波数データを示す。導管に水が充
填されて、流れが無いので、固有周波数は一定である。
しかし、観察された周波数は、振幅変化に応答して相当
量変化する。平均周波数値は、81.41Hzであり、
標準偏差は0.057Hzである。
【0252】図53および図54は、それぞれ、平均値
からの周波数変化と、プロシージャ2600を使用して
生成した補正項とを示す。偏差全体は、極めてよく整合
している。しかし、振幅変化に起因しないさらなる周波
数変化が存在する。図54に示す他の重要な特徴は、上
記のように、誤差項の適当な初期化の結果として、平均
がゼロに近いことである。
【0253】図55および図56は、生の周波数データ
(図55)を、補正関数を適用した結果(図56)と比
較する。平均値周波数に、無視できるシフトがあるが、
標準偏差は、4.4分の1に減少した。図56から、補
正された周波数データの残留構造があることは、明らか
である。観測された周波数に対する影響と周期における
位相の変化に基づいて、更なる分析は更なるノイズ減少
を生ずると予測される。
【0254】図57および図58は、平均周波数への対
応する影響を示す。これは、瞬間的なセンサ電圧周波数
の平均である。平均値周波数がプロセス流体の密度を計
算するために使用されるので、ノイズの減少(ここでは
5.2分の1)は、密度の計算に伝えられる。図59お
よび図60は、ランダムな振幅設定値にさらされている
直径が2インチ(5.08cm)の導管に対する生の平
均周波数と補正された平均周波数とを示す。2インチ
(5.08cm)の流管は、生および補正されたデータ
の両方に対して、1インチ(2.54cm)の流管より
も周波数変化が小さい。ノイズ減少度(noise r
eduction factor)は、4.0である。
【0255】図61および図62は、1インチ流管用の
真のフローデータを有する結果を示す。ランダム設定値
アルゴリズムは、通常の一定の設定値と置換される。こ
の結果、前の例よりも振幅変化が小さくなる。これは、
1.5のより小さなノイズ減少度になる。 d.振幅変調に対する位相測定の補償 再び図50を参照すると、コントローラは、次に位相測
定値を補正して、上記で行われる位相計算を仮定して、
振幅変調を考慮する(ステップ2630)。上記の位相
のフーリエ計算は、振動の振幅が、計算が行われている
データの周期の全体にわたって一定であると仮定してい
る。この領域は、データの周期に亘る振幅の線形変化を
仮定する補正を記載する。
【0256】高次の高調波を無視するとともに、任意の
ゼロオフセットが除去されたと仮定すると、センサ電圧
に対する式は、次式で与えられる。
【0257】
【数93】 但し、λは、時間による相対振幅変化に相当する定数
である。上記のように、積分IおよびIは、次式で
表される。
【0258】
【数94】 これらの積分を評価すると、次式になる。
【0259】
【数95】 これらの式を振幅の計算へ代入し、λに直列に展開す
ると、次式になる。
【0260】
【数96】 λが小さいと仮定し、第1項の後の項を全て無視する
と、これは、次式に示すように簡単になる。
【0261】
【数97】 これは、周期の中間点(t=π/ω)でのSV(t)の
振幅に等しい。従って、振幅計算は、補正のない必要と
される結果を提供する。
【0262】位相計算に対して、真の位相差および周波
数は一定であり、更に、電圧オフセットがない、すなわ
ち、位相値はゼロであるべきである、と仮定する。しか
し、振幅変調の結果として、生の位相データに適用して
振幅変調を補償する補正は、次式である。
【0263】
【数98】 括弧内の値は小さいと仮定すれば、逆タンジェント関数
は無視できる。より精巧な分析は、高次の高調波の影響
を考慮する。センサ電圧を次式として仮定する。
【0264】
【数99】 故に、全ての高調波の振幅は、周期の間に同時に同じ相
対速度で増加し、その結果の積分は、次式で表される。
【0265】
【数100】 更に、正の周期に対しては、
【0266】
【数101】 となり、負の周期に対しては、
【数102】 となる。振幅に対して、これらの式を計算に代入する
と、振幅計算は2次および高次の項においてのみ影響を
受け、その結果、補正は、振幅の1次近似には必要でな
いことが分かる。位相に対して、補正項は、正の周期に
対しては、次式に示すようになり、
【0268】
【数103】 負の周期に対しては、
【0269】
【数104】 となる。これらの補正項は、高次高調波の振幅の有効性
を仮定する。これらは、通常のフーリエ技術を使用して
計算されるが、高調波間の仮定した比を使用してこれら
すべてまたはそのうちのいくつかを近似することも、可
能である。例えば、直径が1インチ(2.54cm)の
導管の1の実施例に対して、振幅の比は、大抵、A
1.0、A=0.01、A=0.005、A
0.001である。
【0270】e.振幅変調補償の位相への適用 シミュレーションは、高次高調波および振幅変調のシミ
ュレーションを含み、デジタル送信機を使用して実行さ
れた。一実施例は、f=80Hz、A(t=0)=
0.3、A=0、A=0、A=0、λ=1e
−5*48kHz(サンプリングレート)=0.476
22を使用し、これは、高い振幅変化率に相当するが、
高次高調波は無い。理論は、−0.02706度の位相
偏差を提案する。1000周期に亘るシミュレーション
において、平均オフセットは、−0.02714度であ
り、標準偏差は2.17e−6である。シミュレーショ
ンと理論との差(シミュレーション誤差の約0.3%)
は、各周期の振幅の線形変化のモデルの仮定に起因して
いるが、シミュレーションは、振幅の指数関数的な変化
を生成する。
【0271】第2の実施例は、第2高調波を含み、f=
80Hz、A(t=0)=0.3、A(t=0)=
0.003、A=0、A=0、λ=−1e−6
48kHz(サンプリングレート)=−0.04762
2のパラメータを有する。本実施例に対し、理論は、正
または負の周期に対して、+2.706e−3、±2.
66%となる位相偏差を予測する。シミュレーションに
おいて、結果は、2.714e−3±2.66%とな
り、これも良く一致する。
【0272】図63乃至図67は、この補正が真の流量
計データをどの程度改善しているかについての例を示
す。図63は、直径が1インチの導管から集められ、一
定と仮定される低流量の、SVからの生の位相データ
を示す。図64は、上記の式を使用して計算された補正
因子を示し、図65は、その結果である補正された位相
を示す。最も明らかな特徴は、補正は位相信号の変化を
増やしたが、図66および図67に示すように、位相差
(すなわち、SV−SV)の標準偏差を1.26分
の1に減らしたことである。従って、特性は改善され
た。何となれば、この補正は、2つの位相間の相関を改
善して、位相差の変化の減少になったからである。この
方法は、他の流れの条件や導管の他のサイズに対して、
同様に機能する。
【0273】f.速度影響に対する位相測定値への補償 位相測定値計算は、速度効果によって影響を受ける。非
常に有効且つ簡単な補正因子(単位、ラジアン)は、次
式に示す形式をとる。
【0274】
【数105】 但し、ΔSV(t)は、振幅の相対変化率でり、次式
で表される。
【0275】
【数106】 但し、tは、ΔSV(t)が測定されている周期の
終了時間であり、tk+1は次の周期の終了時間であ
り、tk−1は前の周期の終了時間である。ΔSVは、
SVの変化率の推定値であり、その絶対値によって拡大
・縮小されて、SVの比例変化率となる。
【0276】図68乃至図72は、この方法を示す。図
68は、上記した振幅変調補正を適用した後に、単一の
センサ(SV)からの生の位相データを示す。図69
は、上記の式を使用して計算された補正因子(単位、
度)を示し、図70は、その結果である補正された位相
を示す。なお、補正された位相の標準偏差は、実際には
生のデータに対して増加している。しかし、対応する計
算が他方のセンサ(SV)に対して行われるとき、2
つの信号の位相の間の負の相関が増加する(−0.8か
ら−0.9まで)。結果として、生の位相測定値(図7
1)に基づいた位相差計算は、補正された位相測定値
(図72)よりもかなりのノイズを有する。図71およ
び図72の比較は、このノイズ減少法の効果を示す。図
72から、周期はおそらく調整が良くないポンプに左右
されるので、プロセス変数が減少し、測定値に多数の周
期が存在することは、明らかである。これのどれも、図
71の未補正の位相差データからは識別できない。
【0277】g.センサレベルノイズ減少の適用 上記の位相雑音低減法を組合せると、図73乃至図84
に図示するように、様々な流れの条件での瞬間的な位相
差測定値がかなり改良された。各グラフは、1インチの
導管に作用するデジタルコリオリ送信機によって、リア
ルタイムで同時に計算された3つの位相差測定値を示
す。中間帯域3600は、簡単な時間−差法を使用して
計算された位相データを示す。最も外側の帯域3605
は、上記のフーリエペースの方法を使用して計算された
位相データを示す。
【0278】フーリエ法は、はるかに多数のデータ、よ
り精巧な分析、より多くの演算効果を使用し、結果とし
て、ノイズが多い計算になる。これは、上記の力学的効
果に対するフーリエ法の感度を原因とする。データ36
10の最も内側の帯域は、センサレベルの雑音低減法を
適用した後の同じフーリエデータを示す。分かるよう
に、各グラフに示される標準偏差値により示されるよう
に、相当なノイズの減少が各場合も発生する。
【0279】図73は、流れが無く、導管が満たされ、
ポンプノイズの無いときの測定値を示す。図74は、流
れが無く、導管が満たされ、ポンプが動作しているとき
の測定値を示す。図75は、導管が空で濡れているとき
の測定値を示す。図76は、流量が低いときの測定値を
示す。図77は、流量が高いときの測定値を示す。図7
8は、高流量で、振動振幅が0.03Vのときの測定値
を示す。図79は、空気混和が小さい低流量の測定値を
示す。図80は、空気混和が大きい低流量の測定値を示
す。図81は、空気混和が小さい高流量の測定値を示
す。図82は、空気混和が大きい高流量の測定値を示
す。図83は、高流量移行への空の測定値を示す。図8
4は、移行を空にする高流量の測定値を示す。
【0280】3.流管レベルの力学のモデル化 力学モデルは、2つの基礎ステージに取り入れられる。
第1ステージにおいて、モデルは、システム識別の方法
を使用してつくられる。流管は、刺激されて、その力学
を明らかにし、一方、本当の質量流および密度の値は、
一定に維持される。流管の反応は、力学モデルを生成す
る際に、測定されて使用される。第2のステージにおい
て、モデルは、標準の流れのデータに適用される。流管
力学の効果の予測は、位相および周波数に対して行われ
る。次に、予測は、観察されたデータから引き算され
て、残留位相および周波数を残す。これは、プロセスに
よる。各ステージを以下に詳細に説明する。
【0281】a.システムの識別 システム識別は、水で満たされ、流れのない流管で始ま
る。振幅は、通常一定に維持されるが、0.05Vおよ
び0.3Vの間のランダムな設定値を割り当てることに
よって変化してもよい。通常は、0.3Vである。結果
であるセンサ電圧を図85に示す。図86および図87
は、それぞれ、計算された対応する位相および周波数値
を示す。これらの値は、周期毎に一回計算される。位相
および周波数の両者は、高度の(high degre
e)「構造」を示す。質量流相当する位相および周波数
が一定であるので、この構造は、流管力学に関係する傾
向がある。真の位相および周波数が一定であることが分
からないときに、この構造を予測する観察可能な変数
を、以下に表す。
【0282】最初に、上記したように、ΔSV(t
は、次式として表される。
【0283】
【数107】 この式は、ΔSVおよびΔSVを測定するために使
用できる。流管の位相は、Δに関係し、これは、(Δ
SV−ΔSV)として定義される。一方、周波数は
Δに関係し、これは、(ΔSV+ΔSV)として
定義される。これらのパラメータを、図88および図8
9に示す。図86を図88に、図87を図89と比較す
ると、Δと位相との関係、Δと周波数との関係が強
いことが分かる。
【0284】流管力学に対する補正は、位相や周波数か
ら適切な予測関数の倍数を引き算することによって得ら
れる。改善された結果は、次式に示す形式のモデルを使
用し
【0285】
【数108】 但し、y(k)は出力(すなわち、位相または周波数)
であり、uは予測関数(すなわち、Δ、またはΔ
である。システム識別の方法は、時間についての効果多
項式にある、次数n、m、係数a,bの値を提示す
る。y(k)の値は、周期毎に計算され、観察された位
相や周波数から引き算されて、残留プロセス値を得る。
【0286】力学補正の無い場合においても、ディジタ
ル流量計が長期に亘り非常に優れた精度を提供すること
を認めることは、重要である。例えば、1ロット200
kgを合計するとき、装置は、0.03%未満の再現性
を容易に達成する。力学のモデル化の目的は、力学精度
を改善することである。このように、生の値および補償
された値は、「変化」や「標準偏差」の減少以外は、類
似した平均値を有するべきである。
【0287】図90および図98は、生の周波数値と補
正された周波数値とを示す。平均値は殆ど同じである
が、標準偏差は3.25分の1に減少している。周波数
の全体の偏差は除去されたが、かなりの「構造」が、残
留ノイズに残っている。この構造は、Δ関数とは無関
係に見える。使用されるモデルは、簡単な1次モデルで
あり、ここでは、(m=n=1)である。
【0288】図91および図99は、対応する位相補正
を示す。平均値への影響はかなり小さいが、標準偏差
は、7.9分の1に減少する。モデルの次数は、n=
2、m=10である。残留ノイズに残っているように見
える構造もある。この構造は、設定値の変化によって、
位相力学の不充分な励振による。より有効な位相識別
は、データ収集中の流管の連続打(セットポイントの変
化は、未だ実行される)によって、流管力学の更なるシ
ミュレーションによって行われる。図92および図10
0は、これらの条件下での補正の効果を示す。図示する
ように、標準偏差は、31分の1に減少する。このより
有効なモデルは、以下の記載において使用される。
【0289】b.フローデータへの応用 識別されたモデルの本当のテストは、新しいデータを提
供する改良である。初めに、多数の観察に注目すること
は、有効である。第一に、例えば10秒以上の間で平均
化された平均位相は、すでにとても正確である。図示し
た実施例において、位相値は、82Hzまたはその近傍
にプロットされている。報告された標準偏差は、10H
zに平均化されたときに示される値のおよそ1/3であ
り、1Hzに平均化されるときに示す値の1/9であ
る。基準として、1インチ流管に関しては、1度の位相
差は、約1kg/sの流量に相当する。
【0290】この方法の予想される利点は、平均精度を
改善するよりは、真のプロセス変化に対してより優れた
力学反応を提供することである。従って、流れがゼロで
はない以下の実施例において、小なる流れのステップ変
化が、補正された位相がよりはっきりとステップ変化を
示すという予測によって、10秒ごとに導かれる。図9
3および図101は、スタートアップ直後に、流れのな
く且つ充填されている流管に適用される補正を示す。ス
タートアップのリングアップ効果の特徴は、生のデータ
では明らかである(図93)、しかし、これは、補正に
より除去され(図101)、故に、標準偏差はデータセ
ット全体に対して23分の1に減少される。補正された
測定値は、ホワイトノイズに似ていて、大抵の流管力学
が取り込まれたことを示している。
【0291】図94および図102は、「排水された」
流管に対する補正を示す。ノイズは、6.5分の1に減
少する。しかし、ノイズに若干の残留構造があるように
見える点に注意すべきである。低流量、中間流量、高流
量に対する方法の効果(図95および図103、図96
および図104、図97および図105)も、10秒ご
との流れのステップ変化毎に示す。各事例において、パ
ターンは同じである。すなわち、補正された平均の流れ
(図103乃至図105)は生の平均の流れ(図95乃
至図97)と同じであるが、力学的ノイズはかなり減少
する。図405において、これは、以前にノイズでおお
い隠されていたステップ変化(図97)の出現になる。
【0292】4.力学的モニタおよび補償技術の拡張 前の章は、力学的動作(センサおよび流管レベルの両方
での流管力学、速度効果、振幅変調によって生じる周波
数および位相ノイズ)の異なる状態をモニタし補償する
ために使用される様々な方法(物理的なモデリング、シ
ステム識別、発見的教授法)を記載した。自然な拡張に
よって、古典的なモデリングおよび識別法と同じよう
に、人工知能、神経網、ファジー論理、遺伝子のアルゴ
リズムのそれらを含む、制御及び計器の実務家に周知の
同様の方法が、計器の力学の性能の状態側面に適用され
る。特に、これらは、流管レベルでの平均周波数および
位相差と同様に、センサレベルでの周波数、振幅、位相
変化のモニタリングおよび補償を含む。何となれば、こ
れらの変化は、測定値間隔の間の時間(測定値間隔が重
ならない)と同様に、各測定値間隔内に起こるからであ
る。
【0293】この方法は、プロセス測定値変化に対する
ノイズ及び力学的反応の両方を減らすとともに改良する
ときに、普通でない。このように、この方法は、流れの
測定値の内容において非常に価値があると見込まれてい
る。 I.空気混和 ディジタル流量計は、導管に空気混和が存在する場合に
は、改善された性能を提供する。空気混和は、質量流量
計によって生成された測定値に相当な負の影響がある導
管のエネルギ損失を生成し、結果として導管が失速す
る。実験によって、ディジタル流量計は、従来のアナロ
グ流量計に対して、空気混和がある場合には性能が実質
的に改善されることを示した。この性能改善は、非常に
広い利得範囲を提供し、負のフィードバックを使用し、
非常に低い振幅レベルで正確に測定値を計算し、振幅の
変化率および流管力学などの力学的効果を補償する計器
の能力から生じ、更に、正確なデジタル振幅制御アルゴ
リズムの計器の使用から生じる。
【0294】ディジタル流量計は、必要な駆動回路利得
が、見かけの流体密度の低下と同時に上がるときに、空
気混和の発生を検出する。次に、ディジタル流量計は、
検出された空気混和に直接反応する。一般に、計器は、
導管を流れる材料の観察密度(すなわち、通常の測定法
から得られる密度測定値)を、材料の空気混和のない周
知の密度と比較することによって、空気混和の有無をモ
ニタする。コントローラは、観察された密度と実際の密
度との間の差に基づいて、空気混和のレベルを測定す
る。次に、コントローラは、質量流の測定値を補正す
る。
【0295】コントローラは、空気混和が存在しない期
間(すなわち、密度の値が安定している期間)での密度
をモニタすることによって、材料の空気混和のない密度
を測定する。あるいは、コントローラが接続されている
制御システムが、接続されて、初期化パラメータとして
空気混和のない密度を提供してもよい。一実施例におい
て、コントローラは、空気混和の影響を考慮する3つの
補正、泡沫作用補正、減衰作用補正、センサ非平衡補正
を使用する。図406乃至図113は、補正プロシージ
ャの効果を示す。
【0296】図106は、測定された密度が異なる質量
流量、空気混和補正のない状態で減少する(すなわち、
空気混和が増加する)ときの、位相測定値の誤差を示
す。示すように、位相誤差は、負であり、振幅は、空気
混和の増加と共に増加する。図107は、結果として生
じた質量流誤差が負であることを示している。ディジタ
ル流量計が、高いレベルの空気混和で作動する点に注目
することは、大切である。対照的に、水平線4000で
示すように、従来のアナログ計器は、低レベルの空気混
和が存在するときに、失速する傾向がある。
【0297】流量計は、低振幅の振動で高駆動電流を流
す十分に大きな駆動回路利得を提供することができない
ときに、失速する。減衰のレベルが、一定の振幅で振動
を維持するために流管によって供給できるもの以上の高
い駆動回路利得を必要とする場合、導管に供給される駆
動エネルギは不十分である。これは、振動振幅の低下に
つながり、これは、最大利得限界値により供給されるよ
り僅かな駆動エネルギになる。破局的な崩壊の結果と、
流管の振動とは、対応する駆動回路利得必要条件が流量
計によって供給されるレベルまで減衰が減少するまで、
不可能である。
【0298】泡効果補正は、無効部分と呼ばれる空気混
和のレベルが増加するにつれて質量流が減少するという
仮定に基づいている。無効な部分と泡沫効果との実際の
関係を予測しようと試みることなく、この補正は、優れ
た理論的な正当化によって、観察された質量流への影響
は、観察された密度への影響と同じであると、仮定す
る。本当の流体密度が公知であるので、泡効果補正は、
同じ割合で質量流量を補正する。この補正は、全ての流
量に対して線形調整である。図108および図109
は、それぞれ、泡沫効果を補正した後の、残留位相およ
び質量流誤差を示す。図示するように、残留誤差は、正
であり、元の誤差よりも振幅がかなり小さい。
【0299】減衰因子補正は、空気混和による導管運動
の減衰を考慮する。一般に、減衰因子の補正は、観察さ
れた位相φobsと、実際の位相φtrueとの次式に
示す関係に基づいている。
【0300】
【数109】 但し、λは減衰係数であり、kは定数である。図110
は、異なる質量流量と異なるレベルの空気混和とに対す
る減衰補正を示す。図111が、減衰補正後の残留位相
誤差を示す。図示するように、位相誤差は、泡沫影響の
補正後に残る位相誤差に対してかなり減少する。
【0301】センサバランス補正は、導管の別々の端部
の間の密度差に基づいている。図114に示すように、
導管の入口と出口との間の圧力低下によって、入口から
出口までの泡の大きさが増加する。材料は、導管の2つ
のループを経由して連続的に流れるので、導管の入口側
(すなわち、第1のセンサ・駆動回路対に隣接する側)
での泡は、導管の出口側(すなわち、第2のセンサ・駆
動回路対に隣接する側)での泡よりも小さい。この泡サ
イズの差は、導管の2つの端部の間の質量及び密度の差
になる。この差は、センサ信号(SVおよびSV
に反映される。従って、センサバランス補正は、2つの
センサ信号の比に基づいている。
【0302】図112は、異なる質量流量と異なるレベ
ルの空気混和とに対するセンサバランス補正を示す。図
113は、センサバランス補正を適用した後の残留位相
誤差を示す。低流量且つ低レベルの空気混和で、位相誤
差は、減衰補正後、残留位相誤差に対して改善される。
他の補正因子も使用できる。例えば、各センサ信号の位
相角がモニタされる。一般に、信号に対する平均位相角
度は、ゼロであるべきである。しかし、平均の位相角度
は、空気混和が増加すると共に増加する傾向がある。従
って、補正因子は、平均位相角度の値に基づいて生成さ
れる。他の補正因子は、導管の温度に基づいている。
【0303】一般に、補正因子の適用は、質量流誤差を
1パーセント以下に保つ傾向がある。さらに、これらの
補正因子は、広範囲にわたる流れおよび空気混和レベル
に対して適用できるように見える。 J.設定値の調整 ディジタル流量計は、導管の振動振幅に対する設定値の
改善された制御部を提供する。アナログ計器において、
フィードバック制御が、導管の振動振幅を、所望のピー
クセンサ電圧(例えば0.3V)に相当する一定レベル
に維持するために使用される。安定した振動振幅は、周
波数および位相測定値の変化の減少につながる。
【0304】一般に、大きな振動振幅が望ましい。何と
なれば、かかる大きな振幅は、測定目的に対して大きな
コリオリ信号を提供するからである。大きな振動振幅
は、導管での高レベルのエネルギ保存になり、これは、
外部振動に対する耐性をより大きくする。駆動回路に供
給される電流の制限のために大きな振動振幅を維持する
ことが不可能な状況が、生じる。例えば、アナログ送信
機の一実施例において、電流は、安全性のために100
mAに制限される。これは、所望の振動振幅を維持する
ために必要な電流の5〜10倍である。しかし、プロセ
ス流体が(例えば、二相の流れを経て)更にかなり減衰
する場合、最適の振幅は、もはや持続できない。
【0305】同様に、後述する2ワイヤ計器などの低電
力流量計は、導管を駆動するために利用できる電力が非
常に少ない。更に、電力レベルは、導管が容量性放電に
よって駆動されるときに変化する。図115を参照する
と、ディジタル流量計のコントローラにより実行される
制御プロシージャ4200は、最大有効電流レベルが与
えられる最高の維持可能な設定値を選択するために使用
される。一般に、プロシージャは、所望の駆動電流出力
が選択される毎に実行され、そして、これは大抵周期毎
に1回である。または、インターリーブされた周期が使
用される場合、半周期毎に1回である。
【0306】コントローラは、設定値をデフォルト値
(例えば0.3V)設定して、センサ電圧(filte
red SV)および駆動電流(filtered
C)のフィルタ処理された式を初期化することによって
スタートする。このプロシージャが実行される毎に、コ
ントローラは、センサ電圧(SV)および駆動電流(D
C)に対する現在の値に基づいて、フィルタ処理された
値を更新する(ステップ4210)。例えば、コントロ
ーラは、99パーセントの(filtered SV)と
1パーセントのSVとの和として、(filtered
SV)に対する新しい値を生成する。
【0307】次に、コントローラは、プロシージャが休
止されて以前の設定値調整のために時間を提供して影響
を及ぼしたか否かを判別する(ステップ4215)。プ
ロシージャの休止は、ゼロを越える値を有する休止周期
カウントによって示される。プロシージャが休止される
場合、コントローラは、周期に対して更なる動作を実行
せず、休止周期カウントを減らす(ステップ422
0)。
【0308】プロシージャが休止しなければ、コントロ
ーラが、フィルタ処理された駆動電流が閾値レベルを越
えるか否かを判別する(ステップ4225)。一実施例
において、閾値レベルは、最大有効電流の95%であ
る。電流が閾値を越える場合、コントローラは設定値を
減らす(ステップ4230)。設定値の変化の後に計器
に対する時間を決めるために、次に、コントローラは、
休止周期カウントを適切な値(例えば100)と等しく
設定することによって、プロシージャの休止を実行する
(ステップ4235)。
【0309】プロシージャが休止しない場合、コントロ
ーラは、フィルタ処理された駆動電流が閾値レベル未満
であるか否かと(ステップ4240)、設定値が最大許
容設定値未満であるか否かと(ステップ4245)を判
別する。一実施例において、閾値レベルは、最大有効電
流の70%に等しい。両方の条件が一致する場合、コン
トローラは、可能な新しい設定値を測定する(ステップ
4250)。一実施例において、コントローラは、新し
い設定値を、(filtered DC)に対する(f
iltered SV)の比で乗算した最大有効電流の
80パーセントとして測定する。設定値の僅かな変化
(すなわち、揺らぎ)を避けるために、次に、コントロ
ーラは、可能な新しい設定値が現在の設定値を相当量だ
け越えるか否かを判別する(ステップ4255)。一実
施例において、可能な新しい設定値は、0.02Vだけ
更に10%だけ電流設定値を越える必要がある。
【0310】可能な新しい設定値が十分に大きい場合、
コントローラは、それが最大許容設定値よりも大きいか
否かを判別する(ステップ4260)。そうであれば、
コントローラは、設定値を最大許容設定値に等しく設定
する(ステップ4265)。そうでなければ、コントロ
ーラは、設定値を可能な新しい設定値に等しく設定する
(ステップ4270)。次に、コントローラは、休止周
期カウントを適切な値と等しく設定することによって、
プロシージャの休止を実行する(ステップ4235)。
【0311】図116乃至図118は、設定値30の調
整プロシージャの動作を示す。図118に示すように、
システムは、0.3Vの設定値で開始する。8秒の動作
で、空気混和によって、導管内の材料の見かけの密度が
低下する(図116)。空気混和を伴って減衰が増加す
ると、駆動電流が増加し(図117)、更に、センサ電
圧のノイズも増加する(図118)。この時、何の変化
もない。何となれば、計器は、所望の設定値を維持する
ことが可能であるからである。
【0312】15秒の動作で、空気混和は増加し、更に
見かけの密度は減少する(図116)。このレベルの空
気混和で、駆動回路電流(図117)は、0.3Vの設
定値を維持するには不充分な最大値に達する。従って、
センサ電圧は、0.26V、すなわち、最大駆動回路電
流が維持できる電圧レベルに低下する(図118)。こ
の状態に応答して、コントローラは、設定値を(約28
秒の動作で)、最大駆動回路電流の生成を必要としない
レベル(0.23V)に調整する。
【0313】約38秒の動作で、空気混和のレベルは減
少し、見かけの密度は増加する(図116)。これによ
って、駆動電流は減少する(図117)。40秒の動作
で、コントローラは、設定値を増やすことによって、こ
の状態に応答する(図118)。空気混和のレベルは減
少し、見かけの密度は、再び約48秒の動作で増加す
る。そして、コントローラは、設定値を0.3Vに増や
すことによって応答する。 K.性能の結果 ディジタル流量計は、従来のアナログ流量計に対して著
しい性能改良を示した。1つの実験において、1ロット
の材料を正確に測定する2つのタイプの計器の能力が調
べられた。いずれの場合も、ロットは、適切な流量計を
経由してタンクに供給され、タンクにおいて、ロットは
計量された。1200および2400ポンド(544.
8及び1089.6kg)のロットに対して、アナログ
計器は、200ポンド(90.8kg)の再現性で50
0ポンド(227kg)の平均偏位を呈した。対照的
に、デジタル計器は、2ポンド(0.908kg)の再
現性で、40ポンド(18.16kg)の平均偏位を呈
した。これは、明らかに相当な改善である。
【0314】いずれの場合も、導管および周囲の配管
は、ロットの開始時には空だった。これは、導管を充填
してバッチを開始することが実用的ではない多数のバッ
チの用途においては、重要である。バッチは、流管が充
填されて、終了した。正のオフセットが予測されるもの
がある。何となれば、流量計は、計量タンクの充填がス
タートする前に、管を満たすために必要な材料を測定し
ているからである。動き出す際の遅延、または空気混和
している流れや低振動振幅によって生じる偏位は、負の
偏位を導入する傾向がある。真のバッチ処理に対して、
最も重要な問題は、測定値の再現性である。
【0315】結果は、アナログ流量計を使用すると、大
きな負の偏位と200ポンドだけの再現性とが存在する
ことを示している。これは、流れの発生後に開始に持っ
ていく時間の長さと、完全な振動振幅が得られるまでの
測定誤差とに起因する。比較すると、ディジタル流量計
は、正の偏位を得、これは、空のパイプを充填すること
に起因し、再現性は2ポンドである。
【0316】他の実験は、2つのタイプの計器の一般的
な測定精度を比較した。図119は、精度と、計器の最
大推奨流量の別々のパーセンテージで2つのタイプの計
器によって生成された測定値の対応する不確実性とを示
す。高流量(すなわち、最大流量の少なくとも25%)
で、アナログ計器は、デジタル計器の0.005%以下
と比較すると、0.15%以下で実際の値に相当する測
定値を生成する。低流量で、アナログ計器の偏位は、デ
ジタル計器の0.25%と比較すると、1.5%のオー
ダである。 L.自己確認性を備えた計器 ディジタル流量計は、自己確認センサを含む制御システ
ムにおいて使用される。このために、ディジタル流量計
は、自己確認計器として実行される。自己確認計器およ
び他のセンサは、「自己確認センサ」と題された米国特
許第5,570,300号に記載されている。
【0317】一般に、自己確認計器は、計器に対して有
効な全ての情報に基づいて、モニタすべきパラメータ
(例えば質量流)の値の最良推定値を提供する。最良推
定値は、一部が非測定データに基づいているので、最良
の推定値は、必ずしも現在の不完全な測定値データによ
って示される値とは一致しない。自己確認計器も、セン
サの動作状況に関する情報と同じように、最良推定値の
信頼性及び不確実性に関する情報を提供する。不確実性
情報は、周知の不確実性分析から導出され、不良が無い
場合においても提供される。
【0318】一般に、自己確認計器は、4つの基礎的な
パラメータ、すなわち、確認された測定値(VMV)、
確認された不確実性(VU)、測定値が生成された状況
の指示(MV状況)、装置状況を提供する。VMVは、
測定されたパラメータの値の計器の最良推定値である。
VUおよびMV状況は、VMVと関連する。計器は、各
測定値に対して、別々のVMV、VU、MV状況を生成
する。装置状況は、計器の動作状況を示す。
【0319】計器は、他の情報も提供する。例えば、制
御システムからの要請で、計器は、計器の状況に関する
詳細な診断情報を提供する。また、測定値が所定限界値
を越えたとき、または越えようとするとき、計器は、ア
ラーム信号を制御システムに送る。異なるアラームレベ
ルは、測定値が所定の値から偏位している程度を示すた
めに使用される。
【0320】VMVおよびVUは、数値である。例え
ば、VMVは、200度且つVUで評価される温度測定
値であり、VMVの不確実性は9度である。この場合、
測定すべき実際温度がVMV近傍の包絡線内に落ちて、
VU(すなわち、191度から209度まで)によって
示される高い可能性(大抵95%)が存在する。コント
ローラは、センサからの基本的なデータに基づいて、V
MVを生成する。第一に、コントローラは、センサから
の信号に基づく生の測定値値(RMV)を導く。一般
に、コントローラが異常を検出しないとき、コントロー
ラは、RMVに公称の信頼性を有して、VMVをRMV
と等しく設定する。コントローラがセンサに異常を検出
するとき、コントローラは、VMVをRMVと等しく設
定しない。その代わりに、コントローラは、コントロー
ラが実際のパラメータのRMVよりは優れた推定値であ
ると考えている値と等しく、VMVを設定する。
【0321】コントローラは、RMVの力学的不確実性
分析の結果である生の不確実性信号(RU)に基づい
て、VUを生成する。コントローラは、各サンプリング
期間中、この不確実性分析を実行する。不確実性分析
は、元々S.J.クライン(S.J.Kline)及び
F.A.マックリントック(F.A.McClinto
ck)により「単一のサンプル実験での不確実性の説明
(Single Sample Experiment
s)」(Mech.Eng.、75、3−8、1953
年)に記載され、広範囲に応用されて、較正の国際基準
の状況を達成した。本質的に、不確実性分析は、測定値
の「質」を表示する。あらゆる測定値が、当然夫知の付
随する誤差を有する。しかし、その誤差の合理的な限度
は、一つの不確実性番号によって表される(ANSI/
ASME PTC 19.1−1985、パート1、測
定値の不確実性(Measurement Uncer
tainty)、機器及び装置(Instrument
s and Apparatus))。
【0322】クライン及びマックリントックによって記
載されるように、観察された測定値Mに対して、M、w
の不確実性は、次式で定義される。
【0323】
【数110】 但し、Mは、あるレベル(大抵は95%)の信頼で真実
(Mtrue)である。この不確実性は、測定値の比
(すなわち、w/M)として相対的に簡単に表現され
る。
【0324】一般に、VUは、理想的な条件下(すなわ
ち、故障のないセンサが制御された研究室環境で作動し
ている)においても、ゼロ以外の値である。これは、セ
ンサにより生成される測定値が完全に確実でなくて、常
に誤差に対する可能性が存在するからである。VMVと
同様に、コントローラが異常を検出しないとき、コント
ローラは、VUをRUと等しく設定する。コントローラ
が、部分的にのみRMVの信頼性に影響を及ぼす欠点を
検出するとき、コントローラは、大抵、欠点の影響を考
慮する新しい不確実性分析を実行し、VUをこの分析の
結果に等しく設定する。コントローラが、RMVは実際
の測定された値に対して何の関係も持たないことを判別
するとき、コントローラは、過去の性能に基づいた値に
VUを値に設定する。
【0325】制御システムがVMVおよびVUを適切に
使用することを確実にするために、MV状況は、それら
の計算方法に関する情報を提供する。コントローラは、
全ての条件下において(センサが停止している時でさ
え)、VMVおよびVUを生成する。制御システムは、
VMVおよびVUが「生」データ、または履歴データの
どちらに基づいているかを知る必要がある。例えば、制
御システムがフィードバック制御においてVMV及びV
Uを使用し、更にセンサが停止している場合、制御シス
テムは、VMVおよびVUが過去の性能に基づいている
ことを知らなければならない。
【0326】MV状況が、異常状態の予想された持続
と、RMVのコントローラの信頼性とに基づいている。
MV状況に対する4つの主要な状態を、表1に従って生
成する。
【0327】
【表5】 クリアMV状況は、RMVがあるプロセス条件に対して
通常の範囲内にあるときに起こる。まぶしいMV状況
は、RMVはかなり異常であるが、異常は短期間である
と予測されることを示す。多くの場合、1つのセンサか
らの信号に突然の変化が存在し、コントローラが、この
変化が、未だ診断されていないセンサの故障、または測
定されるべき変数の急激な変化のいずれによるものかを
明らかにできないとき、コントローラは、MV状況を
「まぶしい」に設定する。「かすむ」MV状況は、RM
Vは異常であるが、測定されるべきパラメータと合理的
に関係していることを示す。例えば、RMVが雑音が多
い信号であるとき、コントローラは、MV状況を「かす
む」に設定する。「盲目」MV状況は、RMVに完全に
信頼性が無く、欠点の持続が予測されることを示す。
【0328】MV状況に対する2つのさらなる状態は、
「未確認」および「安全」である。コントローラがVM
Vを確認しないときに、MV状況は「未確認」である。
コントローラが公称の信頼性を有する冗長な測定値か
ら、VMVが生成されるとき、MV状況は「安全」であ
る。装置状況は、計器の健康を要約している、一般的且
つ離散的な値である。これは、主に、制御システムの欠
点検出および保守ルーチンによって使用される。多くの
場合、装置状況32は、6つの状態のうちの1つにあ
り、6つの状態の各々は、計器毎の異なる動作状況を示
す。これらの状態は、良好(GOOD)、試験(TBS
TING)、疑い(SUSPBCT),損傷(IMPA
IRED),不良(BAD),重要(CRITICA
L)である。GOOD装置状況は、計器が規定状態にあ
ることを意味する。TESTING装置状況は、計器が
自己チェックを実行していることを意味し、この自己チ
ェックは、測定値の質の一時的な低下に対しても責任が
ある。SUSPECT装置状況は、計器は異常な反応を
生成したが、コントローラは詳細な欠陥診断法を持たな
いことを意味する。IMPAIRED装置状況は、計器
が、性能に負の影響を与える診断された欠点に苦しんで
いることを意味する。BAD装置状況は、計器がかなり
誤動作し、保守が必要であることを意味する。最後に、
CRITICAL装置状況は、計器がリークや火事、爆
発等の危険を引き起こしかねない程度まで誤動作するこ
とを意味する。
【0329】図120は、自己確認計器のコントローラ
が、デジタル化されたセンサ信号を処理して、付随する
不確実性および測定値状況によって確認された質量流測
定値と、駆動信号とを生成するプロシージャ4500を
示す。最初に、コントローラは、センサからのデータを
集める(ステップ4505)。このデータを使用して、
コントローラは、センサ信号の周波数を測定する(ステ
ップ4510)。周波数が予測範囲に入る場合(ステッ
プ4515)、コントローラは、センサ信号からゼロオ
フセットを除去して(ステップ4520)、センサ信号
の振幅を測定し(ステップ4525)、位相を測定する
(ステップ4530)。コントローラは、これらの計算
値を使用して駆動信号を生成し(ステップ4535)、
生の質量流測定値および他の測定値を生成する(ステッ
プ4540)。
【0330】周波数が予測範囲に入らない場合(ステッ
プ4515)、コントローラは、失速プロシージャを実
行して(ステップ4545)、導管が失速したか否かを
判別し、それに応じて反応する。失速プロシージャにお
いて、コントローラは、駆動回路の利得を最大にして、
ゼロクロスに対してより広範囲の検索を実行し、導管が
全体で振動しているか否かを判別する。
【0331】導管が正しく振動していない場合(すなわ
ち、導管が振動していない場合、または導管が許容しが
たい高周波数で(例えば、共振周波数の高調波で)振動
している場合)(ステップ4550)、コントローラ
は、例えば駆動回路で方形波を注射することによって、
導管の通常の振動を再スタートせしめようとする(ステ
ップ4555)。振動の再スタートを試みた後に、コン
トローラは、MV状況を「まぶしい」に設定して(ステ
ップ4560)、無意味な生の測定値を生成する(ステ
ップ4565)。導管が正しく振動している場合(ステ
ップ4550)、コントローラは、ゼロオフセットを除
去して(ステップ4520)、上記のごとく進む。
【0332】生の測定値を生成した後(ステップ454
0、4565)、コントローラは、診断を実行し(ステ
ップ4570)、計器が正しく作動しているか否かを判
別する(ステップ4575)。なお、コントローラは、
各周期の間必ずしもこれらの診断を実行しているわけで
はない。次に、コントローラは、不確実性分析を実行し
(ステップ4580)、生の不確実性値を生成する。生
の測定値、診断の結果、他の情報を使用して、コントロ
ーラは、VMV、VU、MV状況と、装置状況とを生成
する(ステップ4585)。その後、コントローラは、
新しいデータセットを収集し、プロシージャを繰り返
す。プロシージャ4500の行程は、順次、または並行
に実行してもよく、更に、順番を変えて実行してもよ
い。
【0333】他の実施例において、空気混和が検出され
る時、質量流の補正は、上記のように適用され、MV状
況は「かすむ」になり、不確実性は、補正法の考えられ
る誤差を反映して増加する。例えば、50%の流速で作
動する流管に対して、通常の動作条件下では、不確実性
は、流速の0.1〜0.2%のオーダである。空気混和
が生じて、上記技術を使用して補正される場合、不確実
性は読取値のおよそ2%に増加する。不確実性値は減少
すべきである。何となれば、空気混和の効果の知識が良
くなり、空気混和を補償する能力が良くなったからであ
る。流量の不確実性が可変であるバッチの状況において
(例えば空に対して一回分で処理する開始や終了時には
高く、または空気混和や空洞の一時的な装入時に)、全
バッチの不確実性は、名目は低い不確実性を有する残り
のバッチに対して高い不確実性の期間の重み付けされた
重要性を反映する。これは、会計や他の測定用途に非常
に役に立つ質の測定基準である。 M.2配線流量計 他の実施例も、考慮する。例えば、図121に示すよう
に、上記方法は、1対の配線4605で双方向通信を実
行する「2配線(two−wire)」コリオリ流量計
4600を操作するために使用される。電力回路461
0は、デジタルコントローラ4615を操作するととも
に、駆動回路4620に給電して導管4625を振動さ
せる電力を受け取る。例えば、電力回路は、コントロー
ラに動作電力を提供する定出力回路4630と、過剰な
電力を使用して充電される駆動コンデンサ4635とを
含む。電力回路は、配線4605から、または第2の配
線対から電力を受け取る。デジタルコントローラは、一
つ以上のセンサ4640から信号を受け取る。
【0334】駆動コンデンサは適切に充電されるとき、
コントローラ4615は、コンデンサ4635を放電せ
しめて導管4625を駆動する。例えば、コントローラ
は、10周期ごとに1回、導管を駆動する。コントロー
ラ4615は、センサ4640から信号を受け取って分
析して質量流測定値を生成し、次に、コントローラは、
配線4605に伝える。
【0335】他の実施例は、記載した特許請求の範囲に
含まれる。
【図面の簡単な説明】
【図1】デジタル質量流量計のブロック図である。
【図2】質量流量計の機械的部品の斜視図である。
【図3】質量流量計の機械的部品の側面図である。
【図4】図1の流量計の動作モードの1つを示す略図で
ある。
【図5】図1の流量計の動作モードの1つを示す略図で
ある。
【図6】図1の流量計の動作モードの1つを示す略図で
ある。
【図7】アナログ制御および測定回路のブロック図であ
る。
【図8】デジタル質量流量計のブロック図である。
【図9】図8の計器の動作を示すフローチャートであ
る。
【図10】センサデータのグラフである。
【図11】センサデータのグラフである。
【図12】時間に対するセンサ電圧の変化を示すグラフ
である。
【図13】時間に対するセンサ電圧の変化を示すグラフ
である。
【図14】曲線はめ込みプロシージャのフローチャート
である。
【図15】位相差を生成するプロシージャのフローチャ
ートである。
【図16】システム始動時の駆動信号を示す。
【図17】システム始動時の駆動信号を示す。
【図18】システム始動時のセンサ電圧を示す。
【図19】システム始動時のセンサ電圧を示す。
【図20】システム始動時の駆動信号を示す。
【図21】システム始動時の駆動信号を示す。
【図22】システム始動時のセンサ電圧を示す。
【図23】システム始動時のセンサ電圧を示す。
【図24】システム始動時の駆動信号を示す。
【図25】システム始動時の駆動信号を示す。
【図26】システム始動時のセンサ電圧を示す。
【図27】システム始動時のセンサ電圧を示す。
【図28】同期変調技術を使用して、センサデータの周
波数、振幅、位相を測定するプロシージャのフローチャ
ートである。
【図29】質量流量計のブロック図である。
【図30】質量流量計のブロック図である。
【図31】図29および図30の計器により実行される
プロシージャのフローチャートである。
【図32】伝達関数のログ・振幅制御を示す。
【図33】根軌跡プロットを示す図である。
【図34】温度に対するアナログ−ディジタルコンバー
タ性能のグラフである。
【図35】温度に対するアナログ−ディジタルコンバー
タ性能のグラフである。
【図36】温度に対するアナログ−ディジタルコンバー
タ性能のグラフである。
【図37】温度に対するアナログ−ディジタルコンバー
タ性能のグラフである。
【図38】位相測定値のグラフである。
【図39】位相測定値のグラフである。
【図40】位相測定値のグラフである。
【図41】位相測定値のグラフである。
【図42】位相測定値のグラフである。
【図43】ゼロオフセットの補償プロシージャのフロー
チャートである。
【図44】位相測定値のグラフである。
【図45】位相測定値のグラフである。
【図46】位相測定値のグラフである。
【図47】位相測定値のグラフである。
【図48】位相測定値のグラフである。
【図49】センサ電圧のグラフである。
【図50】ダイナミックな効果を補償するためのプロシ
ージャのフローチャートである。
【図51】図29のプロシージャの適用を示すグラフで
ある。
【図52】図29のプロシージャの適用を示すグラフで
ある。
【図53】図29のプロシージャの適用を示すグラフで
ある。
【図54】図29のプロシージャの適用を示すグラフで
ある。
【図55】図29のプロシージャの適用を示すグラフで
ある。
【図56】図29のプロシージャの適用を示すグラフで
ある。
【図57】図29のプロシージャの適用を示すグラフで
ある。
【図58】図29のプロシージャの適用を示すグラフで
ある。
【図59】図29のプロシージャの適用を示すグラフで
ある。
【図60】図29のプロシージャの適用を示すグラフで
ある。
【図61】図29のプロシージャの適用を示すグラフで
ある。
【図62】図29のプロシージャの適用を示すグラフで
ある。
【図63】図29のプロシージャの適用を示すグラフで
ある。
【図64】図29のプロシージャの適用を示すグラフで
ある。
【図65】図29のプロシージャの適用を示すグラフで
ある。
【図66】図29のプロシージャの適用を示すグラフで
ある。
【図67】図29のプロシージャの適用を示すグラフで
ある。
【図68】図29のプロシージャの適用を示すグラフで
ある。
【図69】図29のプロシージャの適用を示すグラフで
ある。
【図70】図29のプロシージャの適用を示すグラフで
ある。
【図71】図29のプロシージャの適用を示すグラフで
ある。
【図72】図29のプロシージャの適用を示すグラフで
ある。
【図73】位相測定値を示しているグラフである。
【図74】位相測定値を示しているグラフである。
【図75】位相測定値を示しているグラフである。
【図76】位相測定値を示しているグラフである。
【図77】位相測定値を示しているグラフである。
【図78】位相測定値を示しているグラフである。
【図79】位相測定値を示しているグラフである。
【図80】位相測定値を示しているグラフである。
【図81】位相測定値を示しているグラフである。
【図82】位相測定値を示しているグラフである。
【図83】位相測定値を示しているグラフである。
【図84】位相測定値を示しているグラフである。
【図85】センサ電圧のグラフである。
【図86】図85のセンサ電圧に相当する位相のグラフ
である。
【図87】図85のセンサ電圧に相当する周波数の測定
値である。
【図88】図86および図87の位相および周波数の測
定値に対する補正パラメータのグラフである。
【図89】図86および図87の位相および周波数の測
定値に対する補正パラメータのグラフである。
【図90】生の測定値のグラフである。
【図91】生の測定値のグラフである。
【図92】生の測定値のグラフである。
【図93】生の測定値のグラフである。
【図94】生の測定値のグラフである。
【図95】生の測定値のグラフである。
【図96】生の測定値のグラフである。
【図97】生の測定値のグラフである。
【図98】補正した測定値のグラフである。
【図99】補正した測定値のグラフである。
【図100】補正した測定値のグラフである。
【図101】補正した測定値のグラフである。
【図102】補正した測定値のグラフである。
【図103】補正した測定値のグラフである。
【図104】補正した測定値のグラフである。
【図105】補正した測定値のグラフである。
【図106】空気混和に対する補正を示すグラフであ
る。
【図107】空気混和に対する補正を示すグラフであ
る。
【図108】空気混和に対する補正を示すグラフであ
る。
【図109】空気混和に対する補正を示すグラフであ
る。
【図110】空気混和に対する補正を示すグラフであ
る。
【図111】空気混和に対する補正を示すグラフであ
る。
【図112】空気混和に対する補正を示すグラフであ
る。
【図113】空気混和に対する補正を示すグラフであ
る。
【図114】導管における空気混和の影響を示すブロッ
ク図である。
【図115】設定値の制御プロシージャのフローチャー
トである。
【図116】図114のプロシージャの適用を示すグラ
フである。
【図117】図114のプロシージャの適用を示すグラ
フである。
【図118】図114のプロシージャの適用を示すグラ
フである。
【図119】デジタルおよびアナログ流量計の性能を比
較するグラフである。
【図120】自己確認性を備えた計器の動作を示すフロ
ーチャートである。
【図121】2ワイヤ・デジタル質量流量計のブロック
図である。
【符号の説明】
100 ディジタル流量計 105 制御・測定システム 110、125 センサ 115 駆動回路 120 導管 100 ディジタル流量計 105 制御・測定システム 110、125 センサ 115 駆動回路 120 導管
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成11年3月18日(1999.3.1
8)
【手続補正1】
【補正対象書類名】図面
【補正対象項目名】全図
【補正方法】変更
【補正内容】
【図1】
【図3】
【図2】
【図4】
【図5】
【図6】
【図7】
【図10】
【図11】
【図8】
【図9】
【図12】
【図13】
【図38】
【図14】
【図15】
【図16】
【図28】
【図39】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図40】
【図44】
【図45】
【図46】
【図23】
【図24】
【図25】
【図26】
【図27】
【図32】
【図29】
【図30】
【図31】
【図33】
【図34】
【図35】
【図36】
【図37】
【図41】
【図50】
【図63】
【図42】
【図47】
【図48】
【図51】
【図43】
【図49】
【図52】
【図53】
【図85】
【図54】
【図55】
【図56】
【図57】
【図58】
【図59】
【図60】
【図61】
【図62】
【図64】
【図65】
【図66】
【図67】
【図68】
【図69】
【図70】
【図86】
【図87】
【図71】
【図72】
【図73】
【図74】
【図75】
【図88】
【図89】
【図76】
【図77】
【図90】
【図91】
【図78】
【図79】
【図92】
【図93】
【図80】
【図81】
【図94】
【図95】
【図82】
【図83】
【図96】
【図97】
【図84】
【図98】
【図99】
【図100】
【図101】
【図102】
【図103】
【図104】
【図105】
【図106】
【図114】
【図116】
【図117】
【図118】
【図107】
【図108】
【図109】
【図110】
【図111】
【図112】
【図113】
【図119】
【図115】
【図120】
【図121】
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ディビッド ダブリュ. クラーク イギリス国 オックスフォード 0エック ス3 8エスエックス ヘディントン オ ールドロード 98 (72)発明者 ジェイムズ エイチ. ヴィグノス アメリカ合衆国 マサチューセッツ州 02194 ニードハムハイツ マニングスト リート 129

Claims (105)

    【特許請求の範囲】
  1. 【請求項1】 振動可能な導管と、 導管に接続されて導管に運動を与えるように操作可能な
    駆動回路と、 導管に接続されて導管の運動を検出するように操作可能
    なセンサと、 駆動回路およびセンサの間に接続される制御・測定シス
    テムとからなり、 前記制御・測定システムは、 センサからセンサ信号を受け取り、 ディジタル信号処理を使用してセンサ信号に基づいた駆
    動信号を生成し、 駆動信号を駆動回路への出力し、 センサからの信号に基づいて導管を流れている材料の特
    性の測定値を生成することを特徴とするディジタル流量
    計。
  2. 【請求項2】 更に、導管に接続されて導管の運動を検
    出するように操作可能な第2のセンサを有し、 前記制御・測定システムは、前記第2のセンサに接続さ
    れており、 第2のセンサからの第2のセンサ信号を受け取り、 ディジタル信号処理を使用して第1及び第2のセンサ信
    号に基づいて駆動信号を生成し、 第1および第2のセンサ信号に基づいて導管を流れてい
    る材料の特性の測定値を生成することを特徴とする請求
    項1記載のディジタル流量計。
  3. 【請求項3】 前記制御・測定システムは、2つの駆動
    回路に対して異なる駆動信号を生成することを特徴とす
    る請求項2記載のディジタル流量計。
  4. 【請求項4】 前記制御・測定システムは、異なる周波
    数を有する駆動信号を生成することを特徴とする請求項
    3記載のディジタル流量計。
  5. 【請求項5】 前記制御・測定システムは、異なる振幅
    を有する駆動信号を生成することを特徴とする請求項4
    記載のディジタル流量計。
  6. 【請求項6】 更に、駆動回路に供給される電流を測定
    する回路を有することを特徴とする請求項1記載のディ
    ジタル流量計。
  7. 【請求項7】 前記回路は、 駆動回路と直列に接続された抵抗器と、 前記抵抗器と並列に接続されたアナログ−ディジタルコ
    ンバータと、を有し、 前記アナログ−ディジタルコンバータは、 前記抵抗器の電圧を測定し、 測定された電圧をディジタル値に変換し、 前記ディジタル値を前記制御・測定値システムに供給す
    ることを特徴とする請求項6記載のディジタル流量計。
  8. 【請求項8】 前記導管の入口で第1の圧力を測定する
    ために接続された第1の圧力センサと、 前記導管の出口で第2の圧力を測定するために接続され
    た第2の圧力センサとを更に有することを特徴とする請
    求項1記載のディジタル流量計。
  9. 【請求項9】 アナログ−ディジタルコンバータが更に
    接続され、前記アナログ−ディジタルコンバータは、第
    1の圧力センサおよび第2の圧力センサによって生成さ
    れる信号をディジタル値に変換し、前記ディジタル値を
    前記制御・測定システムに供給することを特徴とする請
    求項8記載のディジタル流量計。
  10. 【請求項10】 前記導管への入り口の第1の温度を測
    定するために接続された第4の温度センサと、 前記導管の出口での第2の温度を測定するために接続さ
    れた第2の温度センサとを更に有することを特徴とする
    請求項8記載のディジタル流量計。
  11. 【請求項11】 前記制御・測定システムは、 第1のセンサ信号の周波数を推定する行程と、 第1のセンサ信号を使用して位相差を計算する行程と、 計算された位相差を使用して測定値を生成する行程と、
    によって、前記特性の測定値を生成することを特徴とす
    る請求項2記載のディジタル流量計。
  12. 【請求項12】 前記制御・測定システムは、センサ信
    号の一つの振幅を調節することによって、センサ信号の
    振幅差を補償することを特徴とする請求項11記載のデ
    ィジタル流量計。
  13. 【請求項13】 前記制御・測定システムは、一方のセ
    ンサ信号の振幅に、一方のセンサ信号の振幅に対する他
    方のセンサ信号の振幅の比を乗算することによって、セ
    ンサ信号の振幅差を補償することを特徴とする請求項1
    2記載のディジタル流量計。
  14. 【請求項14】 前記センサ信号は周期的であり、前記
    制御・測定システムは、前記センサ信号をセット単位で
    処理し、前記セットの各々は、周期的センサ信号の全周
    期に対するデータを含むことを特徴とする請求項1記載
    のディジタル流量計。
  15. 【請求項15】 連続するセットは、周期的センサ信号
    の周期を重畳するためのデータを含むことを特徴とする
    請求項14記載のディジタル流量計。
  16. 【請求項16】 前記制御・測定システムは、前の周期
    の周波数を使用して、周期の終点を評価することを特徴
    とする請求項14記載のディジタル流量計。
  17. 【請求項17】 前記制御・測定システムは、1周期に
    対するデータを分析して、前記周期が更なる処理に値す
    るか否かを判別することを特徴とする請求項14記載の
    ディジタル流量計。
  18. 【請求項18】 前記制御・測定システムは、1の周期
    に対するデータが、前記データの予測動作、すなわち、
    前の周期の1つ以上のパラメータに基づいた予測動作に
    一致しない時、前記周期はさらなる処理に値しないと判
    定することを特徴とする請求項17記載のディジタル流
    量計。
  19. 【請求項19】 前記制御・測定システムは、1の周期
    に対する周波数が前の周期に対する周波数とは閾値量以
    上異なるとき、前記周期はさらなる処理に値しないこと
    を判別することを特徴とする請求項18記載のディジタ
    ル流量計。
  20. 【請求項20】 前記制御・測定システムは、前記周期
    に対する周波数が前の周期に対する周波数と等しい場合
    に生じ得る値に、前記周期でのある点の値を比較するこ
    とによって、前記周期に対する周波数は前の周期に対す
    る周波数とは異なるか否かを判別することを特徴とする
    請求項19記載のディジタル流量計。
  21. 【請求項21】 前記制御・測定システムは、センサ信
    号のゼロクロスを検出し、ゼロクロス間のサンプルを数
    えることによって、センサ信号の周波数を判別すること
    を特徴とする請求項1記載のディジタル流量計。
  22. 【請求項22】 前記制御・測定システムは、反復曲線
    あてはめ法を使用して、センサ信号の周波数を測定する
    ことを特徴とする請求項1記載のディジタル流量計。
  23. 【請求項23】 前記制御・測定システムは、フーリエ
    解析を使用してセンサ信号の振幅を測定し、駆動信号を
    生成する時に、測定された振幅を使用することを特徴と
    する請求項1記載のディジタル流量計。
  24. 【請求項24】 前記制御・測定システムは、センサ信
    号毎のオフセットを測定し、位相偏差を比較して位相差
    を測定することを特徴とする請求項11記載のディジタ
    ル流量計。
  25. 【請求項25】 前記制御・測定システムは、フーリエ
    解析を使用して位相差を測定することを特徴とする請求
    項11記載のディジタル流量計。
  26. 【請求項26】 前記制御・測定システムは、各センサ
    信号に対する周波数、振幅、位相偏差を測定し、センサ
    信号の周波数の平均に対して前記位相偏差を拡大・縮小
    することを特徴とする請求項11記載のディジタル流量
    計。
  27. 【請求項27】 前記制御・測定システムは、複数の方
    法を使用して前記位相差を計算し、前記方法の一つの結
    果を計算された位相差として選択することを特徴とする
    請求項11記載のディジタル流量計。
  28. 【請求項28】 前記制御・測定システムは、前記セン
    サ信号を合成して合成信号を生成し、前記合成信号に基
    づいて駆動信号を生成することを特徴とする請求項2記
    載のディジタル流量計。
  29. 【請求項29】 前記制御・測定システムは、利得を前
    記合成信号に適用することによって、駆動信号を生成す
    ることを特徴とする請求項28記載のディジタル流量
    計。
  30. 【請求項30】 前記制御・測定システムは、大なる利
    得を前記合成信号に適用して導管の運動を開始し、前記
    運動が開始された後、センサ信号の位相および位相に基
    づいた位相及び周波数を有する周期的信号を駆動信号と
    して生成することによって、前記駆動信号を生成するこ
    とを特徴とする請求項29記載のディジタル流量計。
  31. 【請求項31】 前記制御・測定システムは、前記セン
    サ信号を合計することによって、前記センサ信号を合成
    することを特徴とする請求項28記載のディジタル流量
    計。
  32. 【請求項32】 前記制御・測定システムは、前記セン
    サ信号を合計する前に、一方のセンサ信号の振幅に対す
    る他方のセンサ信号の振幅の比を使用して、一方のセン
    サ信号の振幅を訂正することを特徴とする請求項31記
    載のディジタル流量計。
  33. 【請求項33】 前記制御・測定システムは、 駆動信号を生成する第1のモードの信号生成を使用して
    前記導管の運動を開始し、 駆動信号を生成する第2のモードの信号生成を使用して
    前記導管の運動を維持することを特徴とする請求項1記
    載のディジタル流量計。
  34. 【請求項34】 第1のモードの信号生成は所望の特性
    を有する周期的信号の合成から成り、第2のモードの信
    号生成は、センサ信号を含むフィードバックループを使
    用することから成ることを特徴とする請求項33記載の
    ディジタル流量計。
  35. 【請求項35】 前記所望の特性は、導管振動の所望の
    初期周波数であることを特徴とする請求項34記載のデ
    ィジタル流量計。
  36. 【請求項36】 第1のモードの信号生成は、前記セン
    サ信号を含むフィードバックループを使用することから
    成り、第2のモードの信号生成は、所望の特性を有する
    周期的信号の合成から成ることを特徴とする請求項33
    記載のディジタル流量計。
  37. 【請求項37】 所望の特性の周波数および位相は、前
    記センサ信号の周波数および位相に相当することを特徴
    とする請求項36記載のディジタル流量計。
  38. 【請求項38】 前記制御・測定システムは、適応可能
    な周期的駆動信号を生成することを特徴とする請求項1
    記載のディジタル流量計。
  39. 【請求項39】 前記制御・測定システムおよび前記駆
    動回路との間に接続された正及び負直流電源を更に有
    し、 前記制御・測定システムは、センサ信号に基づいた位相
    および周波数を有する間隔で前記電流源のオン・オフを
    切り替えることによって、駆動信号を生成することを特
    徴とする請求項38記載のディジタル流量計。
  40. 【請求項40】 前記制御・測定システムは、センサ信
    号の特性に相当する特性を有するサイン波を合成するこ
    とによって、駆動信号を生成することを特徴とする請求
    項38記載のディジタル流量計。
  41. 【請求項41】 前記制御・測定システムは、センサ信
    号の位相および周波数に相当する位相及び周波数を有す
    るサイン波を合成することによって、駆動信号を生成す
    ることを特徴とする請求項40記載のディジタル流量
    計。
  42. 【請求項42】 前記制御・測定システムは、センサ信
    号の一つ以上の特性に基づいて駆動信号を生成するとき
    に使用する利得を、ディジタル方式で生成することを特
    徴とする請求項1記載のディジタル流量計。
  43. 【請求項43】 前記制御・測定システムは、前記セン
    サ信号の振幅に基づいて前記利得をディジタル方式で生
    成することを特徴とする請求項42記載のディジタル流
    量計。
  44. 【請求項44】 前記制御・測定システムは、PI制御
    アルゴリズムをディジタル方式で実行して、導管振動の
    振幅を調節することを特徴とする請求項42記載のディ
    ジタル流量計。
  45. 【請求項45】 前記駆動回路は、振動運動を前記導管
    に与えるように操作可能であることを特徴とする請求項
    1記載のディジタル流量計。
  46. 【請求項46】 前記制御・測定システムは、導管振動
    の振幅をユーザ制御値に維持するために、センサ符号に
    基づいて、駆動信号をディジタル方式で生成することを
    特徴とする請求項45記載のディジタル流量計。
  47. 【請求項47】 前記制御・測定システムは、 振動の振幅がユーザ制御値を越える時は、負の駆動信号
    を生成し、これによって駆動回路を導管の運動に抵抗せ
    しめ、 振動の振幅がユーザ制御値未満のときは、正の駆動信号
    を生成し、これによって駆動回路が導管に運動を与える
    ことを特徴とする請求項46記載のディジタル流量計。
  48. 【請求項48】 前記制御・測定システムは、 センサ信号に基づいて利得信号を生成するコントローラ
    と、 コントローラに接続されて、利得信号を受け取るととも
    に、前記利得信号に基づいで駆動信号を生成する乗算デ
    ィジタル−アナログ変換器と、を有することを特徴とす
    る請求項1記載のディジタル流量計。
  49. 【請求項49】 前記導管に接続されて前記導管の運動
    を検出するように操作可能な第2のセンサを更に有し、 前記制御・測定システムは、 測定値を生成するコントローラと、 前記第1のセンサおよび前記コントローラの間に接続さ
    れて前記コントローラに第1のディジタルセンサ信号を
    提供する第1のアナログ−ディジタルコンバータと、 前記第2のセンサおよび前記コントローラの間に接続さ
    れて前記コントローラに第2のディジタルセンサ信号を
    提供する第2のアナログーディジタルコンバータと、を
    有することを特徴とする請求項1記載のディジタル流量
    計。
  50. 【請求項50】 前記コントローラは、前記ディジタル
    センサ信号を合成して合成信号を生成し、前記第1およ
    び第2のディジタルセンサ信号に基づいて利得信号を生
    成し、 前記制御・測定システムは、乗算ディジタル−アナログ
    変換器が更に接続され、 前記乗算ディジタル−アナロ
    グ変換器は、コントローラから前記合成信号および利得
    信号を受け取って、前記合成信号および利得信号の積と
    して前記駆動信号を生成することを特徴とする請求項4
    9記載のディジタル流量計。
  51. 【請求項51】 前記制御・測定システムは、選択的に
    負の利得をセンサ信号に適用して導管の運動を減らすこ
    とを特徴とする請求項1記載のディジタル流量計。
  52. 【請求項52】 前記制御・測定システムは、センサ信
    号のゼロオフセットを補償することを特徴とする請求項
    1記載のディジタル流量計。
  53. 【請求項53】 前記ゼロオフセットは、利得変化に起
    因する成分と、利得の非線形性に起因する成分とを含
    み、 前記制御・測定システムは、別々に2つの成分を補償す
    ることを特徴とする請求項52記載のディジタル流量
    計。
  54. 【請求項54】 前記制御・測定システムは、一つ以上
    の補正因子を生成し、前記補正因子を使用してセンサ信
    号を修正することによって、ゼロオフセットを補償する
    ことを特徴とする請求項52記載のディジタル流量計。
  55. 【請求項55】 前記制御・測定システムは、 第1のセンサ信号に対する位相偏差を計算し、 第2のセンサ信号に対する位相偏差を計算することを特
    徴とする請求項2記載のディジタル流量計。
  56. 【請求項56】 位相偏差は、センサ信号のゼロクロス
    点と、センサ信号の基本的周波数に相当するセンサ信号
    の成分に対するゼロ位相の点との間の差として定義され
    ることを特徴とする請求項55記載のディジタル流量
    計。
  57. 【請求項57】 前記制御・測定システムは、計算され
    た位相オフセットを合成して位相差を生成することを特
    徴とする請求項56記載のディジタル流量計。
  58. 【請求項58】 前記制御・測定システムは、 第1のセンサ信号の周波数を推定し、 第1のセンサ信号の周波数とは異なる第2のセンサ信号
    の周波数を推定し、 推定された周波数を使用して、センサ信号間の位相差を
    計算することによって、前記特性の測定値を生成するこ
    とを特徴とする請求項2記載のディジタル流量計。
  59. 【請求項59】 前記センサは速度センサからなり、 前記制御・測定システムは、 センサ信号の周波数、振幅、位相を推定し、 推定した周波数、振幅、位相を補正して、速度センサと
    絶対位置センサとの間の性能差を説明することを特徴と
    する請求項1記載のディジタル流量計。
  60. 【請求項60】 前記制御・測定システムは、 センサ信号の第1のパラメータを推定し、 第2のパラメータの変化率を測定し、 測定された変化率に基づいて推定された第1のパラメー
    タを補正することを特徴とする請求項1記載のディジタ
    ル流量計。
  61. 【請求項61】 前記第4のパラメータは、前記センサ
    信号の周波数からなることを特徴とする請求項60記載
    のディジタル流量計。
  62. 【請求項62】 前記第2のパラメータは、前記導管の
    振動の周波数からなることを特徴とする請求項61記載
    のディジタル流量計。
  63. 【請求項63】 前記第2のパラメータは、前記導管の
    振動の振幅からなることを特徴とする請求項61記載の
    ディジタル流量計。
  64. 【請求項64】 前記第1のパラメータは、前記センサ
    信号の大きさを含むことを特徴とする請求項60記載の
    ディジタル流量計。
  65. 【請求項65】 前記第2のパラメータは、前記導管の
    振動の周波数からなることを特徴とする請求項63記載
    のディジタル流量計。
  66. 【請求項66】 前記第2のパラメータは、前記導管の
    振幅からなることを特徴とする請求項64記載のディジ
    タル流量計。
  67. 【請求項67】 前記第1のパラメータは、前記センサ
    信号の位相からなることを特徴とする請求項60記載の
    ディジタル流量計。
  68. 【請求項68】 前記第2のパラメータは、前記導管の
    振動の周波数からなることを特徴とする請求項67記載
    のディジタル流量計。
  69. 【請求項69】 前記第2のパラメータは、前記導管の
    振動の大きさからなることを特徴とする請求項67記載
    のディジタル流量計。
  70. 【請求項70】 前記制御・測定システムは、 第1のセンサ信号の第1のパラメータを推定し、 第1のセンサ信号の第2のパラメータの変化率を測定
    し、 第1のセンサ信号の第2のパラメータの測定された変化
    率に基づいて、第1のセンサ信号の推定された第1のパ
    ラメータを補正し、 第2のセンサ信号の第1のパラメータを推定し、 第2のセンサ信号の第2のパラメータの変化率を測定
    し、 第2のセンサ信号の第2のパラメータの測定された変化
    率に基づいて、第2のセンサ信号の推定された第1のパ
    ラメータを補正し、 前記制御・測定システムは、他のセンサ信号の推定、測
    定、補正とは独立に、各センサ信号の推定、測定、補正
    を実行することを特徴とする請求項2記載のディジタル
    流量計。
  71. 【請求項71】 前記第1のパラメータは、前記センサ
    信号の周波数からなることを特徴とする請求項70記載
    のディジタル流量計。
  72. 【請求項72】 前記第2のパラメータは、前記導管の
    振動の周波数からなることを特徴とする請求項71記載
    のディジタル流量計。
  73. 【請求項73】 前記第2のパラメータは、前記導管の
    振動の大きさからなることを特徴とする請求項71記載
    のディジタル流量計。
  74. 【請求項74】 前記第4のパラメータは、前記センサ
    信号の大きさからなることを特徴とする請求項70記載
    のディジタル流量計。
  75. 【請求項75】 前記第2のパラメータは、前記導管の
    振動の周波数からなることを特徴とする請求項73記載
    のディジタル流量計。
  76. 【請求項76】 前記第2のパラメータは、前記導管の
    振動の大きさからなることを特徴とする請求項74記載
    のディジタル流量計。
  77. 【請求項77】 前記第1のパラメータは、前記センサ
    信号の位相からなることを特徴とする請求項70記載の
    ディジタル流量計。
  78. 【請求項78】 前記第2のパラメータは、前記導管の
    振動の周波数からなることを特徴とする請求項77記載
    のディジタル流量計。
  79. 【請求項79】 前記第2のパラメータは、前記導管の
    振動の大きさからなることを特徴とする請求項77記載
    のディジタル流量計。
  80. 【請求項80】 前記ディジタル流量計は、質量流量計
    からなり、 導管を流れる材料の特性は、質量流量であることを特徴
    とする請求項1記載のディジタル流量計。
  81. 【請求項81】 前記制御・測定システムは、 初期の質量流量を測定する行程と、 前記導管を流れる材料の見かけの密度を測定する行程
    と、 前記見かけの密度を前記材料の周知の密度と比較して密
    度差を測定する行程と、 前記密度差に基づいて前記初期の質量流量を調節して、
    調節された質量流量を生成する行程と、によって、前記
    導管の空気混和の影響を説明することを特徴とする請求
    項80記載のディジタル流量計。
  82. 【請求項82】 前記制御および測定システムは、調節
    された質量流量を調節して減衰の効果を考慮することに
    よって、前記導管の空気混和の影響を説明することを特
    徴とする請求項81記載のディジタル流量計。
  83. 【請求項83】 前記導管に接続されて前記導管の運動
    を検出するように操作可能な第2のセンサを更に有し、 前記制御・測定システムは、前記第2のセンサに接続さ
    れて、 前記第2のセンサから第2のセンサ信号を受け取り、 第1および第2のセンサ信号の振幅の間の差に基づい
    て、調節された質量流量を調節することによって、前記
    導管の空気混和の影響を説明することを特徴とする請求
    項81記載のディジタル流量計。
  84. 【請求項84】 前記ディジタル流量計は、比重計から
    なり、 前記導管を流れる材料の特性は、材料の密度であること
    を特徴とする請求項1記載のディジタル流量計。
  85. 【請求項85】 前記センサ信号は、アナログ信号であ
    り、 前記制御・測定システムは、アナログ信号をディジタル
    信号に変換するアナログ−ディジタルコンバータから成
    ることを特徴とする請求項1記載のディジタル流量計。
  86. 【請求項86】 振動可能な前記導管は、2つの平らな
    ループから成ることを特徴とする請求項1記載のディジ
    タル流量計。
  87. 【請求項87】 前記センサおよび前記駆動回路は、前
    記ループの間に接続されていることを特徴とする請求項
    86記載のディジタル流量計。
  88. 【請求項88】 第2のセンサと、第2の駆動回路とを
    更に有し、 前記第1の駆動回路およびセンサは、前記ループの第1
    の端部の間に接続され、 前記第2の駆動回路およびセンサは、前記ループの前記
    第1の端部とは対向して位置する第2の端部の間に接続
    されていることを特徴とする請求項87記載のディジタ
    ル流量計。
  89. 【請求項89】 1対のワイヤにて電力を受け取るとと
    もに、電力を前記制御・測定システム及び駆動回路に供
    給する電力回路を更に有し、 前記制御・測定システムは、前記1対のワイヤで、前記
    導管を流れる材料の特性の測定値を伝達するように操作
    可能であることを特徴とする請求項1記載のディジタル
    流量計。
  90. 【請求項90】 前記電力回路は、電力を前記制御・測
    定システムに提供する定出力回路と、2つのワイヤから
    過剰な電力によって充電される駆動コンデンサと、から
    なることを特徴とする請求項89記載のディジタル流量
    計。
  91. 【請求項91】 前記制御・測定システムは、前記駆動
    コンデンサを放電して前記駆動回路に給電することを特
    徴とする請求項90記載のディジタル流量計。
  92. 【請求項92】 前記制御・測定システムは、駆動コン
    デンサの充電レベルをモニタし、前記コンデンサの充電
    レベルが閾値レベルに達した後前記駆動コンデンサを放
    電することを特徴とする請求項91記載のディジタル流
    量計。
  93. 【請求項93】 前記制御・測定システムは、周期的に
    駆動コンデンサを放電することを特徴とする請求項92
    記載のディジタル流量計。
  94. 【請求項94】 前記制御・測定システムは、前記1対
    のワイヤで双方向通信を実行することを特徴とする請求
    項89記載のディジタル流量計。
  95. 【請求項95】 前記センサ信号は周期的信号からな
    り、 前記制御・測定システムは、 周期的信号の期間の間に第1のデータセットを収集し、 第1のデータセットを処理して駆動信号および測定値を
    生成し、 第1のデータセットの処理と同時にセンサ信号の次の期
    間の間に第2のデータセットを収集することを特徴とす
    る請求項1記載のディジタル流量計。
  96. 【請求項96】 前記第1のデータセットに相当する期
    間は、第2のデータセットに相当する期間と重複するこ
    とを特徴とする請求項95記載のディジタル流量計。
  97. 【請求項97】 前記制御・測定システムは、 前記駆動信号を制御してセンサ信号の振幅を一定の設定
    値に維持し、 前記駆動信号が第1の閾値レベルを越える時、一定の設
    定値を減らし、 前記駆動信号が第2の閾値レベル未満であり、前記一定
    の設定値が設定値に対する最大許容値未満のときに、一
    定の設定値を増やすことを特徴とする請求項1記載のデ
    ィジタル流量計。
  98. 【請求項98】 前記第1の閾値レベルは、最大許容駆
    動信号の95%以下であることを特徴とする請求項97
    記載のディジタル流量計。
  99. 【請求項99】 前記制御・測定システムは、測定値に
    ついての不確実性分析を実行することを特徴とする請求
    項1記載のディジタル流量計。
  100. 【請求項100】 制御システムへの接続を更に有し、 前記制御・測定システムは、不確実性分析の結果及び測
    定値を前記制御システムに伝えることを特徴とする請求
    項99記載のディジタル流量計。
  101. 【請求項101】 前記制御・測定システムは、駆動信
    号の位相を調節するディジタル処理を使用して、前記セ
    ンサと前記駆動回路との間に接続された部品と前記セン
    サとに関係する時間遅延を補償することを特徴とする請
    求項1記載のディジタル流量計。
  102. 【請求項102】 振動可能な導管と、 前記導管に接続されて前記導管に運動を与えるように操
    作可能な駆動回路と、 前記導管に接続されて前記導管の運動を検出するととも
    に周期的センサ信号を生成するように操作可能なセンサ
    と、 前記駆動回路と前記センサとの間に接続される制御シス
    テムと、からなり、前記制御システムは、 センサ信号の現在の周期に相当するデータを収集し、 現在の周期の間にデータを処理して、駆動信号と、前記
    導管を流れる材料の特性の測定値とを生成し、 前記駆動回路へ前記駆動信号を出力し、 現在の周期に対するデータの処理と同時にセンサ信号の
    次の周期に相当するデータを収集することを特徴とする
    ディジタル流量計。
  103. 【請求項103】 導管を流れる材料の特性の測定値を
    生成する方法であって、 導管の運動を検出する行程と、 検出された運動に基づいて前記導管を流れる材料の特性
    の測定値を生成する行程と、 ディジタル信号処理を使用して、導管の検出された運動
    に基づいて導管に運動を与える駆動信号を生成する行程
    と、 駆動信号を使用して運動を前記導管に与える行程と、か
    らなることを特徴とする方法。
  104. 【請求項104】 前記導管の運動の検出は、2つの位
    置での運動の検出から成り、 前記測定値は、前記2つの位置で検出された運動に基づ
    いて生成され、 前記駆動信号は、前記2つの位置で検出される運動に基
    づいて生成されることを特徴とする請求項103記載の
    方法。
  105. 【請求項105】 ディジタル信号処理を使用して、前
    記導管の検出された運動に基づいて前記導管に運動を与
    える第2の駆動信号を生成する行程を更に有し、 前記導管へ運動を与える行程は、 前記第1の駆動信号を使用して第1の位置で運動を与え
    る行程と、 前記第2の駆動信号を使用して第2の位置で運動を与え
    る行程と、からなることを特徴とする請求項104記載
    の方法。
JP10378592A 1997-11-26 1998-11-26 ディジタル流量計 Withdrawn JP2000018995A (ja)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US6655497P 1997-11-26 1997-11-26
US09/111,739 US6311136B1 (en) 1997-11-26 1998-07-08 Digital flowmeter
US09/111,739 1998-07-08
US60/066,554 1998-07-08

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2007176131A Division JP2007248479A (ja) 1997-11-26 2007-07-04 ディジタル流量計

Publications (1)

Publication Number Publication Date
JP2000018995A true JP2000018995A (ja) 2000-01-21

Family

ID=26746870

Family Applications (2)

Application Number Title Priority Date Filing Date
JP10378592A Withdrawn JP2000018995A (ja) 1997-11-26 1998-11-26 ディジタル流量計
JP2007176131A Pending JP2007248479A (ja) 1997-11-26 2007-07-04 ディジタル流量計

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2007176131A Pending JP2007248479A (ja) 1997-11-26 2007-07-04 ディジタル流量計

Country Status (4)

Country Link
US (7) US6311136B1 (ja)
EP (1) EP0919793B1 (ja)
JP (2) JP2000018995A (ja)
DE (3) DE69840766D1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008014949A (ja) * 2006-07-04 2008-01-24 Krohne Ag コリオリ質量流量測定装置の作動方法
JP2018508025A (ja) * 2015-03-13 2018-03-22 マイクロ モーション インコーポレイテッド 振動式流量計における信号の温度補償
CN110346006A (zh) * 2019-08-13 2019-10-18 上海一诺仪表有限公司 一种质量流量计变送器和质量流量计

Families Citing this family (220)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102013101369B4 (de) 2013-02-12 2021-02-18 Endress + Hauser Flowtec Ag Coriolis-Massendurchfluß-Meßgerät
US7949495B2 (en) 1996-03-28 2011-05-24 Rosemount, Inc. Process variable transmitter with diagnostics
US7254518B2 (en) * 1996-03-28 2007-08-07 Rosemount Inc. Pressure transmitter with diagnostics
US8290721B2 (en) 1996-03-28 2012-10-16 Rosemount Inc. Flow measurement diagnostics
US6539267B1 (en) 1996-03-28 2003-03-25 Rosemount Inc. Device in a process system for determining statistical parameter
US6017143A (en) 1996-03-28 2000-01-25 Rosemount Inc. Device in a process system for detecting events
US6519546B1 (en) 1996-11-07 2003-02-11 Rosemount Inc. Auto correcting temperature transmitter with resistance based sensor
US6754601B1 (en) 1996-11-07 2004-06-22 Rosemount Inc. Diagnostics for resistive elements of process devices
US6601005B1 (en) * 1996-11-07 2003-07-29 Rosemount Inc. Process device diagnostics using process variable sensor signal
US8467986B2 (en) 1997-11-26 2013-06-18 Invensys Systems, Inc. Drive techniques for a digital flowmeter
US20030216874A1 (en) 2002-03-29 2003-11-20 Henry Manus P. Drive techniques for a digital flowmeter
US6311136B1 (en) * 1997-11-26 2001-10-30 Invensys Systems, Inc. Digital flowmeter
US8447534B2 (en) * 1997-11-26 2013-05-21 Invensys Systems, Inc. Digital flowmeter
US7784360B2 (en) * 1999-11-22 2010-08-31 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
US7404336B2 (en) 2000-03-23 2008-07-29 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
US7124646B2 (en) * 1997-11-26 2006-10-24 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
US6594620B1 (en) * 1998-08-17 2003-07-15 Aspen Technology, Inc. Sensor validation apparatus and method
US6615149B1 (en) 1998-12-10 2003-09-02 Rosemount Inc. Spectral diagnostics in a magnetic flow meter
US6611775B1 (en) 1998-12-10 2003-08-26 Rosemount Inc. Electrode leakage diagnostics in a magnetic flow meter
US7010459B2 (en) * 1999-06-25 2006-03-07 Rosemount Inc. Process device diagnostics using process variable sensor signal
US6505517B1 (en) 1999-07-23 2003-01-14 Rosemount Inc. High accuracy signal processing for magnetic flowmeter
US6701274B1 (en) 1999-08-27 2004-03-02 Rosemount Inc. Prediction of error magnitude in a pressure transmitter
US6543297B1 (en) 1999-09-13 2003-04-08 Rosemount Inc. Process flow plate with temperature measurement feature
US6487507B1 (en) * 1999-10-15 2002-11-26 Micro Motion, Inc. Remote signal conditioner for a Coriolis flowmeter
US6618684B1 (en) * 2000-01-26 2003-09-09 Elster Electricity, Llc System and method for digitally compensating frequency and temperature induced errors in amplitude and phase shift in current sensing of electronic energy meters
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
US6993445B2 (en) * 2001-01-16 2006-01-31 Invensys Systems, Inc. Vortex flowmeter
US6629059B2 (en) 2001-05-14 2003-09-30 Fisher-Rosemount Systems, Inc. Hand held diagnostic and communication device with automatic bus detection
WO2002103327A1 (de) * 2001-06-19 2002-12-27 Endress + Hauser Flowtec Ag Viskositäts-messgerät
US6832179B2 (en) * 2001-06-26 2004-12-14 Invensys Systems, Inc. Evaluating a vortex flow-meter signal
US6772036B2 (en) 2001-08-30 2004-08-03 Fisher-Rosemount Systems, Inc. Control system using process model
US6606917B2 (en) * 2001-11-26 2003-08-19 Emerson Electric Co. High purity coriolis mass flow controller
US20030098069A1 (en) * 2001-11-26 2003-05-29 Sund Wesley E. High purity fluid delivery system
JP3753057B2 (ja) * 2001-12-04 2006-03-08 株式会社日立製作所 気体流量測定装置
JP3669580B2 (ja) * 2002-05-24 2005-07-06 学校法人慶應義塾 超音波流速分布及び流量計
DE10223358A1 (de) * 2002-05-25 2003-12-04 Bosch Gmbh Robert Verfahren und Anordnung zur Erfassung der Bewegung eines Elements
GB2391313B (en) * 2002-07-26 2006-01-18 Abb Ltd Deriving measurements from periodic data
US7040179B2 (en) * 2002-12-06 2006-05-09 Endress+ Hauser Flowtec Ag Process meter
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
US6993396B1 (en) * 2003-03-20 2006-01-31 John Peter Gerry System for determining the health of process control feedback loops according to performance assessment criteria
US6997032B2 (en) * 2003-04-08 2006-02-14 Invensys Systems, Inc. Flowmeter zeroing techniques
US7013740B2 (en) * 2003-05-05 2006-03-21 Invensys Systems, Inc. Two-phase steam measurement system
DE10322763A1 (de) * 2003-05-19 2004-12-09 Helios + Zaschel Gmbh Verfahren und Vorrichtung zur Messung eines Massestroms
CA2470217A1 (en) * 2003-06-06 2004-12-06 Ameriflo, Inc. Lighted fluid flow indication apparatus
US7072775B2 (en) * 2003-06-26 2006-07-04 Invensys Systems, Inc. Viscosity-corrected flowmeter
DE10331126B4 (de) 2003-07-09 2005-09-01 Krohne Ag Coriolis-Massendurchflußmeßgerät und Verfahren zum Betreiben eines Coriolis-Massendurchflußmeßgeräts
US7065455B2 (en) * 2003-08-13 2006-06-20 Invensys Systems, Inc. Correcting frequency in flowtube measurements
AU2003277178B2 (en) * 2003-09-30 2009-04-23 Micro Motion, Inc. Two-wire bus instrument
CN100447535C (zh) * 2003-10-30 2008-12-31 因万西斯系统股份有限公司 流量计的动态响应特性
US7124037B2 (en) * 2003-10-31 2006-10-17 Mpi Incorporated Devices, systems, and methods for monitoring the operation of an injection molding machine
CN100460830C (zh) * 2003-12-12 2009-02-11 罗斯蒙德公司 可调空管功能
US7040181B2 (en) 2004-03-19 2006-05-09 Endress + Hauser Flowtec Ag Coriolis mass measuring device
US20050274200A1 (en) * 2004-05-25 2005-12-15 Henry Manus P Flowmeter batching techniques
US7021121B2 (en) * 2004-05-27 2006-04-04 Asml Holding N.V. Gas gauge proximity sensor with a modulated gas flow
US7302356B2 (en) * 2004-09-15 2007-11-27 Endress + Hauser Flowtec Ag Coriolis flowmeter
DE102004060115A1 (de) * 2004-12-13 2006-06-14 Endress + Hauser Flowtec Ag Meßaufnehmer vom Vibrationstyp
US7350421B2 (en) * 2004-12-13 2008-04-01 Endress + Hauser Flowtec Ag Vibratory measurement transducer
BRPI0519594B1 (pt) * 2004-12-29 2016-09-20 Micro Motion Inc componentes eletrônicos de medidor e método para processar sinais de sensor em um medidor de fluxo
DE102005012505B4 (de) * 2005-02-16 2006-12-07 Krohne Ag Verfahren zum Betreiben eines Massendurchflußmeßgeräts
US7064692B1 (en) * 2005-03-18 2006-06-20 Honeywell International Inc. Solid-state synchro/resolver converter
RU2371677C2 (ru) * 2005-03-29 2009-10-27 Майкро Моушн, Инк. Измерительное электронное устройство и способ для определения жидкой фракции потока в материале газового потока
US8112565B2 (en) 2005-06-08 2012-02-07 Fisher-Rosemount Systems, Inc. Multi-protocol field device interface with automatic bus detection
US7337084B2 (en) 2005-06-21 2008-02-26 Invensys Systems, Inc. Switch-activated zero checking feature for a Coriolis flowmeter
US7313488B2 (en) * 2005-07-11 2007-12-25 Invensys Systems, Inc. Coriolis mode processing techniques
EP1925916A3 (en) 2005-07-11 2011-08-03 Invensys Systems, Inc. Coriolis mode processing techniques
DE102005046319A1 (de) 2005-09-27 2007-03-29 Endress + Hauser Flowtec Ag Verfahren zum Messen eines in einer Rohrleitung strömenden Mediums sowie Meßsystem dafür
US20070068225A1 (en) 2005-09-29 2007-03-29 Brown Gregory C Leak detector for process valve
PL1949047T3 (pl) * 2005-10-18 2015-08-31 Micro Motion Inc Układ elektroniczny miernika i sposoby określania różnicy fazy między pierwszym sygnałem czujnika i drugim sygnałem czujnika przepływomierza
US7505947B2 (en) * 2005-10-20 2009-03-17 International Business Machines Corporation Computer controlled method using genetic algorithms to provide non-deterministic solutions to problems involving physical restraints
US7360452B2 (en) * 2005-12-27 2008-04-22 Endress + Hauser Flowtec Ag In-line measuring devices and method for compensation measurement errors in in-line measuring devices
US7360453B2 (en) * 2005-12-27 2008-04-22 Endress + Hauser Flowtec Ag In-line measuring devices and method for compensation measurement errors in in-line measuring devices
WO2007074055A1 (en) 2005-12-27 2007-07-05 Endress+Hauser Flowtec Ag In-line measuring devices and method for compensating measurement errors in in-line measuring devices
US7532992B2 (en) * 2006-01-20 2009-05-12 Teledyne Isco, Inc. Measuring apparatuses and methods of using them
EP1989517B1 (en) * 2006-02-13 2015-09-02 Invensys Systems, Inc. Compensating for frequency change in flowmeters
EP2027688B1 (en) * 2006-04-28 2018-06-06 Micro Motion, Inc. Bus loop power interface and method
JP2008027263A (ja) * 2006-07-24 2008-02-07 Oval Corp 演算システムでの処理方法、流量変換器、及びコリオリ流量計
US7617055B2 (en) 2006-08-28 2009-11-10 Invensys Systems, Inc. Wet gas measurement
US7953501B2 (en) 2006-09-25 2011-05-31 Fisher-Rosemount Systems, Inc. Industrial process control loop monitor
JP2010505121A (ja) 2006-09-29 2010-02-18 ローズマウント インコーポレイテッド 検証を備える磁気流量計
DE102006058732A1 (de) 2006-12-13 2008-06-26 Abb Ag Verfahren zum Betrieb einer Coriolis-Massendurchflussmesseinrichtung
DE102006062600B4 (de) 2006-12-29 2023-12-21 Endress + Hauser Flowtec Ag Verfahren zum Inbetriebnehmen und/oder Überwachen eines In-Line-Meßgeräts
AT9241U3 (de) * 2007-02-05 2007-12-15 Avl List Gmbh Verfahren und vorrichtung zur kontinuierlichen messung eines dynamischen fluidverbrauchs
US8751171B2 (en) * 2007-03-07 2014-06-10 Invensys Systems, Inc. Coriolis frequency tracking
DE102007021099A1 (de) 2007-05-03 2008-11-13 Endress + Hauser (Deutschland) Ag + Co. Kg Verfahren zum Inbetriebnehmen und/oder Rekonfigurieren eines programmierbaren Feldmeßgeräts
US8289179B2 (en) * 2007-05-25 2012-10-16 Micro Motion, Inc. Vibratory flow meter and method for correcting for entrained gas in a flow material
US8898036B2 (en) 2007-08-06 2014-11-25 Rosemount Inc. Process variable transmitter with acceleration sensor
US7832257B2 (en) 2007-10-05 2010-11-16 Halliburton Energy Services Inc. Determining fluid rheological properties
DE102007058608A1 (de) 2007-12-04 2009-06-10 Endress + Hauser Flowtec Ag Elektrisches Gerät
US7895006B2 (en) * 2007-12-21 2011-02-22 Agilent Technologies, Inc. System and method for determining signal-to-noise ratio (SNR), signal-to-noise and distortion ratio (SINAD) and total harmonic distortion (THD)
DE102008011381A1 (de) * 2008-02-27 2009-09-03 Endress + Hauser Flowtec Ag Verfahren zur Erzeugung eines annähernd sinusförmigen elektrischen Spannungsverlaufs und Verfahren zur Erzeugung mindestens einer Schwingung eines mechanisch schwingfähigen Körpers
DE102008016235A1 (de) 2008-03-27 2009-10-01 Endress + Hauser Flowtec Ag Verfahren zum Betreiben eines auf einer rotierenden Karussell-Abfüllmachine angeordneten Meßgeräts
WO2009134830A1 (en) * 2008-05-01 2009-11-05 Micro Motion, Inc. Vibratory flow meter for determining one or more flow fluid characteristics of a multi-phase flow fluid
DE102008022373A1 (de) 2008-05-06 2009-11-12 Endress + Hauser Flowtec Ag Meßgerät sowie Verfahren zum Überwachen eines Meßgeräts
CN102084224B (zh) * 2008-07-01 2017-04-19 微动公司 用于生成振动测量设备中的驱动信号的系统、方法
DE102008035877A1 (de) 2008-08-01 2010-02-04 Endress + Hauser Flowtec Ag Meßwandler vom Vibrationstyp
DE102008050116A1 (de) 2008-10-06 2010-04-08 Endress + Hauser Flowtec Ag In-Line-Meßgerät
DE102008050113A1 (de) 2008-10-06 2010-04-08 Endress + Hauser Flowtec Ag In-Line-Meßgerät
DE102008050115A1 (de) 2008-10-06 2010-04-08 Endress + Hauser Flowtec Ag In-Line-Meßgerät
DE102008044186A1 (de) 2008-11-28 2010-06-02 Endress + Hauser Flowtec Ag Magneteinrichtung sowie Meßaufnehmer vom Vibrationstyp mit einer solchen Magneteinrichtung
JP4436882B1 (ja) * 2009-02-06 2010-03-24 株式会社オーバル 信号処理方法、信号処理装置、およびコリオリ流量計
JP4436884B1 (ja) * 2009-02-06 2010-03-24 株式会社オーバル 信号処理方法、信号処理装置、およびコリオリ流量計
JP4436883B1 (ja) * 2009-02-06 2010-03-24 株式会社オーバル 信号処理方法、信号処理装置、およびコリオリ流量計
CN102348962B (zh) 2009-03-11 2015-12-02 恩德斯+豪斯流量技术股份有限公司 振动类型测量换能器以及具有该测量换能器的在线测量装置
DE102009027580A1 (de) 2009-07-09 2011-01-13 Endress + Hauser Flowtec Ag Meßaufnehmer vom Vibrationstyp sowie In-line-Meßgerät mit einem solchen Meßaufnehmer
DE102009001472A1 (de) 2009-03-11 2010-09-16 Endress + Hauser Flowtec Ag Meßaufnehmer vom Vibrationstyp sowie In-line-Meßgerät mit einem solchen Meßaufnehmer
DE102009012474A1 (de) 2009-03-12 2010-09-16 Endress + Hauser Flowtec Ag Meßsystem mit einem Messwandler vom Vibrationstyp
US7921734B2 (en) 2009-05-12 2011-04-12 Rosemount Inc. System to detect poor process ground connections
EP2454563A1 (en) * 2009-07-13 2012-05-23 Micro Motion, Inc. Meter electronics and fluid quantification method for a fluid being transferred
CN101608940B (zh) * 2009-07-23 2011-12-21 合肥工业大学 科氏质量流量计正负阶跃交替激励启振方法和系统
DE102009028006A1 (de) 2009-07-24 2011-01-27 Endress + Hauser Flowtec Ag Meßwandler vom Vibrationstyp sowie Meßgerät mit einem solchen Meßwandler
DE102009028007A1 (de) 2009-07-24 2011-01-27 Endress + Hauser Flowtec Ag Meßumwandler vom Vibrationstyp sowie Meßgerät mit einem solchen Meßwandler
CA2770135C (en) * 2009-08-12 2016-06-07 Micro Motion, Inc. Method and apparatus for determining a zero offset in a vibrating flow meter
MX2012001688A (es) * 2009-08-12 2012-03-07 Micro Motion Inc Metodo y aparato para determinar y compensar para un cambio en una diferencial de ajuste a cero de un flujometro vibratorio.
DE102009046839A1 (de) 2009-11-18 2011-05-19 Endress + Hauser Flowtec Ag Meßsystem mit einer zwei parallel durchströmte Meßrohre aufweisenden Rohranordnung sowie Verfahren zu deren Überwachung
DE102009055069A1 (de) 2009-12-21 2011-06-22 Endress + Hauser Flowtec Ag Meßaufnehmer vom Vibrationstyp
DE102010039627A1 (de) 2010-08-20 2012-02-23 Endress + Hauser Flowtec Ag Meßaufnehmer vom Vibrationstyp sowie damit gebildetes Meßsystem
WO2011085852A1 (de) 2009-12-21 2011-07-21 Endress+Hauser Flowtec Ag Messaufnehmer vom vibrationstyp sowie damit gebildetes messsystem
WO2011080172A1 (de) 2009-12-31 2011-07-07 Endress+Hauser Flowtec Ag MEßSYSTEM MIT EINEM MEßWANDLER VOM VIBRATIONSTYP
CA2785933C (en) 2009-12-31 2016-05-24 Endress+Hauser Flowtec Ag Measuring system having a measuring transducer of vibration-type
DE102010000760B4 (de) 2010-01-11 2021-12-23 Endress + Hauser Flowtec Ag Meßsystem mit einem Meßwandler vom Vibrationstyp zum Messen eines statischen Drucks in einem strömenden Medium
DE102010000759A1 (de) 2010-01-11 2011-07-14 Endress + Hauser Flowtec Ag Meßsystem mit einem Meßwandler vom Vibrationstyp
WO2011080171A2 (de) 2009-12-31 2011-07-07 Endress+Hauser Flowtec Ag MEßSYSTEM MIT EINEM MEßWANDLER VOM VIBRATIONSTYP
DE102010000761A1 (de) 2010-01-11 2011-07-28 Endress + Hauser Flowtec Ag Meßsystem mit einem Meßwandler vom Vibrationstyp
JP4694645B1 (ja) * 2010-02-19 2011-06-08 株式会社オーバル 信号処理方法、信号処理装置、及び振動型密度計
CN102859852B (zh) 2010-04-19 2015-11-25 恩德斯+豪斯流量技术股份有限公司 测量变换器的驱动电路及由该驱动电路形成的测量系统
DE202010006553U1 (de) 2010-05-06 2011-10-05 Endress + Hauser Flowtec Ag Elektronisches Meßgerät mit einem Optokoppler
DE102010030924A1 (de) 2010-06-21 2011-12-22 Endress + Hauser Flowtec Ag Elektronik-Gehäuse für ein elektronisches Gerät bzw. damit gebildetes Gerät
DE102010039543A1 (de) 2010-08-19 2012-02-23 Endress + Hauser Flowtec Ag Meßsystem mit einem Meßwandler vom Vibrationstyp
EP2612114B1 (de) 2010-09-02 2019-10-30 Endress+Hauser Flowtec AG Messsystem mit einem messaufnehmer vom vibrationstyp
RU2557409C2 (ru) 2010-09-16 2015-07-20 Эндресс+Хаузер Флоутек Аг Измерительная система для измерения плотности или весовой пропускной способности протекающей в трубопроводе среды
DE102010044179A1 (de) 2010-11-11 2012-05-16 Endress + Hauser Flowtec Ag Meßsystem mit einem Meßwandler von Vibrationstyp
CA2758779C (en) 2010-12-01 2016-01-05 Invensys Systems, Inc. Determining concentrations of components of a mixture
WO2012089431A1 (de) 2010-12-30 2012-07-05 Endress+Hauser Flowtec Ag Messaufnehmer vom vibrationstyp sowie damit gebildetes messsystem
US8593273B2 (en) 2011-02-07 2013-11-26 Infineon Technologies Ag Systems and methods for localization of tire pressure monitoring system wheel modules
CN102128656B (zh) * 2011-02-25 2013-09-04 合肥工业大学 一种微弯型科氏质量流量计数字信号处理方法和系统
US9207670B2 (en) 2011-03-21 2015-12-08 Rosemount Inc. Degrading sensor detection implemented within a transmitter
DE102011006997A1 (de) 2011-04-07 2012-10-11 Endress + Hauser Flowtec Ag Frequenzabgleichsverfahren für eine Rohranordnung
DE102011006971A1 (de) 2011-04-07 2012-10-11 Endress + Hauser Flowtec Ag Meßwandler vom Vibrationstyp sowie Verfahren zu dessen Herstellung
DE102011006919A1 (de) 2011-04-07 2012-10-11 Endress + Hauser Flowtec Ag Verfahren zum Trimmen eines Rohrs
JP5836628B2 (ja) * 2011-04-19 2015-12-24 キヤノン株式会社 制御系の評価装置および評価方法、並びに、プログラム
EP2705334B1 (de) 2011-05-02 2018-01-24 Endress+Hauser Flowtec AG Messaufnehmer vom vibrationstyp sowie damit gebildetes messsystem
GB201108854D0 (en) * 2011-05-26 2011-07-06 Spp Process Technology Systems Uk Ltd Mass flow controller monitoring
DE102011076838A1 (de) 2011-05-31 2012-12-06 Endress + Hauser Flowtec Ag Meßgerät-Elektronik für ein Meßgerät-Gerät sowie damit gebildetes Meßgerät-Gerät
EP2565594A1 (en) 2011-07-11 2013-03-06 Invensys Systems, Inc. Coriolis flowmeter with zero checking feature
US20130080084A1 (en) * 2011-09-28 2013-03-28 John P. Miller Pressure transmitter with diagnostics
US8565967B2 (en) * 2011-12-21 2013-10-22 Infineon Technologies Ag Acceleration detection and angular position determination systems and methods in tire pressure monitoring systems
DE102011089808A1 (de) 2011-12-23 2013-06-27 Endress + Hauser Flowtec Ag Verfahren bzw. Meßsystem zum Ermitteln einer Dichte eines Fluids
RU2579818C1 (ru) 2012-04-03 2016-04-10 Эндресс + Хаузер Флоутек Аг Измерительный преобразователь вибрационного типа, измерительная система для протекающей через трубопровод среды и способ постройки частоты системы труб
DE102012102947B4 (de) 2012-04-03 2023-12-21 Endress + Hauser Flowtec Ag Meßwandler vom Vibrationstyp
JP6166353B2 (ja) * 2012-04-23 2017-07-19 ローズマウント インコーポレイテッド プロセス送信機におけるプロセス変数補償
EP2861941B1 (de) * 2012-06-18 2018-03-21 KROHNE Messtechnik GmbH Verfahren zum betreiben eines resonanzmesssystems
US9052240B2 (en) 2012-06-29 2015-06-09 Rosemount Inc. Industrial process temperature transmitter with sensor stress diagnostics
US9602122B2 (en) 2012-09-28 2017-03-21 Rosemount Inc. Process variable measurement noise diagnostic
DE102012109729A1 (de) 2012-10-12 2014-05-15 Endress + Hauser Flowtec Ag Meßsystem zum Ermitteln eines Volumendruchflusses und/oder einer Volumendurchflußrate eines in einer Rohrleitung strömenden Mediums
US9372107B2 (en) 2012-10-11 2016-06-21 Endress + Hauser Flowtec Ag Measuring system for ascertaining a volume flow and/or a volume flow rate of a medium flowing in a pipeline
US9995666B2 (en) * 2012-10-22 2018-06-12 Rheonics Gmbh Resonant sensors for fluid properties measurement
CN104040301B (zh) * 2012-10-30 2017-12-12 西安东风机电股份有限公司 一种科里奥利流量计数字驱动方法及系统
CN103959019B (zh) 2012-11-19 2017-03-22 因文西斯系统公司 净油气井测试系统
DE102013102708A1 (de) 2013-03-18 2014-09-18 Endress + Hauser Flowtec Ag Meßwandler vom Vibrationstyp sowie damit gebildetes Meßsystem
DE102013102711A1 (de) 2013-03-18 2014-09-18 Endress + Hauser Flowtec Ag Meßwandler vom Vibrationstyp sowie damit gebildetes Meßsystem
DE102013106155A1 (de) 2013-06-13 2014-12-18 Endress + Hauser Flowtec Ag Meßsystem mit einem Druckgerät sowie Verfahren zur Überwachung und/oder Überprüfung eines solchen Druckgeräts
DE102013106157A1 (de) 2013-06-13 2014-12-18 Endress + Hauser Flowtec Ag Meßsystem mit einem Druckgerät sowie Verfahren zur Überwachung und/oder Überprüfung eines solchen Druckgeräts
KR101889831B1 (ko) * 2013-11-14 2018-08-21 마이크로 모우션, 인코포레이티드 코리올리 다이렉트 웰헤드 측정 디바이스들 및 방법들
JP2015105929A (ja) * 2013-12-02 2015-06-08 株式会社東芝 電磁流量計
DE102013113689B4 (de) 2013-12-09 2018-02-01 Endress + Hauser Flowtec Ag Dichte-Meßgerät
WO2015090776A1 (de) 2013-12-20 2015-06-25 Endress+Hauser Flowtec Ag Spule
DE102013114731A1 (de) 2013-12-20 2015-06-25 Endress+Hauser Flowtec Ag Spule
DE102013021915A1 (de) 2013-12-27 2015-07-02 Endress + Hauser Flowtec Ag Meßaufnehmer vom Vibrationstyp
DE102014103430A1 (de) 2014-03-13 2015-09-17 Endress + Hauser Flowtec Ag Wandlervorrichtung sowie damit gebildetes Meßsystem
DE102014103427A1 (de) 2014-03-13 2015-09-17 Endress + Hauser Flowtec Ag Wandlervorrichtung sowie damit gebildetes Meßsystem
AU2014405569C1 (en) * 2014-09-04 2019-06-20 Micro Motion, Inc. Differential flowmeter tool
DE102015103208A1 (de) 2014-10-17 2016-04-21 Endress + Hauser Flowtec Ag Meßsystem zum Messen wenigstens einer Meßgröße eines Fluids sowie Verfahren zum Betreiben eines solchen Meßsystems
US9689736B2 (en) 2014-10-31 2017-06-27 Invensys Systems, Inc. Method to provide a quality measure for meter verification results
DE102014119212A1 (de) * 2014-12-19 2016-06-23 Endress + Hauser Flowtec Ag Messanordnung und Verfahren zum Messen der Dichte von fließfähigen Medien
WO2016107693A1 (de) 2014-12-30 2016-07-07 Endress+Hauser Flowtec Ag Verfahren zum messen einer dichte eines fluids
DE102014019396A1 (de) 2014-12-30 2016-06-30 Endress+Hauser Flowtec Ag Verfahren zum Messen einer Dichte eines Fluids
JP6036864B2 (ja) * 2015-02-05 2016-11-30 横河電機株式会社 測定装置の共振回路
US9863798B2 (en) 2015-02-27 2018-01-09 Schneider Electric Systems Usa, Inc. Systems and methods for multiphase flow metering accounting for dissolved gas
US9664548B2 (en) 2015-03-19 2017-05-30 Invensys Systems, Inc. Testing system for petroleum wells having a fluidic system including a gas leg, a liquid leg, and bypass conduits in communication with multiple multiphase flow metering systems with valves to control fluid flow through the fluidic system
WO2017019024A1 (en) * 2015-07-27 2017-02-02 Micro Motion, Inc. Off-resonance cycling for coriolis flowmeters
US10408647B2 (en) * 2016-02-15 2019-09-10 Olea Networks, Inc. Analysis of pipe systems with sensor devices
US10024707B2 (en) * 2016-02-17 2018-07-17 Schneider Electric Systems Usa, Inc. Electromagnetic flowmeter calibration verification
RU2721312C2 (ru) * 2016-02-26 2020-05-18 Майкро Моушн, Инк. Электронный измеритель для двух или более измерительных узлов
DE102016112599A1 (de) 2016-07-08 2018-01-11 Endress + Hauser Flowtec Ag Meßsystem
DE102016112600A1 (de) 2016-07-08 2018-01-11 Endress + Hauser Flowtec Ag Meßsystem
DE102016114860A1 (de) 2016-08-10 2018-02-15 Endress + Hauser Flowtec Ag Treiberschaltung sowie damit gebildete Umformer-Elektronik bzw. damit gebildetes Meßsystem
WO2018087344A2 (en) 2016-11-11 2018-05-17 Oxford University Innovation Limited Method and system for tracking sinusoidal wave parameters from a received signal that includes noise
DE102017106211A1 (de) 2016-12-29 2018-07-05 Endress+Hauser Flowtec Ag Vibronisches Meßsystem zum Messen einer Massendurchflußrate
WO2018121929A1 (de) 2016-12-29 2018-07-05 Endress+Hauser Flowtec Ag VIBRONISCHES MEßSYSTEM ZUM MESSEN EINER MASSENDURCHFLUßRATE
CN110114641B (zh) 2016-12-29 2021-08-03 恩德斯+豪斯流量技术股份有限公司 用于测量质量流率的电子振动测量系统
CN106643955B (zh) * 2017-02-20 2023-04-07 重庆机床(集团)有限责任公司 流量监视控制器
DE102017118109A1 (de) 2017-08-09 2019-02-14 Endress + Hauser Flowtec Ag Sensorbaugruppe
DE102017121157A1 (de) 2017-08-09 2019-02-14 Endress+Hauser Flowtec Ag Spule sowie Meßwandler mit einer solchen Spule
DE102018102831A1 (de) 2017-10-05 2019-04-11 Endress+Hauser Flowtec Ag Meßwandler für ein vibronisches Meßsystem sowie damit gebildetes vibronisches Meßsystem
DE102017125271A1 (de) * 2017-10-27 2019-05-02 Endress + Hauser Flowtec Ag Massedurchflussmessgerät nach dem Coriolis-Prinzip mit mindestens zwei Messrohrpaaren
DE102017125273A1 (de) 2017-10-27 2019-05-02 Endress + Hauser Flowtec Ag Massedurchflussmessgerät nach dem Coriolis-Prinzip mit mindestens zwei Messrohrpaaren und Verfahren zum Bestimmen des Massedurchflusses
DE102017131199A1 (de) 2017-12-22 2019-06-27 Endress + Hauser Flowtec Ag Coriolis-Massendurchfluß-Meßgerät
US11237550B2 (en) * 2018-03-28 2022-02-01 Honeywell International Inc. Ultrasonic flow meter prognostics with near real-time condition based uncertainty analysis
US11781673B2 (en) * 2018-04-30 2023-10-10 Keto A.I., Inc. Water level control system
DE102018133117A1 (de) 2018-12-20 2020-06-25 Endress+Hauser Flowtec Ag Coriolis-Massendurchfluß-Meßgerät
US20220099543A1 (en) 2018-12-20 2022-03-31 Endress+Hauser Flowtec Ag Coriolis mass flow meter
WO2020126282A1 (de) 2018-12-20 2020-06-25 Endress+Hauser Flowtec Ag CORIOLIS-MASSENDURCHFLUß-MEßGERÄT
US20220057245A1 (en) 2018-12-21 2022-02-24 Endress+Hauser Flowtec Ag Coriolis mass flowmeter with magnetic field detector
US20200209031A1 (en) * 2018-12-26 2020-07-02 Texas Instruments Incorporated Dynamic temperature calibration of ultrasonic transducers
US20220099466A1 (en) * 2019-01-28 2022-03-31 The Texas A&M University System Method and device to measure multiphase flow
CN110470347B (zh) * 2019-09-12 2021-12-17 上海旗浪节能环保科技有限公司 一种用于海洋轮船空调系统的管道实时监测装置
DE102019124709A1 (de) * 2019-09-13 2021-03-18 Endress+Hauser Flowtec Ag Verfahren zum Betreiben eines Messgerätes mit mindestens einem Oszillator und Messgerät zur Durchführung des Verfahrens
DE102019133610A1 (de) 2019-12-09 2021-06-10 Endress + Hauser Flowtec Ag Vibronisches Meßsystem zum Messen eines Massestroms eines fluiden Meßstoff
DE102019009024A1 (de) 2019-12-30 2021-07-01 Endress+Hauser Flowtec Ag Vibronisches Meßsystem
EP4168752A1 (de) 2020-06-18 2023-04-26 Endress + Hauser Flowtec AG VIBRONISCHES MEßSYSTEM
DE102020120054A1 (de) 2020-07-29 2022-02-03 Endress + Hauser Flowtec Ag Verfahren zum Ermitteln einer Meßstoff-Temperatur sowie Meßsystem dafür
DE102020131649A1 (de) 2020-09-03 2022-03-03 Endress + Hauser Flowtec Ag Vibronisches Meßsystem
DE102020127382A1 (de) 2020-10-16 2022-04-21 Endress+Hauser Flowtec Ag Verfahren zum Überprüfen eines vibronischen Meßsystems
DE102021113360A1 (de) 2021-05-21 2022-11-24 Endress + Hauser Flowtec Ag Vibronisches Meßsystem
DE102021123412A1 (de) * 2021-09-09 2023-03-09 Endress+Hauser Flowtec Ag Vibronischer Messaufnehmer zur Massedurchfluss- und Dichtemessung
DE102022112523A1 (de) 2022-05-18 2023-11-23 Endress+Hauser Flowtec Ag Vibronisches Meßsystem
DE102022116111A1 (de) 2022-06-28 2023-12-28 Endress+Hauser Flowtec Ag Vibronisches Meßsystem
DE102022119145A1 (de) 2022-07-29 2024-02-01 Endress+Hauser Flowtec Ag Anschlussschaltung für ein Feldgerät und Feldgerät

Family Cites Families (90)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US29383A (en) * 1860-07-31 Improvement in preserving food
US31450A (en) * 1861-02-19 Improvement in tools used in the manufacture of iron
US5732193A (en) * 1909-01-26 1998-03-24 Aberson; Michael Method and apparatus for behavioristic-format coding of quantitative resource data/distributed automation protocol
US4027662A (en) * 1973-07-11 1977-06-07 Milstein Medical Research Foundation, Inc. Automatic blood pressure recorder
USRE29383E (en) 1974-01-10 1977-09-06 Process Systems, Inc. Digital fluid flow rate measurement or control system
US3956682A (en) * 1974-03-28 1976-05-11 Pennwalt Corporation Two-wire position-to-D.C. current transducer
GB1583545A (en) * 1976-08-04 1981-01-28 Martin Sanchez J Control systems
US4187721A (en) * 1977-07-25 1980-02-12 S & F Associates Method and structure for flow measurement
USRE31450E (en) 1977-07-25 1983-11-29 Micro Motion, Inc. Method and structure for flow measurement
GB2085597B (en) * 1980-10-17 1985-01-30 Redland Automation Ltd Method and apparatus for detemining the mass flow of a fluid
US4422338A (en) 1981-02-17 1983-12-27 Micro Motion, Inc. Method and apparatus for mass flow measurement
US4491025A (en) 1982-11-03 1985-01-01 Micro Motion, Inc. Parallel path Coriolis mass flow rate meter
DE3464189D1 (en) 1984-08-03 1987-07-16 Itt Ind Gmbh Deutsche Videorecorder with a magnetic tape storage medium
US4773257A (en) 1985-06-24 1988-09-27 Chevron Research Company Method and apparatus for testing the outflow from hydrocarbon wells on site
US4688418A (en) 1985-10-17 1987-08-25 Texaco Inc. Method and apparatus for determining mass flow rate and quality in a steam line
JPS62157546A (ja) * 1985-12-31 1987-07-13 Horiba Ltd 自動車排気ガスのモ−ダルマス解析方法
US4879111A (en) 1986-04-17 1989-11-07 Cetus Corporation Treatment of infections with lymphokines
US4823614A (en) * 1986-04-28 1989-04-25 Dahlin Erik B Coriolis-type mass flowmeter
GB8614135D0 (en) 1986-06-10 1986-07-16 Foxboro Co Coriolis mass flowmeters
US4817448A (en) * 1986-09-03 1989-04-04 Micro Motion, Inc. Auto zero circuit for flow meter
DE8712331U1 (de) 1986-09-26 1988-01-28 Flowtec AG, Reinach, Basel Corioliskraft-Massendurchflussmesser
US4852410A (en) 1986-10-03 1989-08-01 Schlumberger Industries, Inc. Omega-shaped, coriolis-type mass flow rate meter
US4911006A (en) 1986-10-03 1990-03-27 Micro Motion Incorporated Custody transfer meter
US4782711A (en) * 1986-10-14 1988-11-08 K-Flow Division Of Kane Steel Co., Inc. Method and apparatus for measuring mass flow
KR960000099B1 (ko) 1986-10-28 1996-01-03 더폭스보로 컴패니 코리올리 유형의 질량유량계
US5271281A (en) 1986-10-28 1993-12-21 The Foxboro Company Coriolis-type mass flowmeter
US5343764A (en) 1986-10-28 1994-09-06 The Foxboro Company Coriolis-type mass flowmeter
US5050439A (en) 1986-10-28 1991-09-24 The Foxboro Company Coriolis-type mass flowmeter circuitry
US4911020A (en) 1986-10-28 1990-03-27 The Foxboro Company Coriolis-type mass flowmeter circuitry
IT1213434B (it) 1986-12-23 1989-12-20 Nuovo Pignone S P A Ind Meccan Procedimento perfezionato per la misura di portate ponderali e relativi dispositivi.
US5027662A (en) 1987-07-15 1991-07-02 Micro Motion, Inc. Accuracy mass flow meter with asymmetry and viscous damping compensation
US5052231A (en) 1988-05-19 1991-10-01 Rheometron Ag Mass flow gauge for flowing media with devices for determination of the Coriolis force
US4852409A (en) * 1988-06-09 1989-08-01 Fischer & Porter Company Signal recovery system for mass flowmeter
US4879911A (en) 1988-07-08 1989-11-14 Micro Motion, Incorporated Coriolis mass flow rate meter having four pulse harmonic rejection
US4852395A (en) * 1988-12-08 1989-08-01 Atlantic Richfield Company Three phase fluid flow measuring system
US5029482A (en) 1989-02-03 1991-07-09 Chevron Research Company Gas/liquid flow measurement using coriolis-based flow meters
US4934196A (en) * 1989-06-02 1990-06-19 Micro Motion, Inc. Coriolis mass flow rate meter having a substantially increased noise immunity
US4996871A (en) 1989-06-02 1991-03-05 Micro Motion, Inc. Coriolis densimeter having substantially increased noise immunity
AU625788B2 (en) * 1989-06-09 1992-07-16 Micro Motion, Inc. Improved stability coriolis mass flow meter
US5054326A (en) 1990-03-05 1991-10-08 The Foxboro Company Density compensator for coriolis-type mass flowmeters
US5259250A (en) 1990-05-14 1993-11-09 Atlantic Richfield Company Multi-phase fluid flow mesurement
US5054313A (en) * 1990-07-17 1991-10-08 National Metal And Refining Company, Ltd. Control circuitry for viscosity sensors
US5497665A (en) 1991-02-05 1996-03-12 Direct Measurement Corporation Coriolis mass flow rate meter having adjustable pressure and density sensitivity
JP2525085B2 (ja) * 1991-02-19 1996-08-14 株式会社クラレ イロンの製造方法
EP0537347B1 (en) * 1991-03-30 1997-01-22 OKADA, Kazuhiro Acceleration sensor having a self test feature
US5228327A (en) 1991-07-11 1993-07-20 Micro Motion, Inc. Technique for determining a mechanical zero value for a coriolis meter
DE69210244T2 (de) * 1991-08-01 1996-10-31 Micro Motion Inc., Boulder, Col. Massendurchflussmesser nach dem Coriolsprinzip
US5295084A (en) * 1991-10-08 1994-03-15 Micromotion, Inc. Vibrating tube densimeter
US5379649A (en) 1991-12-23 1995-01-10 Micro Motion, Inc. Coriolis effect meter using optical fiber sensors
US5218869A (en) 1992-01-14 1993-06-15 Diasonics, Inc. Depth dependent bandpass of ultrasound signals using heterodyne mixing
GB9208704D0 (en) 1992-04-22 1992-06-10 Foxboro Ltd Improvements in and relating to sensor units
US5347874A (en) 1993-01-25 1994-09-20 Micro Motion, Incorporated In-flow coriolis effect mass flowmeter
US5309342A (en) * 1993-03-02 1994-05-03 Cooper Industries, Inc. Recessed lighting fixture
US5774378A (en) 1993-04-21 1998-06-30 The Foxboro Company Self-validating sensors
WO1995010028A1 (en) 1993-10-05 1995-04-13 Atlantic Richfield Company Multiphase flowmeter for measuring flow rates and densities
JP2651891B2 (ja) * 1993-11-05 1997-09-10 株式会社日本開発コンサルタント 円形型金属性帯板張力付与装置
JPH07169833A (ja) * 1993-12-14 1995-07-04 Nec Corp 半導体装置及びその製造方法
US5429002A (en) 1994-05-11 1995-07-04 Schlumberger Industries, Inc. Coriolis-type fluid mass flow rate measurement device and method employing a least-squares algorithm
JP3219122B2 (ja) 1994-07-11 2001-10-15 横河電機株式会社 コリオリ質量流量計
US5469748A (en) 1994-07-20 1995-11-28 Micro Motion, Inc. Noise reduction filter system for a coriolis flowmeter
US5497666A (en) 1994-07-20 1996-03-12 Micro Motion, Inc. Increased sensitivity coriolis effect flowmeter using nodal-proximate sensors
US5594180A (en) 1994-08-12 1997-01-14 Micro Motion, Inc. Method and apparatus for fault detection and correction in Coriolis effect mass flowmeters
EP0698783A1 (de) 1994-08-16 1996-02-28 Endress + Hauser Flowtec AG Auswerte-Elektronik eines Coriolis-Massedurchflussaufnehmers
US5767665A (en) 1994-09-13 1998-06-16 Fuji Electric Co. Ltd. Phase difference measuring apparatus and mass flowmeter thereof
US5555190A (en) 1995-07-12 1996-09-10 Micro Motion, Inc. Method and apparatus for adaptive line enhancement in Coriolis mass flow meter measurement
US5654502A (en) 1995-12-28 1997-08-05 Micro Motion, Inc. Automatic well test system and method of operating the same
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
US5676461A (en) * 1996-03-18 1997-10-14 M. A. Hanna Rubber Compounding A Division Of M. A. Hanna Company Oil injection apparatus and method for polymer processing
US5877954A (en) * 1996-05-03 1999-03-02 Aspen Technology, Inc. Hybrid linear-neural network process control
US5687100A (en) * 1996-07-16 1997-11-11 Micro Motion, Inc. Vibrating tube densimeter
US5734112A (en) * 1996-08-14 1998-03-31 Micro Motion, Inc. Method and apparatus for measuring pressure in a coriolis mass flowmeter
US5804741A (en) 1996-11-08 1998-09-08 Schlumberger Industries, Inc. Digital phase locked loop signal processing for coriolis mass flow meter
US6073495A (en) 1997-03-21 2000-06-13 Endress + Hauser Flowtec Ag Measuring and operating circuit of a coriolis-type mass flow meter
US6185470B1 (en) 1997-11-07 2001-02-06 Mcdonnell Douglas Corporation Neural network predictive control method and system
US6311136B1 (en) * 1997-11-26 2001-10-30 Invensys Systems, Inc. Digital flowmeter
US6092429A (en) 1997-12-04 2000-07-25 Micro Motion, Inc. Driver for oscillating a vibrating conduit
US5969495A (en) * 1997-12-31 1999-10-19 Daewood Heavy Industries Ltd. Accelerator device for electromotive vehicles
US6102846A (en) 1998-02-26 2000-08-15 Eastman Kodak Company System and method of managing a psychological state of an individual using images
DE59907306D1 (de) * 1998-06-12 2003-11-13 Rieter Ag Maschf Fadenchangierung
EP1112532B1 (en) 1998-08-17 2003-04-02 Aspen Technology, Inc. Sensor validation apparatus and method
US6327914B1 (en) * 1998-09-30 2001-12-11 Micro Motion, Inc. Correction of coriolis flowmeter measurements due to multiphase flows
US5969264A (en) * 1998-11-06 1999-10-19 Technology Commercialization Corp. Method and apparatus for total and individual flow measurement of a single-or multi-phase medium
US6301973B1 (en) * 1999-04-30 2001-10-16 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Non-intrusive pressure/multipurpose sensor and method
US6318186B1 (en) * 1999-06-28 2001-11-20 Micro Motion, Inc. Type identification and parameter selection for drive control in a coriolis flowmeter
US6505131B1 (en) * 1999-06-28 2003-01-07 Micro Motion, Inc. Multi-rate digital signal processor for signals from pick-offs on a vibrating conduit
US6318156B1 (en) * 1999-10-28 2001-11-20 Micro Motion, Inc. Multiphase flow measurement system
US6551251B2 (en) * 2000-02-14 2003-04-22 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Passive fetal heart monitoring system
DE01918944T1 (de) * 2000-03-23 2004-10-21 Invensys Systems, Inc., Foxboro Korrektur für eine zweiphasenströmung in einem digitalen durchflussmesser
US6378354B1 (en) 2000-07-21 2002-04-30 Micro Motion, Inc. System for calibrating a drive signal in a coriolis flowmeter to cause the driver to vibrate a conduit in a desired mode of vibration
US6505135B2 (en) * 2001-03-13 2003-01-07 Micro Motion, Inc. Initialization algorithm for drive control in a coriolis flowmeter

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008014949A (ja) * 2006-07-04 2008-01-24 Krohne Ag コリオリ質量流量測定装置の作動方法
JP2018508025A (ja) * 2015-03-13 2018-03-22 マイクロ モーション インコーポレイテッド 振動式流量計における信号の温度補償
CN110346006A (zh) * 2019-08-13 2019-10-18 上海一诺仪表有限公司 一种质量流量计变送器和质量流量计

Also Published As

Publication number Publication date
US20020019710A1 (en) 2002-02-14
EP0919793A2 (en) 1999-06-02
US7136761B2 (en) 2006-11-14
US20100107778A1 (en) 2010-05-06
EP0919793A3 (en) 1999-10-06
US6507791B2 (en) 2003-01-14
US6754594B2 (en) 2004-06-22
DE69840766D1 (de) 2009-06-04
US20040031328A1 (en) 2004-02-19
DE98309694T1 (de) 2004-09-30
US20050209794A1 (en) 2005-09-22
US20020038186A1 (en) 2002-03-28
US8000906B2 (en) 2011-08-16
US7571062B2 (en) 2009-08-04
US6311136B1 (en) 2001-10-30
EP0919793B1 (en) 2009-04-22
DE69841903D1 (de) 2010-10-28
US6917887B2 (en) 2005-07-12
US20070124090A1 (en) 2007-05-31
JP2007248479A (ja) 2007-09-27

Similar Documents

Publication Publication Date Title
JP2000018995A (ja) ディジタル流量計
US9200936B2 (en) Digital flowmeter
US6758102B2 (en) Correcting for two-phase flow in a digital flowmeter
US7865318B2 (en) Meter electronics and methods for verification diagnostics for a flow meter
US7124646B2 (en) Correcting for two-phase flow in a digital flowmeter
US7784360B2 (en) Correcting for two-phase flow in a digital flowmeter
US9014997B2 (en) Drive techniques for a digital flowmeter
EP1484585B1 (en) Coriolis flowmeter with digital control system
RU2453816C2 (ru) Измерение влажного газа
DE FR Coriolisdurchflussmesser mit digitalem Regelsystem Débitmètre à effet Coriolis avec système de contrôle numérique

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060809

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20061102

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20061108

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070208

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20070307

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070604

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20070711

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20070803

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20100107

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20100113